Sociality, density-dependence and microclimates determine the persistence of populations suffering from a novel fungal disease, white-nose syndrome

White-nose syndrome, winter duration, and pre-hibernation climate impact abundance of reproductive female bats

White-nose syndrome (WNS) is an infectious disease that disrupts hibernation in bats, leading to premature exhaustion of fat stores. Though we know WNS does impact reproduction in hibernating female bats, we are unsure how these impacts are exacerbated by local climate factors. We compiled data from four southeastern U.S. states and used generalized linear mixed effects models to compare effects of WNS, pre-hibernation climate variables, and winter duration on the number of reproductive females in species across the range of WNS susceptibility. We predicted we would see a decline in the number of reproductive females in WNS-susceptible species, with the effect exaggerated by longer winter durations and pre-hibernation climate variables that lead to reductions in foraging. We found that the number of reproductive females in WNS-susceptible species was positively correlated with pre-hibernation local climate conditions conducive to foraging; however, WNS-susceptible species experienced an overall decline with the presence of WNS and as winter duration increased. Our long-term dataset provides evidence that pre-hibernation climate, specifically favorable summer weather conditions for foraging, greatly influences the reproduction, regardless of WNS status.

A Palearctic view of a bat fungal disease

The fungal infection causing white-nose disease in hibernating bats in North America has resulted in dramatic population declines of affected species, since the introduction of the causative agent Pseudogymnoascus destructans. The fungus is native to the Palearctic, where it also infects several bat species, yet rarely causes severe pathology or the death of the host. Pseudogymnoascus destructans infects bats during hibernation by invading and digesting the skin tissue, resulting in the disruption of torpor patterns and consequent emaciation. Relations among pathogen, host, and environment are complex, and individuals, populations, and species respond to the fungal pathogen in different ways. For example, the Nearctic Myotis lucifugus responds to infection by mounting a robust immune response, leading to immunopathology often contributing to mortality. In contrast, the Palearctic M. myotis shows no significant immunological response to infection. This lack of a strong response, resulting from the long coevolution between the hosts and the pathogen in the pathogen's native range, likely contributes to survival in tolerant species. After more than 15 years since the initial introduction of the fungus to North America, some of the affected populations are showing signs of recovery, suggesting that the fungus, hosts, or both are undergoing processes that may eventually lead to coexistence. The suggested or implemented management methods of the disease in North America have encompassed, for example, the use of probiotics and fungicides, vaccinations, and modifying the environmental conditions of the hibernation sites to limit the growth of the pathogen, intensity of infection, or the hosts' responses to it. Based on current knowledge from Eurasia, policy makers and conservation managers should refrain from disrupting the ongoing evolutionary processes and adopt a holistic approach to managing the epizootic.

Contribution of host species and pathogen clade to snake fungal disease hotspots in Europe

Infectious diseases are influenced by interactions between host and pathogen, and the number of infected hosts is rarely homogenous across the landscape. Areas with elevated pathogen prevalence can maintain a high force of infection and may indicate areas with disease impacts on host populations. However, isolating the ecological processes that result in increases in infection prevalence and intensity remains a challenge. Here we elucidate the contribution of pathogen clade and host species in disease hotspots caused by Ophidiomyces ophidiicola, the pathogen responsible for snake fungal disease, in 21 species of snakes infected with multiple pathogen strains across 10 countries in Europe. We found isolated areas of disease hotspots in a landscape where infections were otherwise low. O. ophidiicola clade had important effects on transmission, and areas with multiple pathogen clades had higher host infection prevalence. Snake species further influenced infection, with most positive detections coming from species within the Natrix genus. Our results suggest that both host and pathogen identity are essential components contributing to increased pathogen prevalence.

Factors Predicting Apparent Ophidiomycosis in Wild Brown Watersnakes (Nerodia taxispilota)

Ophidiomycosis, also known as snake fungal disease, is caused by Ophidiomyces ophidiicola and is a threat to snake conservation worldwide. Ophidiomycosis has been reported throughout much of the eastern US, and outbreaks have been associated with local population declines of already strained populations. Previous studies report significant variability in ophidiomycosis among species sampled, with higher prevalence typically observed in Nerodia spp. Although ophidiomycosis can lead to morbidity and mortality in affected individuals, little is known about disease dynamics in free-ranging populations. Herein, we examine how individual-specific factors (e.g., life stage [immature, mature], contaminant status, sex, hemograms) may be associated with ophidiomycosis status in the brown watersnake (Nerodia taxispilota). During 2018-19, we sampled 97 N. taxispilota from five locations along the Savannah River in South Carolina and Georgia, US. Ophidiomyces ophidiicola DNA was detected in 66 snakes for a prevalence of 68% (95% confidence interval, 59-77). Mature snakes had a significantly higher risk of apparent ophidiomycosis (skin lesions present and quantitative PCR [qPCR], positive) relative to immature snakes. Snakes classified as having possible (skin lesions present, but qPCR negative) or apparent ophidiomycosis exhibited a relative azurophilia and heterophilia compared with individuals classified as negative (P≤0.037). Nerodia taxispilota in this region appear to have a high prevalence of apparent ophidiomycosis (22%; 95% CI, 14-31), similar to previous reports from the southeastern US. Additional epidemiologic investigations are warranted to further elucidate other individual-specific and environmental factors that may dictate disease risk and outcomes in affected populations.

White-nose syndrome restructures bat skin microbiomes

IMPORTANCE

Inherent complexities in the composition of microbiomes can often preclude investigations of microbe-associated diseases. Instead of single organisms being associated with disease, community characteristics may be more relevant. Longitudinal microbiome studies of the same individual bats as pathogens arrive and infect a population are the ideal experiment but remain logistically challenging; therefore, investigations like our approach that are able to correlate invasive pathogens to alterations within a microbiome may be the next best alternative. The results of this study potentially suggest that microbiome-host interactions may determine the likelihood of infection. However, the contrasting relationship between Pd and the bacterial microbiomes of Myotis lucifugus and Perimyotis subflavus indicate that we are just beginning to understand how the bat microbiome interacts with a fungal invader such as Pd.

Host-pathogen interactions under pressure: A review and meta-analysis of stress-mediated effects on disease dynamics

