Phylogeography in Galaxias maculatus (Jenyns, 1848) along Two Biogeographical Provinces in the Chilean Coast

Contrasting Patterns of Genetic Diversity and Divergence Between Landlocked and Migratory Populations of Fish Galaxias maculatus, Evaluated Through Mitochondrial DNA Sequencing and Nuclear DNA Microsatellites

Galaxias species are interesting biogeographic models due to their distribution and different types of life cycles, with migratory and landlocked populations. To obtain a better understanding of the genetic consequences of the Quaternary glacial cycles in Galaxias maculatus, in this work we compared landlocked and migratory populations collected in areas that were differentially affected by ice advances and retreats. We included nine populations of G. maculatus, four collected from lakes (landlocked) and five from their associated estuaries/rivers (migratory) in three estuary-lake systems across southern Chile. Genetic analyses were performed using the mitochondrial control region and nine microsatellite loci. Genetic diversity measured with both markers was significantly higher in migratory than in landlocked populations across the study area. The levels of genetic differentiation showed higher differentiation among lakes than estuaries. Genetic diversity was higher in migratory populations located in areas that were less impacted by ice during Quaternary glacial processes. These results may be the consequence of recent recolonization of small freshwater bodies following the Last Glacial Maximum (LGM). Finally, the greatest differentiation was observed in populations that were exposed to continental ice advances and retreats during the LGM. Thus, in the present work we corroborate a pattern of differentiation between lakes and estuaries, using mtDNA sequences and microsatellite nuclear markers. This pattern may be due to a combination of biological factors, i.e., resident non-migratory behaviour or landlocking and natal homing-in, as well as geological factors, i.e., Expansion-Contraction Quaternary glacial biogeographic processes.

The right tool for the right question: contrasting biogeographic patterns in the notothenioid fish Harpagifer spp. along the Magellan Province

Molecular-based analysis has become a fundamental tool to understand the role of Quaternary glacial episodes. In the Magellan Province in southern South America, ice covering during the last glacial maximum (20 ka) radically altered the landscape/seascape, speciation rates and distribution of species. For the notothenioid fishes of the genus Harpagifer, in the area are described two nominal species. Nevertheless, this genus recently colonized South America from Antarctica, providing a short time for speciation processes. Combining DNA sequences and genotyping-by-sequencing SNPs, we evaluated the role of Quaternary glaciations over the patterns of genetic structure in Harpagifer across its distribution in the Magellan Province. DNA sequences showed low phylogeographic structure, with shared and dominant haplotypes between nominal species, suggesting a single evolutionary unit. SNPs identified contrastingly two groups in Patagonia and a third well-differentiated group in the Falkland/Malvinas Islands with limited and asymmetric gene flow. Linking the information of different markers allowed us to infer the relevance of postglacial colonization mediated by the general oceanographic circulation patterns. Contrasting rough- and fine-scale genetic patterns highlights the relevance of combined methodologies for species delimitation, which, depending on the question to be addressed, allows discrimination among phylogeographic structure, discarding incipient speciation, and contemporary spatial differentiation processes.

Cenozoic climatic changes drive evolution and dispersal of coastal benthic foraminifera in the Southern Ocean

The Antarctic coastal fauna is characterized by high endemism related to the progressive cooling of Antarctic waters and their isolation by the Antarctic Circumpolar Current. The origin of the Antarctic coastal fauna could involve either colonization from adjoining deep-sea areas or migration through the Drake Passage from sub-Antarctic areas. Here, we tested these hypotheses by comparing the morphology and genetics of benthic foraminifera collected from Antarctica, sub-Antarctic coastal settings in South Georgia, the Falkland Islands and Patagonian fjords. We analyzed four genera (Cassidulina, Globocassidulina, Cassidulinoides, Ehrenbergina) of the family Cassidulinidae that are represented by at least nine species in our samples. Focusing on the genera Globocassidulina and Cassidulinoides, our results showed that the first split between sub-Antarctic and Antarctic lineages took place during the mid-Miocene climate reorganization, probably about 20 to 17 million years ago (Ma). It was followed by a divergence between Antarctic species ~ 10 Ma, probably related to the cooling of deep water and vertical structuring of the water-column, as well as broadening and deepening of the continental shelf. The gene flow across the Drake Passage, as well as between South America and South Georgia, seems to have occurred from the Late Miocene to the Early Pliocene. It appears that climate warming during 7-5 Ma and the migration of the Polar Front breached biogeographic barriers and facilitated inter-species hybridization. The latest radiation coincided with glacial intensification (~ 2 Ma), which accelerated geographic fragmentation of populations, demographic changes, and genetic diversification in Antarctic species. Our results show that the evolution of Antarctic and sub-Antarctic coastal benthic foraminifera was linked to the tectonic and climatic history of the area, but their evolutionary response was not uniform and reflected species-specific ecological adaptations that influenced the dispersal patterns and biogeography of each species in different ways.

