A novel tamanavirus (Flaviviridae) of the European common frog (Rana temporaria) from the UK.

Flavivirids are small, enveloped, positive-sense RNA viruses from the family Flaviviridae with genomes of ~9-13 kb. Metatranscriptomic analyses of metazoan organisms have revealed a diversity of flavivirus-like or flavivirid viral sequences in fish and marine invertebrate groups. However, no flavivirus-like virus has been identified in amphibians. To remedy this, we investigated the virome of the European common frog (Rana temporaria) in the UK, utilizing high-throughput sequencing at six catch locations. De novo assembly revealed a coding-complete virus contig of a novel flavivirid ~11.2 kb in length. The virus encodes a single ORF of 3456 aa and 5' and 3' untranslated regions (UTRs) of 227 and 666 nt, respectively. We named this virus Rana tamanavirus (RaTV), as BLASTp analysis of the polyprotein showed the closest relationships to Tamana bat virus (TABV) and Cyclopterus lumpus virus from Pteronotus parnellii and Cyclopterus lumpus, respectively. Phylogenetic analysis of the RaTV polyprotein compared to Flavivirus and Flavivirus-like members indicated that RaTV was sufficiently divergent and basal to the vertebrate Tamanavirus clade. In addition to the Mitcham strain, partial but divergent RaTV, sharing 95.64-97.39 % pairwise nucleotide identity, were also obtained from the Poole and Deal samples, indicating that RaTV is widespread in UK frog samples. Bioinformatic analyses of predicted secondary structures in the 3'UTR of RaTV showed the presence of an exoribonuclease-resistant RNA (xrRNA) structure standard in flaviviruses and TABV. To examine this biochemically, we conducted an in vitro Xrn1 digestion assay showing that RaTV probably forms a functional Xrn1-resistant xrRNA.

Viruses of Atlantic Bonefish (Albula vulpes) in Florida and the Caribbean show geographic patterns consistent with population declines.

Atlantic Bonefish (Albula vulpes) are economically important due to their popularity with recreational anglers. In the State of Florida, USA, bonefish population numbers declined by approximately 60% between the 1990s and 2015. Habitat loss, water quality impairment, chemical inputs, and other anthropogenic factors have been implicated as causes, but the role of pathogens has been largely overlooked, especially with respect to viruses. We used a metagenomic approach to identify and quantify viruses in the blood of 103 A. vulpes sampled throughout their Western Atlantic range, including populations in Florida that have experienced population declines and populations in Belize, Mexico, Puerto Rico, and The Bahamas that have remained apparently stable. We identified four viruses, all of which are members of families known to infect marine fishes (Flaviviridae, Iflaviridae, Narnaviridae, and Nodaviridae), but all of which were previously undescribed. Bonefish from Florida and Mexico had higher viral richness (numbers of distinct viruses per individual fish) than fish sampled from other areas, and bonefish from the Upper Florida Keys had the highest prevalence of viral infection (proportion of positive fish) than fish sampled from any other location. Bonefish from Florida also had markedly higher viral loads than fish sampled from any other area, both for a novel narnavirus and for all viruses combined. Bonefish viruses may be indicators of environmentally driven physiological and immunological compromise, causes of ill health, or both. Supplementary information: The online version contains supplementary material available at 10.1007/s10641-022-01306-9.

Freshwater Mussels Show Elevated Viral Richness and Intensity during a Mortality Event.

Freshwater mussels (Unionida) are among the world's most imperiled taxa, but the relationship between freshwater mussel mortality events and infectious disease is largely unstudied. We surveyed viromes of a widespread and abundant species (mucket, Actinonaias ligamentina; syn: Ortmanniana ligamentina) experiencing a mortality event of unknown etiology in the Huron River, Michigan, in 2019-2020 and compared them to viromes from mucket in a healthy population in the St. Croix River, Wisconsin and a population from the Clinch River, Virginia and Tennessee, where a mortality event was affecting the congeneric pheasantshell (Actinonaias pectorosa; syn: Ortmanniana pectorosa) population. We identified 38 viruses, most of which were associated with mussels collected during the Huron River mortality event. Viral richness and cumulative viral read depths were significantly higher in moribund mussels from the Huron River than in healthy controls from each of the three populations. Our results demonstrate significant increases in the number and intensity of viral infections for freshwater mussels experiencing mortality events, whereas individuals from healthy populations have a substantially reduced virome comprising a limited number of species at low viral read depths.

