Within- and transgenerational stress legacy effects of ocean acidification on red abalone (Haliotis rufescens) growth and survival.

Understanding the mechanisms by which individual organisms respond and populations adapt to global climate change is a critical challenge. The role of plasticity and acclimation, within and across generations, may be essential given the pace of change. We investigated plasticity across generations and life stages in response to ocean acidification (OA), which poses a growing threat to both wild populations and the sustainable aquaculture of shellfish. Most studies of OA on shellfish focus on acute effects, and less is known regarding the longer term carryover effects that may manifest within or across generations. We assessed these longer term effects in red abalone (Haliotis rufescens) using a multi-generational split-brood experiment. We spawned adults raised in ambient conditions to create offspring that we then exposed to high pCO2 (1180 µatm; simulating OA) or low pCO2 (450 µatm; control or ambient conditions) during the first 3 months of life. We then allowed these animals to reach maturity in ambient common garden conditions for 4 years before returning the adults into high or low pCO2 treatments for 11 months and measuring growth and reproductive potential. Early-life exposure to OA in the F1 generation decreased adult growth rate even after 5 years especially when abalone were re-exposed to OA as adults. Adult but not early-life exposure to OA negatively impacted fecundity. We then exposed the F2 offspring to high or low pCO2 treatments for the first 3 months of life in a fully factorial, split-brood design. We found negative transgenerational effects of parental OA exposure on survival and growth of F2 offspring, in addition to significant direct effects of OA on F2 survival. These results show that the negative impacts of OA can last within and across generations, but that buffering against OA conditions at critical life-history windows can mitigate these effects.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina).

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Local environment drives rapid shifts in composition and phylogenetic clustering of seagrass microbiomes.

Plant microbiomes depend on environmental conditions, stochasticity, host species, and genotype identity. Eelgrass (Zostera marina) is a unique system for plant-microbe interactions as a marine angiosperm growing in a physiologically-challenging environment with anoxic sediment, periodic exposure to air at low tide, and fluctuations in water clarity and flow. We tested the influence of host origin versus environment on eelgrass microbiome composition by transplanting 768 plants among four sites within Bodega Harbor, CA. Over three months following transplantation, we sampled microbial communities monthly on leaves and roots and sequenced the V4-V5 region of the 16S rRNA gene to assess community composition. The main driver of leaf and root microbiome composition was destination site; more modest effects of host origin site did not last longer than one month. Community phylogenetic analyses suggested that environmental filtering structures these communities, but the strength and nature of this filtering varies among sites and over time and roots and leaves show opposing gradients in clustering along a temperature gradient. We demonstrate that local environmental differences create rapid shifts in associated microbial community composition with potential functional implications for rapid host acclimation under shifting environmental conditions.

Genomic responses to parallel temperature gradients in the eelgrass Zostera marina in adjacent bays.

The extent of parallel genomic responses to similar selective pressures depends on a complex array of environmental, demographic, and evolutionary forces. Laboratory experiments with replicated selective pressures yield mixed outcomes under controlled conditions and our understanding of genomic parallelism in the wild is limited to a few well-established systems. Here, we examine genomic signals of selection in the eelgrass Zostera marina across temperature gradients in adjacent embayments. Although we find many genomic regions with signals of selection within each bay there is very little overlap in signals of selection at the SNP level, despite most polymorphisms being shared across bays. We do find overlap at the gene level, potentially suggesting multiple mutational pathways to the same phenotype. Using polygenic models we find that some sets of candidate SNPs are able to predict temperature across both bays, suggesting that small but parallel shifts in allele frequencies may be missed by independent genome scans. Together, these results highlight the continuous rather than binary nature of parallel evolution in polygenic traits and the complexity of evolutionary predictability.

The effect of warming on seagrass wasting disease depends on host genotypic identity and diversity.

