A strategy to assess spillover risk of bat SARS-related coronaviruses in Southeast Asia
=======================================================================================

* Cecilia A. Sánchez
* Hongying Li
* Kendra L. Phelps
* Carlos Zambrana-Torrelio
* Lin-Fa Wang
* Kevin J. Olival
* Peter Daszak

## Abstract

Emerging diseases caused by coronaviruses of likely bat origin (e.g. SARS, MERS, SADS and COVID-19) have disrupted global health and economies for two decades.

Evidence suggests that some bat SARS-related coronaviruses (SARSr-CoVs) could infect people directly, and that their spillover is more frequent than previously recognized. Each zoonotic spillover of a novel virus represents an opportunity for evolutionary adaptation and further spread; therefore, quantifying the extent of this “hidden” spillover may help target prevention programs. We derive biologically realistic range distributions for known bat SARSr-CoV hosts and quantify their overlap with human populations. We then use probabilistic risk assessment and data on human-bat contact, human SARSr-CoV seroprevalence, and antibody duration to estimate that ∼400,000 people (median: ∼50,000) are infected with SARSr-CoVs annually in South and Southeast Asia. These data on the geography and scale of spillover can be used to target surveillance and prevention programs for potential future bat-CoV emergence.

Keywords
*   area of habitat
*   Chiroptera
*   coronavirus
*   human-bat interface
*   probabilistic risk assessment
*   risk mapping
*   sarbecovirus
*   SARSr-CoV
*   spillover
*   COVID-19
*   SARS-CoV-2
*   zoonoses

## Introduction

Emerging coronaviruses (CoVs) of wildlife origin have significantly disrupted global health security and economies during the last two decades1,2. Severe acute respiratory syndrome (SARS) and Middle East respiratory syndrome (MERS) CoVs caused significant human morbidity and mortality in 2002 and 2012 respectively3,4. Swine Acute Diarrheal Syndrome CoV caused substantial mortality in pigs in southern China during 2016 and 20195,6. The emergence of SARS-CoV-2 in 2019 led to the current COVID-19 pandemic that has caused millions of cases and deaths, with economic loss likely to be in the tens of trillions of US dollars2,7. Efforts to increase preparedness and improve surveillance for emerging coronaviruses therefore represent a priority for global health programs8.

Phylogenetic analysis suggests that SARS-CoV, MERS-CoV, SADS-CoV, and SARS-CoV-2 originate within CoV lineages from bat reservoir hosts6,9-11. The initial spillovers of SARS and MERS into human populations are thought to have occurred via intermediate hosts (palm civets and dromedary camels, respectively12,13). However, the role of civets in the emergence of SARS is uncertain, and other bat SARSr-CoVs can directly infect human cells, including airway epithelial cells, and thus have potential to spill over directly from bats to humans14-16. It is uncertain what proportion of bat SARSr-CoVs can infect people either directly or indirectly (via an intermediate or amplifier host). However, serological evidence of prior infection with SARSr-CoVs in communities living near bat populations in China prior to the emergence of COVID-19, including in people who reported no contact with SARSr-CoV intermediate hosts, suggests direct bat-to-human transmission might occur in some regions17,18. Direct bat-to-human spillover events may occur more frequently than has been reported, but go unrecognized because they cause mild symptoms, cause symptoms that are similar to other infections, result in small numbers of cases, or lack sustained chains of human-to-human transmission. However, every wildlife-to-human spillover event represents an opportunity for viral adaptation permitting human-to-human spread19-22. Quantifying the extent of these undetected spillovers could therefore be important to identifying risk of future epidemics or pandemics.

Surveys of bats in China have revealed high diversity of SARSr-CoVs, and often high prevalence (5-10%) in rhinolophid and hipposiderid species that are widely distributed, abundant, resilient to habitat perturbation and are synanthropic (have contact and often interactions with human populations)23-25. Many of the bat species and genera known to harbor these β-CoVs occur in Southeast Asia, a hotspot of bat diversity with 441 species reported, 115 (around a quarter) of which are rhinolophids or hipposiderids26. While some related β-CoVs have been reported outside China (including within Africa), phylogenetic analyses of SARSr-CoVs indicate that regions of south China and parts of neighboring countries (Myanmar, Lao PDR, and Vietnam) act as areas of evolutionary diversification with high diversity of these viruses. However, sampling in China has been far more intense than in nearby Southeast Asian or South Asian countries, and diversity of these viruses is also likely high in these regions23. Furthermore, many of these less well-sampled countries are undergoing dynamic social and environmental changes correlated with zoonotic emergence (e.g. rapid human population growth, movement of rural residents to urban centers, extensive wildlife farming and trade, and rapid land conversion from pristine forested habitats to agricultural monocultures), and thus might represent hitherto unreported hotspots for coronavirus spillover risk23,27-32.

In this study, we use host distribution modeling, ecological, and epidemiological data to estimate the geographic distribution of SARSr-CoV exposure risk, and likely rate of unreported zoonotic spillover in China, South and Southeast Asia. Our results provide the most detailed estimates to date of the distribution and species richness of bat hosts of SARSr-CoVs, and suggest that human exposure to and spillover of SARSr-CoVs may be substantially underestimated. Our approach provides proof of concept for systematic risk assessment of zoonotic spillover, and a strategy to identify key geographic areas that can be prioritized for targeted surveillance of wildlife, livestock, and humans. Given the challenges of identifying the origins of COVID-19 and pathways by which SARS-CoV-2 spilled over to people33,34, this approach may also aid efforts to identify the geographic sites where spillover first occurred.

## Results

We assembled a list of 23 known SARSr-CoV bat host species, mainly members of the *Rhinolophidae* and *Hipposideridae*, occurring in the geographic region of interest (Table S1). We derived an area of habitat (AOH) for each species (Fig. 1, also see rationale for using this approach in Methods), by refining the IUCN geographic range of each species using habitat suitability, elevation limits, and the boundaries of a broadly defined region including parts of South Asia, China and Southeast Asian countries. Removing unsuitable areas of bat occurrence within the IUCN geographic ranges highly improved assessment of species distribution. For example, the reduction in area from the original IUCN geographic range to the more refined AOH for each species ranged from 42% (*Rhinolophus malayanus*) to almost 100% (*R. hipposideros*), with a median of 68% reduction (Fig. S1). Notably, the reduction in area was 55% for *R. rex*, the only species in our list assessed as endangered by the IUCN Red List.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/14/2021.09.09.21263359/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263359/F1)

