Applying landscape metrics to species distribution model predictions to characterize internal range structure and associated changes.

Distributional shifts in species ranges provide critical evidence of ecological responses to climate change. Assessments of climate-driven changes typically focus on broad-scale range shifts (e.g. poleward or upward), with ecological consequences at regional and local scales commonly overlooked. While these changes are informative for species presenting continuous geographic ranges, many species have discontinuous distributions-both natural (e.g. mountain or coastal species) or human-induced (e.g. species inhabiting fragmented landscapes)-where within-range changes can be significant. Here, we use an ecosystem engineer species (Sabellaria alveolata) with a naturally fragmented distribution as a case study to assess climate-driven changes in within-range occupancy across its entire global distribution. To this end, we applied landscape ecology metrics to outputs from species distribution modelling (SDM) in a novel unified framework. SDM predicted a 27.5% overall increase in the area of potentially suitable habitat under RCP 4.5 by 2050, which taken in isolation would have led to the classification of the species as a climate change winner. SDM further revealed that the latitudinal range is predicted to shrink because of decreased habitat suitability in the equatorward part of the range, not compensated by a poleward expansion. The use of landscape ecology metrics provided additional insights by identifying regions that are predicted to become increasingly fragmented in the future, potentially increasing extirpation risk by jeopardising metapopulation dynamics. This increased range fragmentation could have dramatic consequences for ecosystem structure and functioning. Importantly, the proposed framework-which brings together SDM and landscape metrics-can be widely used to study currently overlooked climate-driven changes in species internal range structure, without requiring detailed empirical knowledge of the modelled species. This approach represents an important advancement beyond predictive envelope approaches and could reveal itself as paramount for managers whose spatial scale of action usually ranges from local to regional.

Coastal upwelling generates cryptic temperature refugia.

Understanding the effects of climate-mediated environmental variation on the distribution of organisms is critically important in an era of global change. We used wavelet analysis to quantify the spatiotemporal (co)variation in daily water temperature for predicting the distribution of cryptic refugia across 16 intertidal sites that were characterized as 'no', 'weak' or 'strong' upwelling and spanned 2000 km of the European Atlantic Coast. Sites experiencing weak upwelling exhibited high synchrony in temperature but low levels of co-variability at monthly to weekly timescales, whereas the opposite was true for sites experiencing strong upwelling. This suggests upwelling generates temporal thermal refugia that can promote organismal performance by both supplying colder water that mitigates thermal stress during hot Summer months and ensuring high levels of fine-scale variation in temperature that reduce the duration of thermal extremes. Additionally, pairwise correlograms based on the Pearson-product moment correlation coefficient and wavelet coherence revealed scale dependent trends in temperature fluctuations across space, with a rapid decay in strong upwelling sites at monthly and weekly timescales. This suggests upwelling also generates spatial thermal refugia that can 'rescue' populations from unfavorable conditions at local and regional scales. Overall, this study highlights the importance of identifying cryptic spatiotemporal refugia that emerge from fine-scale environmental variation to map potential patterns of organismal performance in a rapidly changing world.

Fine-scale abundance of rocky shore macroalgae species with distribution limits in NW Iberia in 2020/2021.

Background: Climate change has been increasing at an unprecedented rate in the last decades. Global warming has been causing a variety of impacts in marine ecosystems, including shifts in the geographical ranges of species. The north-western Iberian Peninsula coast is particularly interesting to study distribution shifts as it features a strong latitude thermal gradient, establishing a biogeographical transitional region where several cold- and warm-adapted species have their equatorward or poleward distributions. In the early 2000s, it appeared that, while warm-water species were already responding to warming, cold-water species did not display a coherent response. It is now necessary to gather up-to-date data on the distribution of the same group of species to understand if current patterns of change confirm or deny those observed back then, which may give us important clues about the mechanisms setting species limits in the area. New information: This study provides a fine-scale description of the occurrence of intertidal macroalgae species in the rocky shores of the north-western Iberian coast. Specifically, the spatial distribution and semi-quantitative abundance of 34 native and invasive species were assessed at 70 wave-exposed locations. This included 19 species of cold-water affinity, 10 species of warm-water affinity and five neutral species. When contrasted with historical observations, these new data can be used to quantify and map biodiversity change in the region, as well as help understanding the mechanisms constraining species distributions.

Transcriptomic response of the intertidal limpet Patella vulgata to temperature extremes.

