<figure><img src="image_1.png" /></figure>

STUDY PROTOCOL

# A protocol and a data-based prediction to investigate virus spillover at the wildlife interface in human-dominated and protected habitats in Thailand: The Spillover Interface project

Chuanphot Thinphovong<sup>1*</sup>, Ewan Nordstrom-Schuler<sup>1</sup>, Pipat Soisook<sup>2</sup>, Anamika Kritiyakan<sup>1</sup>, Ronnakrit Chakngean<sup>3</sup>, Sakarin Prapruti<sup>4</sup>, Malee Tanita<sup>5</sup>, Yossapong Paladsing<sup>1</sup>, Phurin Makaew<sup>1</sup>, Awatsaya Pimsai<sup>2</sup>, Abdulloh Samoh<sup>2</sup>, Christophe Mahuzier<sup>6</sup>, Serge Morand<sup>1,7,8</sup>, Kittipong Chaisiri<sup>8</sup>, Waraphon Phimpraphai<sup>9</sup>

1 Faculty of Veterinary Technology, Kasetsart University, Bangkok, Thailand, 2 Princess Maha Chakri Sirindhorn Natural History Museum, Prince of Songkla University, Songkla, Thailand, 3 Nanthaburi Nation Park, Thawangpha, Nan, Thailand, 4 Protected Areas Regional Office 13 (Phrae), Phrae, Thailand, 5 Primary Care Unit (PCU), Saenthong, Thawangpha, Nan, Thailand, 6 Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES Paris)—Centre de Recherche IRD, Montpellier, France, 7 MIVEGEC, CNRS-IRD-MUSE, Montpellier Université, Montpellier, France, 8 Faculty of Tropical Medicine, Department of Helminthology, Mahidol University, Bangkok, Thailand, 9 Faculty of Veterinary Medicine, Kasetsart University, Bangkok, Thailand

* chuanphot.2@gmail.com

<figure><img src="image_2.png" /></figure>

OPEN ACCESS

**Citation:** Thinphovong C, Nordstrom-Schuler E, Soisook P, Kritiyakan A, Chakngean R, Prapruti S, et al. (2024) A protocol and a data-based prediction to investigate virus spillover at the wildlife interface in human-dominated and protected habitats in Thailand: The Spillover Interface project. PLoS ONE 19(1): e0294397. https://doi.org/10.1371/journal.pone.0294397

**Editor:** Nussieba A. Osman, Sudan University of Science and Technology, College of Veterinary Medicine, SUDAN

**Received:** January 31, 2023

**Accepted:** October 27, 2023

**Published:** January 2, 2024

**Copyright:** © 2024 Thinphovong et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** We used an open database compiling the associations between 1,785 virus species (DNA and RNA) and 725 mammalian host species, resulting from the automatic screening of metadata accompanying the nucleotide sequences published between 1950 and 2019? This database is accessible on the Dryad Digital Repository: https://doi.org/10.5061/dryad.p2ngf1vmg The mammal species potentially

## Abstract

The Spillover Interface Project aims at assessing the encounter of wildlife, domestic animals, and humans along a landscape gradient from a protected area to a residential community, through areas of reforestation and agricultural land. Here, we present the protocols of the project that combine virus screening in humans, bats, rodents and dogs with camera trapping, land-use characterization, and network analyses. The project is taking place in the sub-district of Saen Thong (Nan Province, Thailand) in collaboration with local communities, the District Public Health Office, and Nanthaburi National Park. To formulate a predictive hypothesis for the Spillover Interface Project, we assess the wildlife diversity and their viral diversity that could be observed in Saen Thong through a data science analysis approach. Potential mammalian species are estimated using data from the International Union for Conservation of Nature (IUCN) and their associated viral diversity from a published open database. A network analysis approach is used to represent and quantify the transmission of the potential viruses hosted by the mammals present in Saen Thong, according to the IUCN. A total of 57 viruses are expected to be found and shared between 43 host species, including the domestic dog and the human species. By following the protocols presented here, the Spillover Interface Project will collect the data and samples needed to test this data-driven prediction.
present in Saen Thong were assessed using the
open data from IUCN. Shapefiles of terrestrial
mammal species were downloaded from the IUCN
Red List at: https://www.iucnredlist.org/resources/
spatial-data-download.

**Funding:** National Research Council of Thailand, N42A650216, Kittipong Chaisiri, the funders had and will not have a role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing interests:** The authors declare no competing interest.

**Introduction**

For several decades, the number of new emerging diseases as well as the number of epidemics of infectious diseases have continued to increase [1,2]. The vast majority of these diseases are zoonoses, involving domestic and/or wild animals [3]. The emergence of new zoonotic diseases begins with the spillover of a pathogen hosted in an animal reservoir [4], which is essential for maintaining the zoonotic agent [5]. The animal reservoir is assumed to be asymptomatic [6] and living in close relationship with humans [7]. Thus, emerging zoonotic diseases are primarily driven by interspecific interactions described as spillover events [8].

Spillover events are determined by a combination of biological and ecological processes in relation to hosts, pathogens, vectors and environmental conditions. Several studies have investigated the conditions that favor spillover [9,10] and, importantly, where a spillover has a greater chance to occur [6,11,12]. Spillover events may depend on various factors, from habitat change [13], climate variability [14], to social factors [15], which favor the encounter of reservoir and recipient hosts [16]. Habitat changes with agriculture intensification [17,18], livestock expansion [19], urbanization [20,21] and forest conversion [22] enhance the spillover and transmission of zoonotic diseases [23,24], often through generalist and synanthropic species [25]. The various combinations of these factors constitute the encounter interface.