Figure 1. Representative illustration of IUCN range refinement by habitat type and elevation limits to create area of habitat (AOH).
**a** The IUCN range of *Rhinolophus affinis*, outlined in black, is overlaid on a map of terrestrial habitat types, and the habitat map is restricted to only the species range. **b** Only suitable habitat types within the species range are retained. **c** After further restricting habitat by the elevation limits of the species, the remaining habitat represents the AOH of *R. affinis*. The color scheme follows that of 72,73, with the addition of gray used for habitat type 7.1 (caves; carbonate rock outcrops used as proxy). Names of habitat types are found in Fig. 2b and Table S1.

We validated species AOHs using cleaned occurrence records downloaded from the Global Biodiversity Information Facility (GBIF; see Methods for details of data cleaning and validation). After data cleaning, no occurrence records remained for three species (*Hipposideros pratti, R. ferrumequinum*, and *R. hipposideros*), while 1-179 (median = 13) occurrence points remained for all other species (Table S2). In validating each species’ AOH with GBIF occurrence points, we found among species with ≥1 occurrence point, the median percent of points (buffered by 5 km) that overlapped a species’ AOH was 64% (Table S2). Among species with at least 30 occurrence points (*n* = 7), the median overlap was 84%. Overlap was >80% for five species: *R. creaghi* (31/31, 100%), *R. malayanus* (8/8, 100%), *Chaerephon plicatus* (45/53, 85%), *H. galeritus* (37/44, 84%), and *H. larvatus* (41/49, 84%).

The size of individual species AOHs varied widely, but generally, the largest AOHs encompassed the most people (Fig. 2a). For example, the AOH of *R. affinis* was the largest of all species, covering ∼1.9 million km2, and it encompassed ∼132 million people (2nd highest of any species). Two species had AOHs with more limited area but relatively high human population density: *Nyctalus leisleri* (∼262 people/km2) and *R. stheno* (∼226 people/km2; Fig. 2a, Fig. S2). Within species AOHs, forest habitats and ‘carbonate’ (limestone) rock outcrops (used as a proxy for cave distribution) typically comprised the largest proportion of suitable habitat (Fig. 2b). Two species (*R. stheno* and *R. malayanus*) had high proportions (>50%) of artificial habitats, including plantations and arable land.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/14/2021.09.09.21263359/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263359/F2)

Figure 2. Relationships among AOH size, number of people, and habitat proportions.
**a** Scatterplot showing the total number of people living in each AOH versus the total area, for each SARSr-CoV bat host species. A best-fit line was fit through the origin, for which the R2 is displayed. **b** Proportion of each habitat type within species AOHs. Species are listed in order of decreasing AOH size. Color scheme follows that of Jung and colleagues72,73, with the addition of gray for caves (carbonate rock outcrops used as proxy).

The consensus area of all SARSr-CoV bat host species, created by mosaicking the 23 species AOHs, comprised ∼4.5 million km2. We calculated that ∼478 million people live within this consensus area, which covered most of Lao PDR, Cambodia, Thailand, Vietnam, Nepal, Bhutan, peninsular Malaysia, Myanmar, southeast China, and the western islands of Indonesia (Fig. 3a). Bat species distribution was patchier in India, Sri Lanka, East Malaysia, and the Philippines. Species richness ranged from 1-14 species, with the highest richness of SARSr-CoV bat host species in southern China, eastern Myanmar, and northern Lao PDR (Fig. 3a). When we visualized areas with both high host richness and large human populations (‘relative spillover risk’), southern China remained a hotspot, while other areas emerged as important because of their high human population sizes (e.g. Java, parts of northern India; Fig. 3b).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/14/2021.09.09.21263359/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263359/F3)

Figure 3. Hotspot maps of Southeast Asia.
**a** Species richness of SARSr-CoV bat host species in Southeast Asia. This map was created by overlaying AOH maps for all 23 SARSr-CoV bat host species. **b** Bat host species richness * human population count, or “relative spillover risk”. Values were log(x + 1) transformed and then normalized to a 0-1 scale.

After calculating the number of people living in the consensus area of SARSr-CoV bat host species, we further incorporated data from the literature on bat-human contacts, human SARSr-CoV seroprevalence, and SARS antibody duration to estimate the extent of hidden bat-human spillover of SARSr-CoVs in Southeast Asia (see Supplemental Material for details of probabilistic risk assessment methods). We estimated that within the consensus area of SARSr-CoV bat host species, an average of 407,422 people (median: 53,225; range: 1 - 35,632,228) are infected with SARSr-CoVs each year in Southeast Asia (Fig. S3). Sensitivity analyses indicated that the probability of antibody detection given contact with a bat, and the probability of contact with a bat, contributed most to variance in the outcome (Fig. S4). In the process of identifying datasets of value for each step of the probability analysis, we also assessed datasets that are not available or for which quality could be improved, and would likely provide more accurate assessments of risk (Fig. 4).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/09/14/2021.09.09.21263359/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2021/09/14/2021.09.09.21263359/F4)

Figure 4. Conceptual figure depicting key data inputs to better estimate spillover of SARSr-CoVs from bats to humans, and gaps in the available data.

## Discussion

Our paper reports an analytical framework to assess SARSr-CoV spillover risk in a region that includes the site of the first spillover of SARS-CoV, and likely of SARS-CoV-2. By incorporating data on biologically plausible host distributions for all bat species known to harbor SARSr-CoVs in South and Southeast Asia, on human population counts, on behavior of people in the region that drives bat-human contact, and on seroprevalence in people and duration of antibody persistence, we are able to provide detailed maps of potential bat-to-human SARSr-CoV spillover hotspots, as well as estimates of the number of people infected by this viral group annually. Our approach may assist in better identifying regions for targeted surveillance of local communities at risk of spillover, for conducting viral discovery programs to identify novel bat-CoVs, for COVID-19 origins tracing, and for estimating human surveillance targets to identify spillover events earlier and more accurately. All of these are key goals for pandemic preparedness and prevention35-38, and if used to target future surveillance and disease control, may help to reduce the risk of future COVID-like outbreaks.

Our analysis identifies regions in southern China, northeastern Myanmar, Lao PDR, and northern Vietnam as having the highest diversity of SARSr-CoV bat host species. These hotspots of SARSr-CoV bat reservoir host diversity may be particularly fruitful sites for viral discovery of novel SARSr-CoVs, assuming that viral diversity scales with host species diversity39. This finding supports conclusions from prior phylogenetic analyses that particularly diverse SARSr-CoV lineages are found in southern China, and that it may be a center of evolutionary origin for this group of β-CoVs23. It also helps explain the recent identification of multiple strains of SARSr-CoVs in southern China32,40,41 and southeast Asian countries31, despite small sample sizes. It suggests that the less intense sampling in countries bordering southern China has led to an underestimate of the diversity of these viruses there. Given that the bat species known to host the closest relatives of SARS-CoV-2 are found in this region, these maps may also guide efforts to identify the viral clade from which a progenitor virus emerged23,32,34. The map of bat SARSr-CoV host diversity combined with human population density points to hotspots of bat-to-human exposure risk (and potential for human-to-human spread) in southern China and bordering countries, but also in the populous regions of Indonesia. This map may be useful in targeting surveillance to identify SARSr-CoV spillover events in people, including syndromic surveillance for SARS- or COVID-like respiratory disease in communities within these hotspots, as proposed previously18,24,35,39. It may also provide guidance for studies of where initial spillover of the progenitor of SARS-CoV-2 may have occurred, although epidemiological and other data suggest this would most likely have been in China or neighboring countries33,34.

