ORIGINAL RESEARCH article

Front. Cell. Infect. Microbiol., 09 December 2022

Sec. Clinical and Diagnostic Microbiology and Immunology

Volume 12 - 2022 | https://doi.org/10.3389/fcimb.2022.921950

Prevalence of bat viruses associated with land-use change in the Atlantic Forest, Brazil

  • 1. Division of Natural Sciences and Mathematics, Transylvania University, Lexington, KY, United States

  • 2. Durrell Institute of Conservation and Ecology, School of Anthropology and Conservation, University of Kent, Canterbury, United Kingdom

  • 3. Instituto Leônidas e Maria Deane – Fiocruz Amazônia, Manaus, Amazonas, Brazil

  • 4. MRC Unit The Gambia at London School of Hygiene and Tropical Medicine, Fajara, Gambia

  • 5. EcoHealth Alliance, New York, NY, United States

  • 6. Departamento de Recursos Naturais, Faculdade de Ciências Agronomicas, Universidade Estadual Paulista, Botucatu, Brazil

  • 7. Departamento 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

  • 8. Department of Environmental Science and Policy, George Mason University, Fairfax VA, United States

  • 9. Departamento de Microbiologia, Instituto de Ciências Biomédicas-II, Universidade de São Paulo, São Paulo, Brazil

Article metrics

View details

10

Citations

7k

Views

1,7k

Downloads

Abstract

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

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 15). Sequences were segregated into discrete viruses, defined as a viral species, based on distinct monophyletic clustering following Anthony et al. (2013).

Viral FamilyTargetAmplicon sizePrimer namePrimer sequence 5’- 3’Reference
Conventional PCR/Semi-Nested/Nested-PCR reactionAstroviridaeRNA-Dependent RNA Polymerase (RdRp)Round 1 431bpAstr4380FGAYTGGRCNCGNTWYGATGGNACIATAtkins et al., 2009
Astr4811RGGYTTNACCCACATNCCAAA
Round 2 342bpAstr4380F + Astr4722RARNCKRTCATCNCCATA
CoronaviridaeRNA-Dependent RNA Polymerase (RdRp)Round 1 520bpCoV-FWD1CGTTGGIACWAAYBTVCCWYTICARBTRGGQuan et al., 2010
CoV-RVS1GGTCATKATAGCRTCAVMASWWGCNACATG
Round 2 328pbCoV-FWD2GGCWCCWCCHGGNGARCAATT
CoV-RVS2GGWAWCCCCAYTGYTGWAYRTC
HerpesviridaePolymerase (Pol)Round 1 variableDFAgAY TTY gCN AgY YTN TAY CC
ILKTCC Tgg ACA AgC AgC ARN YSg CNM TNA AVan Devanter et al., 1996
KG1gTC TTg CTC ACC AgN TCN ACN CCY TT
Round 2215-315TGVTgT AAC TCg gTg TAY ggN TTY ACN ggN gT
IYGCAC AgA gTC CgT RTC NCC RTA DAT
ParamyxoviridaePolymerase (Pol)Round 1 639bpPAR-F1gAA ggI TAT TgT CAI AAR NTN Tgg ACTong et al., 2008
PAR-RgCT gAA gTT ACI ggI TCI CCD ATR TTN C
Round 2 561bpPAR-R +PAR-F2gTT gCT TCA ATg gTT CAR ggN gAY AA
SYBRGreenHantaviridaeS segment141bpJAN-FCCC TgT Tgg ATC AAC Tgg TTT TgAraujo et al., 2011
JAN-RTgT AAT gTg CTC TTg TTA ACg TCA TCT

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.

Table 1

Predictor variablesDefinition
TreatmentForested versus Deforested
SexMale versus Female
Pregnancy statusYes or No
AgeThree categories including: juveniles, subadults and adults
Genus12 unique genera
Species18 unique species
AbundanceTotal number of individuals captured per species

Description of predictor variables used in the generalized linear models.

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

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

