Tracing Eastern Wolf Origins From Whole-Genome Data in Context of Extensive Hybridization.

Southeastern Canada is inhabited by an amalgam of hybridizing wolf-like canids, raising fundamental questions regarding their taxonomy, origins, and timing of hybridization events. Eastern wolves (Canis lycaon), specifically, have been the subject of significant controversy, being viewed as either a distinct taxonomic entity of conservation concern or a recent hybrid of coyotes (C. latrans) and grey wolves (C. lupus). Mitochondrial DNA analyses show some evidence of eastern wolves being North American evolved canids. In contrast, nuclear genome studies indicate eastern wolves are best described as a hybrid entity, but with unclear timing of hybridization events. To test hypotheses related to these competing findings we sequenced whole genomes of 25 individuals, representative of extant Canadian wolf-like canid types of known origin and levels of contemporary hybridization. Here we present data describing eastern wolves as a distinct taxonomic entity that evolved separately from grey wolves for the past ∼67,000 years with an admixture event with coyotes ∼37,000 years ago. We show that Great Lakes wolves originated as a product of admixture between grey wolves and eastern wolves after the last glaciation (∼8,000 years ago) while eastern coyotes originated as a product of admixture between "western" coyotes and eastern wolves during the last century. Eastern wolf nuclear genomes appear shaped by historical and contemporary gene flow with grey wolves and coyotes, yet evolutionary uniqueness remains among eastern wolves currently inhabiting a restricted range in southeastern Canada.

Transcriptome analysis of smut fungi reveals widespread intergenic transcription and conserved antisense transcript expression.

BACKGROUND: Biotrophic fungal plant pathogens cause billions of dollars in losses to North American crops annually. The model for functional investigation of these fungi is Ustilago maydis. Its 20.5 Mb annotated genome sequence has been an excellent resource for investigating biotrophic plant pathogenesis. Expressed-sequence tag libraries and microarray hybridizations have provided insight regarding the type of transcripts produced by U. maydis but these analyses were not comprehensive and there were insufficient data for transcriptome comparison to other smut fungi. To improve transcriptome annotation and enable comparative analyses, comprehensive strand-specific RNA-seq was performed on cell-types of three related smut species: U. maydis (common smut of corn), Ustilago hordei (covered smut of barley), and Sporisorium reilianum (head smut of corn). RESULTS: In total, >1 billion paired-end sequence reads were obtained from haploid cell, dikaryon and teliospore RNA of U. maydis, haploid cell RNA of U. hordei, and haploid and dikaryon cell RNA of S. reilianum. The sequences were assembled into transfrags using Trinity, and updated gene models were created using PASA and categorized with Cufflinks Cuffcompare. Representative genes that were predicted for the first time with these RNA-seq analyses and genes with novel annotation features were independently assessed by reverse transcriptase PCR. The analyses indicate hundreds more predicted proteins, relative to the previous genome annotation, could be produced by U. maydis from altered transcript forms, and that the number of non-coding RNAs produced, including transcribed intergenic sequences and natural antisense transcripts, approximately equals the number of mRNAs. This high representation of non-coding RNAs appears to be a conserved feature of the smut fungi regardless of whether they have RNA interference machinery. Approximately 50% of the identified NATs were conserved among the smut fungi. CONCLUSIONS: Overall, these analyses revealed: 1) smut genomes encode a number of transcriptional units that is twice the number of annotated protein-coding genes, 2) a small number of intergenic transcripts may encode proteins with characteristics of fungal effectors, 3) the vast majority of intergenic and antisense transcripts do not contain ORFs, 4) a large proportion of the identified antisense transcripts were detected at orthologous loci among the smut fungi, and 5) there is an enrichment of functional categories among orthologous loci that suggests antisense RNAs could have a genome-wide, non-RNAi-mediated, influence on gene expression in smut fungi.

A persistently infecting coronavirus in hibernating Myotis lucifugus, the North American little brown bat.

