key: cord-0429216-19xp903w authors: Guth, Sarah; Hanley, Kathryn; Althouse, Benjamin M.; Boots, Mike title: Ecological processes underlying the emergence of novel enzootic cycles—arboviruses in the neotropics as a case study date: 2020-04-27 journal: bioRxiv DOI: 10.1101/2020.04.24.057430 sha: f74d4ea886475d7d9bb04c582d6fca6d3d9dba26 doc_id: 429216 cord_uid: 19xp903w Pathogens originating from wildlife (zoonoses) pose a significant public health burden, comprising the majority of emerging infectious diseases. Efforts to control and prevent zoonotic disease have traditionally focused on animal-to-human transmission, or “spillover”. However, in the modern era, increasing international mobility and commerce facilitate the spread of infected humans, non-human animals (hereafter animals), and their products worldwide, thereby increasing the risk that zoonoses will be introduced to new geographic areas. Imported zoonoses can potentially ‘spill back’ to infect local wildlife—a danger magnified by urbanization and other anthropogenic pressures that increase contacts between human and wildlife populations. In this way, humans can function as vectors, dispersing zoonoses from their ancestral enzootic systems to establish reservoirs elsewhere in novel animal host populations. Once established, these enzootic cycles are largely unassailable by standard control measures and have the potential to feed human epidemics. Understanding when and why translocated zoonoses establish novel enzootic cycles requires disentangling ecologically complex and stochastic interactions between the zoonosis, the human population, and the natural ecosystem. We address this challenge by delineating potential ecological mechanisms affecting each stage of enzootic establishment—wildlife exposure, enzootic infection, and persistence—applying existing ecological concepts from epidemiology, invasion biology, and population ecology. We ground our study in the neotropics, where four arthropod-borne viruses (arboviruses) of zoonotic origin—yellow fever, dengue, chikungunya, and Zika viruses—have separately been introduced into the human population. This paper is a step towards developing a framework for predicting and preventing novel enzootic cycles in the face of zoonotic translocations. Humans have frequently enabled pathogens to overcome physical barriers to dispersal (1) . 26 The European conquest of the Americas brought Old World diseases to the New World, 27 movement of troops during World War II propagated dengue viruses across the Asia-Pacific 28 region (2), and air travel has provided an international transmission network for emerging 29 infectious diseases (EIDs) such as the 2019 (ongoing) SARS-CoV-2 pandemic (3), the 2002-30 2003 SARS-CoV-1 outbreak (4) and pandemic influenza (5). Today, the majority of 31 pathogens that infect humans are broadly distributed across geographic regions-globalized 32 by human movement and population expansion, particularly during the past century (1) . 33 Animal pathogens have likewise spread globally through anthropogenic channels. The 34 globalization of agriculture has expanded the geographic range of many livestock diseases 35 with major economic repercussions, which continue to disproportionately affect the 36 developing world (6). Domestic and wild animals translocated by humans have introduced 37 their pathogens to new ecosystems, threatening biodiversity conservation-an anthropogenic 38 impact termed "pathogen pollution" (7). In some cases, these invasive animal infections have 39 2 maintained transmission post-emergence in local wildlife, establishing persistent reservoirs 40 that subsequently reseed transmission and thwart control efforts in the original animal host 41 population. Examples include African Swine Fever virus in Eastern Europe, where a novel 42 enzootic cycle of the invasive livestock pathogen in wild boars has prevented disease 43 eradication (8,9); and rabies virus in Africa, where human-mediated dispersal of domestic 44 dogs established wild carnivore reservoirs that now contribute to rabies persistence in both 45 wildlife and human communities (10). 46 Clearly, the global spread of zoonoses poses a unique and critical threat to human 47 health. Novel enzootic cycles occur when zoonoses are introduced to new regions, infect 48 local wildlife (spillback), and persist in local animal host populations (enzootic 49 establishment). Figure 1 provides a diagram of these processes and Table 1 provides 50 definitions of all terms in this paper). Now, more than ever, global conditions are ideal for the 51 generation of novel enzootic cycles. In an increasingly connected world, international trade 52 and travel provide pathways for pathogen introductions, while the recent surge in the 53 emergence and re-emergence of animal pathogens has increased the number of zoonoses 54 poised to exploit those pathways (11). Human population expansion into natural habitat is 55 intensifying contact between humans and animals, creating more opportunities for imported 56 zoonoses to spill back into naïve wildlife populations (12). The probability that these 57 introduced infections persist in animal populations is increasing as human development 58 pushes wildlife into crowded habitat patches and climate change alters transmission 59 conditions (7). 60 In this era of globalization, zoonoses are increasingly being recognized as global 61 threats. The emergence of SARS-CoV-2 in Wuhan, China has since affected 210 countries 62 and territories, causing nearly 2 million cases and 125,000 deaths worldwide as of April 15 th 63 2020 (13). The pandemic has prompted an extraordinary global response-many countries 64 have imposed nationwide lockdowns and closed their borders, nonessential international 65 travel has largely been suspended (14) to build infrastructure for global health security (17); and the WHO has declared recent 70 outbreaks of Ebola, H1N1, and Zika as public health emergencies of international concern 71 (PHEICs) (18). Nevertheless, this dialogue on the globalization of infectious disease 72 continues to conflate zoonoses and human-specific pathogens, often overlooking what makes 73 the spread of zoonoses so uniquely dangerous-the potential for enzootic reservoirs to 74 establish in previously naïve regions. Some zoonoses such as Ebola virus and SARS-CoV-1 75 have remained within the human population after introductions to new regions. Conversely, 76 introductions of yellow fever virus (YFV) in South America, Yersinia pestis (plague) in the 77 Americas ( An emerging body of literature is beginning to discuss the risk that human-to-animal 88 transmission will seed persistent enzootic reservoirs (22-25), but overall, disease emergence 89 and spillover from wildlife continues to dominate the conversation on zoonotic transmission. 90 Our understanding of enzootic establishment is limited by the difficulty in quantifying a 91 process that is both highly stochastic and the product of interactions between multiple 92 systems. In this review, we provide a conceptual framework to begin disentangling this 93 ecological complexity, applying concepts from disease ecology, invasion biology, and 94 population ecology. Cross-species pathogen emergence has previously been compared to 95 species invasions (26) and population ecology used to understand post-introduction 96 persistence (27) in the context of zoonotic spillover in human populations and host shifts 97 within wildlife communities. Adapting this interdisciplinary theory on pathogen emergence 98 to enzootic establishment, we review potential ecological mechanisms affecting the 99 probability that translocated zoonoses emerge in novel enzootic cycles. We discuss the 100 impact of each mechanism on the process of enzootic establishment: (1) local wildlife 101 becomes exposed, (2) the zoonosis successfully infects the novel hosts, and (3) transmission 102 persists indefinitely (28,29). We ground our discussion in the neotropics, where four 103 arboviruses of zoonotic origin-YFV, DENV, CHIKV, and ZIKV-have separately been 104 introduced into the human population. We additionally discuss the utility of modeling 105 approaches, which we illustrate by building a simulation model for our neotropics case study. Our aim is to delineate the ecological processes that shape the outcome of zoonotic 107 translocations as a first step towards developing a framework for predicting and preventing 108 novel enzootic cycles, and therefore we do not discuss other key factors such as immunity, 109 phylogeny, and evolution. preventing the other three arboviruses, as well as any new introductions, from also 139 establishing persistent enzootic reservoirs will be critical for forestalling further human 140 morbidity and mortality from zoonotic transmission in the American tropics. However, the 141 majority of the work on these imported arboviruses has focused on retrospective analysis of 142 the conditions that enabled their introductions-particularly the global invasion of mosquito 143 vectors Ae. aegypti and Ae. albopictus-and the transmission burden in the human population 144 (38-41). Only a few recent papers have begun to discuss the threat of enzootic establishment 145 (22-25,42). We add to this discussion by considering how each ecological mechanism 146 identified in our review may affect the trajectory of DENV, CHIKV, and ZIKV in South 147 America, using YFV as a frame of reference. 148 149 Wildlife exposure 150 Once a translocated zoonosis has established in a new human system, there is an immediate 151 risk that the pathogen spills back into local wildlife populations. The probability of spillback 152 first depends on the rate at which wildlife is exposed, which can be captured by propagule 153 pressure-a concept from invasion biology that represents the number, and temporal and 154 spatial distribution of nonnative individuals introduced to a new system, and a key 155 determinant of invasion success (43,44). 156 Propagule pressure hinges on introduction pathways between a source and recipient 157 population. The propagule pressure of a translocated zoonosis on local wildlife will vary 158 based on the availability of transmission pathways between the human (source) and wildlife 159 (recipient) populations. Borrowing from landscape and movement ecology, Borremans et al. (45) identified permeability-the likelihood that source and recipient hosts, along with the 161 pathogen, enter an ecosystem boundary region-as the ecological basis of pathways available 162 for cross-species pathogen emergence across ecosystem boundaries. With respect to the 163 human-wildlife boundary, permeability for translocated zoonoses will increase with wildlife 164 tolerance of (or preference for) anthropogenically modified landscapes, and human 165 communities' proximity to the edge of a species' habitat and incursions into natural habitat 166 for resource extraction. 167 Host boundary permeability creates opportunities for contacts between infected 168 source and recipient hosts, increasing propagule pressure on the recipient host population. 169 Transmission route determines the type(s) of contact, and thus degree of permeability, 170 required for the translocated zoonosis to cross the human-wildlife interface (45,46). Zoonotic 171 introductions can occur via direct contacts such as bushmeat hunting or the wildlife trade-172 common in the developing world-or via indirect mechanisms such as environmental 173 contamination-e.g., if infected bats leave saliva on forest fruits consumed by humans or 174 shed excreta in the human environment (61,62,66). For directly transmitted zoonoses, 175 propagule pressure on local enzootic systems will require sufficient permeability for humans 176 and wildlife to come into close contact. Conversely, zoonoses that can survive outside their 177 hosts will be less constrained by host boundary permeability. Domestic animals often 178 intersect with both humans and wildlife and thus, have the potential to bridge transmission 179 between the two populations. Vectors can likewise function as bridge hosts; as long as a 180 competent vector is present, vector-borne zoonoses only require some degree of spatial and 181 temporal overlap to transmit between humans and wildlife. 182 Transmission burden in the human population inhabiting the new region-a 183 combination of time since introduction and the number of subsequent human cases-will 184 5 determine the volume of pathogen propagules available to exploit transmission pathways to 185 local wildlife (48). Consistent circulation and a high number of cases in the human 186 population may result in more opportunities to spill back into enzootic transmission cycles, 187 producing greater propagule pressure. The precise propagule pressure that led YFV to invade 188 New World non-human primate populations is unknown, but phylogenetic analyses suggest 189 multiple spillback introductions (49,50). Evidence that YFV reached its current widespread 190 distribution in South America through long-distance, human-mediated dispersal implies that 191 many spillback introductions occurred across a broad spatio-temporal landscape. enzootic establishment, the translocated zoonosis must be able to infect the exposed animals. Infecting novel host species in a novel environment can be described as a niche shift. The 231 6 concept of an ecological niche has many nuanced definitions in ecology, but generally 232 represents the set of abiotic and biotic conditions that allow a species to occupy a particular 233 space within an ecosystem (62). A pathogen niche is defined by its hosts, vectors, 234 ecophysiological requirements, and the many ways in which these parts interact (63,64). Like 235 other species, a translocated zoonosis will have a realized niche-existing transmission 236 cycles-and a fundamental niche-the range of systems the zoonosis could theoretically 237 invade if given the opportunity (63). The probability that the zoonosis undergoes a niche shift 238 to infect novel hosts in its new range depends on the degree of overlap between its 239 fundamental niche and novel environment (62). 240 Host specificity influences the breadth of the pathogen's fundamental niche. 241 Generalists are defined by broad fundamental host ranges (65), which intuitively, will 242 intersect with a wider range of enzootic systems, facilitating shifts to novel hosts (66-68). 243 Alternatively, zoonoses can evolve to invade environments outside of their original 244 fundamental niches (62). Some pathogen types inherently have higher potential for 245 fundamental niche shifts than others. In particular, YFV, DENV, CHIKV, and ZIKV are all 246 single-stranded RNA viruses-a group of pathogens previously shown to be the most likely 247 to shift host species, predisposed to cross-species emergence by high mutation rates (28 The outcome of spillback events depends on the potential for transmission between 306 individuals in the novel animal host population. If animal transmission is limited, spillback 307 might result in an isolated wildlife case, or alternatively, trigger an outbreak that threatens 308 conservation efforts but eventually dies out (94,95). However, above a critical transmission 309 threshold, the translocated zoonosis will persist indefinitely-the final step of successful 310 enzootic establishment. In disease ecology and epidemiology, that critical transmission 311 threshold is represented by the effective reproduction number (R)-the average number of 312 secondary cases generated by a single infected individual in a population of susceptible and 313 non-susceptible hosts. While R > 1, each infected individual will, on average, produce at least 314 one secondary infection, allowing the pathogen to persist. The effective reproduction number 315 is a function of pathogen attributes such as transmissibility and duration of infectiousness, as 316 well as the average rate of contact between infected and susceptible individuals (and vectors 317 if vector-borne) in a given population. 318 The average contact rate between infected and susceptible hosts is largely shaped by 319 population ecology (96). Dynamic demographic rates and structuring of host populations 320 determine the abundance, distribution, and movement of hosts available to translocated 321 zoonoses in enzootic systems. Wildlife populations with high carrying capacities will have 322 larger baseline populations of susceptible hosts. High birth rates replenish the supply of 323 susceptible hosts, inhibiting mortality and conferred immunity from depleting the susceptible 324 population and driving the effective reproduction below the threshold of persistence (97). Spatial structuring of susceptible hosts can either limit or enhance the potential for pathogen 326 persistence. Pathogen dispersal between patches can promote persistence across a 327 metapopulation by buffering against local depletion of susceptible individuals. On the other 328 8 hand, spatial structure can limit contacts between patches, lowering the effective reproduction 329 number and driving an epidemic to extinction (96,98). 330 While some New World monkeys are large species with low birth rates, one of the 331 clade's defining features is its diverse set of remarkably small-bodied, short-lived primates. 332 For example, the pygmy marmoset (Callithrix pygmaea) is a small-bodied New World 333 primate with bimodal annual birth peaks, high twin birth rates, and a natural lifespan of about 334 10 years (99). These smaller-bodied primates also tend to require smaller home ranges and 335 live in higher densities, allowing for greater sympatric species richness. A meta-analysis of 336 New World primate assemblages found that on average, forest sites contained six sympatric 337 species, but this number could reach as high as 14 species, peaking near the equator (100). 338 New World enzootic mosquito species occupy the canopy, preying primarily on primates. 339 These vectors do not appear to demonstrate strong host preferences and thus, likely bridge 340 transmission between sympatric groups of primate species. This spatial structure-vector-341 mediated pathogen dispersal between primate groups (intra-and inter-species)-may have 342 facilitated the enzootic establishment of YFV. It has been hypothesized that enzootic YFV 343 persists within primate metapopulations occupying continuous forest in wandering 344 epizootics, in which transmission continually shifts between subpopulations. We suspect that 345 the demographic rates and structuring of primate host and mosquito vector populations in the 346 New World would, as with YFV, facilitate the persistence of DENV, CHIKV, and ZIKV 347 should they spill back successfully. The ecology of enzootic establishment in the Anthropocene previously non-endemic regions, significantly expanding its range in the New World. 374 375 9 A case for modeling 376 It is challenging to approach the ecological complexity of enzootic establishment through 377 field and experimental studies alone. Mathematical models can however be a useful tool that 378 allows us to integrate existing data and ecological theory to elucidate system dynamics, 379 particularly when data are sparse, as is often the case with enzootic systems (22,104,105) . 380 Many modeling approaches can be applied to a range of questions that allow us to better 381 understand the risk that translocated zoonoses will emerge in novel enzootic cycles. For 382 example, species distribution models, or ecological niche models-statistical approaches that 383 leverage associations between presence-absence information and environmental variables to 384 infer habitat suitability-can be used to identify at-risk systems where high permeability of 385 humans, vectors, and translocated zoonoses create ideal conditions for enzootic establishment 386 (46,106-108). Metapopulation modeling is often an important tool to demonstrate conditions 387 of pathogen persistence in wildlife populations given that they are typically fragmented (96). 388 Next-generation matrix methods are a useful tool to quantify the effective reproduction 389 number of translocated zoonoses in novel enzootic systems to understand the probability of 390 successful invasion (109). Explicit simulations that capture significant amounts of the 391 complexity of systems have been effectively used to compare the impact of interventions in 392 the human and wildlife population (110). We would also argue that simple models can be 393 very useful in helping us understand the system-specific dynamics that influence the 394 ecological processes underlying enzootic establishment. increasing EIP to 10 days, as we might observe with ZIKV, decreased the probability of 433 persistence, whereas shortening EIP to 2 days, as we might observe with CHIKV or at warm 434 temperatures, increased the probability of persistence (Figure 2) The goal of this study is to stimulate research on the emergence of novel enzootic cycles and 452 begin to disentangle the underlying ecology complexity. We have argued here that the 453 establishment of novel enzootic cycles is a pressing threat with the capacity to dramatically 454 alter disease dynamics. The International Task Force for Disease Eradication identifies the 455 existence of an animal reservoir as a barrier to eradicating a disease because enzootic 456 transmission often feeds human epidemics (119). In some cases, enzootic cycles have even 457 contributed to the evolution of pandemic pathogens; for example, pigs have functioned as 458 "mixing vessels" for the evolution of pandemic swine influenza (120). We delineated 459 potential ecological mechanisms at each stage of enzootic establishment, grounding our 460 discussion in the neotropics, where the danger of enzootic establishment is evident in the 461 history of YFV and an ongoing threat given the endemic circulation of DENV, CHIKV, and 462 ZIKV in the human population. Enzootic YFV, which has triggered devastating human 463 epidemics across the neotropics, has expanded its geographic range since its initial 464 establishment (121). There is a real danger that DENV, CHIKV, and ZIKV will also establish 465 persistent enzootic reservoirs in the New World, similarly inhibiting efforts to prevent future 466 human epidemics. Given that enzootic cycles are nearly impossible to control or eradicate, 467 avoiding enzootic establishment will be critical to mitigate the current arbovirus public health 468 emergency in the New World. Moreover, enzootic cycles can thwart efforts to eradicate 469 pathogens in the human population, as has occurred with the carriage of Guinea worm by 470 dogs (122) . In particular, the recent discovery of natural ZIKV (34-36) infection in non-471 human primates in Brazil calls for renewed urgency to understand the potential for enzootic 472 persistence. 473 Spillback events in which humans introduce pathogens into wildlife populations may 474 become common occurrences, as urbanization and anthropogenic pressure on wildlife 475 populations increase opportunities for human-wildlife contact. As a result, understanding the 476 risk that arboviruses will persist in sylvatic cycles after spillback events should be established 477 as a research priority, particularly in the era of SARS-CoV-2. There is significant concern 478 that SARS-CoV-2 could spill back into susceptible wildlife within its expanded geographic 479 range and establish novel enzootic reservoirs, becoming endemic outside of China. The high 480 burden and global distribution of human SARS-CoV-2 transmission places propagule 481 pressure on a wide range of enzootic systems. Bats-the putative origin of SARS-CoV-2 482 (123)-are the second most diverse mammalian group, inhabit every continent except 483 Antarctica, and harbor a large diversity of coronaviruses (124). Many bat species migrate 484 long distances (>1,000 km) (125), increasing their probability of exposure to SARS-CoV-2 485 and capacity to subsequently spread the virus. Although previous work has demonstrated that 486 human-adapted CoV often cannot infect bat cells (126), the cell-surface receptor implicated 487 in SARS-CoV-2 invasion of host cells-angiotensin converting enzyme 2 (ACE2)-has been 488 found to be highly diverse across bat species and in some cases, conducive to human SARS-489 CoV-2 infection (127). ACE2 also appears to be conserved across mammal species, 490 suggesting that the virus may have the potential to establish in a broad range of other 491 mammalian host species (128-130). Furthermore, there is a possibility that human SARS-492 CoV-2 could infect initially unsuitable host species through mutation events (131). 493 To mitigate the risk that SARS-CoV-2 establishes novel enzootic reservoirs and 494 becomes endemic outside of China, it is critical to identify susceptible wildlife species and 495 populations and implement policies that limit their exposure to the virus. The US government 496 has, at present, suspended all bat research to prevent humans from infecting and seeding an 497 enzootic reservoir of SARS-CoV-2 in North American bats (132). However, bat research is 498 still active in many other parts of the world, as are other forms of bat-human contact, though 499 many countries have banned bat bushmeat in response to the pandemic (133). Additional 500 policies may be needed to minimize human contact with other potentially susceptible wildlife 501 populations-species with ACE2 receptors predicted to associate with the SARS-CoV-2 502 spike receptor binding site (RBD), and populations that occupy anthropogenic landscapes 503 (particularly in regions experiencing a high burden of human transmission) and characterized 504 by demographic rates and structuring predicted to facilitate sustained epidemics and enzootic 505 persistence. It is also critical to monitor via surveillance whether these potentially susceptible 506 animals become exposed and begin to demonstrate capacity for between-host transmission. In 507 particular, researchers should monitor domestic animals, which can bridge transmission 508 between human and wildlife populations (45). Notably, efficient SARS-CoV-2 replication 509 and confirmed cases have been reported in cats (129,134) and ACE2 receptors of domestic 510 cattle have been predicted to associate with the SARS-CoV-2 RBD (135). 511 The risk of enzootic persistence depends on a multitude of factors, many of which are 512 unknown or poorly characterized, such as infectivity, number of spillback events, and the 513 transmission conditions needed for persistence. Here we provide a conceptual framework of 514 ecological factors to begin addressing the challenge of predicting this risk. Considering 515 ecological mechanisms is a first step towards developing targeted intervention strategies. We 516 additionally argue for the utility of modeling in detangling ecological complexity, providing a 517 simulation model of arboviral transmission in New World primates and mosquitoes as an 518 example. Our results indicate that although long extrinsic incubation periods (EIPs) can 519 reduce the probability of enzootic persistence, the long mosquito lifespans that are 520 characteristic of tropical New World sylvatic species may negate this effect-suggesting that 521 differences in EIP that we may have expected to be important in determining the translocated 522 arboviruses' relative risk of enzootic establishment are unlikely to be relevant. Overall, our 523 work is important in highlighting the need to be vigilant of imported zoonoses and 524 emphasizing the importance of robust programs to mitigate the risk of spillback events that 525 lead to enzootic persistence. 526 Table 1 . Definitions of terms used in this paper. 527 528 Table 2 populations (enzootic establishment). 553 • Understanding when and why translocated zoonoses establish novel enzootic cycles 554 requires disentangling ecologically complex and stochastic interactions between the 555 zoonosis, the human population, and the natural ecosystem. 556 • Permeability of the human-wildlife interface for humans, wildlife, and the translocated zoonosis creates pathways for propagule pressure. Transmission route determines the type(s) of contact, and thus degree of permeability, required for a translocated zoonosis to cross the humanwildlife interface. Habitat specificity: generalists are more tolerant of human-modified landscapes, increasing wildlife permeability. Natural resource extraction and land conversion increases human permeability. Several New World primate hostsincluding Aotus spp., the most susceptible to YFV-are habitat generalists that often occupy the forest edge. Deforestation, urbanization, and bushmeat hunting encroach on wildlife habitat in South America. Vector-borne transmission may have initially constrained YFV, but now, abundant mosquitoes likely increase human-mediated dispersal and may facilitate wildlife exposure to urban transmission without direct contact with humans. Transmission burden in the human population determines the volume of pathogen propagules available to exploit pathways to enzootic systems. Heightened arboviral burden in the human population likely increases propagule pressure of DENV, ZIKV, and CHIKV on New World enzootic systems. Enzootic transmission of a translocated zoonosis persists indefinitely while R > 1. Demographic rates such as the carrying capacity and birth rate determine the abundance of susceptible hosts. a translocated zoonoses infects exposed animals, undergoing a niche shift. Species diversity creates the host environment available to translocated zoonoses, and facilitates or inhibits invasion depending on the context. When additional species are susceptible and infectious to the translocated zoonosis, increasing species diversity is expected to amplify the probability of enzootic infection. However, when additional species are not competent hosts, there can be a dilution effect. Competitive exclusion can constrain host environments created by biodiversity. Native zoonoses already circulating in an enzootic system can compete for susceptible hosts in an enzootic system, preventing infections of translocated zoonoses. New World primates and mosquitoes are remarkably diverse. Given evidence that these species are competent arboviral hosts, this biodiversity could amplify the probability of enzootic infections of YFV, DENV, CHIKV, and ZIKV. Conserved cellular receptors and transmission dynamics suggest that YFV, DENV, CHIKV, and ZIKV are generalists. The high mutation rates of single-stranded RNA viruses may additionally trigger fundamental niche shifts. There is some evidence of cross-protection between the arboviruses, but cocirculation in urban transmission cycles does not support this hypothesis. Overlap between the fundamental niche of the zoonosis and the novel environment determines the probability that the zoonosis undergoes a niche shift to infect novel hosts in its new range. Pathogen host specificity influences the breadth of the fundamental niche. Generalists are more likely to have the ability to infect and transmit between novel hosts. Host population ecology shapes the average contact rate between infected and susceptible hosts-a key determinant of the effective reproduction number. Population structure determines opportunities to contact susceptible hosts, and in metapopulations, can create pathogen dispersal patterns that either buffer against or promote local extinction. Vector-mediated pathogen dispersal between primate groups (intra-and interspecies) may facilitate the persistence of New World enzootic YFV-a persistence pattern that has been described as wandering epizootics. The New World has a diverse set of remarkably small-bodied, short-lived primates that live in relatively high densities. Our approach builds on Althouse et al. (1) , which explored how the force of infection and 3 primate birthrate affect the probability that ZIKV will establish a sylvatic cycle in the Americas. 4 Expanding their analysis, we explored the effects of the mosquito extrinsic incubation period 5 (EIP) and lifespan, and primate birthrate, and considered DENV, CHIKV, and YFV in addition 6 to ZIKV. We simulate the introduction of a single infected primate into metapopulations of 7 susceptible primate hosts and mosquito vectors, applying the Gillespie stochastic simulation 8 algorithm with the Binomial Tau Leap approximation to the following transition rates for 9 mosquitoes: 10 11 12 The equations divide primate and mosquito metapopulations into primate species 1,...i and 13 mosquito species 1,...j. The model assumes host preference, with mosquito species j biting the 14 corresponding primate species j more often than other primate species-characterized by the on-15 diagonal rates within a contact matrix. Primate-mosquito species pairs are coupled by off-16 diagonal cross-biting rates, calculated as a fixed fraction (10%) of the on-diagonal within-pair 17 biting rates. Given that vertical transmission in the mosquito population is thought to be 18 negligible in enzootic cycles (2), the mosquito transovarial transmission rate is set to zero. Thus, all primates and mosquitoes are born susceptible to arbovirus infection at rates and 20 infected at rates (t), proportional to the number of mosquito bites given or received per day and 21 the probability of successful transmission. We assume that birthrate = 1/lifespan-a conservative 22 estimate given that primates' peak reproductive years occur before the age of mortality (3) and 23 mosquitoes have the capacity to complete multiple gonotrophic cycles within a lifetime (4). 24 Thus, birthrates vary as we explore three mosquito lifespans (7, 14, and 21 days). Transmission 25 probabilities vary seasonally in Equation 14 and 15 due to changes in environmental conditions 26 such as rainfall and temperature (5). After infection, mosquitoes enter the exposed compartment, 27 progressing to the infectious period at one of three fixed rates ! ! (1/2, 1/7, or 1/10 days). 28 Globalization of human infectious 584 disease Urbanization and Globalization: The Unholy Trinity of the 21st 586 Century. Tropical Medicine and Health Pathogenesis, 588 and Control of COVID-19. Viruses Confronting the New Challenge in Travel Medicine: SARS. 590 Spread of a Novel Influenza A 592 (H1N1) Virus via Global Airline Transportation Current drivers and future directions of global livestock 595 disease dynamics Emerging infectious diseases 598 of wildlife: a critical perspective African swine fever in the North 600 Caucasus region and the Russian Federation in years 2007-2012. Virus Research African Swine Fever Status in Europe. Viruses Canine Rabies Elimination in Africa: Dispelling Doubts with Data. PLoS Neglected 606 Tropical Diseases Emerging zoonotic epidemics in the interconnected global community. 608 Veterinary Record Urbanization and Disease Emergence: 610 Dynamics at the Wildlife-Livestock-Human Interface COVID-19) Situation Report -85. World Health 613 Organization Coronavirus Travel Restrictions, Across the Globe. New 615 York Times Saudi Arabia Declares Cease-Fire in Yemen, Citing Fears of 617 Coronavirus WHO and UNICEF to partner on pandemic response through COVID-19 Solidarity 619 Response Fund. WHO News Release The evolution, etiology and eventualities of the global health security 621 regime. Health Policy and Planning Delays in Global Disease Outbreak Responses: Lessons from 623 H1N1, Ebola, and Zika Epidemiology of Human Plague in 625 the United States Kilpatrick AM. Globalization, Land Use, and the Invasion of West Nile Virus Yellow fever in the Americas. Infectious Disease Clinics of 629 North America Potential for Zika 631 Virus to Establish a Sylvatic Transmission Cycle in the Americas. PLOS Neglected Tropical 632 Diseases Urbanization and geographic expansion of zoonotic arboviral diseases: 634 mechanisms and potential strategies for prevention Confronting 637 data sparsity to identify potential sources of Zika virus spillover infection among 638 primates High risk for chikungunya virus to initiate an 640 enzootic sylvatic cycle in the tropical Americas. PLOS Neglected Tropical Diseases Parasites and biological invasions: parallels, interactions, and 643 control Population Dynamics of Pathogens with Multiple Host Species. The American 645 Naturalist Emerging pathogens: the epidemiology and 647 evolution of species jumps Does pathogen plasticity facilitate host shifts? PLOS Pathogens Chikungunya, and Other Emerging 651 Vector-Borne Viral Diseases Fever versus fever: 653 The role of host and vector susceptibility and interspecific competition in shaping the 654 current and future distributions of the sylvatic cycles of dengue virus and yellow fever 655 virus. Infection, Genetics and Evolution Mayaro Virus in Wild Mammals Emerging Infectious Diseases Fever from the forest: 659 prospects for the continued emergence of sylvatic dengue virus and its impact on public 660 health Evidence of 662 natural Zika virus infection in neotropical non-human primates in Brazil Zika 665 Virus in Peridomestic Neotropical Primates, Northeast Brazil. EcoHealth First detection 668 of Zika virus in neotropical primates in Brazil: a possible new reservoir Mosquitoes infected with dengue viruses in Brazil Past and 672 future spread of the arbovirus vectors Aedes aegypti and Aedes albopictus The arboviral burden of 675 disease caused by co-circulation and co-infection of dengue, chikungunya and Zika in 676 the Americas. Travel Medicine and Infectious Disease Climate 678 change influences on global distributions of dengue and chikungunya virus vectors 41. Fauci AS, Morens DM. Zika Virus in the Americas -Yet Another Arbovirus Threat. New 682 England Journal of Medicine Into the 684 woods: Changes in mosquito community composition and presence of key vectors at 685 increasing distances from the urban edge in urban forest parks in Manaus, Brazil. Acta 686 Tropica Species Diversity and Biological Invasions: Relating Local Process to 688 Community Pattern Simberloff D. The Role of Propagule Pressure in Biological Invasions. Annual Review of 690 Ecology, Evolution, and Systematics Cross-species pathogen 692 spillover across ecosystem boundaries: mechanisms and theory Pathways to 696 zoonotic spillover Ecological 698 dynamics of emerging bat virus spillover Dissecting the null model 701 for biological invasions: A meta-analysis of the propagule pressure effect Out of Africa: A Molecular Perspective on the 704 Introduction of Yellow Fever Virus into the Americas 706 Phylodynamics of Yellow Fever Virus in the Americas: new insights into the origin of the 707 2017 Brazilian outbreak Drivers, dynamics, and control of emerging vector-borne 709 zoonotic diseases. The Lancet Assessing habitat utilization by neotropical primates: A new approach. 711 Primates Having your wildlife and eating it too: an analysis of hunting 713 sustainability across tropical ecosystems World Resources Institute, Global Forest Watch (Organization) Human pressure on the Brazilian 716 Amazon forests Gubler DJ. The Global Emergence/Resurgence of Arboviral Diseases As Public Health 718 Problems. Archives of Medical Research The 720 History of Dengue Outbreaks in the Americas. The American Journal of Tropical 721 Medicine and Hygiene Pan American Health 727 Organization (PAHO), World Health Organization (WHO) Epidemiology of Chikungunya 729 in the Americas Pan American Health Organization. Zika Epidemiological Update 732 Current status, challenges and perspectives in the development of vaccines against 733 yellow fever, dengue, Zika and chikungunya viruses Unifying niche shift 735 studies: insights from biological invasions. Trends in Ecology & Evolution The niche 738 reduction approach: an opportunity for optimal control of infectious diseases in low-739 income countries? Vacated niches, competitive release and the community ecology of 741 pathogen eradication Concentration and Specialization: Dynamics of Niche Width in Populations of 744 Organizations Role of evolved host breadth in the initial 746 emergence of an RNA virus: predicting virus emergence Phylogenetic signal in plant pathogen-host range Host and viral 751 traits predict zoonotic spillover from mammals Biodiversity decreases disease 753 through predictable changes in host community competence Frontiers in research on biodiversity and disease Asymmetric competitive suppression between strains 758 of dengue virus Emergence potential of sylvatic dengue virus type 4 in the urban transmission cycle is 761 restrained by vaccination and homotypic immunity A competitive exclusion principle for pathogen virulence Flavivirus 765 Receptors: Diversity, Identity, and Cell Entry Cellular Attachment and Entry Factors for Chikungunya Virus Population biology of emerging and re-emerging pathogens. Trends in 769 Microbiology Evidence for multiple sylvatic transmission cycles during 772 the 2016-2017 yellow fever virus outbreak, Brazil Viral 775 kinetics of primary dengue virus infection in non-human primates: A systematic review 776 and individual pooled analysis 778 Phylogenetic structure and host abundance drive disease pressure in communities. 779 Nature Phylogenetic host specificity and 781 understanding parasite sharing in primates Centrality in primate-parasite networks reveals the 784 potential for the transmission of emerging infectious diseases to humans Host phylogenetic distance drives trends in virus 787 virulence and transmissibility across the animal-human interface Experimental Zika Virus Infection of Neotropical Primates. The American Journal of 790 Tropical Medicine and Hygiene Zika Virus Inoculation in a New World Monkey Model Reproduces Key Features of the 793 Human Infection. Scientific Reports Viral 795 kinetics of primary dengue virus infection in non-human primates: A systematic review 796 and individual pooled analysis An animal model that reflects human disease: the common 798 marmoset (Callithrix jacchus). Current Opinion in Virology Zika Virus Relative to Aedes aegypti. Viruses Yellow fever cases in Asia: primed for an epidemic. 803 International Journal of Infectious Diseases Evaluation of 805 serological cross-reactivity between yellow fever and other flaviviruses. International 806 Journal of Infectious Diseases Does 808 immunity after Zika virus infection cross-protect against dengue? The Lancet Global 809 Health Co-circulation of 811 dengue, chikungunya, and Zika viruses in Colombia from Can Zika virus antibodies 814 cross-protect against dengue virus? The Lancet Global Health Chikungunya virus immunity on Mayaro virus disease and epidemic potential Measles in monkeys: an epidemiological 819 study Context-821 dependent conservation responses to emerging wildlife diseases. Frontiers in Ecology 822 and the Environment Journal of Animal Metapopulation moments: coupling, stochasticity and 824 Ecology 2000, 69, 725±736 persistence The invasion, persistence and spread of infectious diseases 826 within animal and plant communities Metapopulation dynamics of bubonic plague Evolutionary 831 genetics and implications of small size and twinning in callitrichine primates Species coexistence, distribution, and environmental 834 determinants of neotropical primate richness: A community-level zoogeographic 835 analysis. Primate Communities Yellow Fever Remains a Potential Threat to Public 837 Health. Vector-Borne and Zoonotic Diseases Climate change and vector-borne 839 diseases: a regional analysis Habitat 841 connectivity affects specialist species richness more than generalists in veteran trees. 842 Forest Ecology and Management The tortoise or the hare? Impacts of 844 within-host dynamics on transmission success of arthropod-borne viruses Sylvatic Dengue Isolations: A Multi-Host, Multi-Vector SIR Model of Dengue Virus 849 Transmission in Senegal. PLoS Neglected Tropical Diseases Ecological Contexts of Index Cases and 851 Spillover Events of Different Ebolaviruses The impact of human population pressure on 853 flying fox niches and the potential consequences for Hendra virus spillover Mapping the 856 zoonotic niche of Marburg virus disease in Africa. Transactions of The Royal Society of 857 Tropical Medicine and Hygiene Invasion and 859 Persistence of Infectious Agents in Fragmented Host Populations Ecological 862 interventions to prevent and manage zoonotic pathogen spillover Aedes Mosquitoes: Saliva Is Infectious as soon as Two Days after Oral Infection. PLoS 867 ONE Incubation Periods of 869 Yellow Fever Virus. The American Journal of Tropical Medicine and Hygiene The Incubation Periods of Dengue Viruses. Vasilakis N, 872 editor Estimated Incubation Period for Zika Virus 874 Disease. Emerging Infectious Diseases Effect of variability in infection period on the persistence 876 and spatial spread of infectious diseases Dispersal and Survival of Male and Female 879 Aedes albopictus (Diptera: Culicidae) on Réunion Island Release-Recapture Experiments with Canopy Mosquitoes in the Genera Haemagogus 883 and Sabeihes (Diptera: Culicidae) in Brazilian Amazonia 886 et al. Yellow Fever Virus in Haemagogus leucocelaenus and Aedes serratus Mosquitoes Emerging Infectious Diseases International Task Force for Disease Eradication Reverse zoonosis of influenza to swine: new perspectives on 891 the human-animal interface Progress of Sylvan Yellow Fever Wave in Central America: Nicaragua and 893 American Journal of Public Health and the Nations Health The Lancet. Guinea worm disease eradication: a moving target. The Lancet A pneumonia outbreak 898 associated with a new coronavirus of probable bat origin Global Epidemiology of Bat Coronaviruses. Viruses Ecology of bat migration Receptor recognition and cross-species infections of SARS coronavirus Angiotensin-converting enzyme 2 (ACE2) 907 proteins of different bat species confer variable susceptibility to SARS-CoV entry. 908 Archives of Virology Comparative ACE2 variation 910 and primate COVID-19 risk Atlas of ACE2 gene expression in mammals reveals novel 913 insights in transmisson of SARS-Cov-2 [Internet]. Genomics Broad Host 916 Range of SARS-CoV-2 Predicted by Comparative and Structural Analysis of ACE2 in 917 Vertebrates A Genomic Perspective on the Origin and Emergence of SARS-920 CoV-2 advises suspending bat research over concerns coronavirus could 922 infect North American species Potential for Zika 106 Virus to Establish a Sylvatic Transmission Cycle in the Americas. PLOS Neglected Tropical 107 Diseases How important is vertical transmission in mosquitoes for the persistence 109 of dengue? Insights from a mathematical model Reproductive 111 aging patterns in primates reveal that humans are distinct Physiological bases of mosquito ecology The tortoise or the hare? Impacts of within-host 115 dynamics on transmission success of arthropod-borne viruses Effect of variability in infection period on the persistence and 118 spatial spread of infectious diseases Realistic Distributions of Infectious Periods in Epidemic Models: Changing 120 Patterns of Persistence and Dynamics. Theoretical Population Biology Appropriate Models for the Management of Infectious 123 Diseases Discrete stochastic analogs of Erlang epidemic models Holding parameters constant whenever possible allowed us to minimize model complexity and 33target the effects of our parameters of interest. Furthermore, for the majority of these parameters, 34there is not enough empirical data to accurately characterize differences between New World 35 species. Full lists of parameter values used for the simulations presented in Figure 2 of epidemics, and critical population size required for disease persistence in directly transmitted 50 pathogens (6-8). Given the short mosquito lifespan and the resulting interaction between EIP and 51 the length of the infectious period, epidemic dynamics would be more sensitive to the 52 distribution of waiting periods in the mosquito vector than in the primate host. Thus, we explored 53 the effect of modeling the latent period in the mosquito vector (EIP) as an Erlang distribution by 54 splitting the latent period into a series of separate compartments, commonly referred to as a 55 boxcar configuration (9). 56We constructed two models with Erlang-distributed EIPs-a model with 10 exposed 57compartments, and a model with 50 exposed compartments. Figure S1 demonstrates the 58 differences between the probability distributions of the exponentially-distributed EIP relative to 59both Erlang-distributed EIPs. As expected, modeling EIP as less-dispersed Erlang distributions 60 resulted in higher probabilities of longer EIPs, thus decreasing the probability of sylvatic 61persistence. However, the Erlang-distributed models produced the same general trends as the 62 exponentially distributed models, with shorter EIPs and longer mosquito lifespans increasing the 63 probability of sylvatic establishment. As a result, in the interest of minimizing model complexity, 64in our main text, we report results from our exponentially-distributed model. Simulation results 65 from the Erlang-distributed models with respect to differences in the mosquito EIP and lifespan 66 are presented in Figure S2 and S3. periods (such as the EIP) in traditional compartment models follow exponential probability 75 distributions (in blue). In the case of EIP, this assumption is unrealistic because individuals have 76 an equal probability of leaving the incubation compartment regardless of the time since infection. 77In reality, the probability of progressing should be very low immediately after infection and 78highest around the EIP mean-a trajectory better represented by a gamma distribution (in green 79 and purple). A gamma distributed EIP can be constructed in a compartmental model via a boxcar 80 configuration of the latent period (i.e., splitting the latent period up into a series of separate 81 compartments, or boxes). Note how as the number of latent boxes increases (green vs. purple), 82the gamma distribution becomes less dispersed, and more closely centered around the EIP mean. Mosquito birthrate = 1/lifespan and increases from the top to bottom panels, while EIP increases 97 left to right. Here, we constructed an Erlang-distributed EIP by splitting the exposed 98 compartment into 50 separate boxes. Within each panel, the total population size of mosquitoes 99 (in two populations) and primates (in two populations) changes horizontally and vertically, 100respectively. For each parameter set, we simulated the introduction of a single infected primate 101 and subsequent transmission for a three-year period. Blue indicates no simulations establishing, 102whereas red indicates all simulations establishing. Contour lines show 0.25, 0.5, 0.75, and 0.95 103 probability of establishment. 104