Temperature increases due to climate change have affected the distribution and severity of diseases in natural systems, causing outbreaks that can destroy host populations. Host identity, diversity, and the associated microbiome can affect host responses to both infection and temperature, but little is known about how they could function as important mediators of disease in altered thermal environments. We conducted an 8-week warming experiment to test the independent and interactive effects of warming, host genotypic identity, and host genotypic diversity on the prevalence and intensity of infections of seagrass (Zostera marina) by the wasting disease parasite (Labyrinthula zosterae). At elevated temperatures, we found that genotypically diverse host assemblages had reduced infection intensity, but not reduced prevalence, relative to less diverse assemblages. This dilution effect on parasite intensity was the result of both host composition effects as well as emergent properties of biodiversity. In contrast with the benefits of genotypic diversity under warming, diversity actually increased parasite intensity slightly in ambient temperatures. We found mixed support for the hypothesis that a growth-defense trade-off contributed to elevated disease intensity under warming. Changes in the abundance (but not composition) of a few taxa in the host microbiome were correlated with genotype-specific responses to wasting disease infections under warming, consistent with the emerging evidence linking changes in the host microbiome to the outcome of host-parasite interactions. This work emphasizes the context dependence of biodiversity-disease relationships and highlights the potential importance of interactions among biodiversity loss, climate change, and disease outbreaks in a key foundation species.

Detecting clonemate pairs in multicellular diploid clonal species based on a shared heterozygosity index.

Clonal reproduction, the formation of nearly identical individuals via mitosis in the absence of genetic recombination, is a very common reproductive mode across plants, fungi and animals. To detect clonal genetic structure, genetic similarity indices based on shared alleles are widely used, such as the Jaccard index, or identity by state. Here we propose a new pairwise genetic similarity index, the SH index, based on segregating genetic marker loci (typically single nucleotide polymorphisms) that are identically heterozygous for pairs of samples (NSH ). To test our method, we analyse two old seagrass clones (Posidonia australis, estimated to be around 8500 years old; Zostera marina, >750 years old) along with two young Z. marina clones of known age (17 years old). We show that focusing on shared heterozygosity amplifies the power to distinguish sample pairs belonging to different clones compared to methods focusing on all shared alleles. Our proposed workflow can successfully detect clonemates at a location dominated by a single clone. When the collected samples involve two or more clones, the SH index shows a clear gap between clonemate pairs and interclone sample pairs. Ideally NSH should be on the order of approximately ≥3000, a number easily achievable via restriction-site associated DNA (RAD) sequencing or whole-genome resequencing. Another potential application of the SH index is to detect possible parent-descendant pairs under selfing. Our proposed workflow takes advantage of the availability of the larger number of genetic markers in the genomic era, and improves the ability to distinguish clonemates from nonclonemates in multicellular diploid clonal species.

Local adaptation in the transgenerational response to copper pollution in the bryozoan Bugula neritina.

Transgenerational plasticity (TGP)-when a parent or previous generation's environmental experience affects offspring phenotype without involving a genetic change-can be an important mechanism allowing for rapid adaptation. However, despite increasing numbers of empirical examples of TGP, there appears to be considerable variation in its strength and direction, yet limited understanding of what causes this variation. We compared patterns of TGP in response to stress across two populations with high versus low historical levels of stress exposure. Specifically, we expected that exposure to acute stress in the population experiencing historically high levels of stress would result in adaptive TGP or alternatively fixed tolerance (no parental effect), whereas the population with low levels of historical exposure would result in negative parental carryover effects. Using a common sessile marine invertebrate, Bugula neritina, and a split brood design, we exposed parents from both populations to copper or control treatments in the laboratory and then had them brood copper-naïve larvae. We then exposed half of each larval brood to copper and half to control conditions before allowing them to grow to maturity in the field. Maternal copper exposure had a strong negative carryover effect on adult offspring growth and survival in the population without historical exposure, especially when larvae themselves were exposed to copper. We found little to no maternal or offspring treatment effect on adult growth and survival in the population with a history of copper exposure. However, parents from this population produced larger larvae on average and were able to increase the size of their larvae in response to copper exposure, providing a potential mechanism for maintaining fitness and suggesting TGP through maternal provisioning. These results indicate that the ability to adjust offspring phenotype via TGP may be a locally adapted trait and potentially influenced by past patterns of exposure.

A Pleistocene legacy structures variation in modern seagrass ecosystems.

Distribution of Earth's biomes is structured by the match between climate and plant traits, which in turn shape associated communities and ecosystem processes and services. However, that climate-trait match can be disrupted by historical events, with lasting ecosystem impacts. As Earth's environment changes faster than at any time in human history, critical questions are whether and how organismal traits and ecosystems can adjust to altered conditions. We quantified the relative importance of current environmental forcing versus evolutionary history in shaping the growth form (stature and biomass) and associated community of eelgrass (Zostera marina), a widespread foundation plant of marine ecosystems along Northern Hemisphere coastlines, which experienced major shifts in distribution and genetic composition during the Pleistocene. We found that eelgrass stature and biomass retain a legacy of the Pleistocene colonization of the Atlantic from the ancestral Pacific range and of more recent within-basin bottlenecks and genetic differentiation. This evolutionary legacy in turn influences the biomass of associated algae and invertebrates that fuel coastal food webs, with effects comparable to or stronger than effects of current environmental forcing. Such historical lags in phenotypic acclimatization may constrain ecosystem adjustments to rapid anthropogenic climate change, thus altering predictions about the future functioning of ecosystems.