Global warming is challenging wild species in land and water. In the intertidal zone, species are already living at their thermal limits, being vulnerable even to small increases in maximum habitat temperatures. Knowledge of the mechanisms by which many intertidal zone species cope with elevated temperatures is limited. We analysed the molecular thermal stress response of the limpet Patella vulgata under slight and frequent (one-day), and extreme and rare (three-day) warming events. Using RNA-seq to assess differential gene expression among treatments, differing molecular responses were obtained in the two treatments, with more changes in gene expression after the three-day event; with one-third of the differentially expressed transcripts being down-regulated. However, across treatments we observed shifts in gene expression for common aspects of the heat stress response including intra-cellular communication, protein chaperoning, proteolysis and cell cycle arrest. Of the 71,675 transcripts obtained, only 259 were differentially expressed after both heating events. From these, 218 defined the core group (i.e. genes induced by thermal stress with similar expression patterns irrespective of the magnitude of the warming event). The core group was composed of already well-studied genes in heat stress responses in intertidal organisms (e.g. heat shock proteins), but also genes from less explored metabolic pathways, e.g. the ubiquitin system, which were also fundamental regardless of the magnitude of the imposed warming. Moreover, we have also identified 41 signaling genes (i.e. a set of genes responding to both events and with expression patterns specific to the intensity of thermal stress), principally including genes involved in the maintenance of extracellular structure that have previously not been identified as part of the response to thermal stress in intertidal zone organisms. These signaling genes will be useful heat stress molecular biomarkers for monitoring heat stress in natural populations.

Environmental optima for an ecosystem engineer: a multidisciplinary trait-based approach.

A complex interplay of biotic and abiotic factors underpins the distribution of species and operates across different levels of biological organization and life history stages. Understanding ecosystem engineer reproductive traits is critical for comprehending and managing the biodiversity-rich habitats they create. Little is known about how the reproduction of the reef-forming worm, Sabellaria alveolata, varies across environmental gradients. By integrating broad-scale environmental data with in-situ physiological data in the form of biochemical traits, we identified and ranked the drivers of intraspecific reproductive trait variability (ITV). ITV was highest in locations with variable environmental conditions, subjected to fluctuating temperature and hydrodynamic conditions. Our trait selection pointed to poleward sites being the most physiologically stressful, with low numbers of irregularly shaped eggs suggesting potentially reduced reproductive success. Centre-range individuals allocated the most energy to reproduction, with the highest number of intermediate-sized eggs, whilst equatorward sites were the least physiologically stressful, thus confirming the warm-adapted nature of our model organism. Variation in total egg diameter and relative fecundity were influenced by a combination of environmental conditions, which changed depending on the trait and sampling period. An integrated approach involving biochemical and reproductive traits is essential for understanding macro-scale patterns in the face of anthropogenic-induced climate change across environmental and latitudinal gradients.

A comprehensive assessment of the intertidal biodiversity along the Portuguese coast in the early 2000s.

BACKGROUND: The unprecedented rates of current biodiversity loss have motivated a renewed interest in environmental and biodiversity monitoring. The need for sustained monitoring strategies has prompted not only the establisment of new long-term monitoring programmes, but also the rescue of data from historical or otherwise archived sources. Amongst the most valuable datasets are those containing information on intertidal systems, as they are particularly well suited for studying the biological effects of climate change. The Portuguese rocky coast is quite interesting for studying the effects of climate change on the distribution of species due to its geographical orientation, latitudinal patterns in temperature, species richness, species' distribution patterns and availability of historical information. This work aims at providing a comprehensive picture of the distribution and abundance of intertidal macro-invertebrates and macro-algae along the Portuguese rocky coast in the early 2000s. NEW INFORMATION: This study provides a description of the rocky shore intertidal biodiversity of the mainland Portuguese coast in the early 2000s. The spatial distribution and semi-quantitative abundance of a total of 238 taxa were assessed at 49 wave-exposed locations. These data provide a comprehensive baseline against which biodiversity changes can be effectively and objectively evaluated.

Seascape genomics reveals population isolation in the reef-building honeycomb worm, Sabellaria alveolata (L.).

