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.

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.

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.

Chaotic genetic structure and past demographic expansion of the invasive gastropod Tritia neritea in its native range, the Mediterranean Sea.

To better predict population evolution of invasive species in introduced areas it is critical to identify and understand the mechanisms driving genetic diversity and structure in their native range. Here, we combined analyses of the mitochondrial COI gene and 11 microsatellite markers to investigate both past demographic history and contemporaneous genetic structure in the native area of the gastropod Tritia neritea, using Bayesian skyline plots (BSP), multivariate analyses and Bayesian clustering. The BSP framework revealed population expansions, dated after the last glacial maximum. The haplotype network revealed a strong geographic clustering. Multivariate analyses and Bayesian clustering highlighted the strong genetic structure at all scales, between the Black Sea and the Adriatic Sea, but also within basins. Within basins, a random pattern of genetic patchiness was observed, suggesting a superimposition of processes involving natural biological effects (no larval phase and thus limited larval dispersal) and putative anthropogenic transport of specimens. Contrary to the introduced area, no isolation-by-distance patterns were recovered in the Mediterranean or the Black Seas, highlighting different mechanisms at play on both native and introduced areas, triggering unknown consequences for species' evolutionary trajectories. These results of Tritia neritea populations on its native range highlight a mixture of ancient and recent processes, with the effects of paleoclimates and life history traits likely tangled with the effects of human-mediated dispersal.

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.

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.

Population genomics meet Lagrangian simulations: Oceanographic patterns and long larval duration ensure connectivity among Paracentrotus lividus populations in the Adriatic and Ionian seas.

Connectivity between populations influences both their dynamics and the genetic structuring of species. In this study, we explored connectivity patterns of a marine species with long-distance dispersal, the edible common sea urchin Paracentrotus lividus, focusing mainly on the Adriatic-Ionian basins (Central Mediterranean). We applied a multidisciplinary approach integrating population genomics, based on 1,122 single nucleotide polymorphisms (SNPs) obtained from 2b-RAD in 275 samples, with Lagrangian simulations performed with a biophysical model of larval dispersal. We detected genetic homogeneity among eight population samples collected in the focal Adriatic-Ionian area, whereas weak but significant differentiation was found with respect to two samples from the Western Mediterranean (France and Tunisia). This result was not affected by the few putative outlier loci identified in our dataset. Lagrangian simulations found a significant potential for larval exchange among the eight Adriatic-Ionian locations, supporting the hypothesis of connectivity of P. lividus populations in this area. A peculiar pattern emerged from the comparison of our results with those obtained from published P. lividus cytochrome b (cytb) sequences, the latter revealing genetic differentiation in the same geographic area despite a smaller sample size and a lower power to detect differences. The comparison with studies conducted using nuclear markers on other species with similar pelagic larval durations in the same Adriatic-Ionian locations indicates species-specific differences in genetic connectivity patterns and warns against generalizing single-species results to the entire community of rocky shore habitats.

Population genomics of an endemic Mediterranean fish: differentiation by fine scale dispersal and adaptation.

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. For species with large effective population sizes a low number of markers may fail to identify population structure. A solution of this shortcoming can be high-throughput sequencing that allows genotyping thousands of markers on a genome-wide approach while facilitating the detection of genetic structuring shaped by selection. We used Genotyping-by-Sequencing (GBS) on 176 individuals of the endemic East Atlantic peacock wrasse (Symphodus tinca), from 6 locations in the Adriatic and Ionian seas. We obtained a total of 4,155 polymorphic SNPs and we observed two strong barriers to gene flow. The first one differentiated Tremiti Islands, in the northwest, from all the other locations while the second one separated east and south-west localities. Outlier SNPs potentially under positive selection and neutral SNPs both showed similar patterns of structuring, although finer scale differentiation was unveiled with outlier loci. Our results reflect the complexity of population genetic structure and demonstrate that both habitat fragmentation and positive selection are on play. This complexity should be considered in biodiversity assessments of different taxa, including non-model yet ecologically relevant organisms.

Error, Power, and Blind Sentinels: The Statistics of Seagrass Monitoring.

We derive statistical properties of standard methods for monitoring of habitat cover worldwide, and criticize them in the context of mandated seagrass monitoring programs, as exemplified by Posidonia oceanica in the Mediterranean Sea. We report the novel result that cartographic methods with non-trivial classification errors are generally incapable of reliably detecting habitat cover losses less than about 30 to 50%, and the field labor required to increase their precision can be orders of magnitude higher than that required to estimate habitat loss directly in a field campaign. We derive a universal utility threshold of classification error in habitat maps that represents the minimum habitat map accuracy above which direct methods are superior. Widespread government reliance on blind-sentinel methods for monitoring seafloor can obscure the gradual and currently ongoing losses of benthic resources until the time has long passed for meaningful management intervention. We find two classes of methods with very high statistical power for detecting small habitat cover losses: 1) fixed-plot direct methods, which are over 100 times as efficient as direct random-plot methods in a variable habitat mosaic; and 2) remote methods with very low classification error such as geospatial underwater videography, which is an emerging, low-cost, non-destructive method for documenting small changes at millimeter visual resolution. General adoption of these methods and their further development will require a fundamental cultural change in conservation and management bodies towards the recognition and promotion of requirements of minimal statistical power and precision in the development of international goals for monitoring these valuable resources and the ecological services they provide.