key: cord-0813469-9mw1v1nm
authors: González-Salazar, Constantino; Stephens, Christopher R.; Sánchez-Cordero, Víctor
title: Predicting the Potential Role of Non-human Hosts in Zika Virus Maintenance
date: 2017-02-08
journal: Ecohealth
DOI: 10.1007/s10393-017-1206-4
sha: beefb6193302ec5be2c1ae75e0f6247cf070be5d
doc_id: 813469
cord_uid: 9mw1v1nm

Arboviruses are often maintained in complex cycles involving vertebrates such as mammals or birds and blood-feeding mosquitoes. However, the role of wildlife hosts in their emergence or re-emergence in human populations has received little attention. The recent emergence of Zika virus in America, and previous occurrences of chikungunya and dengue, forces us to confront a potential new disease-emergence phenomenon. Using a spatial data mining framework to identify potential biotic interactions, based on the degree of co-occurrence between different species, we identified those mammal species with the highest potential for establishing mammal–vector interactions, considering as principal vector Aedes aegypti. Seven of the top ten identified mammal species with highest potential were bats, with two of them having previously been confirmed as positive hosts for dengue in Mexico. We hope that this will raise interest of Mexican public health authorities and academic institutions to assess the role of wild hosts in the maintenance and spread of arboviruses.

The current outbreak of Zika virus (ZIKV) is the most recent disease threatening human health worldwide (Lucey and Gostin 2016; WHO 2016) . ZIKV is a mosquito-borne virus (arbovirus) belonging to the family Flaviviridae, which was first isolated from a rhesus monkey in the Zika forest of Uganda in 1952 (Dick et al. 1952) . Until 2007, sporadic human disease cases had been reported from countries in Africa and Asia (Musso et al. 2014 ). However, last year the virus was detected in America (Zanluca et al. 2015) , where it has been spreading explosively (ECDC 2015; Petersen et al. 2016) . A great number of cases with neurological and autoimmune complications have been reported from 20 countries in the Americas, mainly in South America (Ventura et al. 2016; WHO 2016) , though a substantial number of human cases have recently been identified in North America (Chen and Hamer 2016; McCarthy 2016) .

Arboviruses are often maintained in complex cycles involving vertebrates and blood-feeding vectors (Blum et al. 2016) and are transmitted by the bite of infected female mosquitoes (Lefèvre et al. 2009 ). Consequently, humans could be at high risk of infections from pathogens circulating within the blood system of other wildlife vertebrate species (Kaddumukasa et al. 2015) . Due to their current relevance to public health, the identification of mosquitoes' wildlife blood-meal sources and their relative importance is important in order to determine the linkages between potential hosts and vectors. The arboviruses ZIKV, dengue (DENV) and chikungunya (CHIKV) have been isolated from several Aedes mosquito species (Klimpel and Mehlhorn 2014) , most notably Aedes aegypti and Aedes albopictus, which are widespread worldwide (Grard et al. 2014; Kraemer et al. 2015) . Although other mosquito species have been identified as having links with these viruses, A. aegypti has generally been considered to be the most relevant in Mexico. However, the role of potential wildlife hosts in maintenance of these arboviruses has received relatively little attention, even though there are results that indicate that they are not necessarily restricted to epizootic cycles in the New World. For instance, although DENV is widely believed to be absent from New World wildlife, all four serotypes of DENV have been found in a wide variety of South American mammals in French Guiana (de Thoisy et al. 2009 ), including bats, rodents and marsupials. Notably, identified hosts were encountered in both DENV endemic areas as well as areas where the disease was nearly absent. Genetic evidence was provided that was consistent with the hypotheses that (1) there exist spillover infections from humans to mammals and (2) there exists an enzootic cycle. In Mexico, six bat species have been identified as positives for DENV (Artibeus jamaicensis, Artibeus lituratus, Carollia brevicauda, Myotis nigricans, Pteronotus parnellii and Natalus stramineus) (Aguilar-Setién et al. 2008; Machain-Williams et al. 2013; Sotomayor-Bonilla et al. 2014) .