Our estimate that an average of ∼400,000 (median ∼50,000) people are infected with SARSr-CoVs each year in Southeast Asia suggests that bat-to-human SARSr-CoV spillover is common in the region, and is undetected by surveillance programs and clinical studies in the majority of cases. While the data suggest significant levels of exposure, many of the diverse viral strains that infect people in the region each year may not be able to replicate well in people, cause illness, or be transmitted sufficiently among people to cause an outbreak. This has been shown in theoretical models of disease emergence21,42, and supports earlier evidence from studies of non-human primate virus spillover43, that cross-species transmission of novel animal-origin viruses is not the rate-limiting step in pandemic viral emergence. However, given the relatively large number of people likely to be infected each year with bat-CoVs, it is plausible that illnesses or clusters of cases due to novel bat-CoV infection occur regularly within the region, and are either not reported, or otherwise missed by clinical surveillance. Evidence of underreporting has been demonstrated for other bat-origin viral infections. For example, targeted syndromic surveillance of encephalitis patients in a small number of clinics in Bangladesh showed that Nipah virus causes outbreaks annually with an overall mortality rate of ∼70%, despite it only recently being reported from the country44. Efforts to increase surveillance for novel SARSr-CoVs (and other bat-origin viruses) in clinical cohorts, particularly using syndromic surveillance, may identify the rate of missed cases, and pre-empt larger scale outbreaks. Estimating the true rate of spillover of previously unknown, potentially zoonotic animal-origin viruses is difficult without serological or genomic surveillance data. For most viruses, the duration of infection in humans is relatively short, e.g. the infectious period for COVID-19 is 2-3 weeks45, and if spillover is rare, PCR surveillance is unlikely to give valuable data on spillover rates due to lack of positives. One exception is for viruses with long infectious periods such as lentiviruses and some retroviruses, and a previous genomic surveillance study of wildlife hunters in Africa was able to find 10/1099 (0.9%) prevalence of non-human primate origin simian foamy viruses46. Because detectable antibodies or T cell responses are more long-lived, serological surveys or activated T cell testing represent more valid strategies to estimate spillover rates.

Our calculation of undetected spillover represents the first published attempt to identify the rate of spillover of novel coronaviruses from bats to people. It relies on a number of input variables and we attempted to account for uncertainty by assigning a probability distribution to each variable, but recognize that further data would likely help refine estimates of SARSr-CoV spillover rates. Perhaps most critical is a lack of information on the role of intermediate hosts in the emergence of SARSr-CoVs. It has been postulated that civets and other commonly farmed and traded mammalian species played a role in the emergence of SARS-CoV by acting as efficient amplifier hosts within wildlife farms and markets12,47-50. SARS-CoV-2 can infect civets, raccoon dogs, and other mammals commonly farmed and traded for food in the region51,52. SARS-CoV-2 has also caused significant outbreaks in animals bred for fur (e.g. mink, raccoon dogs) that have in some cases led to transmission to people53,54. These and other data led an international team to conclude that the most likely pathway for COVID emergence was from bats to people through a farmed mammalian intermediate host33. Accurate data are not available on the number of wildlife farms and potential intermediate hosts bred each year, or on the market systems that supply live animals into cities within China and Southeast Asia. It was estimated that 14 million people were employed in the wildlife farming system within China alone during 2016, in an industry worth $77 billion annually55,56. Additionally, some of these potential intermediate hosts occur naturally in the region (e.g. pangolins, civets), and other livestock species that are extremely common (e.g. pigs, cattle, rabbits) are susceptible to SARSr-CoVs either naturally or experimentally52. Spillover of SARSr-CoVs may therefore be substantially skewed to people who have high exposure to these species, and this would likely have been missed in the relatively small serological surveys upon which our analyses are currently based. Thus, better estimates of the role of farmed and traded intermediate hosts are likely to substantially increase the estimated spillover rates of SARSr-CoVs across the region.

Other data could improve our estimates. These include expanded survey data on the background seroprevalence of bat SARSr-CoVs in people across different populations and geographies, analysis of how specific risk behaviors correlate with seroprevalence, and studies of how serological titers and duration of antibody persistence relate to severity of infection. We estimated bat-human contact using data from previous human-animal contact surveys and ethnographic investigations, but contacts likely vary widely across the region due to different cultural practices and traditions57,58. The type of contact may be critical to assessing the risk of transmission (e.g. hunting and butchering vs. living near a bat colony), and some types of contact may be unrecognized and therefore unreported (e.g. exposure to feces or urine on surfaces). Some of the contact data we identified in our systematic literature search (see Supplemental Material for details) were from populations targeted because of known high risk of bat contact due to occupational or environmental exposure and are likely upwardly biased due to non-random study design. Expanded bat viral survey data would also likely improve our estimates, first by obtaining more contemporary, accurate species presence records, and also by identifying previously unknown SARSr-CoV hosts, improving species-specific estimates of SARSr-CoV prevalence, and clarifying how environmental (e.g. location, season, year) and host trait (e.g. sex, age, body condition) factors affect viral shedding59-62. Finally, it is also possible that taxonomic errors have occurred in the data we analyzed. Taxonomic standards in viral discovery and surveillance studies vary widely, and point to a need for better taxonomic training of field teams, standardized DNA ‘barcoding’ of hosts, collection of voucher specimens, and closer collaboration among disease ecologists, virologists, field biologists and taxonomists63,64.

