Oceanographic connectivity explains the intra-specific diversity of mangrove forests at global scales.

The distribution of mangrove intra-specific biodiversity can be structured by historical demographic processes that enhance or limit effective population sizes. Oceanographic connectivity (OC) may further structure intra-specific biodiversity by preserving or diluting the genetic signatures of historical changes. Despite its relevance for biogeography and evolution, the role of oceanographic connectivity in structuring the distribution of mangrove's genetic diversity has not been addressed at global scale. Here we ask whether connectivity mediated by ocean currents explains the intra-specific diversity of mangroves. A comprehensive dataset of population genetic differentiation was compiled from the literature. Multigenerational connectivity and population centrality indices were estimated with biophysical modeling coupled with network analyses. The variability explained in genetic differentiation was tested with competitive regression models built upon classical isolation-by-distance (IBD) models considering geographic distance. We show that oceanographic connectivity can explain the genetic differentiation of mangrove populations regardless of the species, region, and genetic marker (significant regression models in 95% of cases, with an average R-square of 0.44 ± 0.23 and Person's correlation of 0.65 ± 0.17), systematically improving IBD models. Centrality indices, providing information on important stepping-stone sites between biogeographic regions, were also important in explaining differentiation (R-square improvement of 0.06 ± 0.07, up to 0.42). We further show that ocean currents produce skewed dispersal kernels for mangroves, highlighting the role of rare long-distance dispersal events responsible for historical settlements. Overall, we demonstrate the role of oceanographic connectivity in structuring mangrove intra-specific diversity. Our findings are critical for mangroves' biogeography and evolution, but also for management strategies considering climate change and genetic biodiversity conservation.

A low-latitude species pump: Peripheral isolation, parapatric speciation and mating-system evolution converge in a marine radiation.

Geologically recent radiations can shed light on speciation processes, but incomplete lineage sorting and introgressive gene flow render accurate evolutionary reconstruction and interpretation challenging. Independently evolving metapopulations of low dispersal taxa may provide an additional level of phylogeographic information, given sufficiently broad sampling and genome-wide sequencing. Evolution in the marine brown algal genus Fucus in the south-eastern North Atlantic was shaped by Quaternary climate-driven range shifts. Over this timescale, divergence and speciation occurred against a background of expansion-contraction cycles from multiple refugia, together with mating-system shifts from outcrossing (dioecy) to selfing hermaphroditism. We tested the hypothesis that peripheral isolation of range edge (dioecious) F. vesiculosus led to parapatric speciation and radiation of hermaphrodite lineages. Species tree methods using 876 single-copy nuclear genes and extensive geographic coverage produced conflicting topologies with respect to geographic clades of F. vesiculosus. All methods, however, revealed a new and early diverging hermaphrodite species, Fucus macroguiryi sp. nov. Both the multispecies coalescent and polymorphism-aware models (in contrast to concatenation) support sequential paraphyly in F. vesiculosus resulting from distinct evolutionary processes. Our results support (1) peripheral isolation of the southern F. vesiculosus clade prior to parapatric speciation and radiation of hermaphrodite lineages-a "low-latitude species pump". (2) Directional introgressive gene flow into F. vesiculosus around the present-day secondary contact zone (sympatric-allopatric boundary) between dioecious/hermaphrodite lineages as hermaphrodites expanded northwards, supported by concordance analysis and statistical tests of introgression. (3) Species boundaries in the extensive sympatric range are probably maintained by reproductive system (selfing in hermaphrodites) and reinforcement.

Harnessing positive species interactions as a tool against climate-driven loss of coastal biodiversity.

Habitat-forming species sustain biodiversity and ecosystem functioning in harsh environments through the amelioration of physical stress. Nonetheless, their role in shaping patterns of species distribution under future climate scenarios is generally overlooked. Focusing on coastal systems, we assess how habitat-forming species can influence the ability of stress-sensitive species to exhibit plastic responses, adapt to novel environmental conditions, or track suitable climates. Here, we argue that habitat-former populations could be managed as a nature-based solution against climate-driven loss of biodiversity. Drawing from different ecological and biological disciplines, we identify a series of actions to sustain the resilience of marine habitat-forming species to climate change, as well as their effectiveness and reliability in rescuing stress-sensitive species from increasingly adverse environmental conditions.