Bats are important reservoir hosts for emerging viruses, including coronaviruses that cause diseases in people. Although there have been several studies on the pathogenesis of coronaviruses in humans and surrogate animals, there is little information on the interactions of these viruses with their natural bat hosts. We detected a coronavirus in the intestines of 53/174 hibernating little brown bats (Myotis lucifugus), as well as in the lungs of some of these individuals. Interestingly, the presence of the virus was not accompanied by overt inflammation. Viral RNA amplified from little brown bats in this study appeared to be from two distinct clades. The sequences in clade 1 were very similar to the archived sequence derived from little brown bats and the sequences from clade 2 were more closely related to the archived sequence from big brown bats. This suggests that two closely related coronaviruses may circulate in little brown bats. Sequence variation among coronavirus detected from individual bats suggested that infection occurred prior to hibernation, and that the virus persisted for up to 4 months of hibernation in the laboratory. Based on the sequence of its genome, the coronavirus was placed in the Alphacoronavirus genus, along with some human coronaviruses, bat viruses and the porcine epidemic diarrhoea virus. The detection and identification of an apparently persistent coronavirus in a local bat species creates opportunities to understand the dynamics of coronavirus circulation in bat populations.

Tissue-specific transcriptome characterization for developing tadpoles of the northern leopard frog (Lithobates pipiens).

A potential cause of amphibian population declines are the impacts of environmental degradation on tadpole development. We conducted RNA sequencing on developing northern leopard frog tadpoles and through de novo transcriptome assembly we annotated a large number of open reading frames comparable in number and extent to genes identified in Xenopus. Using our transcriptome, we found transcript level changes between early (Gosner 26-31) and late (Gosner 36-41) stage tadpoles were the greatest in the tail, which is reabsorbed throughout development. There was an up-regulation of immunity genes in both the head and tail of the late tadpoles and a down-regulation of genes associated with the energy pathways of the mitochondria and the production of myosin. Overall, transcript level changes across development were consistent with studies on Xenopus and our findings highlight the broader utility of using RNA-seq to identify genes differentially expressed throughout development and in response to environmental pressures.

Ustilago maydis natural antisense transcript expression alters mRNA stability and pathogenesis.

Ustilago maydis infection of Zea mays leads to the production of thick-walled diploid teliospores that are the dispersal agent for this pathogen. Transcriptome analyses of this model biotrophic basidiomycete fungus identified natural antisense transcripts (NATs) complementary to 247 open reading frames. The U. maydis NAT cDNAs were fully sequenced and annotated. Strand-specific RT-PCR screens confirmed expression and identified NATs preferentially expressed in the teliospore. Targeted screens revealed four U. maydis NATs that are conserved in a related fungus. Expression of NATs in haploid cells, where they are not naturally occurring, resulted in increased steady-state levels of some complementary mRNAs. The expression of one NAT, as-um02151, in haploid cells resulted in a twofold increase in complementary mRNA levels, the formation of sense-antisense double-stranded RNAs, and unchanged Um02151 protein levels. This led to a model for NAT function in the maintenance and expression of stored teliospore mRNAs. In testing this model by deletion of the regulatory region, it was determined that alteration in NAT expression resulted in decreased pathogenesis in both cob and seedling infections. This annotation and functional analysis supports multiple roles for U. maydis NATs in controlling gene expression and influencing pathogenesis.

The impact of elevated sulfur and nitrogen levels on cadmium tolerance in Euglena species.

Heavy metal (HM) pollution threatens human and ecosystem health. Current methods for remediating water contaminated with HMs are expensive and have limited effect. Therefore, bioremediation is being investigated as an environmentally and economically viable alternative. Freshwater protists Euglena gracilis and Euglena mutabilis were investigated for their tolerance to cadmium (Cd). A greater increase in cell numbers under Cd stress was noted for E. mutabilis but only E. gracilis showed an increase in Cd tolerance following pre-treatment with elevated concentrations of S or N. To gain insight regarding the nature of the increased tolerance RNA-sequencing was carried out on E. gracilis. This revealed transcript level changes among pretreated cells, and additional differences among cells exposed to CdCl2. Gene ontology (GO) enrichment analysis reflected changes in S and N metabolism, transmembrane transport, stress response, and physiological processes related to metal binding. Identifying these changes enhances our understanding of how these organisms adapt to HM polluted environments and allows us to target development of future pre-treatments to enhance the use of E. gracilis in bioremediation relating to heavy metals.