BACKGROUND: Under the threat of climate change populations can disperse, acclimatise or evolve in order to avoid fitness loss. In light of this, it is important to understand neutral gene flow patterns as a measure of dispersal potential, but also adaptive genetic variation as a measure of evolutionary potential. In order to assess genetic variation and how this relates to environment in the honeycomb worm (Sabellaria alveolata (L.)), a reef-building polychaete that supports high biodiversity, we carried out RAD sequencing using individuals from along its complete latitudinal range. Patterns of neutral population genetic structure were compared to larval dispersal as predicted by ocean circulation modelling, and outlier analyses and genotype-environment association tests were used to attempt to identify loci under selection in relation to local temperature data. RESULTS: We genotyped 482 filtered SNPs, from 68 individuals across nine sites, 27 of which were identified as outliers using BAYESCAN and ARLEQUIN. All outlier loci were potentially under balancing selection, despite previous evidence of local adaptation in the system. Limited gene flow was observed among reef-sites (FST = 0.28 ± 0.10), in line with the low dispersal potential identified by the larval dispersal models. The North Atlantic reef emerged as a distinct population and this was linked to high local larval retention and the effect of the North Atlantic Current on dispersal. CONCLUSIONS: As an isolated population, with limited potential for natural genetic or demographic augmentation from other reefs, the North Atlantic site warrants conservation attention in order to preserve not only this species, but above all the crucial functional ecological roles that are associated with their bioconstructions. Our study highlights the utility of using seascape genomics to identify populations of conservation concern.

Temperature-related heart rate in water and air and a comparison to other temperature-related measures of performance in the fiddler crab Leptuca pugilator (Bosc 1802).

Performance in poikilotherms is known to be sensitive to temperature, often with a low-sloping increase with temperature to a peak, and a steep decline with increasing temperature past the peak. We complemented past measures of performance by measuring heartbeat rates of the fiddler crab Leptuca pugilator in water and in air as a function of a range of temperatures previously shown to affect other measures of performance. In water over a range of 20-50 °C, heartbeat increased steadily to a peak at 40 °C and then steeply declined to near zero at 50 °C. In air, heartbeat also increased, but to a peak at 35 °C and then with a gentler decline than was found in water. Part of this different response may be due to evaporative water loss, which reduced body temperature in air, and therefore thermal stress, relative to body temperature when crabs were immersed in water. Increased availability of oxygen from air, according to the oxygen and capacity-limited thermal tolerance hypothesis, likely increased aerobic scope past the thermal peak, relative to within water, where oxygen delivery at higher temperatures may have been curtailed. We compared the heart rate performance relations to two previous measures of performance - endurance on a treadmill and sprint speed, both done in air. The peak performance temperature increased in the order: treadmill endurance time, sprint speed, heart rate in air, and heart rate in water, which demonstrates that different performance measures give different perspectives on the relation of thermal tolerance and fitness to temperature. Endurance may therefore be the limiting upper thermal stress factor in male fiddler crabs, when on hot sand flats. Temperature preference, found to be for temperatures <30 °C in air, could be a bet-hedging evolutionary strategy to avoid aerobic scope affecting endurance.

Coastal warming and wind-driven upwelling: A global analysis.

Long-term sea surface temperature (SST) warming trends are far from being homogeneous, especially when coastal and ocean locations are compared. Using data from NOAA's AVHRR OISST, we have analyzed sea surface temperature trends over the period 1982-2015 at around 3500 worldwide coastal points and their oceanic counterparts with a spatial resolution of 0.25 arc-degrees. Significant warming was observed at most locations although with important differences between oceanic and coastal points. This is especially patent for upwelling regions, where 92% of the coastal locations showed lower warming trends than at neighboring ocean locations. This result strongly suggests that upwelling has the potential to buffer the effects of global warming nearshore, with wide oceanographic, climatic, and biogeographic implications.

Drivers of Cape Verde archipelagic endemism in keyhole limpets.

Oceanic archipelagos are the ideal setting for investigating processes that shape species assemblages. Focusing on keyhole limpets, genera Fissurella and Diodora from Cape Verde Islands, we used an integrative approach combining molecular phylogenetics with ocean transport simulations to infer species distribution patterns and analyse connectivity. Dispersal simulations, using pelagic larval duration and ocean currents as proxies, showed a reduced level of connectivity despite short distances between some of the islands. It is suggested that dispersal and persistence driven by patterns of oceanic circulation favouring self-recruitment played a primary role in explaining contemporary species distributions. Mitochondrial and nuclear data revealed the existence of eight Cape Verde endemic lineages, seven within Fissurella, distributed across the archipelago, and one within Diodora restricted to Boavista. The estimated origins for endemic Fissurella and Diodora were 10.2 and 6.7 MY, respectively. Between 9.5 and 4.5 MY, an intense period of volcanism in Boavista might have affected Diodora, preventing its diversification. Having originated earlier, Fissurella might have had more opportunities to disperse to other islands and speciate before those events. Bayesian analyses showed increased diversification rates in Fissurella possibly promoted by low sea levels during Plio-Pleistocene, which further explain differences in species richness between both genera.

