Echoes of a distant time: effects of historical processes on contemporary genetic patterns in Galaxias platei in Patagonia

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.

Genetic Diversity and Population Differentiation of Chinese Lizard Gudgeon (Saurogobio dabryi) in the Upper Yangtze River

Dam construction on the upper Yangtze River has dramatically altered riverine ecosystems and caused habitat fragmentation of fishes, which might influence the genetic structure of fish populations. In this study, we examined the possible genetic effects of dam construction on Chinese lizard gudgeon (Saurogobio dabryi) populations in the upper Yangtze River, China. Seven populations were sampled, and genetic structure was analyzed using single nucleotide polymorphism (SNP) markers through the specific locus amplified fragment sequencing (SLAF-seq) method. The numbers of SNPs were lower in the upstream populations than in the downstream populations. Genetic similarity was increased from downstream to upstream. The upstream populations of S. dabryi might be more vulnerable to genetic drift than those downstream. Structure analysis indicated three distinct genetic groups of S. dabryi in the upper Yangtze River, among which the genetic differentiation values (Fst) were at a high level. The genetic differentiation of S. dabryi exhibited a close correlation with spatial distance. We did not detect a significant correlation between isolation time and genetic differentiation, suggesting that impacts of dams on the genetic structure of S. dabryi can be relatively minimal on a short time scale. The results quantify the genetic diversity and population structure patterns of S. dabryi after habitat fragmentation caused by dams, which will provide a reference for resource protection and management of this species in the upper Yangtze River.

Pleistocene expansion, anthropogenic pressure and ocean currents: Disentangling the past and ongoing evolutionary history of Patella aspera Röding, 1798 in the archipelago of Madeira

AIMS

Rising sea-level following the Last Glacial Maximum lead to fragmentation of coastal limpet populations between islands of the Archipelago of Madeira. This fragmentation is reinforced by recent heavy exploitation reducing effective population size on Madeira Island. We use the limpet P. aspera to understand how the role of processes at different time scales (i.e. changes in the sea level and overexploitation) can influence the genetic composition of an extant species, relating these processes to reproductive phenology and seasonal shifts in ocean currents.

LOCATION

Madeira Island, Porto Santo and Desertas (Archipelago of Madeira, NE Atlantic Ocean).

TAXON

The limpet Patella aspera.

METHODS

Twelve microsatellite genetic markers were used. A power analysis was used to evaluate the power of the microsatellite markers to detect a signal of population differentiation. Long-term past migrations were assessed using a Bayesian Markov Montecarlo approach in the software MIGRATE-n to estimate mutation-scaled migration rates (M = m/μ; m, probability of a lineage immigrating per generation; μ, mutation rate). Two scenarios were evaluated using an Approximate Bayesian Computation (ABC) in the software DIYABC 2.1 (i) Scenario 1: considered a population scenario from a reduced Ne at time t₃ to a higher Ne at time t₂; and (ii) Scenario 2 considering a reduction of Ne from a time t₃ to a time t₂.

RESULTS

Colonization of the archipelago by Portuguese settlers six centuries ago probably led to an important decrease in the genetic diversity of the species (N_e). Contemporary gene flow strongly support a pattern of high asymmetric connectivity explained by the reproductive phenology of the species and spatio-temporal seasonal changes in the ocean currents. Spatio-temporal reconstructions using Bayesian methods, including coalescent and Approximate Bayesian Computation (ABC) approaches, suggest changes in the migration patterns from highly symmetric to highly asymmetric connectivity with subtle population differentiation as consequence of post-glacial maximum sea level rise during the Holocene.

MAIN CONCLUSIONS

Our results suggest that anthropogenic activity could have had serious effects on the genetic diversity of heavily exploited littoral species since the end of the Pleistocene, probably accelerating in recent years.

Intact landscape promotes gene flow and low genetic structuring in the threatened Eastern Massasauga Rattlesnake

Genetic structuring of wild populations is dependent on environmental, ecological, and life-history factors. The specific role environmental context plays in genetic structuring is important to conservation practitioners working with rare species across areas with varying degrees of fragmentation. We investigated fine-scale genetic patterns of the federally threatened Eastern Massasauga Rattlesnake (Sistrurus catenatus) on a relatively undisturbed island in northern Michigan, USA. This species often persists in habitat islands throughout much of its distribution due to extensive habitat loss and distance-limited dispersal. We found that the entire island population exhibited weak genetic structuring with spatially segregated variation in effective migration and genetic diversity. The low level of genetic structuring contrasts with previous studies in the southern part of the species' range at comparable fine scales (~7 km), in which much higher levels of structuring were documented. The island population's genetic structuring more closely resembles that of populations from Ontario, Canada, that occupy similarly intact habitats. Intrapopulation variation in effective migration and genetic diversity likely corresponds to the presence of large inland lakes acting as barriers and more human activity in the southern portion of the island. The observed genetic structuring in this intact landscape suggests that the Eastern Massasauga is capable of sufficient interpatch movements to reduce overall genetic structuring and colonize new habitats. Landscape mosaics with multiple habitat patches and localized barriers (e.g., large water bodies or roads) will promote gene flow and natural colonization for this declining species.

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.

