Oceanographic Conditions Limit the Spread of a Marine Invader along Southern African Shores

Long range gene flow beyond predictions from oceanographic transport in a tropical marine foundation species

The transport of passively dispersed organisms across tropical margins remains poorly understood. Hypotheses of oceanographic transportation potential lack testing with large scale empirical data. To address this gap, we used the seagrass species, Halodule wrightii, which is unique in spanning the entire tropical Atlantic. We tested the hypothesis that genetic differentiation estimated across its large-scale biogeographic range can be predicted by simulated oceanographic transport. The alternative hypothesis posits that dispersal is independent of ocean currents, such as transport by grazers. We compared empirical genetic estimates and modelled predictions of dispersal along the distribution of H. wrightii. We genotyped eight microsatellite loci on 19 populations distributed across Atlantic Africa, Gulf of Mexico, Caribbean, Brazil and developed a biophysical model with high-resolution ocean currents. Genetic data revealed low gene flow and highest differentiation between (1) the Gulf of Mexico and two other regions: (2) Caribbean-Brazil and (3) Atlantic Africa. These two were more genetically similar despite separation by an ocean. The biophysical model indicated low or no probability of passive dispersal among populations and did not match the empirical genetic data. The results support the alternative hypothesis of a role for active dispersal vectors like grazers.

Transcriptome wide analyses reveal intraspecific diversity in thermal stress responses of a dominant habitat-forming species

The impact of climate change on biodiversity has stimulated the need to understand environmental stress responses, particularly for ecosystem engineers whose responses to climate affect large numbers of associated organisms. Distinct species differ substantially in their resilience to thermal stress but there are also within-species variations in thermal tolerance for which the molecular mechanisms underpinning such variation remain largely unclear. Intertidal mussels are well-known for their role as ecosystem engineers. First, we exposed two genetic lineages of the intertidal mussel Perna perna to heat stress treatments in air and water. Next, we ran a high throughput RNA sequencing experiment to identify differences in gene expression between the thermally resilient eastern lineage and the thermally sensitive western lineage. We highlight different thermal tolerances that concord with their distributional ranges. Critically, we also identified lineage-specific patterns of gene expression under heat stress and revealed intraspecific differences in the underlying transcriptional pathways in response to warmer temperatures that are potentially linked to the within-species differences in thermal tolerance. Beyond the species, we show how unravelling within-species variability in mechanistic responses to heat stress promotes a better understanding of global evolutionary trajectories of the species as a whole in response to changing climate.

Thermal stress gradient causes increasingly negative effects towards the range limit of an invasive mussel

Environmental filtering (EF), the abiotic exclusion of species, can have first order, direct effects with cascading consequences for population dynamics, especially at range edges where abiotic conditions are suboptimal. Abiotic stress gradients associated with EF may also drive indirect second order effects, including exacerbating the effects of competitors, disease, and parasites on marginal populations because of suboptimal physiological performance. We predicted a cascade of first and second order EF-associated effects on marginal populations of the invasive mussel Mytilus galloprovincialis, plus a third order effect of EF of increased epibiont load due to second order shell degradation by endoliths. Mussel populations on rocky shores were surveyed across 850 km of the south-southeast coast of South Africa, from the species' warm-edge range limit to sites in the centre of their distribution, to quantify second order (endolithic shell degradation) and third order (number of barnacle epibionts) EF-associated effects as a function of along-shore distance from the range edge. Inshore temperature data were interpolated from the literature. Using in situ temperature logger data, we calculated the effective shore level for several sites by determining the duration of immersion and emersion. Summer and winter inshore water temperatures were linked to distance from the mussel's warm range edge (our proxy for an EF-associated stress gradient), suggesting that seasonality in temperature contributes to first order effects. The gradient in thermal stress clearly affected densities, but its influence on mussel size, shell degradation, and epibiosis was weaker. Relationships among mussel size, shell degradation, and epibiosis were more robust. Larger, older mussels had more degraded shells and more epibionts, with endolithic damage facilitating epibiosis. EF associated with a gradient in thermal stress directly limits the distribution, abundance, and size structure of mussel populations, with important indirect second and third order effects of parasitic disease and epibiont load, respectively.

Weak biodiversity connectivity in the European network of no-take marine protected areas

The need for international cooperation in marine resource management and conservation has been reflected in the increasing number of agreements aiming for effective and well-connected networks of Marine Protected Areas (MPAs). However, the extent to which individual MPAs are connected remains mostly unknown. Here, we use a biophysical model tuned with empirical data on species dispersal ecology to predict connectivity of a vast spectrum of biodiversity in the European network of marine reserves (i.e., no-take MPAs). Our results highlight the correlation between empirical propagule duration data and connectivity potential and show weak network connectivity and strong isolation for major ecological groups, resulting from the lack of direct connectivity corridors between reserves over vast regions. The particularly high isolation predicted for ecosystem structuring species (e.g., corals, sponges, macroalgae and seagrass) might potentially undermine biodiversity conservation efforts if local retention is insufficient and unmanaged populations are at risk. Isolation might also be problematic for populations' persistence in the light of climate change and expected species range shifts. Our findings provide novel insights for management directives, highlighting the location of regions requiring additional marine reserves to function as stepping-stone connectivity corridors.

Historical and contemporary range expansion of an invasive mussel, Semimytlius algosus, in Angola and Namibia despite data scarcity in an infrequently surveyed region

Understanding the spread of invasive species in many regions is difficult because surveys are rare. Here, historical records of the invasive marine mussel, Semimytilus algosus, on the shores of Angola and Namibia are synthesised to re-construct its invasive history. Since this mussel was first discovered in Namibia about 90 years ago, it has spread throughout the western coast of southern Africa. By the late 1960s, the species was well established across a range of 1005 km of coastline in southern Angola and northern Namibia. Although only coarse spatial resolution data are available since the 1990s, the distribution of S. algosus clearly increased substantially over the subsequent decades. Today, the species is distributed over 2785 km of coastline, appearing in southern Namibia in 2014, whence it spread across the border to northern South Africa in 2017, and in northern Angola in 2015. Conspicuously, its current range appears to be relatively contiguous across at least 810 km of shore in southern Angola and throughout Namibia, with isolated, spatially disjunct occurrences towards the southern and northern limits of its distribution. Despite there being few occurrence records that are unevenly distributed spatially and temporally, data for the distributional patterns of S. algosus in Angola and Namibia provide invaluable insights into how marine invasive species spread in developing regions that are infrequently monitored.

A fine-tuned global distribution dataset of marine forests

Species distribution records are a prerequisite to follow climate-induced range shifts across space and time. However, synthesizing information from various sources such as peer-reviewed literature, herbaria, digital repositories and citizen science initiatives is not only costly and time consuming, but also challenging, as data may contain thematic and taxonomic errors and generally lack standardized formats. We address this gap for important marine ecosystem-structuring species of large brown algae and seagrasses. We gathered distribution records from various sources and provide a fine-tuned dataset with ~2.8 million dereplicated records, taxonomically standardized for 682 species, and considering important physiological and biogeographical traits. Specifically, a flagging system was implemented to signal potentially incorrect records reported on land, in regions with limiting light conditions for photosynthesis, and outside the known distribution of species, as inferred from the most recent published literature. We document the procedure and provide a dataset in tabular format based on Darwin Core Standard (DwC), alongside with a set of functions in R language for data management and visualization.

Contemporary climate change hinders hybrid performance of ecologically dominant marine invertebrates

Human activities alter patterns of biodiversity, particularly through species extinctions and range shifts. Two of these activities are human mediated transfer of species and contemporary climate change, and both allow previously isolated genotypes to come into contact and hybridize, potentially altering speciation rates. Hybrids have been shown to survive environmental conditions not tolerated by either parent, suggesting that, under some circumstances, hybrids may be able to expand their ranges and perform well under rapidly changing conditions. However, studies assessing how hybridization influences contemporary range shifts are scarce. We performed crosses on Pyura herdmani and Pyura stolonifera (Chordata, Tunicata), two closely related marine invertebrate species that are ecologically dominant and can hybridize. These sister species live in sympatry along the coasts of southern Africa, but one has a disjunct distribution that includes northern hemisphere sites. We experimentally assessed the performance of hybrid and parental crosses using different temperature regimes, including temperatures predicted under future climate change scenarios. We found that hybrids showed lower performance than parental crosses at the experimental temperatures, suggesting that hybrids are unlikely to expand their ranges to new environments. In turn, we found that the more widespread species performed better at a wide array of temperatures, indicating that this parental species may cope better with future conditions. This study illustrates how offspring fitness may provide key insights to predict range expansions and how contemporary climate change may mediate both the ability of hybrids to expand their ranges and the occurrence of speciation as a result of hybridization.

Rejection of the genetic implications of the "Abundant Centre Hypothesis" in marine mussels

The ‘Abundant-Centre Hypothesis’ is a well-established but controversial hypothesis stating that the abundance of a species is highest at the centre of its range and decreases towards the edges, where conditions are unfavourable. As genetic diversity depends on population size, edge populations are expected to show lower intra-population genetic diversity than core populations, while showing high inter-population genetic divergence. Here, the genetic implications of the Abundant-Centre Hypothesis were tested on two coastal mussels from South Africa that disperse by means of planktonic larvae, the native Perna perna and the invasive Mytilus galloprovincialis. Genetic structure was found within P. perna, which, together with evidence from Lagrangian particle simulations, points to significant reductions in gene flow between sites. Despite this, the expected diversity pattern between centre and edge populations was not found for either species. We conclude that the genetic predictions of the Abundant-Centre Hypothesis are unlikely to be met by high-dispersal species with large population sizes, and may only become evident in species with much lower levels of connectivity.

Biogeographical Patterns of Endolithic Infestation in an Invasive and an Indigenous Intertidal Marine Ecosystem Engineer

By altering the phenotypic properties of their hosts, endolithic parasites can modulate the engineering processes of marine ecosystem engineers. Here, we assessed the biogeographical patterns of species assemblages, prevalence and impact of endolithic parasitism in two mussel species that act as important ecosystem engineers in the southern African intertidal habitat, Perna perna and Mytilus galloprovincialis. We conducted large-scale surveys across three biogeographic regions along the South African coast: the subtropical east coast, dominated by the indigenous mussel, P. perna, the warm temperate south coast, where this species coexists with the invasive Mediterranean mussel, M. galloprovincialis, and the cool temperate west coast dominated by M. galloprovincialis. Infestation increased with mussel size, and in the case of M. galloprovincialis we found a significantly higher infestation in the cool temperate bioregion than the warm temperate region. For P. perna, the prevalence of infestation was higher on the warm temperate than the subtropical region, though the difference was marginally non-significant. On the south coast, there was no significant difference in infestation prevalence between species. Endolith-induced mortality rates through shell collapse mirrored the patterns for prevalence. For P. perna, endolith species assemblages revealed clear grouping by bioregions. Our findings indicate that biogeography affects cyanobacteria species composition, but differences between biogeographic regions in their effects are driven by environmental conditions.

Thermal selection as a driver of marine ecological speciation

Intraspecific genetic structure in widely distributed marine species often mirrors the boundaries between temperature-defined bioregions. This suggests that the same thermal gradients that maintain distinct species assemblages also drive the evolution of new biodiversity. Ecological speciation scenarios are often invoked to explain such patterns, but the fact that adaptation is usually only identified when phylogenetic splits are already evident makes it impossible to rule out the alternative scenario of allopatric speciation with subsequent adaptation. We integrated large-scale genomic and environmental datasets along one of the world's best-defined marine thermal gradients (the South African coastline) to test the hypothesis that incipient ecological speciation is a result of divergence linked to the thermal environment. We identified temperature-associated gene regions in a coastal fish species that is spatially homogeneous throughout several temperature-defined biogeographic regions based on selectively neutral markers. Based on these gene regions, the species is divided into geographically distinct regional populations. Importantly, the ranges of these populations are delimited by the same ecological boundaries that define distinct infraspecific genetic lineages in co-distributed marine species, and biogeographic disjunctions in species assemblages. Our results indicate that temperature-mediated selection represents an early stage of marine ecological speciation in coastal regions that lack physical dispersal barriers.

Managing an invasive corallimorph at Palmyra Atoll National Wildlife Refuge, Line Islands, Central Pacific

In 2007, a phase shift from corals to corallimorpharians (CM) centered around a shipwreck was documented at Palmyra Atoll, Line Islands. Subsequent surveys revealed CM to be overgrowing the reef benthos, including corals and coralline algae, potentially placing coral ecosystems in the atoll at risk. This prompted the U.S. Fish and Wildlife Service, the lead management agency of the atoll, to remove the shipwreck. Subsequent surveys showed reductions in CM around the ship impact site. We explain patterns of spread of the CM in terms of both life history and local currents and show with a pilot study that pulverized bleach may be an effective tool to eradicate CM on a local scale. If applied strategically, particularly in heavily infested (> 66% cover) areas, active intervention such as this could be an effective management tool to reduce CM impact on localized areas and decrease colonization rate of remaining reefs. This is the first documentation of the response of an invasive cnidarian to shipwreck removal. While this was a singular event in Palmyra, the spatial and temporal patterns of this invasion and the eradications lessons described herein, are useful for anticipating and controlling similar situations elsewhere.

Evidence for rangewide panmixia despite multiple barriers to dispersal in a marine mussel

Oceanographic features shape the distributional and genetic patterns of marine species by interrupting or promoting connections among populations. Although general patterns commonly arise, distributional ranges and genetic structure are species-specific and do not always comply with the expected trends. By applying a multimarker genetic approach combined with Lagrangian particle simulations (LPS) we tested the hypothesis that oceanographic features along northeastern Atlantic and Mediterranean shores influence dispersal potential and genetic structure of the intertidal mussel Perna perna. Additionally, by performing environmental niche modelling we assessed the potential and realized niche of P. perna along its entire native distributional range and the environmental factors that best explain its realized distribution. Perna perna showed evidence of panmixia across >4,000 km despite several oceanographic breaking points detected by LPS. This is probably the result of a combination of life history traits, continuous habitat availability and stepping-stone dynamics. Moreover, the niche modelling framework depicted minimum sea surface temperatures (SST) as the major factor shaping P. perna distributional range limits along its native areas. Forthcoming warming SST is expected to further change these limits and allow the species to expand its range polewards though this may be accompanied by retreat from warmer areas.

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.

Habitat continuity and stepping-stone oceanographic distances explain population genetic connectivity of the brown alga Cystoseira amentacea

Effective predictive and management approaches for species occurring in a metapopulation structure require good understanding of interpopulation connectivity. In this study, we ask whether population genetic structure of marine species with fragmented distributions can be predicted by stepping-stone oceanographic transport and habitat continuity, using as model an ecosystem-structuring brown alga, Cystoseira amentacea var. stricta. To answer this question, we analysed the genetic structure and estimated the connectivity of populations along discontinuous rocky habitat patches in southern Italy, using microsatellite markers at multiple scales. In addition, we modelled the effect of rocky habitat continuity and ocean circulation on gene flow by simulating Lagrangian particle dispersal based on ocean surface currents allowing multigenerational stepping-stone dynamics. Populations were highly differentiated, at scales from few metres up to thousands of kilometres. The best possible model fit to explain the genetic results combined current direction, rocky habitat extension and distance along the coast among rocky sites. We conclude that a combination of variable suitable habitat and oceanographic transport is a useful predictor of genetic structure. This relationship provides insight into the mechanisms of dispersal and the role of life-history traits. Our results highlight the importance of spatially explicit modelling of stepping-stone dynamics and oceanographic directional transport coupled with habitat suitability, to better describe and predict marine population structure and differentiation. This study also suggests the appropriate spatial scales for the conservation, restoration and management of species that are increasingly affected by habitat modifications.

High Interannual Variability in Connectivity and Genetic Pool of a Temperate Clingfish Matches Oceanographic Transport Predictions

