Coast-wide recruitment dynamics of Olympia oysters reveal limited synchrony and multiple predictors of failure

Recruitment of a threatened foundation oyster species varies with large and small spatial scales

Understanding how habitat attributes (e.g., patch area and sizes, connectivity) control recruitment and how this is modified by processes operating at larger spatial scales is fundamental to understanding population sustainability and developing successful long-term restoration strategies for marine foundation species-including for globally threatened reef-forming oysters. In two experiments, we assessed the recruitment and energy reserves of oyster recruits onto remnant reefs of the oyster Saccostrea glomerata in estuaries spanning 550 km of coastline in southeastern Australia. In the first experiment, we determined whether recruitment of oysters to settlement plates in three estuaries was correlated with reef attributes within patches (distances to patch edges and surface elevation), whole-patch attributes (shape and size of patches), and landscape attributes (connectivity). We also determined whether environmental factors (e.g., sedimentation and water temperature) explained the differences among recruitment plates. We also tested whether differences in energy reserves of recruits could explain the differences between two of the estuaries (one high- and one low-sedimentation estuary). In the second experiment, across six estuaries (three with nominally high and three with nominally low sedimentation rates), we tested the hypothesis that, at the estuary scale, recruitment and survival were negatively correlated to sedimentation. Overall, total oyster recruitment varied mostly at the scale of estuaries rather than with reef attributes and was negatively correlated with sedimentation. Percentage recruit survival was, however, similar among estuaries, although energy reserves and condition of recruits were lower at a high- compared to a low-sediment estuary. Within each estuary, total oyster recruitment increased with patch area and decreased with increasing tidal height. Our results showed that differences among estuaries have the largest influence on oyster recruitment and recruit health and this may be explained by environmental processes operating at the same scale. While survival was high across all estuaries, growth and reproduction of oysters on remnant reefs may be affected by sublethal effects on the health of recruits in high-sediment estuaries. Thus, restoration programs should consider lethal and sublethal effects of whole-estuary environmental processes when selecting sites and include environmental mitigation actions to maximize recruitment success.

Mapping oysters on the Pacific coast of North America: A coast-wide collaboration to inform enhanced conservation

To conserve coastal foundation species, it is essential to understand patterns of distribution and abundance and how they change over time. We synthesized oyster distribution data across the west coast of North America to develop conservation strategies for the native Olympia oyster (Ostrea lurida), and to characterize populations of the non-native Pacific oyster (Magallana gigas). We designed a user-friendly portal for data entry into ArcGIS Online and collected oyster records from unpublished data submitted by oyster experts and from the published literature. We used the resulting 2,000+ records to examine spatial and temporal patterns and made an interactive web-based map publicly available. Comparing records from pre-2000 vs. post-2000, we found that O. lurida significantly decreased in abundance and distribution, while M. gigas increased significantly. Currently the distribution and abundance of the two species are fairly similar, despite one species being endemic to this region since the Pleistocene, and the other a new introduction. We mapped the networks of sites occupied by oysters based on estimates of larval dispersal distance, and found that these networks were larger in Canada, Washington, and southern California than in other regions. We recommend restoration to enhance O. lurida, particularly within small networks, and to increase abundance where it declined. We also recommend restoring natural biogenic beds on mudflats and sandflats especially in the southern range, where native oysters are currently found most often on riprap and other anthropogenic structures. This project can serve as a model for collaborative mapping projects that inform conservation strategies for imperiled species or habitats.

Rapid enhancement of multiple ecosystem services following the restoration of a coastal foundation species