Predictable Changes in Eelgrass Microbiomes with Increasing Wasting Disease Prevalence across 23° Latitude in the Northeastern Pacific.

Predicting outcomes of marine disease outbreaks presents a challenge in the face of both global and local stressors. Host-associated microbiomes may play important roles in disease dynamics but remain understudied in marine ecosystems. Host-pathogen-microbiome interactions can vary across host ranges, gradients of disease, and temperature; studying these relationships may aid our ability to forecast disease dynamics. Eelgrass, Zostera marina, is impacted by outbreaks of wasting disease caused by the opportunistic pathogen Labyrinthula zosterae. We investigated how Z. marina phyllosphere microbial communities vary with rising wasting disease lesion prevalence and severity relative to plant and meadow characteristics like shoot density, longest leaf length, and temperature across 23° latitude in the Northeastern Pacific. We detected effects of geography (11%) and smaller, but distinct, effects of temperature (30-day max sea surface temperature, 4%) and disease (lesion prevalence, 3%) on microbiome composition. Declines in alpha diversity on asymptomatic tissue occurred with rising wasting disease prevalence within meadows. However, no change in microbiome variability (dispersion) was detected between asymptomatic and symptomatic tissues. Further, we identified members of Cellvibrionaceae, Colwelliaceae, and Granulosicoccaceae on asymptomatic tissue that are predictive of wasting disease prevalence across the geographic range (3,100 kilometers). Functional roles of Colwelliaceae and Granulosicoccaceae are not known. Cellvibrionaceae, degraders of plant cellulose, were also enriched in lesions and adjacent green tissue relative to nonlesioned leaves. Cellvibrionaceae may play important roles in disease progression by degrading host tissues or overwhelming plant immune responses. Thus, inclusion of microbiomes in wasting disease studies may improve our ability to understand variable rates of infection, disease progression, and plant survival. IMPORTANCE The roles of marine microbiomes in disease remain poorly understood due, in part, to the challenging nature of sampling at appropriate spatiotemporal scales and across natural gradients of disease throughout host ranges. This is especially true for marine vascular plants like eelgrass (Zostera marina) that are vital for ecosystem function and biodiversity but are susceptible to rapid decline and die-off from pathogens like eukaryotic slime-mold Labyrinthula zosterae (wasting disease). We link bacterial members of phyllosphere tissues to the prevalence of wasting disease across the broadest geographic range to date for a marine plant microbiome-disease study (3,100 km). We identify Cellvibrionaceae, plant cell wall degraders, enriched (up to 61% relative abundance) within lesion tissue, which suggests this group may be playing important roles in disease progression. These findings suggest inclusion of microbiomes in marine disease studies will improve our ability to predict ecological outcomes of infection across variable landscapes spanning thousands of kilometers.

Disturbance decreases genotypic diversity by reducing colonization: Implications for disturbance-diversity feedbacks.

One objective of eco-evolutionary dynamics is to understand how the interplay between ecology and evolution on contemporary timescales contributes to the maintenance of biodiversity. Disturbance is an ecological process that can alter species diversity through both ecological and evolutionary effects on colonization and extinction dynamics. While analogous mechanisms likely operate among genotypes within a population, empirical evidence demonstrating the relationship between disturbance and genotypic diversity remains limited. We experimentally tested how disturbance altered the colonization (gain) and extinction (loss) of genets within a population of the marine angiosperm Zostera marina (eelgrass). In a 2-year field experiment conducted in northern California, we mimicked grazing disturbance by migratory geese by clipping leaves at varying frequencies during the winter months. Surprisingly, we found the greatest rates of new colonization in the absence of disturbance and that clipping had negligible effects on extinction. We hypothesize that genet extinction was not driven by selective mortality from clipping or from any stochastic loss resulting from the reduced shoot densities in clipped plots. We also hypothesize that increased flowering effort and facilitation within and among clones drove the increased colonization of new genets in the undisturbed treatment. This balance between colonization and extinction resulted in a negative relationship between clipping frequency and net changes in genotypic richness. We interpret our results in light of prior work showing that genotypic diversity increased resistance to grazing disturbance. We suggest that both directions of a feedback between disturbance and diversity occur in this system with consequences for the maintenance of eelgrass genotypic diversity.

