- Research article
Contemporary habitat discontinuity and historic glacial ice drive genetic divergence in Chilean kelp
BMC Evolutionary Biologyvolume 10, Article number: 203 (2010)
South America's western coastline, extending in a near-straight line across some 35 latitudinal degrees, presents an elegant setting for assessing both contemporary and historic influences on cladogenesis in the marine environment. Southern bull-kelp (Durvillaea antarctica) has a broad distribution along much of the Chilean coast. This species represents an ideal model taxon for studies of coastal marine connectivity and of palaeoclimatic effects, as it grows only on exposed rocky coasts and is absent from beaches and ice-affected shores. We expected that, along the central Chilean coast, D. antarctica would show considerable phylogeographic structure as a consequence of the isolating effects of distance and habitat discontinuities. In contrast, we hypothesised that further south - throughout the region affected by the Patagonian Ice Sheet at the Last Glacial Maximum (LGM) - D. antarctica would show relatively little genetic structure, reflecting postglacial recolonisation.
Mitochondrial (COI) and chloroplast (rbcL) DNA analyses of D. antarctica from 24 Chilean localities (164 individuals) revealed two deeply divergent (4.5 - 6.1% for COI, 1.4% for rbcL) clades from the centre and south of the country, with contrasting levels and patterns of genetic structure. Among populations from central Chile (32° - 44°S), substantial phylogeographic structure was evident across small spatial scales, and a significant isolation-by-distance effect was observed. Genetic disjunctions in this region appear to correspond to the presence of long beaches. In contrast to the genetic structure found among central Chilean populations, samples from the southern Chilean Patagonian region (49° - 56°S) were genetically homogeneous and identical to a haplotype recently found throughout the subantarctic region.
Southern (Patagonian) Chile has been recolonised by D. antarctica relatively recently, probably since the LGM. The inferred trans-oceanic ancestry of these Patagonian populations supports the notion that D. antarctica is capable of long-distance dispersal via rafting. In contrast, further north in central Chile, the correspondence of genetic disjunctions in D. antarctica with long beaches indicates that habitat discontinuity drives genetic isolation among established kelp populations. We conclude that rafting facilitates colonisation of unoccupied shores, but has limited potential to enhance gene-flow among established populations. Broadly, this study demonstrates that some taxa may be considered to have either high or low dispersal potential across different temporal and geographic scales.
Oceanic ecosystems have historically been considered relatively open, with high connectivity among populations . Recent studies have, however, demonstrated considerable population structure in a variety of coastal taxa (see [2–5]; and references therein), indicating that connectivity among marine populations is often lower than previously assumed . Levels of connectivity in coastal ecosystems may be affected by a broad range of factors, including oceanography, geological history, habitat continuity, ecology and life history . The relative importance of these factors in structuring marine populations, however, remains largely unknown .
Genetic patterns of isolation-by-distance are common in the marine realm, with neighbouring populations generally more closely related to each other than to those further away, as the chance of successful dispersal decreases with increasing distance [3, 6]. Gene flow among populations with particularly stringent ecological requirements may also be impeded by habitat discontinuities [2, 8]. Extensive stretches of sandy beach might, for example, represent important biogeographic hurdles for rocky-shore marine taxa, even those with a planktonic larval stage [7, 9–12]. The genetic effects of habitat discontinuities should be greatest among populations of taxa with relatively low dispersal potential. Indeed, in several species of marine macroalgae with predicted low dispersal potential, discontinuities in suitable substratum have been linked to reduced gene flow among populations [13–15].
Although most macroalgae have only short-lived gametes or dispersive zygotes, some - particularly highly buoyant species - are hypothesised to be capable of long-distance dispersal via rafting of fertile adults (e.g., Fucus vesiculosus: ; Hormosira banksii: ; Macrocystis pyrifera: ). A broad-scale phylogeographic study of one such species, southern bull-kelp (Durvillaea antarctica), found evidence of recent (postglacial) recolonisation of much of the subantarctic, a feat almost certainly achieved by trans-oceanic rafting of adult specimens . Recolonisation of the subantarctic was inferred from extremely low levels of genetic structure among bull-kelp populations throughout this region, versus high levels of genetic structure in glacial refugia (e.g., New Zealand). Low levels of genetic structure of populations in historically glaciated regions, versus relatively high structure in unglaciated regions, is a signature of postglacial recolonisation that has been demonstrated in a wide range of marine and terrestrial taxa (reviewed by [20–24]). But if bull-kelp is able to disperse over long distances by rafting, how is high genetic structure maintained in refugial regions? Fraser et al.  suggested that the genetic structure observed in refugial areas may be largely a consequence of the density-blocking effects of established populations [20, 25]. In dense kelp populations, any limited available substratum is more likely to be colonised by zygotes from residents rather than those from rare, rafted individuals. Phylogeographic structure within refugial areas may therefore be driven more through isolation-by-distance, or breaks in habitat continuity, than by rafting dispersal. Rafting may, in contrast, enable long-distance dispersal and colonisation of unoccupied habitats.
Chile presents an intriguing system for assessing the relative phylogeographic impacts of isolation-by-distance, habitat discontinuities, and historic glaciations. Central Chile has an almost straight coastline characterised by numerous extensive stretches of beach that isolate rocky-shore habitats ; this coast was largely unaffected by ice during the Pleistocene glaciations . Further south, Chilean Patagonian shores also form isolated habitats of islands and peninsulas, but at the Last Glacial Maximum (LGM), an extensive ice sheet covered most of this region [27, 28] (Fig. 1). This glacial ice probably extended westward to the edge of the continental shelf , and likely extirpated many populations of rocky shore taxa. Such contrasts in geological history between central and Patagonian Chile should be reflected in the phylogeographic structure of taxa occurring across both regions. Durvillaea antarctica currently dominates the rocky shores of Chile from Cape Horn , at around 56°S, up to approximately 32°S . Fraser et al.'s  recent broadscale study of D. antarctica included few samples from central Chile (only seven localities), and none from Chilean Patagonia, and was therefore unable to assess phylogeographic contrasts along the Chilean coast. For the current study, we analysed samples from a total of 24 localities, encompassing the full range of D. antarctica in Chile, to test the predictions that (i) D. antarctica populations in central Chile would show strong phylogeographic structure, and that (ii) D. antarctica populations in southern Chile would show little genetic structure, reflecting postglacial recolonisation.
DNA sequencing of 164 D. antarctica samples from 24 localities along the coast of Chile yielded 14 distinct haplotypes for mitochondrial COI, with 34 variable sites detected across the 629 bp fragment. For a chloroplast marker, rbcL, 78 sequences from nine localities (six in central Chile: Table 1) yielded two haplotypes, with 12 variable sites across the 886 bp fragment. Twenty-seven of the 34 variable COI nucleotide sites were third-codon positions, and transitions accounted for most inferred changes, with transversions inferred at only three sites for COI, and at one site for rbcL. As no genetic diversity for rbcL was observed in 54 samples from six widely-spaced central Chilean localities, obtaining rbcL sequences for all samples for this region was deemed redundant, and only a subset of samples were analysed for this less-informative marker.
Phylogenetic analyses of COI revealed that 13 of the D. antarctica haplotypes in Chile (C-II to C-XIV) were restricted to populations in central Chile, whereas all sampled Chilean Patagonian populations south of 44°S were fixed for a single haplotype (C-I), which was highly divergent from the others (Fig. 2). The widespread Patagonian haplotype, C-I, was previously detected at numerous subantarctic locations throughout the Southern Hemisphere  (Fig. 3). With the exception of the Patagonian haplotype (C-I), all (central) Chilean haplotypes (C-II - C-XIV) were closely related to one another (uncorrected distances 0.2 - 1.3%), together forming a strongly supported clade (Bayesian PP 1.00; Fig. 2a); by contrast, the Patagonian haplotype (C-I) was markedly divergent from all other Chilean haplotypes (C-II - C-XIV) (uncorrected distances 4.5 - 6.1%). The two chloroplast (rbcL) haplotypes detected in Chile were 1.4% divergent (uncorrected distance), with one (R-I) restricted to all sampled Patagonian populations, and the other (R-II) found only in central Chile (Cucao, Punta Loncoyen, Tumbes, Constitución, Montemar and Pichicuy). More broadly, the Patagonian rbcL haplotype R-I has a circum-subantarctic distribution, including the Falkland Islands and South Georgia (not shown; ).
Network analyses of Chilean COI haplotypes recovered two parsimonious networks that could not be joined at ≥90% confidence limits, with one comprising only the Patagonian haplotype C-I. All other Chilean haplotypes (CII-CXIV) formed a single network (Fig. 3 shows these networks in the context of previously-published work ).
The Mantel test revealed a significant relationship between genetic and geographic distance (P = 0.015 for 999 permutations) (Fig. 4). Detailed geographic analyses based on maximum uninterrupted beach length between pairs of adjacent sampling locations in central Chile revealed a clear relationship between rocky shore habitat continuity and genetic connectivity along the central Chilean coast (Fig. 5). Under the definition of genetic disjunction corresponding to Nei's D ≥1.0, no 'disjunct' locality pairs shared any haplotypes. Notably, all phylogeographic disjunctions in central Chilean D. antarctica were associated with extensive (> 20 km) beaches (Fig. 5). The only adjacent locality pair separated by a beach > 20 km long that did not show a genetic disjunction was Lebu/Tirua (locality pair J in Fig. 5), separated by 62 km of uninterrupted beach. Logistic regression (with the binomial response variable of 'genetically disjunct or not', and the continuous fixed effects 'maximum beach length between adjacent localities' and 'coastal distance between adjacent localities') revealed that maximum beach length was a significant predictor of genetic disjunction (β = 0.05; P = 0.035; N = 15), whereas total coastal distance was not (β = 0.01; P = 0.608; N = 15). The outlier formed by Lebu/Tirua (locality pair J, Fig. 5) is, however, likely having a strong effect on the analysis; indeed, removing this locality pair from the analyses greatly strengthens the effect of long beaches in the logistic regression (β = 0.23; P = 0.002; N = 14). With or without the outlier, these results indicate that habitat discontinuity influences genetic connectivity among populations of central Chilean bull-kelp.
Two broad findings were revealed by our phylogeographic research of D. antarctica in Chile. First, substantial phylogeographic structure was observed in central Chile (hypothesis i). Genetic isolation among populations in this region appears driven by habitat discontinuities. Second, and in direct contrast to the pattern observed for central Chile, samples from Chilean Patagonia showed complete genetic homogeneity within and among localities (hypothesis ii). These Patagonian populations were deeply divergent to those from more northern (central Chilean) localities, and appear to have trans-oceanic ancestry.
Habitat discontinuity and population connectivity in central Chile
The substantial phylogeographic structure observed among central Chilean sampling localities of D. antarctica (Fig. 2) suggests a low level of genetic connectivity across this region (hypothesis i), a finding supported by the significant results of the isolation-by-distance test. Furthermore, the strong relationship between long beaches and genetic disjunctions in D. antarctica (Fig. 5) suggests that habitat discontinuity is likely a major factor controlling local genetic connectivity in this rocky-shore macroalga. Population connectivity of other kelp species may be similarly affected by interruptions in habitat: Faugeron et al. , for example, demonstrated genetic differentiation among populations of the kelp Lessonia nigrescens in northern Chile across a habitat discontinuity formed by a mine-waste disposal site. In the English Channel, Billot et al.  found evidence for reduced gene flow among populations of the kelp Laminaria digitata across habitat discontinuities. Patterns of genetic disjunction associated with habitat discontinuity have also been noted for numerous rock-dwelling taxa with predicted high dispersal abilities, including rocky-reef fish [7, 11], rock-pool copepods , and macroalgae (e.g., small-scale patterns of genetic structure in Phyllospora ), suggesting that for marine taxa, as observed for terrestrial plants [30, 31], dispersive life history alone does not necessarily lead to panmixia.
That habitat discontinuity appears to disrupt connectivity among central Chilean populations of D. antarctica is in direct contrast to the broad evidence of long-distance, trans-oceanic dispersal of this species in the subantarctic . Several biological factors may help to explain such a contrast in apparent dispersal ability/effectiveness. First, established and dense kelp populations may resist genetic input from rafted individuals  in a form of density-blocking [20, 25]. In such populations, limited substratum availability and high competition from local D. antarctica zygotes would suppress settlement of zygotes from such migrants. Conversely, shores swept clean by ice  or other large-scale disturbances would offer ample settlement opportunities for zygotes from rafted migrants. Although the genetic effects of density-blocking have yet to be demonstrated experimentally for brown algae, the complete genetic homogeneity of putative postglacially-recolonised populations (this study and ), and the absence of the Patagonian COI haplotype (CI) from central Chilean populations, indicate that such a process is likely to be occurring. Second, although D. antarctica is clearly an effective long-distance disperser, there may be differences in reproductive viability of rafts across latitudinal gradients. Numerous reports suggest that D. antarctica has an extended period of fertility in the subantarctic versus mainland New Zealand or Chile (see references in ), perhaps facilitated by longer optimal reproductive conditions in colder waters. Observations by Macaya et al.  suggest that detached rafts of another kelp, Macrocystis pyrifera, lose reproductive activity at high temperatures, a finding supported by Rothäusler et al. . If D. antarctica rafts are similarly affected by high water temperatures, fewer zygotes from immigrant individuals may be available to colonise central Chile than in the relatively cold waters of Chilean Patagonia (Fig. 1).
Local genetic structure within and among central Chilean D. antarctica populations could additionally be influenced by a range of processes such as seasonal variation and patchy recruitment  or harvesting [34–36]. While the results from Bustamante & Castilla  and Castilla et al.  underline the importance of short-distance dispersal in local population dynamics, coastal islands might serve as stepping stones allowing dispersal across long stretches of sandy beaches. Indeed, the only locality pair separated by a long beach that did not show a genetic disjunction (Lebu-Tirua) is in the immediate vicinity of a large coastal island (Isla Mocha) that may effectively connect the two populations despite the beach barrier. Interestingly, Constitución, flanked by long beaches on both sides, was genetically distinct from adjacent localities but not from some more distant localities (Cartagena, Montemar and Pichicuy). These populations, on the northernmost edge of the range of D. antarctica in Chile, may be sensitive to any changes in water temperature. Changes to the northern range limit of an Australian congeneric, D. potatorum, have been linked to warming waters in recent years . The northern populations of D. antarctica in Chile (Cartagena, Montemar and Pichicuy) may thus be somewhat ephemeral in the long term, and could have been relatively recently colonised by kelp drifting from southern populations (e.g., Constitución) in the north-flowing Humboldt Current.
Major genetic disjunction in southern Chile
Our phylogeographic analyses of Chilean D. antarctica populations indicate a biogeographic break in northern Chilean Patagonia (44° - 49°S), with both mitochondrial and chloroplast DNA phylogenies revealing a deep genetic disjunction between samples from southern Patagonia and those from central Chile (Fig. 2). Although the precise geographic location of this major genetic disjunction cannot be pinpointed with this present dataset, a biogeographic break in the vicinity of the Golfo de Penas/Taitao Peninsula (46.5°S; Fig. 2) has previously been noted [38, 39]. More recently, Lancellotti and Vásquez  suggested this feature may represent the southern boundary of the transitional zone, separating it from the cold-temperate southern waters south of 46°S. The hypothesis of a biogeographic break at Taitao Peninsula has been supported by several recent shallow-marine surveys of Cnidaria in southern Chile [41, 42], with samples showing distinct changes in species composition to the north and south of about 46°S. In a recent study of the density of drifting kelp in channels north of Taitao Peninsula, Hinojosa et al.  suggest that D. antarctica rafts are driven by both surface currents (primarily westward) and wind (primarily eastward) through the fiords. Taitao Peninsula, however, projects across the dense system of channels that could otherwise potentially allow mixing and north-south transport of drifting D. antarctica within the fiords.
The deep genetic disjunction between D. antarctica samples from southern Patagonia and those from central Chile raises the question of the taxonomic status of these two divergent lineages. A recent study of D. antarctica revealed that this taxon comprises two ecologically, morphologically and genetically distinct forms in New Zealand . Given that the genetic disjunction between Patagonian and central Chilean D. antarctica (4.5 - 6.1% for COI) exceeds the divergence observed between the co-occurring, genetically-distinct forms in New Zealand (3.0 - 3.8% for COI: ), it is possible that the distinct Chilean lineages represent reproductively-isolated species, and future studies should assess any physiological, ecological or life history differences between them. Regardless, the phylogeographic contrasts between these lineages are intriguing, with negligible genetic structure observed throughout sampled Patagonian populations (a single haplotype for COI at all five localities across some 1000 km of shoreline), versus high structure further north (eleven COI haplotypes across approximately 1400 km).
Postglacial recolonisation of coastal Chilean Patagonia
The complete genetic homogeneity observed in D. antarctica from all localities south of latitude 44°S (the region covered by the Patagonian Ice Sheet at the LGM: Fig. 1), compared with substantial phylogeographic structure among populations in central Chile for COI (Fig. 2), strongly supports a hypothesis of postglacial recolonisation of Chilean Patagonia (hypothesis ii). Such events are typically marked by low levels of genetic structure in recolonised versus refugial areas (reviewed by [20–24]) although, to date, few studies have been carried out in the Southern Hemisphere . Recolonisation of Patagonia appears to have occurred from a trans-oceanic source, rather than from geographically proximate areas in central Chile. The Patagonian COI haplotype C-I also occurs on numerous subantarctic islands (Fig. 3, and see ), and is more closely related to haplotypes found in the New Zealand subantarctic region than to other Chilean haplotypes , suggesting the New Zealand region as a source, with transport facilitated by the Antarctic Circumpolar Current (Fig. 1) [46–49].
Although genetic evidence of postglacial recolonisation has previously been identified from a variety of Patagonian terrestrial and freshwater taxa, including fish (galaxiids: [50–52]; percichthyids: [53, 54]), lizards , and plants [56–58], this is the first study to reveal evidence of postglacial recolonisation of a Chilean Patagonian marine taxon (but see recent findings on the likely impacts of glacial cycles on marine dispersal by diadromous fish in Patagonia ). Additionally, in comparison with Fraser et al. , who inferred recolonisation of several subantarctic islands by Durvillaea after extirpation by sea ice, this work shows that the glacial (i.e., land-based) ice during the most recent ice age may also have driven kelp population extinctions. If populations of D. antarctica in southern Chile were indeed destroyed by the Patagonian Ice Sheet during the last glacial period, numerous other coastal taxa may have been similarly affected. Intriguingly, a broad-scale phylogeographic study of the rocky-shore gastropod Concholepas concholepas , found no evidence of recent (post-LGM) recolonisation of Chilean Patagonia. Perhaps the habitat of this molluscan species, extending to 40 m below sea level , enabled it to survive throughout the last ice age in regions where more intertidally-restricted taxa were eliminated by ice. Future studies should aim to assess whether a pattern of postglacial recolonisation is common among other subtidal and intertidal species in southern Chile.
Our understanding of evolutionary processes relies largely on our knowledge of the factors that drive and maintain genetic variation and population connectivity. Marine systems can present paradoxical patterns of genetic connectivity, with taxa that apparently have low intrinsic dispersal capacities occasionally showing broader geographic ranges than highly dispersive taxa [61–63]. In such cases, passive dispersal by rafting may explain the broad distributions of otherwise non-dispersive species . The results of the current study suggest that some taxa cannot necessarily be classified as either 'high' or 'low' dispersal organisms, and that mechanisms of dispersal may vary across temporal and geographic scales. Bull-kelp (D. antarctica) is capable of long-distance dispersal via rafting of fertile adults, as evidenced by the trans-oceanic ancestry of populations throughout postglacially-recolonised Chilean Patagonia (and subantarctic islands ). Intriguingly, this dispersal capacity facilitates colonisation of as-yet unoccupied rocky shores, but does not enhance connectivity among established populations. Long-established populations, such as those along the coast of central Chile, seem strongly influenced by the isolating effects of geographic distance and habitat discontinuities, indicating relatively poor inter-population connectivity.
Sample collection, DNA extraction and sequencing
164 samples of D. antarctica were collected from 24 localities along the coast of Chile, from Pichicuy (32°27'37"S) to Cape Horn (55°59'31"S) (Fig. 2, Table 1). These localities encompassed most of the Chilean range of this species, which extends from Los Vilos (31°54'27"S; 60 km north of Pichicuy) to Cape Horn. The data from seven of the localities in Chile (Cucao, Bahia Mansa, Punta Loncoyen, Tumbes, Constitución, Montemar and Pichicuy) were included in previously published analyses  but are here reanalyzed in-depth in the light of extensive additional sampling (106 new samples from 17 new localities). Samples consisted of small (< 2 cm2) pieces of tissue cut from the tip of healthy fronds of attached D. antarctica, and were preserved by desiccation in either 96% ethanol or silica gel beads. DNA extraction, amplification, purification and sequencing were performed as described in Fraser et al. . A 629 bp region of mitochondrial cytochrome c oxidase I (COI) was amplified using primers GazF1 and GazR1  for all samples. In addition, an 886 bp region of chloroplast rbcL was amplified using primers KL2 and KL8  for samples from three localities in Chilean Patagonia, and analysed with previously-published sequences from central Chile (Pichicuy, Montemar, Constitución, Tumbes, Punta Loncoyen and Cucao) (Table 1).
Phylogenies for each DNA dataset were constructed using maximum-likelihood (ML) and Bayesian analyses, and included outgroup sequences from Durvillaea willana (South Island, New Zealand; GenBank accessions: EU918569 for COI and EU918578 for rbcL), Durvillaea potatorum (Tathra, southeastern Australia; GenBank accessions: FJ873092 for COI and FJ872990 for rbcL) and Fucus (GenBank accessions: F. vesiculosus AY494079 for COI and F. gardneri AF195515 for rbcL). To place the new Chilean sequence data in a wider context, published D. antarctica DNA sequences (GenBank accessions FJ550086 - FJ550112) from seven of the localities in central Chile, and numerous localities around the subantarctic and New Zealand, were included in the analyses (Fig. 2). Note that the unique COI and rbcL haplotypes reported from the Crozet Islands by Fraser et al.  were later found to be erroneous, and were in fact identical to the widespread haplotypes C-I and R-I; this error was corrected for the present analyses. ML analyses were performed using PhyML  with a TrN + I  model for both COI and rbcL (proportion of invariable sites estimated, Nst = 6) as selected by the AIC of Modeltest 3.06 . Node support was estimated by bootstrapping , with heuristic analysis of 1000 replicate data sets (Fig. 2). Bayesian posterior probability (PP) values were calculated using MRBAYES 3.1.2  and are shown on the ML phylogenies (Fig. 2). Markov chain Monte Carlo (McMC) searches were performed, with four chains of 5,000,000 generations and trees sampled every 100. The first 10,000 trees sampled were discarded as 'burn-in,' as determined using Tracer v1.4 , while a consensus topology and posterior probability values were calculated with the remaining trees. Unrooted statistical parsimony networks were built using TCS 1.21 . All unique DNA sequences obtained during this study were deposited with GenBank (accession numbers for all Chilean D. antarctica COI and rbcL sequences used in this study: FJ550093-FJ550095; FJ550097; FJ550099; FJ550119; HM103936-HM104173).
Analysis of isolation-by-distance was performed on data from the 19 localities in central Chile (north of 44°S) via the Mantel  test, implemented in GenAlEx6 , using raw (D) average pairwise differences  (calculated using Arlequin version 3.0 ) versus coastal distance (not including deep bays with narrow, < approximately 2 km, openings) among sample localities. Linearised population pairwise FST values could not be used in the Mantel test: several localities were genetically monomorphic for different haplotypes, and pairwise comparisons between such fixed populations gave an FST of 1.0; hence for these locality pairs no linearised FST values (FST/(1 - FST), ) could be calculated.
Approximate measurements of the lengths of all beaches (coastal stretches of uninterrupted sand, gravel or pebbles) between pairs of adjacent rocky coastal localities in central Chile were made using the 'path ruler' tool in Google Earth http://earth.google.com/ to assess the continuity of suitable rocky shore habitat for D. antarctica among localities. As Chilean shores south of about 43°S are almost completely rocky, with few sandy beaches and, moreover, many wide oceanic channels between islands and peninsulas, measurements of beaches were only made for continental Chile and the coast of the Great Island of Chiloé (north of 43°S). The following analysis was performed only on data from sampled localities in this region (Pichichuy, at 32° 20', to Cucao, at 42°40' S). The isolating effects of beach length and total coastal distance between localities in central Chile were assessed by logistic regression. Adjacent locality pairs were first classified binomially as either genetically disjunct (D ≥1.0) or not (D < 1.0). Under this definition, locality pairs not sharing any haplotypes were considered genetically disjunct, whereas those with shared haplotypes were not considered disjunct. Logistic regression was then carried out using R , using the logistf add-on package, with the model: 'genetically disjunct' ~ 'maximum beach length between adjacent localities' + 'coastal distance between adjacent localities'. This analysis enabled assessment of the relative abilities of total coastal distance, versus beach length, to predict whether or not adjacent localities were genetically disjunct.
Caley MJ, Carr MH, Hixon MA, Hughes TP, Jones GP, Menge BA: Recruitment and the local dynamics of open marine populations. Annu Rev Ecol Syst. 1996, 27: 477-500. 10.1146/annurev.ecolsys.27.1.477.
Hellberg ME, Burton RS, Neigel JE, Palumbi SR: Genetic assessment of connectivity among marine populations. Bull Mar Sci. 2002, 70 (1): 273-290.
Palumbi SR: Population genetics, demographic connectivity, and the design of marine reserves. Ecol Appl. 2003, 13 (1): S146-S158. 10.1890/1051-0761(2003)013[0146:PGDCAT]2.0.CO;2.
Levin LA: Recent progress in understanding larval dispersal: new directions and digressions. Integr Comp Biol. 2006, 46 (3): 282-297. 10.1093/icb/icj024.
Cowen RK, Sponaugle S: Larval dispersal and marine population connectivity. Annu Rev Mar Sci. 2009, 1: 443-466. 10.1146/annurev.marine.010908.163757.
Palumbi SR: Genetic divergence, reproductive isolation, and marine speciation. Annu Rev Ecol Syst. 1994, 25: 547-572. 10.1146/annurev.es.25.110194.002555.
Riginos C, Nachman MW: Population subdivision in marine environments: the contributions of biogeography, geographical distance and discontinuous habitat to genetic differentiation in a blennioid fish, Axoclinus nigricaudus. Mol Ecol. 2001, 10 (6): 1439-1453. 10.1046/j.1365-294X.2001.01294.x.
Thiel M, Macaya EC, Acuna E, Arntz WE, Bastias H, Brokordt K, Camus PA, Castilla JC, Castro LR, Cortes M, et al: The Humboldt Current System of northern and central Chile: Oceanographic processes, ecological interactions and socioeconomic feedback. Oceanogr Mar Biol Annu Rev. 2007, 45: 195-344.
Burton RS, Feldman MW: Population genetics of Tigriopus californicus. 2. Differentiation among neighboring populations. Evolution. 1981, 35 (6): 1192-1205. 10.2307/2408132.
Hidas EZ, Costa TL, Ayre DJ, Minchinton TE: Is the species composition of rocky intertidal invertebrates across a biogeographic barrier in south-eastern Australia related to their potential for dispersal?. Marine and Freshwater Research. 2007, 58 (9): 835-842.
Johansson ML, Banks MA, Glunt KD, Hassel-Finnegan HM, Buonaccorsi VP: Influence of habitat discontinuity, geographical distance, and oceanography on fine-scale population genetic structure of copper rockfish (Sebastes caurinus). Mol Ecol. 2008, 17 (13): 3051-3061. 10.1111/j.1365-294X.2008.03814.x.
Coleman MA, Kelaher BP: Connectivity among fragmented populations of a habitat-forming alga, Phyllospora comosa (Phaeophyceae, Fucales) on an urbanised coast. Mar Ecol Prog Ser. 2009, 381: 63-70. 10.3354/meps07977.
Billot C, Engel CR, Rousvoal S, Kloareg B, Valero M: Current patterns, habitat discontinuities and population genetic structure: the case of the kelp Laminaria digitata in the English Channel. Mar Ecol-Prog Ser. 2003, 253: 111-121. 10.3354/meps253111.
Kusumo HT, Pfister CA, Wootton JT: Small-scale genetic structure in the sea palm Postelsia palmaeformis Ruprecht (Phaeophyceae). Mar Biol. 2006, 149 (4): 731-742. 10.1007/s00227-006-0254-z.
Faugeron S, Martinez EA, Correa JA, Billot C: Long-term copper mine waste disposal in northern Chile associated with gene flow disruption of the intertidal kelp Lessonia nigrescens. Mar Ecol Prog Ser. 2005, 288: 129-140. 10.3354/meps288129.
Muhlin JF, Engel CR, Stessel R, Weatherbee RA, Brawley SH: The influence of coastal topography, circulation patterns, and rafting in structuring populations of an intertidal alga. Mol Ecol. 2008, 17 (5): 1198-1210. 10.1111/j.1365-294X.2007.03624.x.
McKenzie PF, Bellgrove A: Dispersal of Hormosira banksii (Phaeophyceae) via detached fragments: Reproductive viability and longevity. J Phycol. 2008, 44 (5): 1108-1115. 10.1111/j.1529-8817.2008.00563.x.
Macaya EC, Boltana S, Hinojosa IA, Macchiavello JE, Valdivia NA, Vasquez NR, Buschmann AH, Vasquez JA, Vega JMA, Thiel M: Presence of sporophylls in floating kelp rafts of Macrocystis spp. (Phaeophyceae) along the Chilean Pacific coast. J Phycol. 2005, 41 (5): 913-922. 10.1111/j.1529-8817.2005.00118.x.
Fraser CI, Nikula R, Spencer HG, Waters JM: Kelp genes reveal effects of subantarctic sea ice during the Last Glacial Maximum. Proc Natl Acad Sci USA. 2009, 106 (9): 3249-3253. 10.1073/pnas.0810635106.
Hewitt GM: Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Linn Soc. 1996, 58 (3): 247-276.
Hewitt G: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405 (6789): 907-913. 10.1038/35016000.
Hewitt GM: Genetic consequences of climatic oscillations in the Quaternary. Philos Trans R Soc London, Ser B. 2004, 359 (1442): 183-195. 10.1098/rstb.2003.1388.
Provan J, Bennett KD: Phylogeographic insights into cryptic glacial refugia. Trends Ecol Evol. 2008, 23 (10): 564-571. 10.1016/j.tree.2008.06.010.
Maggs CA, Castilho R, Foltz D, Henzler C, Jolly MT, Kelly J, Olsen J, Perez KE, Stam W, Vainola R, et al: Evaluating signatures of glacial refugia for North Atlantic benthic marine taxa. Ecology. 2008, 89 (11): S108-S122. 10.1890/08-0257.1.
Hewitt GM: Using molecules to understand the distribution of animal and plant diversity. Evolution: from molecules to ecosystems. Edited by: Moya A, Font E. 2004, Oxford: Oxford University Press, 157-170.
Harrison S: The Pleistocene glaciations of Chile. Quaternary Glaciations: Extent and Chronology. Edited by: Ehlers J, Gibbard PL. 2004, Boston: Elsevier, III:
McCulloch RD, Bentley MJ, Purves RS, Hulton NRJ, Sugden DE, Clapperton CM: Climatic inferences from glacial and palaeoecological evidence at the last glacial termination, southern South America. Journal of Quaternary Science. 2000, 15 (4): 409-417. 10.1002/1099-1417(200005)15:4<409::AID-JQS539>3.0.CO;2-#.
Hulton NRJ, Purves RS, McCulloch RD, Sugden DE, Bentley MJ: The Last Glacial Maximum and deglaciation in southern South America. Quaternary Sci Rev. 2002, 21 (1-3): 233-241. 10.1016/S0277-3791(01)00103-2.
Hay C: Durvillaea (Bory). Biology of Economic Algae. Edited by: Akatsuka I. 1994, The Hague: SPB Academic Publishing, 353-384. p 545
Cain ML, Milligan BG, Strand AE: Long-distance seed dispersal in plant populations. Am J Bot. 2000, 87 (9): 1217-1227. 10.2307/2656714.
Clark JS: Why trees migrate so fast: Confronting theory with dispersal biology and the paleorecord. Am Nat. 1998, 152 (2): 204-224. 10.1086/286162.
Rothäusler E, Gomez I, Hinojosa IA, Karsten U, Tala F, Thiel M: Effect of temperature and grazing on growth and reproduction of floating Macrocystis spp. (Phaeophyceae) along a latitudinal gradient. J Phycol. 2009, 45 (3): 547-559. 10.1111/j.1529-8817.2009.00676.x.
Taylor DI, Schiel DR: Self-replacement and community modification by the southern bull kelp Durvillaea antarctica. Mar Ecol Prog Ser. 2005, 288: 87-102. 10.3354/meps288087.
Castilla JC, Bustamante RH: Human exclusion from rocky intertidal of Las Cruces, central Chile - effects on Durvillaea antarctica (Phaeophyta, Durvilleales). Mar Ecol Prog Ser. 1989, 50 (3): 203-214. 10.3354/meps050203.
Bustamante RH, Castilla JC: Impact of human exploitation on populations of the intertidal southern bull-kelp Durvillaea antarctica (Phaeophyta, Durvilleales) in central Chile. Biol Conserv. 1990, 52 (3): 205-220. 10.1016/0006-3207(90)90126-A.
Castilla JC, Campo MA, Bustamante RH: Recovery of Durvillaea antarctica (Durvilleales) inside and outside Las Cruces marine reserve, Chile. Ecol Appl. 2007, 17 (5): 1511-1522. 10.1890/06-1285.1.
Millar AJK: The Flindersian and Peronian Provinces. Algae of Australia: Introduction. Edited by: McCarthy PM, Orchard AE. 2007, Melbourne: CSIRO Publishing, 554-559.
Pickard GL: Some physical oceanographic features of inlets of Chile. J Fish Res Board Can. 1971, 28 (8): 1077-1106.
Viviani CA: Ecografía del litoral chileno. Studies on Neotropical Fauna and Environment. 1979, 14: 65-123. 10.1080/01650527909360548.
Lancellotti DA, Vasquez JA: Biogeographical patterns of benthic macroinvertebrates in the Southeastern Pacific littoral. J Biogeogr. 1999, 26 (5): 1001-1006. 10.1046/j.1365-2699.1999.00344.x.
Häussermann V: Biodiversity of Chilean sea anemones (Cnidaria: Anthozoa): distribution patterns and zoogeographic implications, including new records for the fjord region. Invest Mar. 2006, 34: 23-35.
Galea HR, Häussermann V, Försterra G: Cnidaria, Hydrozoa: latitudinal distribution of hydroids along the fjords region of southern Chile, with notes on the world distribution of some species. Check List. 2007, 3: 308-320.
Hinojosa IA, Pizarro M, Ramos M, Thiel M: Spatial and temporal distribution of floating kelp in the channels and fjords of southern Chile. Estuar Coast Shelf Sci. 2010, 87 (3): 367-377. 10.1016/j.ecss.2009.12.010.
Fraser CI, Hay CH, Spencer HG, Waters JM: Genetic and morphological analyses of the southern bull kelp Durvillaea antarctica (Phaeophyceae: Durvillaeales) in New Zealand reveal cryptic species. J Phycol. 2009, 45 (2): 436-443. 10.1111/j.1529-8817.2009.00658.x.
Beheregaray LB: Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere. Mol Ecol. 2008, 17 (17): 3754-3774.
Pequeno G, Lamilla J: The littoral fish assemblage of the Desventuradas Islands (Chile) has zoogeographical affinities with the Western Pacific. Global Ecol Biogeogr. 2000, 9 (5): 431-437. 10.1046/j.1365-2699.2000.00207.x.
Smith SDA: Kelp rafts in the Southern Ocean. Global Ecol Biogeogr. 2002, 11 (1): 67-69. 10.1046/j.1466-822X.2001.00259.x.
Donald KM, Kennedy M, Spencer HG: Cladogenesis as the result of long-distance rafting events in South Pacific topshells (Gastropoda, Trochidae). Evolution. 2005, 59 (8): 1701-1711.
Burridge CP, Melendez R, Dyer BS: Multiple origins of the Juan Fernandez kelpfish fauna and evidence for frequent and unidirectional dispersal of cirrhitoid fishes across the south Pacific. Syst Biol. 2006, 55 (4): 566-578. 10.1080/10635150600812585.
Cussac V, Ortubay S, Iglesias G, Milano D, Lattuca ME, Barriga JP, Battini M, Gross M: The distribution of South American galaxiid fishes: the role of biological traits and post-glacial history. J Biogeogr. 2004, 31 (1): 103-121.
Zemlak TS, Habit EM, Walde SJ, Carrea C, Ruzzante DE: Surviving historical Patagonian landscapes and climate: molecular insights from Galaxias maculatus. BMC Evol Biol. 2010, 10: 67-10.1186/1471-2148-10-67.
Zemlak TS, Habit EM, Walde SJ, Battini MA, Adams EDM, Ruzzante DE: Across the southern Andes on fin: glacial refugia, drainage reversals and a secondary contact zone revealed by the phylogeographical signal of Galaxias platei in Patagonia. Mol Ecol. 2008, 17 (23): 5049-5061. 10.1111/j.1365-294X.2008.03987.x.
Ruzzante DE, Walde SJ, Cussac VE, Dalebout ML, Seibert J, Ortubay S, Habit E: Phylogeography of the Percichthyidae (Pisces) in Patagonia: roles of orogeny, glaciation, and volcanism. Mol Ecol. 2006, 15 (10): 2949-2968. 10.1111/j.1365-294X.2006.03010.x.
Ruzzante DE, Walde SJ, Gosse JC, Cussac VE, Habit E, Zemlak TS, Adams EDM: Climate control on ancestral population dynamics: insight from Patagonian fish phylogeography. Mol Ecol. 2008, 17 (9): 2234-2244. 10.1111/j.1365-294X.2008.03738.x.
Victoriano PF, Ortiz JC, Benavides E, Adams BJ, Sites JW: Comparative phylogeography of codistributed species of Chilean Liolaemus (Squamata: Tropiduridae) from the central-southern Andean range. Mol Ecol. 2008, 17 (10): 2397-2416. 10.1111/j.1365-294X.2008.03741.x.
Allnutt TR, Newton AC, Lara A, Premoli A, Armesto JJ, Vergara R, Gardner M: Genetic variation in Fitzroya cupressoides (alerce), a threatened South American conifer. Mol Ecol. 1999, 8 (6): 975-987. 10.1046/j.1365-294x.1999.00650.x.
Pastorino MJ, Gallo LA: Quaternary evolutionary history of Austrocedrus chilensis, a cypress native to the Andean-Patagonian forest. J Biogeogr. 2002, 29 (9): 1167-1178. 10.1046/j.1365-2699.2002.00731.x.
Muellner AN, Tremetsberger K, Stuessy T, Baeza CM: Pleistocene refugia and recolonization routes in the southern Andes: insights from Hypochaeris palustris (Asteraceae, Lactuceae). Mol Ecol. 2005, 14 (1): 203-212. 10.1111/j.1365-294X.2004.02386.x.
Cardenas L, Castilla JC, Viard F: A phylogeographical analysis across three biogeographical provinces of the south-eastern Pacific: the case of the marine gastropod Concholepas concholepas. J Biogeogr. 2009, 36 (5): 969-981. 10.1111/j.1365-2699.2008.02056.x.
Poulin E, Palma AT, Leiva G, Narvaez D, Pacheco R, Navarrete SA, Castilla JC: Avoiding offshore transport of competent larvae during upwelling events: The case of the gastropod Concholepas concholepas in Central Chile. Limnol Oceanogr. 2002, 47 (4): 1248-1255. 10.4319/lo.2002.47.4.1248.
Johannesson K: The paradox of Rockall - why is a brooding gastropod (Littorina saxatilis) more widespread than one having a planktonic larval dispersal stage (L. littorea)?. Mar Biol. 1988, 99 (4): 507-513. 10.1007/BF00392558.
Sponer R, Roy MS: Phylogeographic analysis of the brooding brittle star Amphipholis squamata (Echinodermata) along the coast of New Zealand reveals high cryptic genetic variation and cryptic dispersal potential. Evolution. 2002, 56 (10): 1954-1967.
Waters JM, Roy MS: Out of Africa: The slow train to Australasia. Syst Biol. 2004, 53 (1): 18-24. 10.1080/10635150490264671.
Thiel M, Haye PA: The ecology of rafting in the marine environment. III. Biogeographical and evolutionary consequences. Oceanogr Mar Biol Annu Rev. 2006, 44: 323-429.
Saunders GW: Applying DNA barcoding to red macroalgae: a preliminary appraisal holds promise for future applications. Philos Trans R Soc London, Ser B. 2005, 360 (1462): 1879-1888. 10.1098/rstb.2005.1719.
Lane CE, Mayes C, Druehl LD, Saunders GW: A multi-gene molecular investigation of the kelp (Laminariales, Phaeophyceae) supports substantial taxonomic re-organization. J Phycol. 2006, 42 (2): 493-512. 10.1111/j.1529-8817.2006.00204.x.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.
Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10 (3): 512-526.
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14 (9): 817-818. 10.1093/bioinformatics/14.9.817.
Felsenstein J: Confidence limits on phylogenies - an approach using the bootstrap. Evolution. 1985, 39 (4): 783-791. 10.2307/2408678.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.
Rambaut A, Drummond AJ: Tracer v.1.4 [computer program]. 2007, [http://beast.bio.ed.ac.uk/tracer]
Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9 (10): 1657-1659. 10.1046/j.1365-294x.2000.01020.x.
Mantel N: Detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27: 209-220.
Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6 (1): 288-295. 10.1111/j.1471-8286.2005.01155.x.
Nei M, Li WH: Mathematical model for studying genetic variation in terms of restriction endonucleases. P Natl Acad Sci USA. 1979, 76 (10): 5269-5273. 10.1073/pnas.76.10.5269.
Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: an integrated software for molecular evolutionary genetic analysis. Evolutionary Bioinformatics Online. 2005, 1: 47-50.
Slatkin M: A measure of population subdivision based on microsatellite allele frequencies. Genetics. 1995, 139 (1): 457-462.
R Development Core Team: R: A language and environment for statistical computing v.[computer program]. 2009, [http://www.R-project.org]
Levitus S, Boyer TP: World Ocean Atlas, Temperature. 1994, Washington, DC: US Govt. Printing Office, 4:
Kaiser J, Lamy F, Hebbeln D: A 70-kyr sea surface temperature record off southern Chile (Ocean Drilling Program Site 1233). Paleoceanography. 2005, 20: PA4009-10.1029/2005PA001146.
We are extremely grateful to all who assisted with this research. For field support/sample collection, we thank: Emma Newcombe and Carlos Olavarría (Centro de Estudios del Cuaternario, Punta Arenas, Chile) for advice and collections from Cape Horn and southern Patagonia (Cruceros Australis and Whalesound provided suport for these collections); David Yáñez Jaramillo (Marine Botany Laboratory, Universidad Católica del Norte, Coquimbo, Chile) for assistance with fieldwork; Iván Hinojosa (Universidad Católica del Norte, Coquimbo, Chile), Patricio Díaz and Hector Pavez (Universidad Austral de Chile, Valdivia, Chile) for collections from Patagonia; Erasmo Macaya Horta (Victoria University of Wellington, New Zealand) for facilitating many collections; Thanks to Tania King for laboratory assistance, to Raisa Nikula for advice, and to Shinichi Nakagawa and Catherine Grueber for help with statistical analyses. Comments from Graham Wallis and two anonymous reviewers also helped improve the MS. This work was funded by Marsden contract 07-UOO-099, Department of Zoology and University of Otago Research grants to JMW; a Department of Zoology funding allocation to CIF; and Allan Wilson Centre for Molecular Ecology and Evolution funds to HGS.
C.I.F., J.M.W., M.T. and H.G.S. conceived the ideas; C.I.F. obtained and analysed the data; C.I.F. led the writing; and J.M.W., M.T. and H.G.S assisted with the writing. All authors read and approved the final manuscript.