The global decline of marine foundation species (kelp forests, mangroves, salt marshes, and seagrasses) has contributed to the degradation of the coastal zone and threatens the loss of critical ecosystem services and functions. Restoration of marine foundation species has had variable success, especially for seagrasses, where a majority of restoration efforts have failed. While most seagrass restorations track structural attributes over time, rarely do restorations assess the suite of ecological functions that may be affected by restoration. Here we report on the results of two small-scale experimental seagrass restoration efforts in a central California estuary where we transplanted 117 0.25-m² plots (2,340 shoots) of the seagrass species Zostera marina. We quantified restoration success relative to persistent reference beds, and in comparison to unrestored, unvegetated areas. Within three years, our restored plots expanded ~8,500%, from a total initial area of 29 to 2,513 m² . The restored beds rapidly began to resemble the reference beds in (1) seagrass structural attributes (canopy height, shoot density, biomass), (2) ecological functions (macrofaunal species richness and abundance, epifaunal species richness, nursery function), and (3) biogeochemical functions (modulation of water quality). We also developed a multifunctionality index to assess cumulative functional performance, which revealed restored plots are intermediate between reference and unvegetated habitats, illustrating how rapidly multiple functions recovered over a short time period. Our comprehensive study is one of few published studies to quantify how seagrass restoration can enhance both biological and biogeochemical functions. Our study serves as a model for quantifying ecosystem services associated with the restoration of a foundation species and demonstrates the potential for rapid functional recovery that can be achieved through targeted restoration of fast-growing foundation species under suitable conditions.

Conservation aquaculture as a tool for imperiled marine species: Evaluation of opportunities and risks for Olympia oysters, Ostrea lurida

Conservation aquaculture is becoming an important tool to support the recovery of declining marine species and meet human needs. However, this tool comes with risks as well as rewards, which must be assessed to guide aquaculture activities and recovery efforts. Olympia oysters (Ostrea lurida) provide key ecosystem functions and services along the west coast of North America, but populations have declined to the point of local extinction in some estuaries. Here, we present a species-level, range-wide approach to strategically planning the use of aquaculture to promote recovery of Olympia oysters. We identified 12 benefits of culturing Olympia oysters, including identifying climate-resilient phenotypes that add diversity to growers’ portfolios. We also identified 11 key risks, including potential negative ecological and genetic consequences associated with the transfer of hatchery-raised oysters into wild populations. Informed by these trade-offs, we identified ten priority estuaries where aquaculture is most likely to benefit Olympia oyster recovery. The two highest scoring estuaries have isolated populations with extreme recruitment limitation—issues that can be addressed via aquaculture if hatchery capacity is expanded in priority areas. By integrating social criteria, we evaluated which project types would likely meet the goals of local stakeholders in each estuary. Community restoration was most broadly suited to the priority areas, with limited commercial aquaculture and no current community harvest of the species, although this is a future stakeholder goal. The framework we developed to evaluate aquaculture as a tool to support species recovery is transferable to other systems and species globally; we provide a guide to prioritizing local knowledge and developing recommendations for implementation by using transparent criteria. Our collaborative process engaging diverse stakeholders including managers, scientists, Indigenous Tribal representatives, and shellfish growers can be used elsewhere to seek win-win opportunities to expand conservation aquaculture where benefits are maximized for both people and imperiled species.

Temperature and salinity, not acidification, predict near-future larval growth and larval habitat suitability of Olympia oysters in the Salish Sea

Most invertebrates in the ocean begin their lives with planktonic larval phases that are critical for dispersal and distribution of these species. Larvae are particularly vulnerable to environmental change, so understanding interactive effects of environmental stressors on larval life is essential in predicting population persistence and vulnerability of species. Here, we use a novel experimental approach to rear larvae under interacting gradients of temperature, salinity, and ocean acidification, then model growth rate and duration of Olympia oyster larvae and predict the suitability of habitats for larval survival. We find that temperature and salinity are closely linked to larval growth and larval habitat suitability, but larvae are tolerant to acidification at this scale. We discover that present conditions in the Salish Sea are actually suboptimal for Olympia oyster larvae from populations in the region, and that larvae from these populations might actually benefit from some degree of global ocean change. Our models predict a vast decrease in mean pelagic larval duration by the year 2095, which has the potential to alter population dynamics for this species in future oceans. Additionally, we find that larval tolerance can explain large-scale biogeographic patterns for this species across its range.