Mussel Mass Mortality and the Microbiome: Evidence for Shifts in the Bacterial Microbiome of a Declining Freshwater Bivalve.

Freshwater mussels (Unionida) are suffering mass mortality events worldwide, but the causes remain enigmatic. Here, we describe an analysis of bacterial loads, community structure, and inferred metabolic pathways in the hemolymph of pheasantshells (Actinonaias pectorosa) from the Clinch River, USA, during a multi-year mass mortality event. Bacterial loads were approximately 2 logs higher in moribund mussels (cases) than in apparently healthy mussels (controls). Bacterial communities also differed between cases and controls, with fewer sequence variants (SVs) and higher relative abundances of the proteobacteria Yokenella regensburgei and Aeromonas salmonicida in cases than in controls. Inferred bacterial metabolic pathways demonstrated a predominance of degradation, utilization, and assimilation pathways in cases and a predominance of biosynthesis pathways in controls. Only two SVs correlated with Clinch densovirus 1, a virus previously shown to be strongly associated with mortality in this system: Deinococcota and Actinobacteriota, which were associated with densovirus-positive and densovirus-negative mussels, respectively. Overall, our results suggest that bacterial invasion and shifts in the bacterial microbiome during unionid mass mortality events may result from primary insults such as viral infection or environmental stressors. If so, bacterial communities in mussel hemolymph may be sensitive, if generalized, indicators of declining mussel health.

Soil Reservoir Dynamics of Ophidiomyces ophidiicola, the Causative Agent of Snake Fungal Disease.

Wildlife diseases pose an ever-growing threat to global biodiversity. Understanding how wildlife pathogens are distributed in the environment and the ability of pathogens to form environmental reservoirs is critical to understanding and predicting disease dynamics within host populations. Snake fungal disease (SFD) is an emerging conservation threat to North American snake populations. The causative agent, Ophidiomyces ophidiicola (Oo), is detectable in environmentally derived soils. However, little is known about the distribution of Oo in the environment and the persistence and growth of Oo in soils. Here, we use quantitative PCR to detect Oo in soil samples collected from five snake dens. We compare the detection rates between soils collected from within underground snake hibernacula and associated, adjacent topsoil samples. Additionally, we used microcosm growth assays to assess the growth of Oo in soils and investigate whether the detection and growth of Oo are related to abiotic parameters and microbial communities of soil samples. We found that Oo is significantly more likely to be detected in hibernaculum soils compared to topsoils. We also found that Oo was capable of growth in sterile soil, but no growth occurred in soils with an active microbial community. A number of fungal genera were more abundant in soils that did not permit growth of Oo, versus those that did. Our results suggest that soils may display a high degree of both general and specific suppression of Oo in the environment. Harnessing environmental suppression presents opportunities to mitigate the impacts of SFD in wild snake populations.

Mycobiome Traits Associated with Disease Tolerance Predict Many Western North American Bat Species Will Be Susceptible to White-Nose Syndrome.

White-nose syndrome (WNS), a fungal disease that has caused catastrophic population declines of bats in eastern North America, is rapidly spreading across the continent and now threatens previously unexposed bat species in western North America. The causal agent of WNS, the fungus Pseudogymnoascus destructans, can infect many species of hibernating bats, but susceptibility to WNS varies by host species. We previously reported that certain traits of the skin microbiome, particularly yeast diversity and abundance, of bat species in eastern North America are strongly associated with resistance to WNS. Using these traits, we developed models to predict WNS susceptibility of 13 species of western North American bats. Based on models derived from yeast species diversity, only one bat species, Myotis velifer, was predicted to be WNS resistant (i.e., may develop the disease, but with low mortality rates). We also screened yeasts found on western bats for P. destructans-antagonistic properties by spore germination and growth inhibition/competition assays and found the ability of yeasts to inhibit P. destructans in vitro to be strain specific. Similar to results of inhibition assays performed with yeasts isolated from bats in eastern North America, few yeasts isolated from bats in western North America inhibited P. destructans in vitro. Continued monitoring of western bat populations will serve to validate the accuracy of the mycobiome analysis in predicting WNS susceptibility, document population and susceptibility trends, and identify additional predictors to assess the vulnerability of naive bat populations to WNS. IMPORTANCE White-nose syndrome is one of the most devastating wildlife diseases ever documented. Some bat species are resistant to or tolerant of the disease, and we previously reported that certain traits of the skin mycobiome of bat species in eastern North America are strongly associated with resistance to WNS. Predicting which western bat species will be most susceptible to WNS would be of great value for establishing conservation priorities. Based on models derived from yeast species diversity, only one bat species was predicted to be WNS resistant. High susceptibility to WNS would pose a significant conservation threat to bats in western North America.