Mates Matter: Gametophyte Kinship Recognition and Inbreeding in the Giant Kelp, Macrocystis pyrifera (Laminariales, Phaeophyceae).

Inbreeding, the mating between genetically related individuals, often results in reduced survival and fecundity of offspring, relative to outcrossing. Yet, high inbreeding rates are commonly observed in seaweeds, suggesting compensatory reproductive traits may affect the costs and benefits of the mating system. We experimentally manipulated inbreeding levels in controlled crossing experiments, using gametophytes from 19 populations of Macrocystis pyrifera along its Eastern Pacific coastal distribution (EPC). The objective was to investigate the effects of male-female kinship on female fecundity and fertility, to estimate inbreeding depression in the F1 progeny, and to assess the variability of these effects among different regions and habitats of the EPC. Results revealed that the presence and kinship of males had a significant effect on fecundity and fertility of female gametophytes. Females left alone or in the presence of sibling males express the highest gametophyte size, number, and size of oogonia, suggesting they were able to sense the presence and the identity of their mates before gamete contact. The opposite trend was observed for the production of embryos per female gametes, indicating higher costs of selfing and parthenogenesis than outcrossing on fertility. However, the increased fecundity compensated for the reduced fertility, leading to a stable overall reproductive output. Inbreeding also affected morphological traits of juvenile sporophytes, but not their heatwave tolerance. The male-female kinship effect was stronger in high-latitude populations, suggesting that females from low-latitude marginal populations might have evolved to mate with any male gamete to guarantee reproductive success.

The microbiome of the habitat-forming brown alga Fucus vesiculosus (Phaeophyceae) has similar cross-Atlantic structure that reflects past and present drivers1.

Latitudinal diversity gradients have provided many insights into species differentiation and community processes. In the well-studied intertidal zone, however, little is known about latitudinal diversity in microbiomes associated with habitat-forming hosts. We investigated microbiomes of Fucus vesiculosus because of deep understanding of this model system and its latitudinally large, cross-Atlantic range. Given multiple effects of photoperiod, we predicted that cross-Atlantic microbiomes of the Fucus microbiome would be similar at similar latitudes and correlate with environmental factors. We found that community structure and individual amplicon sequencing variants (ASVs) showed distinctive latitudinal distributions, but alpha diversity did not. Latitudinal differentiation was mostly driven by ASVs that were more abundant in cold temperate to subarctic (e.g., Granulosicoccus_t3260, Burkholderia/Caballeronia/Paraburkholderia_t8371) or warm temperate (Pleurocapsa_t10392) latitudes. Their latitudinal distributions correlated with different humidity, tidal heights, and air/sea temperatures, but rarely with irradiance or photoperiod. Many ASVs in potentially symbiotic genera displayed novel phylogenetic biodiversity with differential distributions among tissues and regions, including closely related ASVs with differing north-south distributions that correlated with Fucus phylogeography. An apparent southern range contraction of F. vesiculosus in the NW Atlantic on the North Carolina coast mimics that recently observed in the NE Atlantic. We suggest cross-Atlantic microbial structure of F. vesiculosus is related to a combination of past (glacial-cycle) and contemporary environmental drivers.

Projected climate changes threaten ancient refugia of kelp forests in the North Atlantic.

Intraspecific genetic variability is critical for species adaptation and evolution and yet it is generally overlooked in projections of the biological consequences of climate change. We ask whether ongoing climate changes can cause the loss of important gene pools from North Atlantic relict kelp forests that persisted over glacial-interglacial cycles. We use ecological niche modelling to predict genetic diversity hotspots for eight species of large brown algae with different thermal tolerances (Arctic to warm temperate), estimated as regions of persistence throughout the Last Glacial Maximum (20,000 YBP), the warmer Mid-Holocene (6,000 YBP), and the present. Changes in the genetic diversity within ancient refugia were projected for the future (year 2100) under two contrasting climate change scenarios (RCP2.6 and RCP8.5). Models predicted distributions that matched empirical distributions in cross-validation, and identified distinct refugia at the low latitude ranges, which largely coincide among species with similar ecological niches. Transferred models into the future projected polewards expansions and substantial range losses in lower latitudes, where richer gene pools are expected (in Nova Scotia and Iberia for cold affinity species and Gibraltar, Alboran, and Morocco for warm-temperate species). These effects were projected for both scenarios but were intensified under the extreme RCP8.5 scenario, with the complete borealization (circum-Arctic colonization) of kelp forests, the redistribution of the biogeographical transitional zones of the North Atlantic, and the erosion of global gene pools across all species. As the geographic distribution of genetic variability is unknown for most marine species, our results represent a baseline for identification of locations potentially rich in unique phylogeographic lineages that are also climatic relics in threat of disappearing.

