Trans-Arctic asymmetries, melting pots and weak species cohesion in the low-dispersal amphiboreal seaweed Fucus distichus

Abstract

Amphiboreal taxa are often composed of vicariant phylogroups and species complexes whose divergence and phylogeographic affinities reflect a shared history of chronic isolation and episodic trans-Arctic dispersal. Ecological filters and shifting selective pressures may also promote selective sweeps, niche shifts and ecological speciation during colonization, but these are seldom considered at biogeographical scales. Here we integrate genetic data and Ecologic Niche Models (ENMs) to investigate the historical biogeography and cohesion of the polymorphic rockweed Fucus distichus throughout its immense amphiboreal range, focusing on trans-Arctic asymmetries, glacial/interglacial dynamics, and integrity of sympatric eco-morphotypes. Populations were sampled throughout the Pacific and the Atlantic, from southern rear-edges to the high-Arctic. They were genotyped for seven microsatellites and an mtDNA spacer, and genetic diversity and structure were assessed from global to local scales. ENMs were used to compare niche divergence and magnitude of post-glacial range shifts in Pacific versus Atlantic sub-ranges. Haplotypic and genotypic data revealed distinct and seemingly isolated Pacific vs Arctic/Atlantic gene-pools, with finer-scale regional sub-structuring pervasive in the Pacific. MtDNA diversity was highly structured and overwhelmingly concentrated in the Pacific. Regionally, Alaska showed the highest intra-population diversity but the lowest levels of endemism. Some sympatric/parapatric ecotypes exhibited distinct genotypic/haplotypic compositions. Strikingly, niche models revealed higher Pacific tolerance to maximum temperatures and predicted a much more consolidated presence in the NE Atlantic. Glacial and modern ranges overlapped extensively in the Pacific, whereas the modern Atlantic range was largely glaciated or emerged during the Last Glacial Maximum. Higher genetic and ecogeographic diversity supports a primary Pacific diversification and secondary Atlantic colonization, also likely reflecting the much larger and more stable climatic refugia in the Pacific. The relic distribution and reduced ecological/morphological plasticity in the NE Atlantic are hypothesized to reflect functional trans-Arctic bottlenecks, recent colonization or competition with congeners. Within the Pacific, Alaska showed signatures of a post-glacial melting pot of eastern and southern populations. Genetic/ecotypic variation was generally not sufficiently discontinuous or consistent to justify recognizing multiple taxonomic entities, but support a separate species in the eastern Pacific, at the southern rear-edge. We predict that layered patterns of phylogeographic structure, incipient speciation and niche differences might be common among widespread low-dispersal amphiboreal taxa.

Introduction

Following the first opening of the Bering Strait at the end of the Miocene some 5 Ma ago, trans-Arctic interchanges have been regularly occurring between Pacific and Atlantic biotas (Lindstrom, 2001; Dodson et al., 2007; Nikula et al., 2007; Cánovas et al., 2011; Hardy et al., 2011; Saunders and McDevit, 2013). Throughout the long Pleistocene glaciations (starting 2.5 Ma), however, the Arctic basin and adjacent continents and shorelines were covered with perennial ice-sheets, ice-shelves and sea-ice, making them uninhabitable to coastal species (Dunton, 1992; Kaufman et al., 2004; Laakkonen et al., 2013; Neiva et al., 2018). Because of lowered sea-levels (down to −130 m; Pico et al., 2020), the Bering Strait (currently at −50 m) was also a terrestrial bridge connecting the Eurasian and American continents (Beringia). During interglacial periods, marine transgressions and seasonally ice-free marine routes briefly resume Pacific/Atlantic oceanographic connections (Dunton, 1992), but for most organisms and for most of the time, the Arctic environment remains a formidable dispersal barrier. Accordingly, most amphiboreal taxa include vicariant phylogroups and species complexes whose systematic and biogeographic affinities reflect such shared history of episodic dispersal against a general background of trans-Arctic isolation (Carr et al., 2011; Halanych et al., 2013; Neiva et al., 2018).

Past the Bering Strait and along marine Arctic routes, key traits (e.g. dispersal ecology, ability to survive extreme Arctic conditions), together with sweepstake events and historical processes, ultimately determine the pool and colonization dynamics of migrants (Hardy et al., 2011). At infra-specific levels, extreme founder events, genetic surfing, priority effects and selective filters are also predicted to result in important trans-Arctic genetic bottlenecks, but these have seldom been assessed. In the absence of detailed fossil records, species origins have been established under the assumption that older ranges contain higher and more deeply structured genetic diversity than recently colonized ranges (Wares and Cunningham, 2001; Bringloe and Saunders, 2018). Differences need careful interpretation, however, as they may also reflect long-term differences in habitat availability and connectivity (e.g. larger climatic refugia, sensu Neiva et al., 2014) between ocean basins.

Many low-dispersal marine species, notably macroalgae, show important intra-specific geographic trait variation, for instance in physiological tolerance to thermal stress (Mohring et al., 2014; Saada et al., 2016; Smolina et al., 2016; King et al., 2018; Liesner et al., 2020), that should also not be overlooked in historical contexts of colonization. If a locally-adapted regional population is the sole source in a long-distance colonization pulse, realized niche/distribution in the establishing range could potentially be a fraction of that in the original range (Holt, 2009). In the longer-term, chronic isolation and different selective pressures and environments may also promote niche shifts through allopatrically evolved environmental tolerances (e.g. thermal adaptation, fundamental niche shifts) or release/introduction to new biotic constraints (e.g. competition or predation, realized niche shifts) (Wiens and Graham, 2005; Mitchell et al., 2006). In a range of invertebrate taxa, correlative models show that (realized) niche differences between eastern Atlantic (putative source) and western Atlantic (putative post-glacial colonization) populations can be negligible (Waltari and Hickerson, 2013). It remains unassessed, however, whether niche conservatism holds for older vicariant taxa, or for highly structured marine species whose “global” niche contains distinct evolutionary lineages and multiple specialized niches among local populations.

Glacial–interglacial cycles also fuel important range shifts and dissections within oceanic basins. At peripheral range margins, populations and their genetic variation are periodically erased during glacial (high-latitude) and warmer interglacial (low-latitude) extremes. These are often associated with asynchronous range expansions (e.g. post-glacial northward colonizations) at opposite distributional boundaries (Assis et al., 2018a; Neiva et al., 2016). At more temperate latitudes, disconnected refugia may allow vicariant subpopulations to persist in the long-term, often producing divergent phylogeographic lineages. During climatic shifts, lineage redistributions lead to secondary contacts, with subsequent outcomes ranging from widespread admixture to stable genetic sectors depending on species dispersal potential and population dynamics. Among low dispersers, genetic breaks can be surprisingly sharp and stable, particularly when large habitat or oceanographic discontinuities, selfing, and strong density-barrier effects are involved (Fraser et al., 2009; Neiva et al., 2012a, 2016, 2018; Assis et al., 2016, 2018b). The stable co-existence of different lineages, particularly when matching distinct morphologic and ecotypic entities, often indicates partial to complete isolation (Tellier et al., 2011; Montecinos et al., 2012; Neiva et al., 2018).

The rockweed Fucus distichus Linnaeus 1767 (Fucales, Ochrophyta) is one of the most widespread, polymorphic and ecologically plastic members of the northern hemisphere amphiboreal marine biota (Powell, 1963; Kucera and Saunders, 2008). This seaweed is the only member of the family Fucaceae occurring on both sides of the Pacific and the Atlantic, and the one extending deepest into the Arctic. Its wide distribution and ecotypic variation have led to the use of multiple taxonomic entities, particularly in the Atlantic (reviewed in Powell, 1957, Powell, 1963). Despite morphological plasticity and even broader range of coastal and estuarine habitats (Kucera and Saunders, 2008; Neiva et al., 2012b), all forms in the NE Pacific were lumped into a single geographical taxon, F. gardneri P.C. Silva 1953. Molecular work has established the monophyly of F. distichus sensu lato and its affinity to F. serratus, but so far has consistently failed to resolve any ecotypic/geographic entity as separate species (Serrão et al., 1999; Coyer et al., 2006, 2011b; Kucera and Saunders, 2008; Cánovas et al., 2011; Laughinghouse et al., 2015). Some degree of phylogeographical structure was apparent though, e.g. discriminating Pacific vs Arctic/Atlantic ranges, with higher diversity suggesting a Pacific origin for the taxa.

Fucus distichus is regarded as a cold-adapted taxon, particularly in the Atlantic. The species undeniably thrives throughout Arctic latitudes, withstanding long periods of freezing low-light conditions and intense ice-scouring (Laughinghouse et al., 2015; Jueterbock et al., 2016). Its distribution however can extend considerably towards more temperate latitudes. In the NW Atlantic it reaches Massachusetts (Ester Serrão, personal observations) and in the NW Pacific it reaches the northern parts of the Sea of Japan (Yamada, 1980; Chernova, 2016), both around ca. 41–42° Lat. By comparison, in the NE Atlantic F. distichus (var. anceps only) occurs as small disconnected populations restricted to very exposed sites only as far south as western Ireland (ca. 52° Lat., Powell, 1957). In the NE Pacific, F. distichus occurs much farther south, almost reaching Point Conception in Central California (ca. 34.5° Lat., present study), suggesting potential niche(s) differences (or niche filling) across disparate oceanic regions.