FAMILY/SpeciesCapturesCorona-Herpes-Hanta-Astro-Paramyxo-
PHYLLOSTOMIDAE
 Artibeus lituratus6610000
 Artibeus fimbriatus3400000
 Artibeus planirostris13063141
 Carollia perspicillata5540000
 Desmodus rotundus100000
 Diaemus youngi200000
 Glossophaga soricina100000
 Phyllostomus hastatus300001
 Sturnira lilium1710000
 Vampyrodes caracciol900000
MOLOSSIDAE
 Molossus molossus300000
 Eumops glaucinus200000
VESPERTILIONIDAE
 Lasiurus blossevillii100000
 Myotis nigricans100000
 Myotis albescens300000
 Myotis riparius300000
 Myotis unidentified A100000
 Myotis unidentified B300000

Total captures of bat species and total viral detections in forested and deforested habitat in the Interior Atlantic Forest.

Figure 3

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.

Table 3

EstimateSEzP value
Intercept-3.2720.492-6.6463e-11
Treatment (forested vs deforested)-1.2740.536-2.3780.017
Sex (male vs female)1.5020.5342.8120.005
Age (juvenile vs adult)1.2470.5962.0940.036
Age (subadult vs adult)-0.5641.067-0.5280.597

Best model of viral detection with different categories of land-use change, sex and age.

Table 4

Model descriptionAIC ΔAICNull d.fResidual devianceResidual d.f.
~ treatment + sex + age6.63033319.24329
~ treatment + sex + age + abundance6.920.2933319.15328
~ treatment + sex + pregnancy + age + abundance8.291.6633319.11327
~ treatment + species + sex + pregnancy + age + abundance31.825.133318.52310

Logistic regression model selection table comparing top four models based on lowest AICc.

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.