Natural antisense transcripts in fungi.

Fungi are models for investigating many eukaryotic molecular processes. The identification of natural antisense transcripts (NATs) in fungi led to the discovery of mechanisms for controlling gene expression through transcriptional interference, chromatin remodelling and dsRNA formation. An overview of these mechanisms and the description of specific NAT functions is provided to give context for a broader discussion of fungal NATs. Transcriptome analyses have revealed a large number of NATs in a divergent group of fungi. The timing of NAT expression suggests roles in core life functions, such as responding to the environment and sexual reproduction. The transcriptome studies also uncover a large number of NATs whose functions remain elusive. These could provide novel control of gene expression, targeted responses to stimuli, or other functions. The goal of this review is provide background for this expanding field of research while highlighting opportunities for future discoveries.

Investigating the Ustilago maydis/Zea mays pathosystem: transcriptional responses and novel functional aspects of a fungal calcineurin regulatory B subunit.

The sustainable control of basidiomycete biotrophic plant pathogenesis requires an understanding of host responses to infection, as well as the identification and functional analysis of fungal genes involved in disease development. The creation and analysis of a suppressive subtractive hybridization (SSH) cDNA library from Ustilago maydis-infected Zea mays seedlings enabled the identification of fungal and plant genes expressed during disease development, and uncovered new insights into the interactions of this model system. Candidate U. maydis pathogenesis genes were identified by using the current SSH cDNA library analysis, and by knowledge generated from previous cDNA microarray and comparative genomic analyses. These identifications were supported by the independent determination of transcript level changes in different cell-types and during pathogenic development. The basidiomycete specific um01632, the highly in planta expressed um03046 (zig1), and the calcineurin regulatory B subunit (um10226, cnb1), were chosen for deletion experiments. um01632 and zig1 mutants showed no difference in morphology and did not have a statistically significant impact on pathogenesis. cnb1 mutants had a distinct cell division phenotype and reduced virulence in seedling assays. Infections with reciprocal wild-type×Δcnb1 haploid strain crosses revealed that the wild-type allele was unable to fully compensate for the lack of a second cnb1 allele. This haploinsufficiency was undetected in other fungal cnb1 mutational analyses. The reported data improves U. maydis genome annotation and expands on the current understanding of pathogenesis genes in this model basidiomycete.

Gene discovery in EST sequences from the wheat leaf rust fungus Puccinia triticina sexual spores, asexual spores and haustoria, compared to other rust and corn smut fungi.

BACKGROUND: Rust fungi are biotrophic basidiomycete plant pathogens that cause major diseases on plants and trees world-wide, affecting agriculture and forestry. Their biotrophic nature precludes many established molecular genetic manipulations and lines of research. The generation of genomic resources for these microbes is leading to novel insights into biology such as interactions with the hosts and guiding directions for breakthrough research in plant pathology. RESULTS: To support gene discovery and gene model verification in the genome of the wheat leaf rust fungus, Puccinia triticina (Pt), we have generated Expressed Sequence Tags (ESTs) by sampling several life cycle stages. We focused on several spore stages and isolated haustorial structures from infected wheat, generating 17,684 ESTs. We produced sequences from both the sexual (pycniospores, aeciospores and teliospores) and asexual (germinated urediniospores) stages of the life cycle. From pycniospores and aeciospores, produced by infecting the alternate host, meadow rue (Thalictrum speciosissimum), 4,869 and 1,292 reads were generated, respectively. We generated 3,703 ESTs from teliospores produced on the senescent primary wheat host. Finally, we generated 6,817 reads from haustoria isolated from infected wheat as well as 1,003 sequences from germinated urediniospores. Along with 25,558 previously generated ESTs, we compiled a database of 13,328 non-redundant sequences (4,506 singlets and 8,822 contigs). Fungal genes were predicted using the EST version of the self-training GeneMarkS algorithm. To refine the EST database, we compared EST sequences by BLASTN to a set of 454 pyrosequencing-generated contigs and Sanger BAC-end sequences derived both from the Pt genome, and to ESTs and genome reads from wheat. A collection of 6,308 fungal genes was identified and compared to sequences of the cereal rusts, Puccinia graminis f. sp. tritici (Pgt) and stripe rust, P. striiformis f. sp. tritici (Pst), and poplar leaf rust Melampsora species, and the corn smut fungus, Ustilago maydis (Um). While extensive homologies were found, many genes appeared novel and species-specific; over 40% of genes did not match any known sequence in existing databases. Focusing on spore stages, direct comparison to Um identified potential functional homologs, possibly allowing heterologous functional analysis in that model fungus. Many potentially secreted protein genes were identified by similarity searches against genes and proteins of Pgt and Melampsora spp., revealing apparent orthologs. CONCLUSIONS: The current set of Pt unigenes contributes to gene discovery in this major cereal pathogen and will be invaluable for gene model verification in the genome sequence.