An encountering interface can be classified as obvious or non-obvious. An obvious encountering interface is a location where interspecific transmission of a pathogen is predictable over time due to the continuous presence of both reservoir and recipient hosts. Wildlife markets or exploited caves for bat guano can be defined as obvious encountering interfaces. In contrast, non-obvious interfaces, such as ecotones related to human activities and rapid land use changes, can be characterized as a transient location with a narrow transmission and spillover opportunity due to short interspecific interactions [26,27]. This interface is difficult to identify because it requires ecological knowledge of the species of interest in order to assess their affinity with different habitats and their responses to environmental factors.

Among the new emerging viral zoonoses, bats harbor a high proportion of viruses [28-30], with 224 species of bats harboring 61 zoonotic diseases [29]. Rodents are also important, with 217 species harboring 66 zoonotic viruses [31]. Among these viruses, coronaviruses (CoVs) have received a great deal of attention [32]. In the last sixteen years, three coronaviruses of bat origin with high epidemic potential have emerged in human populations (SARS-CoV, MERS-CoV, and SARS-CoV-2) and one in domestic animals (SADS-CoV) [33]. Coronaviruses are also well known in veterinary medicine as they infect a wide range of mammalian hosts, such as Transmissible Gastroenteritis Coronavirus (TGEV) of domestic pigs [34], Mouse Hepatitis Virus (MHV) of mice [35], and Canine Coronavirus (CCoV) of dogs [36]. In Malaysia, Canine CoV (CCoV), a novel canine-feline recombinant alphacoronavirus was found in nasopharyngeal swab samples from eight out of 301 patients hospitalized with pneumonia, most of them children living in rural areas [37].

Reverse spillover, called spill-back, can occur from the human population to domestic animals or wildlife, as recently observed for SARS-CoV-2 in several mammal species [38], including white-tailed deer [39]. In Thailand, three dogs out of 35 examined were found infected by SARS-CoV-2 from households with confirmed COVID-19 residing patients [40]. The great diversity of coronaviruses is associated with their evolutionary capacity and ability to jump between species [41], which gives them an exceptional propensity for spillover and emergence [42].

The Spillover Interface Project aims at assessing virus sharing at the interfaces of wildlife, domestic animals, and humans along a gradient from a protected area to a village community and their agricultural land. The project will combine virus screening in various hosts in relation to land-use characterization [43,44].
Here, we present the objectives, design, and protocols used in the Spillover Interface Project. To formulate a predictive hypothesis for the Spillover Interface Project, we develop a predictive hypothesis based on open data to assess the wildlife diversity and associated viral diversity that can be observed in the implementing locality of the project.

# Materials and methods

## A. Spillover Interface Project design

**Project objectives.** The objectives are (1) to investigate the diversity of viruses, and in particular coronaviruses, in humans, domestic dogs, and wildlife (bats and rodents), (2) to identify the potential transmission of viruses among hosts using network analyses, (3) to identify potential interfaces using a land use map, (4) to assess the diversity of wildlife, their habitats and their potential role in virus spillover.

**Location of the study.** Since 2012, our research team has conducted several collaborative studies with local communities and local administrations, such as the District of Public Health Office and the Nanthaburi National Park in the sub-district of Saen Thong (Nan Province, Thailand). The sub-district is divided into two types of landscape, with a lowland agricultural zone (with four villages) close to an urbanized habitat and an upland agricultural zone (four villages) close to the forested area of Nanthaburi National Park [45]. The upland part of the sub-district provides an ideal site for the implementation of the Spillover Interface Project, with a gradient from the protected area of Nanthaburi National Park, with a cave populated by bats, towards the reforestation area, plantations, and agricultural land (Fig 1A and 1B).

<figure><img src="image_3.png" /><figcaption>**Fig 1.** Location of the Spillover Interface project in (A) Saen Thong sub-district (Nan province, Thailand) precisely (B) in the upland part of the subdistrict. A land-use cover describes the different land classes: Multi-specific forests, plantations (rubber, teak, bamboo plantations, orchards), fallows crops (corn, paddy rice, ginger, etc.), urban infrastructural at (C) the level of the sub-district and (D) the upland part of the sub-district.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0294397.g001
**Land use map.** A high-resolution land use map (10 square meters) of Saen Thong subdistrict was developed for a previous project using Copernicus satellite data [46]. Random Forest classifications used the combination of Sentinel-2 optical and Sentinel-1 radar images with the different images acquired in 2019. The classifications cover four categories of lands: (1) unculminated steep mountain slopes with forest (protected areas and community forests), (2) cultivated land mostly on steep slopes with annual crops (such as upland rice and corn), or plantations (such as rubber, teak and bamboo), (3) cultivated lowland with paddy fields, (4) urbanization and infrastructures. All farmers are smallholders, and field sizes are relatively small. The classifications were validated by ground-truth checks [46].

The land use classification makes it possible to differentiate multi-specific and heterogeneous forests (mixed forests), including reforestation areas and community forests, from plantations such as rubber, teak, bamboo or orchards. The land use also distinguishes dominant crops (corn, paddy rice, ginger, etc.), houses, and other infrastructure (Fig 1C and 1D). A terrain elevation model was also developed [46]. This land use map is important (1) to assess the typology of the different habitats of the locality allowing statistical comparisons and (2) to describe the behavior of animals at high resolution, ie. the movement or home range of animals in a small buffer zone around each camera trap or live trap by tags for rodents or pictures for dogs and other wildlife when individual recognition is possible.

**Framework.** A framework was designed for sampling and data collection in order to validate the predictive hypothesis. Rodents will be captured by live traps and bats will be trapped by mist nets and harp traps. Sample collection will include all swabs and guano samples collected for virus screening. Data collection will include qualitative data obtained from questionnaires and focus group discussions with villagers, as well as quantitative data such as weight and body measurements obtained from rodents and bats during the trapping sessions. GPS coordinates will be recorded at individual traps and camera trapping devices. In the field, all information will be collected digitally (Fig 2).