New Insights into the Distribution, Physiology and Life Histories of South American Galaxiid Fishes, and Potential Threats to This Unique Fauna

South American galaxiids occupy both Patagonia and the ichthyogeographic Chilean Province, encompassing glacial Andean deep lakes, shallow plateau lakes, reservoirs, short Pacific rivers and long Atlantic rivers. The total fish fauna includes 29 species, comprising Neotropical fishes (siluriforms and characids), galaxiids, percichthyids, atherinopsids and mugilids, two lamprey species, and several exotic fishes (salmonids, Gambusia spp. and common carp). The family Galaxiidae shares a common ancestry with the Gondwanan temperate fish fauna, played a major role in the post-glacial colonization of Andean lakes and streams, and contributes key species to the food webs. Galaxiid species occupy an enormous latitudinal gradient, show a wide variety of life history patterns and are the southernmost native freshwater fishes of the world. Knowledge of South American galaxiids has improved notably, but new challenges arise due to climate change, biological invasions, damming, aquaculture and contamination. In this changing environment, the future of South American galaxiids should be carefully considered as a legacy of the old Gondwana and a unique attribute of the freshwaters of southern South America.

Invasive species and postglacial colonization: their effects on the genetic diversity of a Patagonian fish

The present distribution of Patagonian species is the result of a complex history involving Quaternary refugial populations, Holocene range expansions and demographic changes occurring during the Anthropocene. Invasive salmonids were introduced in Patagonia during the last century, occupying most rivers and lakes, preying on and competing with native species, including the fish Galaxias platei. Here, we used G. platei as a case study to understand how long-term (i.e. population differentiation during the Holocene) and short-term historical processes (salmonid introductions) affect genetic diversity. Using a suite of microsatellite markers, we found that the number of alleles is negatively correlated with the presence of salmonids (short-term processes), with G. platei populations from lakes with salmonids exhibiting significantly lower genetic diversity than populations from lakes without salmonids. Simulations (100 years backwards) showed that this difference in genetic diversity can be explained by a 99% reduction in population size. Allelic richness and observed heterozygosities were also negatively correlated with the presence of salmonids, but also positively correlated with long-term processes linked to Quaternary glaciations. Our results show how different genetic parameters can help identify processes taking place at different scales and their importance in terms of conservation.

Past, present, and future of a freshwater fish metapopulation in a threatened landscape

It is well documented that hydropower plants can affect the dynamics of fish populations through landscape alterations and the creation of new barriers. Less emphasis has been placed on the examination of the genetic consequences for fish populations of the construction of dams. The relatively few studies that focus on genetics often do not consider colonization history and even fewer tend to use this information for conservation purposes. As a case study, we used a 3-pronged approach to study the influence of historical processes, contemporary landscape features, and potential future anthropogenic changes in landscape on the genetic diversity of a fish metapopulation. Our goal was to identify the metapopulation's main attributes, detect priority areas for conservation, and assess the consequences of the construction of hydropower plants for the persistence of the metapopulation. We used microsatellite markers and coalescent approaches to examine historical colonization processes, traditional population genetics, and simulations of future populations under alternate scenarios of population size reduction and gene flow. Historical gene flow appeared to have declined relatively recently and contemporary populations appeared highly susceptible to changes in landscape. Gene flow is critical for population persistence. We found that hydropower plants could lead to a rapid reduction in number of alleles and to population extirpation 50-80 years after their construction. More generally, our 3-pronged approach for the analyses of empirical genetic data can provide policy makers with information on the potential impacts of landscape changes and thus lead to more robust conservation efforts.

Spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus revealed by museum genomics