Ustilago maydis transcript features identified through full-length cDNA analysis.

Ustilago maydis is the model for investigating basidiomycete biotrophic plant pathogens. To further the annotation of its genome, 12,943 full-length cDNA sequences were used to construct databases for the promoter and untranslated regions of U. maydis genes. A subset of clones was sequenced to determine full cDNA sequences. These and the original ESTs were assembled into contigs representing 3,058, or 45%, of the predicted U. maydis genes. The new sequencing allowed the confirmation of 2,842 gene models, 690 of which contain an intron. The use of full-length cDNA clone sequences ensured that untranslated regions were physically linked to the open reading frames (ORFs), not merely aligned upstream of the start of transcription. Identified sequence features include: (1) over 500 potential short upstream ORFs, (2) 95 gene models that require further annotation, (3) one new potential ORF, (4) varying GC content in different gene regions, (5) a WebLogo motif for the start of translation, (6) the correlation of UTR length with transcript representation in cDNA libraries and with gene function categories, (7) a relationship between natural antisense transcripts and UTR length that differs from that of Saccharomyces cerevisiae, (8) a potential relationship between DNA replication and the control of transcription, and (9) new insights regarding mechanisms for the control of transcription and mRNA maturation in U. maydis.

Genetic structure of immunologically associated candidate genes suggests arctic rabies variants exert differential selection in arctic fox populations.

Patterns of local adaptation can emerge in response to the selective pressures diseases exert on host populations as reflected in increased frequencies of respective, advantageous genotypes. Elucidating patterns of local adaptation enhance our understanding of mechanisms of disease spread and the capacity for species to adapt in context of rapidly changing environments such as the Arctic. Arctic rabies is a lethal disease that largely persists in northern climates and overlaps with the distribution of its natural host, arctic fox. Arctic fox populations display little neutral genetic structure across their North American range, whereas phylogenetically unique arctic rabies variants are restricted in their geographic distributions. It remains unknown if arctic rabies variants impose differential selection upon host populations, nor what role different rabies variants play in the maintenance and spread of this disease. Using a targeted, genotyping-by-sequencing assay, we assessed correlations of arctic fox immunogenetic variation with arctic rabies variants to gain further insight into the epidemiology of this disease. Corroborating past research, we found no neutral genetic structure between sampled regions, but did find moderate immunogenetic structuring between foxes predominated by different arctic rabies variants. FST outliers associated with host immunogenetic structure included SNPs within interleukin and Toll-like receptor coding regions (IL12B, IL5, TLR3 and NFKB1); genes known to mediate host responses to rabies. While these data do not necessarily reflect causation, nor a direct link to arctic rabies, the contrasting genetic structure of immunologically associated candidate genes with neutral loci is suggestive of differential selection and patterns of local adaptation in this system. These data are somewhat unexpected given the long-lived nature and dispersal capacities of arctic fox; traits expected to undermine local adaptation. Overall, these data contribute to our understanding of the co-evolutionary relationships between arctic rabies and their primary host and provide data relevant to the management of this disease.