Our refinement of species ranges has likely produced the most accurate contemporary picture of SARSr-CoV bat host species occurrence to date, and demonstrates substantial range modification for many species, reflecting rapid land use change throughout the region. The AOHs we produced may be useful for targeting surveillance to key species. For example, only a few species AOHs contained a sizable proportion (≥ 25%) of one or more artificial habitats (Fig. 2b): *R. malayanus* (arable land, plantations), *R. stheno* (plantations, arable land), and *N. leisleri* (pastureland). These species may be more likely to come into contact with humans and could represent a particularly important spillover risk considering that global analyses demonstrate zoonotic host diversity increases in areas with anthropogenically modified habitats65. Cave habitats were classified by the IUCN as suitable for nearly all species in our analyses, and carbonate rock outcrops (used here as a proxy for caves) made up a large proportion of species AOHs. Visiting caves, collecting guano, and using bat guano in crop production are likely particularly high-risk activities given these findings66-69. Our validation process indicated good agreement between species AOHs and GBIF occurrence data, suggesting that the maps we generated accurately reflect species presence for occurrence records collected after 1990. Limitations of AOH maps70 include the potential inaccuracy of the IUCN species ranges71, habitat suitability assignments, and elevation limits. Additionally, the map of terrestrial habitat types we used in our analyses72,73 did not include caves, an important habitat type for many bat species. We used carbonate rock outcrop data as a proxy for cave distribution and this could be ground-truthed.

Our analytical framework provides a strategy that has potential for improving preparedness for emerging diseases and pandemic risk. It has produced maps that can be used to conduct more cost-effective field surveys for viral discovery programs, and estimates of spillover rates that can guide targeted human surveillance to identify clusters of cases of a new CoV infection earlier and help prevent spread. It may also give vital guidance for efforts to identify the reservoir hosts of the SARS-CoV-2 progenitor, and the sites of COVID origins or first emergence34. Our analysis pipeline and framework are based on open-source code and can therefore serve as a resource to update and modify spillover risk maps and estimates as new data become available. These could include serological surveys in people using new diagnostic assays that can detect virus-specific neutralizing antibodies to differentiate COVID-19 variants, vaccine strains, and different clades of bat SARSr-CoVs74, or data from ethnographic surveys that identify changes in bat-to-human contact as habitats are increasingly modified or laws to reduce hunting or wildlife trade reduce consumption are enforced. Finally, our framework can be rapidly adapted for spillover risk assessment of other viral groups, such as the HKU-2/SADS-CoV α-CoVs that have recently been found able to infect primary human airway epithelial cells *in vitro*, and therefore pose a heightened spillover risk, or any of the other ∼25 viral families that include known zoonoses39.

## Methods

### 1. Compilation of SARSr-CoV bat host data

We identified all bat species from which molecular evidence of SARSr-CoV infection had been reported, and for which associated sequence confirmation data were available. We supplemented a recently compiled list75 with hosts listed in other more recent publications23,30,31,41,76-78. We removed duplicate records and hosts only identified to genus level. We considered *Rhinolophus paradoxolophus* to be a subspecies of *R. rex*79. We excluded *R. monoceros* because, although recent work has retained it as a distinct species80, no recent assessment of the species has been published by the IUCN Red List of Threatened Species. *Hipposideros pomona* and *H. gentilis* were split in 2018, with *H. pomona* restricted to a small area in southern India while *H. gentilis* is broadly distributed across Southeast Asia81. Although *H. pomona* was listed as a host species in 23 from field sampling in China, we therefore substituted *H. gentilis* in our analyses. We chose to not include several bat species that were predicted to be SARSr-CoV reservoirs using host-trait or species network models82,83, as these predictions remain unvalidated, and vary in reliability across models.

We restricted the compiled list of SARSr-CoV bat hosts to those with geographic distributions either partially or entirely within Southeast Asia, and included the following countries and administrative regions: Bangladesh, Bhutan, Brunei, Cambodia, China, Hong Kong SAR, Macao SAR, India, Indonesia, Lao People’s Democratic Republic, Malaysia, Myanmar, Nepal, the Philippines, Singapore, Sri Lanka, Thailand, Timor-Leste, and Vietnam. See Table S1 for the finalized list of SARSr-CoV bat hosts. Note that while we use the term ‘Southeast Asia’, the region we considered has broader range than is used for political definitions of the region and includes China and parts of South Asia, due to the extensive range of some of the bat host species.

### 2. Calculation and validation of host area of habitat

All analyses were performed in the R statistical environment v. 4.0.3 (R Core Team 2020). To derive a biologically realistic area of habitat (AOH) within our Southeast Asia region for each SARSr-CoV bat host reservoir species, we refined species ranges by incorporating data on habitat suitability and elevation limits70. We used an AOH approach rather than ecological niche modeling (ENM) because ENMs rely on presence data, and this was difficult to find within GBIF, particularly when we excluded older records that may not reflect current distributions. Secondly, ENMs have the potential to overestimate range sizes and are thus not useful for identifying more precisely areas for surveillance, for example. To derive the AOH for each species, we downloaded its geographic range from the IUCN Red List of Threatened Species84 (last update: March 25, 2021), overlaid it onto a raster map of terrestrial habitat types72,73, and selected areas of the map that occurred within its range (Fig. 1a). We then selected areas of suitable habitat (as determined by IUCN) for each species (Fig. 1b, Table S1; see Supplemental Material for further details). Finally, we overlaid elevation data (Shuttle Radar Topography Mission data obtained with the getData function of the “raster” package85) onto each species’ range and habitat map and selected areas that fell within a species’ elevation limits (Fig. 1c; Table S1); these remaining areas represented the species’ AOH. Using the “raster” package, we calculated the size (in km2) of each species’ AOH and compared it to the size of its original IUCN range. Note that all analyses were restricted by the regional boundaries of Southeast Asia as defined above.

We validated the AOH of each species using occurrence data downloaded from the Global Biodiversity Information Facility (GBIF) using the “rgbif” package86,87. We cleaned GBIF data by removing records with inaccurate or imprecise coordinate data (i.e. no coordinates, identical longitude and latitude, coordinates of 0, coordinate uncertainty > 35 km, coordinates in the ocean, country–coordinate mismatch), records from areas outside our geographically defined region, records of absence rather than presence, and records with unreliable data sources (e.g. fossil specimen, literature, unknown). We removed records within 5km of country capitals, within 1km of country/province centroids, and within 100m of biodiversity institutions. We removed records before 1991, as older records tend to have less precise location data88 and could reflect species ranges that have since shifted. Within each species, we removed records with duplicate coordinates. Cleaning was facilitated with the “CoordinateCleaner” package89. For each species, we buffered each of its occurrence points by five kilometers and calculated the percent of buffered points that overlapped the species’ AOH.

### 3. Estimation of bat-human overlap, bat species richness, and habitat proportions

To identify regions where human populations might be exposed directly or indirectly to SARSr-CoV bat hosts, we overlaid human population count data on each host species’ AOH. We used a 1-km resolution raster of 2020 population count data from WorldPop90 and resampled it with bilinear interpolation using the “gdalUtils” package91 so that its extent and resolution matched that of the habitat raster. We calculated the number of people living within the AOH of each species separately, and divided this by the size of each AOH to obtain average human population density values. As living in areas with high diversity of SARSr-CoV bat hosts and large human populations may increase the likelihood of human-bat contact, viral spillover, and subsequent pathogen spread, we also visualized bat species richness and bat species richness multiplied by human population size (which we term “relative spillover risk”) across the consensus map. To examine the relative importance of different habitat types to SARSr-CoV bat hosts, we calculated the proportion (by area) of each habitat type within the AOH of each species.