Carryover effects of temperature and pCO2 across multiple Olympia oyster populations

Predicting how populations will respond to ocean change across generations is critical to effective conservation of marine species. One emerging factor is the influence of parental exposures on offspring phenotype, known as intergenerational carryover effects. Parental exposure may deliver beneficial or detrimental characteristics to offspring that can influence larval recruitment patterns, thus shaping how populations and community structure respond to ocean change. Impacts of adult exposure to elevated winter temperature and pCO₂ on reproduction and offspring viability were examined in the Olympia oyster (Ostrea lurida) using three populations of adult, hatchery-reared O. lurida, plus an additional cohort spawned from one of the populations. Oysters were sequentially exposed to elevated temperature (+4°C, at 10°C), followed by elevated pCO₂ (+2,204 μatm, at 3,045 μatm) during winter months. Male gametes were more developed after elevated temperature exposure and less developed after high pCO₂ exposure, but there was no impact on female gametes or sex ratios. Oysters previously exposed to elevated winter temperature released larvae earlier, regardless of pCO₂ exposure. Those exposed to elevated winter temperature as a sole treatment released more larvae on a daily basis but, when also exposed to high pCO₂ , there was no effect. These combined results indicate that elevated winter temperature accelerates O. lurida spermatogenesis, resulting in earlier larval release and increased production, with elevated pCO₂ exposure negating effects of elevated temperature. Altered recruitment patterns may therefore follow warmer winters due to precocious spawning, but these effects may be masked by coincidental high pCO₂ . Offspring were reared in common conditions for 1 yr, then deployed for 3 months in four estuarine bays with distinct environmental conditions. Offspring of parents exposed to elevated pCO₂ had higher survival rates in two of the four bays. This carryover effect demonstrates that parental conditions can have substantial ecologically relevant impacts that should be considered when predicting impacts of environmental change. Furthermore, Olympia oysters may be more resilient in certain environments when progenitors are pre-conditioned in stressful conditions. Combined with other recent studies, our work suggests that the Olympia may be more equipped than other oysters for the challenge of a changing ocean.

Population structure, genetic connectivity, and adaptation in the Olympia oyster (Ostrea lurida) along the west coast of North America

Effective management of threatened and exploited species requires an understanding of both the genetic connectivity among populations and local adaptation. The Olympia oyster (Ostrea lurida), patchily distributed from Baja California to the central coast of Canada, has a long history of population declines due to anthropogenic stressors. For such coastal marine species, population structure could follow a continuous isolation-by-distance model, contain regional blocks of genetic similarity separated by barriers to gene flow, or be consistent with a null model of no population structure. To distinguish between these hypotheses in O. lurida, 13,424 single nucleotide polymorphisms (SNPs) were used to characterize rangewide population structure, genetic connectivity, and adaptive divergence. Samples were collected across the species range on the west coast of North America, from southern California to Vancouver Island. A conservative approach for detecting putative loci under selection identified 235 SNPs across 129 GBS loci, which were functionally annotated and analyzed separately from the remaining neutral loci. While strong population structure was observed on a regional scale in both neutral and outlier markers, neutral markers had greater power to detect fine-scale structure. Geographic regions of reduced gene flow aligned with known marine biogeographic barriers, such as Cape Mendocino, Monterey Bay, and the currents around Cape Flattery. The outlier loci identified as under putative selection included genes involved in developmental regulation, sensory information processing, energy metabolism, immune response, and muscle contraction. These loci are excellent candidates for future research and may provide targets for genetic monitoring programs. Beyond specific applications for restoration and management of the Olympia oyster, this study lends to the growing body of evidence for both population structure and adaptive differentiation across a range of marine species exhibiting the potential for panmixia. Computational notebooks are available to facilitate reproducibility and future open-sourced research on the population structure of O. lurida.