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

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

What is winter? Modeling spatial variation in bat host traits and hibernation and their implications for overwintering energetics

White-nose syndrome (WNS) has decimated hibernating bat populations across eastern and central North America for over a decade. Disease severity is driven by the interaction between bat characteristics, the cold-loving fungal agent, and the hibernation environment. While we further improve hibernation energetics models, we have yet to examine how spatial heterogeneity in host traits is linked to survival in this disease system. Here, we develop predictive spatial models of body mass for the little brown myotis (Myotis lucifugus) and reassess previous definitions of the duration of hibernation of this species. Using data from published literature, public databases, local experts, and our own fieldwork, we fit a series of generalized linear models with hypothesized abiotic drivers to create distribution-wide predictions of prehibernation body fat and hibernation duration. Our results provide improved estimations of hibernation duration and identify a scaling relationship between body mass and body fat; this relationship allows for the first continuous estimates of prehibernation body mass and fat across the species' distribution. We used these results to inform a hibernation energetic model to create spatially varying fat use estimates for M. lucifugus. These results predict WNS mortality of M. lucifugus populations in western North America may be comparable to the substantial die-off observed in eastern and central populations.

Population genetics revealMyotis keenii(Keen’s myotis) andMyotis evotis(long-eared myotis) to be a single species

Recognizing delineations of gene flow among groups of animals can be challenging but is necessary for conservation and management. Of particular importance is the identification of species boundaries. Several physical and genetic traits have been used with mixed success to distinguish Myotis keenii (Merriam, 1895) (Keen’s myotis) and Myotis evotis (H. Allen, 1864) (long-eared myotis), but it is unclear whether species distinction is biologically warranted. We generated 12–14 microsatellite locus genotypes for 275 long-eared Myotis representing four species — M. keenii, M. evotis, Myotis septentrionalis (Trouessart, 1897) (northern myotis), and Myotis thysanodes Miller, 1897 (fringed myotis) — from across northwestern North America and 23 Myotis lucifugus (Le Conte, 1831) (little brown myotis) as the outgroup. Population genetic analyses revealed four well-defined groups (species): M. septentrionalis, M. thysanodes, M. lucifugus, and a single group comprising M. keenii and M. evotis. We document high rates of gene flow within M. evotis/M. keenii. Cytochrome b gene (mtDNA) sequencing failed to resolve morphologically identifiable species. We highlight the importance of geographically thorough investigation of genetic connectivity (nuclear markers) when assessing taxonomic status of closely related groups. We document a morphometric cline within M. evotis/M. keenii that may in part explain earlier analyses that led to the description of the smaller bodied M. keenii (type locality: Haida Gwaii). We conclude that M. keenii does not qualify as a genetic or biological species.

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.

Local coexistence and genetic isolation of three pollinator species on the same fig tree species

Molecular tools increasingly reveal cryptic lineages and species that were previously unnoticed by traditional taxonomy. The discovery of cryptic species in sympatry prompts the question of how they coexist in the apparent absence of ecological divergence. However, this assumes first that the molecular taxonomy used to identify cryptic lineages delimits species boundaries accurately. This issue is important, because many diversity studies rely heavily or solely on data from mitochondrial DNA sequences for species delimitation, and several factors may lead to poor identification of species boundaries. We used a multilocus population genetics approach to show that three mtDNA-defined cryptic lineages of the fig wasp Pleistodontes imperialis Saunders, which pollinate Port Jackson figs (Ficus rubiginosa) in north-eastern Australia, represent reproductively isolated species. These species coexist locally, with about 13% of figs (where mating occurs) containing wasps from two or three species. However, there was no evidence for gene flow between them. Confirmed cases of coexisting cryptic species provide excellent opportunities for future studies of the ecological and evolutionary forces shaping both species coexistence and fig/pollinator coevolution.

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.

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

Characterizing movement dynamics and spatial aspects of gene flow within a species permits inference on population structuring. As patterns of structuring are products of historical and current demographics and gene flow, assessment of structure through time can yield an understanding of evolutionary dynamics acting on populations that are necessary to inform management. Recent dramatic population declines in hibernating bats in eastern North America from white-nose syndrome have prompted the need for information on movement dynamics for multiple bat species. We characterized population genetic structure of the little brown bat, Myotis lucifugus, at swarming sites in southeastern Canada using 9 nuclear microsatellites and a 292-bp region of the mitochondrial genome. Analyses of F ST, ΦST, and Bayesian clustering (STRUCTURE) found weak levels of genetic structure among swarming sites for the nuclear and mitochondrial genome (Global F ST = 0.001, P < 0.05, Global ΦST = 0.045, P < 0.01, STRUCTURE K = 1) suggesting high contemporary gene flow. Hierarchical AMOVA also suggests little structuring at a regional (provincial) level. Metrics of nuclear genetic structure were not found to differ between males and females suggesting weak asymmetries in gene flow between the sexes. However, a greater degree of mitochondrial structuring does support male-biased dispersal long term. Demographic analyses were consistent with past population growth and suggest a population expansion occurred from approximately 1250 to 12,500 BP, following Pleistocene deglaciation in the region. Our study suggests high gene flow and thus a high degree of connectivity among bats that visit swarming sites whereby mainland areas of the region may be best considered as one large gene pool for management and conservation.

Post-Pleistocene range expansion of the recently imperiled eastern little brown bat (Myotis lucifugus lucifugus) from a single southern refugium

Myotis lucifugus, once among the most widespread and common bats in North America, has been forecast to be extirpated east of the Rockies in as few as 16 years by the spread of white-nose syndrome. Recent genetic research has demonstrated that this species is paraphyletic and part of a broader species complex; however, only one lineage (Myotis lucifugus lucifugus [M. l. lucifugus]) is present in eastern North America. I used molecular tools and niche modeling to validate this and investigate the role that historical biogeography has played in the phylogenetic and population genetic structure of this species to determine if the eastern subspecies represents an evolutionarily distinct population.

To establish the genetic structure within M. l. lucifugus, I densely sampled maternity colonies in Minnesota and sequenced 182 individuals for a portion of cytochrome b. Phylogenetic reconstruction and a haplotype network were used to infer the relationships among mitochondrial haplotypes. Population growth statistics were calculated to determine if there was evidence of significant expansion, and an environmental niche model (ENM) was constructed based on conditions during the last glacial maximum (LGM) to illustrate potential glacial refugia. All individuals derived from a single mitochondrial lineage. Genetic evidence points to population growth starting approximately 18 kya. ENM results show that there was likely a single large southern refugium extending across the southeastern United States and possibly several isolated refugia in western North America. Myotis lucifugus lucifugus likely maintained both a large range and a large population during the peaks of the glacial cycles, and its population appears to have expanded following the retreat of the Laurentide ice sheet. This imperiled lineage likely diverged in isolation from other members of the M. lucifugus/western long-eared Myotis during the Pleistocene.

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.