Analyzing genetic variation through time and space is important to identify key evolutionary and ecological processes in populations. However, using contemporary genetic data to infer the dynamics of genetic diversity may be at risk of a bias, as inferences are performed from a set of extant populations, setting aside unavailable, rare, or now extinct lineages. Here, we took advantage of new developments in next-generation sequencing to analyze the spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus, a steppic Southwestern-Palearctic species. We applied a recently developed hybridization capture (hyRAD) protocol that allows retrieving orthologous sequences even from degraded DNA characteristic of museum specimens. We identified single nucleotide polymorphisms in 68 historical and 51 modern samples in order to (i) unravel the spatial genetic structure across part of the species distribution and (ii) assess the loss of genetic diversity over the past century in Swiss populations. Our results revealed (i) the presence of three potential glacial refugia spread across the European continent and converging spatially in the Alpine area. In addition, and despite a limited population sample size, our results indicate (ii) a loss of allelic richness in contemporary Swiss populations compared to historical populations, whereas levels of expected heterozygosities were not significantly different. This observation is compatible with an increase in the bottleneck magnitude experienced by central European populations of O. decorus following human-mediated land-use change impacting steppic habitats. Our results confirm that application of hyRAD to museum samples produces valuable information to study genetic processes across time and space.

Historical isolation and contemporary gene flow drive population diversity of the brown alga Sargassum thunbergii along the coast of China

BACKGROUND

Long-term survival in isolated marginal seas of the China coast during the late Pleistocene ice ages is widely believed to be an important historical factor contributing to population genetic structure in coastal marine species. Whether or not contemporary factors (e.g. long-distance dispersal via coastal currents) continue to shape diversity gradients in marine organisms with high dispersal capability remains poorly understood. Our aim was to explore how historical and contemporary factors influenced the genetic diversity and distribution of the brown alga Sargassum thunbergii, which can drift on surface water, leading to long-distance dispersal.

RESULTS

We used 11 microsatellites and the plastid RuBisCo spacer to evaluate the genetic diversity of 22 Sargassum thunbergii populations sampled along the China coast. Population structure and differentiation was inferred based on genotype clustering and pairwise F ST and allele-frequency analyses. Integrated genetic analyses revealed two genetic clusters in S. thunbergii that dominated in the Yellow-Bohai Sea (YBS) and East China Sea (ECS) respectively. Higher levels of genetic diversity and variation were detected among populations in the YBS than in the ECS. Bayesian coalescent theory was used to estimate contemporary and historical gene flow. High levels of contemporary gene flow were detected from the YBS (north) to the ECS (south), whereas low levels of historical gene flow occurred between the two regions.

CONCLUSIONS

Our results suggest that the deep genetic divergence in S. thunbergii along the China coast may result from long-term geographic isolation during glacial periods. The dispersal of S. thunbergii driven by coastal currents may facilitate the admixture between southern and northern regimes. Our findings exemplify how both historical and contemporary forces are needed to understand phylogeographical patterns in coastal marine species with long-distance dispersal.

Genomic patterns of diversity and divergence of two introduced salmonid species in Patagonia, South America

Invasive species have become widespread in aquatic environments throughout the world, yet there are few studies that have examined genomic variation of multiple introduced species in newly colonized environments. In this study, we contrast genomic variation in two salmonid species (anadromous Chinook Salmon, Oncorhynchus tshawytscha, 11,579 SNPs and resident Brook Charr Salvelinus fontinalis, 13,522 SNPs) with differing invasion success after introduction to new environments in South America relative to populations from their native range in North America. Estimates of genetic diversity were not significantly different between introduced and source populations for either species, indicative of propagule pressure that has been shown to maintain diversity in founding populations relative to their native range. Introduced populations also demonstrated higher connectivity and gene flow than those in their native range. Evidence for candidate loci under divergent selection was observed, but was limited to specific introduced populations and was not widely evident. Patterns of genomic variation were consistent with general dispersal potential of each species and therefore also the notion that life history variation may contribute to both invasion success and subsequent genetic structure of these two salmonids in Patagonia.

Anthropogenic Stressors Shape Genetic Structure: Insights from a Model Freshwater Population along a Land Use Gradient

Environmental pollution including mutagens from wastewater effluents and discontinuity by man-made barriers are considered typical anthropogenic pressures on microevolutionary processes that are responsible for the loss of biodiversity in aquatic ecosystems. Here, we tested for the effects of wastewater treatment plants (WWTPs), weirs and other stressors on the invertebrate species Gammarus pulex at the population genetic level combining evolutionary ecotoxicology, body burden analysis and testing for exposure to mutagens. Exposure to chemical pollution alone and in combination with the presence of weirs resulted in a depression of allelic richness in native G. pulex populations. Our results suggest that the input of a mutagenic effluent from a WWTP resulted in a strong increase in private alleles over the affected populations. In addition, the presence of weirs along the river disrupted the migration across the river and thus the gene flow between G. pulex upstream and downstream. This study provides strong evidence that the assessment of genetic variation including private alleles together with the contamination of mutagenic and nonmutagenic chemical pollution offers new insights into the regulation of genetic population structure and highlights the relevance of emerging anthropogenic pressures at the genetic level.