Dramatic loss of seagrass habitat under projected climate change in the Mediterranean Sea.

Although climate warming is affecting most marine ecosystems, the Mediterranean is showing earlier impacts. Foundation seagrasses are already experiencing a well-documented regression in the Mediterranean which could be aggravated by climate change. Here, we forecast distributions of two seagrasses and contrast predicted loss with discrete regions identified on the basis of extant genetic diversity. Under the worst-case scenario, Posidonia oceanica might lose 75% of suitable habitat by 2050 and is at risk of functional extinction by 2100, whereas Cymodocea nodosa would lose only 46.5% in that scenario as losses are compensated with gained and stable areas in the Atlantic. Besides, we predict that erosion of present genetic diversity and vicariant processes can happen, as all Mediterranean genetic regions could decrease considerably in extension in future warming scenarios. The functional extinction of Posidonia oceanica would have important ecological impacts and may also lead to the release of the massive carbon stocks these ecosystems stored over millennia.

Cryptic diversity, geographical endemism and allopolyploidy in NE Pacific seaweeds.

BACKGROUND: Molecular markers are revealing a much more diverse and evolutionarily complex picture of marine biodiversity than previously anticipated. Cryptic and/or endemic marine species are continually being found throughout the world oceans, predominantly in inconspicuous tropical groups but also in larger, canopy-forming taxa from well studied temperate regions. Interspecific hybridization has also been found to be prevalent in many marine groups, for instance within dense congeneric assemblages, with introgressive gene-flow being the most common outcome. Here, using a congeneric phylogeographic approach, we investigated two monotypic and geographically complementary sister genera of north-east Pacific intertidal seaweeds (Hesperophycus and Pelvetiopsis), for which preliminary molecular tests revealed unexpected conflicts consistent with unrecognized cryptic diversity and hybridization. RESULTS: The three recovered mtDNA clades did not match a priori species delimitations. H. californicus was congruent, whereas widespread P. limitata encompassed two additional narrow-endemic species from California - P. arborescens (here genetically confirmed) and P. hybrida sp. nov. The congruence between the genotypic clusters and the mtDNA clades was absolute. Fixed heterozygosity was apparent in a high proportion of loci in P. limitata and P. hybrida, with genetic analyses showing that the latter was composed of both H. californicus and P. arborescens genomes. All four inferred species could be distinguished based on their general morphology. CONCLUSIONS: This study confirmed additional diversity and reticulation within NE Pacific Hesperophycus/Pelvetiopsis, including the validity of the much endangered, modern climatic relict P. arborescens, and the identification of a new, stable allopolyploid species (P. hybrida) with clearly discernable ancestry (♀ H. californicus x ♂ P. arborescens), morphology, and geographical distribution. Allopolyploid speciation is otherwise completely unknown in brown seaweeds, and its unique occurrence within this genus (P. limitata possibly representing a second example) remains enigmatic. The taxonomic separation of Hesperophycus and Pelvetiopsis is not supported and the genera should be synonymized; we retain only the latter. The transitional coastline between Point Conception and Monterey Bay represented a diversity hotspot for the genus and the likely sites of extraordinary evolutionary events of allopolyploid speciation at sympatric range contact zones. This study pinpoints how much diversity (and evolutionary processes) potentially remains undiscovered even on a conspicuous seaweed genus from the well-studied Californian intertidal shores let alone in other, less studied marine groups and regions/depths.

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.

Development and characterization of twelve microsatellite markers for Porphyra linearis Greville.