Phylogeography of high Andean killifishes Orestias (Teleostei: Cyprinodontidae) in Caquena and Lauca sub-basins of the Altiplano (Chile): mitochondrial and nuclear analysis of an endangered fish

From the early Miocene, the uplift of the Andes Mountains, intense volcanic activity and the occurrence of successive periods of dryness and humidity would have differentially influenced the modification of Altiplano watersheds, and consequently the evolutionary history of the taxa that live there. We analyzed Orestias populations from the Caquena and Lauca Altiplanic sub-basins of northern Chile to determine their genetic differentiation and relationship to their geographical distribution using mitochondrial (D-loop) and nuclear (microsatellite) molecular markers and to reconstruct its biogeographic history on these sub-basins. The results allowed reconstructing and reevaluating the evolutionary history of the genus in the area; genic diversity and differentiation together with different founding genetic groups suggest that Orestias have been spread homogeneously in the study area and would have experienced local disturbances that promoted isolation and diversification in restricted zones of their distribution.

Assessment of Heavy Metal Contamination in Two Edible Fish Species and Water from North Patagonia Estuary

Estuaries worldwide have been severely degraded and become reservoirs for many types of pollutants, such as heavy metals. This study investigated the levels of Cd, Cu, Mn, Ni, Pb, and Zn in water and whole fish. We sampled 40 juvenile silversides Odontesthes regia and 41 juvenile puye Galaxias maculatus from the Valdivia River estuary, adjacent to the urban area in southern South America (Chile). Samples were analyzed using a flame atomic absorption spectrophotometer. In water samples, metals except Zn were mostly below the detection limits and all metals were below the maximum levels established by local guidelines in this estuary. In whole fish samples, concentrations of Cu, Zn, Pb, Mn, and Cd were significantly higher in puyes than in silversides. Additionally, Zn, Pb, and Mn were correlated to body length and weight in puyes, whereas Cd was correlated to body length in silversides. The mean concentration of heavy metals in silverside and puyes were higher than those reported in the literature. In silversides, all heavy metal levels were below the limits permitted by current legislation (FAO), whereas in puyes Pb and Cd levels were above the recommended maximum level established by international guidelines, therefore putting the human population at risk.

Geographic variation in composition of metazoan parasite infracommunities in Galaxias maculatus Jenyns 1842 (Osmeriformes: Galaxiidae) in southern Chile (38-47° S)

Galaxias maculatus is an abundant freshwater fish species in Chilean continental waters where it plays important ecological functions, yet few parasitological records of this species exist in Chile and all of them cover a very limited geographic range. The objective of this study was to assess large scale geographic variation in composition of parasite infracommunities of Galaxias maculatus. Specifically, parasite infracommunities of this species were compared among 11 locations across 9 degrees of latitude and 3 ecosystem types (lake, river and estuary). Most taxa found had been previously reported in Chile and Argentina. However, this is the first report for Tylodelphys sp. in this host in Chile. Furthermore, the cranial parasite Tylodelphys sp. had the highest overall prevalence and abundance compared to other parasite species. Despite the fact that the abundance of Tylodelphys sp. was not significantly correlated with Fulton’s condition factor of fish, infected fish seem to have a better body condition compared to uninfected ones. The most important source of variation in composition of infracommunities was the sampling location. Furthermore, fish from lakes have a different composition of parasite infracommunities mainly due to higher abundances of Tylodelphys sp.

Historical and Contemporary Diversity of Galaxiids in South America: Biogeographic and Phylogenetic Perspectives

Galaxiid fishes from South America are represented by three genera (Aplochiton, Brachygalaxias and Galaxias) and eight species. Their genetic patterns have been studied over the last two decades to disentangle how historical and contemporary processes influenced their biogeographic distribution and phylogeographic patterns. Here we review and synthesize this body of work. Phylogeographic approaches reveal the important role played by orogeny and the expansion/melting of glacial ice during the Quaternary. Populations retreated to glacial refugia during glacial times and some systems experienced drainage reversals from the Atlantic to the Pacific following deglaciation. Although most species expanded their populations and increased their genetic diversity during the Holocene, the introduction of salmonids and the construction of dams are likely to lead to a decline in genetic diversity for at least some species. An improvement in our understanding of the processes that influenced historical and contemporary diversity patterns among galaxiid and other native fishes in South America is necessary for addressing the cumulative and synergistic impacts of human activity on this unique freshwater fauna.

Contrasting evolutionary responses in two co-distributed species of Galaxias (Pisces, Galaxiidae) in a river from the glaciated range in Southern Chile