Statements

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

  • 1

    AkaikeH. (1973). Information theory and an extension of the maximum likelihood principle (Budapest, Hungary: Akadémiai Kiadó).

  • 2

    AllenT.MurrayK. A.Zambrana-TorrelioC.MorseS. S.RondininiC.Di MarcoM.et al. (2017). Global hotspots and correlates of emerging zoonotic diseases. Nat. Commun.8 (1), 1124. doi: 10.1038/s41467-017-00923-8

  • 3

    AnthonyS. J.EpsteinJ. H.MurrayK. A.Navarrete-MaciasI.Zambrana-TorrelioC. M.SolovyovA.et al. (2013). A strategy to estimate unknown viral diversity in mammals. Mbio4 (5):10281038 doi: 10.1128/mBio.00598-13

  • 4

    AraujoJ.PereiraA.NardiM.HenriquesD.LautenschalagerD.DutraL.et al. (2011). Detection of hantaviruses in Brazilian rodents by SYBR-green-based real-time RT-PCR. Arch. virology.156 (7), 12691274. doi: 10.1007/s00705-011-0968-1

  • 5

    AtkinsA.WellehanJ. F.ChildressA. L.ArcherL. L.FraserW. A.CitinoS. B. (2009). Characterization of an outbreak of astroviral diarrhea in a group of cheetahs (Acinonyx jubatus). Veterinary Microbiol.136 (1), 160165. doi: 10.1016/j.vetmic.2008.10.035

  • 6

    BradleyC. A.GibbsS. E.AltizerS. (2008). Urban land use predicts West Nile virus exposure in songbirds. Ecol. Applications.18 (5), 10831092. doi: 10.1890/07-0822.1

  • 7

    CalisherC. H.ChildsJ. E.FieldH. E.HolmesK. V.SchountzTet al (2006). Bats: important reservoir hosts of emerging viruses. Clinical microbiology reviews19.3, 531545.

  • 8

    CarterG. G.ForssS.PageR. A.RatcliffeJ. M. (2018). Younger vampire bats (Desmodus rotundus) are more likely than adults to explore novel objects. PloS One13 (5), e0196889. doi: 10.1371/journal.pone.0196889

  • 9

    ChaoA.JostL. (2015). Estimating diversity and entropy profiles via discovery rates of new species. Methods Ecol. Evolution.6 (8), 873882. doi: 10.1111/2041-210X.12349

  • 10

    ChasarA.LoiseauC.ValkiūnasG.IezhovaT.SmithT. B.SehgalR. N. (2009). Prevalence and diversity patterns of avian blood parasites in degraded African rainforest habitats. Mol. Ecology.18 (19), 41214133. doi: 10.1111/j.1365-294X.2009.04346.x

  • 11

    CivitelloD. J.CohenJ.FatimaH.HalsteadN. T.LirianoJ.McMahonT. A.et al. (2015). Biodiversity inhibits parasites: broad evidence for the dilution effect. Proc. Natl. Acad. Sci.112 (28), 86678671. doi: 10.1073/pnas.1506279112

  • 12

    de FariaH. H.PiresA. (2006). Parque estadual do morro do diabo: plano de manejo (Santa Cruz do Rio Pardo: Editora Viena), 311.

  • 13

    DelavalM.Charles-DominiqueP. (2006). Edge effects on frugivorous and nectarivorous bat communities in a neotropical primary forest in French Guiana. Revue d'Ecologie, Terre et Vie61.4, 343352.

  • 14

    DuriganG.FrancoG. (2006). “Vegetação,” in Parque estadual do morro do diabo: plano de manejo (Santa Cruz do Rio Pardo: Editora Viena), 111118.

  • 15

    EpsteinJ. H.AnthonyS. J.IslamA.KilpatrickA. M.Ali KhanS.BalkeyM. D.et al. (2020). Nipah virus dynamics in bats and implications for spillover to humans. Proc. Natl. Acad. Sci.117 (46), 2919029201. doi: 10.1073/pnas.2000429117

  • 16

    FrickW. F.KingstonT.FlandersJ. (2020). A review of the major threats and challenges to global bat conservation. Ann. New York Acad. Sci.1469 (1), 525. doi: 10.1111/nyas.14045

  • 17

    GibbR.ReddingD. W.ChinK. Q.DonnellyC. A.BlackburnT. M.NewboldT.et al. (2020). Zoonotic host diversity increases in human-dominated ecosystems. Nature584 (7821), 398402. doi: 10.1038/s41586-020-2562-8

  • 18

    GorresenP. M.WilligM. R. (2004). Landscape responses of bats to habitat fragmentation in Atlantic forest of Paraguay. J. Mammalogy.85 (4), 688697. doi: 10.1644/BWG-125

  • 19

    GottdenkerN. L.StreickerD. G.FaustC. L.CarrollC. (2014). Anthropogenic land use change and infectious diseases: a review of the evidence. Ecohealth11 (4), 619632. doi: 10.1007/s10393-014-0941-z

  • 20

    HandleyC. O.Jr.WilsonD. E.GardnerA. L. (1991). Demography and natural history of the common fruit bat, artibeus jamaicensis, on barro Colorado island, Panama. doi: 10.5479/si.00810282.511

  • 21

    HuB.GeX.WangL.-F.ShiZ. (2015). Bat origin of human coronaviruses. Virol. J.12 (1), 110. doi: 10.1186/s12985-015-0422-1

  • 22

    JánoskaM.VidovszkyM.MolnárV.LiptovszkyM.HarrachB.BenkőM. (2011). Novel adenoviruses and herpesviruses detected in bats. Veterinary J.189 (1), 118121. doi: 10.1016/j.tvjl.2010.06.020

  • 23

    JonesK. E.PatelN.LevyM.StoreygardA.BalkD.GittlemanJ. L.et al. (2008). Global trends in emerging infectious diseases. Nature451, 990993. doi: 10.1038/nature06536

  • 24

    KilpatrickA. M. (2011). Globalization, land use, and the invasion of West Nile virus. Science334 (6054), 323327. doi: 10.1126/science.1201010

  • 25

    KilpatrickA. M.DaszakP.JonesM. J.MarraP. P.KramerL. D. (2006a). Host heterogeneity dominates West Nile virus transmission. Proc. R. Soc. London B: Biol. Sci.273 (1599), 23272333.

  • 26

    KilpatrickA. M.DaszakP.KramerL.MarraP. P.DobsonA. P.HudsonP. J.et al. (2006b). Predicting the transmission of West Nile virus. Am. J. Trop. Med. Hygiene.75 (5), 139.

  • 27

    KilpatrickA. M.JonesM.KramerL. D.MarraP. P.DaszakP. (2005). West Nile Virus vector ecology across an urbanization gradient. Am. J. Trop. Med. Hygiene.73 (6), 307308. doi: 10.3201/eid1103.040364

  • 28

    KingstonT. (2010). Research priorities for bat conservation in southeast Asia: a consensus approach. Biodiversity Conserv.19 (2), 471484. doi: 10.1007/s10531-008-9458-5

  • 29

    LauranceS. G.JonesD.WestcottD.MckeownA.HarringtonG.HilbertD. W. (2013). Habitat fragmentation and ecological traits influence the prevalence of avian blood parasites in a tropical rainforest landscape. PloS One8 (10), e76227. doi: 10.1371/journal.pone.0076227

  • 30

    LoGiudiceK.OstfeldR. S.SchmidtK. A.KeesingF. (2003). The ecology of infectious disease: effects of host diversity and community composition on Lyme disease risk. Proc. Natl. Acad. Sci.100 (2), 567571. doi: 10.1073/pnas.0233733100

  • 31

    LohE. H.Zambrana-TorrelioC.OlivalK. J.BogichT. L.JohnsonC. K.MazetJ. A.et al. (2015). Targeting transmission pathways for emerging zoonotic disease surveillance and control. Vector-Borne Zoonotic Diseases.15 (7), 432437. doi: 10.1089/vbz.2013.1563

  • 32

    MeyerC. F.KalkoE. K. (2008). Assemblage-level responses of phyllostomid bats to tropical forest fragmentation: land-bridge islands as a model system. J. Biogeography.35 (9), 17111726. doi: 10.1111/j.1365-2699.2008.01916.x

  • 33

    MurrayK. A.DaszakP. (2013). Human ecology in pathogenic landscapes: two hypotheses on how land use change drives viral emergence. Curr. Opin. virology.3 (1), 7983. doi: 10.1016/j.coviro.2013.01.006

  • 34

    OlivalK. J.HosseiniP. R.Zambrana-TorrelioC.RossN.BogichT. L.DaszakP. (2017). Host and viral traits predict zoonotic spillover from mammals. Nature546, 646650. doi: 10.1038/nature22975

  • 35

    OstfeldR. S.KeesingF. (2000). Biodiversity and disease risk: the case of Lyme disease. Conserv. Biol.14 (3), 722728. doi: 10.1046/j.1523-1739.2000.99014.x

  • 36

    QuanP.-L.FirthC.StreetC.HenriquezJ. A.PetrosovA.TashmukhamedovaA.et al. (2010). Identification of a severe acute respiratory syndrome coronavirus-like virus in a leaf-nosed bat in Nigeria. MBio1 (4), e00208e00210. doi: 10.1128/mBio.00208-10

  • 37

    RandolphS. E.DobsonA. (2012). Pangloss revisited: a critique of the dilution effect and the biodiversity-buffers-disease paradigm. Parasitology139 (7), 847863. doi: 10.1017/S0031182012000200

  • 38

    ReaserJ. K.HuntB. E.Ruiz-AravenaM.TaborG. M.PatzJ. A.BeckerD. 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

  • 39

    ReisN. R.FregoneziM. N.PeracchiA. L.ShibattaO. A. (2013). Morcegos do brasil: guia de campo (1st edn). (Rio de Janeiro: Technical Books Editora), 225 pp.

  • 40

    ReisN. R.PeracchiA. L.BatistaC. B.de LimaI. P.PereiraA. D. (2017). História natural dos morcegos brasileiros: chave de identificação de espécies (Rio de Janeiro: Technical Books Editora), 416.

  • 41

    RexK.KelmD. H.WiesnerK.KunzT. H.VoigtC. C. (2008). Species richness and structure of three Neotropical bat assemblages. Biol. J. Linn. Society.94 (3), 617629. doi: 10.1111/j.1095-8312.2008.01014.x

  • 42

    RibeiroS.SebaioF.BranquinhoF.MariniM.VagoA.BragaE. (2005). Avian malaria in Brazilian passerine birds: parasitism detected by nested PCR using DNA from stained blood smears. Parasitology130 (03), 261267. doi: 10.1017/S0031182004006596

  • 43

    RulliM. C.SantiniM.HaymanD. T.D’OdoricoP. (2017). The nexus between forest fragmentation in Africa and Ebola virus disease outbreaks. Sci. Rep.7, 41613. doi: 10.1038/srep41613

  • 44

    SalkeldD. J.PadgettK. A.JonesJ. H. (2013). A meta-analysis suggesting that the relationship between biodiversity and risk of zoonotic pathogen transmission is idiosyncratic. Ecol. Letters.16 (5), 679686. doi: 10.1111/ele.12101

  • 45

    SchmidtK. A.OstfeldR. S. (2001). Biodiversity and the dilution effect in disease ecology. Ecology82 (3), 609619. doi: 10.1890/0012-9658(2001)082[0609:BATDEI]2.0.CO;2

  • 46

    SeltmannA.CzirjákCourtiolA.BernardH.StruebigM. J.VoigtC. 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

  • 47

    ShahH. A.HuxleyP.ElmesJ.MurrayK. A. (2019). Agricultural land-uses consistently exacerbate infectious disease risks in southeast Asia. Nat. Commun.10 (1), 113. doi: 10.1038/s41467-019-12333-z

  • 48

    SongX.-P.HansenM. C.StehmanS. V.PotapovP. V.TyukavinaA.VermoteE. F.et al. (2018). Global land change from 1982 to 2016. Nature560 (7720), 639. doi: 10.1038/s41586-018-0411-9

  • 49

    SuzánG.MarcéE.GiermakowskiJ. T.ArmiénB.PascaleJ.MillsJ.et al. (2008). The effect of habitat fragmentation and species diversity loss on hantavirus prevalence in Panama. Ann. New York Acad. Sci.1149 (1), 8083. doi: 10.1196/annals.1428.063

  • 50

    TongS.ChernS.-W. W.LiY.PallanschM. A.AndersonL. J. (2008). Sensitive and broadly reactive reverse transcription-PCR assays to detect novel paramyxoviruses. J. Clin. Microbiol.46 (8), 26522658. doi: 10.1128/JCM.00192-08

  • 51

    UezuA.MetzgerJ. P. (2011). Vanishing bird species in the Atlantic forest: relative importance of landscape configuration, forest structure and species characteristics. Biodiversity Conserv.20 (14), 36273643. doi: 10.1007/s10531-011-0154-5

  • 52

    VanDevanterD. R.WarrenerP.BennettL.SchultzE. R.CoulterS.GarberR. L.et al. (1996). Detection and analysis of diverse herpesviral species by consensus primer PCR. J. Clin. Microbiol.34 (7), 16661671. doi: 10.1128/jcm.34.7.1666-1671.1996

  • 53

    Van VurenP. J.AllamM.WileyM. R.IsmailA.StormN.BirkheadM.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

  • 54

    VenablesW.RipleyB. (2002). Modern applied statistics with s. 4th Edition ed (Berlin: Springer).

  • 55

    VittorA. Y.PanW.GilmanR. H.TielschJ.GlassG.ShieldsT.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.

  • 56

    WijeyaratneP.ArsenaultL. J.MurphyC. (1994). Endemic disease and development: the leishmaniases. Acta Tropica.56 (4), 349364. doi: 10.1016/0001-706X(94)90106-6

  • 57

    WilkinsonG. S.Brunet-RossinniA. (2009). “Methods for age estimation and the study of senescence in bats,” in KunzT. H.ParsonsS. eds Ecological and behavioral methods for the study of bats (Baltimore, Maryland, John Hopkins University Press), 315325.

Summary

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

Volume

12 - 2022

Edited by

Timothy Yong James, University of Michigan, United States

Reviewed by

Kelly A. Speer, Smithsonian Conservation Biology Institute (SI), United States; Jie Cui, Chinese Academy of Sciences (CAS), China

Updates

Copyright

*Correspondence: Elizabeth H. Loh, ; Alessandra Nava,

This article was submitted to Clinical Microbiology, a section of the journal Frontiers in Cellular and Infection Microbiology

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics