BACKGROUND
climate change caused substantial alterations of the landscape and sea level, influencing patterns of species distribution.
pleistocene glaciations, ice-sheet advances and retreats in western patagonia shaped land fragmentation, and the formation of islands and fjords along the pacific coast  <cit> .
throughout the last glacial maximum , ice sheets extended from 56°s up to 35°s along the andes mountain range, and to 41°s in lowland areas and at sea-level in south america  <cit> .
the southern fjords and channels  of southern chile  were covered by an extensive ice sheet.
the pollen record indicates major shifts in most species of plants in this region  <cit> .
as a result of the expected pattern of species range contractions and subsequent expansion after postglacial ice sheet retreat, genetic signatures of low diversity and demographic expansions are expected in these newly colonized areas  <cit> .
this is the case for a vast majority of sigmodontine rodents, either from lowlands, patagonian mountains or tierra de fuego  <cit> .
other studies, however, have suggested the persistence of freshwater species in the region throughout the last glacial period  <cit> .
in the case of the patagonian freshwater fish galaxias platei, distributed along the western side of the andes , survival could have occurred in a southern refuge, possibly due to discontinuities of the ice field  <cit> .
likewise, southern refugia of the freshwater crab aegla alacalufi has been suggested at the el amarillo hot springs, which was possibly left uncovered by glacial ice, but isolated from more northern refugia  <cit> .
thus, it seems populations of freshwater species were not completely extirpated by ice cover in southern areas.
this alternative pattern addresses important questions about the history of this region, the magnitude of the ice cover extension, and the consequences for the phylogeographic pattern of species associated with the sfc system for which there is a complete absence of data.
in chile, the northern glacial limit is also a boundary for major environmental changes , resulting in a major biogeographic transition for marine  <cit>  and freshwater species  <cit> .
this biogeographic boundary is marked by the oceanic divergence between the humboldt and the cape horn current systems, contributing to the relative isolation of the sfc marine fauna  <cit> .
the sfc ecosystem is supplied by marine currents such as those from subantarctic waters.
these become mixed with freshwater from abundant precipitation, river flow, and glacial meltwater  <cit> .
across the biogeographic boundary, hydrographic conditions change in freshwater environments, based on the weather, slope and lithology  <cit> .
two provinces are described for the biogeography of the freshwater fishes: the chilean province, with a southern area of endemism from valdivia river to chiloé island, and the patagonian province, which is restricted to the western watersheds of continental chiloé south to tierra del fuego  <cit> .
different scenarios altered the species assemblages across biogeographical limits and likely generated distinct patterns of genetic diversity.
phylogeographic inferences should help to understand species evolutionary history across marine and freshwater biogeographic breaks.
here we use molecular markers to address the hypothesis of glacial refugia on an aquatic top predator, the huillín or southern river otter, lontra provocax.
l. provocax populations extend across the biogeographic boundary, from exclusively freshwater habitat north of this limit to mostly marine habitat in the sfc .
l. provocax has the smallest geographical range of all otter species  <cit> , being distributed across the andean-patagonian region of southern chile and part of argentina  <cit> .
in continental waters, the species is found from the toltén river basin  to the south of chile  <cit> .
the range extends through a limited area of the andes mountains into nahuel huapi national park and limay river in argentina  <cit> .
south of chiloé island in chile  the species also occurs in marine habitats along the sfc south to 56°s  <cit> .
in the sfc south of taitao peninsula , however, its distribution becomes exclusively marine  <cit> .
as l. provocax is highly dependent on the availability of crustacean prey  <cit> , the absence of the species from continental waters is related to the absence of crustaceans in the oligotrophic waters.
l. provocax is solitary with intrasexual territoriality and an average home range of  <dig>  km in rivers  <cit> .
dispersal is limited, as shown by radio-tracked otters, where the only long-range movement was made by a juvenile male that migrated  <dig> km downstream after release  <cit> .
in freshwater, occurrence of l. provocax is dependent on crustacean distribution, which is strongly influenced by the river slope and altitude  <cit> .
consequently, the distribution of l. provocax along rivers is mostly concentrated below  <dig> m altitude  <cit> , which may limit gene flow across the andes mountain range.
as crustacean and freshwater fishes were able to live in southern chile during the lgm, l. provocax could have had a sufficient amount of food to survive during this period.
however, the occurrence of l. provocax is also dependent on riparian vegetation.
terrestrial plants such as fitzroya cupressoides and hypochaeris palustris were mostly absent along ice sheet cover areas during lgm  <cit> .
thus, the question of l. provocax persistence during lgm on southern chile likely implies a trade off between availability of food and terrestrial habitat.
so far the species has been studied mainly in the rivers and lakes at the northern limit of its distribution.
comparisons between specimens living in freshwater and marine environment are restricted to diet, and consists mainly of crustaceans such as samastacus sp.
and aegla sp.
in rivers and lakes  <cit> , shifting to marine fishes of the genus patagonotothen sp., crustaceans  <cit>  and sea urchins in the marine environment.
although otter species occur mostly in freshwater habitats, the majority of these species have also been recorded in coastal environments  <cit> .
otters distributed along the coast, however, also need access to fresh water for drinking and washing their dense fur to remove accumulating salt and maintain thermo-insulation  <cit> .
fresh water is abundant along the sfc and l. provocax has also been recorded in inland rivers, such as queulat river.
it is important to note that, among all otter species, only two are exclusively adapted to marine habitats, i.e.
the north pacific ocean sea otter  and the chungungo or south pacific marine otter .
lontra felina recently diverged from l. provocax, possibly from populations in the sfc that progressively adapted to the coastal marine habitat  <cit> .
the present study analysed the phylogeographic pattern and the population structure of l. provocax, based on the mitochondrial dna sequences from control region , the nadh dehydrogenase subunit  <dig>  and the cytochrome b gene .
we aimed to infer: i) the demographic processes associated with the lgm south of the pleistocene ice cover limit for l. provocax populations; ii) the evolutionary relationship between freshwater and marine populations; iii) the population structure within each habitat type, among different continental river basins and across the andes mountain range, comparing l. provocax populations from chile and argentina.
RESULTS
a total of  <dig> bp from the mitochondrial dna control region ,  <dig> bp for nd <dig> and  <dig>  bp of the cytochrome b  were sequenced from a total of  <dig> samples.
a total of  <dig> cr haplotypes,  <dig> nd <dig> haplotypes and  <dig> cyt-b haplotypes were found for lontra provocax .
although one cr haplotype did not illustrate a clear geographic separation, nd <dig> and cyt-b haplotypes showed a clear partition according to different environments.
three nd <dig> haplotypes  were exclusively from continental freshwater , whereas haplotype iv was distributed from the chiloé island rivers  to the southern fjords and channels  and nahuel huapi area in argentina.
haplotypes v and vi were exclusively found in sfc .
on the other hand, the cyt-b haplotypes  found in both continental and chiloé island freshwater habitats , were not shared with sfc marine habitat, which showed a high diversity of cyt-b haplotypes .
the partition homogeneity test  indicated statistical congruence between the three mtdna sequences , therefore concatenated sequences were used for data analysis.
seventeen haplotypes and  <dig> polymorphic sites were observed from the concatenated mtdna sequences.
sampling sites are divided by habitat type  and locality numbers , including geographic coordinates, sample size , and concatenated haplotypes found for each location.
the concatenated haplotype a-i-a  was widely distributed in the continental freshwater habitat, except in petrohue where haplotype c-iii-a was observed.
haplotypes a-i-b and b-ii-a were also found only at mahuidanchi river and queule river respectively.
a unique haplotype a-iv-i was found on nahuel huapi area in argentina, not detected on chilean rivers.
all samples from rivers on chiloé island had a unique haplotype a-iv-a, which was not found in continental or marine sites.
in the sfc, each haplotype was restricted to one or few locations.
the samova revealed increasing Φct values when the numbers of groups increased .
petrohue river was separated out in all samova partitions.
the most geographically coherent partition could be defined in  <dig> groups:  most continental rivers and lakes locations including chiloé island ,  continental petrohue river ,  fjords, channels and queulat river  and  southern-most fjords and channels .
high population structure  was found among  <dig> populations and the four groups defined.
no haplotype was shared between the four groups.
we found an overall high haplotype diversity  but low nucleotide diversity  for l. provocax.
the expected pattern of a lower diversity for the glaciated region was not observed, since the haplotype diversity of the southern group  was high and comparable to that of the non-glaciated area, such as the continental freshwater habitat  or the entire northern region .
sample size, polymorphic sites , number of haplotypes  for each of mtdna sequences , haplotype diversity  and nucleotide diversity  for the concatenated sequences cr+nd5+cyt-b by habitat.
the model selected for the three concatenated genes was trn+i, selected based on the lower value of akaike information criterion .
the bayesian phylogenetic analysis  showed two divergent clades including haplotypes described for continental freshwater habitat, whereas haplotypes from chiloé island  and nahuel huapi  belonged to a polytomous group of mixed origins.
a continental freshwater clade is represented by two haplotypes widely distributed along continental rivers and lakes , and another one by two haplotypes  found on queule  and petrohue rivers .
another clade was represented by the majority of haplotypes from the sfc.
one of the sfc haplotypes  was found within the paraphyletic assemblage consisting of freshwater haplotypes.
the divergence between l. felina and l. provocax was  <dig> %, whereas it reached  <dig> % between l. felina and l. longicaudis, and  <dig> % between l. longicaudis and l. provocax.
the mean distance within the l. provocax clade was  <dig> %.
l. provocax mean distance between clades varied from  <dig> % between a freshwater clade  and the sfc clade,  <dig> % among continental freshwater clades , to  <dig> % among the freshwater lineages .
a signature of recent expansion was detected from the sfc clade, as evidenced by the significant fu's test , although tajima was not significant .
the observed mismatch distribution was not significantly different from that expected for both models, however the spatial expansion model fit the sfc-clade data better than demographic expansion model .
in contrast, the median joining network  topology for the entire distribution of the sfc haplotypes, which differs from that expected in a recent population expansion .
mjn indicates high divergence among haplotypes from freshwater habitat and a lack of intermediate haplotypes.
moreover, the bayesian skyline plot did not show any evidence of historical population growth for the populations from the sfc clade, but rather a recent population decrease .
no signature of recent expansion was detected for l. provocax from the distribution in non-glaciated freshwater regions .
discussion
our data evidenced a strong genetic differentiation between continental freshwater and sfc regions.
the limits of the latter correspond not only to a habitat change, but also to a major biogeographic break for marine  <cit>  and freshwater  <cit>  species, and to the northern limit of coastal lgm ice cover  <cit> .
southern marine fjords and channels: glacial survival?
population displacements followed by founder effects due to recolonization of southern areas from a northern refuge would produce a signature of reduced genetic diversity, when compared to northern areas.
the latter pattern was described for southern bull kelp , which has reduced genetic diversity in southern chile with a clear signature of postglacial expansion  <cit> .
in our case, phylogenetic reconstruction for l. provocax indicates that most haplotypes from sfc form a distinctive haplogroup emerging from a basal northern haplotype found in freshwater habitat, which could suggest a postglacial recolonization from the northern freshwater refuge.
however, several results suggest a different scenario: a non-starlike network of haplotypes from sfc, an absence of historical population growth signature from the skyline plot analysis, high haplotype diversity, and highly divergent lineages.
these results are not in agreement with a scenario of post-glacial expansion from northern freshwater populations.
they rather support the hypothesis of a persistence of the species in this region during lgm.
the hypothesis of persistence in glaciated areas is most often rejected by phylogeographic studies.
some studies have, however, shown phylogeographic results strongly supporting demographic persistence in areas supposedly covered by ice  <cit> .
in patagonia, iso-pollen lines indicate the existence of a terrestrial southern refuge  <cit>  along the south-eastern coast of the island of tierra del fuego, east from the beagle channel.
the presence of a southern refuge would explain the present day plant species disjunction between chiloé and magallanes, and the presence of bryophyta species and vascular herbaceous plants not distributed further north of 47-48°s  <cit> .
nevertheless, the unique southern glacial refugium described for terrestrial species was not supported by our data.
single or multiple southern refugia and the persistence of the species within glaciated areas in chilean patagonia were debated for terrestrial  <cit>  and freshwater species  <cit> .
in the case of temperate forest species, patches persisted either at the northern limits of ice cover in chile or on the eastern limits of ice sheet in argentinean patagonia  <cit> .
all these examples point to a post-glacial colonization of western patagonia, i.e.
the chilean side of the andes.
the survival of l. provocax in non-glaciated areas of the eastern side of the andes would have required a recolonization of the species across the andes at multiple sites in order to allow the re-introduction of such a high genetic diversity.
in addition, because the species currently occurs only in the nahuel huapi and limay river area in continental argentina  <cit> , subsequent localised extinctions of l. provocax along most of the eastern refugia would be required to support such a scenario.
lastly, the unique argentinean haplotype  is derived from haplotype a.iv.a present in chiloe island, suggesting that its distribution was likely widespread in the recent past, and that range expansion was more likely from chile to the eastern side of the andes.
thus, the results do not support an origin of l. provocax diversity from eastern patagonia.
on the other hand, glacial refugia on western coast of the andes were suggested for some freshwater species.
mtdna sequences of the fishes galaxias platei and g. maculatus along the western side of the andes mountain  suggest that these species survived in small southern sites due to discontinuities of the ice field  <cit> .
also, exposed portions of the pacific continental shelf could have constituted favourable environment for such aquatic species  <cit> .
similarly, the freshwater crab aegla alacalufi seem to have survived in glaciated areas, at least in a site identified as the el amarillo hot springs, which was possibly left uncovered by glacial ice  <cit> .
whether multiple refugia existed or l. provocax survived all along the western southern patagonia during lgm is still a matter of debate.
the persistence of l. provocax is dependent on terrestrial habitat for dens  and on aquatic habitat for food availability.
such habitat was probably available on the eastern side of the andes, as suggested by the extension of the distribution of fitzroya cupressoides  <cit> , allowing the survival of l. provocax in a large area and therefore allowing the persistence of high genetic diversity.
in the sfc, l. provocax is known to feed not only on crustaceans and sea urchins, but also on intertidal and subtidal fishes   <cit> .
ice scour can eliminate intertidal and shallow water benthos in the southern ocean  <cit> .
in the case of complete elimination of intertidal resources, diet could have been based on species such as the freshwater fish , catfish  or the freshwater crab , which survived in the area during lgm  <cit> , or subtidal organisms.
otter species, such as north american river otter , the eurasian otter  and the sea otter , are distributed throughout extreme cold environments.
l. canadensis inhabiting the marine environments in alaska has access to two major types of prey: intertidal-demersal organisms such as fishes  and crustaceans, and seasonally available schooling pelagic fishes  <cit> .
similarly, l. provocax populations could have shifted their diet according to prey availability and thus persisted during the lgm in the patagonian sfc.
although species such as the eurasian otter  show evidence of a unique glacial refuge and low genetic diversity  <cit> , other mustelids were able to survive during lgm.
gulo gulo, mustela nivalis and mustela erminea show adaptations for survival in pleistocene conditions  <cit> .
southern glacial refugia or adaptation to marine habitat
our data show a monophyletic haplogroup for most haplotypes from the sfc range of l. provocax , distinct from freshwater habitat haplotypes.
such differences between the l. provocax populations inhabiting the freshwater and marine environments suggest either past genetic isolation, and/or restricted gene flow between them at the present time, allowing genetic drift or natural selection to operate.
changes in the l. provocax diet and a greater ability to swim larger distances as required in the sfc could eventually lead to adaptation to the marine environment, if plasticity is limited.
however, patterns of genetic diversity generated by gene surfing during recolonization are similar to those generated by selection and could thus be mistakenly interpreted as adaptive events  <cit> .
similarly to selection and unlike most other demographic effects, gene surfing generally does not affect all loci, and thus seems especially difficult to distinguish from directional selection  <cit> .
all otter species, except the sea otter  and the chungungo , are dependent of freshwater habitat  <cit> .
nevertheless, freshwater sources are abundant in sfc and l. provocax can temporarily return to rivers to access a supply of fresh water.
moreover, all but four freshwater otter species have also been recorded along the coast  <cit>  .
nevertheless, no genetic surveys have been conducted to determine divergence of otter lineages from different environments.
lontra phylogeny based on mtdna markers revealed the recent divergence between l. provocax and l. felina about  <dig>  years ago  with a possible speciation of l. felina from l. provocax living on sfc  <cit> .
this speciation scenario is in agreement with the adaptation hypothesis of l. provocax to the marine habitat in sfc.
conservation implications
although higher haplotype diversity was expected along northern populations due to the persistence of rivers and forests, low haplotype diversity  but high divergence among haplotypes was observed.
our results are concordant with the hypothesis of a recent loss of genetic diversity in freshwater environments due to hunting and habitat destruction.
this is specifically supported by: i) a.i.a haplotype shared among several locations; ii) two highly divergent clades; iii) two divergent haplotypes  in the queule river.
genetic theory predicts that during population bottlenecks low frequency alleles are lost by genetic drift  <cit> .
similarly, l. provocax intermediate haplotypes are not seen on the mjn, suggesting that they were eliminated within the freshwater range.
this pattern is consistent with the history of the l. provocax populations in the region.
indeed, l. provocax populations have been eliminated from the north of its past range.
the northern limit of distribution changed from cauquenes and cachapoal rivers  to tolten river basin   <cit> .
its small geographical range has been strongly impacted by anthropogenic activities resulting in a decline to less than 10% of its former distribution in freshwater habitats.
l. provocax was intensively hunted for its fur; and hunting continued until the 1970's in some southern localities  <cit> .
furthermore, the species activity has been significantly reduced in areas where riparian vegetation was removed or watercourses were disturbed or recently polluted by pulp factories  <cit> .
riparian vegetation significantly influences the presence of crustaceans and consequently the occurrence of l. provocax in the area  <cit> .
the reduction in the distribution of l. provocax led to the classification of the species in chile as ''endangered'' in the northern range between the o'higgins and los lagos regions, in rivers and lakes, and ''insufficiently known'' for the aysén and magallanes regions, where distribution was largely marine  <cit> .
although l. provocax in freshwater habitats mostly occur below  <dig> m altitude  <cit> , the majority of national parks and reserves  in south-central chile  are located above  <dig> m altitude  <cit> , and therefore do not serve the conservation of this species.
our data shows that the survival of the species along sfc during glacial cycles maintained a high diversity along sfc.
large-scale temperate deforestation in chile has progressed from north to south  <cit> .
human populations are concentrated in the central-south region of chile, and are less dense south of chiloé island.
thus, southern l. provocax populations have not been greatly impacted by anthropogenic actions, and have maintained high genetic diversity compared to northern freshwater populations.
l. provocax populations along sfc are, however, barely studied, and increasing human activities in this area are a potential threat to these populations.
CONCLUSIONS
our results evidenced the persistence of a semi-aquatic carnivore species, the huillín, in western patagonia along areas covered by ice sheet during lgm.
marine habitat of the sfc played an important role for l. provocax survival during lgm, probably associated to the survival of other freshwater and marine species that may have represented a persistent food source for the huillín.
therefore, genetic differentiation between northern exclusively freshwater habitat dominated by riparian vegetation and sfc may be explained by some ecological differentiation between both kinds of habitats.
this is an interesting clue for understanding why so many aquatic species seem to have persisted in glaciated areas.

 climate change caused substantial alterations of the landscape and sea level, influencing patterns of species distribution. pleistocene glaciations, ice-sheet advances and retreats in western patagonia shaped land fragmentation, and the formation of islands and fjords along the pacific coast. throughout the last glacial maximum, ice sheets extended from 56°s up to 35°s along the andes mountain range, and to 41°s in lowland areas and at sea-level in south america. the southern fjords and channels of southern chile were covered by an extensive ice sheet. the pollen record indicates major shifts in most species of plants in this region. as a result of the expected pattern of species range contractions and subsequent expansion after postglacial ice sheet retreat, genetic signatures of low diversity and demographic expansions are expected in these newly colonized areas. this is the case for a vast majority of sigmodontine rodents, either from lowlands, patagonian mountains or tierra de fuego. other studies, however, have suggested the persistence of freshwater species in the region throughout the last glacial period. in the case of the patagonian freshwater fish galaxias platei, distributed along the western side of the andes, survival could have occurred in a southern refuge, possibly due to discontinuities of the ice field. likewise, southern refugia of the freshwater crab aegla alacalufi has been suggested at the el amarillo hot springs, which was possibly left uncovered by glacial ice, but isolated from more northern refugia. thus, it seems populations of freshwater species were not completely extirpated by ice cover in southern areas. this alternative pattern addresses important questions about the history of this region, the magnitude of the ice cover extension, and the consequences for the phylogeographic pattern of species associated with the sfc system for which there is a complete absence of data. in chile, the northern glacial limit is also a boundary for major environmental changes, resulting in a major biogeographic transition for marine and freshwater species. this biogeographic boundary is marked by the oceanic divergence between the humboldt and the cape horn current systems, contributing to the relative isolation of the sfc marine fauna. the sfc ecosystem is supplied by marine currents such as those from subantarctic waters. these become mixed with freshwater from abundant precipitation, river flow, and glacial meltwater. across the biogeographic boundary, hydrographic conditions change in freshwater environments, based on the weather, slope and lithology. two provinces are described for the biogeography of the freshwater fishes: the chilean province, with a southern area of endemism from valdivia river to chiloé island, and the patagonian province, which is restricted to the western watersheds of continental chiloé south to tierra del fuego. different scenarios altered the species assemblages across biogeographical limits and likely generated distinct patterns of genetic diversity. phylogeographic inferences should help to understand species evolutionary history across marine and freshwater biogeographic breaks. here we use molecular markers to address the hypothesis of glacial refugia on an aquatic top predator, the huillín or southern river otter, lontra provocax. l. provocax populations extend across the biogeographic boundary, from exclusively freshwater habitat north of this limit to mostly marine habitat in the sfc. l. provocax has the smallest geographical range of all otter species, being distributed across the andean-patagonian region of southern chile and part of argentina. in continental waters, the species is found from the toltén river basin to the south of chile. the range extends through a limited area of the andes mountains into nahuel huapi national park and limay river in argentina. south of chiloé island in chile the species also occurs in marine habitats along the sfc south to 56°s. in the sfc south of taitao peninsula, however, its distribution becomes exclusively marine. as l. provocax is highly dependent on the availability of crustacean prey, the absence of the species from continental waters is related to the absence of crustaceans in the oligotrophic waters. l. provocax is solitary with intrasexual territoriality and an average home range of km in rivers. dispersal is limited, as shown by radio-tracked otters, where the only long-range movement was made by a juvenile male that migrated km downstream after release. in freshwater, occurrence of l. provocax is dependent on crustacean distribution, which is strongly influenced by the river slope and altitude. consequently, the distribution of l. provocax along rivers is mostly concentrated below m altitude, which may limit gene flow across the andes mountain range. as crustacean and freshwater fishes were able to live in southern chile during the lgm, l. provocax could have had a sufficient amount of food to survive during this period. however, the occurrence of l. provocax is also dependent on riparian vegetation. terrestrial plants such as fitzroya cupressoides and hypochaeris palustris were mostly absent along ice sheet cover areas during lgm. thus, the question of l. provocax persistence during lgm on southern chile likely implies a trade off between availability of food and terrestrial habitat. so far the species has been studied mainly in the rivers and lakes at the northern limit of its distribution. comparisons between specimens living in freshwater and marine environment are restricted to diet, and consists mainly of crustaceans such as samastacus sp. and aegla sp. in rivers and lakes, shifting to marine fishes of the genus patagonotothen sp., crustaceans and sea urchins in the marine environment. although otter species occur mostly in freshwater habitats, the majority of these species have also been recorded in coastal environments. otters distributed along the coast, however, also need access to fresh water for drinking and washing their dense fur to remove accumulating salt and maintain thermo-insulation. fresh water is abundant along the sfc and l. provocax has also been recorded in inland rivers, such as queulat river. it is important to note that, among all otter species, only two are exclusively adapted to marine habitats, i.e. the north pacific ocean sea otter and the chungungo or south pacific marine otter. lontra felina recently diverged from l. provocax, possibly from populations in the sfc that progressively adapted to the coastal marine habitat. the present study analysed the phylogeographic pattern and the population structure of l. provocax, based on the mitochondrial dna sequences from control region, the nadh dehydrogenase subunit and the cytochrome b gene. we aimed to infer: i) the demographic processes associated with the lgm south of the pleistocene ice cover limit for l. provocax populations; ii) the evolutionary relationship between freshwater and marine populations; iii) the population structure within each habitat type, among different continental river basins and across the andes mountain range, comparing l. provocax populations from chile and argentina.  a total of bp from the mitochondrial dna control region, bp for nd and bp of the cytochrome b were sequenced from a total of samples. a total of cr haplotypes, nd haplotypes and cyt-b haplotypes were found for lontra provocax. although one cr haplotype did not illustrate a clear geographic separation, nd and cyt-b haplotypes showed a clear partition according to different environments. three nd haplotypes were exclusively from continental freshwater, whereas haplotype iv was distributed from the chiloé island rivers to the southern fjords and channels and nahuel huapi area in argentina. haplotypes v and vi were exclusively found in sfc. on the other hand, the cyt-b haplotypes found in both continental and chiloé island freshwater habitats, were not shared with sfc marine habitat, which showed a high diversity of cyt-b haplotypes. the partition homogeneity test indicated statistical congruence between the three mtdna sequences, therefore concatenated sequences were used for data analysis. seventeen haplotypes and polymorphic sites were observed from the concatenated mtdna sequences. sampling sites are divided by habitat type and locality numbers, including geographic coordinates, sample size, and concatenated haplotypes found for each location. the concatenated haplotype a-i-a was widely distributed in the continental freshwater habitat, except in petrohue where haplotype c-iii-a was observed. haplotypes a-i-b and b-ii-a were also found only at mahuidanchi river and queule river respectively. a unique haplotype a-iv-i was found on nahuel huapi area in argentina, not detected on chilean rivers. all samples from rivers on chiloé island had a unique haplotype a-iv-a, which was not found in continental or marine sites. in the sfc, each haplotype was restricted to one or few locations. the samova revealed increasing Φct values when the numbers of groups increased. petrohue river was separated out in all samova partitions. the most geographically coherent partition could be defined in groups: most continental rivers and lakes locations including chiloé island, continental petrohue river, fjords, channels and queulat river and southern-most fjords and channels. high population structure was found among populations and the four groups defined. no haplotype was shared between the four groups. we found an overall high haplotype diversity but low nucleotide diversity for l. provocax. the expected pattern of a lower diversity for the glaciated region was not observed, since the haplotype diversity of the southern group was high and comparable to that of the non-glaciated area, such as the continental freshwater habitat or the entire northern region. sample size, polymorphic sites, number of haplotypes for each of mtdna sequences, haplotype diversity and nucleotide diversity for the concatenated sequences cr+nd5+cyt-b by habitat. the model selected for the three concatenated genes was trn+i, selected based on the lower value of akaike information criterion. the bayesian phylogenetic analysis showed two divergent clades including haplotypes described for continental freshwater habitat, whereas haplotypes from chiloé island and nahuel huapi belonged to a polytomous group of mixed origins. a continental freshwater clade is represented by two haplotypes widely distributed along continental rivers and lakes, and another one by two haplotypes found on queule and petrohue rivers. another clade was represented by the majority of haplotypes from the sfc. one of the sfc haplotypes was found within the paraphyletic assemblage consisting of freshwater haplotypes. the divergence between l. felina and l. provocax was %, whereas it reached % between l. felina and l. longicaudis, and % between l. longicaudis and l. provocax. the mean distance within the l. provocax clade was %. l. provocax mean distance between clades varied from % between a freshwater clade and the sfc clade, % among continental freshwater clades, to % among the freshwater lineages. a signature of recent expansion was detected from the sfc clade, as evidenced by the significant fu's test, although tajima was not significant. the observed mismatch distribution was not significantly different from that expected for both models, however the spatial expansion model fit the sfc-clade data better than demographic expansion model. in contrast, the median joining network topology for the entire distribution of the sfc haplotypes, which differs from that expected in a recent population expansion. mjn indicates high divergence among haplotypes from freshwater habitat and a lack of intermediate haplotypes. moreover, the bayesian skyline plot did not show any evidence of historical population growth for the populations from the sfc clade, but rather a recent population decrease. no signature of recent expansion was detected for l. provocax from the distribution in non-glaciated freshwater regions. discussion our data evidenced a strong genetic differentiation between continental freshwater and sfc regions. the limits of the latter correspond not only to a habitat change, but also to a major biogeographic break for marine and freshwater species, and to the northern limit of coastal lgm ice cover. southern marine fjords and channels: glacial survival? population displacements followed by founder effects due to recolonization of southern areas from a northern refuge would produce a signature of reduced genetic diversity, when compared to northern areas. the latter pattern was described for southern bull kelp, which has reduced genetic diversity in southern chile with a clear signature of postglacial expansion. in our case, phylogenetic reconstruction for l. provocax indicates that most haplotypes from sfc form a distinctive haplogroup emerging from a basal northern haplotype found in freshwater habitat, which could suggest a postglacial recolonization from the northern freshwater refuge. however, several results suggest a different scenario: a non-starlike network of haplotypes from sfc, an absence of historical population growth signature from the skyline plot analysis, high haplotype diversity, and highly divergent lineages. these results are not in agreement with a scenario of post-glacial expansion from northern freshwater populations. they rather support the hypothesis of a persistence of the species in this region during lgm. the hypothesis of persistence in glaciated areas is most often rejected by phylogeographic studies. some studies have, however, shown phylogeographic results strongly supporting demographic persistence in areas supposedly covered by ice. in patagonia, iso-pollen lines indicate the existence of a terrestrial southern refuge along the south-eastern coast of the island of tierra del fuego, east from the beagle channel. the presence of a southern refuge would explain the present day plant species disjunction between chiloé and magallanes, and the presence of bryophyta species and vascular herbaceous plants not distributed further north of 47-48°s. nevertheless, the unique southern glacial refugium described for terrestrial species was not supported by our data. single or multiple southern refugia and the persistence of the species within glaciated areas in chilean patagonia were debated for terrestrial and freshwater species. in the case of temperate forest species, patches persisted either at the northern limits of ice cover in chile or on the eastern limits of ice sheet in argentinean patagonia. all these examples point to a post-glacial colonization of western patagonia, i.e. the chilean side of the andes. the survival of l. provocax in non-glaciated areas of the eastern side of the andes would have required a recolonization of the species across the andes at multiple sites in order to allow the re-introduction of such a high genetic diversity. in addition, because the species currently occurs only in the nahuel huapi and limay river area in continental argentina, subsequent localised extinctions of l. provocax along most of the eastern refugia would be required to support such a scenario. lastly, the unique argentinean haplotype is derived from haplotype a.iv.a present in chiloe island, suggesting that its distribution was likely widespread in the recent past, and that range expansion was more likely from chile to the eastern side of the andes. thus, the results do not support an origin of l. provocax diversity from eastern patagonia. on the other hand, glacial refugia on western coast of the andes were suggested for some freshwater species. mtdna sequences of the fishes galaxias platei and g. maculatus along the western side of the andes mountain suggest that these species survived in small southern sites due to discontinuities of the ice field. also, exposed portions of the pacific continental shelf could have constituted favourable environment for such aquatic species. similarly, the freshwater crab aegla alacalufi seem to have survived in glaciated areas, at least in a site identified as the el amarillo hot springs, which was possibly left uncovered by glacial ice. whether multiple refugia existed or l. provocax survived all along the western southern patagonia during lgm is still a matter of debate. the persistence of l. provocax is dependent on terrestrial habitat for dens and on aquatic habitat for food availability. such habitat was probably available on the eastern side of the andes, as suggested by the extension of the distribution of fitzroya cupressoides, allowing the survival of l. provocax in a large area and therefore allowing the persistence of high genetic diversity. in the sfc, l. provocax is known to feed not only on crustaceans and sea urchins, but also on intertidal and subtidal fishes. ice scour can eliminate intertidal and shallow water benthos in the southern ocean. in the case of complete elimination of intertidal resources, diet could have been based on species such as the freshwater fish, catfish or the freshwater crab, which survived in the area during lgm, or subtidal organisms. otter species, such as north american river otter, the eurasian otter and the sea otter, are distributed throughout extreme cold environments. l. canadensis inhabiting the marine environments in alaska has access to two major types of prey: intertidal-demersal organisms such as fishes and crustaceans, and seasonally available schooling pelagic fishes. similarly, l. provocax populations could have shifted their diet according to prey availability and thus persisted during the lgm in the patagonian sfc. although species such as the eurasian otter show evidence of a unique glacial refuge and low genetic diversity, other mustelids were able to survive during lgm. gulo gulo, mustela nivalis and mustela erminea show adaptations for survival in pleistocene conditions. southern glacial refugia or adaptation to marine habitat our data show a monophyletic haplogroup for most haplotypes from the sfc range of l. provocax, distinct from freshwater habitat haplotypes. such differences between the l. provocax populations inhabiting the freshwater and marine environments suggest either past genetic isolation, and/or restricted gene flow between them at the present time, allowing genetic drift or natural selection to operate. changes in the l. provocax diet and a greater ability to swim larger distances as required in the sfc could eventually lead to adaptation to the marine environment, if plasticity is limited. however, patterns of genetic diversity generated by gene surfing during recolonization are similar to those generated by selection and could thus be mistakenly interpreted as adaptive events. similarly to selection and unlike most other demographic effects, gene surfing generally does not affect all loci, and thus seems especially difficult to distinguish from directional selection. all otter species, except the sea otter and the chungungo, are dependent of freshwater habitat. nevertheless, freshwater sources are abundant in sfc and l. provocax can temporarily return to rivers to access a supply of fresh water. moreover, all but four freshwater otter species have also been recorded along the coast. nevertheless, no genetic surveys have been conducted to determine divergence of otter lineages from different environments. lontra phylogeny based on mtdna markers revealed the recent divergence between l. provocax and l. felina about years ago with a possible speciation of l. felina from l. provocax living on sfc. this speciation scenario is in agreement with the adaptation hypothesis of l. provocax to the marine habitat in sfc. conservation implications although higher haplotype diversity was expected along northern populations due to the persistence of rivers and forests, low haplotype diversity but high divergence among haplotypes was observed. our results are concordant with the hypothesis of a recent loss of genetic diversity in freshwater environments due to hunting and habitat destruction. this is specifically supported by: i) a.i.a haplotype shared among several locations; ii) two highly divergent clades; iii) two divergent haplotypes in the queule river. genetic theory predicts that during population bottlenecks low frequency alleles are lost by genetic drift. similarly, l. provocax intermediate haplotypes are not seen on the mjn, suggesting that they were eliminated within the freshwater range. this pattern is consistent with the history of the l. provocax populations in the region. indeed, l. provocax populations have been eliminated from the north of its past range. the northern limit of distribution changed from cauquenes and cachapoal rivers to tolten river basin. its small geographical range has been strongly impacted by anthropogenic activities resulting in a decline to less than 10% of its former distribution in freshwater habitats. l. provocax was intensively hunted for its fur; and hunting continued until the 1970's in some southern localities. furthermore, the species activity has been significantly reduced in areas where riparian vegetation was removed or watercourses were disturbed or recently polluted by pulp factories. riparian vegetation significantly influences the presence of crustaceans and consequently the occurrence of l. provocax in the area. the reduction in the distribution of l. provocax led to the classification of the species in chile as''endangered'' in the northern range between the o'higgins and los lagos regions, in rivers and lakes, and''insufficiently known'' for the aysén and magallanes regions, where distribution was largely marine. although l. provocax in freshwater habitats mostly occur below m altitude, the majority of national parks and reserves in south-central chile are located above m altitude, and therefore do not serve the conservation of this species. our data shows that the survival of the species along sfc during glacial cycles maintained a high diversity along sfc. large-scale temperate deforestation in chile has progressed from north to south. human populations are concentrated in the central-south region of chile, and are less dense south of chiloé island. thus, southern l. provocax populations have not been greatly impacted by anthropogenic actions, and have maintained high genetic diversity compared to northern freshwater populations. l. provocax populations along sfc are, however, barely studied, and increasing human activities in this area are a potential threat to these populations.  our results evidenced the persistence of a semi-aquatic carnivore species, the huillín, in western patagonia along areas covered by ice sheet during lgm. marine habitat of the sfc played an important role for l. provocax survival during lgm, probably associated to the survival of other freshwater and marine species that may have represented a persistent food source for the huillín. therefore, genetic differentiation between northern exclusively freshwater habitat dominated by riparian vegetation and sfc may be explained by some ecological differentiation between both kinds of habitats. this is an interesting clue for understanding why so many aquatic species seem to have persisted in glaciated areas. 