The genus Porphyra (and its sister genus Pyropia) contains important red algal species that are cultivated and/or harvested for human consumption, sustaining a billion-dollar aquaculture industry. A vast amount of research has been focused on species of this genus, including studies on genetics and genomics among other areas. Twelve novel microsatellite markers were developed here for Porphyra linearis. Markers were characterized using 32 individuals collected from four natural populations of P. linearis with total heterozygosity varying from 0.098 to 0.916. The number of alleles per locus ranged from 2 to 18. All markers showed cross amplification with Porphyra umbilicalis and/or Porphyra dioica. These polymorphic microsatellite markers are useful for investigating population genetic diversity and differentiation in P. linearis and may become useful for other genetic research on the reproductive biology of this important species.

Genetic diversity of Saccharina latissima (Phaeophyceae) along a salinity gradient in the North Sea-Baltic Sea transition zone.

The North Sea-Baltic Sea transition zone constitutes a boundary area for the kelp species Saccharina latissima due to a strong salinity gradient operating in the area. Furthermore, the existence of S. latissima there, along Danish waters, is fairly patchy as hard bottom is scarce. In this study, patterns of genetic diversity of S. latissima populations were evaluated along the salinity gradient area of Danish waters (here designated brackish) and were compared to reference sites (here designated marine) outside the gradient area, using microsatellite markers. The results showed that the S. latissima populations were structured into two clusters corresponding to brackish versus marine sites, and that gene flow was reduced both between clusters and between populations within clusters. In addition, results provided empirical evidence that marginal populations of S. latissima in the salinity gradient area exhibited a distinct genetic structure when compared to marine ones. Brackish populations were less diverse, more related, and showed increased differentiation over distance compared to marine populations. The isolation of the brackish S. latissima populations within the salinity gradient area of Danish waters in conjunction with their general low genetic diversity makes these populations vulnerable to ongoing environmental and climate change, predicted to result in declining salinity in the Baltic Sea area that may alter the future distribution and performance of S. latissima in the area.

Intraspecific genetic lineages of a marine mussel show behavioural divergence and spatial segregation over a tropical/subtropical biogeographic transition.

BACKGROUND: Intraspecific variability is seen as a central component of biodiversity. We investigated genetic differentiation, contemporary patterns of demographic connectivity and intraspecific variation of adaptive behavioural traits in two lineages of an intertidal mussel (Perna perna) across a tropical/subtropical biogeographic transition. RESULTS: Microsatellite analyses revealed clear genetic differentiation between western (temperate) and eastern (subtropical/tropical) populations, confirming divergence previously detected with mitochondrial (COI) and nuclear (ITS) markers. Gene flow between regions was predominantly east-to-west and was only moderate, with higher heterozygote deficiency where the two lineages co-occur. This can be explained by differential selection and/or oceanographic dynamics acting as a barrier to larval dispersal. Common garden experiments showed that gaping (periodic closure and opening of the shell) and attachment to the substratum differed significantly between the two lineages. Western individuals gaped more and attached less strongly to the substratum than eastern ones. CONCLUSIONS: These behavioural differences are consistent with the geographic and intertidal distributions of each lineage along sharp environmental clines, indicating their strong adaptive significance. We highlight the functional role of diversity below the species level in evolutionary trends and the need to understand this when predicting biodiversity responses to environmental change.

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.

Trade-offs between life-history traits at range-edge and central locations.

The allocation of resources to different life-history traits should represent the best compromise in fitness investment for organisms in their local environment. When resources are limiting, the investment in a specific trait must carry a cost that is expressed in trade-offs with other traits. In this study, the relative investment in the fitness-related traits, growth, reproduction and defence were compared at central and range-edge locations, using the seaweed Ascophyllum nodosum as a model system. Individual growth rates were similar at both sites, whereas edge populations showed a higher relative investment in reproduction (demonstrated by a higher reproductive allocation and extended reproductive periods) when compared to central populations that invested more in defence. These results show the capability of A. nodosum to differentially allocate resources for different traits under different habitat conditions, suggesting that reproduction and defence have different fitness values under the specific living conditions experienced at edge and central locations. However, ongoing climate change may threaten edge populations by increasing the selective pressure on specific traits, forcing these populations to lower the investment in other traits that are also potentially important for population fitness.

Reproductive strategies and isolation-by-demography in a marine clonal plant along an eutrophication gradient.