The biogeography of community assembly: latitude and predation drive variation in community trait distribution in a guild of epifaunal crustaceans.

While considerable evidence exists of biogeographic patterns in the intensity of species interactions, the influence of these patterns on variation in community structure is less clear. Studying how the distributions of traits in communities vary along global gradients can inform how variation in interactions and other factors contribute to the process of community assembly. Using a model selection approach on measures of trait dispersion in crustaceans associated with eelgrass (Zostera marina) spanning 30° of latitude in two oceans, we found that dispersion strongly increased with increasing predation and decreasing latitude. Ocean and epiphyte load appeared as secondary predictors; Pacific communities were more overdispersed while Atlantic communities were more clustered, and increasing epiphytes were associated with increased clustering. By examining how species interactions and environmental filters influence community structure across biogeographic regions, we demonstrate how both latitudinal variation in species interactions and historical contingency shape these responses. Community trait distributions have implications for ecosystem stability and functioning, and integrating large-scale observations of environmental filters, species interactions and traits can help us predict how communities may respond to environmental change.

Sodium molybdate does not inhibit sulfate-reducing bacteria but increases shell growth in the Pacific oyster Magallana gigas.

Recent work on microbe-host interactions has revealed an important nexus between the environment, microbiome, and host fitness. Marine invertebrates that build carbonate skeletons are of particular interest in this regard because of predicted effects of ocean acidification on calcified organisms, and the potential of microbes to buffer these impacts. Here we investigate the role of sulfate-reducing bacteria, a group well known to affect carbonate chemistry, in Pacific oyster (Magallana gigas) shell formation. We reared oyster larvae to 51 days post fertilization and exposed organisms to control and sodium molybdate conditions, the latter of which is thought to inhibit bacterial sulfate reduction. Contrary to expectations, we found that sodium molybdate did not uniformly inhibit sulfate-reducing bacteria in oysters, and oysters exposed to molybdate grew larger shells over the experimental period. Additionally, we show that microbiome composition, host gene expression, and shell size were distinct between treatments earlier in ontogeny, but became more similar by the end of the experiment. Although additional testing is required to fully elucidate the mechanisms, our work provides preliminary evidence that M. gigas is capable of regulating microbiome dysbiosis caused by environmental perturbations, which is reflected in shell development.

Local adaptation in a marine foundation species: Implications for resilience to future global change.

Environmental change is multidimensional, with local anthropogenic stressors and global climate change interacting to differentially impact populations throughout a species' geographic range. Within species, the spatial distribution of phenotypic variation and its causes (i.e., local adaptation or plasticity) will determine species' adaptive capacity to respond to a changing environment. However, comparatively less is known about the spatial scale of adaptive differentiation among populations and how patterns of local adaptation might drive vulnerability to global change stressors. To test whether fine-scale (2-12 km) mosaics of environmental stress can cause adaptive differentiation in a marine foundation species, eelgrass (Zostera marina), we conducted a three-way reciprocal transplant experiment spanning the length of Tomales Bay, CA. Our results revealed strong home-site advantage in growth and survival for all three populations. In subsequent common garden experiments and feeding assays, we showed that countergradients in temperature, light availability, and grazing pressure from an introduced herbivore contribute to differential performance among populations consistent with local adaptation. Our findings highlight how local-scale mosaics in environmental stressors can increase phenotypic variation among neighboring populations, potentially increasing species resilience to future global change. More specifically, we identified a range-center eelgrass population that is pre-adapted to extremely warm temperatures similar to those experienced by low-latitude range-edge populations of eelgrass, demonstrating how reservoirs of heat-tolerant phenotypes may already exist throughout a species range. Future work on predicting species resilience to global change should incorporate potential buffering effects of local-scale population differentiation and promote a phenotypic management approach to species conservation.

Joint effects of patch edges and habitat degradation on faunal predation risk in a widespread marine foundation species.