Incomplete geographic sampling and/or low marker resolution (e.g. IGS/cox1 SSCP genotyping, Coyer et al., 2011b) have so far hindered a full appreciation of the global diversity within the F. distichus complex. The NE Pacific coast, particularly California, has never been adequately surveyed, but is hypothesized to represent a genetically-unique, rear-edge climatic refugium as in related Pelvetiopsis for which it is a center of relic biodiversity and evolutionary novelty (Neiva et al., 2017; Sousa et al., 2019). Regional phylogeographic breaks and “melting pots” resulting from post-glacial extinction/expansion dynamics and secondary contacts are also expected to have formed throughout the North Pacific and potentially in the high Arctic, but have never been adequately investigated. Finally, the co-occurrence of distinct eco-morphotypes in many shorelines (e.g. Pearson and Brawley, 1996) raises inevitable but still unanswered questions regarding potential isolating barriers and incipient speciation.

This study uses phylogeographic (mtDNA, microsatellites) and niche-modelling approaches to investigate the biogeographical history, phylogeographic structure and cohesion of polymorphic F. distichus at multiple spatial scales throughout its immense amphiboreal range. A range of related topics are examined, including 1) trans-Arctic dispersal/vicariance history, including putative origins and modern contacts; 2) potential trans-Arctic genetic/eco-geographic asymmetries; 3) signatures of potential southern hot-spots vs high-latitude melting pots; and 4) the genetic integrity of local (sympatric/parapatric) eco-types.

Materials and methods

Global phylogeography and population structure