### 4. Estimating human exposure to and infection with SARSr-CoVs

We used a probabilistic risk assessment to estimate the extent of SARSr-CoV spillover from bats to humans. We assumed that the number of people infected with SARSr-CoVs by bats each year is equal to 1) the number of people who live in the consensus area of SARSr-CoV bat hosts, multiplied by 2) the probability that a human comes into contact with a bat, multiplied by 3) the probability that a bat-human contact leads to a serologically detectable human infection, multiplied by 4) the probability that serological detection is due to an infection within the previous year. To inform our choice of distribution for each input variable, we gathered data from the literature on human-bat contacts, human SARSr-CoV seroprevalence, and human SARS antibody duration (see the Supplemental Material for details). We accounted for uncertainty and/or variation associated with the input variables by assigning a probability distribution (rather than a single fixed value) to each. We performed Latin hypercube sampling with the “lhs” package92 to generate 400,000 sets of input combinations and calculated the total number of people infected for each set of inputs, creating an output distribution. We then calculated summary statistics (mean, median, range) for this distribution. Finally, we performed a global sensitivity analysis by calculating Sobol sensitivity analyses to understand the relative contribution of each input variable to the outcome.

## Supporting information

Supplemental material [[supplements/263359_file02.docx]](pending:yes)

## Data Availability

Data Availability The datasets analyzed during the current study are available on Zenodo: [https://zenodo.org/record/5251726#.YSWivEApAUE](https://zenodo.org/record/5251726#.YSWivEApAUE). Code Availability Code to replicate the analyses and figures is available on Zenodo: [https://zenodo.org/record/5251726#.YSWivEApAUE](https://zenodo.org/record/5251726#.YSWivEApAUE). 

[https://zenodo.org/record/5251726#.YSWivEApAUE](https://zenodo.org/record/5251726#.YSWivEApAUE) 

## Data Availability

The datasets analyzed during the current study are available on Zenodo: [https://zenodo.org/record/5251726#.YSWivEApAUE](https://zenodo.org/record/5251726#.YSWivEApAUE).

## Code Availability

Code to replicate the analyses and figures is available on Zenodo: [https://zenodo.org/record/5251726#.YSWivEApAUE](https://zenodo.org/record/5251726#.YSWivEApAUE).

## Author Contributions

P.D. conceived of the study, C.A.S., H.L., K.L.P., C.Z.-T., and K.J.O. designed the methods, C.A.S. conducted the analyses, H.L. performed the literature search, C.A.S. and H.L. designed figures and tables, C.A.S. wrote the first draft of the manuscript, and all authors provided comments and edits to the manuscript throughout the writing process. All authors approved the submitted manuscript and agree to be held accountable for their own contributions and ensure that questions related to the accuracy or integrity of any part of the work are appropriately investigated, resolved, and the resolution documented in the literature.

## Competing Interests

P.D. served as a member of the WHO-China joint study on COVID-19 origins, is a current member of the Lancet COVID-19 Commission Taskforce on the Origins and Early Spread of COVID-19 and One Health Solutions to Future Pandemics, and was Chair of the IPBES Pandemics and Biodiversity Workshop. P.D. has made numerous public statements both independently, and as part of these groups, on the likely origins of COVID-19 and the need to assess risk and prioritize targeted surveillance for future disease emergence. L.F.W. serves on multiple committees for WHO, FAO and OIE on COVID-19 including assay and vaccine development and animal models; has ongoing research investigating the origin of COVID-19; and has made statements on this issue to the media.

## Acknowledgments

We acknowledge funding from the National Institute of Allergy and Infectious Diseases of the National Institutes of Health (NIAID-CREID U01AI151797 “EID-SEARCH”), the Samuel Freeman Charitable Trust, Pamela Thye, The Wallace Fund, and an Anonymous Donor c/o Schwab Charitable. L.F.W. is supported by grants from NRF (NRF2012NRF-CRP001-056; NRF2016NRF-NSFC002-013), NMRC (MOH-OFIRG19MAY-0011; COVID19RF-003) and MOE (MOE2019-T2-2-130) of Singapore. The authors thank Cadhla Firth, Noam Ross, and Emma Mendelsohn for useful discussions.

*   Received September 9, 2021.
*   Revision received September 9, 2021.
*   Accepted September 14, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  Lee, J.-W. & McKibbin, W. J. Globalization and Disease: The Case of SARS. Asian Economic Papers 3, 113–131, doi:10.1162/1535351041747932 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1162/1535351041747932&link_type=DOI) 

2.  Cutler, D. M. & Summers, L. H. The COVID-19 Pandemic and the $16 Trillion Virus. JAMA-J. Am. Med. Assoc. 324, 1495–1496, doi:10.1001/jama.2020.19759 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2020.19759&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33044484&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

3.  Peiris, J. S. M., Guan, Y. & Yuen, K. Y. Severe acute respiratory syndrome. Nature Medicine 10, S88–S97, doi:10.1038/nm1143 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm1143&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15577937&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225733900005&link_type=ISI) 

4.  Raj, V. S., Osterhaus, A. D. M. E., Fouchier, R. A. M. & Haagmans, B. L. MERS: emergence of a novel human coronavirus. Current Opinion in Virology 5, 58–62, doi:[https://doi.org/10.1016/j.coviro.2014.01.010](https://doi.org/10.1016/j.coviro.2014.01.010) (2014).
    
    

5.  Zhou, L. et al. The re-emerging of SADS-CoV infection in pig herds in Southern China. Transbound Emerg Dis 66, 2180–2183, doi:10.1111/tbed.13270 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/tbed.13270&link_type=DOI) 

6.  Zhou, P. et al. Fatal swine acute diarrhoea syndrome caused by an HKU2-related coronavirus of bat origin. Nature 556, 255–258, doi:10.1038/s41586-018-0010-9 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0010-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29618817&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

7.  Dong, E., Du, H. & Gardner, L. An interactive web-based dashboard to track COVID-19 in real time. The Lancet Infectious Diseases 20, 533–534, doi:[https://doi.org/10.1016/S1473-3099(20)30120-1](https://doi.org/10.1016/S1473-3099(20)30120-1) (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(20)30120-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

8.  Daszak, P., Keusch, G. T., Phelan, A. L., Johnson, C. K. & Osterholm, M. T. Infectious Disease Threats: A Rebound To Resilience. Health Aff. 40, 204–211, doi:10.1377/hlthaff.2020.01544 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1377/hlthaff.2020.01544&link_type=DOI) 

9.  Li, W. D. et al. Bats are natural reservoirs of SARS-like coronaviruses. Science 310, 676–679, doi:10.1126/science.1118391 (2005).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMTAvNTc0OC82NzYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

10. Anthony, S. et al. Further evidence for bats as the evolutionary source of Middle East respiratory syndrome coronavirus. MBio 8, e00373–00317 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.00373-17&link_type=DOI) 

11. Zhou, P. et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 579, 270–273, doi:10.1038/s41586-020-2012-7 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2012-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32015507&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

12. 1.   James E. Childs, 
    2.   John S. Mackenzie, & 
    3.   Jürgen A. Richt
    
    Wang, L. F. & Eaton, B. T. in Wildlife and Emerging Zoonotic Diseases: The Biology, Circumstances and Consequences of Cross-Species Transmission (eds  James E. Childs,  John S. Mackenzie, &  Jürgen A. Richt) 325–344 (Springer Berlin Heidelberg, 2007).
    
    

13. Dudas, G., Carvalho, L. M., Rambaut, A. & Bedford, T. MERS-CoV spillover at the camel-human interface. eLife 7, e31257, doi:10.7554/eLife.31257 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.31257&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29336306&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

14. Ge, X. Y. et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature 503, 535-+, doi:10.1038/nature12711 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12711&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24172901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000327464200045&link_type=ISI) 

