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.