AbstractMass mortality events are increasing globally in frequency and magnitude, largely as a result of human-induced change. The effects of these mass mortality events, in both the long and short term, are of imminent concern because of their ecosystem impacts. Genomic data can be used to reveal some of the population-level changes associated with mass mortality events. Here, we use reduced-representation sequencing to identify potential short-term genetic impacts of a mass mortality event associated with a sea star wasting outbreak. We tested for changes in the population for genetic differentiation, diversity, and effective population size between pre-sea star wasting and post-sea star wasting populations of Pisaster ochraceus-a species that suffered high sea star wasting-associated mortality (75%-100% at 80% of sites). We detected no significant population-based genetic differentiation over the spatial scale sampled; however, the post-sea star wasting population tended toward more differentiation across sites than the pre-sea star wasting population. Genetic estimates of effective population size did not detectably change, consistent with theoretical expectations; however, rare alleles were lost. While we were unable to detect significant population-based genetic differentiation or changes in effective population size over this short time period, the genetic burden of this mass mortality event may be borne by future generations, unless widespread recruitment mitigates the population decline. Prior results from P. ochraceus indicated that natural selection played a role in altering allele frequencies following this mass mortality event. In addition to the role of selection found in a previous study on the genomic impacts of sea star wasting on P. ochraceus, our current study highlights the potential role the stochastic loss of many individuals plays in altering how genetic variation is structured across the landscape. Future genetic monitoring is needed to determine long-term genetic impacts in this long-lived species. Given the increased frequency of mass mortality events, it is important to implement demographic and genetic monitoring strategies that capture baselines and background dynamics to better contextualize species responses to large perturbations.

Sea star wasting (SSW) disease describes a condition affecting asteroids that resulted in significant Northeastern Pacific population decline following a mass mortality event in 2013. The etiology of SSW is unresolved. We hypothesized that SSW is a sequela of microbial organic matter remineralization near respiratory surfaces, one consequence of which may be limited O₂ availability at the animal-water interface. Microbial assemblages inhabiting tissues and at the asteroid-water interface bore signatures of copiotroph proliferation before SSW onset, followed by the appearance of putatively facultative and strictly anaerobic taxa at the time of lesion genesis and as animals died. SSW lesions were induced in Pisaster ochraceus by enrichment with a variety of organic matter (OM) sources. These results together illustrate that depleted O₂ conditions at the animal-water interface may be established by heterotrophic microbial activity in response to organic matter loading. SSW was also induced by modestly (∼39%) depleted O₂ conditions in aquaria, suggesting that small perturbations in dissolved O₂ may exacerbate the condition. SSW susceptibility between species was significantly and positively correlated with surface rugosity, a key determinant of diffusive boundary layer thickness. Tissues of SSW-affected individuals collected in 2013-2014 bore δ¹⁵N signatures reflecting anaerobic processes, which suggests that this phenomenon may have affected asteroids during mass mortality at the time. The impacts of enhanced microbial activity and subsequent O₂ diffusion limitation may be more pronounced under higher temperatures due to lower O₂ solubility, in more rugose asteroid species due to restricted hydrodynamic flow, and in larger specimens due to their lower surface area to volume ratios which affects diffusive respiratory potential.