No abstract

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.

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 morality 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.