Response of Two Mytilids to a Heatwave: The Complex Interplay of Physiology, Behaviour and Ecological Interactions.

Different combinations of behavioural and physiological responses may play a crucial role in the ecological success of species, notably in the context of biological invasions. The invasive mussel Xenostrobus securis has successfully colonised the inner part of the Galician Rias Baixas (NW Spain), where it co-occurs with the commercially-important mussel Mytilus galloprovincialis. This study investigated the effect of a heatwave on the physiological and behavioural responses in monospecific or mixed aggregations of these species. In a mesocosm experiment, mussels were exposed to simulated tidal cycles and similar temperature conditions to those experienced in the field during a heat-wave that occurred in the summer of 2013, when field robo-mussels registered temperatures up to 44.5°C at low tide. The overall responses to stress differed markedly between the two species. In monospecific aggregations M. galloprovincialis was more vulnerable than X. securis to heat exposure during emersion. However, in mixed aggregations, the presence of the invader was associated with lower mortality in M. galloprovincialis. The greater sensitivity of M. galloprovincialis to heat exposure was reflected in a higher mortality level, greater induction of Hsp70 protein and higher rates of respiration and gaping activity, which were accompanied by a lower heart rate (bradycardia). The findings show that the invader enhanced the physiological performance of M. galloprovincialis, highlighting the importance of species interactions in regulating responses to environmental stress. Understanding the complex interactions between ecological factors and physiological and behavioural responses of closely-related species is essential for predicting the impacts of invasions in the context of future climate change.

Temperature Regimes Impact Coral Assemblages along Environmental Gradients on Lagoonal Reefs in Belize.

Coral reefs are increasingly threatened by global and local anthropogenic stressors such as rising seawater temperature, nutrient enrichment, sedimentation, and overfishing. Although many studies have investigated the impacts of local and global stressors on coral reefs, we still do not fully understand how these stressors influence coral community structure, particularly across environmental gradients on a reef system. Here, we investigate coral community composition across three different temperature and productivity regimes along a nearshore-offshore gradient on lagoonal reefs of the Belize Mesoamerican Barrier Reef System (MBRS). A novel metric was developed using ultra-high-resolution satellite-derived estimates of sea surface temperatures (SST) to classify reefs as exposed to low (lowTP), moderate (modTP), or high (highTP) temperature parameters over 10 years (2003 to 2012). Coral species richness, abundance, diversity, density, and percent cover were lower at highTP sites relative to lowTP and modTP sites, but these coral community traits did not differ significantly between lowTP and modTP sites. Analysis of coral life history strategies revealed that highTP sites were dominated by hardy stress-tolerant and fast-growing weedy coral species, while lowTP and modTP sites consisted of competitive, generalist, weedy, and stress-tolerant coral species. Satellite-derived estimates of Chlorophyll-a (chl-a) were obtained for 13-years (2003-2015) as a proxy for primary production. Chl-a concentrations were highest at highTP sites, medial at modTP sites, and lowest at lowTP sites. Notably, thermal parameters correlated better with coral community traits between site types than productivity, suggesting that temperature (specifically number of days above the thermal bleaching threshold) played a greater role in defining coral community structure than productivity on the MBRS. Dominance of weedy and stress-tolerant genera at highTP sites suggests that corals utilizing these two life history strategies may be better suited to cope with warmer oceans and thus may warrant protective status under climate change.

Equatorial range limits of an intertidal ectotherm are more linked to water than air temperature.

As climate change is expected to impose increasing thermal stress on intertidal organisms, understanding the mechanisms by which body temperatures translate into major biogeographic patterns is of paramount importance. We exposed individuals of the limpet Patella vulgata Linnaeus, 1758, to realistic experimental treatments aimed at disentangling the contribution of water and air temperature for the buildup of thermal stress. Treatments were designed based on temperature data collected at the microhabitat level, from 15 shores along the Atlantic European coast spanning nearly 20° of latitude. Cardiac activity data indicated that thermal stress levels in P. vulgata are directly linked to elevated water temperature, while high air temperature is only stressful if water temperature is also high. In addition, the analysis of the link between population densities and thermal regimes at the studied locations suggests that the occurrence of elevated water temperature may represent a threshold P. vulgata is unable to tolerate. By combining projected temperatures with the temperature threshold identified, we show that climate change will likely result in the westward expansion of the historical distribution gap in the Bay of Biscay (southwest France), and northward contraction of the southern range limit in south Portugal. These findings suggest that even a minor relaxing of the upwelling off northwest Iberia could lead to a dramatic increase in thermal stress, with major consequences for the structure and functioning of the intertidal communities along Iberian rocky shores.

