Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Global patterns in coronavirus diversity</title>
<author><name sortKey="Anthony, Simon J" sort="Anthony, Simon J" uniqKey="Anthony S" first="Simon J." last="Anthony">Simon J. Anthony</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Johnson, Christine K" sort="Johnson, Christine K" uniqKey="Johnson C" first="Christine K." last="Johnson">Christine K. Johnson</name>
<affiliation><nlm:aff id="vex012-aff4">One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Greig, Denise J" sort="Greig, Denise J" uniqKey="Greig D" first="Denise J." last="Greig">Denise J. Greig</name>
<affiliation><nlm:aff id="vex012-aff4">One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kramer, Sarah" sort="Kramer, Sarah" uniqKey="Kramer S" first="Sarah" last="Kramer">Sarah Kramer</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Che, Xiaoyu" sort="Che, Xiaoyu" uniqKey="Che X" first="Xiaoyu" last="Che">Xiaoyu Che</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Wells, Heather" sort="Wells, Heather" uniqKey="Wells H" first="Heather" last="Wells">Heather Wells</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hicks, Allison L" sort="Hicks, Allison L" uniqKey="Hicks A" first="Allison L." last="Hicks">Allison L. Hicks</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Joly, Damien O" sort="Joly, Damien O" uniqKey="Joly D" first="Damien O." last="Joly">Damien O. Joly</name>
<affiliation><nlm:aff id="vex012-aff6">Metabiota, Inc. One Sutter, Suite 600, San Francisco, CA 94104, USA</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="vex012-aff7">Wildlife Conservation Society, New York, NY 10460, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Wolfe, Nathan D" sort="Wolfe, Nathan D" uniqKey="Wolfe N" first="Nathan D." last="Wolfe">Nathan D. Wolfe</name>
<affiliation><nlm:aff id="vex012-aff6">Metabiota, Inc. One Sutter, Suite 600, San Francisco, CA 94104, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Daszak, Peter" sort="Daszak, Peter" uniqKey="Daszak P" first="Peter" last="Daszak">Peter Daszak</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Karesh, William" sort="Karesh, William" uniqKey="Karesh W" first="William" last="Karesh">William Karesh</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Lipkin, W I" sort="Lipkin, W I" uniqKey="Lipkin W" first="W. I." last="Lipkin">W. I. Lipkin</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Morse, Stephen S" sort="Morse, Stephen S" uniqKey="Morse S" first="Stephen S." last="Morse">Stephen S. Morse</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mazet, Jonna A K" sort="Mazet, Jonna A K" uniqKey="Mazet J" first="Jonna A. K." last="Mazet">Jonna A. K. Mazet</name>
<affiliation><nlm:aff id="vex012-aff4">One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Goldstein, Tracey" sort="Goldstein, Tracey" uniqKey="Goldstein T" first="Tracey" last="Goldstein">Tracey Goldstein</name>
<affiliation><nlm:aff id="vex012-aff4">One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">28630747</idno>
<idno type="pmc">5467638</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5467638</idno>
<idno type="RBID">PMC:5467638</idno>
<idno type="doi">10.1093/ve/vex012</idno>
<date when="2017">2017</date>
<idno type="wicri:Area/Pmc/Corpus">001581</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">001581</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Global patterns in coronavirus diversity</title>
<author><name sortKey="Anthony, Simon J" sort="Anthony, Simon J" uniqKey="Anthony S" first="Simon J." last="Anthony">Simon J. Anthony</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Johnson, Christine K" sort="Johnson, Christine K" uniqKey="Johnson C" first="Christine K." last="Johnson">Christine K. Johnson</name>
<affiliation><nlm:aff id="vex012-aff4">One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Greig, Denise J" sort="Greig, Denise J" uniqKey="Greig D" first="Denise J." last="Greig">Denise J. Greig</name>
<affiliation><nlm:aff id="vex012-aff4">One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Kramer, Sarah" sort="Kramer, Sarah" uniqKey="Kramer S" first="Sarah" last="Kramer">Sarah Kramer</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Che, Xiaoyu" sort="Che, Xiaoyu" uniqKey="Che X" first="Xiaoyu" last="Che">Xiaoyu Che</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Wells, Heather" sort="Wells, Heather" uniqKey="Wells H" first="Heather" last="Wells">Heather Wells</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Hicks, Allison L" sort="Hicks, Allison L" uniqKey="Hicks A" first="Allison L." last="Hicks">Allison L. Hicks</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Joly, Damien O" sort="Joly, Damien O" uniqKey="Joly D" first="Damien O." last="Joly">Damien O. Joly</name>
<affiliation><nlm:aff id="vex012-aff6">Metabiota, Inc. One Sutter, Suite 600, San Francisco, CA 94104, USA</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="vex012-aff7">Wildlife Conservation Society, New York, NY 10460, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Wolfe, Nathan D" sort="Wolfe, Nathan D" uniqKey="Wolfe N" first="Nathan D." last="Wolfe">Nathan D. Wolfe</name>
<affiliation><nlm:aff id="vex012-aff6">Metabiota, Inc. One Sutter, Suite 600, San Francisco, CA 94104, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Daszak, Peter" sort="Daszak, Peter" uniqKey="Daszak P" first="Peter" last="Daszak">Peter Daszak</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Karesh, William" sort="Karesh, William" uniqKey="Karesh W" first="William" last="Karesh">William Karesh</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Lipkin, W I" sort="Lipkin, W I" uniqKey="Lipkin W" first="W. I." last="Lipkin">W. I. Lipkin</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Morse, Stephen S" sort="Morse, Stephen S" uniqKey="Morse S" first="Stephen S." last="Morse">Stephen S. Morse</name>
<affiliation><nlm:aff>NONE</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Mazet, Jonna A K" sort="Mazet, Jonna A K" uniqKey="Mazet J" first="Jonna A. K." last="Mazet">Jonna A. K. Mazet</name>
<affiliation><nlm:aff id="vex012-aff4">One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Goldstein, Tracey" sort="Goldstein, Tracey" uniqKey="Goldstein T" first="Tracey" last="Goldstein">Tracey Goldstein</name>
<affiliation><nlm:aff id="vex012-aff4">One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Virus Evolution</title>
<idno type="eISSN">2057-1577</idno>
<imprint><date when="2017">2017</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><title>Abstract</title>
<p>Since the emergence of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and Middle East Respiratory Syndrom Coronavirus (MERS-CoV) it has become increasingly clear that bats are important reservoirs of CoVs. Despite this, only 6% of all CoV sequences in GenBank are from bats. The remaining 94% largely consist of known pathogens of public health or agricultural significance, indicating that current research effort is heavily biased towards describing known diseases rather than the ‘pre-emergent’ diversity in bats. Our study addresses this critical gap, and focuses on resource poor countries where the risk of zoonotic emergence is believed to be highest. We surveyed the diversity of CoVs in multiple host taxa from twenty countries to explore the factors driving viral diversity at a global scale. We identified sequences representing 100 discrete phylogenetic clusters, ninety-one of which were found in bats, and used ecological and epidemiologic analyses to show that patterns of CoV diversity correlate with those of bat diversity. This cements bats as the major evolutionary reservoirs and ecological drivers of CoV diversity. Co-phylogenetic reconciliation analysis was also used to show that host switching has contributed to CoV evolution, and a preliminary analysis suggests that regional variation exists in the dynamics of this process. Overall our study represents a model for exploring global viral diversity and advances our fundamental understanding of CoV biodiversity and the potential risk factors associated with zoonotic emergence.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Anthony, S J" uniqKey="Anthony S">S. J. Anthony</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Anthony, S J" uniqKey="Anthony S">S. J. Anthony</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Anthony, S J" uniqKey="Anthony S">S. J. Anthony</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Azhar, E I" uniqKey="Azhar E">E. I. Azhar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Blackwood, J C" uniqKey="Blackwood J">J. C. Blackwood</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Burnham, K P" uniqKey="Burnham K">K. P. Burnham</name>
</author>
<author><name sortKey="Anderson, D R" uniqKey="Anderson D">D. R. Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Charleston, M A" uniqKey="Charleston M">M. A. Charleston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Conow, C" uniqKey="Conow C">C. Conow</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="Coulibaly N Golo, D" uniqKey="Coulibaly N Golo D">D. Coulibaly-N'Golo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Darriba, D" uniqKey="Darriba D">D. Darriba</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Groot, R J" uniqKey="De Groot R">R. J. de Groot</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Donaldson, E F" uniqKey="Donaldson E">E. F. Donaldson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Drexler, J F" uniqKey="Drexler J">J. F. Drexler</name>
</author>
<author><name sortKey="Corman, V M" uniqKey="Corman V">V. M. Corman</name>
</author>
<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="Edgar, R C" uniqKey="Edgar R">R. C. Edgar</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ellis, V A" uniqKey="Ellis V">V. A. Ellis</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Ge, X Y" uniqKey="Ge X">X. Y. Ge</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hijmans, R J" uniqKey="Hijmans R">R. J. Hijmans</name>
</author>
<author><name sortKey="Williams, E" uniqKey="Williams E">E. Williams</name>
</author>
<author><name sortKey="Vennes, C" uniqKey="Vennes C">C. Vennes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hagberg, A A" uniqKey="Hagberg A">A. A. Hagberg</name>
</author>
<author><name sortKey="Schult, D A" uniqKey="Schult D">D. A. Schult</name>
</author>
<author><name sortKey="Swart, P J" uniqKey="Swart P">P. J. Swart</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hu, B" uniqKey="Hu B">B. Hu</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huang, Y W" uniqKey="Huang Y">Y. W. Huang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Huynh, J" uniqKey="Huynh J">J. Huynh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Irwin, N R" uniqKey="Irwin N">N. R. Irwin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jacomy, M" uniqKey="Jacomy M">M. Jacomy</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jones, K E" uniqKey="Jones K">K. E. Jones</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jost, L" uniqKey="Jost L">L. Jost</name>
</author>
<author><name sortKey="Chao, A" uniqKey="Chao A">A. Chao</name>
</author>
<author><name sortKey="Chazdon, R L" uniqKey="Chazdon R">R. L. Chazdon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Karesh, W B" uniqKey="Karesh W">W. B. Karesh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Keesing, F" uniqKey="Keesing F">F. Keesing</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="King, A M Q" uniqKey="King A">A. M. Q. King</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kreuder Johnson, C" uniqKey="Kreuder Johnson C">C. Kreuder Johnson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Ksiazek, T G" uniqKey="Ksiazek T">T. G. Ksiazek</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lau, S K" uniqKey="Lau S">S. K. Lau</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lau, S K" uniqKey="Lau S">S. K. Lau</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lau, S K" uniqKey="Lau S">S. K. Lau</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lelli, D" uniqKey="Lelli D">D. Lelli</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Li, W" uniqKey="Li W">W. Li</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Liang, K" uniqKey="Liang K">K. Liang</name>
</author>
<author><name sortKey="Zeger, S L" uniqKey="Zeger S">S. L. Zeger</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Maes, P" uniqKey="Maes P">P. Maes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Magurran, A E" uniqKey="Magurran A">A. E. Magurran</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Melade, J" uniqKey="Melade J">J. Melade</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Merkle, D" uniqKey="Merkle D">D. Merkle</name>
</author>
<author><name sortKey="Middendorf, M" uniqKey="Middendorf M">M. Middendorf</name>
</author>
<author><name sortKey="Wieseke, N" uniqKey="Wieseke N">N. Wieseke</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Morse, S S" uniqKey="Morse S">S. S. Morse</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pfefferle, S" uniqKey="Pfefferle S">S. Pfefferle</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Quan, P L" uniqKey="Quan P">P. L. Quan</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Ramsden, C" uniqKey="Ramsden C">C. Ramsden</name>
</author>
<author><name sortKey="Holmes, E C" uniqKey="Holmes E">E. C. Holmes</name>
</author>
<author><name sortKey="Charleston, M A" uniqKey="Charleston M">M. A. Charleston</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="Rubin, D B" uniqKey="Rubin D">D. B. Rubin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sabir, J S" uniqKey="Sabir J">J. S. Sabir</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Streicker, D G" uniqKey="Streicker D">D. G. Streicker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tamura, K" uniqKey="Tamura K">K. Tamura</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tang, X C" uniqKey="Tang X">X. C. Tang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tsoleridis, T" uniqKey="Tsoleridis T">T. Tsoleridis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vijaykrishna, D" uniqKey="Vijaykrishna D">D. Vijaykrishna</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vijgen, L" uniqKey="Vijgen L">L. Vijgen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wacharapluesadee, S" uniqKey="Wacharapluesadee S">S. Wacharapluesadee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, Q" uniqKey="Wang Q">Q. Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, W" uniqKey="Wang W">W. Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Watanabe, S" uniqKey="Watanabe S">S. Watanabe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Woo, P C" uniqKey="Woo P">P. C. Woo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Woo, P C" uniqKey="Woo P">P. C. Woo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Woo, P C" uniqKey="Woo P">P. C. Woo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Woo, P C" uniqKey="Woo P">P. C. Woo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yang, Y" uniqKey="Yang Y">Y. 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-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">28630747</article-id>
<article-id pub-id-type="pmc">5467638</article-id>
<article-id pub-id-type="doi">10.1093/ve/vex012</article-id>
<article-id pub-id-type="publisher-id">vex012</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>Global patterns in coronavirus diversity</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Anthony</surname>
<given-names>Simon J.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff1">1 th</xref>
<xref ref-type="aff" rid="vex012-aff2">2 th</xref>
<xref ref-type="aff" rid="vex012-aff3">3 th</xref>
<xref ref-type="corresp" rid="vex012-cor1"></xref>
<pmc-comment>sja2127@cumc.columbia.edu </pmc-comment>
</contrib>
<contrib contrib-type="author"><name><surname>Johnson</surname>
<given-names>Christine K.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff4">4</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Greig</surname>
<given-names>Denise J.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff4">4</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Kramer</surname>
<given-names>Sarah</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff1">1 th</xref>
<xref ref-type="aff" rid="vex012-aff5">5 th</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Che</surname>
<given-names>Xiaoyu</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff1">1 th</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Wells</surname>
<given-names>Heather</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff1">1 th</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Hicks</surname>
<given-names>Allison L.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff1">1 th</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Joly</surname>
<given-names>Damien O.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff6">6</xref>
<xref ref-type="aff" rid="vex012-aff7">7</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Wolfe</surname>
<given-names>Nathan D.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff6">6</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Daszak</surname>
<given-names>Peter</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff3">3 th</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Karesh</surname>
<given-names>William</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff3">3 th</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Lipkin</surname>
<given-names>W. I.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff1">1 th</xref>
<xref ref-type="aff" rid="vex012-aff2">2 th</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Morse</surname>
<given-names>Stephen S.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff2">2 th</xref>
</contrib>
<contrib contrib-type="author"><collab>PREDICT Consortium</collab>
<xref ref-type="aff" rid="vex012-aff8">8</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Mazet</surname>
<given-names>Jonna A. K.</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff4">4</xref>
<xref ref-type="corresp" rid="vex012-FM1"></xref>
</contrib>
<contrib contrib-type="author"><name><surname>Goldstein</surname>
<given-names>Tracey</given-names>
</name>
<xref ref-type="aff" rid="vex012-aff4">4</xref>
<xref ref-type="corresp" rid="vex012-cor1"></xref>
<xref ref-type="corresp" rid="vex012-FM1"></xref>
<pmc-comment>tgoldstein@ucdavis.edu </pmc-comment>
</contrib>
</contrib-group>
<aff id="vex012-aff1"><label>1</label>
Center for Infection and Immunity, Mailman School of Public Health, Columbia University, 722 West 168<label>th</label>
Street, New York, NY 10032, USA</aff>
<aff id="vex012-aff2"><label>2</label>
Department of Epidemiology, Mailman School of Public Health, Columbia University, 722 West 168<label>th</label>
Street, New York, NY 10032, USA</aff>
<aff id="vex012-aff3"><label>3</label>
EcoHealth Alliance, 460 West 34<label>th</label>
Street, New York, NY 10001, USA</aff>
<aff id="vex012-aff4"><label>4</label>
One Health Institute & Karen C Drayer Wildlife Health Center, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA</aff>
<aff id="vex012-aff5"><label>5</label>
Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, 722 West 168<label>th</label>
Street, New York, NY 10032, USA</aff>
<aff id="vex012-aff6"><label>6</label>
Metabiota, Inc. One Sutter, Suite 600, San Francisco, CA 94104, USA</aff>
<aff id="vex012-aff7"><label>7</label>
Wildlife Conservation Society, New York, NY 10460, USA</aff>
<aff id="vex012-aff8"><label>8</label>
<ext-link ext-link-type="uri" xlink:href="http://www.vetmed.ucdavis.edu/ohi/predict/publications/Authorship.cfm">http://www.vetmed.ucdavis.edu/ohi/predict/publications/Authorship.cfm</ext-link>
</aff>
<author-notes><corresp id="vex012-cor1"><label>*</label>
Corresponding author: E-mail: <email>sja2127@cumc.columbia.edu</email>
; <email>tgoldstein@ucdavis.edu</email>
</corresp>
<corresp id="vex012-FM1"><label>†</label>
Joint senior author.</corresp>
</author-notes>
<pub-date pub-type="collection"><month>1</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="epub" iso-8601-date="2017-06-12"><day>12</day>
<month>6</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>12</day>
<month>6</month>
<year>2017</year>
</pub-date>
<pmc-comment> PMC Release delay is 0 months and 0 days and was based on the . </pmc-comment>
<volume>3</volume>
<issue>1</issue>
<elocation-id>vex012</elocation-id>
<permissions><copyright-statement>© The Author 2017. Published by Oxford University Press.</copyright-statement>
<copyright-year>2017</copyright-year>
<license xlink:href="http://creativecommons.org/licenses/by-nc/4.0/" license-type="cc-by-nc"><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>
<self-uri xlink:href="vex012.pdf"></self-uri>
<abstract><title>Abstract</title>
<p>Since the emergence of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and Middle East Respiratory Syndrom Coronavirus (MERS-CoV) it has become increasingly clear that bats are important reservoirs of CoVs. Despite this, only 6% of all CoV sequences in GenBank are from bats. The remaining 94% largely consist of known pathogens of public health or agricultural significance, indicating that current research effort is heavily biased towards describing known diseases rather than the ‘pre-emergent’ diversity in bats. Our study addresses this critical gap, and focuses on resource poor countries where the risk of zoonotic emergence is believed to be highest. We surveyed the diversity of CoVs in multiple host taxa from twenty countries to explore the factors driving viral diversity at a global scale. We identified sequences representing 100 discrete phylogenetic clusters, ninety-one of which were found in bats, and used ecological and epidemiologic analyses to show that patterns of CoV diversity correlate with those of bat diversity. This cements bats as the major evolutionary reservoirs and ecological drivers of CoV diversity. Co-phylogenetic reconciliation analysis was also used to show that host switching has contributed to CoV evolution, and a preliminary analysis suggests that regional variation exists in the dynamics of this process. Overall our study represents a model for exploring global viral diversity and advances our fundamental understanding of CoV biodiversity and the potential risk factors associated with zoonotic emergence.</p>
</abstract>
<kwd-group><kwd>coronavirus</kwd>
<kwd>viral ecology</kwd>
<kwd>bat</kwd>
<kwd>evolution</kwd>
</kwd-group>
<funding-group><award-group award-type="grant"><funding-source><named-content content-type="funder-name">United States Agency for International Development</named-content>
</funding-source>
</award-group>
</funding-group>
<counts><page-count count="15"></page-count>
</counts>
</article-meta>
</front>
<body><sec><title>1. Introduction</title>
<p>In 2002/3, SARS-Coronavirus (SARS-CoV) emerged in the Guangdong province of southern China (<xref rid="vex012-B16" ref-type="bibr">Drosten et al. 2003</xref>
; <xref rid="vex012-B34" ref-type="bibr">Ksiazek et al. 2003</xref>
). It quickly spread to twenty-seven countries, infecting 8,098 people with 774 deaths, and was declared the first global pandemic of the 21st century. Bats were identified as the reservoir (<xref rid="vex012-B35" ref-type="bibr">Lau et al. 2005</xref>
; <xref rid="vex012-B39" ref-type="bibr">Li et al. 2005</xref>
) and probable source (<xref rid="vex012-B20" ref-type="bibr">Ge et al. 2013</xref>
) of the outbreak. In 2012, the SARS pandemic was followed by MERS-Coronavirus (MERS-CoV), which emerged in the Middle East (<xref rid="vex012-B70" ref-type="bibr">Zaki et al. 2012</xref>
) with 1,782 confirmed human cases and 640 deaths (as of September 2016). Camels were identified as the likely source of human infections (<xref rid="vex012-B52" ref-type="bibr">Reusken et al. 2013</xref>
; <xref rid="vex012-B4" ref-type="bibr">Azhar et al. 2014</xref>
); however, bats were again found to host closely related (MERS-like) viruses and are therefore assumed to be the original evolutionary source (<xref rid="vex012-B65" ref-type="bibr">Woo et al. 2006</xref>
; <xref rid="vex012-B1" ref-type="bibr">Anthony et al. 2012</xref>
; <xref rid="vex012-B9" ref-type="bibr">Corman et al. 2014a</xref>
,<xref rid="vex012-B10" ref-type="bibr">b</xref>
; <xref rid="vex012-B61" ref-type="bibr">Wacharapluesadee et al. 2015</xref>
; <xref rid="vex012-B3" ref-type="bibr">Anthony et al. 2017</xref>
).</p>
<p>Together, these outbreaks have cemented the coronaviridae as a family of zoonotic concern, and stimulated a surge in viral discovery efforts in bats [reviewed by <xref rid="vex012-B15" ref-type="bibr">Drexler et al. (2014)</xref>
]. These efforts appear to show that almost all human CoVs have zoonotic origins or otherwise circulate in animals, including human 229E [bats (<xref rid="vex012-B46" ref-type="bibr">Pfefferle et al. 2009</xref>
); camels (<xref rid="vex012-B54" ref-type="bibr">Sabir et al. 2016</xref>
)], NL63 [bats (<xref rid="vex012-B14" ref-type="bibr">Donaldson et al. 2010</xref>
; <xref rid="vex012-B25" ref-type="bibr">Huynh et al. 2012</xref>
)], and OC43 [cattle (<xref rid="vex012-B60" ref-type="bibr">Vijgen et al. 2005</xref>
)]. Even non-human CoVs such as porcine epidemic diarrhea virus (PEDV) may have emerged by host switching from other animals [bats (<xref rid="vex012-B57" ref-type="bibr">Tang et al. 2006</xref>
; <xref rid="vex012-B24" ref-type="bibr">Huang et al. 2013</xref>
)]. Overall, it seems that CoV diversity in bats is substantial (<xref rid="vex012-B15" ref-type="bibr">Drexler et al. 2014</xref>
), that these viruses are prone to host switching (<xref rid="vex012-B67" ref-type="bibr">Woo et al. 2009</xref>
), and that they are a current, historic, and future threat to public health.</p>
<p>Despite recent efforts, there are several notable gaps in our understanding of CoV diversity. Foremost, our knowledge of CoVs in resource-limited countries is poor (<xref rid="vex012-B15" ref-type="bibr">Drexler et al. 2014</xref>
). This is particularly problematic given that many of these same areas are predicted to be hotspots of disease emergence (<xref rid="vex012-B28" ref-type="bibr">Jones et al. 2008</xref>
). Second, there has been little effort to understand the evolutionary and ecological drivers of CoV diversity on a global scale, or to evaluate host and regional variation in the factors that contribute to the risk of emergence (e.g. host switching). Indeed, most studies to date have been somewhat limited in their geographic scope and have included only small sample sizes with little epidemiologic or ecological context. In short, there is a critical need for a more global perspective on CoV diversity.</p>
<p>In 2009, the PREDICT project was established, in part, to address this need. Focused on strengthening capacity and identifying viruses in wildlife at high-risk interfaces (<xref rid="vex012-B47" ref-type="bibr">PREDICT_Consortium 2014</xref>
), this USAID Emerging Pandemic Threats (EPT) initiative worked with local partners in twenty countries across Latin America, Africa, and Asia over 5 years to better understand the current diversity of CoVs (as well as other viruses) and evaluate the factors that drive this diversity at different scales. Herein, we report a large diversity of CoV sequences (mostly from bats), show that the biogeography of bats has shaped the diversity of CoVs globally, and provide evidence to suggest there could be regional variation in host switching and the risk for zoonotic emergence.</p>
</sec>
<sec sec-type="methods"><title>2. Methods</title>
<sec><title>2.1 Animals and samples</title>
<p>Bats (n = 12,333), rodents (n = 3,387), and non-human primates (n = 3,470) were humanely sampled (capture and release) from twenty ‘hotspot’ countries (<xref rid="vex012-B28" ref-type="bibr">Jones et al. 2008</xref>
), representing central Africa (Cameroon, Gabon, Democratic Republic of Congo, Republic of Congo, Rwanda, Tanzania, Uganda), Latin America (Peru, Bolivia, Brazil, Mexico), and Asia (Bangladesh, Cambodia, China, Indonesia, Laos, Malaysia, Nepal, Thailand, Viet Nam). Samples were also collected from humans (n = 1,124) in seven countries in central Africa and Asia as a pilot to begin to explore the propensity for viral sharing with wildlife. In all twenty countries, wildlife samples were collected from ‘high risk’ interfaces, where direct or indirect contact with humans might promote zoonotic viral transmission. The selected sampling sites included areas of land-use change (deforestation, conversion to agriculture); sites in and around human dwellings; foci of ecotourism; markets restaurants and farms along the animal value chain; and areas where occupational exposure was likely (animal sanctuaries, agricultural activities). When possible, individuals were identified to lowest taxonomic order (genus and species) and assigned to an age class (adult, subadult, neonate) by the field teams. All samples from swabs (e.g. oral, urine, rectal), fluids (e.g. saliva), and tissues were collected into (1) NucliSens<sup>®</sup>
Lysis Buffer (bioMérieux, Inc., Marcy-I’Étoile, France) and (2) viral transport media, and then frozen in the field in liquid nitrogen and transferred to the laboratory for storage at −80 °C. All animal sampling activities were conducted with permissions from local authorities and under the Institutional Animal Care and Use Committee at the University of California, Davis (protocol number: 16048). All human activities were reviewed and approved by the UC Davis Institutional Review Board (IRB), under protocols: 215253 and 432330.</p>
</sec>
<sec><title>2.2 Viral discovery</title>
<p>RNA was extracted from all samples, and cDNA prepared using superscript III (Invitrogen). Two broadly reactive consensus PCR assays targeting non-overlapping fragments of the orf1ab were used to detect both known and novel CoVs (<xref rid="vex012-B48" ref-type="bibr">Quan et al. 2010</xref>
; <xref rid="vex012-B64" ref-type="bibr">Watanabe et al. 2010</xref>
). The first [the ‘Watanabe’ assay (<xref rid="vex012-B64" ref-type="bibr">Watanabe et al. 2010</xref>
)] amplified a ∼434 bp fragment of the RNA-dependent RNA polymerase (RdRp) corresponding to nucleotides (NTs) 15,156–15,589 in the human CoV OC43 genome (NC_005147), while the second [the ‘Quan’ assay (<xref rid="vex012-B48" ref-type="bibr">Quan et al. 2010</xref>
)] amplified a ∼332 bp fragment of a different peptide downstream of the RdRp, corresponding to NTs 18,323–18,654. Amplified products of the expected size were cloned and sequenced (traditional Sanger dideoxy sequencing) according to standard protocols, and sequences edited manually in Geneious Pro (ver 9.1.3, Biomatters, Auckland, NZ). We note that this approach was adopted to facilitate viral discovery in resource-poor settings, which are the target of pandemic preparedness initiatives such as USAID PREDICT, and where techniques such as high throughput sequencing are largely unavailable.</p>
</sec>
<sec><title>2.3 Phylogenetics</title>
<p>Coronavirus sequences from the Quan and Watanabe datasets were aligned with reference sequences collected from GenBank using MUSCLE (<xref rid="vex012-B17" ref-type="bibr">Edgar 2004</xref>
), followed by manual alignment in Se-Al v2.0a11 (<ext-link ext-link-type="uri" xlink:href="http://tree.bio.ed.ac.uk">http://tree.bio.ed.ac.uk</ext-link>
). The best-fitting model of nucleotide substitution for each dataset was selected by Akaike’s Information Criterion (AIC) using jModeltest v2.1.5 (<xref rid="vex012-B12" ref-type="bibr">Darriba et al. 2012</xref>
). For both alignments, the general time reversible model with a gamma distribution of rate heterogeneity was the best-fitting model. Maximum likelihood trees were constructed from each alignment using MEGA v6.06 (<xref rid="vex012-B56" ref-type="bibr">Tamura et al. 2013</xref>
). Histograms of all pairwise sequence identities (%) were plotted, as described previously for hantaviruses (<xref rid="vex012-B41" ref-type="bibr">Maes et al. 2009</xref>
), and used to define cut-offs between viral sequence clusters for analysis (akin to operating taxonomic units).</p>
</sec>
<sec><title>2.4 Diversity measures</title>
<p>The Earth’s surface was divided into grid cells, each 10° latitude by 10° longitude. Thirty-four of these grid cells included areas from which bats were sampled (bats were identified down to the species level in thirty cells), and twenty-seven cells had samples testing positive for CoVs. If the coordinates listed for a particular bat were on the line between grid cells, it was arbitrarily included in the grid cell with higher latitude/longitude. Viral and host species richness for each cell was calculated by determining the total number of unique virus and bat species within each cell. Alpha diversity was calculated using the Shannon index (<xref rid="vex012-B29" ref-type="bibr">Jost et al. 2011</xref>
), which was chosen over the Gini–Simpson index because it does not place disproportionate weight on dominant species. Since we do not know how representative our sampling was, the Shannon index was deemed the best approach. The effective number of viruses and bats in each cell was then correlated using Kendall’s tau, since the values were not normally distributed. Sensitivity analyses were performed by shifting the grid cells in 1° increments and re-calculating richness and effective species numbers.</p>
<p>Matrices of beta diversity between each combination of grid cells were developed using the Jaccard index (<xref rid="vex012-B29" ref-type="bibr">Jost et al. 2011</xref>
), and Mantel tests were used to determine whether the distance between grids was associated with the dissimilarity in virus and bat species. Geographic distance matrices were calculated by finding the centroid of each grid cell and calculating the geodesic distance between each pair using the R package ‘geosphere’(<xref rid="vex012-B21" ref-type="bibr">Geosphere 2015</xref>
) and tests were performed using Spearman’s rank correlation. All analyses of diversity were conducted in R version 2.3.2.</p>
</sec>
<sec><title>2.5 Network model</title>
<p>A presence/absence matrix was constructed to show the distribution of viral sequence clusters across all species. Using the Python package NetworkX (<xref rid="vex012-B22" ref-type="bibr">Hagberg et al. 2008</xref>
), a network model was constructed, connecting bat species to all viral clusters that were identified within that species in our data. The network is made up of thirty-one connected components, eighty-two viral clusters, eighty-five bat species, and 159 edges. Networks were plotted using Gephi using the force-directed algorithm ForceAtlas2 (<xref rid="vex012-B27" ref-type="bibr">Jacomy et al. 2014</xref>
). Specifically, we plotted two bipartite graphs, one where hosts are colored by region and one in which they are colored by family.</p>
</sec>
<sec><title>2.6 Cophylogeny</title>
<p>The Jane (<xref rid="vex012-B8" ref-type="bibr">Conow et al. 2010</xref>
) software tool was used for cophylogenetic reconstructions. This approach applies an <italic>a priori</italic>
‘cost scheme’ to different evolutionary events in order to test the degree of congruence between two trees (virus and host). These ‘events’ include cospeciation, duplication, host switching, and failure to diverge (or sorting) (<xref rid="vex012-B7" ref-type="bibr">Charleston 1998</xref>
; <xref rid="vex012-B8" ref-type="bibr">Conow et al. 2010</xref>
), which are used to create a minimal-cost reconstruction of the evolutionary history between the viruses and hosts. Another cophylogeny program, CoRe-PA (<xref rid="vex012-B44" ref-type="bibr">Merkle et al. 2010</xref>
), uses a parameter-adaptive approach to estimate an appropriate cost scheme without any prior value assignment. For our analysis, we used Jane with the default pre-determined cost scheme where cospeciation is the null hypothesis and, therefore, the ‘cheapest’ event, with the assumption that host switching would theoretically be more costly than genetic drift within a species to which a virus is already adapted. Analysis was repeated with varied cost values in Jane to test for sensitivity to our cost scheme, as well as in CoRe-PA for robustness, and event types and costs from each reconstruction were compared between the two programs. As these changes ultimately did not significantly influence the outcome of the analysis, exact event assignments used in the final analysis were inferred solely from Jane with the default cost scheme.</p>
<p>The analyses were performed using all unique associations between viruses and their respective hosts (limited to bats). Alpha-CoV sequences and beta-CoV sequences were analyzed separately, since each appears to have diversified independently within bats. This gave us a total of four reconstructions (alpha-Quan, alpha-Watanabe, beta-Quan, beta-Watanabe). Only one instance of each virus–host association was included, even if multiple instances of the same association were observed (i.e. our analysis does not account for the frequency of detection). Associations were excluded if the host was not identified to the species level or if the cytochrome B sequence did not exist in GenBank or could not be amplified by PCR. The topology of the virus tree was inferred from the topology of the full tree for consistency, since subsets of sequences often produce different arrangements.</p>
<p>For each reconstruction, the most recent evolutionary event leading to a virus–host association was recorded. Given that sorting is particularly sensitive to missing data (it indicates a loss on the tree when in fact it may simply be an artifact of under-sampling), we have excluded this event from our analysis. We also distinguished between two potential types of host switching, based on an assumption that there could be differences in zoonotic risk between viruses that only move between closely related bats (e.g. those within the same genus) and those that are able to make more substantive jumps into distantly related species (those in other host genera or families). We therefore use ‘host switching’ to refer specifically to viral lineages that have moved into a host belonging to another genus or family, and ‘sharing’ to indicate a viral lineage that has moved into a different host species of the same genus (without speciating). We selected host genus as the cut-off to be consistent with current conventions of taxonomic hierarchy, and verified that the number of species per genus is relatively consistent between regions. We further verified that the mean genetic distance between host species is consistent between regions (i.e. that the taxonomic distance a virus has to navigate when moving between any two species is roughly similar in all three regions).</p>
<p>A binary regression model with logistic link function was used to investigate whether the degree of host switching varied by region. The dependent variable was host switching, and either cospeciation or sharing was the reference group. Region was the independent variable. We used the generalized estimating equations (GEE) (<xref rid="vex012-B40" ref-type="bibr">Liang and Zeger 1986</xref>
) method in order to account for multiple event types within one family and derive a robust estimate of the standard errors. Event data for all four reconstructions were aggregated for the analysis. Where host–virus associations were duplicated, one was removed. If the duplicated event types did not agree [which can occur given the sensitivity of these methods to variation in topology and number of taxa in each tree (<xref rid="vex012-B8" ref-type="bibr">Conow et al. 2010</xref>
)], we randomly selected one of the two possibilities and repeated the randomization 100 times. The GEE estimates from the 100 random iterations were then combined using rules established by <xref rid="vex012-B53" ref-type="bibr">Rubin (1987)</xref>
. The analysis was repeated to explore variation based on bat family (independent variable). Some families only possessed one type of event, which would introduce quasi-complete separation into the GEE model, so they were excluded from the corresponding models. However, the numbers of bat species within these excluded families were very small (<5) compared to the overall number of bat species for the whole study, so we felt this would not impact any overall trends.</p>
</sec>
<sec><title>2.7 Viral discovery effort</title>
<p>The relationship between sampling effort and viral discovery was evaluated using a log-link Poisson regression model, fitted with the count of viruses for each species as the dependent variable and the number of animals sampled as the independent variable. This model is universally applied with count data. The coefficient of the sample number variable is the log of the ratio of the expected number of viruses with (n + 1) samples collected to the expected number of viruses with n samples collected, where n can be any arbitrary positive integer. Using this estimated coefficient, we calculated the estimated numbers of viruses (viral sequence clusters) that can be found with respect to the sampling effort (i.e. all possible numbers of collected animals within a species), and they were then used to plot the estimated curve together with its 95% confidence band. We clarify that we use the term ‘expected’ to indicate the statistical expectation based on the estimate of the regression model. We acknowledge that the model fit is not as good as the model with an additional quadratic term of the independent variable, nonetheless it avoids the problem of over-fitting and realistically reflects the relationship between sampling effort and viral discovery when the number of animals collected is relatively lower.</p>
</sec>
<sec><title>2.8 Factors associated with CoV positivity</title>
<p>Due to the small numbers of positive individuals obtained for rodents, non-human primates, and humans, only bat data were analyzed for factors associated with positivity. A total of 12,333 bats were tested. This included 5,624 males and 5,767 females (922 were not sexed). Among aged bats 7,385 were adults, 1,029 were sub-adults, and 8 were neonates (3,911 were not aged). Binary logistic regression models were used to evaluate significant variables, restricting analyses to species for which ≥50 individuals were tested. Model selection was first performed in R (version 3.2.5) (<xref rid="vex012-B49" ref-type="bibr">R Foundation for Statistical Computing, 2008</xref>
) using the R package ‘gmulti’ and the step function, and the best model determined by AIC. A robust estimation of standard error was calculated by clustering by individual Animal Identification number to account for non-independence of multiple tests from the same animal using STATA 13.1, SE (College Station, TX, USA). The independent variables included in analyses were Specimen Type (blood, feces/rectal swabs, guano, oral/rectal sample, oral/nasal swabs, tissue, urine/urogenital swabs), Host Taxon (family, sub-family, genus), Host Age (adult, sub-adult, neonate), Season (wet, dry), and Interface (animal use, human activity, land use, pristine area). Season was designated as ‘wet’ or ‘dry’ according to month and proximity to the equator for each country. Interfaces were broadly grouped according to the type and intensity of animal contact with people, specifically (1) Land Use (animals sampled in areas with crops, extractive industries, livestock activities), (2) Animal Use (hunting, markets, restaurants, trade, wild animal farms, wildlife managements, zoos/sanctuaries, handling by veterinarian/researchers), (3) Human Activities (ecotourism, in and around human dwellings), and (4) Pristine (where animal human contact was not likely). The model fit for each region was evaluated using the Hosmer–Lemeshow goodness-of-fit test and pseudo R<sup>2</sup>
values.</p>
</sec>
<sec><title>2.9 Virus ‘Hotspot’ Maps</title>
<p>The number of viral sequence clusters (richness) in each sub-clade was summarized by host family (<xref ref-type="supplementary-material" rid="sup1">Supplementary Table S1</xref>
), and chi square tests used to evaluate whether viral sub-clade and host family were independent of each other. Bat species richness maps were then created to infer hotspots of CoV diversity. For each viral clade, spatial distribution data on all bat species belonging to associated families was obtained from the International Union for Conservation of Nature (IUCN). Using ArcGIS version 10 (<xref rid="vex012-B19" ref-type="bibr">ESRI 2011</xref>
), bat species were quantified by counting overlapping polygons. The resulting count data were then converted to a raster and plotted. All sampled bats, as well as all identified sequence clusters within the clade in question, were then plotted over the raster data.</p>
</sec>
</sec>
<sec><title>3. Results</title>
<sec><title>3.1 Global diversity of CoVs</title>
<p>A total of 19,192 animals and humans were assayed for the presence of CoV by consensus PCR (cPCR) (<xref rid="vex012-T1" ref-type="table">Table 1</xref>
). The majority were bats (n = 12,333), representing 282 species from twelve families. Overall, the proportion of CoV positive individuals was 8.6% in bats (n = 1,065/12,333) and 0.2% in non-bats (n = 17/6,859). In other words, over 98% of all positive individuals were bats.
<table-wrap id="vex012-T1" orientation="portrait" position="float"><label>Table 1.</label>
<caption><p> Summary of individuals tested and number positive for at least one coronavirus by host taxa.</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th rowspan="1" colspan="1">Host taxa tested</th>
<th rowspan="1" colspan="1">No. individuals tested</th>
<th rowspan="1" colspan="1">No. individuals positive</th>
<th rowspan="1" colspan="1">No. distinct viruses detected</th>
</tr>
</thead>
<tbody><tr><td rowspan="1" colspan="1">Bats</td>
<td rowspan="1" colspan="1">12,333</td>
<td rowspan="1" colspan="1">1,065</td>
<td rowspan="1" colspan="1">91</td>
</tr>
<tr><td rowspan="1" colspan="1">Non-Human Primates</td>
<td rowspan="1" colspan="1">3,470</td>
<td rowspan="1" colspan="1">4</td>
<td rowspan="1" colspan="1">2</td>
</tr>
<tr><td rowspan="1" colspan="1">Rodents and Shrews</td>
<td rowspan="1" colspan="1">3,387</td>
<td rowspan="1" colspan="1">11</td>
<td rowspan="1" colspan="1">7</td>
</tr>
<tr><td rowspan="1" colspan="1">Humans</td>
<td rowspan="1" colspan="1">1,124</td>
<td rowspan="1" colspan="1">2</td>
<td rowspan="1" colspan="1">2</td>
</tr>
<tr><td rowspan="1" colspan="1">Total</td>
<td rowspan="1" colspan="1">19,192</td>
<td rowspan="1" colspan="1">1,082</td>
<td rowspan="1" colspan="1">100<xref ref-type="fn" rid="tblfn1"><sup>a</sup>
</xref>
</td>
</tr>
</tbody>
</table>
<table-wrap-foot><fn id="tblfn1"><p><sup>a</sup>
Note: Numbers do not total as two viruses were detected in two taxa.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</p>
<p>Partial sequences were obtained from two non-overlapping fragments of the orf1ab gene, yielding 654 sequences for the ‘Quan’ region and 950 for the ‘Watanabe’ region (see ‘Methods’). There was a 27% overlap, where sequences were obtained for both regions from the same sample. The distribution of pairwise sequence identities defined a 90% cut-off between taxa (<xref ref-type="fig" rid="vex012-F1">Fig. 1</xref>
), and the resulting monophyletic groups (in which all sequences had ≥90% identity) used as our operating taxonomic units. Groups that shared less than 90% identity to a known sequence were labeled sequentially as PREDICT_CoV-1, -2, -3 etc; while groups sharing ≥90% identity to a sequence already in GenBank were considered to be strains of a known virus and assigned the same name as the matching sequence (e.g. Kenya_bat_CoV_KY33 or HKU-9). Sequences that shared >90% identity but were found in different hosts were considered part of the same taxonomic unit. Based on these criteria, 100 discrete viral taxa were identified, ninety-one of which were found in bats (<xref rid="vex012-T1" ref-type="table">Table 1</xref>
; <xref ref-type="fig" rid="vex012-F2">Fig. 2</xref>
). Importantly, we make no claims that these groups correspond to ‘species’ as full orf1ab replicase sequences would be required for that (<xref rid="vex012-B32" ref-type="bibr">King et al. 2012</xref>
). Instead we clarify that we are using these partial fragments to cluster sequences into ‘operational taxonomic units’ (analogous to OTU clustering in microbiome research). We also stress that our cut-off was determined based on two discrete regions within the orf1ab separated by >3,000 nucleotides, and that these regions correspond to unique peptides post-translational cleavage.</p>
<fig id="vex012-F1" orientation="portrait" position="float"><label>Figure 1.</label>
<caption><p> Histogram of the relative frequency of pairwise sequence identities used to define the cutoff between operating taxonomic units for all CoV sequences detected. A bimodal distribution was observed and a cutoff of 90% sequence identity used to separate sequences from both the Watanabe (Panel A) and Quan (Panel B) assays into discrete viral taxa.</p>
</caption>
<graphic xlink:href="vex012f1"></graphic>
</fig>
<fig id="vex012-F2" orientation="portrait" position="float"><label>Figure 2.</label>
<caption><p> Maximum likelihood phylogenetic reconstructions for all partial CoV RdRp fragments for both the Quan (Panel A) and Watanabe (Panel B) assays. Sequences are collapsed into clades, representing our operating taxonomic units (sequences sharing ≥90% identity) and number of sequences for each taxon is indicated in parentheses. Representative published sequences from GenBank have been included for comparison (GenBank accession number indicated). Both trees were rooted using the related Breda virus (NC_007447) and all nodes have ≥60% bootstrap support. Pie charts indicate the distribution of viral taxa by host (bat) family, in each virus sub-clade.</p>
</caption>
<graphic xlink:href="vex012f2"></graphic>
</fig>
<p>Of note was the detection of subgroup 2a CoV sequences in bats, humans and non-human primates from four different countries. This sub-clade is largely considered the rodent sub-clade (<xref rid="vex012-B37" ref-type="bibr">Lau et al. 2015</xref>
; <xref rid="vex012-B63" ref-type="bibr">Wang et al. 2015</xref>
; <xref rid="vex012-B58" ref-type="bibr">Tsoleridis et al. 2016</xref>
), but here we detected sequences corresponding to the known virus betacoronavirus-1 in <italic>Pteropus medius</italic>
from Bangladesh, <italic>Pteropus alecto</italic>
from Indonesia, and from humans in China, and sequences corresponding to murine coronavirus in non-human primates from Nepal (these sequences were detected by four different laboratories, and no rodent samples were processed at the same time) (<xref ref-type="fig" rid="vex012-F2">Fig. 2</xref>
). In addition, we report the discovery of several sub-clade 2b and 2c CoV sequence clusters (the SARS and MERS sub-clades, respectively) (<xref ref-type="fig" rid="vex012-F2">Fig. 2</xref>
), and the detection of human coronavirus 229E-like sequences in hipposideros and rhinolophus bats sampled in ROC, Uganda, Cameroon and Gabon, supporting previous suggestions that 229E has zoonotic origins (<xref rid="vex012-B46" ref-type="bibr">Pfefferle et al. 2009</xref>
; <xref rid="vex012-B23" ref-type="bibr">Hu et al. 2015</xref>
). Finally, we highlight the detection of avian-associated infectious bronchitis virus-like sequences (IBV) in bats and a closely related virus (PREDICT_CoV-49) in non-human primates from Bangladesh, as well as porcine epidemic diarrhea virus (PEDV) in bats (<xref ref-type="fig" rid="vex012-F2">Fig. 2</xref>
). Again, we confirm that no livestock samples were processed in any of these laboratories that might have acted as a source of contamination.</p>
</sec>
<sec><title>3.2 Factors driving viral diversity</title>
<p>Viral α-diversity (the ‘effective number’ of species) was significantly correlated with bat α-diversity (tau = 0.325, <italic>P</italic>
= 0.022) (<xref ref-type="fig" rid="vex012-F3">Fig. 3</xref>
), suggesting that more viruses will be found in regions where bat diversity is higher. This association was maintained when viral richness, rather than effective number of species, was used as the diversity index (tau = 0.421, <italic>P</italic>
< 0.001). Viral and bat β-diversity were also significantly correlated (Mantel test; rho = 0.575, <italic>P</italic>
< 0.001). In both cases, diversity differentiated almost entirely into three discrete communities by region (<xref ref-type="fig" rid="vex012-F3">Fig. 3</xref>
). Only Africa and Asia shared any viral sequence clusters (PREDICT_CoV-35 and HKU9; <xref ref-type="fig" rid="vex012-F4">Fig. 4</xref>
), demonstrating that bats have had a strong biogeographic influence on the global ecology and evolution of CoVs.</p>
<fig id="vex012-F3" orientation="portrait" position="float"><label>Figure 3.</label>
<caption><p> Comparison of viral and bat diversity. The earth’s surface was divided into grid cells by latitude and longitude (10 × 10 degree units) for diversity calculations (Panel A). Grid cells where bats were sampled are numbered in each region. Alpha diversity (Shannon H) for virus (Panel B) and host (Panel C) were correlated, indicating that areas of high bat diversity also have high viral diversity. Darker cells indicate higher alpha diversity (i.e. more viral or host taxa) in each grid cell. Beta diversity was also correlated between virus (Panel D) and host (Panel E), and differentiated into three discrete communities by region—Latin America (grid cells 1–10), Africa (grid cells 11–20), and Asia (grid cells 21–34). Shading indicates that either viruses (in red) or hosts (in blue) are shared between two corresponding grid cells, with darker cells indicating higher pairwise similarity.</p>
</caption>
<graphic xlink:href="vex012f3"></graphic>
</fig>
<fig id="vex012-F4" orientation="portrait" position="float"><label>Figure 4.</label>
<caption><p> Network model showing the connection of CoVs and their hosts. Viral sequence clusters (colored grey) are connected to host species, either by region (Panel A) or family (Panel B). Viral and host and communities separate almost entirely by region; only Africa and Asia are connected by two shared viruses (HKU9 and PREDICT_CoV-35) found in species from both continents. Networks also show that viruses appear to be shared by multiple host families in Africa and Asia, while being more restricted to a single family in Latin America.</p>
</caption>
<graphic xlink:href="vex012f4"></graphic>
</fig>
<p>Co-phylogenetic reconciliation analysis was used to investigate the evolutionary mechanisms driving the virus–host associations observed, for example host switching, viral sharing, or co-speciation (see ‘Methods’). Overall, host switching was the dominant mechanism, followed by co-speciation (<xref ref-type="fig" rid="vex012-F5">Fig. 5</xref>
). When analyzed by region, host switching (inter-genus transmission) remained dominant in Africa and Asia, but not in Latin America. Despite fewer host-switching events in Latin America, there was a concomitant increase in virus sharing (intra-genus transmission). This suggests that viruses are still switching hosts in Latin America but preferentially move between closely related species, while in Africa and Asia viruses cross between more distantly related species. To assess the sensitivity of our results to our chosen method (i.e. use of a pre-defined cost scheme in the program Jane), we repeated the analysis for the Quan alpha-CoV and beta-CoV reconstructions using CoRePA, a parameter-adaptive approach that estimates an appropriate cost scheme without any prior value assignment. Comparing the results, we noted 91% agreement between Jane and CoRePA in alpha-CoV reconstruction (thirty-four events, of which thirty-one agreed), and 100% agreement in the beta-CoV reconstructions (twenty-six events, all of which agreed). We further note that the cost scheme calculated in CoRePA was 0.7 for host switching and 0.1 for co-speciation—supporting the default cost scheme used in Jane (i.e. that host switching should be considered more ‘costly’ than co-speciation).</p>
<fig id="vex012-F5" orientation="portrait" position="float"><label>Figure 5.</label>
<caption><p> Relative proportion of evolutionary events leading to observed virus:host associations for CoVs in bats. Cophylogenetic reconstructions were used to identify each event (<xref ref-type="supplementary-material" rid="sup1">Supplementary Fig. S1</xref>
: Panels A–D), and significance evaluated by region. Across all regions, host switching was the dominant evolutionary event (Panel A). When separated by region, host switching remained dominant in Africa (Panel B) and Asia (Panel C), but not in Latin America (Panel D).</p>
</caption>
<graphic xlink:href="vex012f5"></graphic>
</fig>
<p>A bipartite network model connecting viral sequence clusters to their hosts supports these results, illustrating that bat CoVs were connected to several host families in Africa and Asia while being largely restricted to a single bat family in Latin America (<xref ref-type="fig" rid="vex012-F4">Fig. 4</xref>
, Panel B). Further, host switching was nearly four times more likely than virus sharing in Africa, when compared with Latin America (OR: 3.858; <italic>P</italic>
= 0.040). CoVs in Asia were also more likely to host switch than share when compared with Latin America, however the results were not significant (OR: 2.474; <italic>P</italic>
= 0.143). No particular bat family was associated with increased host switching, but this may reflect small sample sizes when the data are summarized by family.</p>
</sec>
<sec><title>3.3 Estimated number of CoVs in bats</title>
<p>A large number of host (bat) species in our study were negative for CoV (n = 197), however none of these species were sampled extensively (<xref ref-type="fig" rid="vex012-F6">Fig. 6</xref>
). We found that all species with sample sizes >110 individuals were positive for one or more CoVs, suggesting that we would have detected CoVs in some of the negative species in our study if sampling effort were increased. Due to the high number of these negative species and their effect on the average number of virus sequence clusters per species, they were excluded from our estimates. By retaining only those species for which >110 individuals were sampled (there were twenty-seven species that qualified), we estimated the average number of CoVs per species to be 2.67 (std = 1.38), thus accounting for the likely scenario that some species will have more viruses and others less. We then extrapolated the average to all 1,200 bat species to estimate a total potential richness of 3,204 CoVs (range = 1,200–6,000 CoVs), most of which have yet to be described.</p>
<fig id="vex012-F6" orientation="portrait" position="float"><label>Figure 6.</label>
<caption><p> Relationship between sampling effort and viral detection among all bat species sampled. A Poisson regression model was used to estimate the expected number of viruses based on the number of animals sampled, by species (each circle indicates a separate species). The 95% confidence intervals are indicated.</p>
</caption>
<graphic xlink:href="vex012f6"></graphic>
</fig>
</sec>
<sec><title>3.4 Factors predicting CoV positivity</title>
<p>In order to refine future surveillance efforts to find the undiscovered diversity of CoVs in bats, we evaluated the factors predicting CoV positivity. For each region, the best model included specimen type, bat family (or in the case of Latin America, subfamily), season, and animal–human interface (<xref rid="vex012-T2" ref-type="table">Table 2</xref>
). Season was not included in the top model for Latin America, but the model with the season included was within two AIC (delta AIC = 1.06) and therefore considered to have as much support as the top model (<xref rid="vex012-B6" ref-type="bibr">Burnham and Anderson 2002</xref>
). Specimen type was highly associated with CoV positivity, and samples containing feces or fecal swabs were significantly more likely to test positive than other specimen types in all regions. Bat family was important in Asia and Africa. Season was also significant, with samples collected during the dry season more likely to test positive in Africa and Asia than those collected during the wet season (albeit with low odds ratios). Age class appears to be important, as sub-adults were more likely to test positive than adults in Africa and Asia. Sex was not related to CoV positivity in Asia or Latin America while males were slightly more likely to be CoV positive in Africa (OR = 1.5, 1.2–1.9 CI, <italic>P</italic>
< 0.001). Broad categories of animal–human interfaces were significantly associated with CoV positivity among bats in all three regions: in Africa and Latin America, bats sampled at the animal use interface category were more likely to be positive, and in Asia, bats sampled at the animal use and human activities interfaces were more likely to be positive than those sampled at other interfaces. The significance of the animal interface in all three regions suggests that practices around animal use may be important for disease transmission.
<table-wrap id="vex012-T2" orientation="portrait" position="float"><label>Table 2.</label>
<caption><p> Variables associated with coronavirus positive bat tests in Africa, Asia, and Latin America. Regional datasets for each logistic regression model consisted of bat species where greater than fifty bats were tested. Age datasets were a subset of the regional datasets as age class was not determined for all individuals. Numbers in bold are statically significant for <italic>P</italic>
< 0.05.</p>
</caption>
<table frame="hsides" rules="groups"><colgroup span="1"><col valign="top" align="left" span="1"></col>
<col valign="top" align="left" span="1"></col>
<col valign="top" align="left" span="1"></col>
<col valign="top" align="char" char="." span="1"></col>
<col valign="top" align="char" char="." span="1"></col>
<col valign="top" align="char" char="." span="1"></col>
<col valign="top" align="char" char="." span="1"></col>
</colgroup>
<thead><tr><th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1"></th>
<th rowspan="1" colspan="1">Odds ratio</th>
<th rowspan="1" colspan="1"><italic>P</italic>
</th>
<th colspan="2" rowspan="1">95% Confidence intervals</th>
</tr>
</thead>
<tbody><tr><td colspan="7" rowspan="1"><bold>Africa</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="2" rowspan="1">Sample type</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Lower</td>
<td rowspan="1" colspan="1">Upper</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Feces/rectal swab</td>
<td rowspan="1" colspan="1"><bold>350.98</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>146.64</bold>
</td>
<td rowspan="1" colspan="1"><bold>840.07</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Oral/nasal swab</td>
<td rowspan="1" colspan="1"><bold>5.05</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>2.03</bold>
</td>
<td rowspan="1" colspan="1"><bold>12.54</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Oral/rectal swab</td>
<td rowspan="1" colspan="1"><bold>30.98</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>13.46</bold>
</td>
<td rowspan="1" colspan="1"><bold>71.30</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Blood</td>
<td rowspan="1" colspan="1">1.65</td>
<td rowspan="1" colspan="1">0.540</td>
<td rowspan="1" colspan="1">0.33</td>
<td rowspan="1" colspan="1">8.09</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Tissue</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Host family</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Hipposideridae</td>
<td rowspan="1" colspan="1">1.22</td>
<td rowspan="1" colspan="1">0.559</td>
<td rowspan="1" colspan="1">0.62</td>
<td rowspan="1" colspan="1">2.41</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Pteropodidae</td>
<td rowspan="1" colspan="1">3.84</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>2.27</bold>
</td>
<td rowspan="1" colspan="1"><bold>6.49</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Molossidae</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Season</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Dry</td>
<td rowspan="1" colspan="1">2.42</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>1.79</bold>
</td>
<td rowspan="1" colspan="1"><bold>3.25</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Wet</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Interface</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Animal use</td>
<td rowspan="1" colspan="1">1.98</td>
<td rowspan="1" colspan="1"><bold>0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>1.35</bold>
</td>
<td rowspan="1" colspan="1"><bold>2.91</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Pristine area</td>
<td rowspan="1" colspan="1">1.39</td>
<td rowspan="1" colspan="1">0.624</td>
<td rowspan="1" colspan="1">0.37</td>
<td rowspan="1" colspan="1">5.29</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Land use</td>
<td rowspan="1" colspan="1">1.66</td>
<td rowspan="1" colspan="1">0.219</td>
<td rowspan="1" colspan="1">0.74</td>
<td rowspan="1" colspan="1">3.74</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Human activity</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Age</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Subadult</td>
<td rowspan="1" colspan="1"><bold>5.91</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>4.28</bold>
</td>
<td rowspan="1" colspan="1"><bold>8.17</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Adult</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td colspan="7" rowspan="1"><bold>Asia</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Sample type</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Feces/rectal swab</td>
<td rowspan="1" colspan="1"><bold>62.80</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>15.44</bold>
</td>
<td rowspan="1" colspan="1"><bold>255.49</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Oral/nasal swab</td>
<td rowspan="1" colspan="1"><bold>13.25</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>3.11</bold>
</td>
<td rowspan="1" colspan="1"><bold>56.38</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Urine/urogenital swab</td>
<td rowspan="1" colspan="1"><bold>46.92</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>10.85</bold>
</td>
<td rowspan="1" colspan="1"><bold>202.80</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Guano</td>
<td rowspan="1" colspan="1"><bold>22.12</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>3.66</bold>
</td>
<td rowspan="1" colspan="1"><bold>133.56</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Tissue, Oral/rectal swab</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Host family</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Emballonuridae</td>
<td rowspan="1" colspan="1">0.31</td>
<td rowspan="1" colspan="1">0.26</td>
<td rowspan="1" colspan="1">0.04</td>
<td rowspan="1" colspan="1">2.37</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Miniopteridae</td>
<td rowspan="1" colspan="1"><bold>9.78</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>5.69</bold>
</td>
<td rowspan="1" colspan="1"><bold>16.81</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Pteropodidae</td>
<td rowspan="1" colspan="1"><bold>2.16</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>1.29</bold>
</td>
<td rowspan="1" colspan="1"><bold>3.62</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Rhinolophidae</td>
<td rowspan="1" colspan="1">1.08</td>
<td rowspan="1" colspan="1">0.82</td>
<td rowspan="1" colspan="1">0.57</td>
<td rowspan="1" colspan="1">2.03</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Vespertilionidae</td>
<td rowspan="1" colspan="1"><bold>3.67</bold>
</td>
<td rowspan="1" colspan="1"><bold><0.001</bold>
</td>
<td rowspan="1" colspan="1"><bold>2.18</bold>
</td>
<td rowspan="1" colspan="1"><bold>6.18</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Hipposideridae</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Season</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Dry</td>
<td rowspan="1" colspan="1"><bold>1.49</bold>
</td>
<td rowspan="1" colspan="1"><bold>0.03</bold>
</td>
<td rowspan="1" colspan="1"><bold>1.03</bold>
</td>
<td rowspan="1" colspan="1"><bold>2.16</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Wet</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Interface</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Animal use</td>
<td rowspan="1" colspan="1"><bold>3.30</bold>
</td>
<td rowspan="1" colspan="1"><bold>0.01</bold>
</td>
<td rowspan="1" colspan="1"><bold>1.29</bold>
</td>
<td rowspan="1" colspan="1"><bold>8.40</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Human activity</td>
<td rowspan="1" colspan="1"><bold>3.48</bold>
</td>
<td rowspan="1" colspan="1"><bold>0.01</bold>
</td>
<td rowspan="1" colspan="1"><bold>1.36</bold>
</td>
<td rowspan="1" colspan="1"><bold>8.95</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Pristine area</td>
<td rowspan="1" colspan="1">1.81</td>
<td rowspan="1" colspan="1">0.35</td>
<td rowspan="1" colspan="1">0.52</td>
<td rowspan="1" colspan="1">6.35</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Land use</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Age</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Subadult</td>
<td rowspan="1" colspan="1"><bold>1.84</bold>
</td>
<td rowspan="1" colspan="1"><bold>0.00</bold>
</td>
<td rowspan="1" colspan="1"><bold>1.26</bold>
</td>
<td rowspan="1" colspan="1"><bold>2.67</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Adult</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td colspan="7" rowspan="1"><bold>Latin America</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Sample type</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Feces/rectal swab</td>
<td rowspan="1" colspan="1"><bold>15.66</bold>
</td>
<td rowspan="1" colspan="1"><bold>0.000</bold>
</td>
<td rowspan="1" colspan="1"><bold>5.02</bold>
</td>
<td rowspan="1" colspan="1"><bold>48.79</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Oral/rectal swab</td>
<td rowspan="1" colspan="1"><bold>27.86</bold>
</td>
<td rowspan="1" colspan="1"><bold>0.000</bold>
</td>
<td rowspan="1" colspan="1"><bold>6.86</bold>
</td>
<td rowspan="1" colspan="1"><bold>113.22</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Oral/nasal swab</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Host subfamily</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Carolliinae</td>
<td rowspan="1" colspan="1">6.95</td>
<td rowspan="1" colspan="1">0.06</td>
<td rowspan="1" colspan="1">0.90</td>
<td rowspan="1" colspan="1">53.76</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Stenodermatinae</td>
<td rowspan="1" colspan="1">5.04</td>
<td rowspan="1" colspan="1">0.12</td>
<td rowspan="1" colspan="1">0.67</td>
<td rowspan="1" colspan="1">37.89</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Glossophaginae</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Interface</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Animal use</td>
<td rowspan="1" colspan="1"><bold>5.25</bold>
</td>
<td rowspan="1" colspan="1"><bold>0.01</bold>
</td>
<td rowspan="1" colspan="1"><bold>1.48</bold>
</td>
<td rowspan="1" colspan="1"><bold>18.65</bold>
</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Human activity</td>
<td rowspan="1" colspan="1">2.73</td>
<td rowspan="1" colspan="1">0.12</td>
<td rowspan="1" colspan="1">0.77</td>
<td rowspan="1" colspan="1">9.69</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Land use</td>
<td rowspan="1" colspan="1">4.22</td>
<td rowspan="1" colspan="1">0.10</td>
<td rowspan="1" colspan="1">0.78</td>
<td rowspan="1" colspan="1">22.95</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Pristine area</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Season</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Dry</td>
<td rowspan="1" colspan="1">1.30</td>
<td rowspan="1" colspan="1">0.52</td>
<td rowspan="1" colspan="1">0.58</td>
<td rowspan="1" colspan="1">2.89</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Wet</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td colspan="6" rowspan="1">Age</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Subadult</td>
<td rowspan="1" colspan="1">1.73</td>
<td rowspan="1" colspan="1">0.15</td>
<td rowspan="1" colspan="1">0.83</td>
<td rowspan="1" colspan="1">3.62</td>
</tr>
<tr><td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1">Adult</td>
<td rowspan="1" colspan="1">1.00</td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
<td rowspan="1" colspan="1"></td>
</tr>
</tbody>
</table>
</table-wrap>
</p>
</sec>
<sec><title>3.5 Future sampling effort</title>
<p>Based on our study, we were able to estimate the sampling effort that would be required to find all CoVs in bats, based on a Poisson regression model of our data (<xref ref-type="fig" rid="vex012-F6">Fig. 6</xref>
). With 154 individuals, we estimate that an average of one CoV will be detected (95% CI 136–177). With 397 individuals, we estimate that up to five CoVs will be detected (95% CI 351–458). As we did not observe any species with greater than five viral sequence clusters, we make the assumption that sampling 397 individuals should capture the full diversity of CoVs in each bat species.</p>
</sec>
<sec><title>3.6 Predicting the distribution of unknown CoVs</title>
<p>Viral sub-clade and host family were not independent of each other (χ<sup>2</sup>
by clade = <italic>P</italic>
< 0.001; <xref ref-type="supplementary-material" rid="sup1">Supplementary Table S1</xref>
), demonstrating that different clades have significant associations with specific bat families. For example, subgroup 2d CoVs were significantly associated with pteropid bats, and were only identified in regions where pteropid bats exist. Likewise, subgroup 2b CoVs were associated with rhinolophid and hipposiderid bats and 2c with vespertilionid bats. To ‘predict’ the potential distribution of unknown CoVs, we plotted the known distribution of bats belonging to these families and assume that ‘hotspots’ of bat diversity infer hotspots of viral diversity for each sub-clade (<xref ref-type="fig" rid="vex012-F7">Fig. 7</xref>
). We note that the alphaCoV genus has been considered as one group, since they do not resolve well into sub-clades (<xref rid="vex012-B13" ref-type="bibr">de Groot et al. 2009</xref>
), and that no 2a map was generated given that only one bat virus from this clade was identified.
</p>
<fig id="vex012-F7" orientation="portrait" position="float"><label>Figure 7.</label>
<caption><p> Viral diversity ‘hotspot’ maps. Panel A shows the potential hotspots for CoVs based on the distribution of bats worldwide. Location of alpha CoV sequences from this study are shown in black and beta CoV sequences in blue, indicating there is no geographical bias based on viral genus (i.e. alpha and beta CoVs are equally likely to be found in all regions). Some viral sub-clades were associated with particular bat families, and the spatial distribution data of all species belonging to these families were plotted to indicate the potential hotspots of viral diversity (richness) for these sub-clades. Panel B indicates the potential distribution of 2b CoVs based on the distribution of rhinolphus and hipposideros bats. Locations of 2b-positive animals identified in this study are indicated in black, and correlate with areas of high species richness (for these families). CoV-positive animals for other sub-clades shown in light blue. Panel C indicates the potential distribution of 2c CoVs based on the distribution of vespertilionid bats. Locations of 2c-positive animals identified in this study are indicated in black. This map suggests there are hotspots of 2c diversity in regions not covered in this study (e.g. Europe). Panel D indicates the potential distribution of 2d viruses based on the distribution of pteropid bats. The map suggests these viruses may have a more limited distribution, compared with viruses of other sub-clades. Locations of 2d-positive animals identified in this study are shown in black.</p>
</caption>
<graphic xlink:href="vex012f7"></graphic>
</fig>
</sec>
</sec>
<sec><title>4. Discussion</title>
<p>The emergence of SARS and MERS has driven a need to understand more about the diversity, ecology and evolution of coronaviruses, particularly at so-called ‘hotspots’ of zoonotic emergence (<xref rid="vex012-B28" ref-type="bibr">Jones et al. 2008</xref>
; <xref rid="vex012-B15" ref-type="bibr">Drexler et al. 2014</xref>
). To this end, we surveyed the diversity of CoVs from twenty countries in Latin America, Africa, and Asia to identify global factors driving viral diversity and to look for regional differences in factors that contribute to the risk of emergence, such as host switching. In total, we identified sequences from 100 discrete phylogenetic clusters, ninety-one of which were found in bats.</p>
<p>Our data suggest that the diversity of bat CoVs has been driven primarily by host ecology. First, viral richness was strongly correlated with bat richness, suggesting that most CoVs will be found in regions where bat diversity is highest. Second, we showed that CoV diversity separates into three distinct communities by region, echoing the distribution of bats and suggesting an ecological dependence on their hosts. And third, we identified particular associations between viral sub-clade and bat family, indicating that CoVs have evolved with (or adapted to) preferred families. Collectively, these data show that the global diversity and distribution of CoVs in bats is non-random and is driven by variation in the biogeography of bats.</p>
<p>Regional variation was also observed in the proportion of host switching events, relative to other evolutionary mechanisms such as co-speciation, duplication, or sharing (see Methods for definitions). Overall, host switching was the dominant mechanism, supporting the general trend observed for CoVs (<xref rid="vex012-B59" ref-type="bibr">Vijaykrishna et al. 2007</xref>
; <xref rid="vex012-B67" ref-type="bibr">Woo et al. 2009</xref>
; <xref rid="vex012-B36" ref-type="bibr">Lau et al. 2012</xref>
) as well as similar findings for other viral families, including paramxyoviruses (<xref rid="vex012-B43" ref-type="bibr">Melade et al. 2016</xref>
), hantaviruses (<xref rid="vex012-B50" ref-type="bibr">Ramsden et al. 2009</xref>
), and arenaviruses (<xref rid="vex012-B11" ref-type="bibr">Coulibaly-N'Golo et al. 2011</xref>
; <xref rid="vex012-B26" ref-type="bibr">Irwin et al. 2012</xref>
). However, when analyzed by region, there were proportionally fewer events in Latin America compared with Africa or Asia. Given that host switching is the first critical step in zoonotic emergence (<xref rid="vex012-B39" ref-type="bibr">Li et al. 2005</xref>
;<xref rid="vex012-B46" ref-type="bibr">Pfefferle et al. 2009</xref>
; <xref rid="vex012-B68" ref-type="bibr">Woo et al. 2012</xref>
; <xref rid="vex012-B4" ref-type="bibr">Azhar et al. 2014</xref>
), we suggest this finding could reflect regional differences in the risk of disease emergence.</p>
<p>It is important to note that we distinguish ‘host switching’ (inter-genus or family switching) from ‘sharing’ (intra-genus switching) in our analysis and that while the proportion of host switches was lower in Latin America, the proportion of sharing events was concomitantly higher. This indicates that CoVs in this region still switch hosts, just like they do in Africa and Asia, but that they preferentially move between closely related species. While we do not yet understand the factors driving this difference, if we assume that the risk of spillover to humans reflects the taxonomic distances that CoVs are able to jump (<xref rid="vex012-B33" ref-type="bibr">Kreuder Johnson et al. 2015</xref>
), our results would indicate a higher risk for zoonotic emergence in Africa and Asia [acknowledging that host switching is not the only important factor driving zoonotic emergence (<xref rid="vex012-B28" ref-type="bibr">Jones et al. 2008</xref>
; <xref rid="vex012-B31" ref-type="bibr">Keesing et al. 2010</xref>
; <xref rid="vex012-B30" ref-type="bibr">Karesh et al. 2012</xref>
)]. While evidence of regional variation in host switching has not been reported previously for viruses (to our knowledge), similar patterns have been observed in hemosporidian parasites (<xref rid="vex012-B18" ref-type="bibr">Ellis et al. 2015</xref>
).</p>
<p>In total, we estimated that there are at least 3,204 CoVs (based on the cluster definition used in this study) in bats. This suggests that although there is still a large number of coronaviruses that remain to be discovered, their detection remains an achievable goal, and we advocate for their discovery given the potential for even greater insights into the global ecology and evolution of these viruses (e.g. further investigation into regional differences in host switching), and the opportunity to elucidate their individual zoonotic potential (<xref rid="vex012-B20" ref-type="bibr">Ge et al. 2013</xref>
; <xref rid="vex012-B62" ref-type="bibr">Wang et al. 2014</xref>
; <xref rid="vex012-B69" ref-type="bibr">Yang et al. 2014</xref>
; <xref rid="vex012-B3" ref-type="bibr">Anthony et al. 2017</xref>
). Building on our study, we suggest that future surveillance efforts should consider the number of samples carefully and include up to 400 individuals (of either sex) per species in order to maximize the chance of detecting all CoVs, and that including less than 154 individuals per species might reduce the chance of finding even one virus and would likely produce a poor return on investment. We also recommend that feces (or rectal swabs) followed by saliva (or oral swabs) should be prioritized if resources are limited and all sample types cannot be processed, and that sampling should be biased towards immature animals and in the dry season.</p>
<p>To predict where this unknown diversity is likely to be, we plotted the known distribution of all bat families associated with each CoV sub-clade. The 2b SARS clade was associated with rhinolophid and hipposiderid bats, so we plotted the distribution of all bat species belonging to these families to generate a map of viral diversity (richness). Based on this map, we would expect the greatest diversity of unknown 2b viruses to be found in South East Asia, which is consistent with previous studies describing 2b viruses in bats, but also in isolated pockets throughout Africa. In contrast, the 2c MERS clade was associated with vespertilionid bats, predicting that 2c diversity will be highest in Mexico, Europe, Central and South East Africa, and parts of South East Asia. Again, this is consistent with the distribution of 2c viruses reported here, and previously in the literature (<xref rid="vex012-B66" ref-type="bibr">Woo et al. 2007</xref>
; <xref rid="vex012-B51" ref-type="bibr">Reusken et al. 2010</xref>
; <xref rid="vex012-B1" ref-type="bibr">Anthony et al. 2012</xref>
; <xref rid="vex012-B38" ref-type="bibr">Lelli et al. 2013</xref>
; <xref rid="vex012-B9" ref-type="bibr">Corman et al. 2014a</xref>
,<xref rid="vex012-B10" ref-type="bibr">b</xref>
; <xref rid="vex012-B3" ref-type="bibr">Anthony et al. 2017</xref>
).</p>
<p>Certain limitations should be considered in the interpretation of our data. First, the diversity of CoV sequences detected is almost certainly biased by our sampling effort. A small number of species were sampled abundantly, thus maximizing our chances of finding a CoV (27/282 species had >110 individuals); however, most were under-sampled (232/282 with ≤50 individuals). While this biases the probability of finding a CoV towards the more abundantly sampled species, we stress that it also reflects the uneven distribution of naturally occurring communities (<xref rid="vex012-B42" ref-type="bibr">Magurran 2011</xref>
), and that targeting the very rare species would be prohibitively expensive. It is also unclear if population size has an effect on the number of viruses present in a species—e.g. do larger populations accommodate more viral diversity? Second our analysis only represents the ‘last’ evolutionary event identified. It does not represent, or account for, the complete evolutionary history of these lineages. Therefore, while a single event has been ascribed to each association (see ‘Methods’), it does not preclude other events having an equal or greater impact in the past. We are limited to using this most recent event type in order to assign a single evolutionary event to each unique association. It should also be noted that events are dependent on the relationships observed in the reconstructed trees and could change as more viruses are added or longer sequences are used. Third, we have only generated partial sequences for analysis, which limits our ability to comment on the specific zoonotic potential of each virus or provide a comprehensive analysis of their genetic histories and taxonomy, for example, estimating the frequency or impact of recombination, which can be an important driver of host switching and disease emergence (<xref rid="vex012-B3" ref-type="bibr">Anthony et al. 2017</xref>
). We certainly support full-genome sequencing wherever possible [and have such efforts underway (<xref rid="vex012-B3" ref-type="bibr">Anthony et al. 2017</xref>
)], however difficulties in virus isolation or sequencing directly from swabs (where material and viral load can be low, yielding variable results) can often preclude extending sequences much beyond the short fragments amplified by consensus PCR (<xref rid="vex012-B15" ref-type="bibr">Drexler et al. 2014</xref>
). Equally, logistic and permitting issues in moving biological samples across borders, and the need to first develop in-country capacity to fully sequence these viruses, have all limited our ability to characterize (most of) these viruses further.</p>
<p>Our ‘state of the art’ is in the unique scope of our study, focusing on a single viral family at a global scale. To achieve this, the USAID PREDICT project first had to establish a strategy for viral discovery in twenty different countries, many of which had no prior capacity to do this work at all. For this reason, we used a simple yet highly implementable approach based on cPCR. Initially, this approach might not appear to be as powerful as high throughput sequencing (HTS) techniques; however, this study has demonstrated the particular and considerable strengths of cPCR as an affordable screening and discovery tool in resource-poor settings, where HTS is still largely unsupported.</p>
<p>Studying viral diversity on a global scale is just one component of a larger strategy to understand the factors driving viral diversity. Ecological and evolutionary mechanisms can contribute differently across scales, and we therefore advocate complimenting the global perspective with studies that focus on entire viral communities (rather than only one viral family) within single individuals or host species (rather than only on a global scale) (<xref rid="vex012-B2" ref-type="bibr">Anthony et al. 2015</xref>
). We propose that it is critical to understand both the component parts of the ‘zoonotic pool’ (<xref rid="vex012-B45" ref-type="bibr">Morse et al. 2012</xref>
) and the nature of their interactions at different scales if we are to move towards better predictive models, and ultimately to reducing zoonotic emergence.</p>
<p>Finally, we offer a specific comment regarding the public health threat posed by bats. While it is tempting to conclude that bats harbor a large number of potentially zoonotic CoVs, most of the putative viruses detected in this study are unlikely to pose any threat to humans—either because they lack the biological pre-requisites to infect human cells or because the ecology of their hosts limits the opportunity for spillover. Studies such as this are intended to advance our understanding of the fundamental biology of viruses, not to create alarm or incite the retaliatory culling of bats. Indeed, such actions often have unanticipated consequences; and can even enhance disease transmission, as seen with rabies in vampire bats (<xref rid="vex012-B55" ref-type="bibr">Streicker et al. 2012</xref>
; <xref rid="vex012-B5" ref-type="bibr">Blackwood et al. 2013</xref>
). Bats are important insectivores, pollinators and seed dispersers and we underscore both their vital ecosystem role and the need to consider any public health interventions carefully.</p>
</sec>
<sec sec-type="supplementary-material"><title>Supplementary data</title>
<p><xref ref-type="supplementary-material" rid="sup1">Supplementary data</xref>
are available at <italic>Virus Evolution</italic>
online.</p>
<p>Conflict of interest: None declared.</p>
</sec>
<sec sec-type="supplementary-material"><title>Supplementary Material</title>
<supplementary-material content-type="local-data" id="sup1"><label>Supplementary Data</label>
<media xlink:href="vex012_Supp.docx"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><title>Acknowledgements</title>
<p>This study was made possible by the generous support of the American people through the United States Agency for International Development (USAID) Emerging Pandemic Threats PREDICT project (cooperative agreement number GHN-A-OO-09-00010-00). We thank the governments of Cameroon, Gabon, Democratic Republic of Congo, Republic of Congo, Rwanda, Tanzania, Uganda, Peru, Bolivia, Brazil, Mexico, Bangladesh, Cambodia, China, Indonesia, Laos, Malaysia, Nepal, Thailand, and Viet Nam for permission to conduct this study, and the field teams and collaborating laboratories that performed sample collection and testing.</p>
</ack>
<ref-list><title>References</title>
<ref id="vex012-B1"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Anthony</surname>
<given-names>S. J.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <article-title>‘Coronaviruses in bats from Mexico’</article-title>
, <source>Journal of General Virology</source>
, <volume>94</volume>
: <fpage>1028</fpage>
–<lpage>38</lpage>
.</mixed-citation>
</ref>
<ref id="vex012-B2"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Anthony</surname>
<given-names>S. J.</given-names>
</name>
</person-group>
, <etal></etal>
(<year>2015</year>
) <article-title>‘Non-random patterns in viral diversity’</article-title>
, <source>Nature Communications</source>
, <volume>6</volume>
: <fpage>8147</fpage>
.</mixed-citation>
</ref>
<ref id="vex012-B3"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Anthony</surname>
<given-names>S. J.</given-names>
</name>
</person-group>
, <etal></etal>
(<year>2017</year>
) <article-title>‘Further evidence for bats as the evolutionary source of MERS coronavirus’</article-title>
, <source>mBio</source>
, <volume>8</volume>
/<issue>2</issue>
. pii: e00373-17. doi: 10.1128/mBio.00373-17.</mixed-citation>
</ref>
<ref id="vex012-B4"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Azhar</surname>
<given-names>E. I.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2014</year>
) <article-title>‘Evidence for camel-to-human transmission of MERS coronavirus’</article-title>
, <source>The New England Journal of Medicine</source>
, <volume>370</volume>
: <fpage>2499</fpage>
–<lpage>505</lpage>
.<pub-id pub-id-type="pmid">24896817</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B5"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Blackwood</surname>
<given-names>J. C.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2013</year>
) <article-title>‘Resolving the roles of immunity, pathogenesis, and immigration for rabies persistence in vampire bats’</article-title>
, <source>Proceedings of the National Academy of Sciences of the United States of America</source>
, <volume>110</volume>
: <fpage>20837</fpage>
–<lpage>42</lpage>
.<pub-id pub-id-type="pmid">24297874</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B6"><mixed-citation publication-type="book"><person-group person-group-type="author"><name name-style="western"><surname>Burnham</surname>
<given-names>K. P.</given-names>
</name>
, <name name-style="western"><surname>Anderson</surname>
<given-names>D. R.</given-names>
</name>
</person-group>
(<year>2002</year>
) <source>Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach</source>
. <publisher-loc>New York</publisher-loc>
: <publisher-name>Springer</publisher-name>
.</mixed-citation>
</ref>
<ref id="vex012-B7"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Charleston</surname>
<given-names>M. A.</given-names>
</name>
</person-group>
(<year>1998</year>
) <article-title>‘Jungles: a new solution to the host/parasite phylogeny reconciliation problem’</article-title>
, <source>Mathematical Biosciences</source>
, <volume>149</volume>
: <fpage>191</fpage>
–<lpage>223</lpage>
.<pub-id pub-id-type="pmid">9621683</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B8"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Conow</surname>
<given-names>C.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2010</year>
) <article-title>‘Jane: a new tool for the cophylogeny reconstruction problem’</article-title>
, <source>Algorithms for Molecular Biology: AMB</source>
, <volume>5</volume>
: <fpage>16</fpage>
.<pub-id pub-id-type="pmid">20181081</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B9"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Corman</surname>
<given-names>V. M.</given-names>
</name>
</person-group>
<etal></etal>
(<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="vex012-B10"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Corman</surname>
<given-names>V. M.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2014b</year>
) <article-title>‘Characterization of a novel betacoronavirus related to middle East respiratory syndrome coronavirus in European hedgehogs’</article-title>
, <source>Journal of Virology</source>
, <volume>88</volume>
: <fpage>717</fpage>
–<lpage>24</lpage>
.<pub-id pub-id-type="pmid">24131722</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B11"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Coulibaly-N'Golo</surname>
<given-names>D.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2011</year>
) <article-title>‘Novel arenavirus sequences in Hylomyscus sp. and Mus (Nannomys) setulosus from Cote d'Ivoire: implications for evolution of arenaviruses in Africa’</article-title>
, <source>PLoS One</source>
, <volume>6</volume>
: <fpage>e20893</fpage>
.<pub-id pub-id-type="pmid">21695269</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B12"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Darriba</surname>
<given-names>D.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <article-title>‘jModelTest 2: more models, new heuristics and parallel computing’</article-title>
, <source>Nature Methods</source>
, <volume>9</volume>
: <fpage>772.</fpage>
</mixed-citation>
</ref>
<ref id="vex012-B13"><mixed-citation publication-type="book"><person-group person-group-type="author"><name name-style="western"><surname>de Groot</surname>
<given-names>R. J.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2009</year>
) <chapter-title>‘Virus taxonomy: classification and nomenclature of viruses’</chapter-title>
, in <person-group person-group-type="editor"><name name-style="western"><surname>King</surname>
<given-names>A. M. Q.</given-names>
</name>
, <name name-style="western"><surname>Adams</surname>
<given-names>M. J.</given-names>
</name>
, <name name-style="western"><surname>Carstens</surname>
<given-names>E. B.</given-names>
</name>
, <name name-style="western"><surname>Lefkkowitch</surname>
<given-names>J.</given-names>
</name>
</person-group>
(eds) <source>Ninth Report of the International Committee for the Taxonomy of Viruses</source>
. <publisher-loc>Amsterdam</publisher-loc>
: <publisher-name>Elsevier</publisher-name>
.</mixed-citation>
</ref>
<ref id="vex012-B14"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Donaldson</surname>
<given-names>E. F.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2010</year>
) <article-title>‘Metagenomic analysis of the viromes of three North American bat species: viral diversity among different bat species that share a common habitat’</article-title>
, <source>Journal of Virology</source>
, <volume>84</volume>
: <fpage>13004</fpage>
–<lpage>18</lpage>
.<pub-id pub-id-type="pmid">20926577</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B15"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Drexler</surname>
<given-names>J. F.</given-names>
</name>
, <name name-style="western"><surname>Corman</surname>
<given-names>V. M.</given-names>
</name>
, <name name-style="western"><surname>Drosten</surname>
<given-names>C.</given-names>
</name>
</person-group>
(<year>2014</year>
) <article-title>‘Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS’</article-title>
, <source>Antiviral Research</source>
, <volume>101</volume>
: <fpage>45</fpage>
–<lpage>56</lpage>
<pub-id pub-id-type="pmid">24184128</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B16"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Drosten</surname>
<given-names>C.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2003</year>
) <article-title>‘Identification of a novel coronavirus in patients with severe acute respiratory syndrome’</article-title>
, <source>The New England Journal of Medicine</source>
, <volume>348</volume>
: <fpage>1967</fpage>
–<lpage>76</lpage>
.<pub-id pub-id-type="pmid">12690091</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B17"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Edgar</surname>
<given-names>R. C.</given-names>
</name>
</person-group>
(<year>2004</year>
) <article-title>‘MUSCLE: multiple sequence alignment with high accuracy and high throughput’</article-title>
, <source>Nucleic Acids Research</source>
, <volume>32</volume>
: <fpage>1792</fpage>
–<lpage>7</lpage>
.<pub-id pub-id-type="pmid">15034147</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B18"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Ellis</surname>
<given-names>V. A.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2015</year>
) <article-title>‘Local host specialization, host-switching, and dispersal shape the regional distributions of avian haemosporidian parasites’</article-title>
, <source>Proceedings of the National Academy of Sciences of the United States of America</source>
, <volume>112</volume>
: <fpage>11294</fpage>
–<lpage>9</lpage>
.<pub-id pub-id-type="pmid">26305975</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B19"><mixed-citation publication-type="other"><collab>ESRI</collab>
. (<year>2011</year>
) <italic>ArcGIS Desktop: Release 10</italic>
. Redlands, CA: Environmental Systems Research Institute.</mixed-citation>
</ref>
<ref id="vex012-B20"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Ge</surname>
<given-names>X. Y.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2013</year>
) <article-title>‘Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor’</article-title>
, <source>Nature</source>
, <volume>503</volume>
: <fpage>535</fpage>
–<lpage>8</lpage>
.<pub-id pub-id-type="pmid">24172901</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B21"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Hijmans</surname>
<given-names>R. J.</given-names>
</name>
, <name name-style="western"><surname>Williams</surname>
<given-names>E.</given-names>
</name>
, <name name-style="western"><surname>Vennes</surname>
<given-names>C.</given-names>
</name>
</person-group>
(<year>2015</year>
) <comment>Geosphere: Spherical Trigonometry v. R Package version 1.5</comment>
.</mixed-citation>
</ref>
<ref id="vex012-B22"><mixed-citation publication-type="book"><person-group person-group-type="author"><name name-style="western"><surname>Hagberg</surname>
<given-names>A. A.</given-names>
</name>
, <name name-style="western"><surname>Schult</surname>
<given-names>D. A.</given-names>
</name>
, <name name-style="western"><surname>Swart</surname>
<given-names>P. J.</given-names>
</name>
</person-group>
(<year>2008</year>
) <chapter-title>‘Exploring network structure, dynamics, and function using NetworkX’</chapter-title>
, in <person-group person-group-type="editor"><name name-style="western"><surname>Varoquaux</surname>
<given-names>G.</given-names>
</name>
, <name name-style="western"><surname>Vaught</surname>
<given-names>T.</given-names>
</name>
, <name name-style="western"><surname>Millman</surname>
<given-names>J.</given-names>
</name>
</person-group>
(eds.) <source>Proceedings of 7th Python in Science Conference (SciPy2008)</source>
, pp. <fpage>11</fpage>
–<lpage>15</lpage>
. <publisher-loc>Pasadena, CA</publisher-loc>
.</mixed-citation>
</ref>
<ref id="vex012-B23"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Hu</surname>
<given-names>B.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2015</year>
) <article-title>‘Bat origin of human coronaviruses’</article-title>
, <source>Virology Journal</source>
, <volume>12</volume>
: <fpage>221.</fpage>
<pub-id pub-id-type="pmid">26689940</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B24"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Huang</surname>
<given-names>Y. W.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2013</year>
) <article-title>‘Origin, evolution, and genotyping of emergent porcine epidemic diarrhea virus strains in the United States’</article-title>
, <source>mBio</source>
, <volume>4</volume>
: <fpage>e00737</fpage>
–<lpage>13</lpage>
.<pub-id pub-id-type="pmid">24129257</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B25"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Huynh</surname>
<given-names>J.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <article-title>‘Evidence supporting a zoonotic origin of human coronavirus strain NL63’</article-title>
, <source>Journal of Virology</source>
, <volume>86</volume>
: <fpage>12816</fpage>
–<lpage>25</lpage>
.<pub-id pub-id-type="pmid">22993147</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B26"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Irwin</surname>
<given-names>N. R.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <article-title>‘Complex patterns of host switching in New World arenaviruses’</article-title>
, <source>Molecular Ecology</source>
, <volume>21</volume>
: <fpage>4137</fpage>
–<lpage>50</lpage>
.<pub-id pub-id-type="pmid">22693963</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B27"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Jacomy</surname>
<given-names>M.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2014</year>
) <article-title>‘ForceAtlas2, a continuous graph layout algorithm for handy network visualization designed for the Gephi software’</article-title>
, <source>PLoS One</source>
, <volume>9</volume>
: <fpage>e98679.</fpage>
<pub-id pub-id-type="pmid">24914678</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B28"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Jones</surname>
<given-names>K. E.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2008</year>
) <article-title>‘Global trends in emerging infectious diseases’</article-title>
, <source>Nature</source>
, <volume>451</volume>
: <fpage>990</fpage>
–<lpage>3</lpage>
.<pub-id pub-id-type="pmid">18288193</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B29"><mixed-citation publication-type="book"><person-group person-group-type="author"><name name-style="western"><surname>Jost</surname>
<given-names>L.</given-names>
</name>
, <name name-style="western"><surname>Chao</surname>
<given-names>A.</given-names>
</name>
, <name name-style="western"><surname>Chazdon</surname>
<given-names>R. L.</given-names>
</name>
</person-group>
(<year>2011</year>
) <chapter-title>‘Compositional similarity and beta diversity’</chapter-title>
in <person-group person-group-type="editor"><name name-style="western"><surname>Magurran</surname>
<given-names>A. E.</given-names>
</name>
, <name name-style="western"><surname>McGill</surname>
<given-names>B. J.</given-names>
</name>
</person-group>
(eds.) <source>Biological Diversity: Frontiers in Measurement and Assessment</source>
. <publisher-loc>Oxford</publisher-loc>
: <publisher-name>Oxford University Press</publisher-name>
.</mixed-citation>
</ref>
<ref id="vex012-B30"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Karesh</surname>
<given-names>W. B.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <article-title>‘Ecology of zoonoses: natural and unnatural histories’</article-title>
, <source>Lancet</source>
, <volume>380</volume>
: <fpage>1936</fpage>
–<lpage>45</lpage>
.<pub-id pub-id-type="pmid">23200502</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B31"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Keesing</surname>
<given-names>F.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2010</year>
) <article-title>‘Impacts of biodiversity on the emergence and transmission of infectious diseases’</article-title>
, <source>Nature</source>
, <volume>468</volume>
: <fpage>647</fpage>
–<lpage>52</lpage>
.<pub-id pub-id-type="pmid">21124449</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B32"><mixed-citation publication-type="book"><person-group person-group-type="author"><name name-style="western"><surname>King</surname>
<given-names>A. M. Q.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <source>Virus Taxonomy: Classification and Nomenclature of Viruses. Ninth Report of the International Committee on Taxonomy of Viruses</source>
. <publisher-name>Elsevier Academic Press</publisher-name>
.</mixed-citation>
</ref>
<ref id="vex012-B33"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Kreuder Johnson</surname>
<given-names>C.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2015</year>
) <article-title>‘Spillover and pandemic properties of zoonotic viruses with high host plasticity’</article-title>
, <source>Scientific Reports</source>
, <volume>5</volume>
: <fpage>14830</fpage>
.<pub-id pub-id-type="pmid">26445169</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B34"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Ksiazek</surname>
<given-names>T. G.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2003</year>
) <article-title>‘A novel coronavirus associated with severe acute respiratory syndrome’</article-title>
, <source>The New England Journal of Medicine</source>
, <volume>348</volume>
: <fpage>1953</fpage>
–<lpage>66</lpage>
.<pub-id pub-id-type="pmid">12690092</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B35"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Lau</surname>
<given-names>S. K.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2005</year>
) <article-title>‘Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats’</article-title>
, <source>Proceedings of the National Academy of Sciences of the United States of America</source>
, <volume>102</volume>
: <fpage>14040</fpage>
–<lpage>5</lpage>
.<pub-id pub-id-type="pmid">16169905</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B36"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Lau</surname>
<given-names>S. K.</given-names>
</name>
</person-group>
, <etal></etal>
(<year>2012</year>
) <article-title>‘Recent transmission of a novel alphacoronavirus, bat coronavirus HKU10, from Leschenault's rousettes to pomona leaf-nosed bats: first evidence of interspecies transmission of coronavirus between bats of different suborders’</article-title>
, <source>Journal of Virology</source>
, <volume>86</volume>
: <fpage>11906</fpage>
–<lpage>18</lpage>
.<pub-id pub-id-type="pmid">22933277</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B37"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Lau</surname>
<given-names>S. K.</given-names>
</name>
</person-group>
, <etal></etal>
(<year>2015</year>
) <article-title>‘Discovery of a novel coronavirus, China Rattus coronavirus HKU24, from Norway rats supports the murine origin of Betacoronavirus 1 and has implications for the ancestor of Betacoronavirus lineage A’</article-title>
, <source>Journal of Virology</source>
, <volume>89</volume>
: <fpage>3076</fpage>
–<lpage>92</lpage>
.<pub-id pub-id-type="pmid">25552712</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B38"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Lelli</surname>
<given-names>D.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2013</year>
) <article-title>‘Detection of coronaviruses in bats of various species in Italy’</article-title>
, <source>Viruses</source>
, <volume>5</volume>
: <fpage>2679</fpage>
–<lpage>89</lpage>
.<pub-id pub-id-type="pmid">24184965</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B39"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Li</surname>
<given-names>W.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2005</year>
) <article-title>‘Bats are natural reservoirs of SARS-like coronaviruses’</article-title>
, <source>Science</source>
, <volume>310</volume>
: <fpage>676</fpage>
–<lpage>9</lpage>
.<pub-id pub-id-type="pmid">16195424</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B40"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Liang</surname>
<given-names>K.</given-names>
</name>
, <name name-style="western"><surname>Zeger</surname>
<given-names>S. L.</given-names>
</name>
</person-group>
(<year>1986</year>
) <article-title>‘Longitudinal data analysis using generalized linear models’</article-title>
, <source>Biometrika</source>
, <volume>73</volume>
: <fpage>13</fpage>
–<lpage>22</lpage>
.</mixed-citation>
</ref>
<ref id="vex012-B41"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Maes</surname>
<given-names>P.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2009</year>
) <article-title>‘A proposal for new criteria for the classification of hantaviruses, based on S and M segment protein sequences’</article-title>
, <source>Infection, Genetics and Evolution: Journal of Molecular Epidemiology and Evolutionary Genetics in Infectious Diseases</source>
, <volume>9</volume>
: <fpage>813</fpage>
–<lpage>20</lpage>
.</mixed-citation>
</ref>
<ref id="vex012-B42"><mixed-citation publication-type="book"><person-group person-group-type="author"><name name-style="western"><surname>Magurran</surname>
<given-names>A. E.</given-names>
</name>
</person-group>
(<year>2011</year>
) <source>Biological Diversity: Frontiers in Measurement and Assessment</source>
. <publisher-loc>Oxford</publisher-loc>
: <publisher-name>Oxford University Press</publisher-name>
.</mixed-citation>
</ref>
<ref id="vex012-B43"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Melade</surname>
<given-names>J.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2016</year>
) <article-title>‘An eco-epidemiological study of Morbilli-related paramyxovirus infection in Madagascar bats reveals host-switching as the dominant macro-evolutionary mechanism’</article-title>
, <source>Scientific Reports</source>
, <volume>6</volume>
: <fpage>23752</fpage>
.<pub-id pub-id-type="pmid">27068130</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B44"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Merkle</surname>
<given-names>D.</given-names>
</name>
, <name name-style="western"><surname>Middendorf</surname>
<given-names>M.</given-names>
</name>
, <name name-style="western"><surname>Wieseke</surname>
<given-names>N.</given-names>
</name>
</person-group>
(<year>2010</year>
) <article-title>‘A parameter-adaptive dynamic programming approach for inferring cophylogenies’</article-title>
, <source>BMC Bioinformatics</source>
, <volume>11 Suppl 1</volume>
: <fpage>S60.</fpage>
<pub-id pub-id-type="pmid">20122236</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B45"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Morse</surname>
<given-names>S. S.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <article-title>‘Prediction and prevention of the next pandemic zoonosis’</article-title>
, <source>The Lancet</source>
, <volume>380</volume>
: <fpage>1956</fpage>
–<lpage>65</lpage>
.</mixed-citation>
</ref>
<ref id="vex012-B46"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Pfefferle</surname>
<given-names>S.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2009</year>
) <article-title>‘Distant relatives of severe acute respiratory syndrome coronavirus and close relatives of human coronavirus 229E in bats, Ghana’</article-title>
, <source>Emerging Infectious Diseases</source>
, <volume>15</volume>
: <fpage>1377</fpage>
–<lpage>84</lpage>
.<pub-id pub-id-type="pmid">19788804</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B47"><mixed-citation publication-type="book"><collab>PREDICT_Consortium</collab>
. (<year>2014</year>
) <source>Reducing Pandemic Risk</source>
. Promoting Global Health. <publisher-name>One Health Institute, University of California Davis</publisher-name>
, <publisher-loc>Davis, CA</publisher-loc>
.</mixed-citation>
</ref>
<ref id="vex012-B48"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Quan</surname>
<given-names>P. L.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2010</year>
) <article-title>‘Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria’</article-title>
, <source>mBio</source>
, <volume>1</volume>
: pii: e00208-10. doi:10.1128/mBio.00208-10.</mixed-citation>
</ref>
<ref id="vex012-B49"><mixed-citation publication-type="book"><collab>R Foundation for Statistical Computing</collab>
. (<year>2008</year>
). <source>R: A Language and Environment for Statistical Computing</source>
. <publisher-loc>Vienna, Austria</publisher-loc>
: <publisher-name>R Foundation for Statistical Computing</publisher-name>
.</mixed-citation>
</ref>
<ref id="vex012-B50"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Ramsden</surname>
<given-names>C.</given-names>
</name>
, <name name-style="western"><surname>Holmes</surname>
<given-names>E. C.</given-names>
</name>
, <name name-style="western"><surname>Charleston</surname>
<given-names>M. A.</given-names>
</name>
</person-group>
(<year>2009</year>
) <article-title>‘Hantavirus evolution in relation to its rodent and insectivore hosts: no evidence for codivergence’</article-title>
, <source>Molecular Biology and Evolution</source>
, <volume>26</volume>
: <fpage>143</fpage>
–<lpage>53</lpage>
.<pub-id pub-id-type="pmid">18922760</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B51"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Reusken</surname>
<given-names>C. B.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2010</year>
) <article-title>‘Circulation of group 2 coronaviruses in a bat species common to urban areas in Western Europe’</article-title>
, <source>Vector Borne Zoonotic Dis</source>
, <volume>10</volume>
: <fpage>785</fpage>
–<lpage>91</lpage>
.<pub-id pub-id-type="pmid">20055576</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B52"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Reusken</surname>
<given-names>C. B.</given-names>
</name>
</person-group>
, <etal></etal>
(<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="vex012-B53"><mixed-citation publication-type="book"><person-group person-group-type="author"><name name-style="western"><surname>Rubin</surname>
<given-names>D. B.</given-names>
</name>
</person-group>
(<year>1987</year>
) <source>Multiple Imputation for Nonresponse in Surveys</source>
. <publisher-loc>New York</publisher-loc>
: <publisher-name>John Wiley and Sons</publisher-name>
.</mixed-citation>
</ref>
<ref id="vex012-B54"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Sabir</surname>
<given-names>J. S.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2016</year>
) <article-title>‘Co-circulation of three camel coronavirus species and recombination of MERS-CoVs in Saudi Arabia’</article-title>
, <source>Science</source>
, <volume>351</volume>
: <fpage>81</fpage>
–<lpage>4</lpage>
.<pub-id pub-id-type="pmid">26678874</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B55"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Streicker</surname>
<given-names>D. G.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <article-title>‘Ecological and anthropogenic drivers of rabies exposure in vampire bats: implications for transmission and control’</article-title>
, <source>Proceedings. Biological Sciences/the Royal Society</source>
, <volume>279</volume>
: <fpage>3384</fpage>
–<lpage>92</lpage>
.</mixed-citation>
</ref>
<ref id="vex012-B56"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Tamura</surname>
<given-names>K.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2013</year>
) <article-title>‘MEGA6: molecular evolutionary genetics analysis version 6.0’</article-title>
, <source>Molecular Biology and Evolution</source>
, <volume>30</volume>
: <fpage>2725</fpage>
–<lpage>9</lpage>
.<pub-id pub-id-type="pmid">24132122</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B57"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Tang</surname>
<given-names>X. C.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2006</year>
) <article-title>‘Prevalence and genetic diversity of coronaviruses in bats from China’</article-title>
, <source>Journal of Virology</source>
, <volume>80</volume>
: <fpage>7481</fpage>
–<lpage>90</lpage>
.<pub-id pub-id-type="pmid">16840328</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B58"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Tsoleridis</surname>
<given-names>T.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2016</year>
) <article-title>‘Discovery of novel alphacoronaviruses in European rodents and shrews’</article-title>
, <source>Viruses</source>
, <volume>8</volume>
: <fpage>84</fpage>
.<pub-id pub-id-type="pmid">27102167</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B59"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Vijaykrishna</surname>
<given-names>D.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2007</year>
) <article-title>‘Evolutionary insights into the ecology of coronaviruses’</article-title>
, <source>Journal of Virology</source>
, <volume>81</volume>
: <fpage>4012</fpage>
–<lpage>20</lpage>
.<pub-id pub-id-type="pmid">17267506</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B60"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Vijgen</surname>
<given-names>L.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2005</year>
) <article-title>‘Complete genomic sequence of human coronavirus OC43: molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event’</article-title>
, <source>Journal of Virology</source>
, <volume>79</volume>
: <fpage>1595</fpage>
–<lpage>604</lpage>
.<pub-id pub-id-type="pmid">15650185</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B61"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Wacharapluesadee</surname>
<given-names>S.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2015</year>
) <article-title>‘Diversity of coronavirus in bats from Eastern Thailand’</article-title>
, <source>Virology Journal</source>
, <volume>12</volume>
: <fpage>57</fpage>
.<pub-id pub-id-type="pmid">25884446</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B62"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Wang</surname>
<given-names>Q.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2014</year>
) <article-title>‘Bat origins of MERS-CoV supported by bat coronavirus HKU4 usage of human receptor CD26’</article-title>
, <source>Cell Host & Microbe</source>
, <volume>16</volume>
: <fpage>328</fpage>
–<lpage>37</lpage>
.<pub-id pub-id-type="pmid">25211075</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B63"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Wang</surname>
<given-names>W.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2015</year>
) <article-title>‘Discovery, diversity and evolution of novel coronaviruses sampled from rodents in China’</article-title>
, <source>Virology</source>
, <volume>474</volume>
: <fpage>19</fpage>
–<lpage>27</lpage>
.<pub-id pub-id-type="pmid">25463600</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B64"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Watanabe</surname>
<given-names>S.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2010</year>
) <article-title>‘Bat Coronaviruses and experimental infection of bats, the Philippines’</article-title>
, <source>Emerging Infectious Diseases</source>
, <volume>16</volume>
: <fpage>1217</fpage>
–<lpage>23</lpage>
.<pub-id pub-id-type="pmid">20678314</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B65"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Woo</surname>
<given-names>P. C.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2006</year>
) <article-title>‘Molecular diversity of coronaviruses in bats’</article-title>
, <source>Virology</source>
, <volume>351</volume>
: <fpage>180</fpage>
–<lpage>7</lpage>
.<pub-id pub-id-type="pmid">16647731</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B66"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Woo</surname>
<given-names>P. C.</given-names>
</name>
</person-group>
, <etal></etal>
(<year>2007</year>
) <article-title>‘Comparative analysis of twelve genomes of three novel group 2c and group 2d coronaviruses reveals unique group and subgroup features’</article-title>
, <source>Journal of Virology</source>
, <volume>81</volume>
: <fpage>1574</fpage>
–<lpage>85</lpage>
.<pub-id pub-id-type="pmid">17121802</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B67"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Woo</surname>
<given-names>P. C.</given-names>
</name>
</person-group>
, <etal></etal>
(<year>2009</year>
) <article-title>‘Coronavirus diversity, phylogeny and interspecies jumping’</article-title>
, <source>Experimental Biology and Medicine</source>
, <volume>234</volume>
: <fpage>1117</fpage>
–<lpage>27</lpage>
.<pub-id pub-id-type="pmid">19546349</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B68"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Woo</surname>
<given-names>P. C.</given-names>
</name>
</person-group>
, <etal></etal>
(<year>2012</year>
) <article-title>‘Genetic relatedness of the novel human group C betacoronavirus to Tylonycteris bat coronavirus HKU4 and Pipistrellus bat coronavirus HKU5’</article-title>
, <source>Emerging Microbes and Infection</source>
, <volume>1</volume>
: <fpage>e35</fpage>
.</mixed-citation>
</ref>
<ref id="vex012-B69"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Yang</surname>
<given-names>Y.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2014</year>
) <article-title>‘Receptor usage and cell entry of bat coronavirus HKU4 provide insight into bat-to-human transmission of MERS coronavirus’</article-title>
, <source>Proceedings of the National Academy of Sciences of the United States of America</source>
, <volume>111</volume>
: <fpage>12516</fpage>
–<lpage>21</lpage>
.<pub-id pub-id-type="pmid">25114257</pub-id>
</mixed-citation>
</ref>
<ref id="vex012-B70"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name name-style="western"><surname>Zaki</surname>
<given-names>A. M.</given-names>
</name>
</person-group>
<etal></etal>
(<year>2012</year>
) <article-title>‘Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia’</article-title>
, <source>The 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/SrasV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 0015818 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 0015818 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Sante |area= SrasV1 |flux= Pmc |étape= Corpus |type= RBID |clé= |texte= }}
![]() | This area was generated with Dilib version V0.6.33. | ![]() |