15. Menachery, V. D. et al. A SARS-like cluster of circulating bat coronaviruses shows potential for human emergence. Nat Med 21, 1508–1513, doi:10.1038/nm.3985 [http://www.nature.com/nm/journal/v21/n12/abs/nm.3985.html#supplementary-information](http://www.nature.com/nm/journal/v21/n12/abs/nm.3985.html#supplementary-information) (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm.3985&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26552008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

16. Menachery, V. D. et al. SARS-like WIV1-CoV poised for human emergence. Proceedings of the National Academy of Sciences 113, 3048, doi:10.1073/pnas.1517719113 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTEzLzExLzMwNDgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

17. Li, H. et al. Human-animal interactions and bat coronavirus spillover potential among rural residents in Southern China. Biosafety and Health 1, 84–90, doi:[https://doi.org/10.1016/j.bsheal.2019.10.004](https://doi.org/10.1016/j.bsheal.2019.10.004) (2019).
    
    

18. Wang, N. et al. Serological evidence of bat SARS-related coronavirus infection in humans, China. Virologica Sinica 33, 104–107, doi:10.1007/s12250-018-0012-7 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12250-018-0012-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29500691&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

19. Parrish, C. R. et al. Cross-Species Virus Transmission and the Emergence of New Epidemic Diseases. Microbiology and Molecular Biology Reviews 72, 457–470, doi:10.1128/mmbr.00004-08 (2008).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoibW1iciI7czo1OiJyZXNpZCI7czo4OiI3Mi8zLzQ1NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA5LzE0LzIwMjEuMDkuMDkuMjEyNjMzNTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

20. Wasik, B. R. et al. Onward transmission of viruses: how do viruses emerge to cause epidemics after spillover? Philosophical transactions of the Royal Society of London. Series B, Biological sciences 374, 20190017–20190017, doi:10.1098/rstb.2019.0017 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.2019.0017&link_type=DOI) 

21. Lloyd-Smith, J. O. et al. Epidemic dynamics at the human-animal interface. Science 326, 1362–1367, doi:10.1126/science.1177345 (2009).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMjYvNTk1OC8xMzYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzM1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

22. Gray, G. C., Robie, E. R., Studstill, C. J. & Nunn, C. L. Mitigating Future Respiratory Virus Pandemics: New Threats and Approaches to Consider. Viruses 13, 637 (2021).
    
    

23. Latinne, A. et al. Origin and cross-species transmission of bat coronaviruses in China. Nat. Commun. 11, 4235, doi:10.1038/s41467-020-17687-3 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-17687-3&link_type=DOI) 

24. Hu, B. et al. Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus. PLoS pathogens 13, doi:10.1371/journal.ppat.1006698 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1006698&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

25. McFarlane, R., Sleigh, A. & McMichael, T. Synanthropy of Wild Mammals as a Determinant of Emerging Infectious Diseases in the Asian-Australasian Region. EcoHealth 9, 24–35, doi:10.1007/s10393-012-0763-9 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10393-012-0763-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22526750&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

26. International Union for Conservation of Nature (IUCN). The IUCN Red List of Threatened Species. Version 2020-1. (2020).
    
    

27. Coker, R. J., Hunter, B. M., Rudge, J. W., Liverani, M. & Hanvoravongchai, P. Emerging infectious diseases in southeast Asia: regional challenges to control. The Lancet 377, 599–609, doi:10.1016/S0140-6736(10)62004-1 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(10)62004-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21269678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000287556700030&link_type=ISI) 

28. Horby, P. W., Pfeiffer, D. & Oshitani, H. Prospects for Emerging Infections in East and Southeast Asia 10 Years after Severe Acute Respiratory Syndrome. Emerging Infectious Disease journal 19, 853, doi:10.3201/eid1906.121783 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1906.121783&link_type=DOI) 

29. Rulli, M. C., D’Odorico, P., Galli, N. & Hayman, D. T. S. Land-use change and the livestock revolution increase the risk of zoonotic coronavirus transmission from rhinolophid bats. Nature Food, doi:10.1038/s43016-021-00285-x (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s43016-021-00285-x&link_type=DOI) 

30. Hul, V. et al. A novel SARS-CoV-2 related coronavirus in bats from Cambodia. bioRxiv, 2021.2001.2026.428212, doi:10.1101/2021.01.26.428212 (2021).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMS4wMS4yNi40MjgyMTJ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA5LzE0LzIwMjEuMDkuMDkuMjEyNjMzNTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

31. Wacharapluesadee, S. et al. Evidence for SARS-CoV-2 related coronaviruses circulating in bats and pangolins in Southeast Asia. Nat. Commun. 12, 972, doi:10.1038/s41467-021-21240-1 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-21240-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33563978&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

32. Zhou, H. et al. Identification of novel bat coronaviruses sheds light on the evolutionary origins of SARS-CoV-2 and related viruses. Cell, doi:[https://doi.org/10.1016/j.cell.2021.06.008](https://doi.org/10.1016/j.cell.2021.06.008) (2021).
    
    

33. Joint WHO-China study. WHO-convened global study of origins of SARS-CoV-2: China part. (2021).
    
    

34. Holmes, E. C. et al. The Origins of SARS-CoV-2: A Critical Review. Cell, doi:10.1016/j.cell.2021.08.017.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2021.08.017&link_type=DOI) 

