Linking avian malaria parasitemia estimates from quantitative PCR and microscopy reveals new infection patterns in Hawai'i

Plasmodium parasites infect thousands of species and provide an exceptional system for studying host-pathogen dynamics, especially for multi-host pathogens. However, understanding these interactions requires an accurate assay of infection. Assessing Plasmodium infections using microscopy on blood smears often misses infections with low parasitemias (the fractions of cells infected), and biases in malaria prevalence estimates will differ among hosts that differ in mean parasitemias. We examined Plasmodium relictum infection and parasitemia using both microscopy of blood smears and quantitative polymerase chain reaction (qPCR) on 299 samples from multiple bird species in Hawai'i and fit models to predict parasitemias from qPCR cycle threshold (Ct) values. We used these models to quantify the extent to which microscopy underestimated infection prevalence and to more accurately estimate infection patterns for each species for a large historical study done by microscopy. We found that most qPCR-positive wild-caught birds in Hawaii had low parasitemias (Ct scores ≥35), which were rarely detected by microscopy. The fraction of infections missed by microscopy differed substantially among eight species due to differences in species' parasitemia levels. Infection prevalence was likely 4-5-fold higher than previous microscopy estimates for three introduced species, including Zosterops japonicus, Hawaii's most abundant forest bird, which had low average parasitemias. In contrast, prevalence was likely only 1.5-2.3-fold higher than previous estimates for Himatione sanguinea and Chlorodrepanis virens, two native species with high average parasitemias. Our results indicate that relative patterns of infection among species differ substantially from those observed in previous microscopy studies, and that differences depend on variation in parasitemias among species. Although microscopy of blood smears is useful for estimating the frequency of different Plasmodium stages and host attributes, more sensitive quantitative methods, including qPCR, are needed to accurately estimate and compare infection prevalence among host species.

White-nose syndrome restructures bat skin microbiomes

IMPORTANCE

Inherent complexities in the composition of microbiomes can often preclude investigations of microbe-associated diseases. Instead of single organisms being associated with disease, community characteristics may be more relevant. Longitudinal microbiome studies of the same individual bats as pathogens arrive and infect a population are the ideal experiment but remain logistically challenging; therefore, investigations like our approach that are able to correlate invasive pathogens to alterations within a microbiome may be the next best alternative. The results of this study potentially suggest that microbiome-host interactions may determine the likelihood of infection. However, the contrasting relationship between Pd and the bacterial microbiomes of Myotis lucifugus and Perimyotis subflavus indicate that we are just beginning to understand how the bat microbiome interacts with a fungal invader such as Pd.

Microbiomes associated with avian malaria survival differ between susceptible Hawaiian honeycreepers and sympatric malaria-resistant introduced birds

Of the estimated 55 Hawaiian honeycreepers (subfamily Carduelinae) only 17 species remain, nine of which the International Union for Conservation of Nature considers endangered. Among the most pressing threats to honeycreeper survival is avian malaria, caused by the introduced blood parasite Plasmodium relictum, which is increasing in distribution in Hawai'i as a result of climate change. Preventing further honeycreeper decline will require innovative conservation strategies that confront malaria from multiple angles. Research on mammals has revealed strong connections between gut microbiome composition and malaria susceptibility, illuminating a potential novel approach to malaria control through the manipulation of gut microbiota. One honeycreeper species, Hawai'i 'amakihi (Chlorodrepanis virens), persists in areas of high malaria prevalence, indicating they have acquired some level of immunity. To investigate if avian host-specific microbes may be associated with malaria survival, we characterized cloacal microbiomes and malaria infection for 174 'amakihi and 172 malaria-resistant warbling white-eyes (Zosterops japonicus) from Hawai'i Island using 16S rRNA gene metabarcoding and quantitative polymerase chain reaction. Neither microbial alpha nor beta diversity covaried with infection, but 149 microbes showed positive associations with malaria survivors. Among these were Escherichia and Lactobacillus spp., which appear to mitigate malaria severity in mammalian hosts, revealing promising candidates for future probiotic research for augmenting malaria immunity in sensitive endangered species.

Host infection and disease-induced mortality modify species contributions to the environmental reservoir

Environmental pathogen reservoirs exist for many globally important diseases and can fuel epidemics, influence pathogen evolution, and increase the threat of host extinction. Species composition can be an important factor that shapes reservoir dynamics and ultimately determines the outcome of a disease outbreak. However, disease-induced mortality can change species communities, indicating that species responsible for environmental reservoir maintenance may change over time. Here we examine the reservoir dynamics of Pseudogymnoascus destructans, the fungal pathogen that causes white-nose syndrome in bats. We quantified changes in pathogen shedding, infection prevalence and intensity, host abundance, and the subsequent propagule pressure imposed by each species over time. We find that highly shedding species are important during pathogen invasion, but contribute less over time to environmental contamination as they also suffer the greatest declines. Less infected species remain more abundant, resulting in equivalent or higher propagule pressure. More broadly, we demonstrate that high infection intensity and subsequent mortality during disease progression can reduce the contributions of high-shedding species to long-term pathogen maintenance.

If You're Not Confused, You're Not Paying Attention: Ochrobactrum Is Not Brucella

Bacteria of the genus Brucella are facultative intracellular parasites that cause brucellosis, a severe animal and human disease. Recently, a group of taxonomists merged the brucellae with the primarily free-living, phylogenetically related Ochrobactrum spp. in the genus Brucella. This change, founded only on global genomic analysis and the fortuitous isolation of some opportunistic Ochrobactrum spp. from medically compromised patients, has been automatically included in culture collections and databases. We argue that clinical and environmental microbiologists should not accept this nomenclature, and we advise against its use because (i) it was presented without in-depth phylogenetic analyses and did not consider alternative taxonomic solutions; (ii) it was launched without the input of experts in brucellosis or Ochrobactrum; (iii) it applies a non-consensus genus concept that disregards taxonomically relevant differences in structure, physiology, population structure, core-pangenome assemblies, genome structure, genomic traits, clinical features, treatment, prevention, diagnosis, genus description rules, and, above all, pathogenicity; and (iv) placing these two bacterial groups in the same genus creates risks for veterinarians, medical doctors, clinical laboratories, health authorities, and legislators who deal with brucellosis, a disease that is particularly relevant in low- and middle-income countries. Based on all this information, we urge microbiologists, bacterial collections, genomic databases, journals, and public health boards to keep the Brucella and Ochrobactrum genera separate to avoid further bewilderment and harm.

Comparison of Brucella abortus population structure based on genotyping methods with different levels of resolution

