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.

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.

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.

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.

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.

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.

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.

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.

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).

Evidence for genetic variation in Natterer's bats (Myotis nattereri) across three regions in Germany but no evidence for co-variation with their associated astroviruses.

BACKGROUND: As bats have recently been described to harbor many different viruses, several studies have investigated the genetic co-variation between viruses and different bat species. However, little is known about the genetic co-variation of viruses and different populations of the same bat species, although such information is needed for an understanding of virus transmission dynamics within a given host species. We hypothesized that if virus transmission between host populations depends on events linked to gene flow in the bats, genetic co-variation should exist between host populations and astroviruses. RESULTS: We used 19 nuclear and one mitochondrial microsatellite loci to analyze the genetic population structure of the Natterer's bat (Myotis nattereri) within and among populations at different geographical scales in Germany. Further, we correlated the observed bat population structure to that of partial astrovirus sequences (323-394 nt fragments of the RNA-dependent RNA polymerase gene) obtained from the same bat populations. Our analyses revealed that the studied bat colonies can be grouped into three distinct genetic clusters, corresponding to the three geographic regions sampled. Furthermore, we observed an overall isolation-by-distance pattern, while no significant pattern was observed within a geographic region. Moreover, we found no correlation between the genetic distances among the bat populations and the astrovirus sequences they harbored. Even though high genetic similarity of some of the astrovirus haplotypes found in several different regions was detected, identical astrovirus haplotypes were not shared between different sampled regions. CONCLUSIONS: The genetic population structure of the bat host suggests that mating sites where several local breeding colonies meet act as stepping-stones for gene flow. Identical astrovirus haplotypes were not shared between different sampled regions suggesting that astroviruses are mostly transmitted among host colonies at the local scale. Nevertheless, high genetic similarity of some of the astrovirus haplotypes found in several different regions implies that occasional transmission across regions with subsequent mutations of the virus haplotypes does occur.

Which temporal resolution to consider when investigating the impact of climatic data on population dynamics? The case of the lesser horseshoe bat (Rhinolophus hipposideros).

Climatic variables are often considered when studying environmental impacts on population dynamics of terrestrial species. However, the temporal resolution considered varies depending on studies, even among studies of the same taxa. Most studies interested in climatic impacts on populations tend to average climatic data across timeframes covering life cycle periods of the organism in question or longer, even though most climatic databases provide at least a monthly resolution. We explored the impact of climatic variables on lesser horseshoe bat (Rhinolophus hipposideros) demography based on count data collected at 94 maternity colonies from 2000 to 2014 in Britanny, France. Meteorological data were considered using different time resolutions (month, life cycle period and year) to investigate their adequacy. Model averaging was used to detect significant predictors for each temporal resolution. Our results show that the finest temporal resolution, e.g. month, was more informative than coarser ones. Precipitation predictors were particularly decisive, with a negative impact on colony sizes when rainfall occurred in October, and a positive impact for June precipitations. Fecundity was influenced by April weather. This highlights the strong impact of climatic conditions during crucial but short time periods on the population dynamics of bats. We demonstrate the importance of choosing an appropriate time resolution and suggest that analogous studies should consider fine-scale temporal resolution (e.g. month) to better grasp the relationship between population dynamics and climatic conditions.

How and why overcome the impediments to resolution: lessons from rhinolophid and hipposiderid bats.

The phylogenetic and taxonomic relationships among the Old World leaf-nosed bats (Hipposideridae) and the closely related horseshoe bats (Rhinolophidae) remain unresolved. In this study, we generated a novel approximately 10-kb molecular data set of 19 nuclear exon and intron gene fragments for 40 bat species to elucidate the phylogenetic relationships within the families Rhinolophidae and Hipposideridae. We estimated divergence times and explored potential reasons for any incongruent phylogenetic signal. We demonstrated the effects of outlier taxa and genes on phylogenetic reconstructions and compared the relative performance of intron and exon data to resolve phylogenetic relationships. Phylogenetic analyses produced a well-resolved phylogeny, supporting the familial status of Hipposideridae and demonstrated the paraphyly of the largest genus, Hipposideros. A fossil-calibrated timetree and biogeographical analyses estimated that Rhinolophidae and Hipposideridae diverged in Africa during the Eocene approximately 42 Ma. The phylogram, the timetree, and a unique retrotransposon insertion supported the elevation of the subtribe Rhinonycterina to family level and which is diagnosed herein. Comparative analysis of diversification rates showed that the speciose genera Rhinolophus and Hipposideros underwent diversification during the Mid-Miocene Climatic Optimum. The intron versus exon analyses demonstrated the improved nodal support provided by introns for our optimal tree, an important finding for large-scale phylogenomic studies, which typically rely on exon data alone. With the recent outbreak of Middle East respiratory syndrome, caused by a novel coronavirus, the study of these species is urgent as they are considered the natural reservoir for emergent severe acute respiratory syndrome (SARS)-like coronaviruses. It has been shown that host phylogeny is the primary factor that determines a virus's persistence, replicative ability, and can act as a predictor of new emerging disease. Therefore, this newly resolved phylogeny can be used to direct future assessments of viral diversity and to elucidate the origin and development of SARS-like coronaviruses in mammals.

Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: Lessons from horseshoe bats (Rhinolophidae: Chiroptera).

Despite many studies illustrating the perils of utilising mitochondrial DNA in phylogenetic studies, it remains one of the most widely used genetic markers for this purpose. Over the last decade, nuclear introns have been proposed as alternative markers for phylogenetic reconstruction. However, the resolution capabilities of mtDNA and nuclear introns have rarely been quantified and compared. In the current study we generated a novel ∼5kb dataset comprising six nuclear introns and a mtDNA fragment. We assessed the relative resolution capabilities of the six intronic fragments with respect to each other, when used in various combinations together, and when compared to the traditionally used mtDNA. We focused on a major clade in the horseshoe bat family (Afro-Palaearctic clade; Rhinolophidae) as our case study. This old, widely distributed and speciose group contains a high level of conserved morphology. This morphological stasis renders the reconstruction of the phylogeny of this group with traditional morphological characters complex. We sampled multiple individuals per species to represent their geographic distributions as best as possible (122 individuals, 24 species, 68 localities). We reconstructed the species phylogeny using several complementary methods (partitioned Maximum Likelihood and Bayesian and Bayesian multispecies-coalescent) and made inferences based on consensus across these methods. We computed pairwise comparisons based on Robinson-Foulds tree distance metric between all Bayesian topologies generated (27,000) for every gene(s) and visualised the tree space using multidimensional scaling (MDS) plots. Using our supported species phylogeny we estimated the ancestral state of key traits of interest within this group, e.g. echolocation peak frequency which has been implicated in speciation. Our results revealed many potential cryptic species within this group, even in taxa where this was not suspected a priori and also found evidence for mtDNA introgression. We demonstrated that by using just two introns one can recover a better supported species tree than when using the mtDNA alone, despite the shorter overall length of the combined introns. Additionally, when combining any single intron with mtDNA, we showed that the result is highly similar to the mtDNA gene tree and far from the true species tree and therefore this approach should be avoided. We caution against the indiscriminate use of mtDNA in phylogenetic studies and advocate for pilot studies to select nuclear introns. The selection of marker type and number is a crucial step that is best based on critical examination of preliminary or previously published data. Based on our findings and previous publications, we recommend the following markers to recover phylogenetic relationships between recently diverged taxa (<20 My) in bats and other mammals: ACOX2, COPS7A, BGN, ROGDI and STAT5A.

Circum-Mediterranean phylogeography of a bat coupled with past environmental niche modeling: A new paradigm for the recolonization of Europe?

The isolation of populations in the Iberian, Italian and Balkan peninsulas during the ice ages define four main paradigms that explain much of the known distribution of intraspecific genetic diversity in Europe. In this study we investigated the phylogeography of a wide-spread bat species, the bent-winged bat, Miniopterus schreibersii around the Mediterranean basin and in the Caucasus. Environmental Niche Modeling (ENM) analysis was applied to predict both the current distribution of the species and its distribution during the last glacial maximum (LGM). The combination of genetics and ENM results suggest that the populations of M. schreibersii in Europe, the Caucasus and Anatolia went extinct during the LGM, and the refugium for the species was a relatively small area to the east of the Levantine Sea, corresponding to the Mediterranean coasts of present-day Syria, Lebanon, Israel, and northeastern and northwestern Egypt. Subsequently the species first repopulated Anatolia, diversified there, and afterwards expanded into the Caucasus, continental Europe and North Africa after the end of the LGM. The fossil record in Iberia and the ENM results indicate continuous presence of Miniopterus in this peninsula that most probably was related to the Maghrebian lineage during the LGM, which did not persist afterwards. Using our results combined with similar findings in previous studies, we propose a new paradigm explaining the general distribution of genetic diversity in Europe involving the recolonization of the continent, with the main contribution from refugial populations in Anatolia and the Middle East. The study shows how genetics and ENM approaches can complement each other in providing a more detailed picture of intraspecific evolution.

Patterns of orofacial clefting in the facial morphology of bats: a possible naturally occurring model of cleft palate.

A normal feature of the facial anatomy of many species of bat is the presence of bony discontinuities or clefts, which bear a remarkable similarity to orofacial clefts that occur in humans as a congenital pathology. These clefts occur in two forms: a midline cleft between the two premaxillae (analogous to the rare midline craniofacial clefts in humans) and bilateral paramedian clefts between the premaxilla and the maxillae (analogous to the typical cleft lip and palate in humans). Here, we describe the distribution of orofacial clefting across major bat clades, exploring the relationship of the different patterns of clefting to feeding mode, development of the vomeronasal organ, development of the nasolacrimal duct and mode of emission of the echolocation call in different bat groups. We also present the results of detailed radiographic and soft tissue dissections of representative examples of the two types of cleft. The midline cleft has arisen independently multiple times in bat phylogeny, whereas the paramedian cleft has arisen once and is a synapomorphy uniting the Rhinolophidae and Hipposideridae. In all cases examined, the bony cleft is filled in by a robust fibrous membrane, continuous with the periosteum of the margins of the cleft. In the paramedian clefts, this membrane splits to enclose the premaxilla but forms a loose fold laterally between the premaxilla and maxilla, allowing the premaxilla and nose-leaf to pivot dorsoventrally in the sagittal plane under the action of facial muscles attached to the nasal cartilages. It is possible that this is a specific adaptation for echolocation and/or aerial insectivory. Given the shared embryological location of orofacial clefts in bats and humans, it is likely that aspects of the developmental control networks that produce cleft lip and palate in humans may also be implicated in the formation of these clefts as a normal feature in some bats. A better understanding of craniofacial development in bats with and without clefts may therefore suggest avenues for research into abnormal craniofacial development in humans.

Will reduced host connectivity curb the spread of a devastating epidemic?

The white-nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans, is threatening the cave-dwelling bat fauna of North America by killing individuals by the thousands in hibernacula each winter since its appearance in New York State less than ten years ago. Epidemiological models predict that WNS will reach the western coast of the USA by 2035, potentially eliminating most populations of susceptible bat species in its path (Frick et al. 2015; O'Regan et al. 2015). These models were built and validated using distributional data from the early years of the epidemic, which spread throughout eastern North America following a route driven by cave density and winter severity (Maher et al. 2012). In this issue of Molecular Ecology, Wilder et al. (2015) refine these findings by showing that connectivity among host populations, as assessed by population genetic markers, is crucial in determining the spread of the pathogen. Because host connectivity is much reduced in the hitherto disease free western half of North America, Wilder et al. make the reassuring prediction that the disease will spread more slowly west of the Great Plains.

A new species of the Miniopterus schreibersii species complex (Chiroptera: Miniopteridae) from the Maghreb Region, North Africa.

We used an integrative approach combining cranio-dental characters, mitochondrial and nuclear data and acoustic data to show the presence in the genus Miniopterus of a cryptic species from the Maghreb region. This species was previously recognised as Miniopterus schreibersii (Kuhl, 1817). Miniopterus maghrebensis sp. nov. can be differentiated from M. schreibersii sensu stricto on the basis of cranial characters and from mitochondrial DNA and microsatellite evidence. Although slight external morphological and acoustic differences were noted between the two species, these criteria alone did not allow reliable species identification from live animals. Based on the specimens identified morphologically and/or genetically, the distribution range of M. maghrebensis sp. nov. extends from northern Morocco to south of the High Atlas Mountains and northern Tunisia. The new cryptic species is found in sympatry with M. schreibersii s.str. near coastal regions of North Africa.

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.