**Research and ethical approvals.** The principles of ethics and responsibility require local community and stakeholder engagement, such as in the Primary Health Care Unit, the District Public Health Office, and Nanthaburi National Park. Several meetings with all stakeholders helped to co-construct the project which made it possible to meet the objectives of human health, animal health, and biodiversity conservation.

Procedures for human health investigation, laboratory investigation, interviews, and questionnaires will be sent for approval by the ethical committee of the Nan Provincial Public Health Office, Ministry of Public Health, Thailand (NAN REC 63-13). All participants who agree to join the study will be asked to read the Participant Information Form, which explaining the objectives, procedures, possible risks, and benefits of the research project.

Similarly, procedures for collection of samples from dogs, laboratory investigation, safety procedures, interviews and questionnaires of dog owners will be sent for approval by the Institutional Animal Care and Use Committee, Kasetsart University (ACKU64-VTN-010).

Rodent and bat species that will be investigated by Spillover Interface Project will be neither on the CITES list nor on the Red List (IUCN). Animals will be treated in accordance with the guidelines of the American Society of Mammalogists, and within the European Union legislation guidelines (Directive 86/609/EEC). Any trapped rodent species listed on CITES will be released without being manipulated. All animals will be released after sampling. Approval notices for trapping and investigation of bats and rodents will be sent for approval by the Institutional Animal Care and Use Committee, Kasetsart University.

Permit has been approved by the Department of National Parks, Wildlife and Plant Conservation (DNP) and the Royal Forest Department, Ministry of Natural Resources and Environment, Thailand. In addition, foreign researchers participating in the research project have also
<figure><img src="image_4.png" /><figcaption>Fig 2. Framework for sampling, swabs collected for virus screening, and data collected from questionnaires and focus groups from villagers, from rodent and bat trappings and from camera trappings (with associated GPS locations for each device). Data will be collected in the field using a data form designed from the Epicollect5 application installed on smart phones.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0294397.g002
been approved by the National Research Council of Thailand (NRCT) with the consent of the DNP.

With community engagement, ethics, and research permits, the investigation will begin with the organization of a meeting involving the researchers, local public health workers, the village leaders, and village volunteers to remind them of the purpose and details of the study. Interviews with villagers will help to describe their interactions with wildlife. Moreover, it will help to select participants who will express close contact with wildlife, either bats or rodents.

## Field implementation and virus screening protocols

**Bat sampling.** Bats will be sampled once inside a cave situated in Nanthaburi National Park and in the village of Santisuk using mist nets and harp traps [47]. Traps will be set before sunset and checked after sunset depending on bat activity time. Measurements will be taken for each captured bat: weight, forearm, hindfoot, tail, ear, head to body, and tibia lengths. Saliva and rectal swabs will be collected, and stored in RNA later and kept at -20°C. For other sessions, only guano will be collected from the cave.