Exposure to solar radiation drives organismal vulnerability to climate: Evidence from an intertidal limpet.

Understanding the physiological abilities of organisms to cope with heat stress is critical for predictions of species' distributions in response to climate change. We investigated physiological responses (respiration and heart beat rate) of the ectotherm limpet Patella vulgata to heat stress events during emersion and the role of seasonal and microclimatic acclimatization for individual thermal tolerance limits. Individuals were collected from 5 microhabitats characterized by different exposure to solar radiation in the high intertidal zone of a semi-exposed rocky shore in winter and summer of 2014. Upper thermal tolerance limits (heat coma temperatures - HCTs, and heart rate Arrhenius break temperatures - ABTs) were determined for individuals from each microhabitat in both seasons under laboratory conditions. While we found a clear seasonal acclimatization, i.e., higher HCTs and ABTs in summer than in winter, we did not find evidence for microhabitat-specific responses that would suggest microclimatic acclimatization. However, operative limpet temperatures derived from in-situ temperature measurements suggest that individuals from sun exposed microhabitats have a much narrower thermal safety margins than those from less exposed surfaces or within crevices. Microhabitat specific thermal safety margins caused by high thermal heterogeneity at small spatial scales and the lack of short term acclimatization will likely shape small scale distribution patterns of intertidal species in response to the predicted increase in the frequency and intensity of heat waves.

Ocean-wide tracking of pelagic sharks reveals extent of overlap with longline fishing hotspots.

Overfishing is arguably the greatest ecological threat facing the oceans, yet catches of many highly migratory fishes including oceanic sharks remain largely unregulated with poor monitoring and data reporting. Oceanic shark conservation is hampered by basic knowledge gaps about where sharks aggregate across population ranges and precisely where they overlap with fishers. Using satellite tracking data from six shark species across the North Atlantic, we show that pelagic sharks occupy predictable habitat hotspots of high space use. Movement modeling showed sharks preferred habitats characterized by strong sea surface-temperature gradients (fronts) over other available habitats. However, simultaneous Global Positioning System (GPS) tracking of the entire Spanish and Portuguese longline-vessel fishing fleets show an 80% overlap of fished areas with hotspots, potentially increasing shark susceptibility to fishing exploitation. Regions of high overlap between oceanic tagged sharks and longliners included the North Atlantic Current/Labrador Current convergence zone and the Mid-Atlantic Ridge southwest of the Azores. In these main regions, and subareas within them, shark/vessel co-occurrence was spatially and temporally persistent between years, highlighting how broadly the fishing exploitation efficiently "tracks" oceanic sharks within their space-use hotspots year-round. Given this intense focus of longliners on shark hotspots, our study argues the need for international catch limits for pelagic sharks and identifies a future role of combining fine-scale fish and vessel telemetry to inform the ocean-scale management of fisheries.

Loss of thermal refugia near equatorial range limits.

This study examines the importance of thermal refugia along the majority of the geographical range of a key intertidal species (Patella vulgata Linnaeus, 1758) on the Atlantic coast of Europe. We asked whether differences between sun-exposed and shaded microhabitats were responsible for differences in physiological stress and ecological performance and examined the availability of refugia near equatorial range limits. Thermal differences between sun-exposed and shaded microhabitats are consistently associated with differences in physiological performance, and the frequency of occurrence of high temperatures is most probably limiting the maximum population densities supported at any given place. Topographical complexity provides thermal refugia throughout most of the distribution range, although towards the equatorial edges the magnitude of the amelioration provided by shaded microhabitats is largely reduced. Importantly, the limiting effects of temperature, rather than being related to latitude, seem to be tightly associated with microsite variability, which therefore is likely to have profound effects on the way local populations (and consequently species) respond to climatic changes.

Understanding complex biogeographic responses to climate change.

Predicting the extent and direction of species' range shifts is a major priority for scientists and resource managers. Seminal studies have fostered the notion that biological systems responding to climate change-impacted variables (e.g., temperature, precipitation) should exhibit poleward range shifts but shifts contrary to that expectation have been frequently reported. Understanding whether those shifts are indeed contrary to climate change predictions involves understanding the most basic mechanisms determining the distribution of species. We assessed the patterns of ecologically relevant temperature metrics (e.g., daily range, min, max) along the European Atlantic coast. Temperature metrics have contrasting geographical patterns and latitude or the grand mean are poor predictors for many of them. Our data suggest that unless the appropriate metrics are analysed, the impact of climate change in even a single metric of a single stressor may lead to range shifts in directions that would otherwise be classified as "contrary to prediction".