A multi-omic investigation of male lower urinary tract symptoms: Potential role for JC virus.

Male lower urinary tract symptoms (LUTS) comprise a common syndrome of aging that negatively impacts quality of life. The etiology of LUTS is multifactorial, involving benign prostatic hyperplasia, smooth muscle and neurologic dysfunction, inflammation, sexually transmitted infections, fibrosis, and potentially dysbiosis, but this aspect remains poorly explored. We investigated whether the presence of infectious agents in urine might be associated with LUTS by combining next-generation DNA sequencing for virus discovery, microbiome analysis for characterization of bacterial communities, and mass spectrometry-based metabolomics. In urine from 29 LUTS cases and 9 controls from Wisconsin, we found a statistically significant association between a diagnosis of LUTS and the presence of JC virus (JCV), a common neurotropic human polyomavirus (Polyomaviridae, Betapolyomavirus) linked to severe neurologic disease in rare cases. This association (based on metagenomics) was not borne out when specific polymerase chain reaction (PCR) testing was applied to this set of samples, likely due to the greater sensitivity of PCR. Interestingly, urine metabolomics analysis identified dysregulation of metabolites associated with key LUTS processes. Microbiome analysis found no evidence of microbial community dysbiosis in LUTS cases, but JCV-positive samples contained more Anaerococcus species, which are involved in polymicrobial infections of the urinary tract. Neither age nor body mass index were significantly associated with the presence of urinary JCV-in the initial group or in an additional, regionally distinct group. These data provide preliminary support the hypothesis that viruses such as JCV may play a role in the development or progression of LUTS, together with other infectious agents and host metabolic responses.

Skin fungal assemblages of bats vary based on susceptibility to white-nose syndrome.

Microbial skin assemblages, including fungal communities, can influence host resistance to infectious diseases. The diversity-invasibility hypothesis predicts that high-diversity communities are less easily invaded than species-poor communities, and thus diverse microbial communities may prevent pathogens from colonizing a host. To explore the hypothesis that host fungal communities mediate resistance to infection by fungal pathogens, we investigated characteristics of bat skin fungal communities as they relate to susceptibility to the emerging disease white-nose syndrome (WNS). Using a culture-based approach, we compared skin fungal assemblage characteristics of 10 bat species that differ in susceptibility to WNS across 10 eastern U.S. states. The fungal assemblages on WNS-susceptible bat species had significantly lower alpha diversity and abundance compared to WNS-resistant species. Overall fungal assemblage structure did not vary based on WNS-susceptibility, but several yeast species were differentially abundant on WNS-resistant bat species. One yeast species inhibited Pseudogymnoascus destructans (Pd), the causative agent on WNS, in vitro under certain conditions, suggesting a possible role in host protection. Further exploration of interactions between Pd and constituents of skin fungal assemblages may prove useful for predicting susceptibility of bat populations to WNS and for developing effective mitigation strategies.

LONG-TERM SURVIVAL OF PSEUDOGYMNOASCUS DESTRUCTANS AT ELEVATED TEMPERATURES.

White-nose syndrome is an emerging fungal disease that has devastated hibernating bat populations across eastern North America. The causal pathogen, Pseudogymnoascus destructans (PD), is a psychrophilic fungus with a known maximal growth temperature of 20 C. Although it is widely speculated that PD is primarily spread between hibernacula by the movement of bats, experimental evidence is lacking to demonstrate that PD can endure temperatures experienced by active bats for periods of time that would facilitate dispersal of viable fungus. We used an in vitro culture-based approach to study the survival of PD conidia on three artificial growth media and bat fur. The fungus was incubated at three temperatures it might realistically be exposed to on nonhibernating bats or in the environment outside of caves and mines (24 C, 30 C, and 37 C). When incubated on artificial media, we found that PD conidia were able to survive for a maximum of 150 d when exposed to temperatures of 24 C, 60 d at 30 C, and 15 d at 37 C. At all temperatures, maximal survival duration was recorded when conidia were incubated on brain-heart infusion agar with 10% volume of sheep (Ovis aries) blood. When incubated on bat fur, viable PD was recovered at 180 d, 60 d, and 5 d when exposed to temperatures of 24 C, 30 C, and 37 C, respectively. Our results suggest that viable PD conidia may be able to survive on or within the bodies of bats, which may facilitate long-distance dispersal. The long-term viability of the fungus on various fomites may differ, and therefore must be assessed for each potential substrate.