35. Dobson, A. P. et al. Ecology and economics for pandemic prevention. Science 369, 379–381, doi:10.1126/science.abc3189 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjkvNjUwMi8zNzkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

36. Petrovan, S. O. et al. Post COVID-19: a solution scan of options for preventing future zoonotic epidemics. Biological Reviews n/a, doi:[https://doi.org/10.1111/brv.12774](https://doi.org/10.1111/brv.12774).
    
    

37. Roche, B. et al. Was the COVID-19 pandemic avoidable? A call for a “solution-oriented” approach in pathogen evolutionary ecology to prevent future outbreaks. Ecology Letters 23, 1557–1560, doi:[https://doi.org/10.1111/ele.13586](https://doi.org/10.1111/ele.13586) (2020).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

38. Naguib, M. M., Ellström, P., Järhult, J. D., Lundkvist, Å. & Olsen, B. Towards pandemic preparedness beyond COVID-19. The Lancet Microbe 1, e185–e186, doi:10.1016/S2666-5247(20)30088-4 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2666-5247(20)30088-4&link_type=DOI) 

39. Carroll, D. et al. The global virome project. Science 359, 872–874 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTkvNjM3OC84NzIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

40. Li, L.-l. et al. A novel SARS-CoV-2 related virus with complex recombination isolated from bats in Yunnan province, China. bioRxiv, 2021.2003.2017.435823, doi:10.1101/2021.03.17.435823 (2021).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMS4wMy4xNy40MzU4MjN2MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA5LzE0LzIwMjEuMDkuMDkuMjEyNjMzNTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

41. Zhou, H. et al. A Novel Bat Coronavirus Closely Related to SARS-CoV-2 Contains Natural Insertions at the S1/S2 Cleavage Site of the Spike Protein. Current Biology 30, 2196-2203.e2193, doi:10.1016/j.cub.2020.05.023 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cub.2020.05.023&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32416074&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

42. Wolfe, N. D., Dunavan, C. P. & Diamond, J. Origins of major human infectious diseases. Nature 447, 279–283, doi:10.1038/nature05775 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature05775&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17507975&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246520300036&link_type=ISI) 

43. Wolfe, N. D. et al. Emergence of unique primate T-lymphotropic viruses among central African bushmeat hunters. Proceedings of the National Academy of Sciences of the United States of America 102, 7994–7999 (2005).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTAyLzIyLzc5OTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

44. Nikolay, B. et al. Transmission of Nipah Virus — 14 Years of Investigations in Bangladesh. New England Journal of Medicine 380, 1804–1814, doi:10.1056/NEJMoa1805376 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1805376&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31067370&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

45. Byrne, A. W. et al. Inferred duration of infectious period of SARS-CoV-2: rapid scoping review and analysis of available evidence for asymptomatic and symptomatic COVID-19 cases. BMJ Open 10, e039856, doi:10.1136/bmjopen-2020-039856 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMjoiMTAvOC9lMDM5ODU2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzM1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

46. Wolfe, N. D. et al. Naturally acquired simian retrovirus infections in central African hunters. Lancet 363, 932–937 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(04)15787-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15043960&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220308600009&link_type=ISI) 

47. Bell, D., Roberton, S. & Hunter, P. R. Animal origins of SARS coronavirus: possible links with the international trade in small carnivores. Philosophical Transactions of the Royal Society of London Series B-Biological Sciences 359, 1107–1114, doi:10.1098/rstb.2004.1492 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.2004.1492&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15306396&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222688200009&link_type=ISI) 

48. He, J. F. et al. Molecular evolution of the SARS coronavirus during the course of the SARS epidemic in China. Science 303, 1666–1669, doi:10.1126/science.1092002 (2004).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMDMvNTY2NC8xNjY2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzM1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

49. Tu, C. et al. Antibodies to SARS coronavirus in civets. Emerg. Infect. Dis. 1, 2244–2248 (2004).
    
    

50. Guan, Y. et al. Isolation and characterization of viruses related to the SARS coronavirus from animals in southern China. Science 302, 276–278 (2003).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMDIvNTY0My8yNzYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

51. Freuling, C. et al. Susceptibility of Raccoon Dogs for Experimental SARS-CoV-2 Infection. Emerging Infectious Disease journal 26, 2982, doi:10.3201/eid2612.203733 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2612.203733&link_type=DOI) 

52. OIE. Infection with SARS-CoV-2 in animals. (OIE, Paris, 2020).
    
    

53. Oreshkova, N. et al. SARS-CoV-2 infection in farmed minks, the Netherlands, April and May 2020. Eurosurveillance 25, 2001005 (2020).
    
    

54. Oude Munnink, B. B. et al. Transmission of SARS-CoV-2 on mink farms between humans and mink and back to humans. Science 371, 172, doi:10.1126/science.abe5901 (2021).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNzEvNjUyNS8xNzIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOS8xNC8yMDIxLjA5LjA5LjIxMjYzMzU5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

55. IPBES et al. Workshop Report on Biodiversity and Pandemics of the Intergovernmental Platform on Biodiversity and Ecosystem Services., (IPBES, Bonn, Germany, 2020).
    
    

56. Chinese Academy of Engineering. Report on Sustainable Development Strategy of China’s Wildlife Farming Industry. (2017).
    
    

57. Mildenstein, T., Tanshi, I. & Racey, P. A. in Bats in the Anthropocene: conservation of bats in a changing world 325–375 (Springer, 2016).
    
    

58. Low, M.-R. et al. Bane or Blessing? Reviewing Cultural Values of Bats across the Asia-Pacific Region. Journal of Ethnobiology 41, 18–34, doi:10.2993/0278-0771-41.1.18 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2993/0278-0771-41.1.18&link_type=DOI) 

59. Wacharapluesadee, S. et al. Longitudinal study of age-specific pattern of coronavirus infection in Lyle’s flying fox (Pteropus lylei) in Thailand. Virology Journal 15, 38, doi:10.1186/s12985-018-0950-6 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12985-018-0950-6&link_type=DOI) 

60. Luo, Y. et al. Longitudinal Surveillance of Betacoronaviruses in Fruit Bats in Yunnan Province, China During 2009–2016. Virologica Sinica 33, 87–95, doi:10.1007/s12250-018-0017-2 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s12250-018-0017-2&link_type=DOI) 

61. Maganga, G. D. et al. Genetic diversity and ecology of coronaviruses hosted by cave-dwelling bats in Gabon. Scientific Reports 10, 7314, doi:10.1038/s41598-020-64159-1 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-64159-1&link_type=DOI) 