Populations of F. distichus were sampled on both sides of the Pacific and the Arctic/Atlantic, including the understudied California rear-edge, the trans-Pacific Aleutian Arc, and the Greenland coast. At each site, tissue samples (16 < n < 24) were, whenever possible, collected every other meter along linear transects. Multiple collections were obtained at specific sites (Maine, California, Alaska) where locally distinct eco-morphotypes and/or discontinuous tidal heights were apparent. Genomic DNA was extracted from silica-dried tissue using the Nucleospin 96 Plant II kit (Macherey-Nagel, Germany). Twenty microsatellite loci previously developed for other Fucus species (Engel et al., 2003; Wallace et al., 2004; Perrin et al., 2007; Coyer et al., 2009) were amplified in a geographically diverse panel of populations to assess cross-amplification success and levels of polymorphism (data not shown). The seven most polymorphic microsatellites were selected and used to generate multi-locus genotypes for all individuals of F. distichus (see Supplementary Table S1 for details). The 23S/trnK mitochondrial intergenic spacer (mtIGS, a marker with good phylogeographic resolution in Fucus; Neiva et al., 2010, 2014; Coyer et al., 2011a; Laughinghouse et al., 2015), was screened in the same populations. Amplified fragments were run in an ABI PRISM 3130xl automated capillary sequencer (Applied Biosystems) at CCMAR, Portugal. Microsatellite alleles were manually scored in STRAND (Veterinary Genetics Laboratory, University of California, Davis; http://www.vgl.ucdavis.edu/STRand) using the 500 LIZ size standard (Applied Biosystems). Sequences were aligned and edited with GENEIOUS 4.8 (Biomatters; http://www.geneious.com).

For microsatellite data, summary statistics of genetic diversity (allele frequency, mean allelic richness [A], unbiased gene diversity [H_E], observed heterozygosity [H_O] and inbreeding coefficients [F_IS]) were calculated with GENETIX 4.05 (Laboratoire Génome, Populations, Interactions, Université de Montpellier II; http://kimura.univ-montp2.fr/genetix) and the R package standArich v1.0 (available at http://alberto-lab.blogspot.pt/p/code.html). Pairwise population differentiation (Jost’s D, Weir and Cockerham’s Θ) was calculated with the R package DiveRsity 1.9 (Keenan et al., 2013). Genetic structure was analyzed at multiple spatial scales with STRUCTURE 2.3 (Pritchard Lab, Stanford University; http://pritchardlab.stanford.edu/structure.html) without prior population assignments. A range of assumed population cluster size (K, set sequentially from 1 to 8–12, depending on the analyses) was run 5 times using a burn-in of 500 000 iterations and a run-length of 1 × 10⁶ iterations, using the ΔK criterion (Evanno et al., 2005) to assist selection of “best” Ks. Structure analyses were complemented with factorial correspondence analyses (FCA) implemented in GENETIX.

The genealogic relationships among mtIGS haplotypes were inferred with Network 5 (Fluxus-Engineering; http://www.fluxus-engineering.com/sharenet.htm) using the median-joining algorithm, and their geographic distribution and frequency mapped. Haplotype (H_hap)and nucleotide (Θ_π) diversities were estimated for each population and selected geographical regions using DNASP 6 (Julio Rozas Lab, Universitat de Barcelona; http://www.ub.edu/dnasp/), the former coding indels as single mutations and the latter discarding indels in pair-wise comparisons only.

Ecological niche modelling

Ecological niche models (ENMs) were developed with Boosted Regression Trees (BRT; Elith et al., 2008), a high-performance machine learning algorithm that automatically fits complex interactions and nonlinear relationships between predictor variables, while reducing the potential of overfitting by forcing monotonicity responses (Hofner et al., 2011) and by optimal parametrization (Elith et al., 2008). The algorithm fitted georeferenced records of occurrence against biologically meaningful predictor variables for macroalgae species (Fragkopoulou et al., 2022), namely, nutrients (as nitrate and phosphate), salinity, sea ice thickness and maximum winter and summer sea surface temperatures.

Occurrence records of F. distichus were collated from the fine-tuned dataset of marine forests (Assis et al., 2020) which integrates data from different sources (e.g. databases, literature and herbaria) and uses a flagging system to exclude potentially erroneous entries (e.g. records outside known biogeographical ranges). Because absence records were unavailable at the scale of the study, pseudo-absences were randomly generated using a 1:1 ratio with presence records, in locations with no presences recorded (Barbet-Massin et al., 2012) and limited to the biogeographical provinces (Spalding et al., 2007) where the species is known to occur, as well as their neighboring provinces. This allowed restricting pseudo-absences to regions where dispersal can potentially occur, a crucial step of ENM (Martins et al., 2021). The negative effects of spatial autocorrelation and sampling bias (e.g. regional differences in density of records) were further reduced by randomly choosing one record only from those at distances considered significantly autocorrelated, as estimated by a correlogram testing the spatial variability of predictors as a function of geographic distance (Boavida et al., 2016). Predictor variables describing sea surface conditions were obtained from Bio-ORACLE v2.1 (Assis et al., 2018c; Tyberghein et al., 2012) for near present-day conditions and for the Last Glacial Maximum (LGM).

A cross-validation (CV) framework using 10-fold latitudinal bands was implemented to tune the optimal parameter combination of BRT and to estimate the performance of the models and their potential transferability in independent data (Assis et al., 2018a; Krause-Jensen et al., 2020; Fragkopoulou et al., 2021). The “grid search” statistical method interactively tested all possible numbers of trees (range 50–1000, step 50), tree complexity (range 1–6, step 1) and learning rate (0.01, 0.005 and 0.001), by fitting occurrence records with all but one latitudinal band at a time. At each step, performance was estimated with the area under the curve (AUC). The best combination of BRT parameters was identified as that resulting in models with higher AUC in CV (Assis et al., 2018a; Krause-Jensen et al., 2020; Fragkopoulou et al., 2021). Overfitting was further controlled by forcing monotonic responses based on the expected biological effect of predictors (positive or negative) on the probability of occurrence (Hofner et al., 2011). Maximum temperatures and ice fitted negative monotonic responses, whereas nutrients and salinity fitted positive responses (e.g. Assis et al., 2018a; Fragkopoulou et al., 2021; Martins et al., 2021; Krause-Jensen et al., 2020).

ENM were performed using two approaches: (1) species modelling considering the entire amphiboreal range and (2) sub-taxon modelling (i.e. within-taxon niche structure; Assis et al., 2016), using distinct models for the Pacific and Arctic/Atlantic ranges where differences in evolutionary history, genetic composition and diversity were apparent (see results and discussion sections). Habitat suitability maps for the present and the LGM were produced with a BRT model using the optimal parameters found in CV. These maps were further reclassified to presence and absence (binomial responses) using a threshold maximizing specificity (true negative rate) and sensitivity (true positive rate; Jiménez-Valverde and Lobo, 2007).

Niche overlap between the Pacific and Arctic/Atlantic ranges was estimated with the Bayesian framework of Swanson et al. (2015). This defines the niche of each range as the 95% probability region of the multivariate environmental space defined by ENM, and determined overlap as the probability of an individual from the Pacific being found in the niche of Arctic/Atlantic, and vice-versa. Additionally, physiological tolerance limits for the Pacific and Arctic/Atlantic ranges were captured from partial dependence functions (Elith et al., 2008).

Results

Genetic analyses

A total of 31 collections of Fucus distichus (590 individuals) were sampled from 16 Pacific and 10 Arctic/Atlantic sites (Table 1), covering most of the species range in the northern hemisphere. Collections from Alaska (Pirates, Kayak) and Oregon (Bob Creek) contained a few F. spiralis L. that were easily recognizable genetically and that were excluded from analyses. Microsatellite polymorphism (f_allele > 0.005) ranged between 2 and 13 alleles per locus, averaging 6.14. A total of 38 mtIGS haplotypes were detected among the 470 individuals sequenced (Table 1; Figure 1; Supplementary Table S2). A subset of 16 haplotypes (≈42%) accounted for almost 90% of the sequences, while another 14 were sampled only once. Several indels of variable size were apparent (Figure 1). Genetic diversity within populations was typically low (A ≤ 3.2; H_E ≤ 0.41, mean n_hap ≈ 2.15), as most populations were monomorphic (or nearly fixed) at several loci (Supplementary Figure S1, Supplementary Table S2). Nearly all populations exhibited significant heterozygote deficits. Avacha Bay (pop. 23), the only populations showing fixed heterozygosity at locus Fg25 (and hence anomalous positive F_IS values), also exhibited significant heterozygote deficit when the putative duplicated locus was removed from analyses (data not shown).

Figure 1

F. distichus exhibited strong phylogeographic structure. MtIGS haplotypes (Genbank accession numbers PP358930–PP359399) formed a complex network of shallow haplotypic radiations corresponding to at least 10 major genetic/geographic discontinuities (Figures 1, 2). Eight of these were endemic to the Pacific and 2 were endemic to the Arctic/Atlantic, i.e. they were never found together. Atlantic haplogroups were genealogically close and they were not more differentiated from Pacific haplogroups than different Pacific haplogroups were between them. Haplogroup DO was widespread throughout both sides of the Atlantic and in the Arctic (encompassing samples originally assigned to F. distichus s.s., F. anceps, F. edentatus and F. evanescens), whereas haplogroup LO was restricted to the cold-temperate NW Atlantic (Maine and Newfoundland, Figure 2A). In Maine (populations 4 and 5), these haplogroups matched two regionally parapatric eco-types – LO corresponding to smaller forms from high-intertidal rock-pools (F. distichus s.s.-like) and DO to larger forms from the low intertidal (F. evanescens-like, Figure 2B, top right). In the NW Pacific and western Aleutians (Adak, pop. 22) only the P, DY (endemic) and Y (shared with Alaska) haplogroups were sampled. All 9 haplotypes detected there were endemic to the region and the majority was restricted to specific populations (Table 1; Supplementary Table S2). The remaining 5 haplogroups were restricted to the NE Pacific (Figure 2A). The LB, DB and DG haplogroups were present in the Mendocinian and Alaska regions, extending into the eastern Aleutians (Unalaska, pop. 21). In the Mendocinian region each haplogroup dominated successively short stretches of coastline, whereas in SE Alaska they co-occurred more extensively (together with the Y haplogroup) at local scales. Remarkably, haplotypic composition among Alaskan populations varied considerably at very local scales, such as among High- vs Low-shore collections at Pirates (Figure 2B, bottom centre). Finally, the R and LG haplogroups were endemic to the Montereyan region. In Pigeon Point (pop. 11), these haplogroups matched a priori “large” and “small” morphological assignments, respectively.

Figure 2

Conservation of Arctic Flora and Fauna (2010), and the bold black bars the known southern range limits. (B) Major microsatellite-based subdivisions in the northern hemisphere (top, K = 2) and in three separate regions of the north Pacific (bottom) according to Structure. Each vertical bar represents the proportions of individual multilocus genotypes assigned to each of K virtual clusters, as illustrated by the different colours, with the corresponding haplogroup on top. Population codes as in Table 1. H and L stand for High- and Low-shore bands, S and L for Small and Large ecotype, D and E for “distichus” and “edentatus” haplotypes. Standardized allelic diversity (A) and nucleotide diversity (Θ_π = pi) are shown below the Structure plots for four selected regions.

Despite the relatively small global allelic diversity, the Structure analyses separated Arctic/Atlantic and Pacific genotypic clusters (Figure 2B, top) with good support (Δk = 2922; Supplementary Figure S2). This high-level subdivision was also recovered by the FCA (Supplementary Figure S3). Within the Pacific there was no support for large-scale regional subdivisions (maximum Δk < 30; Supplementary Figure S2), to some extent due to low microsatellite allelic diversity and resolution. Finer-scale regional sub-structuring was evident nevertheless when different regions were analyzed separately (Figure 2B, bottom). In the NW Pacific and W Aleutians, where most populations exhibited fixed allelic differences in at least one locus (Supplementary Figure S1), recovered genotypic clusters (K = 6, Δk = 254; Supplementary Figure S2) matched closely single populations (Figure 2B, bottom left). The same applied from California to British Columbia, where genotypic clusters (K = 6, Δk = 370; Supplementary Figure S2) again matched to a large degree specific populations or contiguous population sets (Figure 2B, bottom right). Despite local sympatry, “large” and “small” morphotypes from Pigeon Point, characterized by different R and LG haplotypes (see above), also matched distinct genotypic groups. Conversely, in SE Alaska, support for 3-5 genotypic clusters was weak (60 < Δk <90, Supplementary Figure S2), and these did not consistently match distinct haplogroups, shore-heights, or regional populations (Figure 2B, bottom centre). Structure analyses also failed to resolve the two eco-morphotypes/haplogroups from Maine (or any other set of Atlantic populations), even when analyzed alone (Δk < 100, Supplementary Figure S2).

Genetic diversity was much higher in the Pacific than in the Atlantic for both microsatellite (A₍₂₃₄₎ = 6.96 vs 3.86) and haplotypic (n_hap = 29 vs 9, Θ_π × 10⁵ = 857 vs 311) data (Table 1). The Atlantic was also less diverse genetically than any of the Pacific regions considered (Figures 3A, C; Supplementary Table S3). The region between Point Conception (California) and Vancouver (British Columbia), encompassing the Montereyan and Mendocinian biogeographical regions, was (marginally) the most diverse region. Notwithstanding the occurrence of multiple, some quite divergent (e.g. R phylogroup) genetic groups, intra-population diversities were low. Pairwise differentiation, on the other hand, was often large, highlighting the depth of fine-scale spatial structure. Regional diversity in SE Alaska was not obviously larger than in the NW Pacific and western Aleutians, but consistently contained the most diverse populations in the Pacific (Figures 3A–C). Diversity was estimated for high and low shore heights at some sites (Pirates, Point Louise), thus underestimating intra-population (=site) diversity. In any case, small D values indicated relatively weak pair-wise differentiation between these diverse populations (Figure 3D). Haplotypic diversity was also large in this region, despite lacking unique haplogroups.

Figure 3

Ecological niche modelling

We collated 1283 unique occurrence records for the Pacific and 1015 for the Arctic/Atlantic. The ENMs retrieved high accuracy scores in independent data (i.e. in cross-validation procedure). Specifically, the model developed for the entire amphiboreal range achieved an average AUC of 0.829 ± 0.077, while the models for the Pacific and Arctic/Atlantic achieved an average AUC of 0.829 ± 0.077 and 0.848 ± 0.13, respectively. The posterior distribution of the probabilistic niche overlap metric between the Atlantic and Pacific groups was 96.24%, while the probability between the Pacific and Atlantic groups was 93.28% (Supplementary Figure S4). Despite the large niche overlap between ranges, physiological limiting points inferred from partial dependency functions differed between groups (Supplementary Table S4), with the niche of the Pacific showing greater thermal tolerance to maximum winter and summer temperatures (1.61°C and 1.86°C, respectively; Supplementary Table S4; Figure 4).

Figure 4

In the Pacific, ENMs over-estimated the distribution of F. distichus beyond known range limits in Point Conception (California) and Peter the Great Bay (Russia’s Far East), and polewards beyond the Bering Strait into the Chukchi and Beaufort Seas (Figure 5B). In the Arctic/Atlantic, the occurrence of the species was under-estimated in the Canadian Arctic in all models (Figure 5B; Supplementary Figure 5). The major mismatch however was in the NE Atlantic, where Ensemble and Pacific ENMs predicted a much more consolidated presence in the British Isles and contiguous French regions of Brittany and Normandy, as well as an isolated pocket in NW Iberia (Figure 5B; Supplementary Figure S5B). The Atlantic ENM over-predicted the occurrence of F. distichus only as far south as SW England (Cornwall; Supplementary Figure S5A). Hindcast projections supported a massive latitudinal shift in the potential distribution of F. distichus in the Atlantic between the LGM and today (Figure 5A). Suitable climatic conditions in the NW Atlantic would have been restricted to ice-free regions off Virginia to New Jersey and along the then-emerged Grand Banks of Newfoundland. In the NE Atlantic, projections indicated a much more southern and compressed distribution, between western Iberia and the periglacial palaeo-shorelines of the Celtic Sea and western British Isles. Globally, overlap between estimated local LGM and observed present-day distributions was quite restricted, as most of the modern range of F. distichus was either glaciated, emerged, or both during the LGM. According to models, latitudinal shifts were much less extensive in the North Pacific. During full glacial conditions, F. distichus could have spread from the Korean peninsula and Japan in the western Pacific, and Baja California (Mexico) in the NE Pacific, as far north as the paleo-shorelines of the Beringia Land Bridge. This massive range was probably interrupted to some extent between the eastern Aleutians and Vancouver Island (British Columbia), where the Cordilleran ice-sheet reached the coast. Globally, and unlike in the Atlantic, F. distichus still occurs throughout most of its immense glacial range. On the other hand, post-glacial poleward expansions, e.g. past the Bering Strait into the Arctic, have been much more modest.

Figure 5

Discussion

This study represents the most comprehensive phylogeographic survey to date of a widespread amphiboreal taxon. The extensive geographical sampling on both sides of the Pacific and the Arctic/Atlantic, extending to potential climatic refugia in all temperate regions, and the integration of molecular data (complete mtIGS spacer and 7 microsatellites) with niche modelling, revealed remarkable levels of genetic differentiation at all spatial scales examined and uncovered a multitude of historical and ongoing processes contributing to phylogeographic structure in this intertidal alga.

Biogeographic history and trans-Arctic migrations

Globally, F. distichus encompasses multiple phylogeographic lineages and microsatellite genotypic clusters within which shallower genetic discontinuities were also pervasive at more local scales. Nested patterns of structuration are common among fucoid seaweeds and other low-dispersal marine organisms (e.g. kelp, direct-developing invertebrates) and typically result from complex, non-equilibrium demographic/genetic processes (e.g. contraction/expansion and isolation/secondary contact dynamics) most commonly associated with Pleistocene climatic shifts (Marko et al., 2010; Neiva et al., 2016; Assis et al., 2018b; Grant and Chenoweth, 2018, 2021; Grant and Bringloe, 2020). More striking is the unprecedented level of genetic diversity uncovered. F. distichus produced the most complex intra-specific mtIGS network recovered so far among fucoid seaweeds (Neiva et al., 2010, 2014, 2017), with haplotypic diversity comparable to that of an entire Atlantic radiation comprising 5+ species, including F. vesiculosus, F. spiralis, F. guiryi s.l (i.e. F. limitaneus and F. macroguiryi; see Almeida et al., 2022), F. virsoides and a mtDNA-introgressed lineage of F. ceranoides (Supplementary Figure S6). Fucus distichus shares its diplontic life-cycle, hermaphroditic mating system, intertidal habitats and dispersal ecology with most of these seaweeds, and was similarly severely impacted by Pleistocene climatic shifts. The exceptional phylogeographic diversity of F. distichus therefore seems to be best explained by its distinctive amphiboreal range, as large stable ranges (particularly in the Pacific, see below) offer more opportunities for allopatric divergence and lineage persistence. The greater phylogeographic diversity of circum-boreal Saccharina latissima s.l. vs NW Pacific endemic S. japonica supports this prediction (Zhang et al., 2015; Neiva et al., 2018), and additional studies may confirm the generality of the pattern.

Phylogeographic and genotypic diversities were extremely skewed towards the Pacific, supporting the view that F. distichus diversified primarily in the Pacific and that its Arctic/Atlantic range was colonized more recently (e.g. Coyer et al., 2011b). Different haplogroup compositions indicate nevertheless that Atlantic colonization pre-dates the last (de)glaciation, which is consistent with patterns of (nuclear) trans-Arctic divergence (Cánovas et al., 2011). It should be noted that Pacific and Arctic/Atlantic haplogroups were never found co-occurring throughout the Arctic basin, including along known dispersal pathways in the Canadian Arctic and eastern Greenland (Bringloe and Saunders, 2018; Neiva et al., 2018). Unlike the Atlantic side, where F. distichus post-glacial expansions reached the Canadian Arctic and Russia’s Novaya Zemlya archipelago (e.g. Lobus and Udalov, 2021, the latter not genetically confirmed), lack of reports suggest a much more limited post-glacial Pacific expansion past the Bering Strait into the adjacent Laptev, East Siberian and Beauford seas (Wulff et al., 2009; Wilce and Dunton, 2014). Limited availability of intertidal hard substrata, low summer temperatures (<5°C) and reduced salinity concentrations in some Arctic seas may help explain this asymmetry (Jueterbock et al., 2016), but the limiting factors (and actual extent) of post-glacial Arctic colonization from the Pacific should be further clarified.

It is interesting to note that most genera of brown macroalgae with amphiboreal representatives, including several large brown (e.g. Chorda, Alaria and Saccharina) and Arctic-endemic red seaweeds, are better represented (i.e. have more species) in the Pacific than in the Atlantic (Lüning, 1990; Lane et al., 2007; Sasaki and Kawai, 2007; Adey et al., 2008; Neiva et al., 2018; Bringloe et al., 2021), implying a predominant track of Pacific to Atlantic invasion. The biogeography of Fucus, however, does not suggest a primary Pacific origin. The sister genus Pelvetiopsis is composed by 4 NE Pacific-endemic species (Neiva et al., 2017), whereas all 7+ extant congeners of F. distichus are Atlantic endemics, including closely related F. serratus (F. spiralis is regarded as a recent introduction to the NE pacific (Coyer et al., 2006). The split and ensuing vicariant diversification of Fucus and Pelvetiopsis has been hypothesized to result from a trans-Arctic crossing by a common ancestor following the first opening of the Bering Strait at the Miocene/Pliocene transition ca. 5.6 Ma ago (Cánovas et al., 2011; Sousa et al., 2019). The secondary colonization of the Pacific by a F. distichus ancestor represents the most parsimonious hypothesis (Lindstrom, 2001; Cánovas et al., 2011) but is not consensual (Coyer et al., 2006; Laughinghouse et al., 2015). Alternative scenarios include a trans-Arctic split between main Fucus clades I (F. distichus/F. serratus) and II (F. vesiculosus et al.), or an earlier Pacific colonization by a F. distichus/F. serratus ancestor, in both cases involving the Atlantic colonization by (and Pacific extinction of) the F. serratus ancestor. While available data cannot answer this question, it is clear that the inferred colonization of the Atlantic represents only the latest among other trans-Arctic migrations in the evolutionary history of Fucus.

Long-term Pacific vs Atlantic ranges and niche conservatism

Another factor that likely contributed to the greater genetic diversity and differentiation of F. distichus in the North Pacific is the much wider ice-free glacial ranges, and the much greater glacial/interglacial range overlap, when compared with the Arctic/Atlantic. ENMs revealed extensive ice-free habitats suitable for F. distichus in the N Pacific during past major glaciations, from the Korean and Baja California Peninsulas to as far north as the palaeo-shorelines of Beringia. Even in the NE Pacific, where the Cordilleran ice-sheet spread along the coast and potentially created a distributional gap from British Columbia to the eastern Aleutians, suitable conditions existed from Vancouver to Point Conception (and beyond), where the endemic (and seemingly reproductively isolated) lineages R and LG were identified. Globally, and notwithstanding major climatic shifts across the last glacial/interglacial transition, available evidence supports long-term persistence of F. distichus throughout much of the core (and very wide) contemporary Pacific range. By comparison, the core distribution of F. distichus in the Atlantic overlaps extensively with regions that were under the direct influence of the LGM ice-sheets and available for colonization only at the onset of the deglaciation. The existence of a lineage unique to the NW Atlantic is in line with ENMs supporting suitable conditions just south of the Laurentide ice-sheet and along the ice-free palaeo-shorelines of the Grand Banks of Newfoundland (see also Assis et al., 2018a; Bringloe and Saunders, 2018). In the NE Atlantic, ENMs also revealed suitable habitat beyond ice-sheets from north-western British Isles to at least SW Iberia, but this putative southern glacial range is only marginally occupied today. The proximity of the glacial leading edge (just south of ice-sheet limits) and the modern trailing edge (western Ireland) necessarily translates into minimal areas of overlap (i.e. of long-term persistence), where genetic diversity could have accumulated and persisted until today.

Oddly, ENMs reveal suitable present-day climatic conditions throughout much of the British Isles and even Brittany, with Pacific models extending potential distributions as far south as northwest Iberia, i.e. much beyond its core distribution in Norway and its absolute modern range limit in Ireland. Even throughout its southern range in the British Isles, F. distichus populations are exceptionally rare and localized, and exhibit much less morphological and ecological plasticity when compared to other N Atlantic regions (e.g. Powell, 1963). Indeed, the single ecotype present – dwarf F. distichus var. anceps – is restricted to specific and disconnected sites along the western-facing coasts of Scotland and Ireland, subjected to extreme and persistent wave action, where other intertidal fucoids are scarce or absent. The species is ubiquitous throughout equivalent cold-temperate latitudes in the NE Pacific, where it nearly reaches the Southern California Bight warm-temperate eco-region (sensu Spalding et al., 2007). Most Arctic-cold temperate seaweeds with amphiboreal or amphi-Atlantic distributions (including foundational kelps such as Laminaria digitata, Alaria esculenta and Saccharina latissima, intertidal fucoids such as Ascophyllum nodosum and Fucus serratus, and red seaweeds such as Phycodrys rubens and Delesseria sanguinea) are still common in the Brittany-Ireland transition zone, and some even have relic populations in upwelling areas of NW Iberia. Most display unique gene-pools and diversity hotspots in these areas, indicating stable presence across climatic cycles (Wares and Cunningham, 2001; Maggs et al., 2008; Waltari and Hickerson, 2013).

The fragmented distribution of F. distichus in Europe could be explained by ongoing expansion from the NW Atlantic. This hypothesis is unconventional since most post-glacial trans-Atlantic migration has occurred, because of more benign glacial conditions, from East to West, but it explains some of the observations above. The steady expansion of a sheltered form of F. distichus (F. distichus var. evanescens) along the Skagerrak and Kattegat to at least Rostock in the western Baltic Sea (ca. 54° Lat.) since its mid-1890 introduction to Oslofjord (eastern Norway, presumably from northern Europe) is also compatible with this hypothesis, as it confirms suitable thermal regimes and potential for southern range expansions (Coyer et al., 2002; Dietrich and Hendrik, 2017). Atypically high rates of hybridization with F. serratus further suggest that contact between the species in this area is evolutionarily recent (Hoarau et al., 2015). Disentangling among these hypotheses – relic populations left after post-glacial contraction vs ongoing trans-Atlantic (W → E) post-glacial expansion – requires assessing whether genetic variation in the NE Atlantic is to some extent unique or chiefly a subset embedded within NW Atlantic diversity. The NE Atlantic populations, including Irish and Scottish ones, are not particularly remarkable for the studied markers, but more variable markers may be necessary to conclusively answer this question.

Perhaps other factors are playing a role. Whether pre-glacial or more recently, F. distichus may have suffered an important functional bottleneck during the NE Atlantic colonization. Migration corridors along the Canadian Arctic and Greenland require tolerance to extreme Arctic conditions and thus selective filters and gene-surfing along obligatory trans-Arctic pathways could have selected for specialized high-latitude form(s) pre-adapted to colder environments. If so, the regional NE Atlantic “niche” could be a fraction of the global species niche and its “relic” distribution along warmer British Isles (and absence further south) a reflection of it. There, its dependency for extremely exposed habitats could be explained by lower heat/desiccation stress attained in low-shore turbulent environments. The greater thermal tolerance of F. distichus in Pacific vs Atlantic ranges inferred from niche modelling (~1.8°C) is certainly compatible with this hypothesis, but empirical evidence for less thermal tolerance in the NE Atlantic range (or larger thermal tolerance in southern populations of the E Pacific) is missing entirely. Another potential factor to consider is inter-specific competition. F. distichus may be outcompeted by other intertidal fucoids toward temperate latitudes, becoming restricted to extremely exposed habitats (British Isles) or nutrient-enriched harbors (Kattegat and Western Baltic Sea) where other intertidal fucoids (and particularly low-intertidal Fucus serratus, see Jueterbock et al., 2016) are scarce or cannot establish. Intertidal fucoid communities are also represented in the southernmost areas of the NW Atlantic and NE Pacific (Silvetia, Pelvetiopsis), but are comparatively less diverse and saturated than in the temperate NE Atlantic. These contrasting hypotheses differ conceptually – they either assume fundamental or realized niche differences (Figure 6) – but both are rarely considered in marine phylogeographic literature. Comparative physiological/transcriptomic responses and reciprocal transplants, beyond the scope of this study, are likely to bring new insights, and help explain the nature of these apparent range/niche differences between the NE Atlantic and the remaining range.

Figure 6

Pacific refugia and melting pots

In the Pacific, and to a lesser extent in the Atlantic (e.g. Canadian Maritimes), the occurrence of endemic haplotypic radiations provide good indication of glacial range dissections and confirm that some lineages never expanded post-glacially much beyond past refugia. As predicted, the Montereyan region in the southernmost NE Pacific was home to two unique haplogroups (LG and R), the latter the most differentiated within the complex. At higher latitudes, but still beyond the reach of glacial ice, populations tended to be dominated by single haplogroups. For instance, between Vancouver Island and Cape Mendocino, distinct NE Pacific-endemic haplogroups dominated relatively short stretches of shoreline. Populations from Japan and Russia were also dominated by single haplogroups, some endemic to (or nearly so) the NW Pacific (P, DY). Contrastingly, in southern Alaska, NE Pacific haplogroups were locally mixed, also with the trans-Pacific Y haplogroup, resulting in the most haplotypic (and genotypic) diverse populations within F. distichus range. Alaskan populations, even when separately analyzed as high- and low- intertidal bands, also exhibited the largest observed heterozygosity. Ice-sheet reconstructions suggest that the Cordilleran ice-sheet reached the southern Alaskan coastline, presumably creating a distributional gap spanning from the eastern Aleutians to Vancouver Island. This region apparently lacks endemic haplogroups, as the ones co-occurring in the area were also found farther south or farther west. Regional diversity was also not atypically high, at least when comparing to the Montereyan/Mendocinian region, and pair-wise differentiation between populations was quite depressed, indicating weaker genetic structure and lineage sorting when compared to other regions. Finally, local populations were often composed of parapatric bands reminiscent of distinct ecotypes (even if not consistently recovered as separate entities, see next section).

Taken together, these data suggest that the Alaskan Gulf represents a melting pot where unique and more structured southern East and West Pacific sources met. Most Y and DG haplotypes samples were, however, private to this area. It remains a possibility that they are shared with contiguous, unsampled areas beyond the direct influence of the coastal ice-sheet, but if not, they could signal glacial persistence at these higher latitudes as in the kelps Alaria, Hedophyllum nigripes and Saccharina latissima (Grant and Bringloe, 2020; Grant et al., 2020; Grant and Chenoweth, 2021). The mixed evidence is compatible with a more complex scenario were F. distichus survived glacial advances in SE Alaska in a few small pockets of unglaciated shorelines, before invasion from southern and western ranges following ice-retreat – i.e. the area representing both a relic northern refugia and a post-glacial melting pot.

The Aleutian archipelago was poorly sampled but was also expected to bear signatures of recent contact between NW and NE Pacific regions. Central Aleutian Adak showed a more westerner character. This possibly reflects more benign LGM conditions, which presumably allowed survival (or at least earlier colonization) of western vs eastern Aleutians, as the latter was under the influence of the Cordilleran ice-sheet. Additional collections (analysed at a later stage and not included in this study) reveal that eastern Aleutian Islands indeed are genetically closer to the contiguous Alaskan region.

Species cohesion at local vs larger scale

In line with previous studies, molecular data failed to support traditionally recognized morpho-taxa among Atlantic F. distichus, nor, for the most part, subdividing Pacific “Fucus gardneri”. Much of the morphological variation between ecotypes has been shown to be continuous and explained by environmental factors (Powell, 1957; Rice et al., 1985; Sideman and Mathieson, 1985; Ruuskanen and Nappu, 2005), although a few common garden experiments suggest that at least some variation is heritable and not environmentally-induced (Sideman and Mathieson, 1985). In any case, eco-morphological or genetic entities were largely incongruent – some phylogroups included multiple ecotypes, and vice-versa, similar ecotypes belonged to distinct genetic lineages. In other words, at large spatial scales, F. distichus is best described as a single ecologically plastic, morphologically polymorphic and genetically structured species whose variation (genetic, morphologic, etc.) is not sufficiently discontinuous or consistent as to justify recognizing separate genetic or taxonomic entities.

Still, it remains challenging to fully capture and define the biodiversity contained within F. distichus. In some regions clearly distinct ecotypes matching distinct gene-pools co-existed at very local scales as stable sympatric/parapatric entities. This pattern seems to imply (local) isolation mechanisms and is reminiscent of incipient species (Billard et al., 2005; Neiva et al., 2018). Fucus are externally fertilizing broadcast spawners and show broad gametic compatibility, as revealed by natural hybridization/introgression between many species pairs (Coyer et al., 2002; Neiva et al., 2010). Remarkably, hybridization rates between F. distichus and F. serratus are reduced in naturally sympatric ranges when compared to recently established (19^th century) contact zones, implying reinforcement, i.e. the strengthening of pre-zygotic isolating barriers in response to selection against less-fit hybrids (Hoarau et al., 2015). Local build-up of isolation among younger and less differentiated conspecific populations has also been documented in the context of range-shifts, ecological speciation and ring-species (Hopkins, 2013). Many F. distichus haplogroups evolved in allopatry and in some cases may be partially pre-adapted, as suggested in parapatric contexts, to different environmental conditions. Upon secondary contact, their local and regional integrity may be explained by limited gene-flow resulting from high selfing rates, spatial segregation and/or temporal reproductive asynchrony, potentially reinforced by reduced hybrid fitness and/or strong selection for contrasting environmental conditions (Sideman and Mathieson, 1983; Coyer et al., 2002; Billard et al., 2010; Zardi et al., 2011; Hoarau et al., 2015; Monteiro et al., 2016). In Maine, for instance, a high-intertidal tide-pool ecotype (LO haplogroup) reproduces during low-tide when it is isolated from larger, lower-shore (DO haplotype) conspecifics (Pearson and Brawley, 1996), but low microsatellite resolution prevents at present confirming their degree of reproductive isolation. Dwarf and normal-sized intertidal ecotypes from Maine have also been shown to exhibit reproductive asynchrony (Sideman and Mathieson, 1983). In Alaska, high- and low-intertidal bands often showed different haplotypic compositions, but the pattern tended to be site-specific and only weakly correlated with respective nuclear backgrounds, requiring further examination. This area was also the one where populations exhibited the largest observed heterozygosity, implying an important incidence of outcrossing. The most convincing example is perhaps that of Pigeon Point, where haplogroups R (wide-frond type) and LG (thin-frond type) co-occurred sympatrically. Despite overlapping vertical ranges and reproductive period, only a single genetic F1 hybrid was detected (assuming it is not a laboratory contamination during DNA extraction), suggesting nearly complete reproductive isolation. Ongoing analyses of nearby populations support this pattern at a regional scale. This R haplogroup is also the most divergent within the complex, so there is convincing additional ground for suspecting a separate, incipient species.

From a taxonomic perspective, this “paradox” – the existence of putatively isolated genetic entities that are relatively easy to recognize and discriminate at local scales but much harder to define (in terms of distinctive and constant morphological, ecological or genetic characters) at larger spatial scale - is problematic. Dealing with these entities may be even more complicated if reproductive isolation, as seems to be the case in some regions, is contingent and varies geographically depending on the specific genetic entities involved and/or the evolutionary and environmental contexts underlying their regional assembly. Similar studies may show that similar patterns apply to other highly-structured plastic species spanning large geographical areas regularly heavily impacted by glaciations.

Final considerations and future directions

Widespread amphiboreal marine species such as F. distichus are ideal models to investigate how climatic shifts, dispersal barriers, environmental filters and distinct biotic assemblages influence intra-specific genetic structure and niche differentiation across species ranges. Despite the unprecedented effort, a range of open questions remain that could be further addressed in future studies. First, new markers and approaches should be employed to estimate the timing of trans-Arctic migrations and secondary colonization of the Atlantic, and particularly of the NE Atlantic rear-edge, that are difficult to establish with microsatellite and mitochondrial data alone. The same applies with respect to spatially co-existing ecotypes, for which genomic data can provide much more power to clarify the actual degree of divergence and isolation and best taxonomic status, including confirming the putative cryptic species in the eastern Pacific rear-edge. As mentioned above, experimental work also seems necessary to clarify the nature of the apparent and very peculiar niche unfilling in the NE Atlantic. Finally, it is unquestionable that ongoing climatic change will translate into wider areas and windows of ice-free conditions in the near future. Deeper expansions into the high-Arctic seem inevitable and might even benefit from increased shipping along Canadian and Russian Arctic shores (Jueterbock et al., 2016). Vicariant Pacific and Atlantic lineages may eventually meet, creating an excellent biogeographic opportunity to directly investigate the fate of divergent lineages upon secondary contact (e.g. Neiva et al., 2018). These opportunities of research will not only improve our understanding of the complex biogeographical history, trans-Arctic divergence, and species cohesion in Fucus distichus, but of amphiboreal communities in general.

References

• 1

  AdeyW. H.LindstromS. C.HommersandM. H.MullerK. M. (2008). The biogeographic origin of arctic endemic seaweeds: A thermogeographic view. J. Phycology44, 1384–1394. doi: 10.1111/j.1529-8817.2008.00605.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AlmeidaS. C.NeivaJ.SousaF.MartinN.CoxC. J.Melo-FerreiraJ.et al. (2022). A low-latitude species pump: Peripheral isolation, parapatric speciation and mating-system evolution converge in a marine radiation. Mol. Ecol.31, 4797–4817. doi: 10.1111/mec.16623

  • CrossRef
  • Google Scholar

• 3

  AssisJ.AraújoM. B.SerrãoE. A. (2018a). Projected climate changes threaten ancient refugia of kelp forests in the North Atlantic. Global Change Biol.24, e55–e66. doi: 10.1111/gcb.13818

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  AssisJ.CoelhoN. C.LamyT.ValeroM.AlbertoF.SerrãoE. A. (2016). Deep reefs are climatic refugia for genetic diversity of marine forests. J. Biogeography43, 833–844. doi: 10.1111/jbi.12677

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  AssisJ.FragkopoulouE.FradeD.NeivaJ.OliveiraA.AbecasisD.et al. (2020). A fine-tuned global distribution dataset of marine forests. Sci. Data7, 119. doi: 10.1038/s41597-020-0459-x

  • CrossRef
  • Google Scholar

• 6

  AssisJ.SerrãoE. A.CoelhoN. C.TemperaF.ValeroM.AlbertoF. (2018b). Past climate changes and strong oceanographic barriers structured low - latitude genetic relics for the golden kelp Laminaria ochroleuca. J. Biogeography45, 2326–2336. doi: 10.1111/jbi.13425

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  AssisJ.TybergheinL.BoshS.VerbruggenH.SerrãoE. A.De ClerckO. (2018c). Bio-ORACLE v2.0: extending marine data layers for bioclimatic modelling. Global Ecol. Biogeograph. doi: 10.1111/geb.12693

  • CrossRef
  • Google Scholar

• 8

  Barbet-MassinM.JiguetF.AlbertC. H.ThuillerW. (2012). Selecting pseudo-absences for species distribution models: how, where and how many? Methods Ecol. Evol.3, 327–338. doi: 10.1111/j.2041-210X.2011.00172.x

  • CrossRef
  • Google Scholar

• 9

  BillardE.DaguinC.PearsonG. A.SerrãoE. A.EngelC.ValeroM. (2005). Genetic isolation between three closely related taxa: Fucus vesiculosus, F. spiralis and F. ceranoides (Phaophyceae). J. Phycology41, 900–905. doi: 10.1111/j.0022-3646.2005.04221.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BillardE.SerrãoE. A.PearsonG. A.DestombeC.ValeroM. (2010). Fucus vesiculosus and spiralis species complex: a nested model of local adaptation at the shore level. Mar. Ecol. Prog. Ser.405, 163–174. doi: 10.3354/meps08517

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  BoavidaJ.AssisJ.SilvaI.SerrãoE. A. (2016). Overlooked habitat of a vulnerable gorgonian revealed in the Mediterranean and Eastern Atlantic by ecological niche modelling. Sci. Rep.6, 36460. doi: 10.1038/srep36460

  • CrossRef
  • Google Scholar

• 12

  BringloeT. T.SaundersG. W. (2018). Mitochondrial DNA sequence data reveal the origins of postglacial marine macroalgal flora in the Northwest Atlantic. Mar. Ecol. Prog. Ser.589, 45–58. doi: 10.3354/meps12496

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  BringloeT. T.ZaparenkovD.StarkoS.GrantW. S.VieiraC.KawaiH.et al. (2021). Whole-genome sequencing reveals forgotten lineages and recurrent hybridizations within the kelp genus Alaria (Phaeophyceae). J. Phycology57, 1721–1738. doi: 10.1111/jpy.13212

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  CánovasF. G.MotaC. F.SerrãoE. A.PearsonG. A. (2011). Driving south: a multi-gene phylogeny of the brown algal family Fucaceae reveals relationships and recent drivers of a marine radiation. BMC Evol. Biol.11, 371. doi: 10.1186/1471-2148-11-371

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  CarrC. M.HardyS. M.BrownT. M.MacdonaldT. A.HebertP. D. N. (2011). A tri-oceanic perspective: DNA barcoding reveals geographic structure and cryptic diversity in Canadian polychaetes. PloS One6. doi: 10.1371/journal.pone.0022232

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  ChernovaE. N. (2016). The biogeochemical background and trace metal accumulation by brown algae of the genus Fucus in coastal waters of the Sea of Japan, the Sea of Okhotsk, and the White Sea. Russian J. Mar. Biol.42, 87–96. doi: 10.1134/S1063074016010053

  • CrossRef
  • Google Scholar

• 17

  Conservation of Arctic Flora and Fauna. (2010). Arctic biodiversity trends 2010: selected indicators of change. Available online at: https://wedocs.unep.org/20.500.11822/8670.

  • Google Scholar
  • View reference in article

• 18

  CoyerJ. A.HoarauG.BeszteriB.PearsonG.OlsenJ. L. (2009). Expressed sequence tag-derived polymorphic SSR markers for Fucus serratus and amplification in other species of Fucus. Mol. Ecol. Resour.9, 168–170. doi: 10.1111/j.1755-0998.2008.02406.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  CoyerJ. A.HoarauG.CostaJ. F.HogerdijkB.SerrãoE.BillardE.et al. (2011a). Evolution and diversification within the intertidal brown macroalgae Fucus spiralis/F. vesiculosus species complex in the North Atlantic. Mol. Phylogenet. Evol.58, 283–296. doi: 10.1016/j.ympev.2010.11.015

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  CoyerJ. A.HoarauG.Oudot-Le SecqM.-P.StamW. T.OlsenJ. L. (2006). A mtDNA-based phylogeny of the brown algal genus Fucus (Heterokontophyta; Phaeophyta). Mol. Phylogenet. Evol.39, 209–222. doi: 10.1016/j.ympev.2006.01.019

  • CrossRef
  • Google Scholar

• 21

  CoyerJ. A.HoarauG.Van SchaikJ.LuijckxP.OlsenJ. L. (2011b). Trans-Pacific and trans-Arctic pathways of the intertidal macroalga Fucus distichus L. reveal multiple glacial refugia and colonizations from the North Pacific to the North Atlantic. J. Biogeography38, 756–771. doi: 10.1111/jbi.2011.38.issue-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  CoyerJ. A.PetersA. F.HoarauG.StamW. T.OlsenJ. L. (2002). Hybridization of the marine seaweeds, Fucus serratus and Fucus evanescens (Heterokontophyta: Phaeophyceae) in a 100-year-old zone of secondary contact. Proc. R. Soc. B / Biol. Sci.269, 1829–1834. doi: 10.1098/rspb.2002.2093

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  DietrichA.HendrikS. (2017). Fucus evanescens or Fucus edentatus? Taxonomic problems with a non-indigenous species, which arrived at the coast of Mecklenburg-Vorpommern. Rostock. Meeresbiolog. Beitr.27, 117–125.

  • Google Scholar
  • View reference in article

• 24

  DodsonJ. J.TremblayS.ColombaniF.CarscaddenJ. E.LecomteF. (2007). Trans-Arctic dispersals and the evolution of a circumpolar marine fish species complex, the capelin (Mallotus villosus). Mol. Ecol.16, 5030–5043. doi: 10.1111/j.1365-294X.2007.03559.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  DuntonK. H. (1992). Arctic biogeography: the paradox of the marine benthic fauna and flora. Trends Ecol. Evol.7, 183–189. doi: 10.1016/0169-5347(92)90070-R

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  ElithJ.LeathwickJ. R.HastieT. (2008). A working guide to boosted regression trees. J. Anim. Ecol.77, 802–813. doi: 10.1111/j.1365-2656.2008.01390.x

  • CrossRef
  • Google Scholar

• 27

  EngelC. R.BrawleyS. H.EdwardsK. J.SerraoE. (2003). Isolation and cross-species amplification of microsatellite loci from the fucoid seaweeds Fucus vesiculosus, F. serratus and Ascophyllum nodosum (Heterokontophyta, Fucaceae). Mol. Ecol. Notes3, 180–182. doi: 10.1046/j.1471-8286.2003.00390.x

  • CrossRef
  • Google Scholar

• 28

  EvannoG.RegnautS.GoudetJ. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol.14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  FragkopoulouE.SerrãoE. A.HortaP. A.KoerichG.AssisJ. (2021). Bottom trawling threatens future climate refugia of rhodoliths globally. Front. Mar. Sci.7, 1–11. doi: 10.3389/fmars.2020.594537

  • CrossRef
  • Google Scholar

• 30

  FragkopoulouE.SerrãoE. A.De ClerckO.CostelloM. J.AraújoM. B.DuarteC. M.et al. (2022). Global biodiversity patterns of marine forests of brown macroalgae. Global Ecol. Biogeography31, 636–648. doi: 10.1111/geb.13450

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  FraserC. I.NikulaR.SpencerH. G.WatersJ. M. (2009). Kelp genes reveal effects of subantarctic sea ice during the Last Glacial Maximum. Proc. Natl. Acad. Sci. USA106, 3249–3253. doi: 10.1073/pnas.0810635106

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  GrantW. S.BringloeT. T. (2020). Pleistocene ice ages created new evolutionary lineages, but limited speciation in Northeast Pacific Winged Kelp. J. Heredity111, 593–605. doi: 10.1093/jhered/esaa053

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  GrantW. S.ChenowethE. (2021). Phylogeography of sugar kelp: Northern ice-age refugia in the Gulf of Alaska. Ecol. Evol.11, 4670–4687. doi: 10.1002/ece3.7368

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  GrantW. S.LydonA.BringloeT. T. (2020). Phylogeography of split kelp Hedophyllum nigripes: northern ice-age refugia and trans-Arctic dispersal. Polar Biol.43, 1829–1841. doi: 10.1007/s00300-020-02748-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  HalanychK. M.VodotiE. T.SundbergP.DahlgrenT. G. (2013). Phylogeography of the horse mussel Modiolus modiolus. J. Mar. Biol. Assoc. United Kingdom93, 1857–1869. doi: 10.1017/S0025315413000404

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  HardyS. M.CarrC. M.HardmanM.SteinkeD.CorstorphineE.MahC. (2011). Biodiversity and phylogeography of Arctic marine fauna: Insights from molecular tools. Mar. Biodiversity41, 195–210. doi: 10.1007/s12526-010-0056-x

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  HoarauG.CoyerJ. A.GiesbersM. C. W. G.JueterbockA.OlsenJ. L. (2015). Pre-zygotic isolation in the macroalgal genus Fucus from four contact zones a tale of reinforcement? R. Soc. Open Sci.2, 140538. doi: 10.1098/rsos.140538

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  HofnerB.MüllerJ.HothornT. (2011). Monotonicity-constrained species distribution models. Ecology92 (10), 1895–1901. doi: 10.1890/10-2276.1

  • CrossRef
  • Google Scholar

• 39

  HoltR. D. (2009). Bringing the Hutchinsonian niche into the 21st century: Ecological and evolutionary perspectives. Proc. Natl. Acad. Sci.106, 19659–19665. doi: 10.1073/pnas.0905137106

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  HopkinsR. (2013). Minireview reinforcement in plants. New Phytol.197, 1095–1103. doi: 10.1111/nph.12119

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  Jiménez-ValverdeA.LoboJ. M. (2007). Threshold criteria for conversion of probability of species presence to either–or presence–absence. Acta Oecolog.31, 361–369. doi: 10.1016/j.actao.2007.02.001

  • CrossRef
  • Google Scholar

• 42

  JueterbockA.SmolinaI.CoyerJ. A.HoarauG. (2016). The fate of the Arctic seaweed Fucus distichus under climate change: an ecological niche modeling approach. Ecol. Evol.6, 1712–1724. doi: 10.1002/ece3.2001

  • CrossRef
  • Google Scholar

• 43

  KaufmanD. S.AgerT. A.AndersonN. J.AndersonP. M.AndrewsJ. T.BartleinP. J.et al. (2004). Holocene thermal maximum in the western Arctic (0-180°W). Quaternary Sci. Rev.23, 529–560. doi: 10.1016/j.quascirev.2003.09.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  KeenanK.McginnityP.CrossT. F.CrozierW. W.ProdöhlP. A. (2013). DiveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors. Methods Ecol. Evol.4, 782–788. doi: 10.1111/2041-210X.12067

  • CrossRef
  • Google Scholar

• 45

  KingN. G.McKeownN. J.SmaleD. A.MooreP. J. (2018). The importance of phenotypic plasticity and local adaptation in driving intraspecific variability in thermal niches of marine macrophytes. Ecography41, 1469–1484. doi: 10.1111/ecog.03186

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  Krause-JensenD.ArchambaultP.AssisJ.BartschI.BischofK.Filbee-DexterK.et al. (2020). Imprint of climate change on Pan-Arctic marine vegetation. Front. Mar. Sci.7, 1–27. doi: 10.3389/fmars.2020.617324

  • CrossRef
  • Google Scholar

• 47

  KuceraH.SaundersG. W. (2008). Assigning morphological variants of Fucus (Fucales, Phaeophyceae) in Canadian waters to recognized species using DNA barcoding. Botany86, 1065–1079. doi: 10.1139/B08-056

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  LaakkonenH. M.LajusD. L.StrelkovP.VäinöläR. (2013). Phylogeography of amphi-boreal fish: tracing the history of the Pacific herring Clupea pallasii in North-East European seas. BMC Evolutionary Biol.13, 67. doi: 10.1186/1471-2148-13-67

  • CrossRef
  • Google Scholar
  • View reference in article

• 49

  LaneC. E.LindstromS. C.SaundersG. W. (2007). A molecular assessment of northeast Pacific Alaria species (Laminariales, Phaeophyceae) with reference to the utility of DNA barcoding. Mol. Phylogenet. Evol.44, 634–648. doi: 10.1016/j.ympev.2007.03.016

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  LaughinghouseH. D.MüllerK. M.AdeyW. H.LaraY.YoungR.JohnsonG. (2015). Evolution of the northern rockweed, Fucus distichus, in a regime of glacial cycling: implications for benthic algal phylogenetics. PlosOne10, e0143795. doi: 10.1371/journal.pone.0143795

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  LiesnerD.FouqueauL.ValeroM.RoledaM. Y.PearsonG. A.BischofK.et al. (2020). Heat stress responses and population genetics of the kelp Laminaria digitata (Phaeophyceae) across latitudes reveal differentiation among North Atlantic populations. Ecol. Evol.10, 9144–9177. doi: 10.1002/ece3.6569

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  LindstromS. C. (2001). The Bering Strait connection: dispersal and speciation in boreal macroalgae. J. Biogeography28, 243–251. doi: 10.1046/j.1365-2699.2001.00529.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  LobusN. V.UdalovA. A. (2021). The chemical composition of brown algae laminaria digitata (Hudson) J.V. Lamouroux 1893 and Fucus distichus (Linnaeus 1767) from the bays of the novaya zemlya archipelago (the kara sea). Russian J. Mar. Biol.47, 407–412. doi: 10.1134/S1063074021050096

  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  LüningK. (1990). Seaweeds: their environment, biogeography, and ecophysiology (Revised ed) (New York: Wiley).

  • Google Scholar
  • View reference in article

• 55

  MaggsC. A.CastilhoR.FoltzD.HenzlerC.JollyM. T.KellyJ.et al. (2008). Evaluating signatures of glacial refugia for north atlantic benthic marine taxa. Ecology89, 108–122. doi: 10.1890/08-0257.1

  • CrossRef
  • Google Scholar

• 56

  MarkoP. B.HoffmanJ. M.EmmeS. A.McgovernT. M.KeeverC. C.CoxL. N. (2010). The ‘ Expansion – Contraction ‘ model of Pleistocene biogeography : rocky shores suffer a sea change? Mol. Ecol.19 (1), 146–169. doi: 10.1111/mec.2009.19.issue-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  MartinsM. R.AssisJ.AbecasisD. (2021). Biologically meaningful distribution models highlight the benefits of the Paris agreement for demersal fishing targets in the North Atlantic Ocean. Glob. Ecol. Biogeogr.30 (8), 1643–1656. doi: 10.1111/geb.13327

  • CrossRef
  • Google Scholar

• 58

  MitchellC. E.AgawalA. A.BeverJ. D.GilbertG. S.RuthA.KlironomosJ. N.et al. (2006). Biotic interactions and plant invasions. Ecol. Lett.9, 726–740. doi: 10.1111/j.1461-0248.2006.00908.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  MohringM. B.WernbergT.WrightJ. T.ConnellS. D.RussellB. D. (2014). Biogeographic variation in temperature drives performance of kelp gametophytes during warming. Mar. Ecol. Prog. Ser.513, 85–96. doi: 10.3354/meps10916

  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  MontecinosA.BroitmanB. R.FaugeronS.HayeP.TellierF.GuilleminM.-L. (2012). Species replacement along a linear coastal habitat: phylogeography and speciation in the red alga Mazzaella laminarioides along the south east Pacific. BMC Evolutionary Biol.12, 97. doi: 10.1186/1471-2148-12-97

  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  MonteiroC. A.PaulinoC.JacintoR.SerrãoE. A.PearsonG. A. (2016). Temporal windows of reproductive opportunity reinforce species barriers in a marine broadcast spawning assemblage. Sci. Rep.6, 29198. doi: 10.1038/srep29198

  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  NeivaJ.AssisJ.FernandesF.PearsonG. A.SerrãoE. A. (2014). Species distribution models and mitochondrial DNA phylogeography suggest an extensive biogeographical shift in the high-intertidal seaweed Pelvetia canaliculata. J. Biogeography41, 1137–1148. doi: 10.1111/jbi.12278

  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  NeivaJ.HansenG. I.PearsonG. A.VlietM. S.MaggsC. A.SerrãoE. A. (2012a). Fucus cottonii (Fucales, Phaeophyceae) is not a single genetic entity but a convergent salt-marsh morphotype with multiple independent origins. Eur. J. Phycology47, 461–468. doi: 10.1080/09670262.2012.736536

  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  NeivaJ.PaulinoC.NielsenM. M.Krause-JensenD.SaundersG. W.AssisJ.et al. (2018). Glacial vicariance drives phylogeographic diversification in the amphi-boreal kelp Saccharina latissima. Sci. Rep.8, 1112. doi: 10.1038/s41598-018-19620-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  NeivaJ.PearsonG. A.ValeroM.SerrãoE. A. (2010). Surfing the wave on a borrowed board: range expansion and spread of introgressed organellar genomes in the seaweed Fucus ceranoides L. Mol. Ecol.19, 4812–4822. doi: 10.1111/j.1365-294X.2010.04853.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 66

  NeivaJ.PearsonG. A.ValeroM.SerrãoE. A. (2012b). Fine-scale genetic breaks driven by historical range dynamics and ongoing density-barrier effects in the estuarine seaweed Fucus ceranoides L. In BMC Evolutionary Biol.12, 78. doi: 10.1186/1471-2148-12-78

  • CrossRef
  • Google Scholar
  • View reference in article

• 67

  NeivaJ.SerrãoE. A.AndersonL.RaimondiP. T.MartinsN.GouveiaL.et al. (2017). Cryptic diversity, geographical endemism and allopolyploidy in NE Pacific seaweeds. BMC Evolutionary Biol.17, 30. doi: 10.1186/s12862-017-0878-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 68

  NeivaJ.SerrãoE. A.AssisJ.PearsonG. A.CoyerJ. A.OlsenJ. L.et al. (2016). “Climate oscillations, range shifts and phylogeographic patterns of north atlantic fucaceae,” in Seaweed phylogeography - adaptation and evolution of seaweeds under environmental change. Eds. HuZ.-M.FraserC.ed.1. (Springer Netherlands), 279–308.

  • Google Scholar
  • View reference in article

• 69

  NikulaR.StrelkovP.VäinöläR. (2007). Diversity and trans-arctic invasion history of mitochondrial lineages in the North Atlantic Macoma balthica complex (Bivalvia: Tellinidae). Evolution61, 928–941. doi: 10.1111/j.1558-5646.2007.00066.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 70

  PearsonG. A.BrawleyS. H. (1996). Reproductive ecology of Fucus distichus (Phaeophyceae): An intertidal alga with successful external fertilization. Mar. Ecol. Prog. Ser.143, 211–223. doi: 10.3354/meps143211

  • CrossRef
  • Google Scholar
  • View reference in article

• 71

  PerrinC.DaguinC.VlietM.Van DeEngelC. R.PearsonG. A.et al. (2007). Implications of mating system for genetic diversity of sister algal species: Fucus spiralis and Fucus vesiculosus (Heterokontophyta, Phaeophyceae). Eur. J. Phycology42, 219–230. doi: 10.1080/09670260701336554

  • CrossRef
  • Google Scholar
  • View reference in article

• 72

  PicoT.MitrovicaJ. X.MixA. C. (2020). Sea level fingerprinting of the Bering Strait flooding history detects the source of the Younger Dryas climate event. Sci. Adv.6, eaay2935. doi: 10.1126/sciadv.aay2935

  • CrossRef
  • Google Scholar

• 73

  PowellH. T. (1957). Studies in the genus Fucus L. II. Distribution and ecology of forms of Fucus distichus L. Emend. Powell in Britain and Ireland. J. Mar. Biol. Assoc. UK36, 663–693. doi: 10.1017/S0025315400025923

  • CrossRef
  • Google Scholar

• 74

  PowellH. T. (1963). Speciation in the genus Fucus L., and related genera. Systematics Assoc. Publ.5, 63–67.

  • Google Scholar

• 75

  RiceE. L.KenchingtonT. J.ChapmanA. R. O. (1985). Intraspecific geographic-morphological variation patterns in Fucus distichus and F. evanescens. Mar. Biol.88, 207–215. doi: 10.1007/BF00397168

  • CrossRef
  • Google Scholar
  • View reference in article

• 76

  RuuskanenA. T.NappuN. P. (2005). Morphological differences in Fucus gardneri between two shores with equal cartographic exposure values but different levels of wave action. Annales Botanici Fennici42, 27–33.

  • Google Scholar
  • View reference in article

• 77

  SaadaG.NicastroK. R.JacintoR.McQuaidC. D.SerrãoE. A.PearsonG. A.et al. (2016). Taking the heat: distinct vulnerability to thermal stress of central and threatened peripheral lineages of a marine macroalga. Diversity Distributions22, 1060–1068. doi: 10.1111/ddi.12474

  • CrossRef
  • Google Scholar
  • View reference in article

• 78

  SasakiH.KawaiH. (2007). Taxonomic revision of the genus Chorda (Chordaceae, Laminariales) on the basis of sporophyte anatomy and molecular phylogeny. Phycologia46, 10–21. doi: 10.2216/06-06.1

  • CrossRef
  • Google Scholar
  • View reference in article

• 79

  SaundersG. W.McDevitD. C. (2013). DNA barcoding unmasks overlooked diversity improving knowledge on the composition and origins of the Churchill algal flora. BMC Ecol.13, 9. doi: 10.1186/1472-6785-13-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 80

  SerrãoE. A.AliceL. A.BrawleyS. H. (1999). Evolution of the Fucaceae (Phaeophyceae) inferred from nrDNA-ITS. J. Phycology35, 382–394. doi: 10.1046/j.1529-8817.1999.3520382.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 81

  SidemanE. J.MathiesonA. C. (1983). Ecological and genecological distinctions of a high intertidal dwarf form of Fucus distichus (L.) Powell in New England. J. Exp. Mar. Biol. Ecol.72, 171–188. doi: 10.1016/0022-0981(83)90142-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 82

  SidemanE. J.MathiesonA. C. (1985). Morphological variation within and between natural populations of non-tide pool Fucus distichus (Phaeophyta) in New England. J. Phycology21, 250–257. doi: 10.1111/j.0022-3646.1985.00250.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 83

  SmolinaI.KolliasS.JueterbockA.CoyerJ. A.HoarauG. (2016). Variation in thermal stress response in two populations of the brown seaweed, Fucus distichus, from the Arctic and subarctic intertidal. R. Soc. Open Sci.3, 150429. doi: 10.1098/rsos.150429

  • CrossRef
  • Google Scholar
  • View reference in article

• 84

  SousaF.NeivaJ.MartinsN.JacintoR.AndersonL.RaimondiP. T.et al. (2019). Increased evolutionary rates and conserved transcriptional response following allopolyploidisation in brown algae. Evolution73, 59–72. doi: 10.1111/evo.2019.73.issue-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 85

  SpaldingM. D.FoxH. E.AllenG. R.DavidsonN.FerdañaZ. A.FinlaysonM. A. X.et al. (2007). Marine ecoregions of the world: A bioregionalization of coastal and shelf areas. BioScience57, 573–583. doi: 10.1641/B570707

  • CrossRef
  • Google Scholar
  • View reference in article

• 86

  SwansonH. K.LysyM.PowerM.StaskoA. D.JohnsonJ. D.ReistJ. D. (2015). A new probabilistic method for quantifying n-dimensional ecological niches and niche overlap. Ecology96 (2), 318–324. doi: 10.1890/14-0235.1

  • CrossRef
  • Google Scholar

• 87

  TellierF.TapiaJ.FaugeronS.DestombeC.ValeroM. (2011). The Lessonia nigrescens species complex (Laminariales, Phaeophyceae) shows strict parapatry and complete reproductive isolation in a secondary contact zone. J. Phycology47, 894–903. doi: 10.1111/jpy.2011.47.issue-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 88

  TybergheinL.VerbruggenH.PaulyK.TroupinC.MineurF.De ClerckO. (2012). Bio-ORACLE: a global environmental dataset for marine species distribution modelling. Global Ecol. Biogeograph.21, 272–281. doi: 10.1111/j.1466-8238.2011.00656.xv

  • CrossRef
  • Google Scholar

• 89

  WallaceA. L.KleinA. S.MathiesonA. C. (2004). Determining the affinities of salt marsh fucoids using microsatellite markers: evidence of hybridization and introgression between two species of fucus (Phaeophyta) in a maine estuary. J. Phycology40, 1013–1027. doi: 10.1111/j.1529-8817.2004.04085.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 90

  WaltariE.HickersonM. J. (2013). Late Pleistocene species distribution modelling of North Atlantic intertidal invertebrates. J. Biogeography40, 249–260. doi: 10.1111/j.1365-2699.2012.02782.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 91

  WaresJ. P.CunninghamC. W. (2001). Phylogeography and historical ecology of the North Atlantic intertidal. Evolution55, 2455–2469. doi: 10.1111/j.0014-3820.2001.tb00760.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 92

  WiensJ. J.GrahamC. H. (2005). Niche conservatism: integrating evolution, ecology, and conservation biology. Annu. Rev. Ecology Evolution Systematics36, 519–539. doi: 10.1146/annurev.ecolsys.36.102803.095431

  • CrossRef
  • Google Scholar
  • View reference in article

• 93

  WilceR. T.DuntonK. H. (2014). The boulder Patch (North Alaska, Beaufort Sea) and its benthic algal flora. Arctic67, 43–56. doi: 10.14430/arctic4360

  • CrossRef
  • Google Scholar
  • View reference in article

• 94

  WulffA.IkenK.QuartinoM. L.Al-HandalA.WienckeC.ClaytonM. N. (2009). Biodiversity, biogeography and zonation of marine benthic micro- and macroalgae in the Arctic and Antarctic. Botanica Marina52, 491–507. doi: 10.1515/BOT.2009.072

  • CrossRef
  • Google Scholar

• 95

  YamadaI. (1980). Instructions for use Benthic marine algal vegetation along the coasts of Hokkaido , with special reference to the vertical distribution. Journ. Fac. Sci. Hokkaido Univ. Ser. V (Botany)12, 11–98.

  • Google Scholar
  • View reference in article

• 96

  ZardiG. I.NicastroK. R.CanovasF.CostaJ. F.SerrãoE. A.PearsonG. A. (2011). Adaptive traits are maintained on steep selective gradients despite gene flow and hybridization in the intertidal zone. PloS One6, e19402. doi: 10.1371/journal.pone.0019402

  • CrossRef
  • Google Scholar
  • View reference in article

• 97

  ZhangJ.YaoJ.-T.SunZ.-M.FuG.GalaninD.NagasatoC.et al. (2015). Phylogeographic data revealed shallow genetic structure in the kelp Saccharina japonica (Laminariales, Phaeophyta). BMC Evolutionary Biol.15, 237. doi: 10.1186/s12862-015-0517-8

  • CrossRef
  • Google Scholar
  • View reference in article