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.

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.

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.

3D genomics across the tree of life reveals condensin II as a determinant of architecture type.

We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional (3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedly during eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with the absence of condensin II subunits. Moreover, condensin II depletion converts the architecture of the human genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state, centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physical model in which lengthwise compaction of chromosomes by condensin II during mitosis determines chromosome-scale genome architecture, with effects that are retained during the subsequent interphase. This mechanism likely has been conserved since the last common ancestor of all eukaryotes.

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.

The skin microbiome of elasmobranchs follows phylosymbiosis, but in teleost fishes, the microbiomes converge.

BACKGROUND: The vertebrate clade diverged into Chondrichthyes (sharks, rays, and chimeras) and Osteichthyes fishes (bony fishes) approximately 420 mya, with each group accumulating vast anatomical and physiological differences, including skin properties. The skin of Chondrichthyes fishes is covered in dermal denticles, whereas Osteichthyes fishes are covered in scales and are mucous rich. The divergence time among these two fish groups is hypothesized to result in predictable variation among symbionts. Here, using shotgun metagenomics, we test if patterns of diversity in the skin surface microbiome across the two fish clades match predictions made by phylosymbiosis theory. We hypothesize (1) the skin microbiome will be host and clade-specific, (2) evolutionary difference in elasmobranch and teleost will correspond with a concomitant increase in host-microbiome dissimilarity, and (3) the skin structure of the two groups will affect the taxonomic and functional composition of the microbiomes. RESULTS: We show that the taxonomic and functional composition of the microbiomes is host-specific. Teleost fish had lower average microbiome within clade similarity compared to among clade comparison, but their composition is not different among clade in a null based model. Elasmobranch's average similarity within clade was not different than across clade and not different in a null based model of comparison. In the comparison of host distance with microbiome distance, we found that the taxonomic composition of the microbiome was related to host distance for the elasmobranchs, but not the teleost fishes. In comparison, the gene function composition was not related to the host-organism distance for elasmobranchs but was negatively correlated with host distance for teleost fishes. CONCLUSION: Our results show the patterns of phylosymbiosis are not consistent across both fish clades, with the elasmobranchs showing phylosymbiosis, while the teleost fish are not. The discrepancy may be linked to alternative processes underpinning microbiome assemblage, including possible historical host-microbiome evolution of the elasmobranchs and convergent evolution in the teleost which filter specific microbial groups. Our comparison of the microbiomes among fishes represents an investigation into the microbial relationships of the oldest divergence of extant vertebrate hosts and reveals that microbial relationships are not consistent across evolutionary timescales. Video abstract.

Predator type influences the frequency of functional responses to prey in marine habitats.

The functional response of a consumer to a gradient of resource density is a widespread and consistent framework used to quantify the importance of consumption to population dynamics and stability. Within benthic marine ecosystems, both crustaceans and fishes can provide strong top-down pressure on prey populations. Taxon-specific differences in biomechanics or habitat use, among other factors, may lead to variable functional response forms or parameter values (attack rate, handling time). Based on a review of 189 individual functional response fits, we find that these predator guilds differ in their frequency distribution of functional response types, with crustaceans exhibiting nearly double the proportion of sigmoidal, density-dependent functional responses (Holling type III) as predatory fishes. The implications of this finding for prey population stability are significant because type III responses allow prey persistence while type II responses are de-stabilizing and can lead to extinction. Comparing per capita predation rates across diverse taxa can provide integrative insights into predatory effects and the ability of predation to drive community structure.

Experiments reveal limited top-down control of key herbivores in southern California kelp forests.

Predator responses to gradients in prey density have important implications for population regulation and are a potential structuring force for subtidal marine communities, particularly on rocky reefs where herbivorous sea urchins can drive community state shifts. On rocky reefs in southern California where predatory sea otters have been extirpated, top-down control of sea urchins by alternative predators has been hypothesized but rarely tested experimentally. In laboratory feeding assays, predatory spiny lobsters (Panulirus interruptus) demonstrated a saturating functional response to urchin prey, whereby urchin proportional mortality was inversely density-dependent. In field experiments on rocky reefs near San Diego, California, predators (primarily the labrid fish California sheephead, Semicossyphus pulcher) inflicted highly variable mortality on purple urchin (Strongylocentrotus purpuratus) prey across all density levels. However, at low to moderate densities commonly observed within kelp forests, purple urchin mortality increased to a peak at a density of ~11 urchins/m2 . Above that level, at densities typical of urchin barrens, purple urchin mortality was density-independent. When larger red urchins (Mesocentrotus franciscanus) were offered to predators simultaneously with purple urchins, mortality was density-independent. Underwater videography revealed a positive relationship between purple urchin density and both the number and richness of fish predators, but these correlations were not observed when red urchins were present. Our results demonstrate highly variable mortality rates across prey densities in this system and suggest that top-down control of urchins can occur only under limited circumstances. Our findings provide insight into the dynamics of alternate community states observed on rocky reefs.

Eelgrass structural complexity mediates mesograzer herbivory on epiphytic algae.

Structural complexity mediates ecological processes such as predation, competition, and recruitment in marine systems, but relatively little is known about its effects on herbivory. In temperate seagrasses, such as eelgrass (Zostera marina), the primary herbivores are small crustacean and gastropod mesograzers that promote seagrass persistence by preferentially consuming competing epiphytic algae. We used a laboratory grazing experiment, a field colonization experiment, and stable isotope analysis to determine whether one component of eelgrass structural complexity, shoot density, dictates the strength of mesograzer top-down effects on epiphytic algae, and whether this is influenced by mesograzer community composition. Our results suggest that increasing structural complexity shifted eelgrass communities from a bottom-up to a top-down controlled system. In the lab, mesograzers reduced epiphyte standing stock only in high-shoot density experimental communities, though grazing impact varied among different combinations of dominant mesograzer taxa. In our field experiment, epiphyte biomass was inversely correlated with mesograzer density in high but not in low-shoot density eelgrass plots. High-shoot density plots contained lower epiphyte biomass despite housing lower densities of mesograzers, when compared to low-density plots, suggesting potential effects of mesograzer behavior, community composition, or self-shading on epiphyte growth. Our results suggest that structural complexity can strongly influence rates of top-down and bottom-up processes in eelgrass habitat, and should be incorporated into future experiments on the role of herbivores in seagrass ecosystems.

Latitude, temperature, and habitat complexity predict predation pressure in eelgrass beds across the Northern Hemisphere.

Latitudinal gradients in species interactions are widely cited as potential causes or consequences of global patterns of biodiversity. However, mechanistic studies documenting changes in interactions across broad geographic ranges are limited. We surveyed predation intensity on common prey (live amphipods and gastropods) in communities of eelgrass (Zostera marina) at 48 sites across its Northern Hemisphere range, encompassing over 37° of latitude and four continental coastlines. Predation on amphipods declined with latitude on all coasts but declined more strongly along western ocean margins where temperature gradients are steeper. Whereas in situ water temperature at the time of the experiments was uncorrelated with predation, mean annual temperature strongly positively predicted predation, suggesting a more complex mechanism than simply increased metabolic activity at the time of predation. This large-scale biogeographic pattern was modified by local habitat characteristics; predation declined with higher shoot density both among and within sites. Predation rates on gastropods, by contrast, were uniformly low and varied little among sites. The high replication and geographic extent of our study not only provides additional evidence to support biogeographic variation in predation intensity, but also insight into the mechanisms that relate temperature and biogeographic gradients in species interactions.

Structural complexity and fish body size interactively affect habitat optimality.

Predator-prey interactions are strongly influenced by habitat structure, particularly in coastal marine habitats such as seagrasses in which structural complexity (SC) may vary over small spatial scales. For seagrass mesopredators such as juvenile fishes, optimality models predict that fitness will be maximized at levels of SC that enhance foraging but minimize predation risk, both of which are functions of body size. We tested the hypothesis that in eelgrass (Zostera marina) habitat, optimal SC for juvenile giant kelpfish (Heterostichus rostratus), an abundant eelgrass mesopredator in southern California, changes through ontogeny. To do this, we quantified eelgrass SC effects on habitat associations, relative predation risk, and foraging efficiency for three size classes of juvenile giant kelpfish. We found that habitat selection differed with fish size: small fish selected dense eelgrass, whereas larger fish selected sparse eelgrass. Small kelpfish experienced the lowest relative predation risk in dense eelgrass but also had higher foraging efficiency in dense eelgrass, suggesting that dense eelgrass is selected by these fish because it minimizes risk and maximizes potential for growth. Surprisingly, larger kelpfish did not experience lower predation risk than small kelpfish. However, larger kelpfish experienced higher foraging efficiency in sparse eelgrass vs. dense eelgrass, suggesting that they select sparse eelgrass to maximize foraging efficiency. Our study highlights that trade-offs between predation risk and foraging can occur within a single habitat type, that studies should consider how habitat value changes through ontogeny, and that seagrass habitat value may be maximal when within-patch variability in SC is high.

Interactive effects of predator and prey harvest on ecological resilience of rocky reefs.

A major goal of ecosystem-based fisheries management is to prevent fishery-induced shifts in community states. This requires an understanding of ecological resilience: the ability of an ecosystem to return to the same state following a perturbation, which can strongly depend on species interactions across trophic levels. We use a structured model of a temperate rocky reef to explore how multi-trophic level fisheries impact ecological resilience. Increasing fishing mortality of prey (urchins) has a minor effect on equilibrium biomass of kelp, urchins, and spiny lobster predators, but increases resilience by reducing the range of predator harvest rates at which alternative stable states are possible. Size-structured predation on urchins acts as the feedback maintaining each state. Our results demonstrate that the resilience of ecosystems strongly depends on the interactive effects of predator and prey harvest in multi-trophic level fisheries, which are common in marine ecosystems but are unaccounted for by traditional management.

Invasive prey indirectly increase predation on their native competitors.

Ecological theory predicts that invasive prey can interact with native prey directly by competing for shared resources or indirectly by changing the abundance or behavior of shared native predators. However, both the study and management of invasive prey have historically overlooked indirect effects. In southern California estuaries, introduction of the Asian nest mussel Arcuatula senhousia has been linked to profound changes in native bivalve assemblages, but the mechanisms of these interactions remain unclear. We performed three field experiments to assess the mechanisms of competition between Arcuatula and native bivalves, and evaluated the potential for Arcuatula to indirectly mediate native predator-prey dynamics. We found that Arcuatula reduces the diversity, abundance, and size of native bivalve recruits by preemptively exploiting space in surface sediments. When paired with native shallow-dwelling clams (Chione undatella and Laevicardium substriatum), Arcuatula reduces adult survival through overgrowth competition. However, Arcuatula also attracts native predators, causing apparent competition by indirectly increasing predation of native clams, especially for poorly defended species. Therefore, invasive prey can indirectly increase predation rates on native competitors by changing the behavior of shared native predators, but the magnitude of apparent competition strongly depends on the vulnerability of natives to predation. Interestingly, our results indicate that the vulnerability of invasive prey to predation can greatly exacerbate impacts on their native competitors. Our findings suggest that consideration of both direct and indirect effects of invasive prey, as well as native predator-prey relationships, should lead to more effective invasive species management.

Biodiversity mediates top-down control in eelgrass ecosystems: a global comparative-experimental approach.

Nutrient pollution and reduced grazing each can stimulate algal blooms as shown by numerous experiments. But because experiments rarely incorporate natural variation in environmental factors and biodiversity, conditions determining the relative strength of bottom-up and top-down forcing remain unresolved. We factorially added nutrients and reduced grazing at 15 sites across the range of the marine foundation species eelgrass (Zostera marina) to quantify how top-down and bottom-up control interact with natural gradients in biodiversity and environmental forcing. Experiments confirmed modest top-down control of algae, whereas fertilisation had no general effect. Unexpectedly, grazer and algal biomass were better predicted by cross-site variation in grazer and eelgrass diversity than by global environmental gradients. Moreover, these large-scale patterns corresponded strikingly with prior small-scale experiments. Our results link global and local evidence that biodiversity and top-down control strongly influence functioning of threatened seagrass ecosystems, and suggest that biodiversity is comparably important to global change stressors.

Disturbance facilitates the coexistence of antagonistic ecosystem engineers in California estuaries.

Ecological theory predicts that interactions between antagonistic ecosystem engineers can lead to local competitive exclusion, but disturbance can facilitate broader coexistence. However, few empirical studies have tested the potential for disturbance to mediate competition between engineers. We examined the capacity for disturbance and habitat modification to explain the disjunct distributions of two benthic ecosystem engineers, eelgrass Zostera marina and the burrowing ghost shrimp Neotrypaea californiensis, in two California estuaries. Sediment sampling in eelgrass and ghost shrimp patches revealed that ghost shrimp change benthic biogeochemistry over small scales (centimeters) but not patch scales (meters to tens of meters), suggesting a limited capacity for sediment modification to explain species distributions. To determine the relative competitive abilities of engineers, we conducted reciprocal transplantations of ghost shrimp and eelgrass. Local ghost shrimp densities declined rapidly following the addition of eelgrass, and transplanted eelgrass expanded laterally into the surrounding ghost shrimp-dominated areas. When transplanted into eelgrass patches, ghost shrimp failed to persist. Ghost shrimp were also displaced from plots with structural mimics of eelgrass rhizomes and roots, suggesting that autogenic habitat modification by eelgrass is an important mechanism determining ghost shrimp distributions. However, ghost shrimp were able to rapidly colonize experimental disturbances to eelgrass patch edges, which are common in shallow estuaries. We conclude that coexistence in this system is maintained by spatiotemporally asynchronous disturbances and a competition-colonization trade-off: eelgrass is a competitively superior ecosystem engineer, but benthic disturbances permit the coexistence of ghost shrimp at the landscape scale by modulating the availability of space.

Effects of habitat patchiness on American lobster movement across a gradient of predation risk and shelter competition.

The influence of landscape structure on marine ecological processes is receiving increasing attention. However, few studies conducted in coastal marine habitats have evaluated whether the effects of landscape structure on species interactions and organismal behavior are consistent across the range of an organism, over which landscape context and the strength of species interactions typically vary. American lobster (Homarus americanus) juveniles seek refuge from predators within shallow rocky habitat but make short-distance movements to forage outside of shelter. We evaluated how the patchiness of cobble habitat influences juvenile lobster movement by conducting mark-recapture experiments on lobsters placed within patchy and contiguous cobble plots in three regions of New England among which risk of predation and intraspecific shelter competition vary (Rhode Island, mid-coast Maine, and eastern Maine, USA). We also evaluated whether habitat patchiness influenced lobster colonization of plots and whether lobster fidelity to individual shelters corresponds to variability in predator abundance and conspecific density among regions. Cobble patchiness reduced rates of lobster movement in all three regions in 2004 and in two of three regions in 2005, despite large differences in landscape context among regions. Region had much larger effects on lobster colonization than did patchiness, but patchy plots were colonized at higher rates than were contiguous plots where lobster densities were highest. Fidelity to shelter was higher in regions with low conspecific density (Rhode Island and eastern Maine) than in mid-coast Maine where conspecific density is high and where unmarked lobsters often occupied shelters vacated by marked lobsters. Our results indicate that cobble patchiness influences juvenile lobster movement at small scales, but that the effects of patchiness on movement were consistent across much of the range of the American lobster despite strong regional variation in predator abundance and conspecific density.

Biotic resistance to invasion along an estuarine gradient.

Biotic resistance is the ability of native communities to repel the establishment of invasive species. Predation by native species may confer biotic resistance to communities, but the environmental context under which this form of biotic resistance occurs is not well understood. We evaluated several factors that influence the distribution of invasive Asian mussels (Musculista senhousia) in Mission Bay, a southern California estuary containing an extensive eelgrass (Zostera marina) habitat. Asian mussels exhibit a distinct spatial pattern of invasion, with extremely high densities towards the back of Mission Bay (up to 4,000 m(-2)) in contrast with near-complete absence at sites towards the front of the bay. We established that recruits arrived at sites where adult mussels were absent and found that dense eelgrass does not appear to preclude Asian mussel growth and survival. Mussel survival and growth were high in predator-exclusion plots throughout the bay, but mussel survival was low in the front of the bay when plots were open to predators. Additional experiments revealed that consumption by spiny lobsters (Panulirus interruptus) and a gastropod (Pteropurpura festiva) likely are the primary factors responsible for resistance to Asian mussel invasion. However, biotic resistance was dependent on location within the estuary (for both species) and also on the availability of a hard substratum (for P. festiva). Our findings indicate that biotic resistance in the form of predation may be conferred by higher order predators, but that the strength of resistance may strongly vary across estuarine gradients and depend on the nature of the locally available habitat.