In contrast, there are no identified hosts of CHIKV or ZIKV in the New World, although both have multiple mammal hosts in Africa and Asia. In particular, for bats, CHIKV has been isolated in three Chinese bats (Rousettus aegyptiacus, Hipposideros caffer, Chaerephon pumilus) and one from Senegal (Scotophilus sp.). For ZIKV, no natural infections of bats have been documented; however, the bat Myotis lucifugus was infected successfully in laboratory conditions (Reagan et al. 1955) .

Although the question remains as to whether wildlife and/or domestic hosts can maintain DENV, CHIKV and ZIKV in enzootic cycles and, therefore, play a role in its reemergence in human populations, or that infections in human populations can spillover into wildlife hosts, these preliminary results highlight the need to conduct studies focused on identifying and understanding the role of wildlife species in the spread and maintenance of these arboviruses. Thus, we urgently need research on these viruses and the ecologic, entomologic and host determinants of viral maintenance and emergence (Blum et al. 2016) . A further complication comes from the complex evolution of arboviruses (Weaver and Barrett 2004; Turner et al. 2010 ) which, as RNA viruses, are highly mutable, thus allowing for potentially rapid adaptive changes to new hosts and changes in host range.

Additionally, the steady increase in contact between human and wildlife by the invasion and destruction of natural forest, as well as climate and environmental changes, is allowing mosquitoes to be exposed to potentially new wildlife species hosts and pathogens, thus being able to potentially create new disease patterns (Kaddumukasa et al. 2015) . Therefore, a major short-term goal should be to conduct specific, directed studies of potential wildlife hosts involved in arbovirus maintenance. Although host selection by mosquitoes depends on a great number of factors, a first step should be to determine which wildlife species have a high risk of being in contact with the disease vectors in order to optimize resource allocation for these studies. The increased availability of spatial data and recent methodological developments in species distribution modelling (Stephens et al. 2009; González-Salazar et al. 2013 ) allow us to implement spatial analysis that can be used to build predictive models for the presence of emerging diseases. In particular, they allow for the direct incorporation of biotic factors, which are not explicitly accounted for in most niche modelling studies, including inter-specific interactions, as well as abiotic factors. Biotic factors, such as host range, are obviously of crucial importance in the dynamics of zoonoses. The method uses point collection data for potential/known vectors and potential/known hosts and constructs a diagnostic, epsilon (Stephens et al. 2009 ), that measures the statistical significance of the degree of co-occurrence of a pair of species (e.g. a potential vector and a potential host) relative to the null hypothesis that they are uncorrelated. In the present case, the mammal data set used contains 37,297 unique point collections from geo-referenced localities for 427 terrestrial mammals occurring in Mexico (GBIF; www.gbif. org, and CONABIO; www.conabio.gob.mx). For A. aegypti, there were 302 collection points taken from the Sistema Nacional de Informació n sobre la Biodiversidad of the Comisió n Nacional para Uso y Conservació n de la Biodiversidad (www.conabio.gob.mx). Co-occurrences were defined on a uniform grid covering Mexico of 3337 rectangular cells of size 25 km 9 25 km. With these data, a predictive model for the potential distribution of A. aegypti in Mexico was created based only on biotic factors. This approach gives us two principal results: (1) a ranked list of mammals that can potentially be blood-meal sources for the mosquito and therefore potential hosts and (2) a map of potential presence of A. aegypti linked to the presence of potential mammal hosts (Stephens et al. 2009 ).