Numerous genotyping techniques based on different principles and with different costs and levels of resolution are currently available for understanding the transmission dynamics of brucellosis worldwide. We aimed to compare the population structure of the genomes of 53 Brazilian Brucella abortus isolates using eight different genotyping methods: multiple-locus variable-number tandem-repeat analysis (MLVA8, MLVA11, MLVA16), multilocus sequence typing (MLST9, MLST21), core genome MLST (cgMLST) and two techniques based on single nucleotide polymorphism (SNP) detection (parSNP and NASP) from whole genomes. The strains were isolated from six different Brazilian states between 1977 and 2008 and had previously been analyzed using MLVA8, MLVA11, and MLVA16. Their whole genomes were sequenced, assembled, and subjected to MSLT9 MLST21, cgMLST, and SNP analyses. All the genotypes were compared by hierarchical grouping method based on the average distances between the correlation matrices of each technique. MLST9 and MLST21 had the lowest level of resolution, both revealing only four genotypes. MLVA8, MLVA11, and MLVA16 had progressively increasing levels of resolution as more loci were analyzed, identifying 6, 16, and 44 genotypes, respectively. cgMLST showed the highest level of resolution, identifying 45 genotypes, followed by the SNP-based methods, both of which had 44 genotypes. In the assessed population, MLVA was more discriminatory than MLST and was easier and cheaper to perform. SNP techniques and cgMLST provided the highest levels of resolution and the results from the two methods were in close agreement. In conclusion, the choice of genotyping technique can strongly affect one's ability to make meaningful epidemiological conclusions but is dependent on available resources: while the VNTR based techniques are more indicated to high prevalence scenarios, the WGS methods are the ones with the best discriminative power and therefore recommended for outbreaks investigation.

Sex-biased infections scale to population impacts for an emerging wildlife disease

Demographic factors are fundamental in shaping infectious disease dynamics. Aspects of populations that create structure, like age and sex, can affect patterns of transmission, infection intensity and population outcomes. However, studies rarely link these processes from individual to population-scale effects. Moreover, the mechanisms underlying demographic differences in disease are frequently unclear. Here, we explore sex-biased infections for a multi-host fungal disease of bats, white-nose syndrome, and link disease-associated mortality between sexes, the distortion of sex ratios and the potential mechanisms underlying sex differences in infection. We collected data on host traits, infection intensity and survival of five bat species at 42 sites across seven years. We found females were more infected than males for all five species. Females also had lower apparent survival over winter and accounted for a smaller proportion of populations over time. Notably, female-biased infections were evident by early hibernation and likely driven by sex-based differences in autumn mating behaviour. Male bats were more active during autumn which likely reduced replication of the cool-growing fungus. Higher disease impacts in female bats may have cascading effects on bat populations beyond the hibernation season by limiting recruitment and increasing the risk of Allee effects.

Trends in Bacterial Pathogens of Bats: Global Distribution and Knowledge Gaps

Bats have received considerable recent attention for infectious disease research because of their potential to host and transmit viruses, including Ebola, Hendra, Nipah, and multiple coronaviruses. These pathogens are occasionally transmitted from bats to wildlife, livestock, and to humans, directly or through other bridging (intermediate) hosts. Due to their public health relevance, zoonotic viruses are a primary focus of research attention. In contrast, other emerging pathogens of bats, such as bacteria, are vastly understudied despite their ubiquity and diversity. Here, we describe the currently known host ranges and geographic distributional patterns of potentially zoonotic bacterial genera in bats, using published presence-absence data of pathogen occurrence. We identify apparent gaps in our understanding of the distribution of these pathogens on a global scale. The most frequently detected bacterial genera in bats are Bartonella, Leptospira, and Mycoplasma. However, a wide variety of other potentially zoonotic bacterial genera are also occasionally found in bats, such as Anaplasma, Brucella, Borrelia, Coxiella, Ehrlichia, Francisella, Neorickettsia, and Rickettsia. The bat families Phyllostomidae, Vespertilionidae, and Pteropodidae are most frequently reported as hosts of bacterial pathogens; however, the presence of at least one bacterial genus was confirmed in all 15 bat families tested. On a spatial scale, molecular diagnostics of samples from 58 countries and four overseas departments and island states (French Guiana, Mayotte, New Caledonia, and Réunion Island) reported testing for at least one bacterial pathogen in bats. We also identified geographical areas that have been mostly neglected during bacterial pathogen research in bats, such as the Afrotropical region and Southern Asia. Current knowledge on the distribution of potentially zoonotic bacterial genera in bats is strongly biased by research effort towards certain taxonomic groups and geographic regions. Identifying these biases can guide future surveillance efforts, contributing to a better understanding of the ecoepidemiology of zoonotic pathogens in bats.

Shifting effects of host physiological condition following pathogen establishment

Understanding host persistence with emerging pathogens is essential for conserving populations. Hosts may initially survive pathogen invasions through pre-adaptive mechanisms. However, whether pre-adaptive traits are directionally selected to increase in frequency depends on the heritability and environmental dependence of the trait and the costs of trait maintenance. Body condition is likely an important pre-adaptive mechanism aiding in host survival, although can be seasonally variable in wildlife hosts. We used data collected over 7 years on bat body mass, infection and survival to determine the role of host body condition during the invasion and establishment of the emerging disease, white-nose syndrome. We found that when the pathogen first invaded, bats with higher body mass were more likely to survive, but this effect dissipated following the initial epizootic. We also found that heavier bats lost more weight overwinter, but fat loss depended on infection severity. Lastly, we found mixed support that bat mass increased in the population after pathogen arrival; high annual plasticity in individual bat masses may have reduced the potential for directional selection. Overall, our results suggest that some factors that contribute to host survival during pathogen invasion may diminish over time and are potentially replaced by other host adaptations.

First report and whole-genome sequencing of Pseudochrobactrum saccharolyticum in Latin America

The Brucellaceae family comprises microorganisms similar both phenotypically and genotypically, making it difficult to identify the etiological agent of these infections. This study reports the first isolation, identification, and characterization of Pseudochrobactrum saccharolyticum (strain 115) from Latin America. Strain 115 was isolated in 2007 from a bovine in Brazil and was initially classified as Brucella spp. by classical microbiological tests and bcsp31 PCR. The antimicrobial susceptibility of strain 115 was tested against drugs used to treat human brucellosis by minimal inhibitory concentration test. Subsequently, the whole genome of the strain was sequenced, assembled, and characterized. Phylogenetic trees built from 16S rRNA and recA gene sequences enabled the classification of strain 115 as Pseudochrobactrum spp. Phylogenomic analysis using Single Nucleotide Polymorphisms and Average Nucleotide Identity allowed the classification of the strain as P. saccharolyticum. Additionally, a Tetra Correlation Search identified one related genome from the same species, which was compared with strain 115 by analyzing genomic islands. This is the first identification and whole-genome sequence of P. saccharolyticum in Latin America and highlights a challenge in the diagnosis of bovine brucellosis, which could be solved by including the sequencing of 16S rRNA and recA genes in routine diagnostics.

Mechanisms underlying interaction frequencies and robustness in a novel seed dispersal network: lessons for restoration

