Genetic connectivity among swarming sites in the wide ranging and recently declining little brown bat (Myotis lucifugus)

Prehibernation swarming in temperate bats: a critical transition between summer activity and hibernation

In this contribution to Dr. Brock Fenton's Festschrift, we briefly reflect on Dr. Fenton's seminal works examining bat swarming behaviour in Ontario and use these reflections as a launch pad to conduct a global review on autumn swarming in bats, and underlying hypotheses to explain this behaviour. Our review frames the swarming period as a time of critical transitions, during which bats must balance multiple life history trade-offs, and we consider how various intrinsic and extrinsic factors may contribute to inter- and intraspecific differences in autumn behaviour. We discuss the transition away from summer residency, including maternity colony breakup, day roosting, and migration during autumn. We review key life history elements of swarming, including mating behaviours and associated reproductive condition, genetic exchange during swarming, and variation among sexes, ages, and species. Finally, we discuss the behaviours and physiological states of bats transitioning from the swarming period to hibernation. Throughout, we identify common patterns and also exceptions. Over 50 years of research has yielded many insights into autumn swarming, but knowledge gaps remain. Future research focus on a greater diversity of species will reveal general principles underlying the transition from summer active season, through the swarming period, and into winter hibernation.

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.

Similar hibernation physiology in bats across broad geographic ranges