The role of a mechanistic host in maintaining arctic rabies variant distributions: Assessment of functional genetic diversity in Alaskan red fox (Vulpes vulpes).

Populations are exposed to different types and strains of pathogens across heterogeneous landscapes, where local interactions between host and pathogen may present reciprocal selective forces leading to correlated patterns of spatial genetic structure. Understanding these coevolutionary patterns provides insight into mechanisms of disease spread and maintenance. Arctic rabies (AR) is a lethal disease with viral variants that occupy distinct geographic distributions across North America and Europe. Red fox (Vulpes vulpes) are a highly susceptible AR host, whose range overlaps both geographically distinct AR strains and regions where AR is absent. It is unclear if genetic structure exists among red fox populations relative to the presence/absence of AR or the spatial distribution of AR variants. Acquiring these data may enhance our understanding of the role of red fox in AR maintenance/spread and inform disease control strategies. Using a genotyping-by-sequencing assay targeting 116 genomic regions of immunogenetic relevance, we screened for sequence variation among red fox populations from Alaska and an outgroup from Ontario, including areas with different AR variants, and regions where the disease was absent. Presumed neutral SNP data from the assay found negligible levels of neutral genetic structure among Alaskan populations. The immunogenetically-associated data identified 30 outlier SNPs supporting weak to moderate genetic structure between regions with and without AR in Alaska. The outliers included SNPs with the potential to cause missense mutations within several toll-like receptor genes that have been associated with AR outcome. In contrast, there was a lack of genetic structure between regions with different AR variants. Combined, we interpret these data to suggest red fox populations respond differently to the presence of AR, but not AR variants. This research increases our understanding of AR dynamics in the Arctic, where host/disease patterns are undergoing flux in a rapidly changing Arctic landscape, including the continued northward expansion of red fox into regions previously predominated by the arctic fox (Vulpes lagopus).

Transcriptional host-pathogen responses of Pseudogymnoascus destructans and three species of bats with white-nose syndrome.

Understanding how context (e.g., host species, environmental conditions) drives disease susceptibility is an essential goal of disease ecology. We hypothesized that in bat white-nose syndrome (WNS), species-specific host-pathogen interactions may partly explain varying disease outcomes among host species. We characterized bat and pathogen transcriptomes in paired samples of lesion-positive and lesion-negative wing tissue from bats infected with Pseudogymnoascus destructans in three parallel experiments. The first two experiments analyzed samples collected from the susceptible Nearctic Myotis lucifugus and the less-susceptible Nearctic Eptesicus fuscus, following experimental infection and hibernation in captivity under controlled conditions. The third experiment applied the same analyses to paired samples from infected, free-ranging Myotis myotis, a less susceptible, Palearctic species, following natural infection and hibernation (n = 8 sample pairs/species). Gene expression by P. destructans was similar among the three host species despite varying environmental conditions among the three experiments and was similar within each host species between saprophytic contexts (superficial growth on wings) and pathogenic contexts (growth in lesions on the same wings). In contrast, we observed qualitative variation in host response: M. lucifugus and M. myotis exhibited systemic responses to infection, while E. fuscus up-regulated a remarkably localized response. Our results suggest potential phylogenetic determinants of response to WNS and can inform further studies of context-dependent host-pathogen interactions.

Bioinformatic identification of Ustilago maydis meiosis genes.

In the corn smut pathogen, Ustilago maydis, meiosis and teliospore germination are temporally linked. We review teliospore dormancy and germination in U. maydis and present an overview of meiosis in basidiomycetes. The relevant available expressed sequence tag data is discussed, the databases used in reciprocal best hit blastp analysis are presented and potential U. maydis meiosis genes are identified. The implications of identifying these genes are discussed and hypotheses are presented regarding the control of meiosis in U. maydis.

Development of a genotype-by-sequencing immunogenetic assay as exemplified by screening for variation in red fox with and without endemic rabies exposure.