As human-caused extinctions and invasions accumulate across the planet, understanding the processes governing ecological functions mediated by species interactions, and anticipating the effect of species loss on such functions become increasingly urgent. In seed dispersal networks, the mechanisms that influence interaction frequencies may also influence the capacity of a species to switch to alternative partners (rewiring), influencing network robustness. Studying seed dispersal interactions in novel ecosystems on O'ahu island, Hawai'i, we test whether the same mechanisms defining interaction frequencies can regulate rewiring and increase network robustness to simulated species extinctions. We found that spatial and temporal overlaps were the primary mechanisms underlying interaction frequencies, and the loss of the more connected species affected networks to a greater extent. Further, rewiring increased network robustness, and morphological matching and spatial and temporal overlaps between partners were more influential on network robustness than species abundances. We argue that to achieve self-sustaining ecosystems, restoration initiatives can consider optimal morphological matching and spatial and temporal overlaps between consumers and resources to maximize chances of native plant dispersal. Specifically, restoration initiatives may benefit from replacing invasive species with native species possessing characteristics that promote frequent interactions and increase the probability of rewiring (such as long fruiting periods, small seeds and broad distributions).

Spatial and temporal variation in New Hampshire bat diets

Insectivorous bats consume a diverse array of arthropod prey, with diets varying by bat species, sampling location, and season. North American bat diets remain incompletely described, which is concerning at a time when many bat and insect populations appear to be declining. Understanding the variability in foraging is thus an essential component for effective bat conservation. To comprehensively evaluate local foraging, we assessed the spatial and temporal variability in prey consumed by the little brown bat, Myotis lucifugus, in New Hampshire, USA. We collected bat guano samples from 20 sites over 2 years and analyzed sequence data for 899 of these samples using a molecular metabarcoding approach targeting the cytochrome oxidase I subunit (COI) gene. Some prey items were broadly shared across locations and sampling dates, with the most frequently detected arthropod orders broadly similar to previous morphological and molecular analyses; at least one representative sequence variant was assigned to Coleoptera in 92% of samples, with other frequently detected orders including Diptera (73%), Lepidoptera (65%), Trichoptera (38%), and Ephemeroptera (32%). More specifically, two turf and forest pests were routinely detected: white grubs in the genus Phyllophaga (50%), and the Asiatic Garden beetle, Maladera castanea (36%). Despite the prevalence of a few taxa shared among many samples and distinct seasonal peaks in consumption of specific arthropods, diet composition varied both temporally and spatially. However, species richness did not strongly vary indicating consumption of a broad diversity of taxa throughout the summer. These data characterize little brown bats as flexible foragers adept at consuming a broad array of locally available prey resources.

Phenotypic homogeneity of emetic Bacillus cereus isolates in China

Emetic Bacillus cereus strains produce a potent cereulide cytotoxin, which can cause acute and fatal cases of food poisoning. We isolated 18 emetic B. cereus strains from a food poisoning event, and from clinical and non-random food surveillance in China and phenotypic characteristics of haemolysis, starch hydrolysis, salicin fermentation, gelatin liquefaction, cytotoxicity, and susceptibility to antibiotics were assessed. All isolates were positive for haemolysis and gelatin liquefaction, and negative for starch hydrolysis and salicin fermentation. Their haemolytic potentials were intermediate to Bacillus anthracis and B. cereus ATCC 14579 (a non-emetic strain). All isolates were cytotoxic to CHO, Hep-2, and Vero cells, and were sensitive to ampicillin. The homogeneous phenotypes of emetic isolates from China are similar to the corresponding traits of European and Japanese isolates that have been characterized, suggesting highly similar phenotypes of emetic B. cereus worldwide.

Borreliosis Transmission from Ticks Associated with Desert Tortoise Burrows: Examples of Tick-Borne Relapsing Fever in the Mojave Desert

Ticks transmit pathogens and parasitize wildlife in turn causing zoonotic diseases in many ecosystems. Argasid ticks, such as Ornithodoros spp., harbor and transmit Borrelia spp., resulting in tick-borne relapsing fever (TBRF) in people. In the western United States, TBRF is typically associated with the bite of an infected Ornithodoros hermsi tick found in habitats at high elevations (>1500 ft). This report describes the first TBRF cases in people in the Mojave Desert (Clark County, NV). Individuals documented in these case studies were exposed to Ornithodoros ticks during excavation of soil burrows associated with Mojave Desert tortoises (Gopherus agassizii), with bacteria from one of the human case's blood sample genetically matching to Borrelia turicatae as determined by quantitative PCR and sequencing. Our findings should serve as a precaution to individuals working with tortoises or animal burrows, or those in contact with Ornithodoros ticks in this region.

St. Louis Encephalitis Virus in the Southwestern United States: A Phylogeographic Case for a Multi-Variant Introduction Event

Since the reemergence of St. Louis Encephalitis (SLE) Virus (SLEV) in the Southwest United States, identified during the 2015 outbreak in Arizona, SLEV has been seasonally detected within Culex spp. populations throughout the Southwest United States. Previous work revealed the 2015 outbreak was caused by an importation of SLEV genotype III, which had only been detected previously in Argentina. However, little is known about when the importation occurred or the transmission and genetic dynamics since its arrival into the Southwest. In this study, we sought to determine whether the annual detection of SLEV in the Southwest is due to enzootic cycling or new importations. To address this question, we analyzed 174 SLEV genomes (142 sequenced as part of this study) using Bayesian phylogenetic analyses to estimate the date of arrival into the American Southwest and characterize the underlying population structure of SLEV. Phylogenetic clustering showed that SLEV variants circulating in Maricopa and Riverside counties form two distinct populations with little evidence of inter-county transmission since the onset of the outbreak. Alternatively, it appears that in 2019, Yuma and Clark counties experienced annual importations of SLEV that originated in Riverside and Maricopa counties. Finally, the earliest representatives of SLEV genotype III in the Southwest form a polytomy that includes both California and Arizona samples. We propose that the initial outbreak most likely resulted from the importation of a population of SLEV genotype III variants, perhaps in multiple birds, possibly multiple species, migrating north in 2013, rather than a single variant introduced by one bird.

First Case of Brucellosis Caused by an Amphibian-type Brucella

We report the first case of brucellosis caused by an isolate whose genome is identical that of a frog isolate from Texas, demonstrating the zoonotic potential of amphibian-type Brucella. Importantly, with such atypical Brucella, correct diagnosis cannot be performed using routine serological tests or identification methods.

Emerging diversity and ongoing expansion of the genus Brucella

Remarkable genetic diversity and breadth of host species has been uncovered in the Brucella genus over the past decade, fundamentally changing our concept of what it means to be a Brucella. From ocean fishes and marine mammals, to pond dwelling amphibians, forest foxes, desert rodents, and cave-dwelling bats, Brucella have revealed a variety of previously unknown niches. Classical microbiological techniques have been able to help us classify many of these new strains but at times have limited our ability to see the true relationships among or within species. The closest relatives of Brucella are soil bacteria and the adaptations of Brucella spp. to live intracellularly suggest that the genus has evolved to live in vertebrate hosts. Several recently discovered species appear to have phenotypes that are intermediate between soil bacteria and core Brucella, suggesting that they may represent ancestral traits that were subsequently lost in the traditional species. Remarkably, the broad relationships among Brucella species using a variety of sequence and fragment-based approaches have been upheld when using comparative genomics with whole genomes. Nonetheless, genomes are required for fine-scale resolution of many of the relationships and for understanding the evolutionary history of the genus. We expect that the coming decades will reveal many more hosts and previously unknown diversity in a wide range of environments.