Species with broad geographic ranges may experience varied environmental conditions throughout their range leading to local adaptation. Variation among populations reflects potential adaptability or plasticity, with implications for populations impacted by disease, climate change, and other anthropogenic influences. However, behavior may counteract divergent selection among populations. We studied intraspecific variation in hibernation physiology of Myotis lucifugus (little brown myotis) and Corynorhinus townsendii (Townsend's big-eared bat), two species of bats with large geographic ranges. We studied M. lucifugus at three hibernacula which spanned a latitudinal gradient of 1500 km, and C. townsendii from 6 hibernacula spread across 1200 km latitude and 1200 km longitude. We found no difference in torpid metabolic rate among populations of either species, nor was there a difference in the effect of ambient temperature among sites. Evaporative water loss was similar among populations of both species, with the exception of one C. townsendii pairwise site difference and one M. lucifugus site that differed from the others. We suggest the general lack of geographic variation is a consequence of behavioral microhabitat selection. As volant animals, bats can travel relatively long distances in search of preferred microclimates for hibernation. Despite dramatic macroclimate differences among populations, hibernating bats are able to find preferred microclimate conditions within their range, resulting in similar selection pressures among populations spread across wide geographic ranges.

Landscape Genetic Connectivity and Evidence for Recombination in the North American Population of the White-Nose Syndrome Pathogen, Pseudogymnoascus destructans

White-Nose Syndrome is an ongoing fungal epizootic caused by epidermal infections of the fungus, Pseudogymnoascus destructans (P. destructans), affecting hibernating bat species in North America. Emerging early in 2006 in New York State, infections of P. destructans have spread to 38 US States and seven Canadian Provinces. Since then, clonal isolates of P. destructans have accumulated genotypic and phenotypic variations in North America. Using microsatellite and single nucleotide polymorphism markers, we investigated the population structure and genetic relationships among P. destructans isolates from diverse regions in North America to understand its pattern of spread, and to test hypotheses about factors that contribute to transmission. We found limited support for genetic isolation of P. destructans populations by geographic distance, and instead identified evidence for gene flow among geographic regions. Interestingly, allelic association tests revealed evidence for recombination in the North American P. destructans population. Our landscape genetic analyses revealed that the population structure of P. destructans in North America was significantly influenced by anthropogenic impacts on the landscape. Our results have important implications for understanding the mechanism(s) of P. destructans spread.

Population structure, gene flow and relatedness of Natterer’s bats in Northern England

There have been significant declines in population numbers of many bat species in the United Kingdom, including Natterer’s batsMyotis nattereri, over the last century, largely due to anthropogenic changes. The philopatry, which temperate-zone bats often exhibit to their natal landscapes, in combination with anthropogenic threats, can lead to fragmentation, isolation and sub-division of populations. This may result in bottlenecks and declines in genetic diversity. Multi-scaled research is required to disentangle how the variation in the physical traits of bat species (e.g. affecting flight), as well as their social and behavioural traits (e.g. community size, migration, breeding systems), may affect the genetic health of populations and provide a potential buffer against fragmentation. We used microsatellite markers to characterise the genetic diversity and population structure present in Natterer’s bat colonies to determine whether summer roosting bat colonies were spatially differentiated or part of a meta-population. Analyses of population structure and measures of genetic relatedness suggest spatially differentiated populations of bats exhibit long term site fidelity to summer roosting sites, whilst high genetic diversity at sites indicates gene exchange occurs via swarming sites. Natterer’s bats in northern England may travel greater distances to swarming sites than has been previously documented.

Population Connectivity Predicts Vulnerability to White-Nose Syndrome in the Chilean Myotis (Myotis chiloensis) - A Genomics Approach

Despite its peculiar distribution, the biology of the southernmost bat species in the world, the Chilean myotis (Myotis chiloensis), has garnered little attention so far. The species has a north-south distribution of c. 2800 km, mostly on the eastern side of the Andes mountain range. Use of extended torpor occurs in the southernmost portion of the range, putting the species at risk of bat white-nose syndrome, a fungal disease responsible for massive population declines in North American bats. Here, we examined how geographic distance and topology would be reflected in the population structure of M. chiloensis along the majority of its range using a double digestion RAD-seq method. We sampled 66 individuals across the species range and discovered pronounced isolation-by-distance. Furthermore, and surprisingly, we found higher degrees of heterozygosity in the southernmost populations compared to the north. A coalescence analysis revealed that our populations may still not have reached secondary contact after the Last Glacial Maximum. As for the potential spread of pathogens, such as the fungus causing WNS, connectivity among populations was noticeably low, especially between the southern hibernatory populations in the Magallanes and Tierra del Fuego, and more northerly populations. This suggests the probability of geographic spread of the disease from the north through bat-to-bat contact to susceptible populations is low. The study presents a rare case of defined population structure in a bat species and warrants further research on the underlying factors contributing to this. See the graphical abstract here. https://doi.org/10.25387/g3.12173385

Genome-Wide Changes in Genetic Diversity in a Population of Myotis lucifugus Affected by White-Nose Syndrome

Novel pathogens can cause massive declines in populations, and even extirpation of hosts. But disease can also act as a selective pressure on survivors, driving the evolution of resistance or tolerance. Bat white-nose syndrome (WNS) is a rapidly spreading wildlife disease in North America. The fungus causing the disease invades skin tissues of hibernating bats, resulting in disruption of hibernation behavior, premature energy depletion, and subsequent death. We used whole-genome sequencing to investigate changes in allele frequencies within a population of Myotis lucifugus in eastern North America to search for genetic resistance to WNS. Our results show low F_ST values within the population across time, i.e., prior to WNS (Pre-WNS) compared to the population that has survived WNS (Post-WNS). However, when dividing the population with a geographical cut-off between the states of Pennsylvania and New York, a sharp increase in values on scaffold GL429776 is evident in the Post-WNS samples. Genes present in the diverged area are associated with thermoregulation and promotion of brown fat production. Thus, although WNS may not have subjected the entire M. lucifugus population to selective pressure, it may have selected for specific alleles in Pennsylvania through decreased gene flow within the population. However, the persistence of remnant sub-populations in the aftermath of WNS is likely due to multiple factors in bat life history.

Observed Resiliency of Little Brown Myotis to Long-Term White-Nose Syndrome Exposure

White-nose syndrome (WNS) is a disease that has killed millions of bats in eastern North America and has steadily been spreading across the continent. Little brown myotis Myotis lucifugus populations have experienced extensive declines; however, some localized populations have remained resilient, with bats surviving multiple years past initial WNS exposure. These persistent populations may be critical to species recovery, and understanding mechanisms leading to this long-term survival and persistence may provide insight into overall bat and disease management. We monitored a maternity colony of little brown myotis on Fort Drum Military Installation in northern New York between 2006 and 2017 to determine basic demographic parameters and find evidence of what may be leading to resiliency and persistence at this site. Total colony size declined by approximately 88% from 2008 to 2010 due primarily to impacts of WNS. Counts of all adults returning to the colony stabilized during 2010–2014 (mean = 94, range 84–101) and increased after 2014 (mean = 132, range = 108–166). We captured 727 little brown myotis (575 females, 152 males) and banded 534 individuals (389 females, 145 males) at the colony. The majority of sampled bats showed evidence of recent past WNS infection and exposure to Pseudogymnoascus destructans, and we documented pervasive presence and limited viability of the fungus within the colony's main roosting structure. We recaptured 98 individually marked females in years after initial banding, and some individuals survived at least 6 y. Ninety-one percent of all adult females, 93% of recaptured bats, and 90% of 1-y-old females (i.e., bats recaptured the first year after initial capture as juveniles) showed evidence of reproduction during the monitoring period. Using mark–recapture models, we estimated annual survival rates of juvenile and adult little brown myotis during 2009–2016 and examined whether reproductive condition or evidence of recent infection of WNS had any effect on survival. Annual survival rates were similar between juveniles and adults, but highly variable, ranging from 41.0 to 86.5%. Models indicated that neither evidence of recent past exposure to WNS nor reproductive status were related to survival. No one parameter stood out as being responsible for this colony's continued existence, and it is likely that many interwoven factors were responsible for the observed resilience. Although relatively high reproductive effort from all females (i.e., both1-y-old and >1-y-old ) and intermittently suitable survival rates have led to the continued persistence of, and population increases in, this summer colony, mortality from WNS and inherently low reproductive potential still seemed to be limiting population growth. Until there is a better understanding of this overall potential resiliency in little brown myotis, we recommend considering minimizing disturbance and direct human involvement within these persisting populations to allow whatever natural recovery that may be occurring to evolve uninterrupted.

Prelude to a panzootic: Gene flow and immunogenetic variation in northern little brown myotis vulnerable to bat white-nose syndrome

The fungus that causes bat white-nose syndrome (WNS) recently leaped from eastern North America to the Pacific Coast. The pathogen’s spread is associated with the genetic population structure of a host ( Myotis lucifugus). To understand the fine-scale neutral and immunogenetic variation among northern populations of M. lucifugus, we sampled 1142 individuals across the species’ northern range. We used genotypes at 11 microsatellite loci to reveal the genetic structure of, and directional gene flow among, populations to predict the likely future spread of the pathogen in the northwest and to estimate effective population size ( Ne). We also pyrosequenced the DRB1-like exon 2 of the class II major histocompatibility complex (MHC) in 160 individuals to explore immunogenetic selection by WNS. We identified three major neutral genetic clusters: Eastern, Montane Cordillera (and adjacent sampling areas), and Haida Gwaii, with admixture at intermediate areas and significant substructure west of the prairies. Estimates of Ne were unexpectedly low (289–16 000). Haida Gwaii may provide temporary refuge from WNS, but the western mountain ranges are not barriers to its dispersal in M. lucifugus and are unlikely to slow its spread. Our major histocompatibility complex (MHC) data suggest potential selection by WNS on the MHC, but gene duplication limited the immunogenetic analyses.

Novel Alphacoronaviruses and Paramyxoviruses Cocirculate with Type 1 and Severe Acute Respiratory System (SARS)-Related Betacoronaviruses in Synanthropic Bats of Luxembourg

Several infectious disease outbreaks with high mortality in humans have been attributed to viruses that are thought to have evolved from bat viruses. In this study from Luxembourg, the genetic diversity and epidemiology of paramyxoviruses and coronaviruses shed by the bat species Rhinolophus ferrumequinum and Myotis emarginatus were evaluated. Feces collection (n = 624) was performed longitudinally in a mixed-species colony in 2015 and 2016. In addition, feces (n = 254) were collected cross-sectionally from six Myotis emarginatus colonies in 2016. By use of degenerate primers in a nested format, overall prevalences of 1.1% (10/878) and 4.9% (43/878) were determined for paramyxoviruses and coronaviruses. Sequences of the partial RNA-dependent RNA polymerase and spike glycoprotein genes of coronaviruses, as well as sequences of the partial L gene of paramyxoviruses, were obtained. Novel paramyxovirus and Alphacoronavirus strains were identified in different Myotis emarginatus colonies, and severe acute respiratory syndrome (SARS)-related Betacoronavirus strains were shed by Rhinolophus ferrumequinum Logistic regression revealed that the level of Alphacoronavirus shedding was highest in July (odds ratio, 2.8; P < 0.01), probably due to periparturient stress. Phylogenetic analyses point to close virus-host coevolution, and the high genetic similarity of the study strains suggests that the Myotis emarginatus colonies in Luxembourg are socially connected. Most interestingly, we show that bats also host Betacoronavirus1 strains. The high similarity of the spike gene sequences of these viruses with mammalian Betacoronavirus 1 strains may be of concern. Both the SARS-related and Betacoronavirus 1 strains detected in bats in Luxembourg may cross the species barrier after a host adaptation process.IMPORTANCE Bats are a natural reservoir of a number of zoonotic pathogens. Several severe outbreaks in humans (e.g., a Nipah virus outbreak in Malaysia in 1998, and the almost global spread of severe acute respiratory syndrome in 2003) have been caused by bat-borne viruses that were transmitted to humans mostly after virus adaptation (e.g., in intermediate animal hosts). Despite the indigenousness of bat species that host viruses with suspected zoonotic potential and despite the zoonotic transmission of European bat 1 lyssavirus in Luxembourg, knowledge about the diversity and epidemiology of bat viruses remains limited in this country. Moreover, in contrast to other European countries, bat viruses are currently not included in the national surveillance activities of this land-locked country. We suggest that this gap in disease surveillance should be addressed, since we show here that synanthropic bats host viruses that may be able to cross the species barrier.

Range-wide genetic structure and demographic history in the bat ectoparasite Cimex adjunctus

Background

Evolutionary histories of parasite and host populations are intimately linked such that their spatial genetic structures may be correlated. While these processes have been relatively well studied in specialist parasites and their hosts, less is known about the ecological and evolutionary consequences of relationships between generalist ectoparasites and their hosts. The aim of this study was to investigate the genetic structure and demographic history of a bat ectoparasite, Cimex adjunctus, whose host affinity is weak but the biology of the potential hosts have been well studied. This ectoparasite has been hypothesized to rely on its hosts for dispersal due to its low inherent dispersal potential. Here we describe genetic diversity and demographic history in C. adjunctus through most of its range in North America. We investigated variation at the cytochrome c oxidase 1 mitochondrial gene and nine microsatellite markers, and tested the prediction that genetic diversity in C. adjunctus is spatially structured. We also tested the prediction that demographic history in C. adjunctus is characterized by range and demographic expansion as a consequence of post-Pleistocene climate warming.

Results

We found stronger spatial structuring of genetic diversity in C. adjunctus than has been quantified in two of its hosts, but contrast in amount of variation explained by host association with different genetic markers (i.e., nuclear vs mitochondrial DNA). Also, C. adjunctus’ history is not primarily characterized by demographic and range expansion, as is the case with two of its key hosts.

Conclusions

Our study shows different patterns of genetic structure and demographic history in C. adjunctus than have been detected in two of its key hosts. Our results suggest an effect of a loose parasite-host relationship and anti-parasitism strategies on genetic structure and post-Pleistocene recovery of population size.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0839-1) contains supplementary material, which is available to authorized users.

Population genetic structure of a common host predicts the spread of white-nose syndrome, an emerging infectious disease in bats

Landscape complexity influences patterns of animal dispersal, which in turn may affect both gene flow and the spread of pathogens. White-nose syndrome (WNS) is an introduced fungal disease that has spread rapidly throughout eastern North America, causing massive mortality in bat populations. We tested for a relationship between the population genetic structure of the most common host, the little brown myotis (Myotis lucifugus), and the geographic spread of WNS to date by evaluating logistic regression models of WNS risk among hibernating colonies in eastern North America. We hypothesized that risk of WNS to susceptible host colonies should increase with both geographic proximity and genetic similarity, reflecting historical connectivity, to infected colonies. Consistent with this hypothesis, inclusion of genetic distance between infected and susceptible colonies significantly improved models of disease spread, capturing heterogeneity in the spatial expansion of WNS despite low levels of genetic differentiation among eastern populations. Expanding our genetic analysis to the continental range of little brown myotis reveals strongly contrasting patterns of population structure between eastern and western North America. Genetic structure increases markedly moving westward into the northern Great Plains, beyond the current distribution of WNS. In western North America, genetic differentiation of geographically proximate populations often exceeds levels observed across the entire eastern region, suggesting infrequent and/or locally restricted dispersal, and thus relatively limited opportunities for pathogen introduction in western North America. Taken together, our analyses suggest a possibly slower future rate of spread of the WNS pathogen, at least as mediated by little brown myotis.

Range-Wide Genetic Analysis of Little Brown Bat (Myotis lucifugus) Populations: Estimating the Risk of Spread of White-Nose Syndrome

The little brown bat (Myotis lucifugus) is one of the most widespread bat species in North America and is experiencing severe population declines because of an emerging fungal disease, white-nose syndrome (WNS). To manage and conserve this species effectively it is important to understand patterns of gene flow and population connectivity to identify possible barriers to disease transmission. However, little is known about the population genetic structure of little brown bats, and to date, no studies have investigated population structure across their entire range. We examined mitochondrial DNA and nuclear microsatellites in 637 little brown bats (including all currently recognized subspecific lineages) from 29 locations across North America, to assess levels of genetic variation and population differentiation across the range of the species, including areas affected by WNS and those currently unaffected. We identified considerable spatial variation in patterns of female dispersal and significant genetic variation between populations in eastern versus western portions of the range. Overall levels of nuclear genetic differentiation were low, and there is no evidence for any major barriers to gene flow across their range. However, patterns of mtDNA differentiation are highly variable, with high ΦST values between most sample pairs (including between all western samples, between western and eastern samples, and between some eastern samples), while low mitochondrial differentiation was observed within two groups of samples found in central and eastern regions of North America. Furthermore, the Alaskan population was highly differentiated from all others, and western populations were characterized by isolation by distance while eastern populations were not. These data raise the possibility that the current patterns of spread of WNS observed in eastern North America may not apply to the entire range and that there may be broad-scale spatial variation in the risk of WNS transmission and occurrence if the disease continues to spread west.

