Bat white-nose disease fungus diversity in time and space.

White-nose disease (WND), caused by the psychrophilic fungus Pseudogymnoascusdestructans, represents one of the greatest threats for North American hibernating bats. Research on molecular data has significantly advanced our knowledge of various aspects of the disease, yet more studies are needed regarding patterns of P.destructans genetic diversity distribution. In the present study, we investigate three sites within the native range of the fungus in detail: two natural hibernacula (karst caves) in Bulgaria, south-eastern Europe and one artificial hibernaculum (disused cellar) in Germany, northern Europe, where we conducted intensive surveys between 2014 and 2019. Using 18 microsatellite and two mating type markers, we describe how P.destructans genetic diversity is distributed between and within sites, the latter including differentiation across years and seasons of sampling; across sampling locations within the site; and between bats and hibernaculum walls. We found significant genetic differentiation between hibernacula, but we could not detect any significant differentiation within hibernacula, based on the variables examined. This indicates that most of the pathogen's movement occurs within sites. Genotypic richness of P.destructans varied between sites within the same order of magnitude, being approximately two times higher in the natural caves (Bulgaria) compared to the disused cellar (Germany). Within all sites, the pathogen's genotypic richness was higher in samples collected from hibernaculum walls than in samples collected from bats, which corresponds with the hypothesis that hibernacula walls represent the environmental reservoir of the fungus. Multiple pathogen genotypes were commonly isolated from a single bat (i.e. from the same swab sample) in all study sites, which might be important to consider when studying disease progression.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation.

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals.

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Robust evidence for bats as reservoir hosts is lacking in most African virus studies: a review and call to optimize sampling and conserve bats.

Africa experiences frequent emerging disease outbreaks among humans, with bats often proposed as zoonotic pathogen hosts. We comprehensively reviewed virus-bat findings from papers published between 1978 and 2020 to evaluate the evidence that African bats are reservoir and/or bridging hosts for viruses that cause human disease. We present data from 162 papers (of 1322) with original findings on (1) numbers and species of bats sampled across bat families and the continent, (2) how bats were selected for study inclusion, (3) if bats were terminally sampled, (4) what types of ecological data, if any, were recorded and (5) which viruses were detected and with what methodology. We propose a scheme for evaluating presumed virus-host relationships by evidence type and quality, using the contrasting available evidence for Orthoebolavirus versus Orthomarburgvirus as an example. We review the wording in abstracts and discussions of all 162 papers, identifying key framing terms, how these refer to findings, and how they might contribute to people's beliefs about bats. We discuss the impact of scientific research communication on public perception and emphasize the need for strategies that minimize human-bat conflict and support bat conservation. Finally, we make recommendations for best practices that will improve virological study metadata.

Bat pluripotent stem cells reveal unusual entanglement between host and viruses.

Bats are distinctive among mammals due to their ability to fly, use laryngeal echolocation, and tolerate viruses. However, there are currently no reliable cellular models for studying bat biology or their response to viral infections. Here, we created induced pluripotent stem cells (iPSCs) from two species of bats: the wild greater horseshoe bat (Rhinolophus ferrumequinum) and the greater mouse-eared bat (Myotis myotis). The iPSCs from both bat species showed similar characteristics and had a gene expression profile resembling that of cells attacked by viruses. They also had a high number of endogenous viral sequences, particularly retroviruses. These results suggest that bats have evolved mechanisms to tolerate a large load of viral sequences and may have a more intertwined relationship with viruses than previously thought. Further study of bat iPSCs and their differentiated progeny will provide insights into bat biology, virus host relationships, and the molecular basis of bats' special traits.

Development and Application of Loop-Mediated Isothermal Amplification (LAMP) Assays for Rapid Diagnosis of the Bat White-Nose Disease Fungus Pseudogymnoascus destructans.

Pseudogymnoascus destructans (= Geomyces destructans) is a psychrophilic filamentous fungus that causes White-Nose Disease (WND; the disease associated with White-Nose Syndrome, WNS) in hibernating bats. The disease has caused considerable reductions in bat populations in the USA and Canada since 2006. Identification and detection of the pathogen in pure cultures and environmental samples is routinely based on qPCR or PCR after DNA isolation and purification. Rapid and specific direct detection of the fungus in the field would strongly improve prompt surveillance, and support control measures. Based on the genes coding for ATP citrate lyase1 (acl1) and the 28S-18S ribosomal RNA intergenic spacer (IGS) in P. destructans, two independent LAMP assays were developed for the rapid and sensitive diagnosis of the fungus. Both assays could discriminate P. destructans from 159 tested species of filamentous fungi and yeasts. Sensitivity of the assays was 2.1 picogram per reaction (pg/rxn) and 21 femtogram per reaction (fg/rxn) for the acl1 and IGS based assays, respectively. Moreover, both assays also work with spores and mycelia of P. destructans that are directly added to the master mix without prior DNA extraction. For field-diagnostics, we developed and tested a field-applicable version of the IGS-based LAMP assay. Lastly, we also developed a protocol for preparation of fungal spores and mycelia from swabs and tape liftings of contaminated surfaces or infected bats. This protocol in combination with the highly sensitive IGS-based LAMP-assay enabled sensitive detection of P. destructans from various sources.