Life-history traits are among the most important factors affecting population abundance and genetic diversity of species. Here, we analysed the genetic patterns of two Galaxias species with different life-history traits to investigate how these biological differences impacted their evolution in the Valdivia River basin, Southern Chile. We analysed mitochondrial DNA (mtDNA) sequences from 225 individuals of Galaxias maculatus and 136 of G. platei to compare patterns of genetic diversity, structure and demographic growth across the basin. Galaxias maculatus presented higher genetic diversity and higher genetic structure than G. platei. Demographic analyses showed G. maculatus kept a higher population size over time, with a signal of demographic expansion in the last 250 kyr. Whereas Galaxias platei, exhibited lower, but constant population size over time. Furthermore, haplotype networks revealed higher lineage diversity in G. maculatus with a tendency to occupy different areas of the basin. Coalescent simulations ruled out that genetic differences between species could be explained by stochastic processes (genetic drift), suggesting species-specific biological differences as responsible for the observed genetic differences. We discuss how differences in life-history traits and past glaciations interact to shape the evolutionary history of the two Galaxias species.

Multiple drainage reversal episodes and glacial refugia in a Patagonian fish revealed by sequenced microsatellites

The rise of the southern Andes and the Quaternary glacial cycles influenced the landscape of Patagonia, affecting the phylogeographic and biogeographic patterns of its flora and fauna. Here, we examine the phylogeography of the freshwater fish, Percichthys trucha, using 53 sequenced microsatellite DNA markers. Fish (n= 835) were collected from 16 river systems (46 locations) spanning the species range on both sides of the Andes. Eleven watersheds drain to the Pacific, five of which are trans-Andean (headwaters east of Andes). The remaining five drainages empty into the Atlantic. Three analytical approaches (neighbour-joining tree, hierarchical AMOVAs, Structure) revealed evidence of historic drainage reversals: fish from four of the five trans-Andean systems (Puelo, Futalaufquen/Yelcho, Baker, Pascua) exhibited greater genetic similarity with Atlantic draining systems than with Pacific systems with headwaters west of Andes. Present-day drainage (Pacific versus Atlantic) explained only 5% of total genetic variance, while ancestral drainage explained nearly 27% of total variance. Thus, the phylogeographic structure of P. trucha is consistent with episodes of drainage reversal in multiple systems and suggests a major role for deglaciation in the genetic and indeed the geographical distribution of P. trucha in Patagonia. The study emphasizes the significant role of historical processes in the current pattern of genetic diversity and differentiation in a fish from a southern temperate region.

Phenotypic and genetic differentiation between diadromous and landlocked puyen Galaxias maculatus

This study reports the phenotypic and genetic differences between individuals of puyen Galaxias maculatus from two sites in the same river basin in Tierra del Fuego National Park, southern South America. Individuals from the two sampling sites presented morphometric and genetic differences. The morphometric differences indicated that individuals from Laguna Negra (LN) were short and more robust and had large eyes, whereas those from Arroyo Negro (AN) were thin and elongated and had small eyes. Genetic differences showed that AN individuals had a greater genetic structuration and an older demographic history than LN individuals. The results of this study affirmed that the individuals from the two sampling sites belong to different populations with a high degree of isolation. The demographic history could indicate that the individuals of G. maculatus which migrated to northern areas during the last glaciation settled in the Beagle Channel after its formation. The LN population could have originated after the retreat of the glaciers, migrating from AN.

HaCeD-Seq: a Novel Method for Reliable and Easy Estimation About the Fish Population Using Haplotype Count from eDNA

It is common to count the numbers of specified fish in the field after speciation of captured fish according to their morphology and to subject these counts to appropriate statistical analyses. In recent years, a non-invasive method to estimate the abundance of a particular fish species using environmental DNA (eDNA) has been developed. However, it is still difficult to determine accurate numbers of fish species using such method. We predict that the estimation of individuals of certain fish species in the field is more accurate and easier by using haplotypes of DNA in the fast evolutionary region. Therefore, we focused on the regulatory region (D-loop) in mitochondrial DNA, which is known to have a high genetic variation at the intraspecific level of the targeting eel. We investigated haplotype diversity in eel at first and then determined the number of D-loop haplotypes contained in their exfoliated cells in breeding water. Finally, we developed a novel analytical method, HaCeD-Seq, to estimate the number of individuals based on the abovementioned data.

Local persistence of Mann's soft-haired mouse Abrothrix manni (Rodentia, Sigmodontinae) during Quaternary glaciations in southern Chile