Mobility and infectiousness in the spatial spread of an emerging fungal pathogen

Emerging infectious diseases can have devastating effects on host communities, causing population collapse and species extinctions. The timing of novel pathogen arrival into naïve species communities can have consequential effects that shape the trajectory of epidemics through populations. Pathogen introductions are often presumed to occur when hosts are highly mobile. However, spread patterns can be influenced by a multitude of other factors including host body condition and infectiousness. White-nose syndrome (WNS) is a seasonal emerging infectious disease of bats, which is caused by the fungal pathogen Pseudogymnoascus destructans. Within-site transmission of P. destructans primarily occurs over winter; however, the influence of bat mobility and infectiousness on the seasonal timing of pathogen spread to new populations is unknown. We combined data on host population dynamics and pathogen transmission from 22 bat communities to investigate the timing of pathogen arrival and the consequences of varying pathogen arrival times on disease impacts. We found that midwinter arrival of the fungus predominated spread patterns, suggesting that bats were most likely to spread P. destructans when they are highly infectious, but have reduced mobility. In communities where P. destructans was detected in early winter, one species suffered higher fungal burdens and experienced more severe declines than at sites where the pathogen was detected later in the winter, suggesting that the timing of pathogen introduction had consequential effects for some bat communities. We also found evidence of source-sink population dynamics over winter, suggesting some movement among sites occurs during hibernation, even though bats at northern latitudes were thought to be fairly immobile during this period. Winter emergence behaviour symptomatic of white-nose syndrome may further exacerbate these winter bat movements to uninfected areas. Our results suggest that low infectiousness during host migration may have reduced the rate of expansion of this deadly pathogen, and that elevated infectiousness during winter plays a key role in seasonal transmission. Furthermore, our results highlight the importance of both accurate estimation of the timing of pathogen spread and the consequences of varying arrival times to prevent and mitigate the effects of infectious diseases.

Lord of the Diptera (and Moths and a Spider): Molecular Diet Analyses and Foraging Ecology of Indiana Bats in Illinois

Effective management of endangered or threatened wildlife requires an understanding of how foraging habitats are used by those populations. Molecular diet analysis of fecal samples offers a cost-effective and non-invasive method to investigate how diets of wild populations vary with respect to spatial and temporal factors. For the federally endangered Indiana bat (Myotis sodalis), documenting its preferred food sources can provide critical information to promote effective conservation of this federally endangered species. Using cytochrome oxidase I amplicon sequence data from Indiana bat guano samples collected at two roosting areas in Cypress Creek National Wildlife Refuge, we found that dipteran taxa (i.e., flies) associated with riparian habitats were the most frequently detected taxon and represented the majority of the sequence diversity among the arthropods sampled. A select few arthropods from other taxa—especially spiders—are also likely important to Indiana bat diets in this refuge. A supervised learning analysis of diet components suggest only a small fraction of the frequently detected taxa are important contributors to spatial and temporal variation. Overall, these data depict the Indiana bat as a generalist consumer whose diet includes some prey items associated with particular seasonal or spatial components, along with other taxa repeatedly consumed throughout the entire foraging season. These molecular diet analyses suggest that protecting foraging resources specifically associated with the riparian habitat of Cypress Creek National Wildlife Refuge is essential to promote effective Indiana bat conservation.

Herd immunity drives the epidemic fadeout of avian cholera in Arctic-nesting seabirds

Avian cholera, caused by the bacterium Pasteurella multocida, is a common and important infectious disease of wild birds in North America. Between 2005 and 2012, avian cholera caused annual mortality of widely varying magnitudes in Northern common eiders (Somateria mollissima borealis) breeding at the largest colony in the Canadian Arctic, Mitivik Island, Nunavut. Although herd immunity, in which a large proportion of the population acquires immunity to the disease, has been suggested to play a role in epidemic fadeout, immunological studies exploring this hypothesis have been missing. We investigated the role of three potential drivers of fadeout of avian cholera in eiders, including immunity, prevalence of infection, and colony size. Each potential driver was examined in relation to the annual real-time reproductive number (R_t) of P. multocida, previously calculated for eiders at Mitivik Island. Each year, colony size was estimated and eiders were closely monitored, and evaluated for infection and serological status. We demonstrate that acquired immunity approximated using antibody titers to P. multocida in both sexes was likely a key driver for the epidemic fadeout. This study exemplifies the importance of herd immunity in influencing the dynamics and fadeout of epidemics in a wildlife population.

Evolutionary history and current distribution of the West Mediterranean lineage of Brucella melitensis in Italy

Ovine and caprine brucellosis, caused by Brucella melitensis, is one of the world’s most widespread zoonoses and is a major cause of economic losses in domestic ruminant production. In Italy, the disease remains endemic in several southern provinces, despite an ongoing brucellosis eradication programme. In this study, we used whole-genome sequencing to detail the genetic diversity of circulating strains, and to examine the origins of the predominant sub-lineages of B. melitensis in Italy. We reconstructed a global phylogeny of B. melitensis, strengthened by 339 new whole-genome sequences, from Italian isolates collected from 2011 to 2018 as part of a national livestock surveillance programme. All Italian strains belonged to the West Mediterranean lineage, which further divided into two major clades that diverged roughly between the 5th and 7th centuries. We observed that Sicily serves as a brucellosis burden hotspot, giving rise to several distinct sub-lineages. More than 20 putative outbreak clusters of ovine and caprine brucellosis were identified, several of which persisted over the 8 year survey period despite an aggressive brucellosis eradication campaign. While the outbreaks in Central and Northern Italy were generally associated with introductions of single clones of B. melitensis and their subsequent dissemination within neighbouring territories, we observed weak geographical segregation of genotypes in the southern regions. Biovar determination, recommended in routine analysis of all Brucella strains by the World Organisation for Animal Health (OIE), could not discriminate among the four main global clades. This demonstrates a need for updating the guidelines used for monitoring B. melitensis transmission and spread, both at the national and international level, and to include whole-genome-based typing as the principal method for identification and tracing of brucellosis outbreaks.

A total crapshoot? Evaluating bioinformatic decisions in animal diet metabarcoding analyses

