Research article | Open | Published:
Towards a model of postglacial biogeography in shallow marine species along the Patagonian Province: lessons from the limpet Nacella magellanica(Gmelin, 1791)
BMC Evolutionary Biologyvolume 12, Article number: 139 (2012)
Patagonia extends for more than 84,000 km of irregular coasts is an area especially apt to evaluate how historic and contemporary processes influence the distribution and connectivity of shallow marine benthic organisms. The true limpet Nacella magellanica has a wide distribution in this province and represents a suitable model to infer the Quaternary glacial legacy on marine benthic organisms. This species inhabits ice-free rocky ecosystems, has a narrow bathymetric range and consequently should have been severely affected by recurrent glacial cycles during the Quaternary. We performed phylogeographic and demographic analyses of N. magellanica from 14 localities along its distribution in Pacific Patagonia, Atlantic Patagonia, and the Falkland/Malvinas Islands.
Mitochondrial (COI) DNA analyses of 357 individuals of N. magellanica revealed an absence of genetic differentiation in the species with a single genetic unit along Pacific Patagonia. However, we detected significant genetic differences among three main groups named Pacific Patagonia, Atlantic Patagonia and Falkland/Malvinas Islands. Migration rate estimations indicated asymmetrical gene flow, primarily from Pacific Patagonia to Atlantic Patagonia (Nem=2.21) and the Falkland/Malvinas Islands (Nem=16.6). Demographic reconstruction in Pacific Patagonia suggests a recent recolonization process (< 10 ka) supported by neutrality tests, mismatch distribution and the median-joining haplotype genealogy.
Absence of genetic structure, a single dominant haplotype, lack of correlation between geographic and genetic distance, high estimated migration rates and the signal of recent demographic growth represent a large body of evidence supporting the hypothesis of rapid postglacial expansion in this species in Pacific Patagonia. This expansion could have been sustained by larval dispersal following the main current system in this area. Lower levels of genetic diversity in inland sea areas suggest that fjords and channels represent the areas most recently colonized by the species. Hence recolonization seems to follow a west to east direction to areas that were progressively deglaciated. Significant genetic differences among Pacific, Atlantic and Falkland/Malvinas Islands populations may be also explained through disparities in their respective glaciological and geological histories. The Falkland/Malvinas Islands, more than representing a glacial refugium for the species, seems to constitute a sink area considering the strong asymmetric gene flow detected from Pacific to Atlantic sectors. These results suggest that historical and contemporary processes represent the main factors shaping the modern biogeography of most shallow marine benthic invertebrates inhabiting the Patagonian Province.
Climatic changes are considered as one of the main factors regulating abundance, composition, and distribution of species at different temporal and spatial scales [1–3]. Direct historical evidence from fossil records indicates that many terrestrial species underwent rapid latitudinal shifts during the Quaternary glacial period and especially after the Last Glacial Maximum (LGM) around 23 to 18 ka [4–6]. Paleontological and palynological records from the Northern Hemisphere, along with biogeographic evidence, provided the empirical basis for the Expansion-Contraction (EC) model of Pleistocene biogeography  which describes the response of populations and species to climatic changes [3, 4, 6, 8, 9]. Under a basic EC model, cool-temperate species from the Northern Hemisphere survived the LGM at lower latitude refugia and then recolonized higher latitudes through range expansion [6, 7].
The application of molecular-based approaches in population genetics has provided new insights into the history of many species and helped us to further understand the impact of the Quaternary glacial cycles on patterns of genetic variation and structure [2, 3, 6, 7]. However, most examples come from studies of Northern Hemisphere biota mainly of terrestrial species [10, 11]. In marine ecosystems, interglacial period deposits are rich in fossils but glacial records are in most cases unavailable due to the Holocene rise in sea level . Phylogeographical studies in non-tropical areas of the Southern Hemisphere are scarce , but during the last decade more data has been accumulated in different southern South American groups [13–20]. Periodic global cooling during Quaternary glacial cycles (1.8 Ma – 10 ka) generated shifts in climate, landscape and sea level. For instance, during the LGM an ice sheet about 1800 km long covered the west slope of the Andes from 35°S to almost 56°S [21–24]. Much of the Atlantic side of Patagonia and northeastern Tierra del Fuego remained unglaciated through the late Pleistocene . These glacial changes in Patagonia led to regional isolations and local extinctions, shaping the current patterns of species diversity in temperate areas of southern South America [13, 16, 25, 26]. Genetic evidence of postglacial recolonization has been found in several Patagonian groups including galaxiid [16, 26, 27] and percichthyd fishes [13, 14, 28], lizards [29–31], amphibians , mammals [11, 20, 33–35] and plants [15, 18, 36]. These studies have provided conflicting results, indicating either postglacial colonization from restricted glacial refugia [26, 32, 36], recolonization from geographically distant ice-free regions [15, 34], or local persistence through glacial cycles [16, 28, 31, 33, 36]. Few genetic studies have examined the effect of the Quaternary glacial cycles in marine organisms of southern South America and most of these were restricted to Pacific sectors of Patagonia [15, 17, 37–39]. Moreover, due to logistic problems especially in the hard-to-access fiordal region of Chilean Patagonia, most of these studies present unbalanced sampling, only including localities from easy-access areas which represent the northern (Reloncaví Archipelago and Chiloé Island) and the southern extremes of Patagonian species distributions (Magellan Strait and Tierra del Fuego).
The true limpet genus Nacella (Patellogastropoda: Nacellidae) includes 15 nominal species distributed in different biogeographical provinces of the Southern Ocean . Along the Patagonian coast, Nacella represent a dominant group of benthic macro-invertebrates, especially in the marine rocky ecosystems [41–44]. Based on morphological characters, at least eight nominal species of the genus have been described in this region (Nacella chiloensis, N. deaurata, N. delicatissima, N. flammea, N. fuegiensis, N. magellanica, N. mytilina, and N. venosa; [40, 43]. On the basis of species richness, Powell  considered Patagonia as the center of origin and diversification of Nacella, from where it expanded eastward through the West Wind Drift (WWD). However, this assumption has been recently rejected by phylogenetic reconstructions showing that the Patagonian group of Nacella is the most derived one and diversified no more than 2.0 Ma . Molecular and morphological comparisons of Patagonian species suggest that the number of nominal species in Nacella was overestimated [17, 46]. For instance, González-Wevar et al. using COI sequences and geometric morphometrics in seven sympatric nominal species recognized only four units of Nacella in Patagonia. In spite of the absence of reciprocal monophyly , morphological, genetic and habitat preference differentiation among congeners are maintained even in sympatry. Considering these results, the diversification of Nacella in Patagonia includes four Evolutionarily Significant Units (ESUs): N. deaurata, N. flammea, N. mytilina and N. magellanica.
Nacella magellanica (Gmelin, 1791) exhibits the widest distribution in Patagonia, extending from Puerto Montt in the Pacific (42°S) to the Buenos Aires Province in the Atlantic (35°–40°S), including the Strait of Magellan, Cape Horn, Tierra del Fuego and the Falkland/Malvinas Islands [40, 43]. This species is the most abundant and conspicuous limpet in intertidal and shallow subtidal areas of Patagonia . It has been also reported in the Beagle Channel as an organism associated with holdfasts of the macroalga Macrocystis pyrifera. As in other nacellid limpets, N. magellanica is a broadcast spawner that reproduces during austral spring  but its free-living larval duration is still unknown. A recent phylogeographic study of the species in Atlantic Patagonia identified an absence of genetic structure and a very recent geographic-demographic expansion (~ 11 ka) . Nevertheless, there is still an absence of knowledge about the patterns of genetic diversity, structure and connectivity of the species in Pacific Patagonia.
The southern tip of South America constitutes an interesting system to evaluate the relative effects of habitat discontinuities, oceanography, and glaciological history in marine benthic organisms with limited autonomous motility that exhibit some degree of larval dispersal. The presence of an extensive ice sheet during the LGM in this region likely eradicated many populations of shallow-water marine benthic organisms. Considering the current distribution of N. magellanica, its narrow bathymetric range and its high abundance along both sides of Patagonia, this species constitutes a suitable model to infer how historical and contemporary climatic events shaped the patterns of population genetic diversity and structure. We analyzed samples from a total of 14 localities encompassing most of the species range in Pacific Patagonia and Tierra del Fuego, as well as two population from Atlantic Patagonia and individuals from the Falkland/Malvinas Islands. We aimed to test the hypothesis that (i) N. magellanica in the Pacific sector of Patagonia represents a post-glacial recolonization from restricted glacial refugia in the northern limit of its distribution, or alternatively, (ii) N. magellanica persisted unaffected through glacial cycles in this area. Also, we aimed to determine if glacial-interglacial periods promoted genetic differentiation or even divergence between Atlantic and Pacific populations. Finally, considering the glaciological history of Patagonia and the pattern of genetic diversity and structure in the species it will be possible to evaluate the role of the Falkland/Malvinas Islands in the phylogeography of the species as a source or sink area.
We analyzed a total of 357 individuals; the COI sequence data set consisted of 671 nucleotide positions coding 223 aminoacids. As expected for coding regions, no indels or stop codons were detected, sequences were not saturated at any position and no amino acid substitution was detected using the invertebrate mitochondrial table . In the whole data set, Nacella magellanica exhibited intermediate levels of genetic diversity with 58 polymorphic characters (8.6%); 37 of them (5.5%) were parsimoniously informative. As previously estimated for nacellids [17, 51], N. magellanica sequences were A-T rich (61.6%) compared to mean G-C content (39.4%). Haplotype diversity (H) varied from 0.370 (Costa Channel) to 0.872 (Falkland/Malvinas Islands). The number of haplotypes and polymorphic sites per locality ranged from 4 (Costa Channel) to 14 (Puerto Montt) and from 5 (Costa Channel) to 13 (Puerto Montt), respectively (Table 1). However, rarefaction analysis of the number of haplotypes using PAST  showed that most of the variations were the result of different sampling sizes, particularly in the case of Puerto Montt. The average number of nucleotide differences (Π) and the nucleotide diversity (π) were low in most of the localities with the exception of the Falkland/Malvinas Islands (Table 1).
The median-joining network depicted from the COI data set exhibited 56 different haplotypes (Figure 1). In the Pacific sector we observed a typical star-like topology in which the central haplotype (H1) was the most frequent (> 50%) and widely distributed. As proposed by Posada & Crandall  this haplotype should represent the most ancestral one, whereas the most derived ones are related to it with a maximum branch length of twelve mutational steps (H55 and H52). Two haplotypes H3 and H42, located no more than two mutational steps away from H1, were present in several localities and showed intermediate frequencies (H3=9.5% and H42=5.3%; Additional files 1 and 2). We did not detect an association between haplotype identity and geographical locality in Pacific Patagonia. Several haplotypes (H2, H8, H12, H15, H26, and H30) were present in more than five individuals belonging to different localities. The remaining haplotypes occurred at low frequencies and we identified 29 singletons. However, in Puerto Deseado H41 and H42 were the most common haplotypes, while H1 (the most frequent in Pacific Patagonia) was present in a single individual (Figure 1). We found some degree of haplotype similarity of haplotype frequencies among localities from the southern tip of the Pacific and Atlantic Patagonia. For instance, none of the dominant haplotypes in Puerto Deseado (H41 and H42) were found in the Reloncavi Fjord, Chonos Archipelago, and Pacific localities in the Strait of Magellan but these haplotypes were observed in Cape Horn sites and in Possession Bay at the eastern mouth of the Strait of Magellan. Finally, in the Falkland/Malvinas Islands we detected two groups of haplotypes. The first includes haplotypes closely related to (H54 and H56) and even shared with (H1 and H26) Patagonian diversity, while the second consists of private haplotypes, separated from the others by approximately 10 mutational steps (Figure 1).
SAMOVA analyses in the whole data set recovered two maximally differentiated groups explaining 51.73% of the total variation. These groups were 1) the Falkland/Malvinas Islands and 2) the rest of the localities along Pacific and Atlantic Patagonia. In a second SAMOVA analysis excluding the Falkland/Malvinas Island sample it recognized two maximally differentiated groups explaining 37.63% of the total variation named: a) Pacific Patagonia (including Reloncaví Fjord, Chonos Archipelago, the Strait of Magellan and Cape Horn localities), and b) Puerto Deseado (Atlantic Patagonia). Finally, in a third SAMOVA analysis including only Pacific and Cape Horn localities it did not recover significant spatial structure in the species. Evidence of this are the levels of variance among groups that explained just between 3.3% and 4.7% of the total variation, while within localities differences represented 95% to 97.5%. According to this, we recognized three main groups named: a) Pacific Patagonia including Reloncaví Fjord, Chonos Archipelago, Strait of Magellan and Cape Horn localities, b) Atlantic Patagonia including Puerto Deseado and c) the Falkland/Malvinas Islands. General differentiation coefficient measured over 14 populations of N. magellanica was low, especially taking into account average GST=0.126 and NST=0.190. Even when some pairwise comparisons between Pacific Patagonia localities showed significant levels of genetic structure (Table 2), none of them were statistically significant after Bonferroni correction. The permutation test indicated that NST is significantly higher than GST (P < 0.05), pointing to a phylogeographical structure for N. magellanica mtDNA haplotypes between the Falkland/Malvinas Islands and the rest of the localities.
Migration rates among SAMOVA’s defined groups (Pacific Patagonia, Atlantic Patagonia and Falkland/Malvinas Islands) showed clear evidence of asymmetrical gene flow. The total number of immigrants per generation (Nem) from the Pacific to the Atlantic was 16.6, while from the Pacific to the Falkland/Malvinas Islands gene flow was lower, Nem=2.21 (Figure 2). In contrast, gene flow from the Atlantic to the Pacific was low (Nem=0.09) and from the Atlantic to the Falkland/Malvinas Islands even lower (Nem=0.0004). Similarly, the migration rate from the Falkland/Malvinas Islands to the Pacific was 0.26 and from the Falkland/Malvinas Islands to Atlantic Patagonia was extremely low, 0.0002 (Figure 2). We detected a small but significant correlation between genetic and geographic distances when all the analyzed localities were included (r=0.37; P < 0.001). This result is expected, considering that most of the significant pairwise comparisons (GST=69% and NST=71%) included Puerto Deseado (Atlantic Patagonia) and Falkland/Malvinas Islands (Table 2). Considering this, we performed a new Mantel test including only Pacific Patagonia localities and the analysis did not detect significant correlation (r=0.14; P=0.10) between geographic and genetic distance from the Reloncaví Fjord to Cape Horn.
Tajima’s D and Fu’s F S neutrality tests showed contrasting results among the three defined groups in N. magellanica of Patagonia. These indices were negative and highly significant in Pacific Patagonia, while in Atlantic Patagonia they were negative but not significant. Finally, in the Falkland/Malvinas Islands both indices were positive and not significant, pointing to different demographic histories among the analyzed sectors. Similarly, the distribution of pairwise differences varied considerably among the recognized genetic groups in N. magellanica. For instance, the mismatch distribution in Pacific Patagonia was L-shaped (Figure 2), in Atlantic Patagonia it showed a bimodal pattern (Figure 2) and in the Falkland/Malvinas Islands it had a multimodal distribution (Figure 2).
Recent molecular studies recognized an error in the use of substitution rates inferred from phylogenetic analyses in studies at the population level [54–57]. It has been demonstrated in different groups of organisms that short-term mutation rates may be tenfold higher than long-term rates [56–58]. Including this tenfold correction to the specific molecular rate estimated for nacellids  and used in another study in the species , the start of the expansion in N. magellanica under a sudden growth model occurred ~ 6.3 ka. The Bayesian skyline plot analysis indicates that the most common recent ancestor of the current N. magellanica’s diversity occurred ~ 24.1 ka while population expansion began around 9 ka (Figure 3).
Understanding how ecosystems and species respond to climate change has become a major focus of ecology and conservation biology [59–62]. One of the central premises of biogeography is that climate exerts a dominant control over the distribution of species [63, 64]. Evidence from the fossil record [65, 66] and from reconstructions based on molecular data [59, 67, 68] have demonstrated that changing climate generates a profound influence on species’ range expansion and contraction, as well as in their patterns of genetic diversity and structure [3, 6, 7, 10, 11, 68, 69, 70, 71]. In the particular case of Nacella magellanica it may be possible to ascribe the observed patterns of genetic diversity and structure to drastic demographic effects of the glacial cycles on the species in its distribution in Patagonia. For instance, COI diversity in the species are lower than those observed in temperate patellogastropods [72–74] but higher than in its Antarctic relative, N. concinna. These results agree with molecular studies in Northern Hemisphere biota where the impact of the Quaternary glacial cycles exerted stronger effects in the demography of marine benthic populations at higher latitudes and particularly in polar regions .
Genetic homogeneity and recent population expansion in N. magellanicain Pacific Patagonia
According to Camus , the Chilean coast of Pacific Patagonia can be considered as a major insular system that includes many islands, gulfs, peninsulas, fjords, and channels that generate a very complex landscape as a result of the marked climatic changes during the Quaternary. Despite such a complex landscape, we found a single genetic unit in N. magellanica from the Reloncaví Fjord (41.5°S) to Cape Horn (55.9°S), an area that includes ~ 1300 km in a straight line but about 84,000 km of irregular coasts [77–80]. For marine benthic organisms, duration of planktonic larval stages is expected to correlate with dispersal ability [81, 82]. Regretfully, there is no direct information about larval duration in the species, but it is expected that its development should be similar to the Antarctic limpet, with a free-swimming planktonic period for 1 to 2 months in the water column [83, 84]. Consequently, it could be expected that the N. magellanica larval period extends for at least four weeks, considering the effect of temperature on development and metabolism [85, 86]. Gene flow mediated by larval dispersal may have been enhanced by oceanographic conditions in the area and constitutes a suitable explanation for the low levels of genetic diversity and the high degree homogeneity in the N. magellanica populations from Pacific Patagonia. Even more, considering the observations of N. magellanica in holdfasts of M. pyrifera, rafting could also constitute an important dispersal mechanism, particularly for the colonization of geographically distant areas.
The general genetic pattern of N. magellanica along its distribution in Pacific Patagonia indicates low levels of nucleotide diversity and number of nucleotide differences. On one hand, the lowest values were found at localities that should have been severely ice-impacted during the LGM including channels and fjords such as Costa Channel, Puerto Aguirre, and Concoto Island. On the other hand, higher diversity levels were found at northern localities (Metri and Puerto Montt), more oceanic areas (Serrano Channel and London Island), the Strait of Magellan and Cape Horn. In spite of these slight diversity differences among localities, theory predicts that large population sizes should maintain high levels of genetic variability, because genetic drift is low and the mutational rate is high. General molecular diversity indices estimated in N. magellanica in the Pacific (θ k =4.84; θ S =2.93; θ H =2.41; θπ=1.77) would be sustained by effective sizes (Ne) between 133,750 and 372,000 individuals. These estimations are smaller by far than the expected population sizes in the species, considering the high densities reported [41, 42, 44, 87–89]. Low levels of genetic diversity together with dominant haplotypes widely distributed are consistent with the hypothesis of a recent range expansion [90, 91]and high levels of migration [3,70]. Moreover, significant negative Tajima’s D and Fu’s F S indices, together with a unimodal mismatch distribution in Pacific Patagonia are the result of an excess of low frequency haplotypes, commonly explained by recent demographic processes.
Traditional genetic models of glacial refugia and routes of recolonization include the prediction of low genetic diversity in formerly glaciated areas with a small number of haplotypes dominating disproportionally large areas, and high diversity in glacial refugia [2, 7, 71]. Based on the patterns of genetic structure in N. magellanica, the hypothesis of persistence of the species in multiple glacial refugia along Pacific Patagonia followed by expansion from surviving populations is most unlikely. If these periglacial populations experienced strong bottlenecks during the LGM, they may exhibit low genetic diversity as expected in recolonized areas with no refugia, but should have more endemic diversity than recently recolonized areas . Even in the presence of high levels of gene flow under a multiple in situ refugia hypothesis, it is expected to find more than one haplotype exhibiting high frequency and each of these haplotypes could derive from a putative glacial refuge. However, in the case of N. magellanica, coupled with the low levels of nucleotide diversity we found an absence of genetic differentiation along Pacific Patagonia with just one dominant haplotype (H1). The lack of structure in a large geographical area with a single dominant haplotype, the absence of correlation between geographic and genetic distance and the evidence of recent demographic growth support the hypothesis of a recent expansion in the species, possibly mediated by its indirect development. Such larval-mediated postglacial recolonization processes in the Northern Hemisphere have been frequently recognized in marine organisms [70, 92, 93]. In contrast, few studies have established the importance of the developmental mode in the phylogeographic patterns of marine invertebrates in Patagonia. Absence of genetic structure, as observed in N. magellanica in Pacific Patagonia, have been also recognized in other marine organisms with indirect development, including the mytilid Mytilus edulis, fishes like Eleginops maclovinus and Sebastes oculatus, and in macroalgae including Durvillaea antartica and Macrocystis pyrifera. These studies contrast with the results obtained in the direct developer Acanthina monodon that exhibits marked differentiation between northern and southern Pacific Patagonia localities . The patterns of genetic structure observed in different groups of marine organisms across Patagonia further support the importance of the developmental mode and the prevailing directions of currents and winds.
According to Hein et al. the timing of the LGM extent and the onset of deglaciation occurred broadly synchronously throughout Patagonia. In northern areas of Pacific Patagonia the final ice advance is dated about 17.9 ka [97, 98] and warming began at 17.5 ka . Similarly, around the Strait of Magellan the final ice advance occurred prior to ca. 17 ka  while a major and rapid warming period occurred between 14–10 ka [23, 100, 101]. Pollen statrigraphic studies in the Beagle Channel and Tierra del Fuego suggest that the disappearance of ice in that sector occurred ~ 11.6 ka [96, 101, 102]. Based on our estimations, population expansion in N. magellanica would have occurred ~ 6.3 ka under a sudden growth model and ~ 9.0 ka under the Bayesian Skyline Plot approximation. Estimated dates of population expansion in N. magellanica are consistent with previous analysis in the species  and with thermal records of warmer conditions in Patagonia. Also, paleontological studies on postglacial mollusk faunas in the northern coast of the Beagle Channel suggest that major expansion of taxa occurred after the glaciers receded fully (~ 10 ka). Under the relatively warmer conditions of the middle Holocene (5.0 to 4.0 ka), the fossil record indicates a process of diversification of several mollusk taxa including Nacella[103–105].
Genetic differentiation among Pacific Patagonia, Atlantic Patagonia and the Falkland/Malvinas Islands
In spite that N. magellanica constitute a single panmictic unit along Pacific Patagonia, we detected clear differences among three main areas in Patagonia. Marked levels of genetic structure in the species among these areas may be explained by differences in their respective glaciological histories, rocky shore availability and the prevailing currents and winds among them. First, ice-shelf advances and retreats differentially affected Pacific and Atlantic Patagonia during Pleistocene glacial cycles. Pacific Patagonia was almost completely covered by ice during the LGM and shallow marine habitats should have been severely affected. In contrast Atlantic Patagonia was only affected over the piedmont areas to the east and to the current submarine platform south of Río Gallegos [96, 106, 107]. Sea levels changes during the LGM might have differentially affected Atlantic and Pacific populations. Atlantic and Falkland Island populations might have had even less rocky substrated than today, with exposure of the shelf, and would also have moved north and eastward during the LGM following the shorelines. Geomorphologic evidence in West Falkland/Malvinas suggests that during the Cenozoic the largest glaciers were no more than 2.7 km long . According to this scenario, the Pacific population would have been more severely hampered during LGM, as suggested by significant negative Tajima and Fu’s tests and an L-shaped mismatch distribution. In contrast, populations from Atlantic Patagonia and Falkland/Malvinas Islands did not show deviation from the mutation-drift equilibrium model and also exhibited multimodal mismatch distributions, in agreement with the expectation for more stable populations (Figure 4). The results observed in Puerto Deseado did not match with those recently published by de Aranzamendi et al. that detected signal of recent demographic expansion along Atlantic Patagonia. In this respect, contrasted demographic signals detected here between Pacific and Atlantic populations must be taken cautiously considering the differences in sample size . However, the main signal detected here in the whole Pacific data set is also observed in each locality from the Reloncaví Fjord to Cape Horn and sample size in each one of them is comparable to the ones of Puerto Deseado and the Falkland/Malvinas Islands (Table 1). Second, the Atlantic coast includes less rocky shore ecosystems than Pacific Patagonia where more suitable rocky habitat are available for the species. These differences between Pacific and Atlantic Patagonia in terms of habitat availability, has been also recognized in comparative biodiversity studies in the southern tip of South America . Off southern Chile the Cape Horn Current (CHC) flows southward around Cape Horn while the Malvinas-Falkland Current (M-FC) is a northward running branch of the former that moves about up to 28°S following the continental shelf margin [110–112]. The minor Patagonian Coastal Current (PCC) influences the Atlantic Patagonia coast and moves northward up to 38°S [110–112]. Moreover, oceanographic fronts such as the Atlantic Patagonian cold estuarine front on the eastern margin of the Strait of Magellan has been described as an oceanographic barrier between Pacific Patagonia and the southern Atlantic coast [113–115]. As expected under the general circulation pattern in this region, we found asymmetrical gene flow in the species among Pacific Patagonia, Atlantic Patagonia and the Falkland/Malvinas Islands. Migration rates to Pacific Patagonia from the Atlantic population and from the Falkland/Malvinas Islands were very low. Similarly, migration rates between these islands and Atlantic Patagonia were almost negligible. In this respect, it is proble that Atlantic Patagonia is continuously receiving haplotypes from the Pacific. In fact, haplotypes from Puerto Deseado are shared with Pacific Patagonia localities and especially with those from Tierra del Fuego, the Beagle Channel, and the eastern mouth of the Strait of Magellan Additional file 2. Considering the migration estimation from the Atlantic to the Pacific and the prevailing direction of the oceanic fronts and winds, the hypothesis of postglacial recolonization from Atlantic Patagonia to the Pacific is most unlikely. Similar patterns of genetic discontinuities between Pacific and Atlantic localities in Patagonia have been found in Enteroctopus megalocyathus and in the scallop Aequipecten tehuelchus. Moreover, de Aranzamendi et al. found significant genetic differentiation between N. magellanica populations from Tierra del Fuego and northern localities such as Golfo San José and Golfo Nuevo.
In the Falkland/Malvinas Islands, even when we included only 13 individuals we detected the highest levels of genetic diversity. Furthermore, this population was characterized by positive Tajima’ D and Fu’s F S indices, a multimodal mismatch distribution and an expanded genealogy (Figure 4). According to the Quaternary genetic model [2, 3, 7, 71], the Falkland/Malvinas Islands could be considered as a glacial refugium, considering the higher level of genetic diversity and the presence of endemic haplotypes (H53, H52 and H55) clearly differentiated from the Patagonian diversity. Moreover, these islands have been previously proposed as refugial areas for several plant species during the LGM [21, 108, 118] and as relevant area for conservation .
However, even if our data support the persistence of N. magellanica in Falkland/Malvinas Islands during the LGM, they do not support a scenario of posterior recolonization from Falkland/Malvinas Islands to Atlantic and/or Pacific Patagonia. Considering our migration rate estimations, most of the gene flow in N. magellanica is derived from Pacific Patagonia to the other areas. The particular case of the Falkland/Malvinas Islands seems to represent a sink area where private surviving haplotypes are mixed together with recently arrived ones from Pacific Patagonia.
Historical factors and life-history traits such as its indirect development play a main role in the connectivity of N. magellanica. In concert with the high dispersal potential of the species, we detected a rapid postglacial recolonization process in a very complex landscape likely related the deglaciation process along Pacific Patagonia. In contrast to the model of Pleistocene biogeography, where higher levels of genetic diversity are expected at lower latitudes, in N. magellanica we did not detect a clear relationship between latitude and genetic diversity. The absence of evidence of a progressive southward recolonization through recurrent founder effects may be the result of the synchronous deglaciation process along Pacific Patagonia  together with high dispersal capacities. In contrast, lower genetic diversity detected in the inland sea, characterized by fjords and channels, could indicate that these areas represent those most recently recolonized by N. magellanica. In this region, the timing of recolonization would therefore have followed a west to east trend, contrasting to the usual north–south model of Pleistocene biogeography in South America. At the same time, this study gives further evidence for the role of the major current systems among different areas of Patagonia through the existence of an asymmetrical pattern of gene flow from Pacific Patagonia to Atlantic Patagonia and the Falkland/Malvinas Islands following the CHC, the M-FC and the PCC. According to our results, N. magellanica persisted in the Falkland/Malvinas Islands during the Quaternary glacial cycles and therefore represents a relict population. However, the pattern of genetic diversity strongly suggests that this population did not participate to the postglacial recolonization of southern South America. Considering oceanic and atmospheric circulation in the province and the pattern of gene flow we found among Patagonian sectors, the Falkland/Malvinas Islands seem to represent a sink area where recently arrived and endemic haplotypes coexist. The main pattern of genetic diversity and structure in N. magellanica appears to be the result of the combination of the impact of the last glacial period in Pacific Patagonia and the prevailing oceanographic circulation, together with life-history traits like its indirect development and its narrow bathymetric range. These historical and contemporary processes may constitute important factors in shaping the modern biogeography of most shallow marine benthic invertebrates inhabiting the Patagonian Province.
Future research in N. magellanica will include a broader sampling effort along the Atlantic coast and the use of recently developed fast-evolving markers  in order to corroborate the observed pattern of genetic structure. Finally, more studies on other species of Nacella as well as other marine benthic taxa are required in order to provide a better understanding of the historical and recent processes governing the patterns of genetic structure and connectivity in southern South America. This information will provide an empirical framework in order to generalize the postglacial biogeographic model proposed here for the limpet N. magellanica.
Sample collection, DNA extraction, PCR amplification and sequencing
Individuals were collected between 2007–2011 from the intertidal zone in 14 localities (Figure 5). Along Pacific Patagonia we included two localities in the northern limit of the species distribution in the Reloncaví Fjord (41.5°S), four localities from the Chonos Archipelago (44° - 46°S), two localities in the Strait of Magellan (52° - 53°S) and three localities from the Cape Horn Biosphere Reserve (54° - 55°S; Figure 5). We also included in the analyses two localities in Atlantic Patagonia, Puerto Deseado (47°45’ S; 65°52’ W) and Possession Bay (52°17’ S; 68°57’ W) in the eastern mouth of the Strait of Magellan. Finally, we included in the analyses 13 individuals from the Falkland/Malvinas Islands (51°41’ S; 57°50 W; Figure 5). Specimens were identified based on shell morphology, sculpture, height, and coloration  and with the help of diagnostic external characteristic of the species such as coloration of the foot muscle and the mantle tentacle . Individuals were fixed in ethanol (95%) and whole DNA was extracted from the mantle using the salting-out method described by . PCR amplifications of a partial fragment of the mtDNA gene Cytochrome c Oxidase Subunit I (COI) were performed using the universal primers described by Folmer et al. and PCR conditions were done following . Double stranded PCR products were purified using the QIAquick Gel Extraction Kit (QUIAGEN), and sequenced in both directions using an Automatic Sequencer ABI3730 x 1 at Macrogen, Inc. (Seoul, Korea).
DNA chromatograms were manually edited using Proseq v. 2.91  and aligned with ClustalW . COI sequences were translated to amino acids to check for sequencing errors and/or the presence of pseudogenes with MEGA 5.0 . We performed a DNA saturation analysis following Roe & Sperling  to evaluate the levels of saturation changes along the N. magellanica COI data set. New COI sequences have been submitted to GenBank database (Accession Numbers: JX262742 – JX262797).
We estimated the levels of polymorphism in N. magellanica using standard diversity indices, haplotype number (k), the number of segregating sites (S) and haplotypic diversity (H) for each locality and for the whole COI data set using DnaSP v.5.00.07 . We also estimated average pairwise sequence differences (Π) and nucleotide diversity (π). Population parameters (Tajima’s D and Fu’s F S ) were calculated for all populations using DnaSP and Arlequin v.3.11 .
Genetic differentiation was determined in two ways following [129, 130] using mean pair-wise differences (NST) and through their haplotype frequencies (GST) in Arlequin. We performed permutation tests (25,000 random iterations) of both coefficients to confirm statistical differences among the analyzed localities. Moreover, both parameters were compared with Permut (http://www.pierroton.inra.fr/genetics/labo/Software/) using an analytical test. We tested whether NST >> GST by comparison of the NST values measured directly with those obtained after 1000 random permutation of haplotype identities . Using SAMOVA v.1 (Spatial Analysis of Molecular Variance)  we defined the number and composition of geographically homogeneous, maximally differentiated groups of localities. This method aims to maximize the proportion of total genetic variance due to differences among groups minimizing the variance portion among population within groups. Once these groups were defined, we estimated the levels of migration among them using a Markov Monte Carlo coalescent genealogy sampler implemented in LAMARC v.2.1.8 (Likelihood Analysis with Metropolis Algorithm using Random Coalescence) . This approximation allows to estimate migration levels among the recognized groups of N. magellanica and at the same time to test whether the migration was symmetric or asymmetric among them. We examined the significance of the correlation between genetic divergence measured as Slatkin’s linearized FST [PhiST/(1 - PhiST)] and geographical distance between localities using a Mantel test implemented in Arlequin; associated probabilities were estimated with 25,000 permutations.
We reconstructed genealogical relationships for N. magellanica using median-joining haplotype networks in Network v.4.6 (http://www.fluxus-engineering.com) . To estimate past population dynamics in the species within Pacific Patagonia we applied two methods. First, we used the sudden population growth model  which rests on the assumption that population growth and decline events leave characteristic signatures in the distribution of nucleotide site differences between pairs of individuals. We constructed the distribution of pairwise differences (mismatch distribution) in N. magellanica to determine whether N. magellanica has undergone sudden population growth. We compared the distribution of pairwise differences in N. magellanica with expectations of a sudden expansion model. Three main parameters were estimated: i) the date of growth/decline (τ=2μt) measured in units of 1/2 μ generations where t=time in years and μ=mutation rate per sequence per generation, initial population size (Θ i ) before the population growth/decline and a final theta (Θ f ) after population growth/decline. These demographic expansion parameters were determined using a nonlinear least squares approach implemented in Arlequin . The goodness of fit between the observed and expected mismatch distributions was tested using a parametric bootstrap approach that uses the sum of squared deviations as a statistic test implemented in Arlequin. Second, we used a Bayesian skyline plot method implemented in BEAST v. 1.6 , which detects demographic signatures from nucleotide sequences that are not readily described by simple demographic models [137, 138]. We analyzed the data set under an uncorrelated lognormal relaxed molecular clock model using an evolutionary rate of 1.1% per million years estimated for COI in nacellids , using the GTR + G + I model previously estimated with MrModeltest v.2.3 (http://www.abc.se/~nylander/) and a piecewise constant Bayesian skyline model with 10 groups. Before choosing this model we performed Bayesian Skyline Plot analyses using N. magellanica’s COI data set with three different models: an uncorrelated lognormal relaxed clock, an uncorrelated exponential relaxed clock, and a strict molecular clock. Estimated bayes factors among these models strongly supported the molecular clock hypothesis. We ran the analyses for 350 × 106 generations, making sure that the effective sampling sizes for each statistic were at least 1500. Convergence was examined in Tracer v. 1.5 (http://beast.bio.ed.ac.uk/Tracer) .
Clark J, Fastie C, Hurtt G, Jackson S, Johnson C, King G, Lewis M, Lynch J, Pacala S, Prentice C, Schupp E, Webb T, Wyckoff P: Reid's paradox of rapid plant migration: Dispersal theory and interpretation of paleoecological records. Bioscience. 1998, 48: 13-24. 10.2307/1313224.
Hewitt G: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405: 907-913. 10.1038/35016000.
Marko P, Hoffman J, Emme S, McGovern T, Keever C, Cox N: The “Expansion-Contraction” model of Pleistocene biogeography: rocky shores suffer a sea change?. Mol Ecol. 2010, 19: 146-169. 10.1111/j.1365-294X.2009.04417.x.
Bennett K, Tzedakis P, Willis K: Quaternary refugia of North European trees. J Biogeogr. 1991, 18: 103-115. 10.2307/2845248.
Webb T, Bartlein P: Global changes during the last 3 million years: climatic controls and biotic responses. Annu Rev Ecol Evol Syst. 1992, 23: 141-173.
Hewitt G: Genetic consequences of climatic oscillations in the Quaternary. Philosophical Transactions: Biological Sciences. 2004, 359: 183-195. 10.1098/rstb.2003.1388.
Provan J, Bennett K: Phylogeographic insights into cryptic glacial refugia. Trends in Ecology & Evolution. 2008, 23: 564-571. 10.1016/j.tree.2008.06.010.
Williams J, Summers R, Webb T: Applying plant functional types to construct biome maps from eastern North American pollen data. Quat Sci Rev. 1998, 17: 607-627. 10.1016/S0277-3791(98)00014-6.
Hewitt G: Some genetic consequences of ice ages, and their role, in divergence and speciation. Biol J Linn Soc. 1996, 58: 247-276.
Lessa E, Cook J, Patton J: Genetic footprints of demographic expansion in North America, but not Amazonia, during the late Quaternary. Proceedings of the National Academy of Science, USA. 2003, 100: 10331-10334. 10.1073/pnas.1730921100.
Lessa E, D'Elía G, Pardiñas U: Genetic footprints of late Quaternary climate change in the diversity of Patagonian-Fueguian rodents. Mol Ecol. 2010, 19: 3031-3037. 10.1111/j.1365-294X.2010.04734.x.
Beheregaray L: Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere. Mol Ecol. 2008, 17: 3754-3774.
Ruzzante D, Walde S, Cussac V, Dalebout M, Seibert J, Ortubay S, Habit E: Phylogeography of the Percichthyidae (Pisces) in Patagonia: roles of orogeny, glaciation, and volcanism. Mol Ecol. 2006, 15: 2949-2968. 10.1111/j.1365-294X.2006.03010.x.
Ruzzante D, Walde S, Gosse J, Cussac V, Habit E, Zemlak T, Adams E: Climate control on ancestral population dynamics: insight from Patagonian fish phylogeography. Mol Ecol. 2008, 17: 2234-2244. 10.1111/j.1365-294X.2008.03738.x.
Fraser C, Thiel M, Spencer H, Waters J: Contemporary habitat discontinuity and historic glacial ice drive genetic divergence in Chilean kelp. BMC Evol Biol. 2010, 10: 203-10.1186/1471-2148-10-203.
Zemlak T, Habit E, Walde S, Carrea C, Ruzzante D: Surviving historical Patagonian landscapes and climate: molecular insights from Galaxias maculatus. BMC Evol Biol. 2010, 10: 67-10.1186/1471-2148-10-67.
González-Wevar CA, Nakano T, Cañete JI, Poulin E: Concerted genetic, morphological and ecological diversification in Nacella limpets in the Magellanic Province. Mol Ecol. 2011, 20: 1936-1951. 10.1111/j.1365-294X.2011.05065.x.
Sérsic A, Cosacov A, Cocucci A, Johnson L, Pozner R, Avila L, Sites J, Morando M: Emerging phylogeographical patterns of plants and terrestrial vertebrates from Patagonia. Biol J Linn Soc. 2011, 103: 475-494. 10.1111/j.1095-8312.2011.01656.x.
Ceballos SG, Lessa EP, Victorio MF, Fernández DA: Phylogeography of the sub-Antarctic notothenioid fish Eleginops maclovinus: evidence of population expansion. Mar Biol. 2012, 159: 499-505. 10.1007/s00227-011-1830-4.
Palma R, Boric-Bargetto D, Torres-Pérez F, Hernández C, Yates T: Glaciation effects on the phylogeographic structure of Oligoryzomys longicaudatus (Rodentia: Sigmodontinae) in the Southern Andes. PLoS One. 2012, 7: e32206-10.1371/journal.pone.0032206.
Clapperton C: Quaternary Geology and Geomorphology of South America. 1993, Elsevier Amsterdam, NL: Elsevier Science, 796-
McCulloch R, Bentley M, Purves R, Hulton N, Sugden D, Clapperton C: Climatic inferences from glacial and paleoecological evidence at the last glacial termination, southern South America. J Quat Sci. 2000, 15: 409-417. 10.1002/1099-1417(200005)15:4<409::AID-JQS539>3.0.CO;2-#.
Hulton N, Purves P, McCulloch R, Sugden D, Bentley M: The last glacial maximum and deglaciation in southern South America. Quat Sci Rev. 2002, 21: 233-241. 10.1016/S0277-3791(01)00103-2.
RABASSA J: Late Cenozoic glaciation in Patagonia and Tierra del Fuego. Developments in Quaternary Sciences. 2008, 11: 151-205.
Kim I, Phillips C, Monjeau J, Birney E, Noack K, Pumo D, Sikes R, Dole J: Habitat islands, genetic diversity, and gene flow in a Patagonian rodent. Mol Ecol. 1998, 7: 667-678. 10.1046/j.1365-294x.1998.00369.x.
Smith M, Kelt D, Patton J: Testing models of diversification in mice in the Abrothrix olivaceus/xanthorhinus complex in Chile and Argentina. Mol Ecol. 2001, 10: 397-405. 10.1046/j.1365-294x.2001.01183.x.
Cussac V, Ortubay S, Iglesias G, Milano D, Lattuca M, Barriga J, Battini M, Gross M: The distribution of South American galaxiid fishes: the role of biological traits and post-glacial history. J Biogeogr. 2003, 31: 103-121.
Zemlak T, Habit E, Walde S, Battini M, Adams E, Ruzzante D: 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: 5049-5061. 10.1111/j.1365-294X.2008.03987.x.
Victoriano P, Ortiz J, Benavides E, Adams B, Sites J: Comparative phylogeography of codistributed species of Chilean Liolaemus (Squamata: Tropiduridae) from the central-southern Andean range. Mol Ecol. 2008, 17: 2397-2416. 10.1111/j.1365-294X.2008.03741.x.
Breitman M, Avila L, Sites J, Morando M: Lizards from the end of the world: Phylogenetic relationships of the Liolaemus lineomaculatus section (Squamata: Iguania: Liolaemini). Molecular Phylogenetics and Evolution. 2011, 59: 364-376. 10.1016/j.ympev.2011.02.008.
Vidal M, Moreno P, Poulin E: Genetic diversity and insular colonization of Liolaemus pictus (Squamata, Liolaeminae) in north-western Patagonia. Austral Ecology. 2011, 37: 67-77.
Nuñez J, Wood N, Rabanal F, Fontanella F, Sites J: Amphibian phylogeography in the Antipodes: Refugia and postglacial colonization explain mitochondrial haplotype distribution in the Patagonian frog Eupsophus calcaratus (Cycloramphidae). Molecular Phylogenetics and Evolution. 2011, 58: 343-352. 10.1016/j.ympev.2010.11.026.
Himes C, Gallardo M, Kenagy G: Historical biogeography and post-glacial recolonization of South American temperate rain forest by the relictual marsupial Dromiciops gliroides. J Biogeogr. 2008, 35: 1415-1424. 10.1111/j.1365-2699.2008.01895.x.
Marín J, Spotorno A, González B, Bonacic C, Wheeler J, Casey C, Bruford M, Palma R, Poulin E: Mitochondrial DNA variation and systematics of the guanaco (Lama guanicoe, Artiodactyla: Camelidae). J Mammal. 2008, 89: 1-13. 10.1644/07-MAMM-F-268.1.
Ufj P, Teta P, D'Elía G, Lessa E: The evolutionary history of sigmodontine rodents in Patagonia and Tierra del Fuego. Biol J Linn Soc. 2011, 103: 495-513. 10.1111/j.1095-8312.2011.01679.x.
Jakob S, Martinez-Meyer E, Blattner F: Phylogeographic analyses and paleodistribution modeling indicate Pleistocene In Situ survival of Hordeum species (Poaceae) in Southern Patagonia without genetic or spatial restriction. Mol Biol Evol. 2009, 26: 907-923. 10.1093/molbev/msp012.
Cárdenas L, Castilla J, 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: 969-981. 10.1111/j.1365-2699.2008.02056.x.
Macaya E, Zuccarello G: Genetic structure of the giant kelp Macrocystis pyrifera along the southeastern Pacific. Mar Ecol Prog Ser. 2010, 420: 103-112.
Sánchez R, Sepúlveda R, Brante A, Cárdenas L: Spatial pattern of genetic and morphological diversity in the direct developer Acanthina monodon (Gastropoda: Mollusca). Mar Ecol Prog Ser. 2011, 434: 121-131.
Powell A: The Patellid limpets of the World (Patellidae). Indo-Pacific Mollusca. Abbot. 1973, Greenville: The Department of Mollusks, 75-206. 3
Gúzman L: Patrón de distribución espacial y densidad de Nacella magellanica (Gmelin, 1791) en el intermareal del Estrecho de Magallanes (Mollusca, Gastropoda). Anales del Instituto de la Patagonia. 1978, 9: 205-219.
Ríos C, Mutschke E: Community structure of intertidal boulder-cobble fields in the Straits of Magellan, Chile*. Sci Mar. 1999, 63: 193-201.
Valdovinos C, Rüth M: Nacellidae limpets of southern South America: taxonomy and distribution. Rev Chil Hist Nat. 2005, 78: 497-517.
Bazterrica M, Silliman B, Hidalgo F, Crain C, Bertness M: Limpet grazing on a physically stressful Patagonian rocky shore. J Exp Mar Biol Ecol. 2007, 353: 22-34. 10.1016/j.jembe.2007.08.018.
González-Wevar C, Nakano T, Cañete J, Poulin E: Molecular phylogeny and historical biogeography of Nacella (Patellogastropoda: Nacellidae) in the Southern Ocean. Molecular Phylogenetics and Evolution. 2010, 56: 115-124. 10.1016/j.ympev.2010.02.001.
De Aranzamendi M, Gardenal C, Martin J, Bastida R: Limpets of the genus Nacella (Patellogastropoda) from the Southwestern Atlantic: species identification based on molecular data. Journal of Molluscan Studies. 2009, 75: 241-251. 10.1093/mollus/eyp025.
Morriconi E: Reproductive biology of the limpet Nacella (P.) deaurata (Gmelin, 1791) in Bahía Lapataia (Beagle Channel). Sci Mar. 1999, 63: 417-426.
Adami M, Gordillo S: Structure and dynamics of the biota associated with Macrocystis pyrifera (Phaeophyta) from the Beagle Channel, Tierra del Fuego. Sci Mar. 1999, 63: 183-191.
De Aranzamendi M, Bastidas R, Gardenal C: Different evolutionary histories in two sympatric limpets of the genus Nacella (Patellogastropoda) in the South-western Atlantic coast. Mar Biol. 2011, 158: 2405-2418. 10.1007/s00227-011-1742-3.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
González-Wevar CA, David B, Poulin E: Phylogeography and demographic inference in Nacella (Patinigera) concinna (Strebel, 1908) in the western Antarctic Peninsula. Deep-Sea Res II Top Stud Oceanogr. 2011, 58: 220-229. 10.1016/j.dsr2.2010.05.026.
Hammer O, Harper D, Ryan P: PAST: Paleontological statistics software package for education and data analysis. Palaeontol Electron. 2001, 4: 1-9.
Posada D, Crandall K: Intraspecific gene genealogies: trees grafting into networks. Trends in Ecology & Evolution. 2001, 16: 37-45. 10.1016/S0169-5347(00)02026-7.
Ho S, Phillips M, Cooper A, Drummond A: Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol Biol Evol. 2005, 22: 1561-1568. 10.1093/molbev/msi145.
Ho S, Larson G: Molecular clocks: when times are a-changin'. Trends in Genetics. 2006, 22: 79-83. 10.1016/j.tig.2005.11.006.
Ho S, Shapiro B, Phillips M, Cooper A, Drummond A: Evidence for time dependency of molecular rates. Syst Biol. 2007, 56: 515-522. 10.1080/10635150701435401.
Ho S, Lanfear R, Bromham L, Phillips M, Soubrier J, Rodrigo A, Cooper A: Time-dependent rates of molecular evolution. Mol Ecol. 2011, 20: 3087-3101. 10.1111/j.1365-294X.2011.05178.x.
Burridge C, Craw D, Fletcher D, Waters J: Geological dates and molecular rates: fish DNA sheds light on time dependency. Mol Biol Evol. 2008, 25: 624-633. 10.1093/molbev/msm271.
McCarty J: Ecological consequences of recent climate change. Conserv Biol. 2001, 15: 320-331. 10.1046/j.1523-1739.2001.015002320.x.
Davis M, Shaw R, Etterson J: Evolutionary responses to climate change. Ecology. 2005, 86: 1704-1714. 10.1890/03-0788.
Parmesan C: Ecological and evolutionary responses to recent climate change. The Annual Reviews of Ecology, Evolution, and Systematics. 2006, 37: 1-35. 10.1146/annurev.ecolsys.37.091305.110035.
Marshall J, Blair J, Peters D, Okin G, Rango A, Williams M: Predicting and understanding ecosystem responses to climate change at continental scales. Front Ecol Environ. 2008, 6: 273-280. 10.1890/070165.
Guisan A, Thuiller W: Predicting species distribution: offering more than simple habitat models. Ecol Lett. 2005, 8: 993-1009. 10.1111/j.1461-0248.2005.00792.x.
Lesbarrères D: Post-glacial phylogeography: New insight into an old story: the post-glacial recolonization of European biota. Heredity. 2008, 102: 213-213.
Huntley B: How plants respond to climate change: migration rates, individualism and the consequences for the plant communities. Ann Bot. 1991, 67: 15-22.
Davis M, Shaw R: Range shifts and adaptive responses to Quaternary climate change. Science. 2001, 292: 673-679. 10.1126/science.292.5517.673.
Hughes L: Biological consequences of global warming: is the signal already apparent?. Trends in Ecology & Evolution. 2000, 15: 56-61. 10.1016/S0169-5347(99)01764-4.
Lima F, Ribeiro P, Queiroz N, Hawkins S, Santos A: Do distributional shifts of northern and southern species of algae match the warming pattern?. Glob Chang Biol. 2007, 13: 2592-2604. 10.1111/j.1365-2486.2007.01451.x.
Pearson R, Dawson T: Predicting the impacts of climate change on the distribution of species: are bioclimate envelope models useful?. Glob Ecol Biogeogr. 2003, 12: 361-371. 10.1046/j.1466-822X.2003.00042.x.
Marko P: 'What‘s larvae got to do with it?’ Disparate patterns of post-glacial population structure in two benthic marine gastropods with identical dispersal potential. Mol Ecol. 2004, 13: 597-611. 10.1046/j.1365-294X.2004.02096.x.
Maggs C, Castihlo R, Foltz D, Henzler C, Jolly M, Kelly J, Olsen J, Perez K, Stam W, Vainöla R, Viard F, Wares J: Evaluating signatures of glacial refugia for North Atlantic benthic marine taxa. Ecology. 2008, 89: S108-S122. 10.1890/08-0257.1.
Nakano T, Ozawa T: Systematic revision of Patelloida pygmaea (Dunker, 1860) (Gastropoda: Lottiidae), with a description of a new species. Journal of Molluscan Studies. 2005, 71: 357-370. 10.1093/mollus/eyi039.
Bird C, Holland B, Bowen B, Toonen R: Contrasting phylogeography in three endemic Hawaiian limpets (Cellana spp.) with similar life histories. Mol Ecol. 2007, 16: 3173-3186. 10.1111/j.1365-294X.2007.03385.x.
Sá-Pinto A, Branco M, Sayanda D, Alexandrino P: Patterns of colonization, evolution and gene flow in species of the genus Patella in the Macaronesian Islands. Mol Ecol. 2007, 17: 519-532. 10.1111/j.1365-294X.2007.03563.x.
Clarke A, Crame J: Evolutionary dynamics at high latitudes: speciation and extinction in polar marine faunas. Philosophical Transactions of the Royal Society B: Biological Sciences. 2010, 365: 3655-3666. 10.1098/rstb.2010.0270.
Camus P: Biogeografía marina de Chile continental. Rev Chil Hist Nat. 2001, 74: 587-617.
Silva N, Palma S: The CIMAR Program in the austral Chilean channels and fjords. Progress in the. Edited by: Silva N, Palma S. 2012, Comité Oceanográfico Nacional, 11-15.
Miloslavich P, Klein E, Díaz J, Hernández C, Bigatti G, Campos L, Artigas F, Castillo J, Penchaszadeh P, Neill P, Carranza A, Retana M, Díaz de Astarloa J, Lewis M, Yorio P, Piriz M, Rodríguez D, Yoneshigue-Valentin Y, Gamboa L, Martín A: Marine Biodiversity in the Atlantic and Pacific Coasts of South America: Knowledge and Gaps. PLoS One. 2011, 6: e14631-10.1371/journal.pone.0014631.
Pantoja S, Iriarte J, Daneri G: Oceanography of the Chilean Patagonia. Cont Shelf Res. 2011, 31: 149-153. 10.1016/j.csr.2010.10.013.
Rozzi R, Armesto J, Gutiérrez F, Linkens G, Anderson C, Poole A, Moses K, Hargrove E, Mansilla A, Kennedy J, Wilson M, Jax K, Jones C, Callicott J, Arroyo M: Integrating ecology and environmental ethics: Earth stewardship in the southern end of the Americas. Bioscience. 2012, 62: 226-236.
Ayre D, Minchinton T, PERRIN C: Does life history predict past and current connectivity for rocky intertidal invertebrates across a marine biogeographic barrier?. Mol Ecol. 2009, 18: 1887-1903. 10.1111/j.1365-294X.2009.04127.x.
Paulay G, Meyer C: Dispersal and divergence across the greatest ocean region: Do larvae matter?. Integr Comp Biol. 2006, 46: 269-281. 10.1093/icb/icj027.
Bowden D, Clarke A, Peck L, Barnes D: Antarctic sessile marine benthos: colonisation and growth on artificial substrata over three years. Mar Ecol Prog Ser. 2006, 316: 1-16.
Bowden D, Clarke A, Peck L: Seasonal variation in the diversity and abundance of pelagic larvae of Antarctic marine invertebrates. Mar Biol. 2009, 156: 2033-2047. 10.1007/s00227-009-1235-9.
Green B, Fisher R: Temperature influences swimming speed, growth and larval duration in coral reef fish larvae. J Exp Mar Biol Ecol. 2004, 299: 115-132. 10.1016/j.jembe.2003.09.001.
O'Connor M, Bruno J, Gaines S, Halpern B, Lester S, Kinlan B, Weiss J: Temperature control larval dispersal and the implications for marine ecology, evolution, and conservation. Proceedings of the National Academy of Science, USA. 2007, 104: 1266-1271. 10.1073/pnas.0603422104.
Gúzman L, Ríos C: Age and growth of the Subantarctic limpet Nacella (Patinigera) magellanica magellanica (Gmelin,1791) from the Strait of Magellan, Chile. Veliger. 1987, 30: 159-166.
Ríos C, Mutschke E, Morrison E: Biodiversidad bentónica sublitoral en el estrecho de Magallanes, Chile. Rev Biol Mar Oceanogr. 2003, 38: 1-12.
Bertness M, Crain C, Silliman B, Bazterrica M, Reyna M, Hidalgo F, Farina J: The community structure of western Atlantic Patagonian rocky shores. Ecol Monogr. 2006, 76: 439-460. 10.1890/0012-9615(2006)076[0439:TCSOWA]2.0.CO;2.
Slatkin M, Hudson R: Pairwise comparisons of mitochondrial DNA sequences in stable and exponentially growing populations. Genetics. 1991, 129: 555-562.
Rogers A, Harpending H: Population growth makes waves. Mol Biol Evol. 1992, 9: 552-569.
Hickerson M, Ross J: Post-glacial population history and genetic structure of the northern clingfish (Gobbiesox meandricus), revealed from mtDNA analysis. Mar Biol. 2001, 138: 1-13. 10.1007/s002270000437.
Flight P, O'Brien M, Schmidt P, Rand D: Genetic Structure and the North American Postglacial Expansion of the Barnacle, Semibalanus balanoides. J Hered. 2012, 103: 153-165. 10.1093/jhered/esr083.
Toro J, Ojeda J, Vergara A: The genetic structure of Mytilus chilensis (Hupe 1854) populations along the Chilean coast based on RAPDs analysis. Aquac Res. 2004, 35: 1466-1471. 10.1111/j.1365-2109.2004.01172.x.
Nuñez J, González M, Pérez-Losada M: Testing species boundaries between Atlantic and Pacific lineages of the Patagonian rockfish Sebastes oculatus (Teleostei: Scorpaenidae) through mitochondrial DNA sequences. Rev Biol Mar Oceanogr. 2010, 45: 565-573.
Hein A, Hulton N, Dunai T, Sugden D, Kaplan M, Xu S: The chronology of the Last Glacial Maximum and deglacial events in central Argentine Patagonia. Quat Sci Rev. 2010, 29: 1212-1227. 10.1016/j.quascirev.2010.01.020.
Heusser C, Heusser L, Lowell T: Paleoecology of the Southern chilean lake district-Isla Grande de Chiloe during the Middle-Late Llanquihue glaciation and deglaciation. Geografiska Annaler. Series A, Physical Geography. 1999, 81: 231-284. 10.1111/j.0435-3676.1999.00058.x.
Moreno P, Lowell T, Jacobson G, Denton G: Abrupt vegetation and climate changes during the last termination in the Chilean lake district: a case study from Canal de la Puntilla (41 S). Geogr Ann. 1999, 81: 285-311. 10.1111/j.0435-3676.1999.00059.x.
Moreno P, León A: Abrupt vegetation changes during the last glacial to Holocene transition in mid-latitude South America. J Quat Sci. 2003, 18: 787-800. 10.1002/jqs.801.
McCulloch R, Fogwill C, Sugden D, Bentley M, Kubik P: Chronology of the last glaciation in Central Strait of Magellan and Bahía Inútil, Southernmost South America. Geogr Ann. 2005, 87: 1-24.
Heusser C: Deglacial paleoclimate of the Americansector of the Southern Ocean: Late Glacial–Holocene records from the latitude of Canal Beagle (55°S), Argentine Tierra del Fuego. Palaeogeogr Palaeoclimatol Palaeoecol. 1998, 141: 277-301. 10.1016/S0031-0182(98)00053-4.
De Pol-Holz R, Ulloa O, Dezileau L, Kaiser J, Lamy F, Hebbeln D: Melting of the Patagonian Ice Sheet and deglacial perturbations of the nitrogen cycle in the eastern South Pacific. Geophys Res Lett. 2006, 33: L04704-
Gordillo S, Coronato A, Rabassa J: Quaternary molluscan faunas from the islands of Tierra del Fuego after the Last Glacial Maximum. Sci Mar. 2005, 69: 337-348.
Gordillo S, Rabassa J, Coronato A: Paleoecology and paleobiogeographic pattern of mid-Holocene mollusks from the Beagle Channel (southern Tierra del Fuego, Argentina). Revista Geológica de Chile. 2008, 35: 321-333.
Cárdenas J, Gordillo S: Paleoenvironmental interpretation of late Quaternary molluscan assemblages from southern South America: A taphonomic comparison between the Strait of Magellan and the Beagle Channel. Andean Geology. 2009, 36: 81-93.
Rabassa J, Coronato A, Salemme M: Chronology of the Late Cenozoic Patagonian glaciations and their correlation with biostratigraphic units of the Pampean region (Argentina). Journal of South American Earth Science. 2005, 20: 81-103. 10.1016/j.jsames.2005.07.004.
Rabassa J, Coronato A, Martínez O: Late Cenozoic glaciations in Patagonia and Tierra del Fuego: an updated review. Biol J Linn Soc. 2011, 103: 316-335. 10.1111/j.1095-8312.2011.01681.x.
Clapperton C, Sugden D: The Maximum extent of glaciers in part of west Falkland. J Glaciol. 1976, 17: 73-77.
Ramos-Onsins SE, Rozas J: Statistical properties of new neutrality tests against population growth. Mol Biol Evol. 2002, 19: 2092-2100. 10.1093/oxfordjournals.molbev.a004034.
Bastida R, Roux A, Martínez D: Benthic communities of the Argentine continental shelf. Oceanol Acta. 1992, 15: 687-698.
Aguirre M, Hlebzsebitch J, Delatorre F: Late Cenozoic Invertebrate Paleontology, with emphasis on Mollusks. Development in Quaternary Science 11, The Late Cenozoic of Patagonia and Tierra del Fuego. Edited by: Rabassa J. 2008, Elsevier, 285-326. 1
Aguirre M, Richiano S, Sirch Y: Palaeoenvironments and palaeoclimates of the Quaternary molluscan faunas from the coastal area of Bahía Vera-Camarones (Chubut, Patagonia). Palaeogeogr Palaeoclimatol Palaeoecol. 2006, 229: 251-286. 10.1016/j.palaeo.2005.06.025.
Hinojosa I, Pizarro M, Ramos M, Thiel M: Spatial and temporal distribution of floating kelp in the channels and fjords of southern Chile. Estuarine Coastal Shelf Sci. 2010, 87: 367-377. 10.1016/j.ecss.2009.12.010.
Acha E, Mianzan H, Guerrero R, Favero M, Bava J: Marine fronts at the continental shelves of austral South America. J Mar Syst. 2004, 44: 83-105. 10.1016/j.jmarsys.2003.09.005.
Alemany D, Acha E, Iribarne O: The relationship between marine fronts and fish diversity in the Patagonian Shelf Large Marine Ecosystem. J Biogeogr. 2009, 36: 2111-2124. 10.1111/j.1365-2699.2009.02148.x.
Pardo-Gandarillas M: Aproximaciones filogeográficas para evaluar patrones de conectividad y de historia demográfica en dos especies de pulpos, habitantes de distintas provincias biogeográficas de Sudamérica. Ph.D Thesis. 2012, Universidad de Chile: Ecology Department
Real L, Julio N, Gardenal N, Ciocco N: Genetic variability of Tehuelche scallop, Aequipecten tehuelchus, populations from the Patagonian coasts (Argentina). J Mar Biol Assoc UK. 2004, 84: 235-238. 10.1017/S0025315404009099h.
Morrone JI: Island evolutionary biogeography: analysis of the weevils (Coleoptera: Curculinodae) of the Falkland Islands (Islas Malvinas). J Biogeogr. 2011, 38: 2078-2090. 10.1111/j.1365-2699.2011.02553.x.
Papadopoulou A, Jones AG, Hammond PM, Vogler AP: DNA taxonomy and phylogeography of beetles of the Falkland Islands (Islas Malvinas). Molecular Phylogenetics and Evolution. 2009, 53: 935-947. 10.1016/j.ympev.2009.08.027.
Pöhlmann K, Held C: Isolation and characterization of eight polymorphic microsatellite markers from South American limpets of the Nacella species complex. Conservation Genetic Resources. 2011, 3: 673-676. 10.1007/s12686-011-9430-7.
Aljanabi S, Martinez I: Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Res. 1997, 25: 4692-4693. 10.1093/nar/25.22.4692.
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994, 3: 294-299.
Filatov D: PROSEQ: A software for preparation and evolutionary analysis of DNA sequence data sets. Molecular Ecology Notes. 2002, 2: 621-624. 10.1046/j.1471-8286.2002.00313.x.
Thomson J, Higgins D, Gibson T: CLUSTALW: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Kumar S, Nei M, Dudley J, Tamura K: MEGA: A biologist-centric software for evolutionary analysis of DNA and protein sequences. Briefings in Bioinformatics. 2008, 9: 299-306. 10.1093/bib/bbn017.
Roe A, Sperling F: Patterns of evolution of mitochondrial cytochrome c oxidase I and II DNA and implications for DNA barcoding. Molecular Phylogenetics and Evolution. 2007, 44: 325-345. 10.1016/j.ympev.2006.12.005.
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.
Excoffier L, Laval G, Schneider S: Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005, 1: 37-50.
Pons O, Chaouche K: Estimation, variance and optimal sampling of gene diversity II. Diploid locus. Theoretical and Applied Genetics. 1995, 91: 122-130.
Pons O, Petit R: Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics. 1996, 144: 1237-1245.
Burban C, Petit R, Carcreef E, Jactel H: Rangewide variation of the maritime pine bast scale Matsucoccus feytaudi Duc. (Homoptera: Matsucoccidae) in relation to the genetic structure of its host. Mol Ecol. 1999, 8: 1593-1602. 10.1046/j.1365-294x.1999.00739.x.
Dupanloup I, Schneider S, Excoffier L: A simulated annealing approach to define the genetic structure of population. Mol Ecol. 2002, 11: 2571-2581. 10.1046/j.1365-294X.2002.01650.x.
Kuhner M: LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. 2006, 22: 768-770. 10.1093/bioinformatics/btk051.
Bandelt H, Forster P, Röhl A: Median-Joining Networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48. 10.1093/oxfordjournals.molbev.a026036.
Schneider S, Excoffier L: Estimation of past demographic parameters from a distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics. 1999, 152: 1879-1889.
Drummond A, Rambaut A: BEAST: Bayesian Evolutionary Analysis by Sampling Trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.
Ho S, Shapiro B: Skyline-plot methods for estimating demographic history from nucleotide sequences. Mol Ecol Resour. 2011, 11: 423-434. 10.1111/j.1755-0998.2011.02988.x.
Drummond A, Suchard M, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.
We thank especially the following persons and institutions for help with field work, data analyses and for contributing with specimens to this study: Ceridwen Fraser, César Cárdenas, Luis Vargas Chacoff, Alejandro Riedemann, Angie Díaz, Leyla Cárdenas, Cristián Ibañez, María Cecilia Pardo-Gandarillas, Paola Grendi, Francisco Villarroel; National Museum of Nature and Science, Tokyo, National Museum of Natural History Chile (MNHN), Martin Gusinde Museum, Puerto Williams, Chile.
This study was supported by the following Grants, Institutions and Projects: Thesis support grants INACH B_01_07 and Conicyt 24090009, Conicyt PhD Grant N° D-21060218 and 24090009, IDEAWILD (C.G-W), Grant-in-Aid for JSPS Fellows No. 207024 and the Japan Society for Promotion of Science (T.N.). Projects P05-002 ICM and PFB 023 (Institute of Ecology and Biodiversity, Universidad de Chile) and INACH 02–02, 13–05 and ECOS C06B02 to E.P., A.M. and C.G-W.; FIP Project 2005–44 (A.M. and J.I.C.). Research Programs PR-F2-01CNR-10 and 0273, Universidad de Magallanes to J.I.C. Thanks are also due to international programmes such as CAML, EBA-SCAR and PROSUL–Brazil for encouraging and supporting Antarctic and Subantarctic research in evolution.
The authors declare that they have no competing interests.
CG-W and E.P generated the ideas and design of the study, and all the authors took part in sample collections. CG-W produced the molecular data, and CG-W and E.P were responsible for data analyses. CG-W drafted the original version of the manuscript, and CG-W and E.P produced subsequent versions. All authors participated in a critical review of the manuscript and approved the final version for submission.
Electronic supplementary material
Additional file 1: A) Haplotype network including 357 Nacella magellanica mtDNA COI sequences. Each haplotype is represented by a circle and its size is proportional to its frequency in the whole data set. mv=median vector (theoretical haplotype that has not been collected but should exist). B) Haplotype Network Puerto Deseado; C) Haplotype Network Falkland/Malvinas Islands. (EPS 381 KB)
About this article
- Cape horn current
- Last glacial maximum
- Post-glacial recolonization
- Expansion-contraction model
- Nacella magellanica
- Larval dispersal
- Asymmetric gene flow
- Patagonian Province
- Falkland/Malvinas Islands