Of 427 potential vector-host pair interactions, we show in Table 1 the top 25% (108) of most statistically significant positive associations based on our co-occurrence index, epsilon. This ranked mammal list is a first-pass predictive model for the most important potential bloodmeal sources for A. aegypti; consequently, they are, at this level of description, the main candidates to be positive for DENV, ZIKV or CHIKV. Of course, such a statistical association does not necessarily prove that there is a direct ''causal'' interaction between these taxa. Neither does it allow for a more detailed description of the potential host role, i.e., as a competent host, dead-end host, alternative host. However, previous studies using and validating this method led to the prediction and subsequent discovery of 22 previously unknown mammals hosts for Leishmaniasis (Stephens et al. 2016) . Additionally, evaluating our mammal list against the known results for DENV in Mexico, five of the six bat species positive to DENV (Aguilar-Setién et al. 2008 ) are in the top 50 highest ranked species (P < 10 -4 ). It is notable that 7 of the top 10 species are bats, with two being previously identified as positive for the presence of dengue, and, interestingly, four of them (Glossophaga soricina, A. jamaicensis, A. lituratus and Sturnira lilium) are identified hosts of Leishmania (L.) mexicana and Trypanosoma cruzi (Villegas-García and Santillán-Alarcó n 2001; Berzunza-Cruz et al. 2015; Ló pez-Cancino et al. 2015) . Within the highest ranked species are all major groups of mammals, with many of them being synanthropic species (Table 1) (Estrada et al. 1994; Medellín et al. 2000; Ruiz-Piña and Cruz-Reyes 2002; Martínez-Hernández et al. 2014; Hennessy et al. 2015; Ló pez-Cancino et al. 2015; Cruz-Salazar et al. 2016) . Additionally, the genus Artibeus, Sturnira and Carollia have all been identified as hosts for other pathogens in fragmented landscapes, including coronavirus, pegivirus, hepacivirus, rabies and leptospira Quan et al. 2013; Chávez et al. 2015) .

To identify those Mexican states with a high probability of significant mosquito-mammal interaction, we modelled the potential distribution of A. aegypti based on mammal point collection data using a score function (see Stephens et al. 2009) (Fig. 1a) . This map shows areas where vector and potential mammal hosts are most likely to both be present and, therefore, potentially interact. As stated, an important potential use of this information is to target further investigation by concentrating field studies and/or surveillance in those areas. Note that, our list is not restricted to bat species, including several other groups of mammals, for instance rodents, marsupials, carnivores, etc., indicating that field studies should have the capacity to collect a wide, representative set of potential mammal hosts.

To test model accuracy as a function of score, we grouped Mexican states by score quartiles. The 4th quartile corresponds to the 25% of states with the highest score values, the 3rd quartile to the next 25% of states with highest score values, etc. This allows us to establish predictability profiles across the different score quartiles for our biotic model. We randomly selected 30% of mosquito records (92 points) to calculate for each score quartile the percentage of associated A. aegypti records. Forty percentages of points were observed in the top quartile, and this percentage significantly decreased towards the lower quartiles (X 2 = 38.64, P < 0.00001). Taking as a null hypothesis that A. aegypti are distributed randomly in Mexico, the outcomes showed that biotic model predicted accurately the most important potential areas for A. aegypti presence.