Metabarcoding studies provide a powerful approach to estimate the diversity and abundance of organisms in mixed communities in nature. While strategies exist for optimizing sample and sequence library preparation, best practices for bioinformatic processing of amplicon sequence data are lacking in animal diet studies. Here we evaluate how decisions made in core bioinformatic processes, including sequence filtering, database design, and classification, can influence animal metabarcoding results. We show that denoising methods have lower error rates compared to traditional clustering methods, although these differences are largely mitigated by removing low-abundance sequence variants. We also found that available reference datasets from GenBank and BOLD for the animal marker gene cytochrome oxidase I (COI) can be complementary, and we discuss methods to improve existing databases to include versioned releases. Taxonomic classification methods can dramatically affect results. For example, the commonly used Barcode of Life Database (BOLD) Classification API assigned fewer names to samples from order through species levels using both a mock community and bat guano samples compared to all other classifiers (vsearch-SINTAX and q2-feature-classifier's BLAST + LCA, VSEARCH + LCA, and Naive Bayes classifiers). The lack of consensus on bioinformatics best practices limits comparisons among studies and may introduce biases. Our work suggests that biological mock communities offer a useful standard to evaluate the myriad computational decisions impacting animal metabarcoding accuracy. Further, these comparisons highlight the need for continual evaluations as new tools are adopted to ensure that the inferences drawn reflect meaningful biology instead of digital artifacts.

Ambiguities in using telomere length for age determination in two North American bat species

The age of an animal, determined by time (chronological age) as well as genetic and environmental factors (biological age), influences the likelihood of mortality and reproduction and thus the animal’s contribution to population growth. For many long-lived species, such as bats, a lack of external and morphological indicators has made determining age a challenge, leading researchers to examine genetic markers of age for application to demographic studies. One widely studied biomarker of age is telomere length, which has been related both to chronological and biological age across taxa, but only recently has begun to be studied in bats. We assessed telomere length from the DNA of known-age and minimum known-age individuals of two bat species using a quantitative PCR assay. We determined that telomere length was quadratically related to chronological age in big brown bats (Eptesicus fuscus), although it had little predictive power for accurate age determination of unknown-age individuals. The relationship was different in little brown bats (Myotis lucifugus), where telomere length instead was correlated with biological age, apparently due to infection and wing damage associated with white-nose syndrome. Furthermore, we showed that wing biopsies currently are a better tissue source for studying telomere length in bats than guano and buccal swabs; the results from the latter group were more variable and potentially influenced by storage time. Refinement of collection and assessment methods for different non-lethally collected tissues will be important for longitudinal sampling to better understand telomere dynamics in these long-lived species. Although further work is needed to develop a biomarker capable of determining chronological age in bats, our results suggest that biological age, as reflected in telomere length, may be influenced by extrinsic stressors such as disease.

Persistence of Brucella abortus lineages revealed by genomic characterization and phylodynamic analysis

Brucellosis, caused by Brucella abortus, is a major disease of cattle and humans worldwide distributed. Eradication and control of the disease has been difficult in Central and South America, Central Asia, the Mediterranean and the Middle East. Epidemiological strategies combined with phylogenetic methods provide the high-resolution power needed to study relationships between surveillance data and pathogen population dynamics, using genetic diversity and spatiotemporal distributions. This information is crucial for prevention and control of disease spreading at a local and worldwide level. In Costa Rica (CR), the disease was first reported at the beginning of the 20th century and has not been controlled despite many efforts. We characterized 188 B. abortus isolates from CR recovered from cattle, humans and water buffalo, from 2003 to 2018, and whole genome sequencing (WGS) was performed in 95 of them. They were also assessed based on geographic origin, date of introduction, and phylogenetic associations in a worldwide and national context. Our results show circulation of five B. abortus lineages (I to V) in CR, phylogenetically related to isolates from the United States, United Kingdom, and South America. Lineage I was dominant and probably introduced at the end of the 19th century. Lineage II, represented by a single isolate from a water buffalo, clustered with a Colombian sample, and was likely introduced after 1845. Lineages III and IV were likely introduced during the early 2000s. Fourteen isolates from humans were found within the same lineage (lineage I) regardless of their geographic origin within the country. The main CR lineages, introduced more than 100 years ago, are widely spread throughout the country, in contrast to new introductions that seemed to be more geographically restricted. Following the brucellosis prevalence and the farming practices of several middle- and low-income countries, similar scenarios could be found in other regions worldwide.

Pathogen to commensal? Longitudinal within-host population dynamics, evolution, and adaptation during a chronic >16-year Burkholderia pseudomallei infection

Although acute melioidosis is the most common outcome of Burkholderia pseudomallei infection, we have documented a case, P314, where disease severity lessened with time, and the pathogen evolved towards a commensal relationship with the host. In the current study, we used whole-genome sequencing to monitor this long-term symbiotic relationship to better understand B. pseudomallei persistence in P314’s sputum despite intensive initial therapeutic regimens. We collected and sequenced 118 B. pseudomallei isolates from P314’s airways over a >16-year period, and also sampled the patient’s home environment, recovering six closely related B. pseudomallei isolates from the household water system. Using comparative genomics, we identified 126 SNPs in the core genome of the 124 isolates or 162 SNPs/indels when the accessory genome was included. The core SNPs were used to construct a phylogenetic tree, which demonstrated a close relationship between environmental and clinical isolates and detailed within-host evolutionary patterns. The phylogeny had little homoplasy, consistent with a strictly clonal mode of genetic inheritance. Repeated sampling revealed evidence of genetic diversification, but frequent extinctions left only one successful lineage through the first four years and two lineages after that. Overall, the evolution of this population is nonadaptive and best explained by genetic drift. However, some genetic and phenotypic changes are consistent with in situ adaptation. Using a mouse model, P314 isolates caused greatly reduced morbidity and mortality compared to the environmental isolates. Additionally, potentially adaptive phenotypes emerged and included differences in the O-antigen, capsular polysaccharide, motility, and colony morphology. The >13-year co-existence of two long-lived lineages presents interesting hypotheses that can be tested in future studies to provide additional insights into selective pressures, niche differentiation, and microbial adaptation. This unusual melioidosis case presents a rare example of the evolutionary progression towards commensalism by a highly virulent pathogen within a single human host.

Pathogens frequently jump between different hosts, and associated adaptation may lead to the emergence of new infectious agents. Such host-jumping evolution is witnessed through endpoint analyses but these cannot capture genetic changes in lineages that have gone extinct. In this study, we have identified and monitored an example of the evolution of a bacterium often deadly to its mammalian host, in an unprecedented case whereby disease lessened through time and the pathogen became a part of the commensal human flora. We used genomic analyses to characterize more than 16 years of this evolutionary process and the stepwise mutations that control pathogen interactions with the patient. Soon after infection, mutational changes occurred that allowed the bacterium to remain in the airways without causing disease. This shift towards avirulence was determined based on clinical data and virulence testing in an animal model. In addition, mutations occurred that contributed to the persistence of the bacteria in the patient's lungs. Finally, we found evidence for the evolutionary emergence and persistence of two distinct lineages of the bacterium over the last 13 years, presenting interesting questions about niche utilization. Bacteria are ubiquitous in the human body and almost all are beneficial or benign. In this study, we document the evolutionary conversion of a normally deadly bacterium towards a commensal.

Genomic epizootiology of a Brucella abortus outbreak in Northern Ireland (1997-2012)

BACKGROUND

In the recent past (1997-2012), Northern Ireland in the United Kingdom suffered an outbreak of Brucella abortus, which at its height affected over 200 cattle herds. Initially, isolates were characterized using multi-locus variable number tandem repeats analysis (MLVA). While informative in this setting, hyper-variability in some loci limited the resolution necessary to infer fine-scale disease transmission networks. Consequently, we applied whole-genome sequencing to isolates from this outbreak to evaluate higher resolution markers for disease epizootiology.

RESULTS

Phylogenetic analysis revealed that the B. abortus outbreak in Northern Ireland was caused by two distinct pathogen lineages. One contained isolates consistent with the 1997-2012 outbreak being linked to a previous endemic infection thought eradicated. The dominant second lineage exhibited little genetic diversity throughout the recrudescent outbreak, with limited population sub-structure evident. This finding was inconsistent with prior MLVA molecular characterizations that suggested the presence of seven clonal complexes. Spatio-temporal modeling revealed a significant association of pairwise SNP differences between isolates and geographic distances. However, effect sizes were very small due to reduced pathogen diversity.

CONCLUSIONS

Genome sequence data suggested that hyper-variability in some MLVA loci contributed to an overestimate of pathogen diversity in the most recent outbreak. The low diversity observed in our genomic dataset made it inappropriate to apply phylodynamic methods to these data. We conclude that maintaining data repositories of genome sequence data will be invaluable for source attribution/epizootiological inference should recrudescence ever re-occur. However genomic epizootiological methods may have limited utility in some settings, such as when applied to recrudescent/re-emergent infections of slowly-evolving bacterial pathogens.

Structure, spatial dynamics, and stability of novel seed dispersal mutualistic networks in Hawaiʻi

Increasing rates of human-caused species invasions and extinctions may reshape communities and modify the structure, dynamics, and stability of species interactions. To investigate how such changes affect communities, we performed multiscale analyses of seed dispersal networks on Oʻahu, Hawaiʻi. Networks consisted exclusively of novel interactions, were largely dominated by introduced species, and exhibited specialized and modular structure at local and regional scales, despite high interaction dissimilarity across communities. Furthermore, the structure and stability of the novel networks were similar to native-dominated communities worldwide. Our findings suggest that shared evolutionary history is not a necessary process for the emergence of complex network structure, and interaction patterns may be highly conserved, regardless of species identity and environment. Introduced species can quickly become well integrated into novel networks, making restoration of native ecosystems more challenging than previously thought.

Genomic sequencing is required for identification of tuberculosis transmission in Hawaii

Background

Tuberculosis (TB) caused an estimated 1.4 million deaths and 10.4 million new cases globally in 2015. TB rates in the United States continue to steadily decline, yet rates in the State of Hawaii are perennially among the highest in the nation due to a continuous influx of immigrants from the Western Pacific and Asia. TB in Hawaii is composed of a unique distribution of genetic lineages, with the Beijing and Manila families of Mycobacterium tuberculosis (Mtb) comprising over two-thirds of TB cases. Standard fingerprinting methods (spoligotyping plus 24-loci Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats [MIRU-VNTR] fingerprinting) perform poorly when used to identify actual transmission clusters composed of isolates from these two families. Those typing methods typically group isolates from these families into large clusters of non-linked isolates with identical fingerprints. Next-generation whole-genome sequencing (WGS) provides a new tool for molecular epidemiology that can resolve clusters of isolates with identical spoligotyping and MIRU-VNTR fingerprints.

Methods

We performed WGS and SNP analysis and evaluated epidemiological data to investigate 19 apparent TB transmission clusters in Hawaii from 2003 to 2017 in order to assess WGS’ ability to resolve putative Mtb clusters from the Beijing and Manila families. This project additionally investigated MIRU-VNTR allele prevalence to determine why standard Mtb fingerprinting fails to usefully distinguish actual transmission clusters from these two Mtb families.

Results

WGS excluded transmission events in seven of these putative clusters, confirmed transmission in eight, and identified both transmission-linked and non-linked isolates in four. For epidemiologically identified clusters, while the sensitivity of MIRU-VNTR fingerprinting for identifying actual transmission clusters was found to be 100%, its specificity was only 28.6% relative to WGS. We identified that the Beijing and Manila families’ significantly lower Shannon evenness of MIRU-VNTR allele distributions than lineage 4 was the cause of standard fingerprinting’s poor performance when identifying transmission in Beijing and Manila family clusters.

Conclusions

This study demonstrated that WGS is necessary for epidemiological investigation of TB in Hawaii and the Pacific.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3502-1) contains supplementary material, which is available to authorized users.

