- 1Division of Natural Sciences and Mathematics, Transylvania University, Lexington, KY, United States
- 2Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, University of Kent, Canterbury, United Kingdom
- 3Instituto Leônidas e Maria Deane – Fiocruz Amazônia, Manaus, Amazonas, Brazil
- 4MRC Unit The Gambia at London School of Hygiene and Tropical Medicine, Fajara, Gambia
- 5EcoHealth Alliance, New York, NY, United States
- 6Departamento de Recursos Naturais, Faculdade de Ciências Agronomicas, Universidade Estadual Paulista, Botucatu, Brazil
- 7Departamento de Medicina Veterinária Preventiva e Saúde Animal da Faculdade de Medicina Veterinária e Zootecnia da Universidade de São Paulo, São Paulo, Brazil
- 8Department of Environmental Science and Policy, George Mason University, Fairfax VA, United States
- 9Departamento de Microbiologia, Instituto de Ciências Biomédicas-II, Universidade de São Paulo, São Paulo, Brazil
Introduction: Bats are critical to maintaining healthy ecosystems and many species are threatened primarily due to global habitat loss. Bats are also important hosts of a range of viruses, several of which have had significant impacts on global public health. The emergence of these viruses has been associated with land-use change and decreased host species richness. Yet, few studies have assessed how bat communities and the viruses they host alter with land-use change, particularly in highly biodiverse sites.
Methods: In this study, we investigate the effects of deforestation on bat host species richness and diversity, and viral prevalence and richness across five forested sites and three nearby deforested sites in the interior Atlantic Forest of southern Brazil. Nested-PCR and qPCR were used to amplify and detect viral genetic sequence from six viral families (corona-, adeno-, herpes-, hanta-, paramyxo-, and astro-viridae) in 944 blood, saliva and rectal samples collected from 335 bats.
Results: We found that deforested sites had a less diverse bat community than forested sites, but higher viral prevalence and richness after controlling for confounding factors. Viral detection was more likely in juvenile males located in deforested sites. Interestingly, we also found a significant effect of host bat species on viral prevalence indicating that viral taxa were detected more frequently in some species than others. In particular, viruses from the Coronaviridae family were detected more frequently in generalist species compared to specialist species.
Discussion: Our findings suggest that deforestation may drive changes in the ecosystem which reduce bat host diversity while increasing the abundance of generalist species which host a wider range of viruses.
Introduction
Emerging viruses with wildlife origins are a significant threat to global health (e.g. Ebolaviruses, SARS and MERS coronaviruses) (Jones et al., 2008). Analyses of recent emerging infectious disease (EID) events show that anthropogenic changes including land-use change (e.g. habitat degradation, deforestation, forest fragmentation), intensification of food production, and global trade and travel are key factors in disease emergence (Loh et al., 2015; Allen et al., 2017; Rulli et al., 2017; Reaser et al., 2022). Further, nearly one-third of all EIDs, and a higher proportion of zoonoses, are associated with land-use change specifically (Loh et al., 2015). This suggests that increasing and/or novel interactions among hosts, vectors and pathogens following land-use change are significant contributors to disease emergence.
In tropical and subtropical environments, the pace of land-use change is unprecedented and continues to increase globally as demand for natural resources grows (Song et al., 2018). Bats are globally threatened, with 15% of bat species being listed as threatened or vulnerable, and habitat loss in the tropics is a major driver of population declines (Frick et al., 2020). Land-use change has also been associated with the emergence of many recent zoonotic diseases (Gibb et al., 2020). Yet, the relationship between land-use change and disease emergence is poorly understood. Recent studies have hypothesized that land-use change may increase the risk of disease emergence through more frequent human-animal interactions, or by influencing pathogen diversity, either directly by changing pathogen prevalence and/or diversity, or indirectly via impacts on host assemblages (Bradley et al., 2008; Vittor et al., 2009; Murray and Daszak, 2013; Rulli et al., 2017). However, mechanistic studies have tended to focus on how abundance and prevalence of specific pathogens, or their vectors and hosts, vary over the landscape (Ostfeld and Keesing, 2000; LoGiudice et al., 2003; Kilpatrick et al., 2006a; Kilpatrick et al., 2006b; Bradley et al., 2008; Vittor et al., 2009). Others have used meta-analyses to try to identify generality and mechanisms involved (Salkeld et al., 2013; Gottdenker et al., 2014; Civitello et al., 2015). Few empirical studies have taken a community approach to examine how viral assemblages in host communities vary with land-use change.
In this study, we investigate the effects of deforestation on bat host abundance and diversity, and viral prevalence and richness. We work with bats because they are diverse, abundant, and geographically widespread (Rex et al., 2008), comprising species from nearly every trophic level, with wide differences in their dispersal abilities (Kingston, 2010). Further, some of their life history traits and characteristics (e.g. diet, ability to fly, torpor and hibernations, and roosting behaviors) make them suitable hosts of viruses and other pathogens (Calisher et al., 2006) and many bat species are strongly impacted by land-use changes. Bats are important hosts of pathogens that have had significant impact on public health (e.g. Ebola, SARS, MERS, rabies). They also harbor the highest proportion of zoonotic viruses of any mammal order (Jones et al., 2008; Olival et al., 2017), as well as significant emerging diseases of people, livestock and wildlife. Finally, while no bat viruses have emerged from the Atlantic Forest to our knowledge, we chose this region as our study site because of the high biodiversity it contains and the large-scale deforestation it has undergone. Our study focuses on three questions: (1) Does bat abundance and diversity differ in forested versus non-forested areas? (2) Does viral prevalence differ between bat communities in forested versus deforested areas? (3) What biological and ecological factors determine the likelihood of viral detection?
Materials and methods
Ethical statement
This study was carried out with animal handling permits issued from the Brazilian Ministry of the Environment (#33078-4). Animal handling ethics approval was provided by the University of California, Davis (#16048). Bat handling followed strict personal protection and biosafety requirements and short capture times to minimize stress on individual animals. All captured individuals were released at the point of capture.
Study site
Morro do Diabo State Park (municipality of Teodoro Sampaio, São Paulo state, Brazil, Figure 1) is located in the Pontal do Paranapanema region and contains the largest preserved area of interior Atlantic Forest in Sao Paulo State. The park covers an area of 33,845 ha (Durigan and Franco, 2006) and is comprised of mesophytic semideciduous forest and a small area of Cerrado (savanna-like vegetation). The climate is characterized as subtropical, with dry winters and wet summers. Mean annual temperature is 22°C, and annual rainfall ranges between 1100 and 1300 mm (de Faria and Pires, 2006). The matrix around the park is comprised of 63 small properties of agrarian reform settlements, as well as pasture (~60%) and agriculture (~15%), and forest fragments ranging from 2 to 2000 ha in area, most of which are privately owned (Uezu and Metzger, 2011). The forested study sites were chosen to control for similar characteristics including elevation, vegetation structure and rainfall. We sampled bats and viruses at five intact forested sites (>200ha) and three nearby deforested sites, located 3-5km away and defined as areas where more than 20% of forest cover has been removed and converted from the original forest to agrarian reform settlements.
Figure 1 The study area, in and around Morro do Diabo State Park, located in the Pontal do Paranapanema region, Sao Paulo state. Bat surveys (yellow and red circles) were undertaken across forested (n=5) and deforested sites (n=3).
Bat capture and sample collection
Bats were sampled during April to November of 2014. At each site we sampled a 100m x 100m grid using eight horizontal mist nets (9m x 3m), one canopy mist net (6m x 3m) and one harp trap (1.5m x 1.5m). At least 150 individuals were collected from deforested sites and from forested sites each. Additional sampling effort was required in the forested areas due to lower catch rates. Bats were captured for a period of five consecutive nights at each site, totaling 2040 m2/270 hours capture effort across all sites. Mist nets were opened at sunset and remained open for six hours. Nets were checked at 30-minute intervals and bats processed immediately.
Samples were collected from bats with no clinical or neurological symptoms and in good body condition, defined as mass divided by forearm length which has been validated in temperate bats (Wilkinson and Brunet-Rossinni, 2009). All animals were released immediately following processing. Blood, saliva, and rectal swabs were collected from each captured animal, with feces and urine opportunistically collected. All samples were placed in cryovials containing 200 ml of Viral Transportation Media (VTM) and stored in liquid nitrogen in a dry shipper while in the field, then transferred to -80C freezers at the Institute of Biomedical Sciences at the University of Sao Paulo. External morphological measurements (including forearm/radius length, body length, head length) were collected by a bat taxonomist to assist in species identification. Several bat identification keys from the region were also used for reference (Reis et al., 2013; Reis et al., 2017). Sex and age were determined by the presence of scrotal testes and well developed teats (Wilkinson and Brunet-Rossinni, 2009). Before release, each individual was marked with a non-toxic pen to determine the rate of within-trip recapture. This was used to ensure that the same bat was not re-sampled within sampling trips.
Viral detection
Total nucleic acid was extracted from all samples using the EasyMag (bioMérieux, Inc.) platform, and cDNA synthesis performed using SuperScript III first-strand synthesis supermix (Invitrogen), all according to the manufacturer’s instructions. Viral discovery was performed using nested-PCR assays targeting coronaviruses (Quan et al., 2010), astroviruses (Atkins et al., 2009), paramyxoviruses (Tong et al., 2008), and herpesviruses (VanDevanter et al., 1996), while real-time PCR was used to target hantaviruses (Araujo et al., 2011) described below (Box 1). PCR results were visualized on a 2% agarose gel and Sanger sequencing was performed using ABI3100 (Applied Biosystems) equipment and BigDye Terminator v3.1 Cycle Sequencing Kit at the Institute of Biomedical Sciences II at the University of Sao Paulo. Sequences were analysed and edited using Geneious (version 6.0.3). Sequences were aligned with ClustalW and MUSCLE, and phylogenetic trees (see Text S1) constructed with neighbor-joining (p-distance, pairwise deletion, 1,000 bootstraps), maximum-likelihood (1,000 bootstraps), and Bayesian (GTR+I - Mr Bayes) algorithms. In Mr. Bayes, we discarded the first 25% of trees as burn-in, and used the remaining trees to estimate the posterior probability value (PP) of 0.7. The chains ran for 2,000,000 cycles (mcmc ngen = 2,000,000). Trees were reconstructed with unconstrained branch lengths and unrooted. In MEGA 7 (macOS available in: https://www.megasoftware.net/) we used Maximum Likelihood with heuristic search and GTR+gamma+I algorithm. For the ML tree, we conducted 1,000 fast bootstrap ML replicates to assess the support values of internal nodes and visualized the trees in FigTree software version 1.4.4 with Midpoint Root (available in: http://tree.bio.ed.ac.uk/software/figtree/)(Supplementary Figures 1–5). Sequences were segregated into discrete viruses, defined as a viral species, based on distinct monophyletic clustering following Anthony et al. (2013).
Data analysis
Statistical analysis was performed using R 3.5.1, with ggplot2 for graphing. To compare estimated bat species diversity between forested and deforested sites, we calculated abundance-based diversity profiles with Hill numbers (effective number of species) using the iNEXT package based on the parameter q (Chao and Jost, 2015). This parameter controls the relative emphasis placed on rare or common species. In addition to providing information on species richness, this diversity profile estimator also accounts for species abundances to differing degrees. With increasing order q, the weight of dominant species increases in the calculation of species diversity. We used three widely used species diversity measures: Species richness (number of observed species; q=0), Shannon diversity (number of typical species; q=1) and Simpson diversity (number of most common species; q=2). We then applied a bootstrap method (1,000 bootstraps) using observed detections to obtain approximate variances of the proposed profiles and to construct the associated confidence intervals. These estimations take into account the effect of undetected species in samples. Estimated viral diversity could not be explored using these methods due to limited sample sizes. However we compared viral species richness and overall viral prevalence across treatments using a Fisher’s Exact Test. To account for the uneven number of captures per bat species, we used Bartels rank test of randomness to determine whether viruses were randomly distributed among bat host species by examining whether viral prevalence significantly differed among species. Due to low detection rates in other viral families, our analysis was limited to the coronaviridae family.
We use a Generalized Linear model (GLM) of viral detection with a logit link function. We use “viral detection” as the response variable in our model based on the presence or absence of a viral detection for each individual bat. After testing for collinearity among the response variables, no variables were excluded based on their variance inflation factor (VIF) scores. Seven variables were selected for the final analysis. Definitions of the variables used are given in Table 1. In a “stepwise backwards-selection”, factors were eliminated from the full model in an iterative process based on the Akaike information criterion (AIC) (Akaike, 1973) with the stepAIC function of the MASS package (Venables and Ripley, 2002) in the statistical software R 3.5.1.
Results
Bat diversity
We recorded 18 bat species from three families (Phyllostomidae, Molossidae, and Vespertilionidae) and five dietary guilds (frugivorous, insectivorous, nectarivorous, sanguivorous, omnivourous) from 335 mist-net captures. No bats were captured using the vertical canopy net or harp trap. After accounting for sampling effort, capture rates were similar between forested (n = 163 captures) and deforested (n = 172 captures) sites (Paired t-test, t = 1.883, p = 0.081). Bat species richness in the deforested sites (n=11 species) was slightly higher than the forested sites (n = 9 species); however, this difference was not significant as indicated by the empirical diversity profiles that show overlap between the 95% confidence intervals at q = 0 (Figure 2). In contrast, at q>1 (i.e. measures of diversity that incorporate abundance information) the forested sites were found to be more diverse than deforested sites. When correcting for the bias introduced by the non-detection of species in the samples, bat diversity was reduced in deforested sites; species richness was slightly higher in forested areas (n=15 species) compared to deforested areas (n = 11 species). However, for q > 1.25, this difference in community diversity was statistically significant, as reflected by the two non-overlapped confidence intervals.
Figure 2 Estimated diversity profiles for bat species data in forested (green line) and deforested (red line) sites for q between 0 and 3 with 95% confidence interval (shaded areas based on a bootstrap method of 1000 replications). The numbers show the estimated diversities for q = 0, 1, 2 and 3.
Viral prevalence and richness
Overall, a total of 22 individual bats from three families (Phyllostomidae, Molossidae, Vespertilionidae) were positive for 13 viral species in the following viral families: Astroviridae, Coronaviridae, Hantaviridae, Herpesviridae, and Paramyxoviridae, with a combined viral prevalence of 6.6% (22/335) (Table 2). None of the samples were positive for adenoviruses, despite previous studies documenting their presence in other bat species (Jánoska et al., 2011; Van Vuren et al., 2018). Only one individual bat yielded more than one viral species - a coinfection by a coronavirus and herpesvirus was found in Artibeus planirostris. Viral species were not evenly distributed among bat species, with all detected viruses coming from just five of 18 sampled bat species (Figure 3). After accounting for the number of captures per bat species by looking at viral prevalence as opposed to number of positive detections, we found a significant effect of host bat species on viral prevalence indicating that viral taxa were detected more frequently in some species than others. In particular, viruses from the Coronaviridae family were detected more frequently in generalist species compared to specialist species (P<0.01, Bartel’s Rank Test). Viral prevalence also differed among viral families; Coronaviridae had the highest prevalence of 3.6%, followed by Astroviridae (1.2%), Paramyxoviridae (0.6%), Herpesviridae (0.9%) and Hantaviridae (0.3%).
Table 2 Total captures of bat species and total viral detections in forested and deforested habitat in the Interior Atlantic Forest.
Figure 3 Overall viral prevalence by bat species with sample sizes. Error bars represent the 95% Clopper-Pearson binomial confidence intervals. Species codes: ARTLIT, Artibeus lituratus; ARTPLA, Artibeus planirostris; CARPER, Carollia perspicillata; PHYHAS, Phyllostomus hastatus; STULIL, Sturnira lillium.
Treatment (forested vs. deforested), sex and age were the only significant predictors of overall viral detection (presence/absence) (P < 0.05, df = 1) (Table 3). This result was supported by the logistic regression model with the lowest AIC value (Table 4), which demonstrates that the odds of a positive viral detection decreases in forested habitat.
With all viral families combined, viral prevalence in deforested sites (9.3%) was significantly higher than in forested sites (3.68%) (P<0.05, Fisher’s Exact Test). Deforested sites also had higher viral richness (n = 13 unique viral taxa) compared to forested sites (n = 2).
Discussion
In the Atlantic Forest of Brazil, higher bat host diversity is not associated with higher viral prevalence or richness. Despite lower bat host species richness in deforested areas, viral richness and prevalence is significantly higher. This result does not appear to be associated with the abundance of bat hosts, which was not significantly different (based on mist-net capture frequency) in deforested versus forested areas.
Bats are known to harbor a wide diversity of viruses, and have received growing attention due to their role in the emergence of several recent infectious disease outbreaks (e.g. Severe Acute Respiratory Syndrome (SARS), Middle East Respiratory Syndrome, Nipah virus) (Quan et al., 2010; Hu et al., 2015; Epstein et al., 2020). While studies exploring viral diversity in bat host species have increased, few studies have assessed how bat communities and the viruses they host alter with land-use change, particularly in highly biodiverse sites. Overall, this study identified 13 unique viral taxa from four viral families known to infect humans. We found that different viral families were not evenly distributed within different bat host species and between habitats (forested vs deforested). Specifically, viruses from Coronoviridae were primarily found in species considered to be generalists, including Artibeus planirostris, Carollia perspicillata, Artibeus literatus, and Sturnira lillium. We suggest that such differences in virus prevalence could be related to viral ecology (i.e., their ability to infect host cells and to persist and replicate) and to the ecology and behavior of the bat hosts in a given habitat. Specifically, we found that viral detection is more likely in juvenile, male bats. Indeed, previous studies have shown that the behavior between many species of young male and female bats differ considerably, with young males immediately leaving the maternity roost once they are weaned, while females continue to forage with their mothers. This difference in behavior could result in younger males having a greater frequency of contacts with new host species or with shared food resources that increase their exposure to potential pathogens. For example, younger vampire bats appear to have higher exposure to pathogens such as rabies virus because younger male bats are more exploratory and are more likely to feed on novel hosts (Carter et al., 2018).
Previous studies of bats have demonstrated that even moderate forest disturbance can result in an increase of certain generalist species that can successfully adapt to human-modified landscapes (Delaval and Charles-Dominique, 2006; Meyer and Kalko, 2008). The strategies they employ such as greater dispersal ability, and the ability to exploit a variety of resources, allow these species to tolerate a wide range of habitats, leading to higher colonization rates throughout human-modified landscapes. In our study area, A. planirostris and A. fimbriatus were the two species most commonly captured in deforested sites. Both species are large-bodied frugivores, which feed heavily on figs in the canopy (Handley et al., 1991). In many tropical landscapes, figs are not regularly available throughout the year, thus Artibeus species are more likely to occupy disturbed landscapes which provide a variety of food resources (Gorresen and Willig, 2004). Artibeus spp. bats accounted for 73% (n=16/22) of all viruses detected, after accounting for the number of captures per bat species.
Here, we show that deforested sites support higher viral richness despite lower bat species richness. While we did not measure disease risk directly, we hypothesize that humans living close to forest edges disturbed by deforestation may be particularly exposed to zoonotic infections not only because of the higher likelihood for humans to be in contact with disease reservoirs, but also because of the higher viral richness found in deforested areas. Yet, previous studies examining the link between land-use change and disease in have been equivocal (Randolph and Dobson, 2012; Salkeld et al., 2013; Civitello et al., 2015; Rulli et al., 2017). Some studies of single-pathogens (e.g. West Nile virus, Hantavirus, the Lyme disease pathogen Borellia burgdorferi) in multi-host systems have found that higher pathogen prevalence is associated with decreased continuous forest area (LoGiudice et al., 2003; Suzán et al., 2008; Kilpatrick, 2011). Results from a recent meta-analysis from studies in Southeast Asia shows that people who live or work on agricultural land are more likely to be infected with zoonotic diseases (Shah et al., 2019), and in West and Central Africa, previous research shows that the index cases of Ebola virus outbreaks (i.e. spillover cases from wildlife reservoirs) occurred mostly in areas of forest fragmentation and deforestation (Rulli et al., 2017). Further, urbanized and agricultural areas that have undergone deforestation have been associated with higher rates of disease transmission of West Nile Virus in the United States, increased risk of malaria in Peru (Vittor et al., 2009), Leishmaniasis in Costa Rica (Wijeyaratne et al., 1994) and hantavirus in Panama (Suzán et al., 2008), in part because of changes in host and vector abundance in human-modified areas.
Our study provides further evidence from a multi-pathogen, multi-host species system that deforestation can increase viral prevalence and richness in bat hosts. However, studies of Plasmodium infections in Australia (Laurance et al., 2013), Cameroon (Chasar et al., 2009), and Brazil (Ribeiro et al., 2005) found a positive correlation between continuous forest area and pathogen prevalence. In Sabah, Seltmann et al., 2017 found that reduced body mass in bats in logged forests was associated with chronic stress and impaired health status for some species of bats. Interestingly, this did not translate into an increase in coronavirus and astrovirus detection rates among more disturbed sites (Seltmann et al., 2017), perhaps due to the extent of disturbance. Unlike our system, which is more than 30 years post-fragmentation and fully converted, the Sabah study sites are still ongoing active deforestation and fragmentation, which may result in delays in species’ responses.
Our study examines some of the complexities in the relationship among deforestation, viral prevalence and host and viral community assemblages by addressing how viral richness and prevalence in bat hosts varies with land-use change. Our findings suggest that deforestation can increase the abundance of generalist species that, in our case, host the majority of viruses detected. From a theoretical point of view, the dilution effect hypothesis explores how the decrease of biodiversity may increase the amplification of zoonotic diseases. It suggests that high species diversity in a community can reduce infectious disease risk, provided that hosts differ in competency for transmitting a pathogen (Schmidt and Ostfeld, 2001). While this study does not test the ‘dilution effect’ as laid out for Lyme disease and other single pathogen systems (LoGiudice et al., 2003; Kilpatrick et al., 2005), these findings provide further evidence that anthropogenic land use change can in some cases, lead to increased abundance of reservoirs that harbor a higher diversity and prevalence of potential pathogens. As pressures on the environment continue to grow, further research is needed on viral and host ecology and how they are structured across varying landscapes.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics statement
This study was carried out with animal handling permits from the Brazilian Ministry of the Environment (#33078-4). Animal handling ethics approval was provided by the University of California, Davis (#16048) and by the Ethics Committee of the Faculty of Veterinary Medicine and Zootechnics at the University of Sao Paulo.
Author contributions
All authors contributed to the article and approved the submitted version.
Funding
This study was made possible by the United States Agency for International Development (USAID) Emerging Pandemic Threats PREDICT program (cooperative agreement no., GHN-A-OO-09-00010-00). The contents are the responsibility of the authors and do not necessarily reflect the views of USAID or the United States Government. This work was also supported by The American Association of University Women (Dissertation Fellowship) and by the International Development Research Centre (project no. 106150-001).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.921950/full#supplementary-material
Supplementary Figure 1 | Paramyxovirus Maximum Likelihood Phylogenetic Tree. Genetic analysis of 558 nucleotide partial L gene. Tree reconstructed by MEGA7 with heuristic search, Neighbor-Joining “NJ” algorithm and Model GTR+gamma+I. The principal node values superior to 70% represent 1,000 bootstrap replicates.
Supplementary Figure 2 | Coronavirus Maximum Likelihood Phylogenetic Tree. Genetic analysis of 394 nucleotide partial RdRp gene. Tree reconstructed by MEGA7 with heuristic search, Neighbor-Joining “NJ” algorithm and Model GTR+gamma+I. The principal node values superior to 70% represent 1,000 bootstrap replicates.
Supplementary Figure 3 | Astrovirus Maximum Likelihood Phylogenetic Tree. Genetic analysis of 369 nucleotide partial RdRp gene. Tree reconstructed by MEGA7 with heuristic search, Neighbor-Joining “NJ” algorithm and Model GTR+gamma+I. The principal node values superior to 70% represent 1,000 bootstrap replicates.
Supplementary Figure 4 | Herpesvirus Maximum Likelihood Phylogenetic Tree. Genetic analysis of 189 nucleotide partial Polymerase (pol) gene. Tree reconstructed by MEGA7 with heuristic search, Neighbor-Joining “NJ” algorithm and Model GTR+gamma+I. The principal node values superior to 70% represent 1,000 bootstrap replicates.
Supplementary Figure 5 | Hantavirus Maximum Likelihood Phylogenetic Tree. Genetic analysis of 89 nucleotide partial S segment. Tree reconstructed by MEGA7 with heuristic search, Neighbor-Joining “NJ” algorithm and Model GTR+gamma+I. The principal node values superior to 70% represent 1,000 bootstrap replicates.
References
Akaike, H. (1973). Information theory and an extension of the maximum likelihood principle (Budapest, Hungary: Akadémiai Kiadó).
Allen, T., Murray, K. A., Zambrana-Torrelio, C., Morse, S. S., Rondinini, C., Di Marco, M., et al. (2017). Global hotspots and correlates of emerging zoonotic diseases. Nat. Commun. 8 (1), 1124. doi: 10.1038/s41467-017-00923-8
Anthony, S. J., Epstein, J. H., Murray, K. A., Navarrete-Macias, I., Zambrana-Torrelio, C. M., Solovyov, A., et al. (2013). A strategy to estimate unknown viral diversity in mammals. Mbio 4 (5):1028–1038 doi: 10.1128/mBio.00598-13
Araujo, J., Pereira, A., Nardi, M., Henriques, D., Lautenschalager, D., Dutra, L., et al. (2011). Detection of hantaviruses in Brazilian rodents by SYBR-green-based real-time RT-PCR. Arch. virology. 156 (7), 1269–1274. doi: 10.1007/s00705-011-0968-1
Atkins, A., Wellehan, J. F., Childress, A. L., Archer, L. L., Fraser, W. A., Citino, S. B. (2009). Characterization of an outbreak of astroviral diarrhea in a group of cheetahs (Acinonyx jubatus). Veterinary Microbiol. 136 (1), 160–165. doi: 10.1016/j.vetmic.2008.10.035
Bradley, C. A., Gibbs, S. E., Altizer, S. (2008). Urban land use predicts West Nile virus exposure in songbirds. Ecol. Applications. 18 (5), 1083–1092. doi: 10.1890/07-0822.1
Calisher, C. H., Childs, J. E., Field, H. E., Holmes, K. V., Schountz, T, et al (2006). Bats: important reservoir hosts of emerging viruses. Clinical microbiology reviews 19.3, 531–545.
Carter, G. G., Forss, S., Page, R. A., Ratcliffe, J. M. (2018). Younger vampire bats (Desmodus rotundus) are more likely than adults to explore novel objects. PloS One 13 (5), e0196889. doi: 10.1371/journal.pone.0196889
Chao, A., Jost, L. (2015). Estimating diversity and entropy profiles via discovery rates of new species. Methods Ecol. Evolution. 6 (8), 873–882. doi: 10.1111/2041-210X.12349
Chasar, A., Loiseau, C., Valkiūnas, G., Iezhova, T., Smith, T. B., Sehgal, R. N. (2009). Prevalence and diversity patterns of avian blood parasites in degraded African rainforest habitats. Mol. Ecology. 18 (19), 4121–4133. doi: 10.1111/j.1365-294X.2009.04346.x
Civitello, D. J., Cohen, J., Fatima, H., Halstead, N. T., Liriano, J., McMahon, T. A., et al. (2015). Biodiversity inhibits parasites: broad evidence for the dilution effect. Proc. Natl. Acad. Sci. 112 (28), 8667–8671. doi: 10.1073/pnas.1506279112
de Faria, H. H., Pires, A. (2006). Parque estadual do morro do diabo: plano de manejo (Santa Cruz do Rio Pardo: Editora Viena), 311.
Delaval, M., Charles-Dominique, P. (2006). Edge effects on frugivorous and nectarivorous bat communities in a neotropical primary forest in French Guiana. Revue d'Ecologie, Terre et Vie 61.4, 343–352.
Durigan, G., Franco, G. (2006). “Vegetação,” in Parque estadual do morro do diabo: plano de manejo (Santa Cruz do Rio Pardo: Editora Viena), 111–118.
Epstein, J. H., Anthony, S. J., Islam, A., Kilpatrick, A. M., Ali Khan, S., Balkey, M. D., et al. (2020). Nipah virus dynamics in bats and implications for spillover to humans. Proc. Natl. Acad. Sci. 117 (46), 29190–29201. doi: 10.1073/pnas.2000429117
Frick, W. F., Kingston, T., Flanders, J. (2020). A review of the major threats and challenges to global bat conservation. Ann. New York Acad. Sci. 1469 (1), 5–25. doi: 10.1111/nyas.14045
Gibb, R., Redding, D. W., Chin, K. Q., Donnelly, C. A., Blackburn, T. M., Newbold, T., et al. (2020). Zoonotic host diversity increases in human-dominated ecosystems. Nature 584 (7821), 398–402. doi: 10.1038/s41586-020-2562-8
Gorresen, P. M., Willig, M. R. (2004). Landscape responses of bats to habitat fragmentation in Atlantic forest of Paraguay. J. Mammalogy. 85 (4), 688–697. doi: 10.1644/BWG-125
Gottdenker, N. L., Streicker, D. G., Faust, C. L., Carroll, C. (2014). Anthropogenic land use change and infectious diseases: a review of the evidence. Ecohealth 11 (4), 619–632. doi: 10.1007/s10393-014-0941-z
Handley, C. O., Jr., Wilson, D. E., Gardner, A. L. (1991). Demography and natural history of the common fruit bat, artibeus jamaicensis, on barro Colorado island, Panama. doi: 10.5479/si.00810282.511
Hu, B., Ge, X., Wang, L.-F., Shi, Z. (2015). Bat origin of human coronaviruses. Virol. J. 12 (1), 1–10. doi: 10.1186/s12985-015-0422-1
Jánoska, M., Vidovszky, M., Molnár, V., Liptovszky, M., Harrach, B., Benkő, M. (2011). Novel adenoviruses and herpesviruses detected in bats. Veterinary J. 189 (1), 118–121. doi: 10.1016/j.tvjl.2010.06.020
Jones, K. E., Patel, N., Levy, M., Storeygard, A., Balk, D., Gittleman, J. L., et al. (2008). Global trends in emerging infectious diseases. Nature 451, 990–993. doi: 10.1038/nature06536
Kilpatrick, A. M. (2011). Globalization, land use, and the invasion of West Nile virus. Science 334 (6054), 323–327. doi: 10.1126/science.1201010
Kilpatrick, A. M., Daszak, P., Jones, M. J., Marra, P. P., Kramer, L. D. (2006a). Host heterogeneity dominates West Nile virus transmission. Proc. R. Soc. London B: Biol. Sci. 273 (1599), 2327–2333.
Kilpatrick, A. M., Daszak, P., Kramer, L., Marra, P. P., Dobson, A. P., Hudson, P. J., et al. (2006b). Predicting the transmission of West Nile virus. Am. J. Trop. Med. Hygiene. 75 (5), 139.
Kilpatrick, A. M., Jones, M., Kramer, L. D., Marra, P. P., Daszak, P. (2005). West Nile Virus vector ecology across an urbanization gradient. Am. J. Trop. Med. Hygiene. 73 (6), 307–308. doi: 10.3201/eid1103.040364
Kingston, T. (2010). Research priorities for bat conservation in southeast Asia: a consensus approach. Biodiversity Conserv. 19 (2), 471–484. doi: 10.1007/s10531-008-9458-5
Laurance, S. G., Jones, D., Westcott, D., Mckeown, A., Harrington, G., Hilbert, D. W. (2013). Habitat fragmentation and ecological traits influence the prevalence of avian blood parasites in a tropical rainforest landscape. PloS One 8 (10), e76227. doi: 10.1371/journal.pone.0076227
LoGiudice, K., Ostfeld, R. S., Schmidt, K. A., Keesing, F. (2003). The ecology of infectious disease: effects of host diversity and community composition on Lyme disease risk. Proc. Natl. Acad. Sci. 100 (2), 567–571. doi: 10.1073/pnas.0233733100
Loh, E. H., Zambrana-Torrelio, C., Olival, K. J., Bogich, T. L., Johnson, C. K., Mazet, J. A., et al. (2015). Targeting transmission pathways for emerging zoonotic disease surveillance and control. Vector-Borne Zoonotic Diseases. 15 (7), 432–437. doi: 10.1089/vbz.2013.1563
Meyer, C. F., Kalko, E. K. (2008). Assemblage-level responses of phyllostomid bats to tropical forest fragmentation: land-bridge islands as a model system. J. Biogeography. 35 (9), 1711–1726. doi: 10.1111/j.1365-2699.2008.01916.x
Murray, K. A., Daszak, P. (2013). Human ecology in pathogenic landscapes: two hypotheses on how land use change drives viral emergence. Curr. Opin. virology. 3 (1), 79–83. doi: 10.1016/j.coviro.2013.01.006
Olival, K. J., Hosseini, P. R., Zambrana-Torrelio, C., Ross, N., Bogich, T. L., Daszak, P. (2017). Host and viral traits predict zoonotic spillover from mammals. Nature 546, 646–650. doi: 10.1038/nature22975
Ostfeld, R. S., Keesing, F. (2000). Biodiversity and disease risk: the case of Lyme disease. Conserv. Biol. 14 (3), 722–728. doi: 10.1046/j.1523-1739.2000.99014.x
Quan, P.-L., Firth, C., Street, C., Henriquez, J. A., Petrosov, A., Tashmukhamedova, A., et al. (2010). Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria. MBio 1 (4), e00208–e00210. doi: 10.1128/mBio.00208-10
Randolph, S. E., Dobson, A. (2012). Pangloss revisited: a critique of the dilution effect and the biodiversity-buffers-disease paradigm. Parasitology 139 (7), 847–863. doi: 10.1017/S0031182012000200
Reaser, J. K., Hunt, B. E., Ruiz-Aravena, M., Tabor, G. M., Patz, J. A., Becker, D. J., et al. (2022). Fostering landscape immunity to protect human health: A science-based rationale for shifting conservation policy paradigms. Conserv. Lett., e12869. doi: 10.1111/conl.12869
Reis, N. R., Fregonezi, M. N., Peracchi, A. L., Shibatta, O. A. (2013). Morcegos do brasil: guia de campo (1st edn). (Rio de Janeiro: Technical Books Editora), 225 pp.
Reis, N. R., Peracchi, A. L., Batista, C. B., de Lima, I. P., Pereira, A. D. (2017). História natural dos morcegos brasileiros: chave de identificação de espécies (Rio de Janeiro: Technical Books Editora), 416.
Rex, K., Kelm, D. H., Wiesner, K., Kunz, T. H., Voigt, C. C. (2008). Species richness and structure of three Neotropical bat assemblages. Biol. J. Linn. Society. 94 (3), 617–629. doi: 10.1111/j.1095-8312.2008.01014.x
Ribeiro, S., Sebaio, F., Branquinho, F., Marini, M., Vago, A., Braga, E. (2005). Avian malaria in Brazilian passerine birds: parasitism detected by nested PCR using DNA from stained blood smears. Parasitology 130 (03), 261–267. doi: 10.1017/S0031182004006596
Rulli, M. C., Santini, M., Hayman, D. T., D’Odorico, P. (2017). The nexus between forest fragmentation in Africa and Ebola virus disease outbreaks. Sci. Rep. 7, 41613. doi: 10.1038/srep41613
Salkeld, D. J., Padgett, K. A., Jones, J. H. (2013). A meta-analysis suggesting that the relationship between biodiversity and risk of zoonotic pathogen transmission is idiosyncratic. Ecol. Letters. 16 (5), 679–686. doi: 10.1111/ele.12101
Schmidt, K. A., Ostfeld, R. S. (2001). Biodiversity and the dilution effect in disease ecology. Ecology 82 (3), 609–619. doi: 10.1890/0012-9658(2001)082[0609:BATDEI]2.0.CO;2
Seltmann, A., Czirják, GÁ, Courtiol, A., Bernard, H., Struebig, M. J., Voigt, C. C. (2017). Habitat disturbance results in chronic stress and impaired health status in forest-dwelling paleotropical bats. Conserv. Physiol. 5 (1). doi: 10.1093/conphys/cox020
Shah, H. A., Huxley, P., Elmes, J., Murray, K. A. (2019). Agricultural land-uses consistently exacerbate infectious disease risks in southeast Asia. Nat. Commun. 10 (1), 1–13. doi: 10.1038/s41467-019-12333-z
Song, X.-P., Hansen, M. C., Stehman, S. V., Potapov, P. V., Tyukavina, A., Vermote, E. F., et al. (2018). Global land change from 1982 to 2016. Nature 560 (7720), 639. doi: 10.1038/s41586-018-0411-9
Suzán, G., Marcé, E., Giermakowski, J. T., Armién, B., Pascale, J., Mills, J., et al. (2008). The effect of habitat fragmentation and species diversity loss on hantavirus prevalence in Panama. Ann. New York Acad. Sci. 1149 (1), 80–83. doi: 10.1196/annals.1428.063
Tong, S., Chern, S.-W. W., Li, Y., Pallansch, M. A., Anderson, L. J. (2008). Sensitive and broadly reactive reverse transcription-PCR assays to detect novel paramyxoviruses. J. Clin. Microbiol. 46 (8), 2652–2658. doi: 10.1128/JCM.00192-08
Uezu, A., Metzger, J. P. (2011). Vanishing bird species in the Atlantic forest: relative importance of landscape configuration, forest structure and species characteristics. Biodiversity Conserv. 20 (14), 3627–3643. doi: 10.1007/s10531-011-0154-5
VanDevanter, D. R., Warrener, P., Bennett, L., Schultz, E. R., Coulter, S., Garber, R. L., et al. (1996). Detection and analysis of diverse herpesviral species by consensus primer PCR. J. Clin. Microbiol. 34 (7), 1666–1671. doi: 10.1128/jcm.34.7.1666-1671.1996
Van Vuren, P. J., Allam, M., Wiley, M. R., Ismail, A., Storm, N., Birkhead, M., et al. (2018). A novel adenovirus isolated from the Egyptian fruit bat in south Africa is closely related to recent isolates from China. Sci. Rep. 8 (1), 9584. doi: 10.1038/s41598-018-27836-w
Venables, W., Ripley, B. (2002). Modern applied statistics with s. 4th Edition ed (Berlin: Springer).
Vittor, A. Y., Pan, W., Gilman, R. H., Tielsch, J., Glass, G., Shields, T., et al. (2009). Linking deforestation to malaria in the Amazon: characterization of the breeding habitat of the principal malaria vector, anopheles darlingi. Am. J. Trop. Med. Hyg. 81 (1), 5.
Wijeyaratne, P., Arsenault, L. J., Murphy, C. (1994). Endemic disease and development: the leishmaniases. Acta Tropica. 56 (4), 349–364. doi: 10.1016/0001-706X(94)90106-6
Keywords: viral richness, diversity, bat host, deforestation, land-use change, viral prevalence
Citation: Loh EH, Nava A, Murray KA, Olival KJ, Guimarães M, Shimabukuro J, Zambrana-Torrelio C, Fonseca FR, de Oliveira DBL, Campos ACdA, Durigon EL, Ferreira F, Struebig MJ and Daszak P (2022) Prevalence of bat viruses associated with land-use change in the Atlantic Forest, Brazil. Front. Cell. Infect. Microbiol. 12:921950. doi: 10.3389/fcimb.2022.921950
Received: 16 April 2022; Accepted: 28 October 2022;
Published: 09 December 2022.
Edited by:
Timothy Yong James, University of Michigan, United StatesReviewed by:
Kelly A. Speer, Smithsonian Conservation Biology Institute (SI), United StatesJie Cui, Chinese Academy of Sciences (CAS), China
Copyright © 2022 Loh, Nava, Murray, Olival, Guimarães, Shimabukuro, Zambrana-Torrelio, Fonseca, de Oliveira, Campos, Durigon, Ferreira, Struebig and Daszak. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elizabeth H. Loh, eloh@transy.edu; Alessandra Nava, navaveterinaria@gmail.com