Revealing the extent of sea otter impacts on bivalve prey through multi-trophic monitoring and mechanistic models.

Sea otters are apex predators that can exert considerable influence over the nearshore communities they occupy. Since facing near extinction in the early 1900s, sea otters are making a remarkable recovery in Southeast Alaska, particularly in Glacier Bay, the largest protected tidewater glacier fjord in the world. The expansion of sea otters across Glacier Bay offers both a challenge to monitoring and stewardship and an unprecedented opportunity to study the top-down effect of a novel apex predator across a diverse and productive ecosystem. Our goal was to integrate monitoring data across trophic levels, space, and time to quantify and map the predator-prey interaction between sea otters and butter clams Saxidomus gigantea, one of the dominant large bivalves in Glacier Bay and a favoured prey of sea otters. We developed a spatially-referenced mechanistic differential equation model of butter clam dynamics that combined both environmental drivers of local population growth and estimates of otter abundance from aerial survey data. We embedded this model in a Bayesian statistical framework and fit it to clam survey data from 43 intertidal and subtidal sites across Glacier Bay. Prior to substantial sea otter expansion, we found that butter clam density was structured by an environmental gradient driven by distance from glacier (represented by latitude) and a quadratic effect of current speed. Estimates of sea otter attack rate revealed spatial heterogeneity in sea otter impacts and a negative relationship with local shoreline complexity. Sea otter exploitation of productive butter clam habitat substantially reduced the abundance and altered the distribution of butter clams across Glacier Bay, with potential cascading consequences for nearshore community structure and function. Spatial variation in estimated sea otter predation processes further suggests that community context and local environmental conditions mediate the top-down influence of sea otters on a given prey. Overall, our framework provides high-resolution insights about the interaction among components of this food web and could be applied to a variety of other systems involving invasive species, epidemiology or migration.

Genomic analyses reveal range-wide devastation of sea otter populations.

The genetic consequences of species-wide declines are rarely quantified because the timing and extent of the decline varies across the species' range. The sea otter (Enhydra lutris) is a unique model in this regard. Their dramatic decline from thousands to fewer than 100 individuals per population occurred range-wide and nearly simultaneously due to the 18th-19th century fur trade. Consequently, each sea otter population represents an independent natural experiment of recovery after extreme population decline. We designed sequence capture probes for 50 Mb of sea otter exonic and neutral genomic regions. We sequenced 107 sea otters from five populations that span the species range to high coverage (18-76×) and three historical Californian samples from ~1500 and ~200 years ago to low coverage (1.5-3.5×). We observe distinct population structure and find that sea otters in California are the last survivors of a divergent lineage isolated for thousands of years and therefore warrant special conservation concern. We detect signals of extreme population decline in every surviving sea otter population and use this demographic history to design forward-in-time simulations of coding sequence. Our simulations indicate that this decline could lower the fitness of recovering populations for generations. However, the simulations also demonstrate how historically low effective population sizes prior to the fur trade may have mitigated the effects of population decline on genetic health. Our comprehensive approach shows how demographic inference from genomic data, coupled with simulations, allows assessment of extinction risk and different models of recovery.

Ecosystem response persists after a prolonged marine heatwave.

Some of the longest and most comprehensive marine ecosystem monitoring programs were established in the Gulf of Alaska following the environmental disaster of the Exxon Valdez oil spill over 30 years ago. These monitoring programs have been successful in assessing recovery from oil spill impacts, and their continuation decades later has now provided an unparalleled assessment of ecosystem responses to another newly emerging global threat, marine heatwaves. The 2014-2016 northeast Pacific marine heatwave (PMH) in the Gulf of Alaska was the longest lasting heatwave globally over the past decade, with some cooling, but also continued warm conditions through 2019. Our analysis of 187 time series from primary production to commercial fisheries and nearshore intertidal to offshore oceanic domains demonstrate abrupt changes across trophic levels, with many responses persisting up to at least 5 years after the onset of the heatwave. Furthermore, our suite of metrics showed novel community-level groupings relative to at least a decade prior to the heatwave. Given anticipated increases in marine heatwaves under current climate projections, it remains uncertain when or if the Gulf of Alaska ecosystem will return to a pre-PMH state.

Spatial epidemiological patterns suggest mechanisms of land-sea transmission for Sarcocystis neurona in a coastal marine mammal.