Pathogens are recognized as major drivers of local adaptation in wildlife systems. By determining which gene variants are favored in local interactions among populations with and without disease, spatially explicit adaptive responses to pathogens can be elucidated. Much of our current understanding of host responses to disease comes from a small number of genes associated with an immune response. High-throughput sequencing (HTS) technologies, such as genotype-by-sequencing (GBS), facilitate expanded explorations of genomic variation among populations. Hybridization-based GBS techniques can be leveraged in systems not well characterized for specific variants associated with disease outcome to "capture" specific genes and regulatory regions known to influence expression and disease outcome. We developed a multiplexed, sequence capture assay for red foxes to simultaneously assess ~300-kbp of genomic sequence from 116 adaptive, intrinsic, and innate immunity genes of predicted adaptive significance and their putative upstream regulatory regions along with 23 neutral microsatellite regions to control for demographic effects. The assay was applied to 45 fox DNA samples from Alaska, where three arctic rabies strains are geographically restricted and endemic to coastal tundra regions, yet absent from the boreal interior. The assay provided 61.5% on-target enrichment with relatively even sequence coverage across all targeted loci and samples (mean = 50×), which allowed us to elucidate genetic variation across introns, exons, and potential regulatory regions (4,819 SNPs). Challenges remained in accurately describing microsatellite variation using this technique; however, longer-read HTS technologies should overcome these issues. We used these data to conduct preliminary analyses and detected genetic structure in a subset of red fox immune-related genes between regions with and without endemic arctic rabies. This assay provides a template to assess immunogenetic variation in wildlife disease systems.

Growth medium and incubation temperature alter the Pseudogymnoascus destructans transcriptome: implications in identifying virulence factors.

Pseudogymnoascus destructans is the causal agent of bat white-nose syndrome (WNS), which is devastating some North American bat populations. Previous transcriptome studies provided insight regarding the molecular mechanisms involved in WNS; however, it is unclear how different environmental parameters could influence pathogenicity. This information could be useful in developing management strategies to mitigate the negative impacts of P. destructans on bats. We cultured three P. destructans isolates from Atlantic Canada on two growth media (potato dextrose agar and Sabouraud dextrose agar) that differ in their nitrogen source, and at two separate incubation temperatures (4 C and 15 C) that approximate the temperature range of bat hibernacula during the winter and a temperature within its optimal mycelial growth range. We conducted RNA sequencing to determine transcript levels in each sample and performed differential gene expression (DGE) analyses to test the influence of growth medium and incubation temperature on gene expression. We also compared our in vitro results with previous RNA-sequencing data sets generated from P. destructans growing on the wings of a susceptible host, Myotis lucifugus. Our findings point to a critical role for substrate and incubation temperature in influencing the P. destructans transcriptome. DGE analyses suggested that growth medium plays a larger role than temperature in determining P. destructans gene expression and that although the psychrophilic fungus responds to different nitrogen sources, it may have evolved for continued growth at a broad range of low temperatures. Further, our data suggest that down-regulation of the RNA-interference pathway and increased fatty acid metabolism are involved in the P. destructans-bat interaction. Finally, we speculate that to reduce the activation of host defense responses, P. destructans minimizes changes in the expression of genes encoding secreted proteins during bat colonization.

Environmentally persistent pathogens present unique challenges for studies of host-pathogen interactions: Reply to Field (2018).

Linked article: https://doi.org/10.1002/ece3.4034.

White-nose syndrome is associated with increased replication of a naturally persisting coronaviruses in bats.

