Phylogeography of the Patagonian otter Lontra provocax: adaptive divergence to marine habitat or signature of southern glacial refugia?
© Vianna et al; licensee BioMed Central Ltd. 2011
Received: 1 September 2010
Accepted: 28 February 2011
Published: 28 February 2011
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 (Lontra provocax).
We examined mtDNA sequences of the control region, ND5 and Cytochrome-b (2151 bp in total) in 75 samples of L. provocax from 21 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 (SFC) region (43°38' to 53°08'S), 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 (41°S latitude). This limit also corresponds to a strong environmental change, which might have spurred L. provocax differentiation between the two environments.
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 [21–23]. In continental waters, the species is found from the Toltén River basin (39°S latitude) 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 [24–26]. South of Chiloé Island in Chile (42°S) the species also occurs in marine habitats along the SFC south to 56°S . In the SFC south of Taitao Peninsula (46°S latitude), however, its distribution becomes exclusively marine [28, 27]. As L. provocax is highly dependent on the availability of crustacean prey [29, 30], 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 11.3 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 46 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 300 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 [34, 35]. 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 [24, 29], 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 (Enhydra lutris) and the chungungo or south Pacific marine otter (Lontra felina). 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 (CR), the NADH dehydrogenase subunit 5 (ND5) and the cytochrome b gene (Cyt-b). 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.
Sampling sites of Lontra provocax analysed in this study.
Huilio River, Chile
rivers and lakes
Queule River, Chile
Mahuidanche River, Chile
Lingue River, Chile
Cua cua River, Chile
Riñihue Lake, Chile
Traful Lake, Argentina
Nahuel Huapi Lake, Argentina
Petrohue River, Chile
Darwin, Chiloé Island
Chepu River, Chiloé Island
Seno Magdalena, Magdalena Island
Valle Marta, Magdalena Island
Puyuhuapi Channel, Magdalena Island
Madre de Dios Island
Around Puerto Natales
Around Punta Arenas
D.V.D; H.IV.F; G.VI.E
The concatenated haplotype A-I-A (for CR-ND5-Cyt-b respectively) 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 (K = 2-6; Φ CT = 0.55-0.72). Petrohue River was separated out in all SAMOVA partitions. The most geographically coherent partition could be defined in 4 groups: (i) most continental rivers and lakes locations including Chiloé Island (location 1 to 11), (ii) continental Petrohue river (location 9), (iii) fjords, channels and Queulat River (locations 11 to 18) and (iv) southern-most fjords and channels (locations 19 to 21). High population structure (Φ ST = 0.85, p < 0.0001) was found among 21 populations and the four groups defined. No haplotype was shared between the four groups.
Genetic diversity of mtDNA sequences for Lontra provocax.
RC+ND5+Cyt-b - Genetic Diversity
Continental rivers and lakes
0.001306 +/- 0.000791
Southern Fjords and Channels
0.7220 +/- 0.0531
0.001060 +/- 0.000663
0.8775 +/- 0.0195
0.001610 +/- 0.000922
The divergence between L. felina and L. provocax was 1.5%, whereas it reached 5.2% between L. felina and L. longicaudis, and 4.6% between L. longicaudis and L. provocax. The mean distance within the L. provocax clade was 0.3%. L. provocax mean distance between clades varied from 0.5% between a freshwater clade (B.II.A+C.III.A) and the SFC clade, 0.4% among continental freshwater clades (A.I.A.+A.I.B and B.II.A.+CI.II.A), to 0.1% among the freshwater lineages (A.I.A.+A.I.B, A.IV.A, A.IV.I, A.IV.G).
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 [2, 3].
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 (Durvillaea antarctica), 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 [7–9, 34, 35, 41–46] and freshwater species [10–13, 47, 49]. 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 [34, 35, 42, 44, 48]. 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 [24–26], 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 (A.IV.I) 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 (39°-49°S) suggest that these species survived in small southern sites due to discontinuities of the ice field [12, 49]. 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 (shelter) 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 (46% of its diet) . 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 (Galaxias platei and G. maculatus), catfish (Trichomycterus areolatus) or the freshwater crab (Aegla alacalufi), which survived in the area during LGM [12, 13, 47, 49], or subtidal organisms. Otter species, such as North American river otter (Lontra candensis), the Eurasian otter (Lutra lutra) and the sea otter (Enhydra lutris), 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 (Cottidae, Hexagrammidae) 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 (Lutra lutra) 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 (43°S to 53°S), 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 (Enhydra lutris) and the chungungo (Lontra felina), 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  (Lontra canadensis, Lutra lutra and Lontra longicaudis, among others). 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 883,000 years ago (95% HPD: 0.16-1.89 mya) 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.
Although higher haplotype diversity was expected along northern populations due to the persistence of rivers and forests, low haplotype diversity (compared to SFC) 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 (A.I.A and B.II.A) 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 (34°S) to Tolten River basin (39°S) . 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 [57, 58]. Riparian vegetation significantly influences the presence of crustaceans and consequently the occurrence of L. provocax in the area [30, 32, 57].
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 300 m altitude , the majority of National Parks and Reserves (> 90%) in south-central Chile (35.6° to 41.3°S) are located above 600 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.
Study area, sample collection and DNA extraction
A total of 57 feces and 18 tissue samples were collected. These included blood from captured animals, muscle from carcasses of animals that died of natural causes, and pelts confiscated by authorities due to illegal hunting. Feces and muscle tissue were preserved in pure ethanol. The sampling range included the full range of the species (Figure 1) from the northern limit in Chile, the Huilio River tributary of Tolten river (38°58'S) to southern Patagonia (53°08'S, Table 1, Figure 1, 2). In addition we sample the only known population of the eastern side of the Andes: the Nahuel Huapi Area. A total of 21 locations were sampled, nine in continental rivers and lakes (CRL), two in rivers from Chiloé Island (CI) and ten marine locations in the fjords and channels (SFC).
Mitochondrial DNA control region (CR), NADH dehydrogenase subunit 5 (ND5) gen and cythocrome b (Cyt-b) gen were amplified using primers described by : LfCR-2F and LfCR-1R; ND5-DF1 and LfND5-R; and LfCYTB-1F; LfCYTB-1F; LfCYTB-2F and LfCYTB-2R.
PCR reactions were carried out following . PCR amplification was confirmed by electrophoresis of products with ethidium bromide in 0.8% agarose gels and visualization under UV light. Amplicons were purified using QIAquick PCR purification kit (Qiagen) and sequenced using amplification primers by Macrogen Inc., Seoul, Korea. All sequences have been deposited in the GenBank [GenBank: GQ843803-GQ843824, and HM997011-HM997017].
The sequences were aligned and mutations were confirmed by eye according to the chromatogram using Proseq ver. 2.91 . All sequences were aligned and haplotypes were identified using ClustalX ver. 1.83 .
Spatial analysis of molecular variance (SAMOVA) was implemented by SAMOVA ver. 1.0  to define groups based on the geographic distribution of the genetic diversity. SAMOVA were performed for 21 locations testing from 2 to 7 groups, each of which with 100 initial conditions. The groups of populations geographically homogeneous are defined by maximizing Φ CT values (among group variance) and minimizing Φ SC values (among populations within group variance). Haplotype (h) and nucleotide diversity (π) were calculated using Arlequin program ver. 3.0  for all data set and for the different environments.
Deviations from a neutral Wright-Fisher model were performed by calculating Tajima's D and Fu's Fs statistics [65, 66]. We tested the demographic  and spatial expansion  models by calculating the sum of squared differences (SSD) between the observed and an estimated mismatch distribution obtained by 1,000 bootstrap. The P-value of the SSD statistic was calculated as the proportion of simulated cases that show a SSD value distinctive from the original. Calculations were performed in ARLEQUIN, using 1,000 bootstrap to evaluate significance. To estimate the shape of population growth through time for the individuals distributed along the ice sheet coverage area during the LGM, we constructed Bayesian skyline plots implemented in BEAST v 1.4.8 . The appropriate model of nucleotide substitution was HKY+I determined using ModelTest ver. 3.06 . Five million iterations were performed, of which the model parameters were sampled every 1000 iterations. Throughout our analysis, we assumed a within-lineage per site mutation rate of 6%Ma. Demographic plots for each analysis were visualized using Tracer v1.0.1 .
We applied the Partition Homogeineity test (10,000 permutation) to assess the congruence of the evolution rates among CR, ND5 and Cytb using PAUP ver. 4.0b8 . The evolutionary relationship among concatenated CR+ND5+Cyt-b haplotypes was investigated by a Median Joining Network using Network ver. 22.214.171.124 , and phylogenetic reconstructions based on Bayesian (BA) methods. Four divergent haplotypes of Lontra felina , were incorporated in the phylogenetic reconstruction, whereas a L. longicaudis haplotype from the Amazon (this manuscript) was used as an outgroup. The substitution model of DNA evolution was selected based on AIC using Modeltest ver. 3.06 . BA was performed by MrBayes ver. 3.1.2  using the general type of the best fit model parameters defined for the data set, in which four independent analyses were run with four chains each, for six million generations and then sampled at intervals of 1,000 generations. The first 25% of sampled trees were discarded to ensure stabilization and the remaining used to compute a consensus tree. The split frequency was below 0.004, confirming that sampling was from the posterior probability distribution. Mean distance between clades and species was calculated using Mega v.3.1  using p-distance.
This work was supported by Universidad Andrés Bello-DI-06-06/R, Rufford Small Grant for Nature Conservation, Earthwatch Institute and FONDECYT 1100139. Vianna was supported by a CONICYT Doctoral Fellowship, CONICYT Thesis Project AT-23070034.. Special thanks to René Monsalves, Attia Zerega, Juan Carlos Marín, Gerardo Porro, Carla Pozzi, Javier Lucotti who helped with sample collection. Samples from southern Patagonia were collected by Servicio Agricola y Ganadero (SAG) after illegal hunting. Florance Tellier, Andrés Parada and Emma Newcombe helped with analysis and English. All Chilean samples were collected according to permits: Subsecretaria de Pesca (686-2006; 1588-2009; 1228-2009).
- Denton GH, Lowell TV, Moreno PI, Andersen BG, Schlüchter C: Geomorphology, stratigraphy, and radiocarbon chronology of Llanquihue drift in the area of the Southern Lake District, Seno de Reloncaví and Isla de Chiloé, Chile. Geographyrafiska annaler. 1999, 81: 167-229. 10.1111/j.0435-3676.1999.00057.x.View ArticleGoogle Scholar
- Clapperton CC: Quaternary geology and geomorphology of South America. 1993, Amsterdam: ElsevierGoogle Scholar
- McCulloch RD, Bentley MJ, Purves RS, Hulton NRJ, Sugden DE, Clapperton CM: Climatic inferences from glacial and paleoecological evidence at the last glacial termination, southern South America. J Quat Science. 2000, 15: 409-417. 10.1002/1099-1417(200005)15:4<409::AID-JQS539>3.0.CO;2-#.View ArticleGoogle Scholar
- Villagrán C, Moreno P, Villa R: Antecedentes palinológicos acerca de la historia cuaternaria de los bosques chilenos. Ecología de los bosques nativos de Chile. Edited by: Armesto JJ, Villagrán C, Arroyo K. 1995, Santiago: Editorial Universitaria, 51-69.Google Scholar
- Hewitt GM: The genetic legacy of the Quaternary ice ages. Nature. 2000, 405: 907-913. 10.1038/35016000.View ArticlePubMedGoogle Scholar
- Provan J, Bennett KD: Phylogeographic insights into cryptic glacial refugia. Trends Ecol Evol. 2008, 23: 564-571. 10.1016/j.tree.2008.06.010.View ArticlePubMedGoogle Scholar
- Palma ER, Rivera-Milla E, Salazar-Bravo J, Torres-Perez F, Pardinas UFJ, Marquet PA, Spotorno AE, Meynard AP, Yates TL: Phylogeography of Oligoryzomys longicaudatus (Rodentia: Sigmodontinae) in temperate South America. J Mamm. 2005, 86: 191-200. 10.1644/1545-1542(2005)086<0191:POOLRS>2.0.CO;2.View ArticleGoogle Scholar
- Belmar-Lucero S, Godoy P, Ferrés M, Vial P, Palma RE: Range expansion of Olygoryzomys longicaudatus (Rodentia, Sigmodontinae) in Patagonian Chile, and first record of Hantavirus in the region. Revista Chilena de Historia Natural. 2009, 82: 265-275. 10.4067/S0716-078X2009000200008.View ArticleGoogle Scholar
- Lessa EP, D'Elia G, Pardiñas UFJ: 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.View ArticlePubMedGoogle Scholar
- Ruzzante DE, Walde SJ, Cussac VE, Dalebout ML, Seibert J, Ortubay S, Habit E: Phylogeography of the Percichthyidae (Pisces) in Patagonia: roles of orogeny, glaciation, and volcanism. Mol Ecol. 2006, 15: 2949-2968. 10.1111/j.1365-294X.2006.03010.x.View ArticlePubMedGoogle Scholar
- Ruzzante DE, Walde SJ, Gosse JC, Cussac VE, Habit E, Zemlak TS, Adams EDM: Climate control on ancestral population dynamics: insight from Patagonian fish phylogeography. Mol Ecol. 2008, 17: 2234-2244. 10.1111/j.1365-294X.2008.03738.x.View ArticlePubMedGoogle Scholar
- Zemlak TS, Habit EM, Walde SJ, Battini MA, Adams EDM, Ruzzante DE: Across the southern Andes on fin: glacial refugia, drainage reversals and a secondary contact zone revealed by the phylogeographical signal of Galaxias platei in Patagonia. Mol Ecol. 2008, 17: 5049-5061. 10.1111/j.1365-294X.2008.03987.x.View ArticlePubMedGoogle Scholar
- Xu JW, Pérez-Losada M, Jara CG, Crandall KA: Pleistocene glaciation leaves deep signature on the freshwater crab Aegla alacalufi in Chilean Patagonia. Mol Ecol. 2009, 18: 904-918. 10.1111/j.1365-294X.2008.04070.x.View ArticlePubMedGoogle Scholar
- Camus PA: Biogeografía marina de Chile continental. Revista Chilena de Historia Natural. 2001, 74: 587-617. 10.4067/S0716-078X2001000300008.View ArticleGoogle Scholar
- Dyer B: Systematic review and biogeography of the freshwater fishes of Chile. Estudios Oceanológicos. 2000, 19: 77-98.Google Scholar
- Castilla JC, Navarrete SA, Lubchenco J: Southeastern Pacific coastal environments: main features, largescale perturbations, and global climate change. Earth system responses to global change. Edited by: Mooney H, Fuentes E, Kronberg B. 1993, USA: Academic Press, 167-188.Google Scholar
- Valdovinos C, Navarrete S, Marquet P: Mollusk species diversity in the Southeastern Pacific: why are the most species towards the pole?. Ecography. 2003, 26: 139-144. 10.1034/j.1600-0587.2003.03349.x.View ArticleGoogle Scholar
- Valdenegro A, Silva N: Caracterización oceanográfica física y química de la zona de canales y fiordos australes de Chile entre el Estrecho de Magallanes y Cabo de Hornos (Cimar 3 Fiordos). Ciencia y Tecnología del Mar. 2003, 26: 19-60.Google Scholar
- Niemeyer H, Cerceda P: Geografía de Chile. Hidrografía. 1984, Santiago: Instituto Geográfico Militar de ChileGoogle Scholar
- Kruuk H: Otters: Ecology, Behaviour and Conservation. 2006, USA: Oxford University PressView ArticleGoogle Scholar
- Medina G: Conservation and status of Lutra provocax in Chile. Pacific Conserv Biol. 1996, 2: 414-419.Google Scholar
- Larivière S: Lontra provocax. Mammalian Species. 1999, 610: 1-4.Google Scholar
- Aued MB, Chehebar C, Porro G, Macdonald DW, Cassini MH: Environmental correlates of the distribution of Southern River Otters Lontra provocax. Oryx. 2003, 37: 413-421. 10.1017/S0030605303000772.View ArticleGoogle Scholar
- Chehébar C: A survey of the southern river otter Lutra provocax Thomas in Nahuel Huapi National Park, Argentina. Biol Conserv. 1985, 32: 299-307.View ArticleGoogle Scholar
- Chehébar CE, Gallur A, Giannico G, Gotteli MD, Yorio P: A survey of the southern river otter Lutra provocax in Lanín, Puelo and Los Alerces National Parks, Argentina, and evaluation of its conservation status. Biol Conserv. 1986, 38: 293-304.View ArticleGoogle Scholar
- Cassini MH, Fasola L, Chéhebar C, Macdonald DW: Scale-dependent analysis of an otter-crustacean system in Argentinean Patagonia. Naturwissenschaften. 2009, 96: 593-599. 10.1007/s00114-009-0512-2.View ArticlePubMedGoogle Scholar
- Sielfeld W: Características del hábitat de Lutra felina (Molina) y L. provocax Thomas (Carnivora, Mustelidae) en Fuego-Patagonia. Investigaciones en Ciencia y Tecnologia, Serie: Ciencias del Mar. 1990, 1: 30-36.Google Scholar
- Sielfeld W: Sobreposición de nicho y patrones de distribución de Lutra felina y L. provocax (Mustelidae: Carnivora) en el medio marino de Sudamérica austral. Anales del Museo de Historia Natural de Valparaíso. 1989, 20: 103-108.Google Scholar
- Medina G: Seasonal variation and changes in the diet of southern river otter in different freshwater habitats in Chile. Acta Theriologica. 1998, 43: 285-292.View ArticleGoogle Scholar
- Sepúlveda MA, Bartheld JL, Meynard C, Benavides M, Astorga C, Parra D, Medina-Vogel G: Landscape features and crustacean prey as predictors of the Southern river otter distribution in Chile. Anim Conserv. 2009, 12: 522-530.View ArticleGoogle Scholar
- Sepúlveda MA, Bartheld JL, Monsalve R, Gomez V, Medina-Vogel G: Habitat use and spatial behaviour of the endangered Southern river otter (Lontra provocax) in riparian habitats of Chile: Conservation implications. Biol Conserv. 2007, 140: 329-338.View ArticleGoogle Scholar
- Parra DA: Determinación de la distribución de crustáceos y clasificación de la estructura del hábitat en la cuenca del río Toltén y su relación con la presencia-ausencia del huillín (Lontra provocax). Undergraduation thesis. 2006, Universidad Iberoamerica de Ciencias y Tecnología, Medicina VeterinariaGoogle Scholar
- Astorga C, Benavides M, Sepúlveda M, Bartheld JL, Medina-Vogel G: Variables de paisaje y su relación con al distribución del huillin en las cuencas del río Toltén y Queule. El Huillín Lontra provocax: Investigaciones sobre una nutria patagónica amenazada de extinción. Edited by: Cassini MH, Sepulveda M. 2006, Buenos Aires: Fundación Organización Profauna, 83-87.Google Scholar
- Premoli AC, Kitzberger T, Veblen TT: Isozyme variation and recent biogeographical history of the long-lived conifer Fitzroya cupressoides. J Biogeography. 2000, 27: 251-260. 10.1046/j.1365-2699.2000.00402.x.View ArticleGoogle Scholar
- Muellner AN, Tremetsberger K, Stuessy T, Baeza CM: Pleistocene refugia and recolonization routes in the southern Andes: insights from Hypochaeris palustris (Astraceae, Lactuceae). Mol Ecol. 2005, 14: 203-212. 10.1111/j.1365-294X.2004.02386.x.View ArticlePubMedGoogle Scholar
- Sielfeld W, Castilla JC: Estado de conservación y conocimiento de las nutrias en Chile. Estudios Oceanológicos. 1999, 18: 69-79.Google Scholar
- Vianna JA, Ayerdi P, Medina-Vogel G, Mangel JC, Zeballos H, Apaza M, Faugeron S: Phylogeography of the marine otter (Lontra felina): historical and contemporary factors determining its distribution. J Heredity. 2010, 101: 676-689. 10.1093/jhered/esq088.View ArticleGoogle Scholar
- Fraser CI, Thiel M, Spencer HG, Waters JM: 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Rowe KC, Heske EJ, Brown PW, Paige KN: Surviving the ice: Northern refugia and postglacial colonization. Proc Natl Acad Sci USA. 2004, 101: 10355-10359. 10.1073/pnas.0401338101.View ArticlePubMedPubMed CentralGoogle Scholar
- He S: A Checklist of the mooses of Chile. Journal of the Hattori Botanical Laboratory. 1998, 85: 103-189.Google Scholar
- Kim I, Phillips CJ, Monjeau JA, Birney EC, Noack K, Pumo DE, Sikes RS, Dole JA: 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.View ArticlePubMedGoogle Scholar
- Allnutt TR, Newton AC, Lara A, Premoli A, Armesto JJ, Vergara R, Gardner M: Genetic variation in Fitzroya cupressoides (alerce), a threatened South American conifer. Mol Ecol. 1999, 8: 975-987. 10.1046/j.1365-294x.1999.00650.x.View ArticlePubMedGoogle Scholar
- Smith MF, Kelt DA, Patton JL: 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.View ArticlePubMedGoogle Scholar
- Pastorino MJ, Galloa LA, Hattemerb HH: Genetic variation in natural populations of Austrocedrus chilensis, a cypress of the Andean-Patagonian Forest. Biochemical Systematics and Ecology. 2004, 32: 993-1008. 10.1016/j.bse.2004.03.002.View ArticleGoogle Scholar
- Himes CMT, Gallardo MH, Kenagy GJ: Historical biogeography and post-glacial recolonization of South American temperate rain forest by the relictual marsupial Dromiciops gliroides. J Biogeography. 2008, 35: 1415-1424. 10.1111/j.1365-2699.2008.01895.x.View ArticleGoogle Scholar
- González-Ittig RE, Rossi-Fraire HJ, Cantoni GE, Herrero ER, Benedetti R, Gallardo MH, Gardenal CN: Population genetic structure of long-tailed pygmy rice rats (Oligoryzomys longicaudatus) from Argentina and Chile based on the mitochondrial control region. Canadian J Zoology. 2010, 88: 23-35.View ArticleGoogle Scholar
- Unmack PJ, Bennin AP, Habit EM, Victoriano PF, Johnson JB: Impact of ocean barriers, topography, and glaciation on the phylogeography of the catfish Trichomycterus areolatus (Teleostei: Trichomycteridae) in Chile. Biol J Linn Soc. 2009, 97: 876-892. 10.1111/j.1095-8312.2009.01224.x.View ArticleGoogle Scholar
- Armesto JJ, Rozzi R, Smith-Ramírez C, Arroyo MTK: Conservation Targets in South American Temperate Forests. Science. 1995, 282: 1271-1272. 10.1126/science.282.5392.1271.View ArticleGoogle Scholar
- Zemlak TS, Habit EM, Walde SJ, Carrera C, Ruzzante DE: Surviving historical Patagonian landscapes and climate: molecular insights from Galaxias maculates. BMC Evol Biol. 2010, 10: 67-10.1186/1471-2148-10-67.View ArticlePubMedPubMed CentralGoogle Scholar
- Pugh PJA, Davenport J: Colonisation vs. disturbance: The effects of sustained ice-scouring on intertidal communities. Journal of Experimental Marine Biology and Ecology. 1997, 210: 1-21. 10.1016/S0022-0981(96)02711-6.View ArticleGoogle Scholar
- Blundell GM, Ben-David M, Groves P, Bowyer RT, Geffen E: Characteristics of sex-biased dispersal and gene flow in coastal river otters: implications for natural recolonization of extirpated populations. Mol Ecol. 2002, 11: 289-303. 10.1046/j.0962-1083.2001.01440.x.View ArticlePubMedGoogle Scholar
- Cassens I, Tiedemann R, Suchentrunk F: Mitochondrial DNA variation in the European otter (Lutra lutra) and the use of Spatial Autocorrelation Analysis in conservation. J Heredity. 2000, 91: 31-35. 10.1093/jhered/91.1.31.View ArticleGoogle Scholar
- Sommer R, Benecke N: Late- and Post-Glacial history of the Mustelidae in Europe. Mammal Review. 2004, 34: 249-284. 10.1111/j.1365-2907.2004.00043.x.View ArticleGoogle Scholar
- Klopfstein S, Currat M, Excoffier L: The Fate of Mutations Surfing on the Wave of a Range Expansion. Mol Biol Evol. 2006, 23: 482-490. 10.1093/molbev/msj057.View ArticlePubMedGoogle Scholar
- Excoffier L, Ray N: Surfing during population expansions promotes genetic revolutions and structuration. Trends Ecol Evol. 2008, 23: 347-351. 10.1016/j.tree.2008.04.004.View ArticlePubMedGoogle Scholar
- Nei M, Maruyama T, Chakraborty R: The bottleneck effect and genetic variability in populations. Evolution. 1975, 29: 1-10. 10.2307/2407137.View ArticleGoogle Scholar
- Medina-Vogel G, Kaufman VS, Monsalve R, Gomez V: The influence of riparian vegetation, woody debris, stream morphology and human activity on the use of rivers by southern river otters in Lontra provocax in Chile. Oryx. 2003, 37: 422-430. 10.1017/S0030605303000784.View ArticleGoogle Scholar
- Medina-Vogel G, González-Lagos C: Habitat use and diet of endangered southern river otter Lontra provocax in a predominantly palustrine wetland in Chile. Wildlife Biology. 2008, 14: 211-220. 10.2981/0909-6396(2008)14[211:HUADOE]2.0.CO;2.View ArticleGoogle Scholar
- República de Chile: Clasificación de especies según conservación. Decreto Supremo 151/06. 2007, Santiago: Ministerio Secretaría General de la PresidenciaGoogle Scholar
- Smith-Ramírez C: The Chilean coastal range: a vanishing center of biodiversity and endemism in South American temperate rainforests. Biodiversity and Conservation. 2004, 13: 373-393.View ArticleGoogle Scholar
- Filatov DA: ProSeq: A software for preparation and evolutionary analysis of DNA sequence data sets. Mol Ecol Notes. 2002, 2: 621-624. 10.1046/j.1471-8286.2002.00313.x.View ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research. 1997, 24: 4876-4882. 10.1093/nar/25.24.4876.View ArticleGoogle Scholar
- Dupanloup I, Schneider S, Excoffier L: A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002, 11: 2571-2581. 10.1046/j.1365-294X.2002.01650.x.View ArticlePubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005, 1: 47-50.Google Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMedPubMed CentralGoogle Scholar
- Fu YX: Statistical Tests of Neutrality of Mutations against Population Growth, Hitchhiking and Background Selection. Genetics. 1997, 147: 915-925.PubMedPubMed CentralGoogle Scholar
- Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9: 552-569.PubMedGoogle Scholar
- Excoffier L: Patterns of DNA sequence diversity and genetic structure after a range expansion: lessons from the infinite-island model. Mol Ecol. 2004, 13: 853-864. 10.1046/j.1365-294X.2003.02004.x.View ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian Evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 213-10.1186/1471-2148-7-214.View ArticleGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics Application Note. 1998, 14: 817-818.View ArticleGoogle Scholar
- Swofford DL: PAUP*, Phylogenetic Analysis Using Parsimony (*and Other Methods). 2000, Sunderland: Sinauer AssociatesGoogle Scholar
- Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.View ArticlePubMedGoogle Scholar
- Huelsenbeck JL, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA 3.1: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Briefings in Bioinformatics. 2004, 5: 150-163. 10.1093/bib/5.2.150.View ArticlePubMedGoogle Scholar
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.