Sarcocystis neurona was recognised as an important cause of mortality in southern sea otters (Enhydra lutris nereis) after an outbreak in April 2004 and has since been detected in many marine mammal species in the Northeast Pacific Ocean. Risk of S. neurona exposure in sea otters is associated with consumption of clams and soft-sediment prey and is temporally associated with runoff events. We examined the spatial distribution of S. neurona exposure risk based on serum antibody testing and assessed risk factors for exposure in animals from California, Washington, British Columbia and Alaska. Significant spatial clustering of seropositive animals was observed in California and Washington, compared with British Columbia and Alaska. Adult males were at greatest risk for exposure to S. neurona, and there were strong associations with terrestrial features (wetlands, cropland, high human housing-unit density). In California, habitats containing soft sediment exhibited greater risk than hard substrate or kelp beds. Consuming a diet rich in clams was also associated with increased exposure risk. These findings suggest a transmission pathway analogous to that described for Toxoplasma gondii, with infectious stages traveling in freshwater runoff and being concentrated in particular locations by marine habitat features, ocean physical processes, and invertebrate bioconcentration.

Defining the risk landscape in the context of pathogen pollution: Toxoplasma gondii in sea otters along the Pacific Rim.

Pathogens entering the marine environment as pollutants exhibit a spatial signature driven by their transport mechanisms. The sea otter (Enhydra lutris), a marine animal which lives much of its life within sight of land, presents a unique opportunity to understand land-sea pathogen transmission. Using a dataset on Toxoplasma gondii prevalence across sea otter range from Alaska to California, we found that the dominant drivers of infection risk vary depending upon the spatial scale of analysis. At the population level, regions with high T. gondii prevalence had higher human population density and a greater proportion of human-dominated land uses, suggesting a strong role for population density of the felid definitive host of this parasite. This relationship persisted when a subset of data were analysed at the individual level: large-scale patterns in sea otter T. gondii infection prevalence were largely explained by individual exposure to areas of high human housing unit density, and other landscape features associated with anthropogenic land use, such as impervious surfaces and cropping land. These results contrast with the small-scale, within-region analysis, in which age, sex and prey choice accounted for most of the variation in infection risk, and terrestrial environmental features provided little variation to help in explaining observed patterns. These results underscore the importance of spatial scale in study design when quantifying both individual-level risk factors and landscape-scale variation in infection risk.

Size, growth, and density data for shallow-water sea urchins from Mexico to the Aleutian Islands, Alaska, 1956-2016.

Size, growth, and density have been studied for North American Pacific coast sea urchins Strongylocentrotus purpuratus, S. droebachiensis, S. polyacanthus, Mesocentrotus (Strongylocentrotus) franciscanus, Lytechinus pictus, Centrostephanus coronatus, and Arbacia stellata by various workers at diverse sites and for varying lengths of time from 1956 to present. Numerous peer-reviewed publications have used some of these data but some data have appeared only in graduate theses or the gray literature. There also are data that have never appeared outside original data sheets. Motivation for studies has included fisheries management and environmental monitoring of sewer and power plant outfalls as well as changes associated with disease epidemics. Studies also have focused on kelp restoration, community effects of sea otters, basic sea urchin biology, and monitoring. The data sets presented here are a historical record of size, density, and growth for a common group of marine invertebrates in intertidal and nearshore environments that can be used to test hypotheses concerning future changes associated with fisheries practices, shifts of predator distributions, climate and ecosystem changes, and ocean acidification along the Pacific Coast of North America and islands of the north Pacific. No copyright restrictions apply. Please credit this paper when using the data.

The interaction of intraspecific competition and habitat on individual diet specialization: a near range-wide examination of sea otters.