Genetic diversity in clonal organisms includes two distinct components, (i) the diversity of genotypes or clones (i.e. genotypic richness) in a population and (ii) that of the alleles (i.e. allelic and gene diversity within populations, and differentiation between populations). We investigated how population differentiation and genotypic components are associated across a gradient of eutrophication in a clonal marine plant. To that end, we combined direct measurements of sexual allocation (i.e. flower and seed counts) and genotypic analyses, which are used as an estimator of effective sexual reproduction across multiple generations. Genetic differentiation across sites was also modelled according to a hypothesis here defined as isolation-by-demography, in which we use population-specific factors, genotypic richness and eutrophication that are hypothesized to affect the source-sink dynamics and thus influence the genetic differentiation between a pair of populations. Eutrophic populations exhibited lower genotypic richness, in agreement with lower direct measurements of sexual allocation and contemporaneous gene flow. Genetic differentiation, while not explained by distance, was best predicted by genotypic richness and habitat quality. A multiple regression model using these two predictors was considered the best model (R(2) = 0.43). In this study, the relationship between environment and effective sexual-asexual balance is not simply (linearly) predicted by direct measurements of sexual allocation. Our results indicate that population-specific factors and the isolation-by-demography model should be used more often to understand genetic differentiation.

Disentangling the Influence of Mutation and Migration in Clonal Seagrasses Using the Genetic Diversity Spectrum for Microsatellites.

The recurrent lack of isolation by distance reported at regional scale in seagrass species was recently suggested to stem from stochastic events of large-scale dispersal. We explored the usefulness of phylogenetic information contained in microsatellite loci to test this hypothesis by using the Genetic Diversity Spectrum (GDS) on databases containing, respectively, 7 and 9 microsatellites genotypes for 1541 sampling units of Posidonia oceanica and 1647 of Cymodocea nodosa The simultaneous increase of microsatellite and geographic distances that emerges reveals a coherent pattern of isolation by distance in contrast to the chaotic pattern previously described using allele frequencies, in particular, for the long-lived P. oceanica These results suggest that the lack of isolation by distance, rather than the resulting from rare events of large-scale dispersal, reflects at least for some species a stronger influence of mutation over migration at the scale of the distribution range. The global distribution of genetic polymorphism may, therefore, result predominantly from ancient events of step-by-step (re)colonization followed by local recruitment and clonal growth, rather than contemporary gene flow. The analysis of GDS appears useful to unravel the evolutionary forces influencing the dynamics and evolution at distinct temporal and spatial scales by accounting for phylogenetic information borne by microsatellites, under an appropriate mutation model. This finding adds nuance to the generalization of the influence of large-scale dispersal on the dynamics of seagrasses.

Shift happens: trailing edge contraction associated with recent warming trends threatens a distinct genetic lineage in the marine macroalga Fucus vesiculosus.

BACKGROUND: Significant effects of recent global climate change have already been observed in a variety of ecosystems, with evidence for shifts in species ranges, but rarely have such consequences been related to the changes in the species genetic pool. The stretch of Atlantic coast between North Africa and North Iberia is ideal for studying the relationship between species distribution and climate change as it includes the distributional limits of a considerable number of both cold- and warm-water species.We compared temporal changes in distribution of the canopy-forming alga Fucus vesiculosus with historical sea surface temperature (SST) patterns to draw links between range shifts and contemporary climate change. Moreover, we genetically characterized with microsatellite markers previously sampled extinct and extant populations in order to estimate resulting cryptic genetic erosion. RESULTS: Over the past 30 years, a geographic contraction of the southern range edge of this species has occurred, with a northward latitudinal shift of approximately 1,250 km. Additionally, a more restricted distributional decline was recorded in the Bay of Biscay. Coastal SST warming data over the last three decades revealed a significant increase in temperature along most of the studied coastline, averaging 0.214°C/decade. Importantly, the analysis of existing and extinct population samples clearly distinguished two genetically different groups, a northern and a southern clade. Because of the range contraction, the southern group is currently represented by very few extant populations. This southern edge range shift is thus causing the loss of a distinct component of the species genetic background. CONCLUSIONS: We reveal a climate-correlated diversity loss below the species level, a process that could render the species more vulnerable to future environmental changes and affect its evolutionary potential. This is a remarkable case of genetic uniqueness of a vanishing cryptic genetic clade (southern clade).