Human activities have increased the intensity and frequency of natural stressors and created novel stressors, altering host-pathogen interactions and changing the risk of emerging infectious diseases. Despite the ubiquity of such anthropogenic impacts, predicting the directionality of outcomes has proven challenging. Here, we conduct a review and meta-analysis to determine the primary mechanisms through which stressors affect host-pathogen interactions and to evaluate the impacts stress has on host fitness (survival and fecundity) and pathogen infectivity (prevalence and intensity). We assessed 891 effect sizes from 71 host species (representing seven taxonomic groups) and 78 parasite taxa from 98 studies. We found that infected and uninfected hosts had similar sensitivity to stressors and that responses varied according to stressor type. Specifically, limited resources compromised host fecundity and decreased pathogen intensity, while abiotic environmental stressors (e.g., temperature and salinity) decreased host survivorship and increased pathogen intensity, and pollution increased mortality but decreased pathogen prevalence. We then used our meta-analysis results to develop susceptible-infected theoretical models to illustrate scenarios where infection rates are expected to increase or decrease in response to resource limitations or environmental stress gradients. Our results carry implications for conservation and disease emergence and reveal areas for future work.

Contact-independent exposure to Rhodococcus rhodochrous DAP96253 volatiles does not improve the survival rate of Myotis lucifugus (little brown bats) affected by White-nose Syndrome

Since the emergence of White-nose Syndrome, a fungal disease in bats, caused by Pseudogymnoascus destructans, hibernating populations of little brown bats (Myotis lucifugus) have declined by 70-90% within P. destructans positive hibernacula. To reduce the impact of White-nose Syndrome to North American little brown bat populations we evaluated if exposure to volatile organic compounds produced by induced cells from Rhodococcus rhodochrous strain DAP96253 could improve the overwinter survival of bats infected by P. destructans. Two simultaneous field treatment trials were conducted at natural hibernacula located in Rockcastle and Breckinridge counties, Kentucky, USA. A combined total of 120 little brown bats were randomly divided into control groups (n = 60) which were not exposed to volatile organic compounds and treatment groups (n = 60) which were exposed to volatile organic compounds produced by non-growth, fermented cell paste composed of R. rhodochrous strain DAP96253 cells. Cox proportional hazard models revealed a significant decreased survival at the Rockcastle field trial site but not the Breckinridge field site. At the Breckinridge hibernacula, overwinter survival for both treatment and control groups were 60%. At the Rockcastle hibernacula, Kaplan-Meier survival curves indicated significantly increased overwinter survival of bats in the control group (43% survived) compared to the treatment group (20% survived). Although complete inhibition of P. destructans by volatile organic compounds produced by induced R. rhodochrous strain DAP96253 cells was observed in vitro studies, our results suggest that these volatile organic compounds do not inhibit P. destructans in situ and may promote P. destructans growth.

BatTool: projecting bat populations facing multiple stressors using a demographic model