Human activities degrade and fragment coastal marine habitats, reducing their structural complexity and making habitat edges a prevalent seascape feature. Though habitat edges frequently are implicated in reduced faunal survival and biodiversity, results of experiments on edge effects have been inconsistent, calling for a mechanistic approach to the study of edges that explicitly includes indirect and interactive effects of habitat alteration at multiple scales across biogeographic gradients. We used an experimental network spanning 17 eelgrass (Zostera marina) sites across the Atlantic and Pacific oceans and the Mediterranean Sea to determine (1) if eelgrass edges consistently increase faunal predation risk, (2) whether edge effects on predation risk are altered by habitat degradation (shoot thinning), and (3) whether variation in the strength of edge effects among sites can be explained by biogeographical variability in covarying eelgrass habitat features. Contrary to expectations, at most sites, predation risk for tethered crustaceans (crabs or shrimps) was lower along patch edges than in patch interiors, regardless of the extent of habitat degradation. However, the extent to which edges reduced predation risk, compared to the patch interior, was correlated with the extent to which edges supported higher eelgrass structural complexity and prey biomass compared to patch interiors. This suggests an indirect component to edge effects in which the impact of edge proximity on predation risk is mediated by the effect of edges on other key biotic factors. Our results suggest that studies on edge effects should consider structural characteristics of patch edges, which may vary geographically, and multiple ways that humans degrade habitats.

Climate drives the geography of marine consumption by changing predator communities.

The global distribution of primary production and consumption by humans (fisheries) is well-documented, but we have no map linking the central ecological process of consumption within food webs to temperature and other ecological drivers. Using standardized assays that span 105° of latitude on four continents, we show that rates of bait consumption by generalist predators in shallow marine ecosystems are tightly linked to both temperature and the composition of consumer assemblages. Unexpectedly, rates of consumption peaked at midlatitudes (25 to 35°) in both Northern and Southern Hemispheres across both seagrass and unvegetated sediment habitats. This pattern contrasts with terrestrial systems, where biotic interactions reportedly weaken away from the equator, but it parallels an emerging pattern of a subtropical peak in marine biodiversity. The higher consumption at midlatitudes was closely related to the type of consumers present, which explained rates of consumption better than consumer density, biomass, species diversity, or habitat. Indeed, the apparent effect of temperature on consumption was mostly driven by temperature-associated turnover in consumer community composition. Our findings reinforce the key influence of climate warming on altered species composition and highlight its implications for the functioning of Earth's ecosystems.

Previous exposure mediates the response of eelgrass to future warming via clonal transgenerational plasticity.

Mortality and shifts in species distributions are among the most obvious consequences of extreme climatic events. However, the sublethal effects of an extreme event can have persistent impacts throughout an individual's lifetime and into future generations via within-generation and transgenerational phenotypic plasticity. These changes can either confer resilience or increase susceptibility to subsequent stressful events, with impacts on population, community, and potentially ecosystem processes. Here, we show how a simulated extreme warming event causes persistent changes in the morphology and growth of a foundation species (eelgrass, Zostera marina) across multiple clonal generations and multiple years. The effect of previous parental exposure to warming increased aboveground biomass, shoot length, and aboveground-belowground biomass ratios while also greatly decreasing leaf growth rates. Long-term increases in aboveground-belowground biomass ratios could indicate an adaptive clonal transgenerational response to warmer climates that reduces the burden of increased respiration in belowground biomass. These transgenerational responses were likely decoupled from clonal parent provisioning as rhizome size of clonal offspring was standardized at planting and rhizome starch reserves were not impacted by warming treatments. Future investigations into potential epigenetic mechanisms underpinning such clonal transgenerational plasticity will be necessary to understand the resilience of asexual foundation species to repeated extreme climatic events.

Human-induced reductions in fish predator boldness decrease their predation rates in kelp forests.

Humans have restructured food webs and ecosystems by depleting biomass, reducing size structure and altering traits of consumers. However, few studies have examined the ecological impacts of human-induced trait changes across large spatial and temporal scales and species assemblages. We compared behavioural traits and predation rates by predatory fishes on standard squid prey in protected areas of different protection levels and ages, and found that predation rates were 6.5 times greater at old, no-take (greater than 40 years) relative to new, predominantly partial-take areas (approx. 8 years), even accounting for differences in predatory fish abundance, body size and composition across sites. Individual fishes in old protected areas consumed prey at nearly twice the rate of fishes of the same species and size at new protected areas. Predatory fish exhibited on average 50% longer flight initiation distance and lower willingness to forage at new protected areas, which partially explains lower foraging rates at new relative to old protected areas. Our experiments demonstrate that humans can effect changes in functionally important behavioural traits of predator guilds at large (30 km) spatial scales within managed areas, which require protection for multiple generations of predators to recover bold phenotypes and predation rates, even as abundance rebounds.