Outbreaks of an Emerging Viral Disease Covary With Differences in the Composition of the Skin Microbiome of a Wild United Kingdom Amphibian.

There is growing appreciation of the important role of commensal microbes in ensuring the normal function and health of their hosts, including determining how hosts respond to pathogens. A range of infectious diseases are threatening amphibians worldwide, and evidence is accumulating that the host-associated bacteria that comprise the microbiome may be key in mediating interactions between amphibian hosts and infectious pathogens. We used 16S rRNA amplicon sequencing to quantify the skin microbial community structure of over 200 individual wild adult European common frogs (Rana temporaria), from ten populations with contrasting history of the lethal disease ranavirosis, caused by emerging viral pathogens belonging to the genus Ranavirus. All populations had similar species richness irrespective of disease history, but populations that have experienced historical outbreaks of ranavirosis have a distinct skin microbiome structure (beta diversity) when compared to sites where no outbreaks of the disease have occurred. At the individual level, neither age, body length, nor sex of the frog could predict the structure of the skin microbiota. Our data potentially support the hypothesis that variation among individuals in skin microbiome structure drive differences in susceptibility to infection and lethal outbreaks of disease. More generally, our results suggest that population-level processes are more important for driving differences in microbiome structure than variation among individuals within populations in key life history traits such as age and body size.

An emerging viral pathogen truncates population age structure in a European amphibian and may reduce population viability.

Infectious diseases can alter the demography of their host populations, reducing their viability even in the absence of mass mortality. Amphibians are the most threatened group of vertebrates globally, and emerging infectious diseases play a large role in their continued population declines. Viruses belonging to the genus Ranavirus are responsible for one of the deadliest and most widespread of these diseases. To date, no work has used individual level data to investigate how ranaviruses affect population demographic structure. We used skeletochronology and morphology to evaluate the impact of ranaviruses on the age structure of populations of the European common frog (Rana temporaria) in the UK. We compared ecologically similar populations that differed most notably in their historical presence or absence of ranavirosis (the acute syndrome caused by ranavirus infection). Our results suggest that ranavirosis may truncate the age structure of R. temporaria populations. One potential explanation for such a shift might be increased adult mortality and subsequent shifts in the life history of younger age classes that increase reproductive output earlier in life. Additionally, we constructed population projection models which indicated that such increased adult mortality could heighten the vulnerability of frog populations to stochastic environmental challenges.

A novel approach to wildlife transcriptomics provides evidence of disease-mediated differential expression and changes to the microbiome of amphibian populations.

Ranaviruses are responsible for a lethal, emerging infectious disease in amphibians and threaten their populations throughout the world. Despite this, little is known about how amphibian populations respond to ranaviral infection. In the United Kingdom, ranaviruses impact the common frog (Rana temporaria). Extensive public engagement in the study of ranaviruses in the UK has led to the formation of a unique system of field sites containing frog populations of known ranaviral disease history. Within this unique natural field system, we used RNA sequencing (RNA-Seq) to compare the gene expression profiles of R. temporaria populations with a history of ranaviral disease and those without. We have applied a RNA read-filtering protocol that incorporates Bloom filters, previously used in clinical settings, to limit the potential for contamination that comes with the use of RNA-Seq in nonlaboratory systems. We have identified a suite of 407 transcripts that are differentially expressed between populations of different ranaviral disease history. This suite contains genes with functions related to immunity, development, protein transport and olfactory reception among others. A large proportion of potential noncoding RNA transcripts present in our differentially expressed set provide first evidence of a possible role for long noncoding RNA (lncRNA) in amphibian response to viruses. Our read-filtering approach also removed significantly more bacterial reads from libraries generated from positive disease history populations. Subsequent analysis revealed these bacterial read sets to represent distinct communities of bacterial species, which is suggestive of an interaction between ranavirus and the host microbiome in the wild.