Contrasting patterns from two invasion fronts suggest a niche shift of an invasive predator of native bees.

Background: The accuracy of predictions of invasive species ranges is dependent on niche similarity between invasive and native populations and on our ability to identify the niche characteristics. With this work we aimed to compare the niche dynamics of two genetically related invasive populations of Vespa velutina (an effective predator of honeybees and wild pollinators), in two distinct climatic regions, one in central Europe and another one in the north-western Iberian Peninsula, and hence to identify uninvaded regions susceptible to invasion. Methods: Niche dynamics and shifts of V. velutina were assessed by comparing the environmental niches of the native and of the two invasive populations, using climatic, topographic and land use variables. We also ran reciprocal distribution models using different algorithms and records from both native and invasive ranges to compare model predictions and estimate which regions are at a greater risk of being invaded. Results: An apparent niche shift was detected in the population of the NW of Iberian Peninsula, where the species is living under environmental conditions different from the native niche. In central Europe, large suitable areas remain unoccupied. The fact that both invasive populations are well established, despite occupying environmentally distinct regions indicates that V. velutina has a high ability to successfully invade different environmental envelopes from those existing in its native range. For example, in north-western Iberian Peninsula the species is now thriving out of its native niche limits. Moreover, the large extent of still unoccupied environmental space with similar conditions to those used by the species in its native range suggests that there is still a large area of central and eastern Europe that can be potentially invaded by the species.

Population genetics as a tool to elucidate pathogen reservoirs: Lessons from Pseudogymnoascus destructans, the causative agent of White-Nose disease in bats.

Emerging infectious diseases pose a major threat to human, animal, and plant health. The risk of species-extinctions increases when pathogens can survive in the absence of the host. Environmental reservoirs can facilitate this. However, identifying such reservoirs and modes of infection is often highly challenging. In this study, we investigated the presence and nature of an environmental reservoir for the ascomycete fungus Pseudogymnoascus destructans, the causative agent of White-Nose disease. Using 18 microsatellite markers, we determined the genotypic differentiation between 1497 P. destructans isolates collected from nine closely situated underground sites where bats hibernate (i.e., hibernacula) in Northeastern Germany. This approach was unique in that it ensured that every isolate and resulting multilocus genotype was not only present, but also viable and therefore theoretically capable of infecting a bat. The distinct distribution of multilocus genotypes across hibernacula demonstrates that each hibernaculum has an essentially unique fungal population. This would be expected if bats become infected in their hibernaculum (i.e., the site they spend winter in to hibernate) rather than in other sites visited before they start hibernating. In one hibernaculum, both the walls and the hibernating bats were sampled at regular intervals over five consecutive winter seasons (1062 isolates), revealing higher genotypic richness on walls compared to bats and a stable frequency of multilocus genotypes over multiple winters. This clearly implicates hibernacula walls as the main environmental reservoir of the pathogen, from which bats become reinfected annually during the autumn.

Genetic diversity in a long-lived mammal is explained by the past's demographic shadow and current connectivity.

Within-species genetic diversity is crucial for the persistence and integrity of populations and ecosystems. Conservation actions require an understanding of factors influencing genetic diversity, especially in the context of global change. Both population size and connectivity are factors greatly influencing genetic diversity; the relative importance of these factors can, however, change through time. Hence, quantifying the degree to which population size or genetic connectivity are shaping genetic diversity, and at which ecological time scale (past or present), is challenging, yet essential for the development of efficient conservation strategies. In this study, we estimated the genetic diversity of 42 colonies of Rhinolophus hipposideros, a long-lived mammal vulnerable to global change, sampling locations spanning its continental northern range. Here, we present an integrative approach that disentangles and quantifies the contribution of different connectivity measures in addition to contemporary colony size and historic bottlenecks in shaping genetic diversity. In our study, the best model explained 64% of the variation in genetic diversity. It included historic bottlenecks, contemporary colony size, connectivity and a negative interaction between the latter two. Contemporary connectivity explained most genetic diversity when considering a 65 km radius around the focal colonies, emphasizing the large geographic scale at which the positive impact of connectivity on genetic diversity is most profound and hence, the minimum scale at which conservation should be planned. Our results highlight that the relative importance of the two main factors shaping genetic diversity varies through time, emphasizing the relevance of disentangling them to ensure appropriate conservation strategies.