Considering the correlation between current human cases of ZIKV in Mexico with the predicted distribution by state of A. aegypti using our biotic model we see that the vast majority of autochthonous cases have occurred in states with high values of scores (i.e. high probability to find significant mosquito-mammal interactions). For instance, Chiapas and Guerrero reported 500 and 558 cases, respectively, corresponding to 56% of the 1884 cases confirmed at present (http://www.epidemiologia.salud.gob.mx/). We believe there is a significant risk that synanthropic mammals, such as bats, rodents or marsupials, could be infected thus favouring virus maintenance. We have restricted our analysis, conservatively, to A. aegypti, as vector as it has been identified positively for all three viruses-DENV, CHIKV and ZIKV. However, our methodology can be applied to consider the relation between potential mammal hosts and any other potential vector or combination thereof.

In distinction to important zoonosis like Leishmaniasis and Chagas, for which the role of wild hosts is better understood (De Almeida Curi et al. 2006; Herrera 2010) , for DENV, CHIKV and ZIKV and other arbovirus, our understanding of the role of non-human hosts, either as a spillover from an epizootic cycle or via the potential existence of enzootic cycles, is very poor. However, given the potentially important epidemiological consequences of either of these scenarios, it is important to explore these possibilities. Our results are a first step in identifying ecological characteristics of those regions with the highest risk of a potentially significant role for non-human hosts. We hope that this will raise the interest of Mexican public health authorities and academic institutions in considering the risk arising from potential non-human hosts of ZIKV. Obviously our methodology is also applicable to prediction of potential non-human hosts for ZIKV in other countries, the only requirement being the relevant collection data analogous to that used in the present study. Finally, we emphasise again that the scope of the model of the present paper is to serve as a focus for future studies and show that potentially useful information can be gleaned from the method, which, at this level, is not capable of predicting detailed elements such as potential host competency.

Dengue virus in Mexican bats

Coronaviruses in bats from Mexico

Leishmania (L.) mexicana infected bats in Mexico: novel potential reservoirs

Zika virus in the Americas-yet another arbovirus threat

Viral diversity of bat communities in human-dominated landscapes in Mexico

Zika virus: rapid spread in the Western hemisphere

Genetic diversity of Didelphis virginiana related to different levels of disturbance in the Highlands and the Central Depression regions of Chiapas

Serologic evidence of Leishmania infection in free-ranging wild and domestic canids around a Brazilian National Park

Dengue infection in neotropical forest mammals

Zika virus (I). Isolations and serological specificity

Rapid risk assessment: Zika virus epidemic in the Americas: potential association with microcephaly and Guillain-Barré syndrome

Non flying mammals and landscape changes in the tropical rain forest region of Los Tuxtlas

Comparing the relative contributions of biotic and abiotic factors as mediators of species' distributions

Zika virus in Gabon (Central Africa)-2007: a new threat from Aedes albopictus

Elucidation of population connectivity in synanthropic mesopredators: using genes to define relevant spatial scales for management of raccoons and virginia opossums

Una revisió n sobre reservorios de Trypanosoma (Schizotrypanum) cruzi (Chagas, 1909), agente etioló gico de la Enfermedad de Chagas

High proportion of mosquito vectors in Zika forest, Uganda, feeding on humans has implications for the spread of new arbovirus pathogens

Bats (Chiroptera) as Vectors of Diseases and Parasites

Beyond nature and nurture: phenotypic plasticity in blood-feeding behavior of Anopheles gambiae s.s. when humans are not readily accessible

Landscape ecology of Trypanosoma cruzi in the southern Yucatan Peninsula

The emerging Zika pandemic: enhancing preparedness

Serologic evidence of flavivirus infection in bats in the Yucatan Peninsula of Mexico

Zika virus was transmitted by sexual contact in Texas, health officials report

Bat diversity and abundance as indicators of disturbance in neotropical rainforest

Rapid spread of emerging Zika virus in the Pacific area

Unexpected and rapid spread of Zika virus in the Americas-implications for public health preparedness for mass gatherings at the 2016 Brazil Olympic Games

Bats are a major natural reservoir for hepaciviruses and pegiviruses

Effect of Zika virus and Bwamba virus in the cave bat (Myotis lucifugus)

The opossum Didelphis virginiana as a synanthropic reservoir of Trypanosoma cruzi in Dzidzilché

Short report: dengue virus in bats from southeastern Mexico

Can you judge a disease host by the company it keeps? Predicting disease hosts and their relative importance: a case study for Leishmaniasis

Using biotic interaction networks for prediction in biodiversity and emerging diseases

Role of evolved host breadth in the initial emergence of an RNA virus

Zika virus in Brazil and macular atrophy in a child with microcephaly

Sylvatic focus of American trypanosomiasis in the State of

Transmission cycles, host range, evolution and emergence of arboviral disease

Report Zika Situation Neurological Syndrome and Congenital Anomalies

First report of autochthonous transmission of Zika virus in Brazil

Mosquito data were kindly provided by the National Commission for the Knowledge and Use of Biodiversity (Comision Nacional para el Conocimiento y Uso de la Biodiversidad, CONABIO) of Mexico. We deeply thank two anonymous referees for their significant comments and suggestions to improve this manuscript. We are grateful for financial support from PAPIIT-UNAM grant number IG200217. CGS is grateful for financial support from the Centro de Ciencias de la Complejidad and the Laboratorio Nacional de Ciencias de la Complejidad.