A portrait of a sucker using landscape genetics: how colonization and life history undermine the idealized dendritic metapopulation

Dendritic metapopulations have been attributed unique properties by in silico studies, including an elevated genetic diversity relative to a panmictic population of equal total size. These predictions have not been rigorously tested in nature, nor has there been full consideration of the interacting effects among contemporary landscape features, colonization history and life history traits of the target species. We tested for the effects of dendritic structure as well as the relative importance of life history, environmental barriers and historical colonization on the neutral genetic structure of a longnose sucker (Catostomus catostomus) metapopulation in the Kogaluk watershed of northern Labrador, Canada. Samples were collected from eight lakes, genotyped with 17 microsatellites, and aged using opercula. Lakes varied in differentiation, historical and contemporary connectivity, and life history traits. Isolation by distance was detected only by removing two highly genetically differentiated lakes, suggesting a lack of migration-drift equilibrium and the lingering influence of historical factors on genetic structure. Bayesian analyses supported colonization via the Kogaluk's headwaters. The historical concentration of genetic diversity in headwaters inferred by this result was supported by high historical and contemporary effective sizes of the headwater lake, T-Bone. Alternatively, reduced allelic richness in headwaters confirmed the dendritic structure's influence on gene flow, but this did not translate to an elevated metapopulation effective size. A lack of equilibrium and upstream migration may have dampened the effects of dendritic structure. We suggest that interacting historical and contemporary factors prevent the achievement of the idealized traits of a dendritic metapopulation in nature.

Phylogeographic model selection leads to insight into the evolutionary history of four-eyed frogs

Phylogeographic research investigates biodiversity at the interface between populations and species, in a temporal and geographic context. Phylogeography has benefited from analytical approaches that allow empiricists to estimate parameters of interest from the genetic data (e.g., θ = 4Neμ, population divergence, gene flow), and the widespread availability of genomic data allow such parameters to be estimated with greater precision. However, the actual inferences made by phylogeographers remain dependent on qualitative interpretations derived from these parameters' values and as such may be subject to overinterpretation and confirmation bias. Here we argue in favor of using an objective approach to phylogeographic inference that proceeds by calculating the probability of multiple demographic models given the data and the subsequent ranking of these models using information theory. We illustrate this approach by investigating the diversification of two sister species of four-eyed frogs of northeastern Brazil using single nucleotide polymorphisms obtained via restriction-associated digest sequencing. We estimate the composite likelihood of the observed data given nine demographic models and then rank these models using Akaike information criterion. We demonstrate that estimating parameters under a model that is a poor fit to the data is likely to produce values that lead to spurious phylogeographic inferences. Our results strongly imply that identifying which parameters to estimate from a given system is a key step in the process of phylogeographic inference and is at least as important as being able to generate precise estimates of these parameters. They also illustrate that the incorporation of model uncertainty should be a component of phylogeographic hypothesis tests.

Non-Antarctic notothenioids: Past phylogenetic history and contemporary phylogeographic implications in the face of environmental changes

The non-Antarctic Notothenioidei families, Bovichtidae, Pseudaphritidae and Eleginopsidae, diverged early from the main notothenioid lineage. They are important in clarifying the early evolutionary processes that triggered notothenioid evolution in the Antarctic. The early-diverged group represents 8% of all notothenioid species and never established themselves on the Antarctic shelf. Most attention has been paid to the Antarctic notothenioids and their limited physiological tolerance to climate change and increased temperatures. In this review, we discuss key life history traits that are characteristic of the non-Antarctic early-diverged notothenioid taxa as well as the genetic resources and population differentiation information available for this group. We emphasise the population fitness and dynamics of these species and indicate how resource management and conservation of the group can be strengthened through an integrative approach. Both Antarctic waters and the non-Antarctic regions face rapid temperature rises combined with strong anthropogenic exploitation. While it is expected that early-diverged notothenioid species may have physiological advantages over high Antarctic species, it is difficult to predict how climate changes might alter the geographic range, behaviour, phenology and ultimately genetic variability of these species. It is possible, however, that their high degree of endemism and dependence on local environmental specificities to complete their life cycles might enhance their vulnerability.