Contrasting Genetic Structure and Diversity of Galaxias maculatus (Jenyns, 1848) Along the Chilean Coast: Stock Identification for Fishery Management

Cryptic diversity within two widespread diadromous freshwater fishes (Teleostei: Galaxiidae)

Identification of taxonomically cryptic species is essential for the effective conservation of biodiversity. Freshwater-limited organisms tend to be genetically isolated by drainage boundaries, and thus may be expected to show substantial cryptic phylogenetic and taxonomic diversity. By comparison, populations of diadromous taxa, that migrate between freshwater and marine environments, are expected to show less genetic differentiation. Here we test for cryptic diversity in Australasian populations (both diadromous and non-diadromous) of two widespread Southern Hemisphere fish species, Galaxias brevipinnis and Galaxias maculatus. Both mtDNA and nuclear markers reveal putative cryptic species within these taxa. The substantial diversity detected within G. brevipinnis may be explained by its strong climbing ability which allows it to form isolated inland populations. In island populations, G. brevipinnis similarly show deeper genetic divergence than those of G. maculatus, which may be explained by the greater abundance of G. maculatus larvae in the sea allowing more ongoing dispersal. Our study highlights that even widespread, 'high-dispersal' species can harbour substantial cryptic diversity and therefore warrant increased taxonomic and conservation attention.

Species tree analyses and speciation-based species delimitation support new species in the relict catfish family Diplomystidae and provide insights on recent glacial history in Patagonia

Diplomystidae is an early-diverged family of freshwater catfish endemic to southern South America. We have recently collected five juvenile specimens belonging to this family from the Bueno River Basin, a basin which the only previous record was a single juvenile specimen collected in 1996. This finding confirms the distribution of the family further South in northern Patagonia, but poses new questions about the origin of this population in an area with a strong glacial history. We used phylogenetic analyses to evaluate three different hypotheses that could explain the origin of this population in the basin. First, the population could have originated in Atlantic basins (East of the Andes) and dispersed to the Bueno Basin after the Last Glacial Maximum (LGM) via river reversals, as it has been proposed for other population of Diplomystes as well as for other freshwater species from Patagonia. Second, the population could have originated in the geographically close Valdivia Basin (West of the Andes) and dispersed south to its current location in the Bueno Basin. Third, regardless of its geographic origin (West or East of the Andes), the Bueno Basin population could have a longer history in the basin, surviving in situ through the LGM. In addition, we conducted species delimitation analyses using a recently developed method that uses a protracted model of speciation. Our goal was to test the species status of the Bueno Basin population along with another controversial population in Central Chile (Biobío Basin), which appeared highly divergent in previous studies with mtDNA. The phylogenetic analyses showed that the population from the Bueno Basin is more related to Atlantic than to Pacific lineages, although with a deep divergence that predated the LGM, supporting in situ survival rather than postglacial dispersal. In addition, these analyses also showed that the species D. nahuelbutaensis is polyphyletic, supporting the need for a taxonomic reevaluation. The species delimitation analyses supported two new species which are described using molecular diagnostic characters: Diplomystes arratiae sp. nov. from the Biobío, Carampangue, and Laraquete basins, maintaining D. nahuelbutaensis valid only for the Imperial Basin, and Diplomystes habitae sp. nov. from the Bueno Basin. This study greatly increases the number of species within both the family Diplomystidae and Patagonia, and contributes substantially to the knowledge of the evolution of southern South American freshwater biodiversity during its glacial history. Given the important contribution to the phylogenetic diversity of the family, we recommend a high conservation priority for both new species. Finally, this study highlights an exemplary scenario where species descriptions based only on DNA data are particularly valuable, bringing additional elements to the ongoing debate on DNA-based taxonomy.

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.

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.

Quaternary ice sheets and sea level regression drove divergence in a marine gastropod along Eastern and Western coasts of South America

The southern coastline of South America is a remarkable area to evaluate how Quaternary glacial processes impacted the demography of the near-shore marine biota. Here we present new phylogeographic analyses in the pulmonate Siphonaria lessonii across its distribution, from northern Chile in the Pacific to Uruguay in the Atlantic. Contrary to our expectations, populations from the southwestern Atlantic, an area that was less impacted by ice during glacial maxima, showed low genetic diversity and evidence of recent expansion, similar to the patterns recorded in this study across heavily ice-impacted areas in the Pacific Magellan margin. We propose that Atlantic and Pacific shallow marine hard-substrate benthic species were both affected during the Quaternary in South America, but by different processes. At higher latitudes of the southeast Pacific, ice-scouring drastically affected S. lessonii populations compared to non-glaciated areas along the Chile-Peru province where the species was resilient. In the southwest Atlantic, S. lessonii populations would have been dramatically impacted by the reduction of near-shore rocky habitat availability as a consequence of glacio-eustatic movements. The increase of gravelly and rocky shore substrates in the southwest Atlantic supports a hypothesis of glacial refugia from where the species recolonized lower latitudes across the Atlantic and Pacific margins. Our results suggest that current patterns of genetic diversity and structure in near-shore marine benthic species do not solely depend on the impact of Quaternary glacial ice expansions but also on the availability of suitable habitats and life-history traits, including developmental mode, bathymetry and the likelihood of dispersal by rafting.

Otolith microchemistry and diadromy in Patagonian river fishes

Coastal habitats in Chile are hypothesized to support a number of diadromous fish species. The objective of this study was to document migratory life histories of native galaxiids and introduced salmonids from a wide latitudinal range in Chilean Patagonia (39–48°S). Otolith microchemistry data were analysed using a recursive partitioning approach to test for diadromy. Based on annular analysis of Sr:Ca ratios, a diadromous life history was suggested for populations of native Aplochiton taeniatus, A. marinus, and Galaxias maculatus. Lifetime residency in freshwater was suggested for populations of A. zebra and G. platei. Among introduced salmonids, populations of Oncorhynchus tshawytscha and O. kisutch exhibited patterns consistent with anadromy, whereas the screened population of O. mykiss appeared restricted to freshwater. Salmo trutta exhibited variable patterns suggesting freshwater residency and possibly anadromy in one case. The capacity and geographic scope of hydropower development is increasing and may disrupt migratory routes of diadromous fishes. Identification of diadromous species is a critical first step for preventing their loss due to hydropower development.

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.

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

Major geologic and climatic changes during the Quaternary exerted a major role in shaping past and contemporary distribution of genetic diversity and structure of aquatic organisms in southern South America. In fact, the northern glacial limit along the Pacific coast, an area of major environmental changes in terms of topography, currents, and water salinity, represents a major biogeographic transition for marine and freshwater species. We used mitochondrial DNA sequences (D-loop) to investigate the consequences of Quaternary glacial cycles over the pattern of genetic diversity and structure of G. maculatus (Pisces: Galaxiidae) along two biogeographical provinces in the Chilean coast. Extreme levels of genetic diversity and strong phylogeographic structure characterize the species suggesting a low amount of influence of the last glacial cycle over its demography. However, we recognized contrasting patterns of genetic diversity and structure between main biogeographical areas here analyzed. Along the Intermediate Area (38°-41° S) each estuarine population constitutes a different unit. In contrast, Magellanic populations (43°-53° S) exhibited low levels of genetic differentiation. Contrasting patterns of genetic diversity and structure recorded in the species between the analyzed biogeographic areas are consistent with the marked differences in abiotic factors (i.e., different coastal configurations, Quaternary glacial histories, and oceanographic regimes) and to inherent characteristics of the species (i.e., salt-tolerance, physiology, and reproductive behavior).