Heterothermy and antifungal responses in bats.

Hibernation, a period where bats have suppressed immunity and low body temperatures, provides the psychrophilic fungus Pseudogymnoascus destructans the opportunity to colonise bat skin, leading to severe disease in susceptible species. Innate immunity, which requires less energy and may remain more active during torpor, can control infections with local inflammation in some bat species that are resistant to infection. If infection is not controlled before emergence from hibernation, ineffective adaptive immune mechanisms are activated, including incomplete Th1, ineffective Th2, and variable Th17 responses. The Th17 and neutrophil responses, normally beneficial antifungal mechanisms, appear to be sources of immunopathology for susceptible bat species, because they are hyperactivated after return to homeothermy. Non-susceptible species show both well-balanced and suppressed immune responses both during and after hibernation.

Determinants of defence strategies of a hibernating European bat species towards the fungal pathogen Pseudogymnoascus destructans.

Pseudogymnoascus destructans (Pd), the causative agent of white-nose syndrome in North America, has decimated bat populations within a decade. The fungus impacts bats during hibernation when physiological functions, including immune responses, are down-regulated. Studies have shown that Pd is native to Europe, where it is not associated with mass mortalities. Moreover, genomic and proteomic studies indicated that European bats may have evolved an effective immune defence, which is lacking in North American bats. However, it is still unclear which defence strategy enables European bats to cope with the pathogen. Here, we analyzed selected physiological and immunological parameters in torpid, Pd infected European greater mouse-eared bats (Myotis myotis) showing three different levels of infection (asymptomatic, mild and severe symptoms). From a subset of the studied bats we tracked skin temperatures during one month of hibernation. Contrasting North American bats, arousal patterns remained unaffected by Pd infections in M. myotis. In general, heavier M. myotis aroused more often from hibernation and showed less severe disease symptoms than lean individuals; most likely because heavy bats were capable of reducing the Pd load more effectively than lean individuals. In the blood of severely infected bats, we found higher gene expression levels of an inflammatory cytokine (IL-1ß), but lower levels of an acute phase protein (haptoglobin), reactive oxygen metabolites (ROMs) and plasma non-enzymatic antioxidant capacity (OXY) compared to conspecifics with lower levels of infection. We conclude that M. myotis, and possibly also other European bat species, tolerate Pd infections during torpor by using selected acute phase response parameters at baseline levels, yet without arousing from torpor and without synthesizing additional immune molecules.

Timescale and colony-dependent relationships between environmental conditions and plasma oxidative markers in a long-lived bat species.

To increase the applicability and success of physiological approaches in conservation plans, conservation physiology should be based on ecologically relevant relationships between physiological markers and environmental variation that can only be obtained from wild populations. Given their integrative and multifaceted aspects, markers of oxidative status have recently been considered in conservation physiology, but still need to be validated across environmental conditions and locations. Here, we examined whether inter-annual variation in two oxidative markers, plasma antioxidant capacity and plasma hydroperoxides, followed inter-annual variation in temperature anomalies and associated vegetation changes in four colonies of long-lived greater mouse-eared bats (Myotis myotis) monitored over five consecutive years. We found that the plasma antioxidant capacity of bats decreased while plasma hydroperoxide concentrations increased with increasing temperature anomalies occurring in the two weeks before blood sampling. Moreover, the antioxidant defences of these bats reflected vegetation indices, which themselves reflected the thermal conditions experienced by bats in their foraging habitat. Variation in oxidative markers therefore appears to be due to variation in thermoregulatory costs and to indirect changes in foraging costs. Overall, these results validate the use of markers of oxidative status in conservation physiology to monitor thermal perturbations recently experienced by animals in their natural habitat. However, even though oxidative markers varied in the same direction in all four bat colonies across years, the amplitude of their response differed. If these different physiological responses reflect different performances (e.g. productivity, survival rate) between colonies, this implies that, if necessary, conservation measures may need to be applied at the local scale.

Six reference-quality genomes reveal evolution of bat adaptations.