Effect of oceanic straits on gene flow in the recently endangered little brown bat (Myotis lucifugus) in maritime Canada: implications for the spread of white-nose syndrome

White-nose syndrome is rapidly spreading in eastern North America, causing mass mortality of hibernating bats. We characterized levels of genetic diversity and population structure of the little brown bat (Myotis lucifugus (Le Conte, 1831)) in eastern Canada to infer the extent to which oceanic straits may be barriers to movement. To quantify metrics of gene flow and infer movement dynamics, we genotyped 679 M. lucifugus at nine nuclear microsatellites (nDNA) and sequenced a portion of the mitochondrial DNA (mtDNA). We found high levels of genetic diversity and little population structure, with ≈13-fold higher differentiation of mtDNA than nDNA markers, suggesting that structuring patterns largely result from female philopatry. Discriminant analysis of principle components suggested that the subtle underlying structure was not concordant with sampling site. Regional differentiation (FST, Dest, Mantel test residuals) is mostly consistent with genetic isolation by distance. However, samples from Newfoundland showed genetic differentiation over and above the effects of distance, lower levels of genetic diversity, and less genetic connectivity with other sampled regions. Despite this, oceanic straits in the Gulf of Saint Lawrence do not appear to create an impenetrable barrier to movement, therefore it may be possible for white-nose syndrome to spread to Newfoundland.

Population Genetic Structure Within and among Seasonal Site Types in the Little Brown Bat (Myotis lucifugus) and the Northern Long-Eared Bat (M. septentrionalis)

During late summer and early autumn, temperate bats migrate from their summering sites to swarming sites, where mating likely occurs. However, the extent to which individuals of a single summering site migrate to the same swarming site, and vice versa, is not known. We examined the migratory connectivity between summering and swarming sites in two temperate, North American, bat species, the little brown bat (Myotis lucifugus) and the northern long-eared bat (Myotis septentrionalis). Using mitochondrial and microsatellite DNA markers, we examined population structuring within and among summering and swarming sites. Both species exhibited moderate degrees of mitochondrial DNA differentiation (little brown bat: FST(SUMMER) = 0.093, FST(SWARMING) = 0.052; northern long-eared bat: FST(SUMMER) = 0.117, FST(SWARMING) = 0.043) and little microsatellite DNA differentiation among summering and among swarming sites[corrected]. Haplotype diversity was significantly higher at swarming sites than summering sites, supporting the idea that swarming sites are comprised of individuals from various summering sites. Further, pairwise analyses suggest that swarming sites are not necessarily comprised of only individuals from the most proximal summering colonies.

Carbon (δ13C) and Nitrogen (δ15N) Stable Isotope Signatures in Bat Fur Indicate Swarming Sites Have Catchment Areas for Bats from Different Summering Areas

Migratory patterns of bats are not well understood and traditional methods to study this, like capture-mark-recapture, may not provide enough detail unless there are many records. Stable isotope profiles of many animal species have been used to make inferences about migration. Each year Myotis lucifugus and M. septentrionalis migrate from summering roosts to swarming caves and mines in the fall, but the pattern of movement between them is not well understood. In this study, fur δ13C and δ15N values of 305 M. lucifugus and 200 M. septentrionalis were analyzed to make inferences about migration patterns between summering areas and swarming sites in Nova Scotia, Canada. We expected that there would be greater variability in δ13C and δ15N among individuals at swarming sites because it was believed that these sites are used by individuals originating from many summering areas. There was extensive overlap in the standard ellipse area, corrected for small sample sizes (SEAc), of bats at swarming sites and much less overlap in SEAc among groups sampled at summering areas. Meaningful inference could not be made on M. septentrionalis because their low variation in SEAc may have been the result of sampling only 3 summering areas. However, for M. lucifugus, swarming sites had larger SEAc than summering areas and predictive discriminant analysis assigned swarming bats to multiple summering areas, supporting the contention that swarming bats are mixed aggregations of bats from several summering areas. Together, these data support the contention that swarming sites have catchment areas for bats from multiple summering areas and it is likely that the catchment areas for swarming sites overlap. These data suggest that δ13C and δ15N profiling of bat fur offer some potential to make inferences about regional migration in bats.