Bats provide ecologically and agriculturally important ecosystem services but are currently experiencing population declines caused by multiple environmental stressors, including mortality from white-nose syndrome and wind energy development. Analyses of the current and future health and viability of these species may support conservation management decision making. Demographic modeling provides a quantitative tool for decision makers and conservation managers to make more informed decisions, but widespread adoption of these tools can be limited because of the complexity of the mathematical, statistical, and computational components involved in implementing these models. In this work, we provide an exposition of the BatTool R package, detailing the primary components of the matrix projection model, a publicly accessible graphical user interface ( https://rconnect.usgs.gov/battool ) facilitating user-defined scenario analyses, and its intended uses and limitations (Wiens et al., US Geol Surv Data Release 2022; Wiens et al., US Geol Surv Softw Release 2022). We present a case study involving wind energy permitting, weighing the effects of potential mortality caused by a hypothetical wind energy facility on the projected abundance of four imperiled bat species in the Midwestern United States.

Habitat suitability mapping using logistic regression analysis of long-term bioacoustic bat survey dataset in the Cassadaga Creek watershed (USA)

Bat species show global ecological importance, yet their numbers are declining worldwide. Understanding bat-habitat interactions is crucial in terms of developing effective conservation plans. In an effort to model bat habitat suitability in the Cassadaga Creek watershed, long-term bioacoustic bat data (spanning 2009-2020) was compiled, georeferenced and statistically analyzed using logistic regression techniques. In total, 1600 bat occurrence records from five species of bat (559 Eptesicus fuscus, 560 Lasionycteris noctivagans, 143 Lasiurus borealis, 260 Lasiurus cinereus, and 78 Myotis lucifugus) were paired with pseudo-absence points to study the relationship between bat calling behavior and land cover. All bats but Myotis lucifugus had a statistically significant relationship with forested land cover, and all bats had negative interactions with agricultural habitats. Geospatial data was coupled with the statistical output to create maps of habitat suitability and echolocation calling density. This work provides a model that can be employed worldwide to evaluate bat habitat needs or patterns in echolocation behavior. Future research will incorporate a more recently collected dataset that is of greater geographic diversity with a larger number of environmental variables in the species distribution model.

Host infection and disease-induced mortality modify species contributions to the environmental reservoir

Environmental pathogen reservoirs exist for many globally important diseases and can fuel epidemics, influence pathogen evolution, and increase the threat of host extinction. Species composition can be an important factor that shapes reservoir dynamics and ultimately determines the outcome of a disease outbreak. However, disease-induced mortality can change species communities, indicating that species responsible for environmental reservoir maintenance may change over time. Here we examine the reservoir dynamics of Pseudogymnoascus destructans, the fungal pathogen that causes white-nose syndrome in bats. We quantified changes in pathogen shedding, infection prevalence and intensity, host abundance, and the subsequent propagule pressure imposed by each species over time. We find that highly shedding species are important during pathogen invasion, but contribute less over time to environmental contamination as they also suffer the greatest declines. Less infected species remain more abundant, resulting in equivalent or higher propagule pressure. More broadly, we demonstrate that high infection intensity and subsequent mortality during disease progression can reduce the contributions of high-shedding species to long-term pathogen maintenance.

Divergent population responses following salamander mass mortalities and declines driven by the emerging pathogen Batrachochytrium salamandrivorans

Understanding wildlife responses to novel threats is vital in counteracting biodiversity loss. The emerging pathogen Batrachochytrium salamandrivorans (Bsal) causes dramatic declines in European salamander populations, and is considered an imminent threat to global amphibian biodiversity. However, real-life disease outcomes remain largely uncharacterized. We performed a multidisciplinary assessment of the longer-term impacts of Bsal on highly susceptible fire salamander (Salamandra salamandra) populations, by comparing four of the earliest known outbreak sites to uninfected sites. Based on large-scale monitoring efforts, we found population persistence in strongly reduced abundances to over a decade after Bsal invasion, but also the extinction of an initially small-sized population. In turn, we found that host responses varied, and Bsal detection remained low, within surviving populations. Demographic analyses indicated an ongoing scarcity of large reproductive adults with potential for recruitment failure, while spatial comparisons indicated a population remnant persisting within aberrant habitat. Additionally, we detected no early signs of severe genetic deterioration, yet nor of increased host resistance. Beyond offering additional context to Bsal-driven salamander declines, results highlight how the impacts of emerging hypervirulent pathogens can be unpredictable and vary across different levels of biological complexity, and how limited pathogen detectability after population declines may complicate surveillance efforts.

Mumps epidemic dynamics in the United States before vaccination (1923-1932)

Mumps is a vaccine-preventable, reemerging, and highly transmissible infectious disease. Widespread vaccination dramatically reduced cases; however, case counts have been increasing over the past 20 years. To provide a quantitative overview of historical mumps dynamics that can act as baseline information to help identify causes of mumps reemergence, we analyzed timeseries of cases reported from 1923 to 1932 in the United States. During that time, 239,230 mumps cases were reported in 70 cities. Larger cities reported annual epidemics and smaller cities reported intermittent, sporadic outbreaks. The critical community size above which transmission continuously occurred was likely between 365,583 and 781,188 individuals but could range as high as 3,376,438 individuals. Mumps cases increased as city size increased, suggesting density-dependent transmission. Using a density-dependent SEIR model, we calculated a mean effective reproductive number (Re) of 1.2. Re varied by city and over time, with periodic high values that could characterize short periods of very high transmission known as superspreading events. Case counts most often peaked in March, with higher-than-average transmission from December through April and showed a correlation with weekly births. While certain city pairs in Midwestern states had synchronous outbreaks, most outbreaks were less synchronous and not driven by distance between cities. This work demonstrates the importance of long-term infectious disease surveillance data and will inform future studies on mumps reemergence and control.

COULD WHITE-NOSE SYNDROME MANIFEST DIFFERENTLY IN MYOTIS LUCIFUGUS IN WESTERN VERSUS EASTERN REGIONS OF NORTH AMERICA? A REVIEW OF FACTORS

White-nose syndrome (WNS) has notably affected the abundance of Myotis lucifugus (little brown myotis) in North America. Thus far, substantial mortality has been restricted to the eastern part of the continent where the cause of WNS, the invasive fungus Pseudogymnoascus destructans, has infected bats since 2006. To date, the state of Washington is the only area in the Western US or Canada (the Rocky Mountains and further west in North America) with confirmed cases of WNS in bats, and there the disease has spread more slowly than it did in Eastern North America. Here, we review differences between M. lucifugus in western and eastern parts of the continent that may affect transmission, spread, and severity of WNS in the West and highlight important gaps in knowledge. We explore the hypothesis that western M. lucifugus may respond differently to WNS on the basis of different hibernation strategies, habitat use, and greater genetic structure. To document the effect of WNS on M. lucifugus in the West most effectively, we recommend focusing on maternity roosts for strategic disease surveillance and monitoring abundance. We further recommend continuing the challenging work of identifying hibernation and swarming sites to better understand the microclimates, microbial communities, and role in disease transmission of these sites, as well as the ecology and hibernation physiology of bats in noncavernous hibernacula.

Sex-biased infections scale to population impacts for an emerging wildlife disease

Demographic factors are fundamental in shaping infectious disease dynamics. Aspects of populations that create structure, like age and sex, can affect patterns of transmission, infection intensity and population outcomes. However, studies rarely link these processes from individual to population-scale effects. Moreover, the mechanisms underlying demographic differences in disease are frequently unclear. Here, we explore sex-biased infections for a multi-host fungal disease of bats, white-nose syndrome, and link disease-associated mortality between sexes, the distortion of sex ratios and the potential mechanisms underlying sex differences in infection. We collected data on host traits, infection intensity and survival of five bat species at 42 sites across seven years. We found females were more infected than males for all five species. Females also had lower apparent survival over winter and accounted for a smaller proportion of populations over time. Notably, female-biased infections were evident by early hibernation and likely driven by sex-based differences in autumn mating behaviour. Male bats were more active during autumn which likely reduced replication of the cool-growing fungus. Higher disease impacts in female bats may have cascading effects on bat populations beyond the hibernation season by limiting recruitment and increasing the risk of Allee effects.

Environmental transmission of Pseudogymnoascus destructans to hibernating little brown bats

Pathogens with persistent environmental stages can have devastating effects on wildlife communities. White-nose syndrome (WNS), caused by the fungus Pseudogymnoascus destructans, has caused widespread declines in bat populations of North America. In 2009, during the early stages of the WNS investigation and before molecular techniques had been developed to readily detect P. destructans in environmental samples, we initiated this study to assess whether P. destructans can persist in the hibernaculum environment in the absence of its conclusive bat host and cause infections in naive bats. We transferred little brown bats (Myotis lucifugus) from an unaffected winter colony in northwest Wisconsin to two P. destructans contaminated hibernacula in Vermont where native bats had been excluded. Infection with P. destructans was apparent on some bats within 8 weeks following the introduction of unexposed bats to these environments, and mortality from WNS was confirmed by histopathology at both sites 14 weeks following introduction. These results indicate that environmental exposure to P. destructans is sufficient to cause the infection and mortality associated with WNS in naive bats, which increases the probability of winter colony extirpation and complicates conservation efforts.

Shifting effects of host physiological condition following pathogen establishment

Understanding host persistence with emerging pathogens is essential for conserving populations. Hosts may initially survive pathogen invasions through pre-adaptive mechanisms. However, whether pre-adaptive traits are directionally selected to increase in frequency depends on the heritability and environmental dependence of the trait and the costs of trait maintenance. Body condition is likely an important pre-adaptive mechanism aiding in host survival, although can be seasonally variable in wildlife hosts. We used data collected over 7 years on bat body mass, infection and survival to determine the role of host body condition during the invasion and establishment of the emerging disease, white-nose syndrome. We found that when the pathogen first invaded, bats with higher body mass were more likely to survive, but this effect dissipated following the initial epizootic. We also found that heavier bats lost more weight overwinter, but fat loss depended on infection severity. Lastly, we found mixed support that bat mass increased in the population after pathogen arrival; high annual plasticity in individual bat masses may have reduced the potential for directional selection. Overall, our results suggest that some factors that contribute to host survival during pathogen invasion may diminish over time and are potentially replaced by other host adaptations.

Higher white-nose syndrome fungal isolate yields from UV-guided wing biopsies compared with skin swabs and optimal culture media

BACKGROUND

North American bat populations have suffered severe declines over the last decade due to the Pseudogymnoascus destructans fungus infection. The skin disease associated with this causative agent, known as white-nose syndrome (WNS), is specific to bats hibernating in temperate regions. As cultured fungal isolates are required for epidemiological and phylogeographical studies, the purpose of the present work was to compare the efficacy and reliability of different culture approaches based on either skin swabs or wing membrane tissue biopsies for obtaining viable fungal isolates of P. destructans.

RESULTS

In total, we collected and analysed 69 fungal and 65 bacterial skin swabs and 51 wing membrane tissue biopsies from three bat species in the Czech Republic, Poland and the Republic of Armenia. From these, we obtained 12 viable P. destructans culture isolates.

CONCLUSIONS

Our results indicated that the efficacy of cultures based on wing membrane biopsies were significantly higher. Cultivable samples tended to be based on collections from bats with lower body surface temperature and higher counts of UV-visualised lesions. While cultures based on both skin swabs and wing membrane tissue biopsies can be utilised for monitoring and surveillance of P. destructans in bat populations, wing membrane biopsies guided by UV light for skin lesions proved higher efficacy. Interactions between bacteria on the host's skin also appear to play an important role.

Synchronicity of viral shedding in molossid bat maternity colonies

Infection dynamics in vertebrates are driven by biological and ecological processes. For bats, population structure and reproductive cycles have major effects on RNA virus transmission. On Reunion Island, previous studies have shown that parturition of pregnant females and aggregation of juvenile Reunion free-tailed bats (Mormopterus francoismoutoui) are associated with major increase in the prevalence of bats shedding RNA viruses. The synchronicity of such shedding pulses, however, is yet to be assessed between viruses but also maternity colonies. Based on 3422 fresh faeces collected every 2-5 weeks during four consecutive birthing seasons, we report the prevalence of bats shedding astroviruses (AstVs), coronaviruses (CoVs) and paramyxoviruses (PMVs) in two maternity colonies on Reunion Island. We found that the proportion of bats shedding viruses is highly influenced by sampling collection periods, and therefore by the evolution of the population age structure. We highlight that virus shedding patterns are consistent among years and colonies for CoVs and to a lesser extent for PMVs, but not for AstVs. We also report that 1% of bats harbour co-infections, with two but not three of the viruses, and most co-infections were due to CoVs and PMVs.

U.S. National Park visitor perceptions and behavioral intentions towards actions to prevent white-nose syndrome

In the United States, the discovery and spread of white-nose syndrome (WNS) has drastically changed how bats and caves are managed. The U.S. National Park Service has been instrumental in the national response to WNS, as it manages extensive cave resources and has a close relationship with the public. However, managers lack information on visitor support for disease prevention measures designed to slow the spread of WNS and minimize human disturbance of vulnerable bat populations. This study utilized the Theory of Planned Behavior to determine how visitor attitudes, subjective norms, and perceived behavioral controls influenced their behavior regarding WNS preventive actions, including participation in educational programming on bats, wearing clothes or shoes in caves that have not been contaminated with the fungus that causes WNS, walking over decontamination mats, and complying with cave closures. During summer of 2019, data were collected using an on-site survey of 1365 visitors to eight U.S. national park units: Oregon Caves, Lava Beds, Carlsbad Caverns, El Malpais, Wind Cave, Jewel Cave, Mammoth Cave, and Cumberland Gap. Visitors were willing to participate in all preventative actions addressed in the survey (77.7%-96.7%). Visitors expressed that engaging in these actions was very desirable (36.0%-65.6%), and their decision to engage in these actions was most strongly influenced by park staff (39.2%-68.8%) or signage (35.5%-61.9%). Attitudes and subjective norms were positive predictors of behavioral intentions for all measures. Perceived behavioral control was not a direct predictor for behavioral intent, but its interaction with attitudes and subjective norms had a moderating influence on intention to comply with multiple WNS preventive actions. With the continued spread of WNS and emergence of other threats to bats, understanding visitor behavioral intent and underlying factors will facilitate successful implementation of preventive actions that are publicly supported and promote conservation of bat populations in U.S. national parks.

Re-examining extreme sleep duration in bats: implications for sleep phylogeny, ecology, and function

Bats, quoted as sleeping for up to 20 h a day, are an often used example of extreme sleep duration amongst mammals. Given that duration has historically been one of the primary metrics featured in comparative studies of sleep, it is important that species specific sleep durations are well founded. Here, we re-examined the evidence for the characterization of bats as extreme sleepers and discuss whether it provides a useful representation of the sleep behavior of Chiroptera. Although there are a wealth of activity data to suggest that the diurnal cycle of bats is dominated by rest, estimates of sleep time generated from electrophysiological analyses suggest considerable interspecific variation, ranging from 83% to a more moderate 61% of the 24 h day spent asleep. Temperature-dependent changes in the duration and electroencephalographic profile of sleep suggest that bats represent a unique model for investigating the relationship between sleep and torpor. Further sources of intra-specific variation in sleep duration, including the impact of artificial laboratory environments and sleep intensity, remain unexplored. Future studies conducted in naturalistic environments, using larger sample sizes and relying on a pre-determined set of defining criteria will undoubtedly provide novel insights into sleep in bats and other species.

Functional Redundancy in Bat Microbial Assemblage in the Presence of the White Nose Pathogen

Understanding how host-associated microbial assemblages respond to pathogen invasion has implications for host health. Until recently, most investigations have focused on understanding the taxonomic composition of these assemblages. However, recent studies have suggested that microbial assemblage taxonomic composition is decoupled from its function, with assemblages being taxonomically varied but functionally constrained. The objective of this investigation was to understand how the Tri-colored bat, Perimyotis subflavus cutaneous microbial assemblage responds to fungal pathogen invasion within a functional context. We hypothesized that at a broad scale (e.g., KEGG pathways), there will be no difference in the functional assemblages between the white nose pathogen, Pseudogymnoascus destructans, positive and negative bats; and this pattern will be driven by the functional redundancy of bacterial taxa. At finer scales (e.g., gene models), we postulate differences in function attributed to interactions between bacteria and P. destructans, resulting in the production of antifungal metabolites. To test this, we used a combination of shotgun metagenomic and amplicon sequencing to characterize the bat cutaneous microbial assemblage in the presence/absence of P. destructans. Results showed that while there was a shift in taxonomic assemblage composition between P. destructans positive and negative bats, there was little overall difference in microbial function. Functional redundancy across bacterial taxa was clear at a broad-scale; however, both redundancy and variation in bacterial capability related to defense against pathogens was evident at finer scales. While functionality of the microbial assemblage was largely conserved in relation to P. destructans, the roles of particular functional pathways in resistance to fungal pathogens require further attention.

Emergence activity at hibernacula differs among four bat species affected by white-nose syndrome

Prior to the introduction of white‐nose syndrome (WNS) to North America, temperate bats were thought to remain within hibernacula throughout most of the winter. However, recent research has shown that bats in the southeastern United States emerge regularly from hibernation and are active on the landscape, regardless of their WNS status. The relationship between winter activity and susceptibility to WNS has yet to be explored but warrants attention, as it may enable managers to implement targeted management for WNS‐affected species. We investigated this relationship by implanting 1346 passive integrated transponder (PIT) tags in four species that vary in their susceptibility to WNS. Based on PIT‐tag detections, three species entered hibernation from late October to early November. Bats were active at hibernacula entrances on days when midpoint temperatures ranged from −1.94 to 22.78°C (mean midpoint temperature = 8.70 ± 0.33°C). Eastern small‐footed bats (Myotis leibii), a species with low susceptibility to WNS, were active throughout winter, with a significant decrease in activity in mid‐hibernation (December 16 to February 15). Tricolored bats (Perimyotis subflavus), a species that is highly susceptible to WNS, exhibited an increase in activity beginning in mid‐hibernation and extending through late hibernation (February 16 to March 31). Indiana bats (M. sodalis), a species determined to have a medium–high susceptibility to WNS, remained on the landscape into early hibernation (November 1 to December 15), after which we did not record any again until the latter portion of mid‐hibernation. Finally, gray bats (M. grisescens), another species with low susceptibility to WNS, maintained low but regular levels of activity throughout winter. Given these results, we determined that emergence activity from hibernacula during winter is highly variable among bat species and our data will assist wildlife managers to make informed decisions regarding the timing of implementation of species‐specific conservation actions.

Changes in hibernating tricolored bat (Perimyotis subflavus) roosting behavior in response to white-nose syndrome

Understanding animals' behavioral and physiological responses to pathogenic diseases is critical for management and conservation. One such disease, white‐nose syndrome (WNS), has greatly affected bat populations throughout eastern North America leading to significant population declines in several species. Although tricolored bat (Perimyotis subflavus) populations have experienced significant declines, little research has been conducted on their responses to the disease, particularly in the southeastern United States. Our objective was to document changes in tricolored bat roost site use after the appearance of WNS in a hibernaculum in the southeastern U.S. and relate these to microsite temperatures, ambient conditions, and population trends. We censused a tricolored bat hibernaculum in northwestern South Carolina, USA, once each year between February 26 and March 2, 2014–2021, and recorded species, section of the tunnel, distance from the entrance, and wall temperature next to each bat. The number of tricolored bats in the hibernaculum dropped by 90.3% during the first 3 years after the arrival of WNS. However, numbers stabilized and slightly increased from 2018 to 2021. Prior to the arrival of WNS, 95.6% of tricolored bats roosted in the back portion of the tunnel that was the warmest. After the arrival of WNS, we observed a significant increase in the proportion of bats using the front, colder portions of the tunnel, particularly during the period of population stabilization and increase. Roost temperatures of bats were also positively associated with February external temperatures. Our results suggest that greater use of the colder sections of the tunnel by tricolored bats could have led to increased survival due to slower growth rates of the fungus that causes WNS in colder temperatures or decreased energetic costs associated with colder hibernation temperatures. Thus, management actions that provide cold hibernacula may be an option for long‐term management of hibernacula, particularly in southern regions.

Great diversity of KSα sequences from bat-associated microbiota suggests novel sources of uncharacterized natural products

Polyketide synthases (PKSs) are multidomain enzymes in microorganisms that synthesize complex, bioactive molecules. PKS II systems are iterative, containing only a single representative of each domain: ketosynthase alpha (KS[Formula: see text]), ketosynthase beta and the acyl carrier protein. Any gene encoding for one of these domains is representative of an entire PKS II biosynthetic gene cluster (BGC). Bat skin surfaces represent an extreme environment prolific in Actinobacteria that may constitute a source for bioactive molecule discovery. KS[Formula: see text] sequences were obtained from culturable bacteria from bats in the southwestern United States. From 467 bat bacterial isolates, we detected 215 (46%) had KS[Formula: see text] sequences. Sequencing yielded 210 operational taxonomic units, and phylogenetic placement found 45 (21%) shared <85% homology to characterized metabolites. Additionally, 16 Actinobacteria genomes from the bat microbiome were analyzed for biosynthetic capacity. A range of 69-93% of the BGCs were novel suggesting the bat microbiome may contain valuable uncharacterized natural products. Documenting and characterizing these are important in understanding the susceptibility of bats to emerging infectious diseases, such as white-nose syndrome. Also noteworthy was the relationship between KS [Formula: see text] homology and total BGC novelty within each fully sequenced strain. We propose amplification and detection of KS[Formula: see text] could predict a strain's global biosynthetic capacity.

Winter torpor expression varies in four bat species with differential susceptibility to white-nose syndrome

Studies examining the overwintering behaviors of North American hibernating bats are limited to a handful of species. We deployed temperature-sensitive transmitters on four species of bat that exhibit differences in their susceptibility to white nose syndrome (WNS; Myotis grisescens, M. leibii, M. sodalis, and Perimyotis subflavus) to determine if these differences are correlated with behavior exhibited during hibernation (i.e., torpor expression and arousal frequency). Mean torpor skin temperature (T_sk) and torpor bout duration varied significantly among species (P ≤ 0.024), but arousal T_sk and duration did not (P ≥ 0.057). One of the species with low susceptibility to WNS, M. leibii, had significantly shorter torpor bout durations (37.67 ± 26.89 h) than M. sodalis (260.67 ± 41.33 h), the species with medium susceptibility to WNS. Myotis leibii also had significantly higher torpor T_sk (18.57 °C ± 0.20) than M. grisescens (13.33 °C ± 0.60), a second species with low WNS susceptibility. The high susceptibility species, Perimyotis subflavus, exhibited low torpor T_sk (14.42 °C ± 0.36) but short torpor bouts (72.36 ± 32.16 h). We demonstrate that the four cavernicolous species examined exhibit a wide range in torpid skin temperature and torpor bout duration. Information from this study may improve WNS management in multispecies hibernacula or individual species management by providing insight into how some species may differ in their techniques for overwinter survival.

Cooling of bat hibernacula to mitigate white-nose syndrome

White-nose syndrome (WNS) is a fungal disease that has caused precipitous declines in several North American bat species, creating an urgent need for conservation. We examined how microclimates and other characteristics of hibernacula have affected bat populations following WNS-associated declines and evaluated whether cooling of warm, little-used hibernacula could benefit bats. During the period following mass mortality (2013-2020), we conducted 191 winter surveys of 25 unmanipulated hibernacula and 6 manipulated hibernacula across Pennsylvania (USA). We joined these data with additional datasets on historical (pre-WNS) bat counts and on the spatial distribution of underground sites. We used generalized linear mixed models and model selection to identify factors affecting bat populations. Winter counts of Myotis lucifugus were higher and increased over time in colder hibernacula (those with midwinter temperatures of 3-6 °C) compared with warmer (7-11 °C) hibernacula. Counts of Eptesicus fuscus, Myotis leibii, and Myotis septentrionalis were likewise higher in colder hibernacula (temperature effects = -0.73 [SE 0.15], -0.51 [0.18], and -0.97 [0.28], respectively). Populations of M. lucifugus and M. septentrionalis increased most over time in hibernacula surrounded by more nearby sites, whereas Eptesicus fuscus counts remained high where they had been high before WNS onset (pre-WNS high count effect = 0.59 [0.22]). Winter counts of M. leibii were higher in hibernacula with high vapor pressure deficits (VPDs) (particularly over 0.1 kPa) compared with sites with lower VPDs (VPD effect = 15.3 [4.6]). Counts of M. lucifugus and E. fuscus also appeared higher where VPD was higher. In contrast, Perimyotis subflavus counts increased over time in relatively warm hibernacula and were unaffected by VPD. Where we manipulated hibernacula, we achieved cooling of on average 2.1 °C. At manipulated hibernacula, counts of M. lucifugus and P. subflavus increased over time (years since manipulation effect = 0.70 [0.28] and 0.51 [0.15], respectively). Further, there were more E. fuscus where cooling was greatest (temperature difference effect = -0.46 [SE 0.11]), and there was some evidence there were more P. subflavus in hibernacula sections that remained warm after manipulation. These data show bats are responding effectively to WNS through habitat selection. In M. lucifugus, M. septentrionalis, and possibly P. subflavus, this response is ongoing, with bats increasingly aggregating at suitable hibernacula, whereas E. fuscus remain in previously favored sites. Our results suggest that cooling warm sites receiving little use by bats is a viable strategy for combating WNS.

Host traits and environment interact to determine persistence of bat populations impacted by white-nose syndrome

Emerging infectious diseases have resulted in severe population declines across diverse taxa. In some instances, despite attributes associated with high extinction risk, disease emergence and host declines are followed by host stabilisation for unknown reasons. While host, pathogen, and the environment are recognised as important factors that interact to determine host-pathogen coexistence, they are often considered independently. Here, we use a translocation experiment to disentangle the role of host traits and environmental conditions in driving the persistence of remnant bat populations a decade after they declined 70-99% due to white-nose syndrome and subsequently stabilised. While survival was significantly higher than during the initial epidemic within all sites, protection from severe disease only existed within a narrow environmental space, suggesting host traits conducive to surviving disease are highly environmentally dependent. Ultimately, population persistence following pathogen invasion is the product of host-pathogen interactions that vary across a patchwork of environments.

Experimental inoculation trial to determine the effects of temperature and humidity on White-nose Syndrome in hibernating bats

Disease results from interactions among the host, pathogen, and environment. Inoculation trials can quantify interactions among these players and explain aspects of disease ecology to inform management in variable and dynamic natural environments. White-nose Syndrome, a disease caused by the fungal pathogen, Pseudogymnoascus destructans (Pd), has caused severe population declines of several bat species in North America. We conducted the first experimental infection trial on the tri-colored bat, Perimyotis subflavus, to test the effect of temperature and humidity on disease severity. We also tested the effects of temperature and humidity on fungal growth and persistence on substrates. Unexpectedly, only 37% (35/95) of bats experimentally inoculated with Pd at the start of the experiment showed any infection response or disease symptoms after 83 days of captive hibernation. There was no evidence that temperature or humidity influenced infection response. Temperature had a strong effect on fungal growth on media plates, but the influence of humidity was more variable and uncertain. Designing laboratory studies to maximize research outcomes would be beneficial given the high costs of such efforts and potential for unexpected outcomes. Understanding the influence of microclimates on host-pathogen interactions remains an important consideration for managing wildlife diseases, particularly in variable environments.

Contaminant exposure as an additional stressor to bats affected by white-nose syndrome: current evidence and knowledge gaps

Bats are exposed to numerous threats including pollution and emerging diseases. In North America, the fungal disease white-nose syndrome (WNS) has caused declines in many bat species. While the mechanisms of WNS have received considerable research attention, possible influences of contaminants have not. Herein, we review what is known about contaminant exposure and toxicity for four species whose populations have been severely affected by WNS (Myotis sodalis, M. septentrionalis, M. lucifugus, and Perimyotis subflavus) and identify temporal and spatial data gaps. We determine that there is limited information about the effects of contaminants on bats, and many compounds that have been detected in these bat species have yet to be evaluated for toxicity. The four species examined were exposed to a wide variety of contaminants; however, large spatial and knowledge gaps limit our ability to evaluate if contaminants contribute to species-level declines and if contaminant exposure exacerbates infection by WNS.

Skin pH varies among bat species and seasons and between wild and captive bats

Skin is a key aspect of the immune system in the defence against pathogens. Skin pH regulates the activity of enzymes produced both by hosts and by microbes on host skin, thus implicating pH in disease susceptibility. Skin pH varies inter- and intra-specifically and is influenced by a variety of intrinsic and extrinsic variables. Increased skin alkalinity is associated with a predisposition to cutaneous infections in humans and dogs, and inter-specific and inter-individual variation in skin pH is implicated in differential susceptibility to some skin diseases. The cutaneous pH of bats has not been characterized but is postulated to play a role in susceptibility to white-nose syndrome (WNS), a fungal infection that has decimated several Nearctic bat species. We used non-invasive probes to measure the pH of bat flight membranes in five species with differing susceptibility to WNS. Skin pH ranged from 4.67 to 8.59 and varied among bat species, geographic locations, body parts, age classes, sexes and seasons. Wild Eptesicus fuscus were consistently more acidic than wild Myotis lucifugus, Myotis leibii and Perimyotis subflavus. Juvenile bats had more acidic skin than adults during maternity season but did not differ during swarming. Male M. lucifugus were more acidic than females during maternity season, yet this trend reversed during swarming. Bat skin was more acidic in summer compared to winter, a pattern also reported in humans. Skin pH was more acidic in captive than wild E. fuscus, suggesting environmental impacts on skin pH. The pH of roosting substrates affects skin pH in captive bats and may partially explain seasonal patterns in wild bats that use different roost types across seasons. Future research on the influence of pH on microbial pathogenic factors and skin barrier function may provide valuable insights on new therapeutic targets for treating bat skin conditions.

Counterintuitive scaling between population abundance and local density: Implications for modelling transmission of infectious diseases in bat populations

Models of host-pathogen interactions help to explain infection dynamics in wildlife populations and to predict and mitigate the risk of zoonotic spillover. Insights from models inherently depend on the way contacts between hosts are modelled, and crucially, how transmission scales with animal density. Bats are important reservoirs of zoonotic disease and are among the most gregarious of all mammals. Their population structures can be highly heterogeneous, underpinned by ecological processes across different scales, complicating assumptions regarding the nature of contacts and transmission. Although models commonly parameterise transmission using metrics of total abundance, whether this is an ecologically representative approximation of host-pathogen interactions is not routinely evaluated. We collected a 13-month dataset of tree-roosting Pteropus spp. from 2,522 spatially referenced trees across eight roosts to empirically evaluate the relationship between total roost abundance and tree-level measures of abundance and density-the scale most likely to be relevant for virus transmission. We also evaluate whether roost features at different scales (roost level, subplot level, tree level) are predictive of these local density dynamics. Roost-level features were not representative of tree-level abundance (bats per tree) or tree-level density (bats per m² or m³ ), with roost-level models explaining minimal variation in tree-level measures. Total roost abundance itself was either not a significant predictor (tree-level 3D density) or only weakly predictive (tree-level abundance). This indicates that basic measures, such as total abundance of bats in a roost, may not provide adequate approximations for population dynamics at scales relevant for transmission, and that alternative measures are needed to compare transmission potential between roosts. From the best candidate models, the strongest predictor of local population structure was tree density within roosts, where roosts with low tree density had a higher abundance but lower density of bats (more spacing between bats) per tree. Together, these data highlight unpredictable and counterintuitive relationships between total abundance and local density. More nuanced modelling of transmission, spread and spillover from bats likely requires alternative approaches to integrating contact structure in host-pathogen models, rather than simply modifying the transmission function.

Interspecific variation in evaporative water loss and temperature response, but not metabolic rate, among hibernating bats

Hibernation is widespread among mammals in a variety of environmental contexts. However, few experimental studies consider interspecific comparisons, which may provide insight into general patterns of hibernation strategies. We studied 13 species of free-living bats, including populations spread over thousands of kilometers and diverse habitats. We measured torpid metabolic rate (TMR) and evaporative water loss (two key parameters for understanding hibernation energetics) across a range of temperatures. There was no difference in minimum TMR among species (i.e., all species achieved similarly low torpid metabolic rate) but the temperature associated with minimum TMR varied among species. The minimum defended temperature (temperature below which TMR increased) varied from 8 °C to < 2 °C among species. Conversely, evaporative water loss varied among species, with species clustered in two groups representing high and low evaporative water loss. Notably, species that have suffered population declines due to white-nose syndrome fall in the high evaporative water loss group and less affected species in the low evaporative water loss group. Documenting general patterns of physiological diversity, and associated ecological implications, contributes to broader understanding of biodiversity, and may help predict which species are at greater risk of environmental and anthropogenic stressors.

The scope and severity of white-nose syndrome on hibernating bats in North America

Assessing the scope and severity of threats is necessary for evaluating impacts on populations to inform conservation planning. Quantitative threat assessment often requires monitoring programs that provide reliable data over relevant spatial and temporal scales, yet such programs can be difficult to justify until there is an apparent stressor. Leveraging efforts of wildlife management agencies to record winter counts of hibernating bats, we collated data for 5 species from over 200 sites across 27 U.S. states and 2 Canadian provinces from 1995 to 2018 to determine the impact of white-nose syndrome (WNS), a deadly disease of hibernating bats. We estimated declines of winter counts of bat colonies at sites where the invasive fungus that causes WNS (Pseudogymnoascus destructans) had been detected to assess the threat impact of WNS. Three species undergoing species status assessment by the U.S. Fish and Wildlife Service (Myotis septentrionalis, Myotis lucifugus, and Perimyotis subflavus) declined by more than 90%, which warrants classifying the severity of the WNS threat as extreme based on criteria used by NatureServe. The scope of the WNS threat as defined by NatureServe criteria was large (36% of Myotis lucifugus range) to pervasive (79% of Myotis septentrionalis range) for these species. Declines for 2 other species (Myotis sodalis and Eptesicus fuscus) were less severe but still qualified as moderate to serious based on NatureServe criteria. Data-sharing across jurisdictions provided a comprehensive evaluation of scope and severity of the threat of WNS and indicated regional differences that can inform response efforts at international, national, and state or provincial jurisdictions. We assessed the threat impact of an emerging infectious disease by uniting monitoring efforts across jurisdictional boundaries and demonstrated the importance of coordinated monitoring programs, such as the North American Bat Monitoring Program (NABat), for data-driven conservation assessments and planning.

Bacterial community dynamics on bats and the implications for pathogen resistance

The bats skin microbiota plays an important role in reducing pathogen infection, including the deadly fungal pathogen Pseudogymnoascus destructans, the causative agent of white-nose syndrome. However, the dynamic of skin bacterial communities response to environmental perturbations remains poorly described. We characterized skin bacterial community over time and space in Rhinolophus ferrumequinum, a species with high resistance to the infection with P. destructans. We collected environmental covariate data to determine what factors influenced changes in community structure. We observed significant temporal and spatial shifts in the skin bacterial community, which was mainly associated with variation in operational taxonomic units. The skin bacterial community differed by the environmental microbial reservoirs and was most influenced by host body condition, bat roosting temperature and geographic distance between sites, but was not influenced by pathogen infection. Furthermore, the skin microbiota was enriched in particular taxa with antifungal abilities, such as Enterococcus, Burkholderia, Flavobacterium, Pseudomonas, Corynebacterium and Rhodococcus. And specific strains of Pseudomonas, Corynebacterium and Rhodococcus even inhibited P. destructans growth. Our findings provide new insights in characterizing the variation in bacterial communities can inform us about the processes of driving community assembly and predict the host's ability to resist or survive pathogen infection.

Long-term changes in occurrence, relative abundance, and reproductive fitness of bat species in relation to arrival of White-nose Syndrome in West Virginia, USA

White-nose syndrome (WNS) is a disease caused by the fungus Pseudogymnoascus destructans which has resulted in the deaths of millions of bats across eastern North America. To date, hibernacula counts have been the predominant means of tracking the spread and impact of this disease on bat populations. However, an understanding of the impacts of WNS on demographic parameters outside the winter season is critical to conservation and recovery of bat populations impacted by this disease. We used long-term monitoring data to examine WNS-related impacts to summer populations in West Virginia, where WNS has been documented since 2009. Using capture data from 290 mist-net sites surveyed from 2003 to 2019 on the Monongahela National Forest, we estimated temporal patterns in presence and relative abundance for each bat species. For species that exhibited a population-level response to WNS, we investigated post-WNS changes in adult female reproductive state and body mass. Myotis lucifugus (little brown bat), M. septentrionalis (northern long-eared bat), and Perimyotis subflavus (tri-colored bat) all showed significant decreases in presence and relative abundance during and following the introduction of WNS, while Eptesicus fuscus (big brown bat) and Lasiurus borealis (eastern red bat) responded positively during the WNS invasion. Probability of being reproductively active was not significantly different for any species, though a shift to earlier reproduction was estimated for E. fuscus and M. septentrionalis. For some species, body mass appeared to be influenced by the WNS invasion, but the response differed by species and reproductive state. Results suggest that continued long-term monitoring studies, additional research into impacts of this disease on the fitness of WNS survivors, and a focus on providing optimal nonwintering habitat may be valuable strategies for assessing and promoting recovery of WNS-affected bat populations.

The Milieu Souterrain Superficiel as hibernation habitat for bats: implications for white-nose syndrome

Recent studies have revealed that western populations of little brown bats (Myotis lucifugus) in North America exhibit different hibernation behavior than their eastern counterparts. Understanding these differences is essential for assessing the risk white-nose syndrome (WNS) poses to western bat populations. We used acoustic monitoring and radiotelemetry to study the overwintering behavior of little brown bats near Juneau, Alaska during 2011-2014. Our objectives were to identify the structures they use for hibernation, measure the microclimates within those structures, and determine the timing of immergence and emergence and the length of the hibernation season. We radiotracked 10 little brown bats to underground hibernacula dispersed along two ridge systems. All hibernacula were ≤ 24.2 km from where the bats were captured. Eight bats hibernated in the "Milieu Souterrain Superficiel" (MSS), a network of air-filled underground voids between the rock fragments found in scree (talus) deposits. Two bats hibernated in holes in the soil beneath the root system of a tree or stump (rootball). At least two hibernacula in the MSS were reused in subsequent years. Average MSS and rootball temperatures were warmer and more stable than ambient temperature and were well below the optimal growth range of the fungus that causes WNS. Temperatures in the MSS dropped below freezing, but MSS temperatures increased with depth, indicating bats could avoid subfreezing temperatures by moving deeper into the MSS. Relative humidity (RH) approached 100% in the MSS and under rootballs and was more stable than ambient RH, which also was high, but dropped substantially during periods of extreme cold. Acoustic monitoring revealed that bats hibernated by late October and began emerging by the second week of April; estimates of minimum length of the hibernation season ranged from 156 to 190 days. The cold temperatures, dispersed nature of the hibernacula, and close proximity of hibernacula to summering areas may slow the spread and reduce the impacts of WNS on local populations of little brown bats.

Winter roosting ecology of tricolored bats (Perimyotis subflavus) in trees and bridges

Tricolored bats (Perimyotis subflavus) that roost in subterranean hibernacula have experienced precipitous declines from white-nose syndrome (WNS); however, understudied populations also use during winter non-subterranean roosts such as tree cavities, bridges, and foliage. Our objectives were to determine winter roost use by tricolored bats in an area devoid of subterranean roosts, determine roost microclimates to relate them to growth requirements of the fungal causal agent of WNS, and determine habitat factors influencing winter tree selection. From November to March 2017–2019, we used radiotelemetry to track 15 bats to their day roosts in the upper Coastal Plain of South Carolina and recorded microclimates in accessible tree cavities and bridges. We also characterized habitat and tree characteristics of 24 used trees and 153 random, available trees and used discrete choice models to determine selection. Roost structures included I-beam bridges, cavities in live trees, and foliage. Bridges were warmer and less humid than cavities. Roost temperatures often were amenable to fungal growth (&lt; 19.5°C) but fluctuated widely depending on ambient temperatures. Bats used bridges on colder days (8.7°C ± 5.0 SD) and trees on warmer days (11.3°C ± 5.4). Bats selected low-decay trees closer to streams in areas with high canopy closure and cavity abundance. Bats also appeared to favor hardwood forests and avoid pine forests. Our results suggest that access to multiple roost microclimates might be important for tricolored bats during winter, and forest management practices that retain live trees near streams and foster cavity formation in hardwood forests likely will benefit this species. Our results also suggest tricolored bats using bridge and tree roosts might be less susceptible to WNS than bats using subterranean hibernaculum roosts. Thus, forests in areas without subterranean hibernacula in the southeastern United States that support bats during winter might represent important refugia from WNS for multiple species.

Temperature alone is insufficient to understand hibernation energetics

Energy conservation has long been a focal point in hibernation research. A long-standing assumption is that ambient temperature (Ta) largely defines the rate of energy expenditure because of well-known relationships between Ta, metabolic rate and frequency of arousal from torpor. Body condition and humidity also affect energy expenditure but are usually considered secondary factors. We held tricolored bats (Perimyotis subflavus) in captivity under multiple environmental conditions to directly compare the importance of Ta, fat mass and humidity for hibernation energy expenditure. Fat mass was the best predictor of female mass loss, followed by Ta and humidity. However, males had less fat and adopted a more energetically conservative hibernation strategy. Our results demonstrate that understanding the evolution of behavior, physiology and ecology of hibernation requires disentangling the relative contributions of multiple drivers of hibernation energetics, and that Ta is not always the most important factor driving energy expenditure.