**Rodent sampling.** The rodent trapping is designed to cover a gradient from the village to the bat cave situated in the Nanthaburi National Park area. Rodent trapping will then cover several habitat types, including agricultural crops, plantations, and the reforestation area at the interface between plantations and the National Park. Live-traps, locally made (dimension: 14 X 14 X 27 cm) will be used. A total of 210 traps will be set: 40 traps in and around the village, 60 traps in agricultural land and plantations, 100 traps in the reforestation area, and 10 traps inside the cave. Traps will be set in lines of five traps separated by a distance of 10 meters. The different lines will be set to evenly cover the gradient. Trapping will be conducted over a period of four nights for each session, which corresponds to a total of 840 trap nights per trapping session. The traps will be placed in the same positions for each session by using plastic plant labels as notification tags on the ground. Pictures, habitat descriptions and rodent trap coordinates available for the CERoPath Project (https://healthdeep.shinyapps.io/Small_mammals_CERoPath/) will be used for rodent species identification [44].

Rodents will be anesthetized using a mixture of 20% v/v isoflurane in propylene glycol. Two ml of mixed isoflurane will be placed on cotton or gauze in a transparent plastic box chamber (32 X 22 X 15 cm). The individual animal will be introduced into the chamber, and respiratory rate will be observed after anesthetizing for one minute in the case of a mouse-size rodent and for two minutes for a rat-size rodent [48,49]. The animal will be removed from the anesthetic chamber, its weight will be recorded, and pictures will be taken for further body measurements. Saliva and rectal swabs will be taken, and stored in RNA later and kept at -20°C. A microchip (Passive Integrated Transponder: PIT Tag) of 1.4 X 8 mm will be injected subcutaneously for individual permanent identification and for capture-mark-recapture (Absonutrix Sale Company, ISO11784/785). Captured animals will be read with a microchip reader (Absonutrix Sale Company, mini scanner ISO FDX-B 134.2 KHz) to evaluate re-trapping or new trapping. Rodents will be carefully monitored before releasing them at their trap location captures. Any re-trapped individual during a session will be released without manipulation.

**Dog sampling.** Most dogs roam freely in the upland villages of Saen Thong and may go to agricultural or reforestation areas with their owners. Free-roaming dogs will be selected with the help of the village leader, who will invite dog owners to freely join the research project. The participant dog owners will take their dogs to the village hall, where we will collect nasopharyngeal swabs, conjunctival, swabs and rectal swab with the consent of their owners. Nasopharyngeal, conjunctival and rectal swabs that will be placed in cryotube contained 0.5 µl of RNA
later, will be stored in RNA and kept at -20°C. Conjunctival swabs will be used to screen asymptomatic [50] as well as symptomatic cases of COVID-19 [51].

**Virus screening.** For RNA extraction, the QIAamp Viral RNA Mini Kit (Qiagen) will be used following the manufacturer's instructions. The extracted RNA will be measured for concentration, purity and quality checking using Nanodrop [41]. Concerning CoVs, extracted RNA will be used for reverse transcriptase (RT) reactions using SuperScript III Reverse Transcriptase (Life tech) and hexamers with primers and semi-nested PCR protocols following Gouilh et al. [43]. We will keep extracted DNA and RNA extraction for further viral analyses [52].

**Camera trapping.** Camera trapping helps to investigate wildlife diversity along the gradient from the protected area, especially species visiting the bat cave, to the reforestation area and agricultural land. A total of 40 camera traps will be positioned following a land use map illustrating forest or reforested areas, rodent trapping lines, and suggestions from the rangers of Nanthaburi National Park and the village volunteers based on their ecological knowledge. The camera traps will be left in place for a year. The camera traps will be checked every two months. Batteries will be changed, and pictures taken by the camera traps will be transferred to a laptop computer. Later, the pictures will be sorted, and the identification of species will be assessed by a consensus of experts. The sorted pictures by species and by camera traps will be analyzed using the 'camtrapR' package [53] implemented in R, which will allow the exploration of the spatio-temporal activities of animals, including roaming dogs.

## B. Predictive hypothesis

**Potential wildlife species and viruses.** We developed a predictive modeling that helps formulate hypotheses to be tested by the Spillover Interface Project, such as which mammal species and which viruses can be found in the study locality. To describe the potential association between mammal species and virus, we estimated the potential presence of mammal species in the study site using data from the International Union for Conservation of Nature (IUCN). Shapefiles of terrestrial mammal species were downloaded from the IUCN Red List (https://www.iucnredlist.org/resources/spatial-data-download) providing geographical distribution for each animal species. The IUCN Red List the status of each mammal species and their CITES status. The list of mammal species from the orders Carnivora, Cetartiodactyla, Chiroptera, Pholidata, Rodentia, and Scandentia, that are potentially present in Saen Thong was then extracted using the land use map presented above. Animal species were identified on the basis of their overlap distribution, given by the IUCN shapefiles, with Saen Thong area given by the land-use map [46]. Concerning the viral diversity, we used an open database of nucleotide sequences published between 1950 and 2019 compiling the associations between 1,785 virus species (DNA and RNA) with 725 mammalian host species that resulted from the automatic screening of the accompanying metadata [54,55].

**Network analyses.** We used a network analysis approach to represent and quantify the ecology of virus transmission between different hosts [56], as it allowed representation and mathematical quantification of the importance of a given host in the transmission of viruses. Viral diversity and potential transmission to humans were assessed using network architectures and associated indices. In bipartite networks, only host-virus interactions are considered, representing how host species and virus species are associated. Viruses and hosts represent nodes of the network, and edges represent the observed interactions between nodes in the network. In a unipartite network, nodes represent the hosts and edges the viruses shared between each node (i.e. host).

Modularity in bipartite and unipartite networks of shared viruses help assess host communities sharing similar virus transmission and then potential risks of transmission of other
viruses between species [57]. Modularity measures the strength of division of a network into modules, i.e. groups or clusters of hosts sharing common viruses.

Network centrality indices provided useful information on the relative importance of a given host to the entire network [58]. A host occupying a highly central position (i.e., having a high centrality value) links different hosts clustered into subgroups within the network and participates strongly in the transmission among a large community of hosts. Hosts with high values of centrality in networks can help target key reservoirs for viral disease surveillance [54,56,57].

We first linked the list of mammal species extracted from the IUCN, adding the domestic dog (*Canis lupus familiaris*) and the human species (*Homo sapiens*), to the dataset of virus species from mammals [47]. Using the information linking each host species with their virus, we obtained bipartite networks and unipartite projections in which each node was a host species. Modules were identified for bipartite and unipartite networks of shared viruses between mammal species. We used bipartite network analysis, with nodes describing the mammal species interacting with nodes describing the viruses using the 'bipartite' package [58,59] implemented in R [60]. We used the function 'computeModules' of the package 'bipartite' to compute modules using the modularity algorithm of Dormann & Strau [61]. We projected these bipartite networks onto unipartite networks using the 'tnet' package [62] implemented in R. A unipartite network represented relative interactions amongst hosts through the sharing of virus species. Each host within the network played a different role in virus sharing relative to all other nodes, which was examined using its centrality measurement. We used the function 'cluster_louvain' implemented in the package 'igraph' [63] to identify the modularity structure of the unipartite networks. This function is based on a multilevel modularity optimization algorithm [64]. A central node (i.e. a host with a high centrality value) was the one that was highly connected to other nodes and thus was supposed to have a greater transmission potential for virus species. We calculated the eigenvalue centrality (EC) with the 'evcent' function from the package 'igraph' [63].

## Results

The predictive modeling used network analysis to explore the association between mammal species and their viruses. A total of 130 mammal species are potentially present in Saen Thong sub-district according to the IUCN Red List, although we removed several species known for their absence in the area (such as *Mus musculus* and *Rattus rattus*). The species were checked and confirmed using the CERoPath database and the data collected by the Smart Patrol Project of Nanthaburi National Park. This list of mammal species, also including the domestic dog (*Canis lupus familiaris*) and the human species (*Homo sapiens*), was linked to the dataset of virus species from mammals. A total of 57 viruses were found to be shared between 43 host species (see S1 and S2 Tables), including the domestic dog and the human species (Fig 3A), with the human species hosting the highest number of viruses followed by the domestic dog, the Indochinese rhesus macaque (*Macaca mulatta*), and several bat species. Bats showed the highest diversity of viruses shared among the whole host communities, followed by human and non-human primates, carnivores, rodents, and cetartiodactylids (Cervidae) (Fig 3B).

Bipartite networks and unipartite projections were obtained using the information linking each host species with their virus. Modules identified for bipartite and unipartite networks differed by their numbers of modules (Figs 4 and 5). Using the bipartite network, five modules were identified (Fig 4), such as the one grouping the human species, domestic dog,
<figure><img src="image_5.png" /><figcaption>**Fig 3.** Diversity of potential viruses hosted in species potentially present in Saen Thong with (A) the number of viruses per mammal species including the human species and the domestic dog and (B) the number of viruses per mammal orders. Data on the potential presence of wild mammals were extracted from the IUCN Red List and data on the potential viruses from an open database.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0294397.g003

Indochinese rhesus macaque, and Sambar deer (*Rusa unicolor*) with ten virus species. The bat species were separated into three modules, and the rodent species into two modules.

Four modules were identified using the unipartite network. The first module grouped the human species with domestic dog, Indochinese rhesus macaque, Sambar deer, and two species of rodents (*Rattus exulans*, *Rattus nitidus*). The bat species were separated into three modules, with one comprising two rodent species (*Rattus tanezumi*, *Rattus losea*), and one comprising three carnivore species (*Paguma larvata*, *Paradoxurus hermaphroditus*, *Viverricula indica*) (Fig 5).

In terms of host centrality in the network of virus sharing, the primates showed the highest centrality because of the human species, followed by carnivores due to the domestic dog, the deer species (Cetartiodactylida), the rodents, and the bats (Fig 6).

Discussion

Pathogen spillover events depend on various factors that favor the encounter of reservoir hosts and recipient hosts. Encountering interfaces is not always easy to characterize, mostly because of the lack of ecological knowledge about the species of interest. The protocols of the Spillover Interface Project will combine pathogen screening in humans, domestic dogs, bats, and
<figure><img src="image_6.png" /><figcaption>**Fig 4.** Bipartite network of potentially shared virus species among wild mammal species potentially present in Saen Thong and the human species and the domestic dog. Five modules were identified with one the one grouping humans, domestic dog, Indochinese rhesus macaque, and the Sambar deer with ten viruses. Bat species were separated in three modules, and rodent species in two modules.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0294397.g004

rodents with extensive camera trapping, land use characterization, and network analyses in order to assess the encountering of wildlife, domestic animals and humans.

Using a data-based approach, we have made predictions about the wildlife diversity and viral diversity that can be observed in Saen Thong. The use of IUCN and virus open databases helped to predict the potential diversity of virus-mammal interactions in the locality. Using network analysis, we showed the importance of bat species as major reservoirs of potential zoonotic viruses. However, carnivores, along with domestic dogs, appeared central in the virus sharing among all mammal species [57,65]. The predictive analysis using open access data and network analysis showed the importance of several species. Synanthropic rodents and the domestic dog appeared central in the network of sharing viruses among humans and wildlife species that are potentially present at the location site. This data-based prediction confirmed
<figure><img src="image_7.png" /><figcaption>**Fig 5.** (A) Unipartite network with modules differentiated by colors of shared virus species among mammal species including humans and domestic dog. The links among hosts (nodes) of the unipartite network depict shared viruses (vertices were placed according to the Fruchterman–Reingold algorithm) with thickness of links proportional to number of viruses shared and size of vertices proportional to the degree centrality of hosts. (B) Clustering representation of the unipartite network with the different clusters differentiated by colors. Four modules were identified with one grouping the human species with the domestic dog, the Indochinese rhesus macaque, the Sambar deer and two species of rodents (*Rattus exulans*, *Rattus nitidus*). The bat species were separated in three modules, with one comprising two rodent species (*Rattus tanezumi*, *Rattus losea*), and one comprising three carnivore species (Paguma larvata, *Paradoxurus hermaphroditus*, *Viverricula indica*).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0294397.g005

the importance of investigating not only bats and rodents but also the domestic dog, which are then the targeted species of the Spillover Interface Project.

The results of this data-based prediction will be interrogated with the data gathered by the Spillover Interface Project. First, with the data on the diversity and co-occurrence of wildlife species in the different habitats gathered by the camera and live traps. The land-use map will allow the precise identification of habitat use and shared by each species trapped or pictured, and potential interfaces. Second, the screening of viruses will permit the investigation of their diversity and the number of species that individual viral species are able to infect such as small terrestrial mammals, bats, dogs, and importantly, across habitats. Last, the trapping design will standardize the sampling effort to assess the presence and abundance of the investigated small mammals and domestic dogs in each habitat.

## Conclusion

We presented the results of a predictive modeling approach illustrating the potential presence and association of hosts and their viruses in Saen Thong using network analysis. The protocols and framework of the Spillover Interface Project were designed to test the results of the predictive modeling in the Spillover Interface Project, first starting with coronaviruses.
<figure><img src="image_8.png" /><figcaption>**Fig 6.** Box and whiskers plots of host centrality measures for sharing virus. Notched boxes for orders of mammals correspond to 95% confidence interval (horizontal lines inside the boxes are the medians of the centrality values, grey dots are means of centrality values, standard deviation are represented by vertical lines).</figcaption></figure>

https://doi.org/10.1371/journal.pone.0294397.g006

**Supporting information**

S1 Table. The list of mammal species in Saen Thong, Thawangpha, Nan Thailand that was reported in the International Union for Conservation of Nature (IUCN).
(CSV)

S2 Table. The list of viruses which link to mammal species in Saen Thong, Thawangpha, Nan Thailand. Their associated viruses from a published open database.
(CSV)

**Acknowledgments**

This project will benefit from the support and help of the staff of Nanthaburi National Park and local villagers for survey of forest area and camera trap setting. We also thank Phutita Wongwaiyut and Pachara Prompnopwong from Princess Maha Chakri Sirindhorn Natural History Museum, Prince of Songkla University.
## Author Contributions

**Conceptualization:** Chuanphot Thinphovong, Serge Morand, Kittipong Chaisiri, Waraphon Phimpraphai.

**Data curation:** Serge Morand.

**Funding acquisition:** Kittipong Chaisiri.

**Methodology:** Chuanphot Thinphovong, Pipat Soisook, Anamika Kritiyakan, Ronnakrit Chakngean, Sakarin Prapruti, Malee Tanita, Yossapong Paladsing, Phurin Makaew, Awatsaya Pimsai, Abdulloh Samoh, Christophe Mahuzier, Serge Morand, Kittipong Chaisiri, Waraphon Phimpraphai.

**Visualization:** Serge Morand.

**Writing – original draft:** Chuanphot Thinphovong, Ewan Nordstrom-Schuler, Serge Morand.

## References

1. Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL et al. Global trends in emerging infectious diseases. Nature. 2008; 451: 990–993. https://doi.org/10.1038/nature06536 PMID: 18288193

2. Smith KF, Goldberg M, Rosenthal SCarlson LChen J Chen C et al. Global rise in human infectious disease outbreaks. J R Soc Interface. 2014; 11: 20140950. https://doi.org/10.1098/rsif.2014.0950 PMID: 25401184

3. Wourc'h G, Moutou F, Morand S, Jourdain E. Zoonoses, Diseases that Link Animals and Humans. 2022; Versailles: Quae.

4. Daszak P, Cunningham A, Hyatt AD. Emerging infectious diseases of wildlife-threats to biodiversity and human health. Science. 2000; 287: 443–449. https://doi.org/10.1126/science.287.5452.443 PMID: 10642539

5. Mollentze N, Streicker DG. Viral zoonotic risk is homogenous among taxonomic orders of mammalian and avian reservoir hosts. Proc Natl Acad Sci USA. 2020; 117: 9423–9430. https://doi.org/10.1073/pnas.1919176117 PMID: 32284401

6. Plowright RK, Parrish CR, McCallum H, Hudson PJ, Ko AI, Graham ALel at. Pathways to zoonotic spillover. Nat Rev Microbiol. 2017; 15: 502–510. https://doi.org/10.1038/nrmicro.2017.45.

7. Maia C, Campino L. Can domestic cats be considered reservoir hosts of zoonotic leishmaniasis?. Trends Parasitol/ 2011; 27: 341–344. https://doi.org/10.1016/j.pt.2011.03.008 PMID: 21570915

8. Power AG, Mitchell CE. Pathogen spillover in disease epidemics. Am Nat. 2004; 164 Suppl 5: S79–89. https://doi.org/10.1086/424610 PMID: 15540144

9. Olson SH, Benedum CM, Mekaru SR, Preston ND, Mazet JA, Joly DO, et al. Drivers of emerging infectious disease events as a framework for digital detection. Emer Infect Dis. 2015; 21:1285–92. https://doi.org/10.3201/eid2108.141156 PMID: 26196106

10. Borremans B, Faust C, Manlove KR, Sokolow SH, Lloyd-Smith JO. Cross-species pathogen spillover across ecosystem boundaries: mechanisms and theory. Phil Trans R Soc B. 2019; 374: 20180344. https://doi.org/10.1098/rstb.2018.0344 PMID: 31401953

11. Becker DJ, Washburne AD, Faust CL, Mordecai EA, Plowright RK. The problem of scale in the prediction and management of pathogen spillover. Phil Trans R Soc B. 2019; 374: 20190224. https://doi.org/10.1098/rstb.2019.0224 PMID: 31401958

12. Plowright RK, Reaser JK, Locke H, Woodley SJ, Patz JA, Becker DJ, et al. Land use-induced spillover: a call to action to safeguard environmental, animal, and human health. Lancet Planet Health. 2021; 5: e237–e245. https://doi.org/10.1016/S2542-5196(21)00031-0 PMID: 33684341

13. Murray KA, Daszak P. Human ecology in pathogenic landscapes: two hypotheses on how land use change drives viral emergence. Curr Opin Virol. 2013; 3: 79–83. https://doi.org/10.1016/j.coviro.2013.01.006 PMID: 23415415

14. Latinne A, Morand S, Climate anomalies and spillover of bat-borne viral diseases in the Asia-Pacific region and Arabian Peninsula. Viruses. 2022; 14: 1100. https://doi.org/10.3390/v14051100
15. Weiss R, McMichael A. Social and environmental risk factors in the emergence of infectious diseases. Nature Med. 2004; 10: S70-S76. https://doi.org/10.1038/nm1150 PMID: 15577934

16. Holmes EC, Rambaut A, Andersen KG. Pandemics: spend on surveillance, not prediction. Nature 2018; 558: 180-182. https://doi.org/10.1038/d41586-018-05373-w PMID: 29880819

17. Shah HA, Huxley P, Elmes J, Murray KA. Agricultural land-uses consistently exacerbate infectious disease risks in Southeast Asia. Nature Comm. 2019; 10: 4299. https://doi.org/10.1038/s41467-019-12333-z PMID: 31541099

18. Morand S. The role of agriculture in human infectious disease outbreaks. CABI Review. 2022, https://doi.org/10.1079/cabireviews202217060.

19. Morand S. Emerging diseases, livestock expansion and biodiversity loss are positively related at global scale. Bio Cons. 2020; 248: 108707. https://doi.org/10.1016/j.biocon.2020.108707 PMID: 32834060

20. Hassell JM, Begon M, Ward MJ, Fèvre EM. Urbanization and disease emergence: Dynamics at the wildlife-livestock-human interface. Trends Ecol Evol. 2017; 32: 55-67. https://doi.org/10.1016/j.tree.2016.09.012 PMID: 28029378

21. Dobigny G, Morand S. Zoonotic emergence at the animal-environment-human interface: the forgotten urban socio-ecosystems. Peer Community Journal. 2022; 2; e79. https://peercommunityjournal.org/articles/10.24072/pcjournal.206/.

22. Morand S, Lajaunie C. Outbreaks of vector-borne and zoonotic diseases are associated with changes in forest cover and oil palm expansion at global scale. Front Vet Sci. 2021; 8: 661063. https://doi.org/10.3389/fvets.2021.661063 PMID: 33842581

23. Faust CL, Dobson AP, Gottdenker N, Bloomfield LSP, McCallum HI, Gillespie TR, et al. Null expectations for disease dynamics in shrinking habitat: dilution or amplification? Phil Trans R Soc B. 2017; 372: 20160173. https://doi.org/10.1098/rstb.2016.0173 PMID: 28438921

24. Rulli MC D'Odorico P Galli N, Hayman DTS. Land-use change and the livestock revolution increase the risk of zoonotic coronavirus transmission from rhinolophid bats. Nature Food. 2021; 2: 409-16. https://doi.org/10.1038/s43016-021-00285-x PMID: 37118224

25. Gibb R, Redding DW, Chin KQ, Donnelly CA, Blackburn TM, Newbold T, et al. Zoonotic host diversity increases in human-dominated ecosystems. Nature. 2020; 584: 398-402. https://doi.org/10.1038/s41586-020-2562-8 PMID: 32759999

26. Despommier D, Ellis BR, Wilcox. The Role of Ecotones in Emerging Infectious Diseases. Ecohealth. 2006; 3; 281-289. https://doi.org/10.1007/s10393-006-0063-3.

27. Borremans B, Faust C, Manlove KR, Sokolow SH, Lloyd-Smith JO. Cross-species pathogen spillover across ecosystem boundaries: mechanisms and theory. Phil Trans R Soc London B. 2019; 374: 20180344. https://doi.org/10.1098/rstb.2018.0344 PMID: 31401953

28. Calisher CH, Childs JE, Field HE, Holmes KV, Schountz T. Bats: important reservoir hosts of emerging viruses. Clin Microbiol Rev. 2006; 19: 531-545. https://doi.org/10.1128/CMR.00017-06 PMID: 16847084

29. Luis AD, Hayman DTS, O'Shea TJ, Cryan PM, Gilbert AT, Pulliam JRC et al. A comparison of bats and rodents as reservoirs of zoonotic viruses: are bats special? Proc Biol Sci. 2013; 280: 20122753. https://doi.org/10.1098/rspb.2012.2753 PMID: 23378666

30. Olival K, Hosseini P, Zambrana-Torrelio C, Ross N, Bogich TL, Daszak P. Host and viral traits predict zoonotic spillover from mammals. Nature. 2017; 546: 646-650. https://doi.org/10.1038/nature22975 PMID: 28636590

31. Han BA, Schmidt JP, Bowden SE, Drake JM. Rodent reservoirs of future zoonotic diseases. Proc Natl Acad Sci USA. 2015; 112: 7039-7044. https://doi.org/10.1073/pnas.1501598112 PMID: 26038558

32. Corman VM, Muth D, Niemeyer D, Drosten C. Hosts and sources of endemic human coronaviruses. Adv Virus Res. 2018; 100: 163-188. https://doi.org/10.1016/bs.aivir.2018.01.001 PMID: 29551135

33. Li W, Shi Z, Yu M, Ren W, Smith C, Epstein JH et al. Bats are natural reservoirs of SARS-like coronaviruses. Science. 2005; 310: 676-679. https://doi.org/10.1126/science.1118391 PMID: 16195424

34. Liu Q, Gerdts V. Transmissible Gastroenteritis Virus of pigs and Porcine Epidemic Diarrhea Virus (Coronaviridae). Encyclopedia of Virology. 2021; 850-853. https://doi.org/10.1016/B978-0-12-809633-8.20928-X).

35. Körner RW, Majjouti M. Alcazar MA, Mahabir E. Of mice and men: The Coronavirus MHV and mouse models as a translational approach to understand SARS-CoV-2. Viruses. 2020; 12: 880. https://doi.org/10.3390/v12080880 PMID: 32806708

36. Buonavoglia C, Decaro N, Martella V, Elia G, Campolo M, Desario C, et al. Canine coronavirus highly pathogenic for dogs. Emerg Infect Dis. 2006; 12: 492-494. https://doi.org/10.3201/eid1203.050839 PMID: 16704791
37. Vlasova AN, Diaz A, Damtie D, Xiu L, Toh TH, Lee JS, a, et al. Novel Canine Coronavirus isolated from a hospitalized patient with pneumonia in East Malaysia. Clin Infect Dis. 2022; 74: 446–454. https://doi.org/10.1093/cid/ciab456 PMID: 34013321

38. Fagre AC, Cohen LE, Eskew EA, Farrell M, Glennon E, et al. Assessing the risk of human-to-wildlife pathogen transmission for conservation and public health. Ecol Lett. 2022; 25: 1534–1549. https://doi.org/10.1111/ele.14003 PMID: 35318793

39. Hale VL, Dennis PM, McBride DS, Nolting JM, Madden C, Huey D, et al. SARS-CoV-2 infection in free-ranging white-tailed deer. Nature. 2022; 602: 481–486. https://doi.org/10.1038/s41586-021-04353-x PMID: 34942632

40. Jairak W, Charoenkul K, Chamsai E, Udom K, Chaiyawong S, Bunpapong N, et al. First cases of SARS-CoV-2 infection in dogs and cats in Thailand. Transbound Emerg Dis. 2022; 69: e979–e991. https://doi.org/10.1111/tbed.14383 PMID: 34739748

41. Ar Gouilh M., Puechmaille SJ, Diancourt L, Vandenbogaert M, Serra-Cobo J, Lopez Roïg, M, et al. SARS-CoV related Betacoronavirus and diverse Alphacoronavirus members found in western old-world. Virology. 2018; 517: 88–97. https://doi.org/10.1016/j.virol.2018.01.014 PMID: 29482919

42. Wang LF, Anderson DE. Viruses in bats and potential spillover to animals and humans. Curr Opin Virol. 2019; 34: 79–89. https://doi.org/10.1016/j.coviro.2018.12.007 PMID: 30665189

43. Mar Gouilh, Puechmaille SJ Gonzalez JP, Teeling E, Kittayapong P, Manuguerra JC. SARS-Coronavirus ancestor's foot-prints in South-East Asian bat colonies and the refuge theory. Infect Genet Evol. 2011; 1: 1690–1702. https://doi.org/10.1016/j.meegid.2011.06.021.

44. Morand S, Blasdell K, Bordes F, Buchy P, Carcy B, Chaisiri K, et al. Changing landscapes of Southeast Asia and rodent-borne diseases: decreased diversity but increased transmission risks. Ecol Appl. 2019, 29: e01886. https://doi.org/10.1002/eap.1886 PMID: 30986339

45. Chaisiri K, Tanganuchitcharnchai A, Kritiyakan A, Thinphovong C. Tanita M, Morand S, et al. Risk factors analysis for neglected human rickettsioses in rural communities in Nan province, Thailand: A community-based observational study along a landscape gradient. PLoS Negl Trop Dis. 2022; 16: e0010256. https://doi.org/10.1371/journal.pntd.0010256 PMID: 35320277

46. Mahuzier C, Morand S, Chaisiri K, De Rouw A, Soulileuth B, Thinphovong C, et al. Random Forest land cover classifications of Sentinel satellite images in 2019, Saen Thong, Thailand. DataSuds. https://doi.org/10.23708/GENR6J.

47. Morrison K. BaTML Factsheet: Harp trapping guidance notes for bat workers. BaTML Publications. 2005. https://batability.co.uk/wp-content/uploads/2018/11/BaTML_Factsheet_Organising_Bat_WalksTalks-1.pdf.

48. AVMA. Guidelines for the Euthanasia of Animals: 2020 Edition. https://www.avma.org/resources-tools/avma-policies/avma-guidelines-euthanasia-animals.

49. Navarro KL, Huss M, Smith JC, Sharp P, Marx JO, Pacharinsak C. Mouse anesthesia: The art and science. ILAR J. 2021; 238–273. https://doi.org/10.1093/ilar/ilab016 PMID: 34180990

50. IDEXX SARS-CoV-2 (COVID-19) RealPCR Test.: //www.idexx.com/en/veterinary/reference-laboratories/idexx-realpcr-tests/idexx-sars-2-covid-19-realpcr-test/.

51. Gijs M, Veugen JMJ, Wolffs PFG, Savelkoul PHM, Tas J, van Bussel BCT, et al. In-Depth investigation of conjunctival swabs and tear fluid of symptomatic COVID-19 patients, an observational cohort study. Transl Vis Sci Technol. 2021; 10: 32. https://doi.org/10.1167/tvst.10.12.32 PMID: 34673905

52. Wu Z, Han Y, Liu B, Li H, Zhu G, Latinne A, et al. (2021) Decoding the RNA viromes in rodent lungs provides new insight into the origin and evolutionary patterns of rodent-borne pathogens in Mainland Southeast Asia. Microbiome. 9: 18. https://doi.org/10.1186/s40168-020-00965-z PMID: 33478588

53. Niedballa J, Sollmann R, Courtiol A, Wilting A. camtrapR: an R package for efficient camera trap data management. Methods Ecol Evol. 2016; 7: 1457–1462. https://doi.org/10.1111/2041-210X.12600.

54. Wells K, Morand S, Wardeh M, Baylis M. Distinct spread of DNA and RNA viruses among mammals amid prominent role of domestic species. Glob Ecol Biogeo. 2020a; 29: 470–481. https://doi.org/10.1111/geb.13045 PMID: 32336945

55. Wells K, Morand S, Wardeh M, Baylis M. Data from: Distinct spread of DNA and RNA viruses among mammals amid prominent role of domestic species, 2020b; v3, Dryad, Dataset. Dryad Digital Repository: https://doi.org/10.5061/dryad.p2ngf1vmg.

56. Morand S, Chaisiri K, Kritiyakan A, Kumlert R. Disease Ecology of Rickettsial Species: A Data Science Approach. Trop Med Infect Dis. 2020; 5: 56. https://doi.org/10.3390/tropicalmed5020064.

57. Bordes F, Morand S, Pilosof S, Claude J, Krasnov B, Cosson JF, et al. Habitat fragmentation alters the properties of a host-parasite network: Rodents and their helminths in South-East Asia. J Anim Ecol. 2015; 84: 1253–1263. https://doi.org/10.1111/1365-2656.12368 PMID: 25777342
58. Dormann CF, Fruend J, Bluethgen N, Gruber B. Indices, graphs and null models: analyzing bipartite ecological networks. The Open Ecology Journal. 2009; 2: 7–24.

59. Dormann CF, Gruber B, Fruend J. Introducing the bipartite Package: Analysing Ecological Networks. R news. 2008; 8/2: 8–11.

60. R Development Core Team. The R Project for Statistical Computing, R version 4.2.0 2022 Available online: https://www.r-project.org (accessed on 2023-08-01).

61. Dormann C.F. & Strau R. 2013 Detecting modules in quantitative bipartite networks: the QuaBiMo algorithm. ArXiv. 1304, 3218.

62. Opsahl T. Structure and Evolution of Weighted Networks. University of London (Queen Mary College): London, UK. 2009; pp. 104–122.

63. Csardi G, Nepusz T. The igraph software package for complex network research. Int J Complex Syst. 2006; 1695.

64. Blondel VD, Guillaume JL, Lambiotte R, Lefebvre Étienne. Fast unfolding of community hierarchies in large networks. J Stat Mech. 2008; 10: 10008.

65. Morand S, McIntyre KM, Baylis M (2014) Domesticated animals and human infectious diseases of zoo-notic origins: domestication time matters. Infect Genet Evol. 24:76–81. https://doi.org/10.1016/j.meegid.2014.02.013 PMID: 24642136