Evolution at a different pace: distinctive phylogenetic patterns of cone snails from two ancient oceanic archipelagos.

Ancient oceanic archipelagos of similar geological age are expected to accrue comparable numbers of endemic lineages with identical life history strategies, especially if the islands exhibit analogous habitats. We tested this hypothesis using marine snails of the genus Conus from the Atlantic archipelagos of Cape Verde and Canary Islands. Together with Azores and Madeira, these archipelagos comprise the Macaronesia biogeographic region and differ remarkably in the diversity of this group. More than 50 endemic Conus species have been described from Cape Verde, whereas prior to this study, only two nonendemic species, including a putative species complex, were thought to occur in the Canary Islands. We combined molecular phylogenetic data and geometric morphometrics with bathymetric and paleoclimatic reconstructions to understand the contrasting diversification patterns found in these regions. Our results suggest that species diversity is even lower than previously thought in the Canary Islands, with the putative species complex corresponding to a single species, Conus guanche. One explanation for the enormous disparity in Conus diversity is that the amount of available habitat may differ, or may have differed in the past due to eustatic (global) sea level changes. Historical bathymetric data, however, indicated that sea level fluctuations since the Miocene have had a similar impact on the available habitat area in both Cape Verde and Canary archipelagos and therefore do not explain this disparity. We suggest that recurrent gene flow between the Canary Islands and West Africa, habitat losses due to intense volcanic activity in combination with unsuccessful colonization of new Conus species from more diverse regions, were all determinant in shaping diversity patterns within the Canarian archipelago. Worldwide Conus species diversity follows the well-established pattern of latitudinal increase of species richness from the poles towards the tropics. However, the eastern Atlantic revealed a striking pattern with two main peaks of Conus species richness in the subtropical area and decreasing diversities toward the tropical western African coast. A Random Forests model using 12 oceanographic variables suggested that sea surface temperature is the main determinant of Conus diversity either at continental scales (eastern Atlantic coast) or in a broader context (worldwide). Other factors such as availability of suitable habitat and reduced salinity due to the influx of large rivers in the tropical area also play an important role in shaping Conus diversity patterns in the western coast of Africa.

Thermal adaptation and clinal mitochondrial DNA variation of European anchovy.

Natural populations of widely distributed organisms often exhibit genetic clinal variation over their geographical ranges. The European anchovy, Engraulis encrasicolus, illustrates this by displaying a two-clade mitochondrial structure clinally arranged along the eastern Atlantic. One clade has low frequencies at higher latitudes, whereas the other has an anti-tropical distribution, with frequencies decreasing towards the tropics. The distribution pattern of these clades has been explained as a consequence of secondary contact after an ancient geographical isolation. However, it is not unlikely that selection acts on mitochondria whose genes are involved in relevant oxidative phosphorylation processes. In this study, we performed selection tests on a fragment of 1044 bp of the mitochondrial cytochrome b gene using 455 individuals from 18 locations. We also tested correlations of six environmental features: temperature, salinity, apparent oxygen utilization and nutrient concentrations of phosphate, nitrate and silicate, on a compilation of mitochondrial clade frequencies from 66 sampling sites comprising 2776 specimens from previously published studies. Positive selection in a single codon was detected predominantly (99%) in the anti-tropical clade and temperature was the most relevant environmental predictor, contributing with 59% of the variance in the geographical distribution of clade frequencies. These findings strongly suggest that temperature is shaping the contemporary distribution of mitochondrial DNA clade frequencies in the European anchovy.

Three decades of high-resolution coastal sea surface temperatures reveal more than warming.

Understanding and forecasting current and future consequences of coastal warming require a fine-scale assessment of the near-shore temperature changes. Here we show that despite the fact that 71% of the world's coastlines are significantly warming, rates of change have been highly heterogeneous both spatially and seasonally. We demonstrate that 46% of the coastlines have experienced a significant decrease in the frequency of extremely cold events, while extremely hot days are becoming more common in 38% of the area. Also, we show that the onset of the warm season is significantly advancing earlier in the year in 36% of the temperate coastal regions. More importantly, it is now possible to analyse local patterns within the global context, which is useful for a broad array of scientific fields, policy makers and general public.