The quantification of individuality is a common research theme in the fields of population, community, and evolutionary ecology. The potential for individuality to arise is likely context-dependent, and the influence of habitat characteristics on its prevalence has received less attention than intraspecific competition. We examined individual diet specialization in 16 sea otter (Enhydra lutris) populations from southern California to the Aleutian Islands in Alaska. Because population histories, relative densities, and habitat characteristics vary widely among sites, we could examine the effects of intraspecific competition and habitat on the prevalence of individual diet specialization. Using observed diet data, we classified half of our sites as rocky substrate habitats and the other half containing a mixture of rocky and unconsolidated (soft) sediment substrates. We used stable isotope data to quantify population- and individual-level diet variation. Among rocky substrate sites, the slope [±standard error (SE)] of the positive significant relationship between the within-individual component (WIC) and total isotopic niche width (TINW) was shallow (0.23 ± 0.07) and negatively correlated with sea otter density. In contrast, the slope of the positive WIC/TINW relationship for populations inhabiting mixed substrate habitats was much higher (0.53 ± 0.14), suggesting a low degree of individuality, irrespective of intraspecific competition. Our results show that the potential for individuality to occur as a result of increasing intraspecific competition is context-dependent and that habitat characteristics, which ultimately influence prey diversity, relative abundance, and the range of skillsets required for efficient prey procurement, are important in determining when and where individual diet specialization occurs in nature.

Structure and mechanism of diet specialisation: testing models of individual variation in resource use with sea otters.

Studies of consumer-resource interactions suggest that individual diet specialisation is empirically widespread and theoretically important to the organisation and dynamics of populations and communities. We used weighted networks to analyze the resource use by sea otters, testing three alternative models for how individual diet specialisation may arise. As expected, individual specialisation was absent when otter density was low, but increased at high-otter density. A high-density emergence of nested resource-use networks was consistent with the model assuming individuals share preference ranks. However, a density-dependent emergence of a non-nested modular network for 'core' resources was more consistent with the 'competitive refuge' model. Individuals from different diet modules showed predictable variation in rank-order prey preferences and handling times of core resources, further supporting the competitive refuge model. Our findings support a hierarchical organisation of diet specialisation and suggest individual use of core and marginal resources may be driven by different selective pressures.

Gene transcription in sea otters (Enhydra lutris); development of a diagnostic tool for sea otter and ecosystem health.

Gene transcription analysis for diagnosing or monitoring wildlife health requires the ability to distinguish pathophysiological change from natural variation. Herein, we describe methodology for the development of quantitative real-time polymerase chain reaction (qPCR) assays to measure differential transcript levels of multiple immune function genes in the sea otter (Enhydra lutris); sea otter-specific qPCR primer sequences for the genes of interest are defined. We establish a 'reference' range of transcripts for each gene in a group of clinically healthy captive and free-ranging sea otters. The 10 genes of interest represent multiple physiological systems that play a role in immuno-modulation, inflammation, cell protection, tumour suppression, cellular stress response, xenobiotic metabolizing enzymes, antioxidant enzymes and cell-cell adhesion. The cycle threshold (C(T)) measures for most genes were normally distributed; the complement cytolysis inhibitor was the exception. The relative enumeration of multiple gene transcripts in simple peripheral blood samples expands the diagnostic capability currently available to assess the health of sea otters in situ and provides a better understanding of the state of their environment.

Assessment of clinical pathology and pathogen exposure in sea otters (Enhydra lutris) bordering the threatened population in Alaska.

Northern sea otter (Enhydra lutris kenyoni) abundance has decreased dramatically over portions of southwest Alaska, USA, since the mid-1980s, and this stock is currently listed as threatened under the Endangered Species Act. In contrast, adjacent populations in south central Alaska, USA, and Russia have been stable to increasing during the same period. Sea otters bordering the area classified in the recent decline were live-captured during 2004-2006 at Bering Island, Russia, and the Kodiak Archipelago, Alaska, USA, to evaluate differences in general health and current exposure status to marine and terrestrial pathogens. Although body condition was lower in animals captured at Bering Island, Russia, than it was at Kodiak, USA, clinical pathology values did not reveal differences in general health between the two regions. Low prevalences of antibodies (<5%) were found in Kodiak, USA, and on Bering Island, Russia, to Toxoplasma gondii, Sarcocystis neurona, and Leptospira interrogans. Exposure to phocine herpesvirus-1 was found in both Kodiak, USA (15.2%), and Bering Island, Russia (2.3%). Antibodies to Brucella spp. were found in 28% of the otters tested on Bering Island, Russia, compared with only 2.7% of the samples from Kodiak, USA. Prevalence of exposure to Phocine distemper virus (PDV) was 41% in Kodiak, USA, but 0% on Bering Island, Russia. Archived sera from southwest and south-central Alaska dating back to 1989 were negative for PDV, indicating exposure occurred in sea otters in Kodiak, USA, in recent years. Because PDV can be highly pathogenic in naïve and susceptible marine mammal populations, tissues should be examined to explore the contribution of this virus to otter deaths. Our results reveal an increase in exposure to pathogens in sea otters in Kodiak, Alaska, USA, since the 1990 s.