Spillover of viruses from bats to other animals may be associated with increased contact between them, as well as increased shedding of viruses by bats. Here, we tested the prediction that little brown bats (Myotis lucifugus) co-infected with the M. lucifugus coronavirus (Myl-CoV) and with Pseudogymnoascus destructans (Pd), the fungus that causes bat white-nose syndrome (WNS), exhibit different disease severity, viral shedding and molecular responses than bats infected with only Myl-CoV or only P. destructans. We took advantage of the natural persistence of Myl-CoV in bats that were experimentally inoculated with P. destructans in a previous study. Here, we show that the intestines of virus-infected bats that were also infected with fungus contained on average 60-fold more viral RNA than bats with virus alone. Increased viral RNA in the intestines correlated with the severity of fungus-related pathology. Additionally, the intestines of bats infected with fungus exhibited different expression of mitogen-activated protein kinase pathway and cytokine related transcripts, irrespective of viral presence. Levels of coronavirus antibodies were also higher in fungal-infected bats. Our results suggest that the systemic effects of WNS may down-regulate anti-viral responses in bats persistently infected with M. lucifugus coronavirus and increase the potential of virus shedding.

Profiling the immunome of little brown myotis provides a yardstick for measuring the genetic response to white-nose syndrome.

White-nose syndrome (WNS) has devastated populations of hibernating bats in eastern North America, leading to emergency conservation listings for several species including the previously ubiquitous little brown myotis (Myotis lucifugus). However, some bat populations near the epicenter of the WNS panzootic appear to be stabilizing after initial precipitous declines, which could reflect a selective immunogenetic sweep. To investigate the hypothesis that WNS exerts significant selection on the immunome of affected bat populations, we developed a novel, high-throughput sequence capture assay targeting 138 adaptive, intrinsic, and innate immunity genes of putative adaptive significance, as well as their respective regulatory regions (~370 kbp of genomic sequence/individual). We used the assay to explore baseline immunogenetic variation in M. lucifugus and to investigate whether particular immune genes/variants are associated with WNS susceptibility. We also used our assay to detect 1,038 putatively neutral single nucleotide polymorphisms and characterize contemporary population structure, providing context for the identification of local immunogenetic adaptation. Sequence capture provided a cost-effective, "all-in-one" assay to test for neutral genetic and immunogenetic structure and revealed fine-scale, baseline immunogenetic differentiation between sampling sites <600 km apart. We identified functional immunogenetic variants in M. lucifugus associated with WNS susceptibility. This study lays the foundations for future investigations of rangewide immunogenetic adaptation to WNS in M. lucifugus and provides a blueprint for studies of evolutionary rescue in other host-pathogen systems.

The other white-nose syndrome transcriptome: Tolerant and susceptible hosts respond differently to the pathogen Pseudogymnoascus destructans.

Mitigation of emerging infectious diseases that threaten global biodiversity requires an understanding of critical host and pathogen responses to infection. For multihost pathogens where pathogen virulence or host susceptibility is variable, host-pathogen interactions in tolerant species may identify potential avenues for adaptive evolution in recently exposed, susceptible hosts. For example, the fungus Pseudogymnoascus destructans causes white-nose syndrome (WNS) in hibernating bats and is responsible for catastrophic declines in some species in North America, where it was recently introduced. Bats in Europe and Asia, where the pathogen is endemic, are only mildly affected. Different environmental conditions among Nearctic and Palearctic hibernacula have been proposed as an explanation for variable disease outcomes, but this hypothesis has not been experimentally tested. We report the first controlled, experimental investigation of response to P. destructans in a tolerant, European species of bat (the greater mouse-eared bat, Myotis myotis). We compared body condition, disease outcomes and gene expression in control (sham-exposed) and exposed M. myotis that hibernated under controlled environmental conditions following treatment. Tolerant M. myotis experienced extremely limited fungal growth and did not exhibit symptoms of WNS. However, we detected no differential expression of genes associated with immune response in exposed bats, indicating that immune response does not drive tolerance of P. destructans in late hibernation. Variable responses to P. destructans among bat species cannot be attributed solely to environmental or ecological factors. Instead, our results implicate coevolution with the pathogen, and highlight the dynamic nature of the "white-nose syndrome transcriptome." Interspecific variation in response to exposure by the host (and possibly pathogen) emphasizes the importance of context in studies of the bat-WNS system, and robust characterization of genetic responses to exposure in various hosts and the pathogen should precede any attempts to use particular bat species as generalizable "model hosts."