Cryptic connections illuminate pathogen transmission within community networks

Understanding host interactions that lead to pathogen transmission is fundamental to the prediction and control of epidemics^1-5. Although the majority of transmissions often occurs within social groups^6-9, the contribution of connections that bridge groups and species to pathogen dynamics is poorly understood^10-12. These cryptic connections-which are often indirect or infrequent-provide transmission routes between otherwise disconnected individuals and may have a key role in large-scale outbreaks that span multiple populations or species. Here we quantify the importance of cryptic connections in disease dynamics by simultaneously characterizing social networks and tracing transmission dynamics of surrogate-pathogen epidemics through eight communities of bats. We then compared these data to the invasion of the fungal pathogen that causes white-nose syndrome, a recently emerged disease that is devastating North American bat populations^13-15. We found that cryptic connections increased links between individuals and between species by an order of magnitude. Individuals were connected, on average, to less than two per cent of the population through direct contact and to only six per cent through shared groups. However, tracing surrogate-pathogen dynamics showed that each individual was connected to nearly fifteen per cent of the population, and revealed widespread transmission between solitarily roosting individuals as well as extensive contacts among species. Connections estimated from surrogate-pathogen epidemics, which include cryptic connections, explained three times as much variation in the transmission of the fungus that causes white-nose syndrome as did connections based on shared groups. These findings show how cryptic connections facilitate the community-wide spread of pathogens and can lead to explosive epidemics.

Pathogen dynamics during invasion and establishment of white-nose syndrome explain mechanisms of host persistence

Disease dynamics during pathogen invasion and establishment determine the impacts of disease on host populations and determine the mechanisms of host persistence. Temporal progression of prevalence and infection intensity illustrate whether tolerance, resistance, reduced transmission, or demographic compensation allow initially declining populations to persist. We measured infection dynamics of the fungal pathogen Pseudogymnoascus destructans that causes white-nose syndrome in bats by estimating pathogen prevalence and load in seven bat species at 167 hibernacula over a decade as the pathogen invaded, became established, and some host populations stabilized. Fungal loads increased rapidly and prevalence rose to nearly 100% at most sites within 2 yr of invasion in six of seven species. Prevalence and loads did not decline over time despite huge reductions in colony sizes, likely due to an extensive environmental reservoir. However, there was substantial variation in fungal load among sites with persisting colonies, suggesting that both tolerance and resistance developed at different sites in the same species. In contrast, one species disappeared from hibernacula within 3 yr of pathogen invasion. Variable host responses to pathogen invasion require different management strategies to prevent disease-induced extinction and to facilitate evolution of tolerance or resistance in persisting populations.

Whole genome SNP analysis suggests unique virulence factor differences of the Beijing and Manila families of Mycobacterium tuberculosis found in Hawaii