Quaternary climatic oscillations have impacted Patagonian sigmodontine fauna, leaving traceable genetic footprints. In southern Chile, changes in the landscape included transitions to different vegetation formations as well as the extension of ice sheets. In this study, we focus on the Valdivian forest endemic and recently described sigmodontine species Abrothrix manni. We aim to assess the genetic structure of this species, testing for the existence of intraspecific lineages, and inferring the recent demographic history of the species. Analyses were based on the first 801 bp of the mitochondrial gene Cytocrhome-b from 49 individuals of A. manni collected at 10 localities that covers most part of its geographic distribution. Genealogical analyses recovered two main intraspecific lineages that are geographically segregated and present an intermediate site of secondary contact. Historical demography shows signal of recent population decrease. Based on these results, we proposed that current genetic diversity of A. manni differentiated in at least two distinct refugial areas in southern Chile. This scenario, in addition to be unique among those uncovered for the so far studied Valdivian forest rodents, is noteworthy because of the reduced geographic scale inhabited by the species.

Historical dimensions of population structure in a continuously distributed marine species: The case of the endemic Chilean dolphin

The complementarity of historical and contemporary processes contributes to understanding the genetic structure of continuously distributed marine species with high dispersal capabilities. Cephalorhynchus eutropia, has a continuous coastal distribution with strong genetic differentiation identified by nuclear DNA markers. We explored the historical dimension of this genetic differentiation between northern and southern populations to evaluate phylogeographic structure. Additionally, we conducted mtDNA and microsatellite analyses to detect past and recent demographic changes. The southern population was characterized by lower genetic diversity with a signal of population expansion, likely associated with ice retreat and habitat extension after the Last Glacial Maximum (LGM). In contrast, structure within the northern population was more consistent with stable historical population size. Approximate Bayesian Computation analyses suggested that during the LGM, C. eutropia persisted in the northern area; while the south was colonized by dispersal ~11,000 years ago followed by population expansion. This study shows that Chilean dolphin population structure is consistent with predictions from the Expansion-Contraction biogeographic model, with a poleward post-glacial shift revealed in current genetic structure. The results also confirm the validity of the population units previously identified, demonstrating their historical origin and highlighting the utility of integrating genetic markers with different temporal scale resolutions.

Genetics, Gene Flow, and Glaciation: The Case of the South American Limpet Nacella mytilina

Glacial episodes of the Quaternary, and particularly the Last Glacial Maximum (LGM) drastically altered the distribution of the Southern-Hemisphere biota, principally at higher latitudes. The irregular coastline of Patagonia expanding for more than 84.000 km constitutes a remarkable area to evaluate the effect of Quaternary landscape and seascape shifts over the demography of near-shore marine benthic organisms. Few studies describing the biogeographic responses of marine species to the LGM have been conducted in Patagonia, but existing data from coastal marine species have demonstrated marked genetic signatures of post-LGM recolonization and expansion. The kelp-dweller limpet Nacella mytilina is broadly distributed along the southern tip of South America and at the Falkland/Malvinas Islands. Considering its distribution, abundance, and narrow bathymetry, N. mytilina represents an appropriate model to infer how historical and contemporary processes affected the distribution of intraspecific genetic diversity and structure along the southern tip of South America. At the same time, it will be possible to determine how life history traits and the ecology of the species are responsible for the current pattern of gene flow and connectivity across the study area. We conducted phylogeographic and demographic inference analyses in N. mytilina from 12 localities along Pacific Patagonia (PP) and one population from the Falkland/Malvinas Islands (FI). Analyses of the mitochondrial gene COI in 300 individuals of N. mytilina revealed low levels of genetic polymorphism and the absence of genetic differentiation along PP. In contrast, FI showed a strong and significant differentiation from Pacific Patagonian populations. Higher levels of genetic diversity were also recorded in the FI population, together with a more expanded genealogy supporting the hypothesis of glacial persistence of the species in these islands. Haplotype genealogy, and mismatch analyses in the FI population recognized an older and more complex demographic history than in PP. Demographic reconstructions along PP suggest a post-LGM expansion process (7.5 ka), also supported by neutrality tests, mismatch distribution and maximum parsimony haplotype genealogies. Migration rate estimations showed evidence of asymmetrical gene flow from PP to FI. The absence of genetic differentiation, the presence of a single dominant haplotype, high estimated migration rates, and marked signal of recent demographic growth, support the hypothesis of rapid post-glacial expansion in N. mytilina along PP. This expansion could have been sustained by larval and rafting-mediated dispersal of adults from northernmost populations following the Cape Horn Current System. Marked genetic differentiation between PP and FI could be explained through differences in their respective glacial histories. During the LGM, Pacific Patagonia (PP) was almost fully covered by the Patagonian Ice Sheet, while sheet coverage in the FI ice was restricted to small cirques and valleys. As previously recorded in the sister-species N. magellanica, the FI rather than represent a classical glacial refugium for N. mytilina, seems to represent a sink area and/or a secondary contact zone. Accordingly, historical and contemporary processes, contrasting glacial histories between the analyzed sectors, as well as life history traits constitute the main factors explaining the current biogeographical patterns of most shallow Patagonian marine benthic organisms.