Towards a model of postglacial biogeography in shallow marine species along the Patagonian Province: lessons from the limpet Nacella magellanica(Gmelin, 1791)
© González-Wevar et al.; licensee BioMed Central Ltd. 2012
Received: 26 April 2012
Accepted: 25 July 2012
Published: 7 August 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.
KeywordsQuaternary Cape horn current Last glacial maximum Post-glacial recolonization Expansion-contraction model Nacella magellanica Larval dispersal Asymmetric gene flow Patagonian Province Falkland/Malvinas Islands
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.
Number of individuals per locality, their respective diversity indices and neutrality tests results based on mtDNA (COI) sequences
Fu’s F S
G ST (below diagonal) and N ST (above diagonal) pairwise comparisons among the analyzed sites in Patagonia
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).
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 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
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) .
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.
- 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.PubMed
- 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.PubMed
- 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.PubMed Central
- 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.PubMed
- Beheregaray L: Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere. Mol Ecol. 2008, 17: 3754-3774.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMedPubMed Central
- 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.PubMedPubMed Central
- 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.PubMed
- 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.PubMedPubMed Central
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMed
- 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.PubMedPubMed Central
- 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.PubMed
- 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.PubMed
- 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.PubMedPubMed Central
- Rogers A, Harpending H: Population growth makes waves. Mol Biol Evol. 1992, 9: 552-569.PubMed
- 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.PubMedPubMed Central
- 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.PubMed
- 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.PubMed
- 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.PubMedPubMed Central
- 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.PubMed
- 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.PubMedPubMed Central
- 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.PubMed
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.PubMed
- 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.PubMed
- Pons O, Petit R: Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics. 1996, 144: 1237-1245.PubMedPubMed Central
- 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.PubMed
- 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.PubMed
- Kuhner M: LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. 2006, 22: 768-770. 10.1093/bioinformatics/btk051.PubMed
- 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.PubMed
- 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.PubMedPubMed Central
- 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.PubMed
- 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.PubMedPubMed Central
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.