here, we further test this hypothesis based on phylogeographic inferences from a semi-aquatic species that is a top predator of river and marine fauna, the huillín or southern river otter. a number of studies have described the extension of ice cover in western patagonia during the last glacial maximum, providing evidence of a complete cover of terrestrial habitat from 41°s to 56°s and two main refugia, one in south-eastern tierra del fuego and the other north of the chiloé island. our data support the hypothesis of the persistence of l. provocax in western patagonia, south of the ice sheet limit, during last glacial maximum. we examined mtdna sequences of the control region, nd and cytochrome-b in samples of l. provocax from locations in river and marine habitats. a highly diverse clade was represented by haplotypes from the marine habitat of the southern fjords and channels region, whereas only one of these haplotypes was paraphyletic and associated with northern river haplotypes. phylogenetic analysis illustrates two main divergent clades for l. provocax in continental freshwater habitat. this limit also corresponds to a strong environmental change, which might have spurred l. provocax differentiation between the two environments. however, recent evidence of high genetic diversity in patagonian river species suggests the existence of aquatic refugia in this region.
 a number of studies have described the extension of ice cover in western patagonia during the last glacial maximum, providing evidence of a complete cover of terrestrial habitat from 41°s to 56°s and two main refugia, one in south-eastern tierra del fuego and the other north of the chiloé island. however, recent evidence of high genetic diversity in patagonian river species suggests the existence of aquatic refugia in this region. here, we further test this hypothesis based on phylogeographic inferences from a semi-aquatic species that is a top predator of river and marine fauna, the huillín or southern river otter.  we examined mtdna sequences of the control region, nd and cytochrome-b in samples of l. provocax from locations in river and marine habitats. phylogenetic analysis illustrates two main divergent clades for l. provocax in continental freshwater habitat. a highly diverse clade was represented by haplotypes from the marine habitat of the southern fjords and channels region, whereas only one of these haplotypes was paraphyletic and associated with northern river haplotypes.  our data support the hypothesis of the persistence of l. provocax in western patagonia, south of the ice sheet limit, during last glacial maximum. this limit also corresponds to a strong environmental change, which might have spurred l. provocax differentiation between the two environments. 