PCB exposure in sea otters and harlequin ducks in relation to history of contamination by the Exxon Valdez oil spill.

Exposure to contaminants other than petroleum hydrocarbons could confound interpretation of Exxon Valdez oil spill effects on biota at Prince William Sound, Alaska. Hence, we investigated polychlorinated biphenyls (PCBs) in blood of sea otters and harlequin ducks sampled during 1998. PCB concentrations characterized by lower chlorinated congeners were highest in sea otters from the unoiled area, whereas concentrations were similar among harlequin ducks from the oiled and unoiled area. Blood enzymes often elevated by xenobiotics were not related to PCB concentrations in sea otters. Only sea otters from the unoiled area had estimated risk from PCBs, and PCB composition or concentrations did not correspond to reported lower measures of population performance in sea otters or harlequin ducks from the oiled area. PCBs probably did not influence limited sea otter or harlequin duck recovery in the oiled area a decade after the spill.

Using ecological function to develop recovery criteria for depleted species: sea otters and kelp forests in the Aleutian archipelago.

Recovery criteria for depleted species or populations normally are based on demographic measures, the goal being to maintain enough individuals over a sufficiently large area to assure a socially tolerable risk of future extinction. Such demographically based recovery criteria may be insufficient to restore the functional roles of strongly interacting species. We explored the idea of developing a recovery criterion for sea otters (Enhydra lutris) in the Aleutian archipelago on the basis of their keystone role in kelp forest ecosystems. We surveyed sea otters and rocky reef habitats at 34 island-time combinations. The system nearly always existed in either a kelp-dominated or deforested phase state, which was predictable from sea otter density. We used a resampling analysis of these data to show that the phase state at any particular island can be determined at 95% probability of correct classification with information from as few as six sites. When sea otter population status (and thus the phase state of the kelp forest) was allowed to vary randomly among islands, just 15 islands had to be sampled to estimate the true proportion that were kelp dominated (within 10%) with 90% confidence. We conclude that kelp forest phase state is a more appropriate, sensitive, and cost-effective measure of sea otter recovery than the more traditional demographically based metrics, and we suggest that similar approaches have broad potential utility in establishing recovery criteria for depleted populations of other functionally important species.

Incorporating diverse data and realistic complexity into demographic estimation procedures for sea otters.

Reliable information on historical and current population dynamics is central to understanding patterns of growth and decline in animal populations. We developed a maximum likelihood-based analysis to estimate spatial and temporal trends in age/sex-specific survival rates for the threatened southern sea otter (Enhydra lutris nereis), using annual population censuses and the age structure of salvaged carcass collections. We evaluated a wide range of possible spatial and temporal effects and used model averaging to incorporate model uncertainty into the resulting estimates of key vital rates and their variances. We compared these results to current demographic parameters estimated in a telemetry-based study conducted between 2001 and 2004. These results show that survival has decreased substantially from the early 1990s to the present and is generally lowest in the north-central portion of the population's range. The greatest temporal decrease in survival was for adult females, and variation in the survival of this age/sex class is primarily responsible for regulating population growth and driving population trends. Our results can be used to focus future research on southern sea otters by highlighting the life history stages and mortality factors most relevant to conservation. More broadly, we have illustrated how the powerful and relatively straightforward tools of information-theoretic-based model fitting can be used to sort through and parameterize quite complex demographic modeling frameworks.

Long-term ecosystem response to the Exxon Valdez oil spill.

The ecosystem response to the 1989 spill of oil from the Exxon Valdez into Prince William Sound, Alaska, shows that current practices for assessing ecological risks of oil in the oceans and, by extension, other toxic sources should be changed. Previously, it was assumed that impacts to populations derive almost exclusively from acute mortality. However, in the Alaskan coastal ecosystem, unexpected persistence of toxic subsurface oil and chronic exposures, even at sublethal levels, have continued to affect wildlife. Delayed population reductions and cascades of indirect effects postponed recovery. Development of ecosystem-based toxicology is required to understand and ultimately predict chronic, delayed, and indirect long-term risks and impacts.