While tuberculosis (TB) remains a global disease, the WHO estimates that 62% of the incident TB cases in 2016 occurred in the WHO South-East Asia and Western Pacific regions. TB in the Pacific is composed predominantly of two genetic families of Mycobacterium tuberculosis (Mtb): Beijing and Manila. The Manila family is historically under-studied relative to the families that comprise the majority of TB in Europe and North America (e.g. lineage 4), and it remains unclear why this lineage has persisted in Filipino populations despite the predominance of more globally successful Mtb lineages in most of the world. The Beijing family is of particular interest as it is increasingly associated with drug resistance throughout the world. Both of these lineages are important to the State of Hawaii, where they comprise over two-thirds of TB cases. Here, we performed whole genome sequencing on 82 Beijing family, Manila family, and outgroup clinical Mtb isolates from Hawaii to identify lineage-specific SNPs (SNPs found in all isolates from their respective families, and exclusively in those families) in established virulence factor genes. Six non-silent lineage-specific virulence factor SNPs were found in the Beijing family, including mutations in alternative sigma factor sigG and polyketide synthases pks5 and pks7. The Manila family displayed more than eleven non-silent lineage-specific and characteristic virulence factor mutations, including in genes coding for MCE-family protein Mce1B, two mutations in fatty-acid-AMP ligase FadD26, and virulence-regulating transcriptional regulator VirS. This study further identified an ancient clade that shared some virulence factor mutations with the Manila family, and investigated the relationship of those and other “Manila-like” spoligotypes to the Manila family with this SNP dataset. This work identified a set of virulence genes that are worth pursuing to determine potential differences in transmission or virulence displayed by these two Mtb families.

White-Nose Syndrome Fungus in a 1918 Bat Specimen from France

White-nose syndrome, first diagnosed in North America in 2006, causes mass deaths among bats in North America. We found the causative fungus, Pseudogymnoascusdestructans, in a 1918 sample collected in Europe, where bats have now adapted to the fungus. These results are consistent with a Eurasian origin of the pathogen.

African Lineage Brucella melitensis Isolates from Omani Livestock

Brucellosis is a common livestock disease in the Middle East and North Africa, but remains poorly described in the region both genetically and epidemiologically. Traditionally found in goats and sheep, Brucella melitensis is increasingly recognized as infecting camels. Most studies of brucellosis in camels to date have focused on serological surveys, providing only limited understanding of the molecular epidemiology of circulating strains. We genotyped B. melitensis isolates from Omani camels using whole genome SNP assays and VNTRs to provide context for regional brucellosis cases. We identified a lineage of B. melitensis circulating in camels as well as in goats, sheep, and cattle in Oman. This lineage is genetically distinct from most genotypes from the Arabian Peninsula and from isolates from much of the rest of the Middle East. We then developed diagnostic assays that rapidly identify strains from this lineage. In analyses of genotypes from throughout the region, Omani isolates were genetically most closely related to strains from brucellosis cases in humans and livestock in North Africa. Our findings suggest an African origin for B. melitensis in Oman that has likely occurred through the trade of infected livestock. Moreover, African lineages of B. melitensis appear to be undersampled and consequently are underrepresented in genetic databases for Brucella. As we begin to more fully understand global genomic diversity of B. melitensis, finding and characterizing these unique but widespread lineages is essential. We predict that increased sampling of humans and livestock in Africa will reveal little known diversity in this important zoonotic pathogen.

Extreme sensitivity to ultraviolet light in the fungal pathogen causing white-nose syndrome of bats

Bat white-nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans, has decimated North American hibernating bats since its emergence in 2006. Here, we utilize comparative genomics to examine the evolutionary history of this pathogen in comparison to six closely related nonpathogenic species. P. destructans displays a large reduction in carbohydrate-utilizing enzymes (CAZymes) and in the predicted secretome (~50%), and an increase in lineage-specific genes. The pathogen has lost a key enzyme, UVE1, in the alternate excision repair (AER) pathway, which is known to contribute to repair of DNA lesions induced by ultraviolet (UV) light. Consistent with a nonfunctional AER pathway, P. destructans is extremely sensitive to UV light, as well as the DNA alkylating agent methyl methanesulfonate (MMS). The differential susceptibility of P. destructans to UV light in comparison to other hibernacula-inhabiting fungi represents a potential “Achilles’ heel” of P. destructans that might be exploited for treatment of bats with WNS.

White-nose syndrome, caused by the fungus Pseudogymnoascus destructans, is decimating North American bats. Here, Palmer et al. use comparative genomics to examine the evolutionary history of this pathogen, and show that it has lost a crucial DNA repair enzyme and is extremely sensitive to UV light.

Resistance in persisting bat populations after white-nose syndrome invasion

Increases in anthropogenic movement have led to a rise in pathogen introductions and the emergence of infectious diseases in naive host communities worldwide. We combined empirical data and mathematical models to examine changes in disease dynamics in little brown bat (Myotis lucifugus) populations following the introduction of the emerging fungal pathogen Pseudogymnoascus destructans, which causes the disease white-nose syndrome. We found that infection intensity was much lower in persisting populations than in declining populations where the fungus has recently invaded. Fitted models indicate that this is most consistent with a reduction in the growth rate of the pathogen when fungal loads become high. The data are inconsistent with the evolution of tolerance or an overall reduced pathogen growth rate that might be caused by environmental factors. The existence of resistance in some persisting populations of little brown bats offers a glimmer of hope that a precipitously declining species will persist in the face of this deadly pathogen.This article is part of the themed issue 'Human influences on evolution, and the ecological and societal consequences'.

Drivers of variation in species impacts for a multi-host fungal disease of bats

Disease can play an important role in structuring species communities because the effects of disease vary among hosts; some species are driven towards extinction, while others suffer relatively little impact. Why disease impacts vary among host species remains poorly understood for most multi-host pathogens, and factors allowing less-susceptible species to persist could be useful in conserving highly affected species. White-nose syndrome (WNS), an emerging fungal disease of bats, has decimated some species while sympatric and closely related species have experienced little effect. We analysed data on infection prevalence, fungal loads and environmental factors to determine how variation in infection among sympatric host species influenced the severity of WNS population impacts. Intense transmission resulted in almost uniformly high prevalence in all species. By contrast, fungal loads varied over 3 orders of magnitude among species, and explained 98% of the variation among species in disease impacts. Fungal loads increased with hibernating roosting temperatures, with bats roosting at warmer temperatures having higher fungal loads and suffering greater WNS impacts. We also found evidence of a threshold fungal load, above which the probability of mortality may increase sharply, and this threshold was similar for multiple species. This study demonstrates how differences in behavioural traits among species-in this case microclimate preferences-that may have been previously adaptive can be deleterious after the introduction of a new pathogen. Management to reduce pathogen loads rather than exposure may be an effective way of reducing disease impact and preventing species extinctions.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Brucella spp. of amphibians comprise genomically diverse motile strains competent for replication in macrophages and survival in mammalian hosts