Bats possess extraordinary adaptations, including flight, echolocation, extreme longevity and unique immunity. High-quality genomes are crucial for understanding the molecular basis and evolution of these traits. Here we incorporated long-read sequencing and state-of-the-art scaffolding protocols1 to generate, to our knowledge, the first reference-quality genomes of six bat species (Rhinolophus ferrumequinum, Rousettus aegyptiacus, Phyllostomus discolor, Myotis myotis, Pipistrellus kuhlii and Molossus molossus). We integrated gene projections from our 'Tool to infer Orthologs from Genome Alignments' (TOGA) software with de novo and homology gene predictions as well as short- and long-read transcriptomics to generate highly complete gene annotations. To resolve the phylogenetic position of bats within Laurasiatheria, we applied several phylogenetic methods to comprehensive sets of orthologous protein-coding and noncoding regions of the genome, and identified a basal origin for bats within Scrotifera. Our genome-wide screens revealed positive selection on hearing-related genes in the ancestral branch of bats, which is indicative of laryngeal echolocation being an ancestral trait in this clade. We found selection and loss of immunity-related genes (including pro-inflammatory NF-κB regulators) and expansions of anti-viral APOBEC3 genes, which highlights molecular mechanisms that may contribute to the exceptional immunity of bats. Genomic integrations of diverse viruses provide a genomic record of historical tolerance to viral infection in bats. Finally, we found and experimentally validated bat-specific variation in microRNAs, which may regulate bat-specific gene-expression programs. Our reference-quality bat genomes provide the resources required to uncover and validate the genomic basis of adaptations of bats, and stimulate new avenues of research that are directly relevant to human health and disease1.

Seasonal patterns of Pseudogymnoascus destructans germination indicate host-pathogen coevolution.

Emerging infectious diseases rank among the most important threats to human and wildlife health. A comprehensive understanding of the mode of infection and presence of potential reservoirs is critical for the development of effective counter strategies. Fungal pathogens can remain viable in environmental reservoirs for extended periods of time before infecting susceptible individuals. This may be the case for Pseudogymnoascus destructans (Pd), the causative agent of bat white-nose disease. Owing to its cold-loving nature, this fungal pathogen only grows on bats during hibernation, when their body temperature is reduced. Bats only spend part of their life cycle in hibernation and do not typically show signs of infection in summer, raising the question of whether Pd remains viable in hibernacula during this period (roughly six months). If so, this could facilitate the re-infection of bats when they return to the sites the following winter. In a laboratory experiment, we determined the germination rate of Pd spores kept under constant conditions on a wall-like substrate, over the course of two years. Results showed that the seasonal pattern in Pd germination mirrored the life cycle of the bats, with an increased germination rate at times when hibernating bats would naturally be present and lower germination rates during their absence. We suggest that Pd is dependent on the presence of hibernating bats and has therefore coupled its germination rate to host availability. Furthermore, we demonstrate that Pd spores survive extended periods of host absence and can remain viable for at least two years. There is, however, a strong decrease in spore viability between the first and second years (98%). Pd viability for at least two years on a solid mineral-based substrate establishes the potential for environmental reservoirs in hibernacula walls and has strong implications for the efficacy of certain management strategies (e.g. bat culling).

Drivers of longitudinal telomere dynamics in a long-lived bat species, Myotis myotis.

Age-related telomere shortening is considered a hallmark of the ageing process. However, a recent cross-sectional ageing study of relative telomere length (rTL) in bats failed to detect a relationship between rTL and age in the long-lived genus Myotis (M. myotis and M. bechsteinii), suggesting some other factors are responsible for driving telomere dynamics in these species. Here, we test if longitudinal rTL data show signatures of age-associated telomere attrition in M. myotis and differentiate which intrinsic or extrinsic factors are likely to drive telomere length dynamics. Using quantitative polymerase chain reaction, rTL was measured in 504 samples from a marked population, from Brittany, France, captured between 2013 and 2016. These represent 174 individuals with an age range of 0 to 7+ years. We find no significant relationship between rTL and age (p = .762), but demonstrate that within-individual rTL is highly variable from year to year. To investigate the heritability of rTL, a population pedigree (n = 1744) was constructed from genotype data generated from a 16-microsatellite multiplex, designed from an initial, low-coverage, Illumina genome for M. myotis. Heritability was estimated in a Bayesian, mixed model framework, and showed that little of the observed variance in rTL is heritable (h2  = 0.01-0.06). Rather, correlations of first differences, correlating yearly changes in telomere length and weather variables, demonstrate that, during the spring transition, average temperature, minimum temperature, rainfall and windspeed correlate with changes in longitudinal telomere dynamics. As such, rTL may represent a useful biomarker to quantify the physiological impact of various environmental stressors in bats.

Longitudinal comparative transcriptomics reveals unique mechanisms underlying extended healthspan in bats.