62. Epstein, J. H. et al. Nipah virus dynamics in bats and implications for spillover to humans. Proceedings of the National Academy of Sciences 117, 29190, doi:10.1073/pnas.2000429117 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzQ2LzI5MTkwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDkvMTQvMjAyMS4wOS4wOS4yMTI2MzM1OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

63. Thompson, C. W. et al. Preserve a Voucher Specimen! The Critical Need for Integrating Natural History Collections in Infectious Disease Studies. mBio 12, e02698–02620, doi:10.1128/mBio.02698-20 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.02698-20&link_type=DOI) 

64. Phelps, K. et al. Bat Research Networks and Viral Surveillance: Gaps and Opportunities in Western Asia. Viruses 11, doi:10.3390/v11030240 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/v11030240&link_type=DOI) 

65. Gibb, R. et al. Zoonotic host diversity increases in human-dominated ecosystems. Nature, 1–5 (2020).
    
    

66. Robertson, K. et al. Rabies-Related Knowledge and Practices Among Persons At Risk of Bat Exposures in Thailand. Plos Neglect. Trop. Dis. 5, e1054, doi:10.1371/journal.pntd.0001054 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0001054&link_type=DOI) 

67. Wacharapluesadee, S. et al. Group C Betacoronavirus in Bat Guano Fertilizer, Thailand. Emerging Infectious Disease journal 19, 1349, doi:10.3201/eid1908.130119 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1908.130119&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23880503&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

68. Suwannarong, K. et al. Risk factors for bat contact and consumption behaviors in Thailand; a quantitative study. BMC Public Health 20, 841, doi:10.1186/s12889-020-08968-z (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12889-020-08968-z&link_type=DOI) 

69. Valitutto, M. T. et al. Detection of novel coronaviruses in bats in Myanmar. PLOS ONE 15, e0230802, doi:10.1371/journal.pone.0230802 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0230802&link_type=DOI) 

70. Brooks, T. M. et al. Measuring Terrestrial Area of Habitat (AOH) and Its Utility for the IUCN Red List. Trends Ecol. Evol. 34, 977–986, doi:10.1016/j.tree.2019.06.009 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tree.2019.06.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31324345&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

71. Herkt, K. M. B., Skidmore, A. K. & Fahr, J. Macroecological conclusions based on IUCN expert maps: A call for caution. Global Ecology and Biogeography 26, 930–941, doi:[https://doi.org/10.1111/geb.12601](https://doi.org/10.1111/geb.12601) (2017).
    
    

72. Jung, M. et al. A global map of terrestrial habitat types. Scientific Data 7, 256, doi:10.1038/s41597-020-00599-8 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41597-020-00599-8&link_type=DOI) 

73. Jung, M. et al. (Zenodo, 2020).
    
    

74. Tan, C. W. et al. A SARS-CoV-2 surrogate virus neutralization test based on antibody-mediated blockage of ACE2–spike protein–protein interaction. Nature biotechnology 38, 1073–1078, doi:10.1038/s41587-020-0631-z (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41587-020-0631-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32704169&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F09%2F14%2F2021.09.09.21263359.atom) 

75. Olival, K. J. et al. Possibility for reverse zoonotic transmission of SARS-CoV-2 to free-ranging wildlife: A case study of bats. PLOS Pathogens 16, e1008758, doi:10.1371/journal.ppat.1008758 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1008758&link_type=DOI) 

76. Anthony, S. J. et al. Global patterns in coronavirus diversity. Virus Evolution 3, doi:10.1093/ve/vex012 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ve/vex012&link_type=DOI) 

77. Anthony, S. J. et al. Global patterns in coronavirus diversity. Virus Evolution 3, vex012–vex012, doi:10.1093/ve/vex012 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ve/vex012&link_type=DOI) 

78. Murakami, S. et al. Detection and Characterization of Bat Sarbecovirus Phylogenetically Related to SARS-CoV-2, Japan. Emerging Infectious Disease journal 26, 3025, doi:10.3201/eid2612.203386 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2612.203386&link_type=DOI) 

79. Zhang, L. et al. Multilocus phylogeny and species delimitation within the philippinensis group (Chiroptera: Rhinolophidae). Zoologica Scripta 47, 655–672, doi:[https://doi.org/10.1111/zsc.12308](https://doi.org/10.1111/zsc.12308) (2018).
    
    

80. Wilson, D. E. & Mittermeier, R. A. Handbook of the Mammals of the World. Vol. 9. Bats. (Lynx Edicions, 2019).
    
    

81. Srinivasulu, B. & Srinivasulu, C. In plain sight: Bacular and noseleaf morphology supports distinct specific status of Roundleaf Bats Hipposideros pomona Andersen, 1918 and Hipposideros gentilis Andersen, 1918 (Chiroptera: Hipposideridae). Journal of Threatened Taxa 10, 12018–12026, doi:10.11609/jott.4111.10.8.12018-12026 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.11609/jott.4111.10.8.12018-12026&link_type=DOI) 

82. Becker, D. J. et al. Predicting wildlife hosts of betacoronaviruses for SARS-CoV-2 sampling prioritization: a modeling study. bioRxiv, 2020.2005.2022.111344, doi:10.1101/2020.05.22.111344 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wNS4yMi4xMTEzNDR2NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA5LzE0LzIwMjEuMDkuMDkuMjEyNjMzNTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

83. Wardeh, M., Baylis, M. & Blagrove, M. S. C. Predicting mammalian hosts in which novel coronaviruses can be generated. Nat. Commun. 12, 780, doi:10.1038/s41467-021-21034-5 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-21034-5&link_type=DOI) 

84. IUCN. (2021).
    
    

85. raster: Geographic Data Analysis and Modeling. R package version 2.9-23. (2019).
    
    

86. Chamberlain, S. & Boettiger, C. (PeerJ Preprints, 2017).
    
    

87. rgbif: Interface to the Global Biodiversity Information Facility API. R package version 3.5.2 (2021).
    
    

88. Feng, X. et al. A checklist for maximizing reproducibility of ecological niche models. Nature Ecology & Evolution 3, 1382–1395, doi:10.1038/s41559-019-0972-5 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41559-019-0972-5&link_type=DOI) 

89. Zizka, A. et al. CoordinateCleaner: Standardized cleaning of occurrence records from biological collection databases. Methods Ecol. Evol. 10, 744–751, doi:[https://doi.org/10.1111/2041-210X.13152](https://doi.org/10.1111/2041-210X.13152) (2019).
    
    

90. WorldPop. (2018).
    
    

91. gdalUtils: Wrappers for the Geospatial Data Abstraction Library (GDAL) Utilities. R package version 2.0.3.2. (2020).
    
    

92. lhs: Latin Hypercube Samples. R package version 1.0.1. (2019).