Adults of most marine benthic and demersal fish are site-attached, with the dispersal of their larval stages ensuring connectivity among populations. In this study we aimed to infer spatial and temporal variation in population connectivity and dispersal of a marine fish species, using genetic tools and comparing these with oceanographic transport. We focused on an intertidal rocky reef fish species, the shore clingfish Lepadogaster lepadogaster, along the southwest Iberian Peninsula, in 2011 and 2012. We predicted high levels of self-recruitment and distinct populations, due to short pelagic larval duration and because all its developmental stages have previously been found near adult habitats. Genetic analyses based on microsatellites countered our prediction and a biophysical dispersal model showed that oceanographic transport was a good explanation for the patterns observed. Adult sub-populations separated by up to 300 km of coastline displayed no genetic differentiation, revealing a single connected population with larvae potentially dispersing long distances over hundreds of km. Despite this, parentage analysis performed on recruits from one focal site within the Marine Park of Arrábida (Portugal), revealed self-recruitment levels of 2.5% and 7.7% in 2011 and 2012, respectively, suggesting that both long- and short-distance dispersal play an important role in the replenishment of these populations. Population differentiation and patterns of dispersal, which were highly variable between years, could be linked to the variability inherent in local oceanographic processes. Overall, our measures of connectivity based on genetic and oceanographic data highlight the relevance of long-distance dispersal in determining the degree of connectivity, even in species with short pelagic larval durations.

Future climate change is predicted to shift long-term persistence zones in the cold-temperate kelp Laminaria hyperborea

Global climate change is shifting species distributions worldwide. At rear edges (warmer, low latitude range margins), the consequences of small variations in environmental conditions can be magnified, producing large negative effects on species ranges. A major outcome of shifts in distributions that only recently received attention is the potential to reduce the levels of intra-specific diversity and consequently the global evolutionary and adaptive capacity of species to face novel disturbances. This is particularly important for low dispersal marine species, such as kelps, that generally retain high and unique genetic diversity at rear ranges resulting from long-term persistence, while ranges shifts during climatic glacial/interglacial cycles. Using ecological niche modelling, we (1) infer the major environmental forces shaping the distribution of a cold-temperate kelp, Laminaria hyperborea (Gunnerus) Foslie, and we (2) predict the effect of past climate changes in shaping regions of long-term persistence (i.e., climatic refugia), where this species might hypothetically harbour higher genetic diversity given the absence of bottlenecks and local extinctions over the long term. We further (3) assessed the consequences of future climate for the fate of L. hyperborea using different scenarios of greenhouse gas emissions (RCP 2.6 and RCP 8.5). Results show NW Iberia, SW Ireland and W English Channel, Faroe Islands and S Iceland, as regions where L. hyperborea may have persisted during past climate extremes until present day. All predictions for the future showed expansions to northern territories coupled with the significant loss of suitable habitats at low latitude range margins, where long-term persistence was inferred (e.g., NW Iberia). This pattern was particularly evident in the most agressive scenario of climate change (RCP 8.5), likely driving major biodiversity loss, changes in ecosystem functioning and the impoverishment of the global gene pool of L. hyperborea. Because no genetic baseline is currently available for this species, our results may represent a first step in informing conservation and mitigation strategies.

Seascape drivers of Macrocystis pyrifera population genetic structure in the northeast Pacific

At small spatial and temporal scales, genetic differentiation is largely controlled by constraints on gene flow, while genetic diversity across a species' distribution is shaped on longer temporal and spatial scales. We assess the hypothesis that oceanographic transport and other seascape features explain different scales of genetic structure of giant kelp, Macrocystis pyrifera. We followed a hierarchical approach to perform a microsatellite-based analysis of genetic differentiation in Macrocystis across its distribution in the northeast Pacific. We used seascape genetic approaches to identify large-scale biogeographic population clusters and investigate whether they could be explained by oceanographic transport and other environmental drivers. We then modelled population genetic differentiation within clusters as a function of oceanographic transport and other environmental factors. Five geographic clusters were identified: Alaska/Canada, central California, continental Santa Barbara, California Channel Islands and mainland southern California/Baja California peninsula. The strongest break occurred between central and southern California, with mainland Santa Barbara sites forming a transition zone between the two. Breaks between clusters corresponded approximately to previously identified biogeographic breaks, but were not solely explained by oceanographic transport. An isolation-by-environment (IBE) pattern was observed where the northern and southern Channel Islands clustered together, but not with closer mainland sites, despite the greater distance between them. The strongest environmental association with this IBE pattern was observed with light extinction coefficient, which extends suitable habitat to deeper areas. Within clusters, we found support for previous results showing that oceanographic connectivity plays an important role in the population genetic structure of Macrocystis in the Northern hemisphere.