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.

Freshwater mussels of the order Unionida are among the most endangered animal groups globally, but the causes of their population decline are often enigmatic, with little known about the role of disease. In 2018, we collected wild adult pheasantshell (Actinonaias pectorosa) and mucket (Actinonaias ligamentina) during an epidemiologic survey investigating an ongoing mussel mass mortality event in the Clinch River, Virginia and Tennessee, USA. Histopathology and transmission electron microscopy showed a novel microsporidian parasite primarily infecting the ovary of pheasantshell. Sequencing of the small subunit rRNA gene produced a 1333 bp sequence with the greatest similarity to Pseudonosema cristatellae (AF484694.1; 86.36%; e-value = 0), a microsporidium infecting the freshwater bryozoan (Cristatella mucedo). Microsporidia were observed in 65% (17/26) of the examined female pheasantshell (A. pectorosa) and in no (0/2) female muckets (A. ligamentina) and occurred at mortality and non-mortality sites. Our findings indicate that a novel parasite, Microsporidium clinchi n. sp., is present in pheasantshell in the Clinch River, and while likely not a cause of mass mortality, could reduce fecundity and recruitment in this declining population and threaten the success of reintroductions. Surveillance of M. clinchi n. sp. and evaluation of broodstock and their progeny for microsporidia would therefore be prudent.

Freshwater mussels (order Unionida) play a key role in freshwater systems as ecosystem engineers and indicators of aquatic ecosystem health. The fauna is globally imperilled due to a diversity of suspected factors; however, causes for many population declines and mortality events remain unconfirmed due partly to limited health assessment tools. Mussel-monitoring activities often rely on population-level measurements, such as abundance and age structure, which reflect delayed responses to environmental conditions. Measures of organismal health would enable preemptive detection of declining condition before population-level effects manifest. Metabolomic analysis can identify shifts in biochemical pathways in response to stressors and changing environmental conditions; however, interpretation of the results requires information on inherent variability of metabolite concentrations in mussel populations. We targeted metabolites in the haemolymph of two common mussels, Lampsilis cardium and Lampsilis siliquoidea, from three Indiana streams (USA) using ultra-high-performance liquid chromatography combined with quadrupole time-of-flight mass spectroscopy. The influence of species, stream and sex on metabolite variability was examined with distance-based redundancy analysis. Metabolite variability was most influenced by species, followed by site and sex. Inter- and intraspecies metabolite variability among sexes was less distinct than differences among locations. We further categorized metabolites by occurrence and variability in mussel populations. Metabolites with high occurrence (Categories 1 and 2) included those indicative of energy status (catabolism versus anabolism; arginine, proline, carnitine, nicotinic acid, pantothenic acid), oxidative stress (proline, glutamine, glutamate) and protein metabolism (thymidine, cytidine, inosine). Metabolites with lower occurrence (Category 3) are constituents of assorted metabolic pathways and can be important biomarkers with additional temporal sampling to characterize their variability. These data provide a reference for future temporal (before/after) monitoring and for studies of stressor–metabolite linkages in freshwater mussels.

Biologists monitoring freshwater mussel (order Unionida) populations rely on behavioral, often subjective, signs to identify moribund (“sick”) or stressed mussels, such as gaping valves and slow response to probing, and they lack clinical indicators to support a diagnosis. As part of a multi-year study to investigate causes of reoccurring mortality of pheasantshell (Ortmanniana pectorosa; synonym Actinonaias pectorosa) in the Clinch River, Virginia and Tennessee, USA, we analyzed the hemolymph metabolome of a subset of mussels from the 2018 sampling period. Mussels at the mortality sites were diagnosed in the field as affected (case) or unaffected (control) based on behavioral and physical signs. Hemolymph was collected in the field by non-lethal methods from the anterior adductor muscle for analysis. We used ultra-high-performance liquid chromatography with quadrupole time-of-flight mass spectroscopy to detect targeted and untargeted metabolites in hemolymph and compared metabolomic profiles by field assessment of clinical status. Targeted biomarker analysis found 13 metabolites associated with field assessments of clinical status. Of these, increased gamma-linolenic acid and N-methyl-l-alanine were most indicative of case mussels, while adenine and inosine were the best indicators of control mussels. Five pathways in the targeted analysis differed by clinical status; two of these, purine metabolism and glycerophospholipid metabolism, were also indicated in the untargeted analysis. In the untargeted nalysis, 22 metabolic pathways were associated with clinical status. Many of the impacted pathways in the case group were catabolic processes, such as degradation of amino acids and fatty acids. Hierarchical clustering analysis matched clinical status in 72% (18 of 25) of mussels, with control mussels more frequently (5 of 16) not matching clinical status. Our study demonstrated that metabolomic analysis of hemolymph is suitable for assessing mussel condition and complements field-based indicators of health.