Range map data of marine ecosystem structuring species under global climate change.

Data on contemporary and future geographical distributions of marine species are crucial for guiding conservation and management policies in face of climate change. However, available distributional patterns have overlooked key ecosystem structuring species, despite their numerous ecological and socioeconomic services. Future range estimates are mostly available for few species at regional scales, and often rely on the outdated Representative Concentration Pathway scenarios of climate change, hindering global biodiversity estimates within the framework of current international climate policies. Here, we provide range maps for 980 marine structuring species of seagrasses, kelps, fucoids, and cold-water corals under present-day conditions (from 2010 to 2020) and future scenarios (from 2090 to 2100) spanning from low carbon emission scenarios aligned with the goals of the Paris Agreement (Shared Socioeconomic Pathway 1-1.9), to higher emissions under reduced mitigation strategies (SSP3-7.0 and SSP5-8.5). These models were developed using state-of-the-art and advanced machine learning algorithms linking the most comprehensive and quality-controlled datasets of occurrence records with high-resolution, biologically relevant predictor variables. By integrating the best aspects of species distribution modelling over key ecosystem structuring species, our datasets hold the potential to enhance the ability to inform strategic and effective conservation policy, ultimately supporting the resilience of ocean ecosystems.

Cryptic diversity in southern African kelp.

The southern coast of Africa is one of the few places in the world where water temperatures are predicted to cool in the future. This endemism-rich coastline is home to two sister species of kelps of the genus Ecklonia maxima and Ecklonia radiata, each associated with specific thermal niches, and occuring primarily on opposite sides of the southern tip of Africa. Historical distribution records indicate that E. maxima has recently shifted its distribution ~ 70 km eastward, to sites where only E. radiata was previously reported. The contact of sister species with contrasting thermal affinities and the occurrence of mixed morphologies raised the hypothesis that hybridization might be occurring in this contact zone. Here we describe the genetic structure of the genus Ecklonia along the southern coast of Africa and investigate potential hybridization and cryptic diversity using a combination of nuclear microsatellites and mitochondrial markers. We found that both species have geographically discrete genetic clusters, consistent with expected phylogeographic breaks along this coastline. In addition, depth-isolated populations were found to harbor unique genetic diversity, including a third Ecklonia lineage. Mito-nuclear discordance and high genetic divergence in the contact zones suggest multiple hybridization events between Ecklonia species. Discordance between morphological and molecular identification suggests the potential influence of abiotic factors leading to convergent phenotypes in the contact zones. Our results highlight an example of cryptic diversity and hybridization driven by contact between two closely related keystone species with contrasting thermal affinities.

Looking into the black box: simulating the role of self-fertilization and mortality in the genetic structure of Macrocystis pyrifera.

Patterns of spatial genetic structure (SGS), typically estimated by genotyping adults, integrate migration over multiple generations and measure the effective gene flow of populations. SGS results can be compared with direct ecological studies of dispersal or mating system to gain additional insights. When mismatches occur, simulations can be used to illuminate the causes of these mismatches. Here, we report a SGS and simulation-based study of self-fertilization in Macrocystis pyrifera, the giant kelp. We found that SGS is weaker than expected in M. pyrifera and used computer simulations to identify selfing and early mortality rates for which the individual heterozygosity distribution fits that of the observed data. Only one (of three) population showed both elevated kinship in the smallest distance class and a significant negative slope between kinship and geographical distance. All simulations had poor fit to the observed data unless mortality due to inbreeding depression was imposed. This mortality could only be imposed for selfing, as these were the only simulations to show an excess of homozygous individuals relative to the observed data. Thus, the expected data consistently achieved nonsignificant differences from the observed data only under models of selfing with mortality, with best fits between 32% and 42% selfing. Inbreeding depression ranged from 0.70 to 0.73. The results suggest that density-dependent mortality of early life stages is a significant force in structuring Macrocystis populations, with few highly homozygous individuals surviving. The success of these results should help to validate simulation approaches even in data-poor systems, as a means to estimate otherwise difficult-to-measure life cycle parameters.