Bats are the longest-lived mammals, given their body size. However, the underlying molecular mechanisms of their extended healthspans are poorly understood. To address this question we carried out an eight-year longitudinal study of ageing in long-lived bats (Myotis myotis). We deep-sequenced ~1.7 trillion base pairs of RNA from 150 blood samples collected from known aged bats to ascertain the age-related transcriptomic shifts and potential microRNA-directed regulation that occurred. We also compared ageing transcriptomic profiles between bats and other mammals by analysis of 298 longitudinal RNA sequencing datasets. Bats did not show the same transcriptomic changes with age as commonly observed in humans and other mammals, but rather exhibited a unique, age-related gene expression pattern associated with DNA repair, autophagy, immunity and tumour suppression that may drive their extended healthspans. We show that bats have naturally evolved transcriptomic signatures that are known to extend lifespan in model organisms, and identify novel genes not yet implicated in healthy ageing. We further show that bats' longevity profiles are partially regulated by microRNA, thus providing novel regulatory targets and pathways for future ageing intervention studies. These results further disentangle the ageing process by highlighting which ageing pathways contribute most to healthy ageing in mammals.

Considering adaptive genetic variation in climate change vulnerability assessment reduces species range loss projections.

Local adaptations can determine the potential of populations to respond to environmental changes, yet adaptive genetic variation is commonly ignored in models forecasting species vulnerability and biogeographical shifts under future climate change. Here we integrate genomic and ecological modeling approaches to identify genetic adaptations associated with climate in two cryptic forest bats. We then incorporate this information directly into forecasts of range changes under future climate change and assessment of population persistence through the spread of climate-adaptive genetic variation (evolutionary rescue potential). Considering climate-adaptive potential reduced range loss projections, suggesting that failure to account for intraspecific variability can result in overestimation of future losses. On the other hand, range overlap between species was projected to increase, indicating that interspecific competition is likely to play an important role in limiting species' future ranges. We show that although evolutionary rescue is possible, it depends on a population's adaptive capacity and connectivity. Hence, we stress the importance of incorporating genomic data and landscape connectivity in climate change vulnerability assessments and conservation management.

Phenotypic plasticity closely linked to climate at origin and resulting in increased mortality under warming and frost stress in a common grass.

Phenotypic plasticity is important for species responses to global change and species coexistence. Phenotypic plasticity differs among species and traits and changes across environments. Here, we investigated phenotypic plasticity of the widespread grass Arrhenatherum elatius in response to winter warming and frost stress by comparing phenotypic plasticity of 11 geographically and environmentally distinct populations of this species to phenotypic plasticity of populations of different species originating from a single environment. The variation in phenotypic plasticity was similar for populations of a single species from different locations compared to populations of functionally and taxonomically diverse species from one environment for the studied traits (leaf biomass production and root integrity after frost) across three indices of phenotypic plasticity (RDPI, PIN, slope of reaction norm). Phenotypic plasticity was not associated with neutral genetic diversity but closely linked to the climate of the populations' origin. Populations originating from warmer and more variable climates showed higher phenotypic plasticity. This indicates that phenotypic plasticity can itself be considered as a trait subject to local adaptation to climate. Finally, our data emphasize that high phenotypic plasticity is not per se positive for adaptation to climate change, as differences in stress responses are resulting in high phenotypic plasticity as expressed by common plasticity indices, which is likely to be related to increased mortality under stress in more plastic populations.

Population level mitogenomics of long-lived bats reveals dynamic heteroplasmy and challenges the Free Radical Theory of Ageing.

Bats are the only mammals capable of true, powered flight, which drives an extremely high metabolic rate. The "Free Radical Theory of Ageing" (FTRA) posits that a high metabolic rate causes mitochondrial heteroplasmy and the progressive ageing phenotype. Contrary to this, bats are the longest-lived order of mammals given their small size and high metabolic rate. To investigate if bats exhibit increased mitochondrial heteroplasmy with age, we performed targeted, deep sequencing of mitogenomes and measured point heteroplasmy in wild, long lived Myotis myotis. Blood was sampled from 195 individuals, aged between <1 and at 6+ years old, and whole mitochondria deep-sequenced, with a subset sampled over multiple years. The majority of heteroplasmies were at a low frequency and were transitions. Oxidative mutations were present in only a small number of individuals, suggesting local oxidative stress events. Cohort data showed no significant increase in heteroplasmy with age, while longitudinal data from recaptured individuals showed heteroplasmy is dynamic, and does not increase uniformly over time. We show that bats do not suffer from the predicted, inevitable increase in heteroplasmy as posited by the FRTA, instead heteroplasmy was found to be dynamic, questioning its presumed role as a primary driver of ageing.