Twenty-one small Gram-negative motile coccobacilli were isolated from 15 systemically diseased African bullfrogs (Pyxicephalus edulis), and were initially identified as Ochrobactrum anthropi by standard microbiological identification systems. Phylogenetic reconstructions using combined molecular analyses and comparative whole genome analysis of the most diverse of the bullfrog strains verified affiliation with the genus Brucella and placed the isolates in a cluster containing B. inopinata and the other non-classical Brucella species but also revealed significant genetic differences within the group. Four representative but molecularly and phenotypically diverse strains were used for in vitro and in vivo infection experiments. All readily multiplied in macrophage-like murine J774-cells, and their overall intramacrophagic growth rate was comparable to that of B. inopinata BO1 and slightly higher than that of B. microti CCM 4915. In the BALB/c murine model of infection these strains replicated in both spleen and liver, but were less efficient than B. suis 1330. Some strains survived in the mammalian host for up to 12 weeks. The heterogeneity of these novel strains hampers a single species description but their phenotypic and genetic features suggest that they represent an evolutionary link between a soil-associated ancestor and the mammalian host-adapted pathogenic Brucella species.

Host persistence or extinction from emerging infectious disease: insights from white-nose syndrome in endemic and invading regions

Predicting species' fates following the introduction of a novel pathogen is a significant and growing problem in conservation. Comparing disease dynamics between introduced and endemic regions can offer insight into which naive hosts will persist or go extinct, with disease acting as a filter on host communities. We examined four hypothesized mechanisms for host-pathogen persistence by comparing host infection patterns and environmental reservoirs for Pseudogymnoascus destructans (the causative agent of white-nose syndrome) in Asia, an endemic region, and North America, where the pathogen has recently invaded. Although colony sizes of bats and hibernacula temperatures were very similar, both infection prevalence and fungal loads were much lower on bats and in the environment in Asia than North America. These results indicate that transmission intensity and pathogen growth are lower in Asia, likely due to higher host resistance to pathogen growth in this endemic region, and not due to host tolerance, lower transmission due to smaller populations, or lower environmentally driven pathogen growth rate. Disease filtering also appears to be favouring initially resistant species in North America. More broadly, determining the mechanisms allowing species persistence in endemic regions can help identify species at greater risk of extinction in introduced regions, and determine the consequences for disease dynamics and host-pathogen coevolution.

Deconstructing the Bat Skin Microbiome: Influences of the Host and the Environment

Bats are geographically widespread and play an important role in many ecosystems, but relatively little is known about the ecology of their associated microbial communities and the role microbial taxa play in bat health, development, and evolution. Moreover, few vertebrate animal skin microbiomes have been comprehensively assessed, and thus characterizing the bat skin microbiome will yield valuable insight into the variability of vertebrate skin microbiomes as a whole. The recent emergence of the skin fungal disease white-nose syndrome highlights the potentially important role bat skin microbial communities could play in bat health. Understanding the determinant of bat skin microbial communities could provide insight into important factors allowing individuals to persist with disease. We collected skin swabs from a total of 11 bat species from the eastern United States (n = 45) and Colorado (n = 119), as well as environmental samples (n = 38) from a subset of sites, and used 16S rRNA marker gene sequencing to observe bacterial communities. In addition, we conducted a literature survey to compare the skin microbiome across vertebrate groups, including the bats presented in this study. Host species, region, and site were all significant predictors of the variability across bat skin bacterial communities. Many bacterial taxa were found both on bats and in the environment. However, some bacterial taxa had consistently greater relative abundances on bat skin relative to their environments. Bats shared many of their abundant taxa with other vertebrates, but also hosted unique bacterial lineages such as the class Thermoleophilia (Actinobacteria). A strong effect of site on the bat skin microbiome indicates that the environment very strongly influences what bacteria are present on bat skin. Bat skin microbiomes are largely composed of site-specific microbiota, but there do appear to be important host-specific taxa. How this translates to differences in host-microbial interactions and bat health remains an important knowledge gap, but this work suggests that habitat variability is very important. We identify some bacterial groups that are more consistent on bats despite site differences, and these may be important ones to study in terms of their function as potential core microbiome members.

VNTR diversity in Yersinia pestis isolates from an animal challenge study reveals the potential for in vitro mutations during laboratory cultivation

Underlying mutation rates and other evolutionary forces shape the population structure of bacteria in nature. Although easily overlooked, similar forces are at work in the laboratory and may influence observed mutations. Here, we investigated tissue samples and Yersinia pestis isolates from a rodent laboratory challenge with strain CO92 using whole genome sequencing and multi-locus variable-number tandem repeat (VNTR) analysis (MLVA). We identified six VNTR mutations that were found to have occurred in vitro during laboratory cultivation rather than in vivo during the rodent challenge. In contrast, no single nucleotide polymorphism (SNP) mutations were observed, either in vivo or in vitro. These results were consistent with previously published mutation rates and the calculated number of Y. pestis generations that occurred during the in vitro versus the in vivo portions of the experiment. When genotyping disease outbreaks, the potential for in vitro mutations should be considered, particularly when highly variable genetic markers such as VNTRs are used.

Genomics reveals historic and contemporary transmission dynamics of a bacterial disease among wildlife and livestock

Whole-genome sequencing has provided fundamental insights into infectious disease epidemiology, but has rarely been used for examining transmission dynamics of a bacterial pathogen in wildlife. In the Greater Yellowstone Ecosystem (GYE), outbreaks of brucellosis have increased in cattle along with rising seroprevalence in elk. Here we use a genomic approach to examine Brucella abortus evolution, cross-species transmission and spatial spread in the GYE. We find that brucellosis was introduced into wildlife in this region at least five times. The diffusion rate varies among Brucella lineages (∼3 to 8 km per year) and over time. We also estimate 12 host transitions from bison to elk, and 5 from elk to bison. Our results support the notion that free-ranging elk are currently a self-sustaining brucellosis reservoir and the source of livestock infections, and that control measures in bison are unlikely to affect the dynamics of unrelated strains circulating in nearby elk populations.

White-Nose Syndrome Disease Severity and a Comparison of Diagnostic Methods

White-nose syndrome is caused by the fungus Pseudogymnoascus destructans and has killed millions of hibernating bats in North America but the pathophysiology of the disease remains poorly understood. Our objectives were to (1) assess non-destructive diagnostic methods for P. destructans infection compared to histopathology, the current gold-standard, and (2) to evaluate potential metrics of disease severity. We used data from three captive inoculation experiments involving 181 little brown bats (Myotis lucifugus) to compare histopathology, quantitative PCR (qPCR), and ultraviolet fluorescence as diagnostic methods of P. destructans infection. To assess disease severity, we considered two histology metrics (wing area with fungal hyphae, area of dermal necrosis), P. destructans fungal load (qPCR), ultraviolet fluorescence, and blood chemistry (hematocrit, sodium, glucose, pCO2, and bicarbonate). Quantitative PCR was most effective for early detection of P. destructans, while all three methods were comparable in severe infections. Correlations among hyphae and necrosis scores, qPCR, ultraviolet fluorescence, blood chemistry, and hibernation duration indicate a multi-stage pattern of disease. Disruptions of homeostasis occurred rapidly in late hibernation. Our results provide valuable information about the use of non-destructive techniques for monitoring, and provide novel insight into the pathophysiology of white-nose syndrome, with implications for developing and implementing potential mitigation strategies.