Evolutionary history and population genetics of a cyprinid fish (Iberochondrostoma olisiponensis) endangered by introgression from a more abundant relative

Fish diversification at the pace of geomorphological changes: evolutionary history of western Iberian Leuciscinae (Teleostei: Leuciscidae) inferred from multilocus sequence data

The evolutionary history of western Iberian Leuciscinae, obligatory freshwater fish, is directly linked to the evolution of the hydrographic network of the Iberian Peninsula after its isolation from the rest of Europe, which involved dramatic rearrangements such as the transition from endorheic lakes to open basins draining to the Atlantic. Previous phylogenetic research on western Iberian leuciscines, using mainly mitochondrial DNA and more recently one or two nuclear genes, has found contradictory results and there remain many unresolved issues regarding species relationships, taxonomy, and evolutionary history. Moreover, there is a lack of integration between phylogenetic and divergence time estimates and information on the timing of geomorphological changes and paleobasin rearrangements in the Iberian Peninsula. This study presents the first comprehensive fossil-calibrated multilocus coalescent species tree of western Iberian Leuciscinae (including 14 species of Achondrostoma, Iberochondrostoma, Pseudochondrostoma and Squalius endemic to the Iberian Peninsula, seven of which endemic to Portugal) based on seven nuclear genes, and from which we infer their biogeographic history by comparing divergence time estimates to known dated geological events. The phylogenetic pattern suggests slow-paced evolution of leuciscines during the Early-Middle Miocene endorheic phase of the main Iberian river basins, with the shift to exorheism in the late Neogene-Quaternary allowing westward dispersals that resulted in many cladogenetic events and a high rate of endemism in western Iberia. The results of this study also: (i) confirm the paraphyly of S. pyrenaicus with respect to S. carolitertii, and thus the possible presence of a new taxon in the Portuguese Tagus currently assigned to S. pyrenaicus; (ii) support the taxonomic separation of the Guadiana and Sado populations of S. pyrenaicus; (iii) show the need for further population sampling and taxonomic research to clarify the phylogenetic status of A. arcasii from the Minho basin and of the I. lusitanicum populations in the Sado and Tagus basins; and (iv) indicate that A. occidentale, I. olisiponensis and P. duriensis are the most ancient lineages within their respective genera.

Broad-scale sampling of primary freshwater fish populations reveals the role of intrinsic traits, inter-basin connectivity, drainage area and latitude on shaping contemporary patterns of genetic diversity

Background. Worldwide predictions suggest that up to 75% of the freshwater fish species occurring in rivers with reduced discharge could be extinct by 2070 due to the combined effect of climate change and water abstraction. The Mediterranean region is considered to be a hotspot of freshwater fish diversity but also one of the regions where the effects of climate change will be more severe. Iberian cyprinids are currently highly endangered, with over 68% of the species raising some level of conservation concern.

Methods. During the FISHATLAS project, the Portuguese hydrographical network was extensively covered (all the 34 river basins and 47 sub-basins) in order to contribute with valuable data on the genetic diversity distribution patterns of native cyprinid species. A total of 188 populations belonging to 16 cyprinid species of Squalius, Luciobarbus, Achondrostoma, Iberochondrostoma, Anaecypris and Pseudochondrostoma were characterized, for a total of 3,678 cytochrome b gene sequences.

Results. When the genetic diversity of these populations was mapped, it highlighted differences among populations from the same species and between species with identical distribution areas. Factors shaping the contemporary patterns of genetic diversity were explored and the results revealed the role of latitude, inter-basin connectivity, migratory behaviour, species maximum size, species range and other species intrinsic traits in determining the genetic diversity of sampled populations. Contrastingly, drainage area and hydrological regime (permanent vs. temporary) seem to have no significant effect on genetic diversity. Species intrinsic traits, maximum size attained, inter-basin connectivity and latitude explained over 30% of the haplotype diversity variance and, generally, the levels of diversity were significantly higher for smaller sized species, from connected and southerly river basins.

Discussion. Targeting multiple co-distributed species of primary freshwater fish allowed us to assess the relative role of historical versus contemporary factors affecting genetic diversity. Since different patterns were detected for species with identical distribution areas we postulate that contemporary determinants of genetic diversity (species’ intrinsic traits and landscape features) must have played a more significant role than historical factors. Implications for conservation in a context of climate change and highly disturbed habitats are detailed, namely the need to focus management and conservation actions on intraspecific genetic data and to frequently conduct combined genetic and demographic surveys.

Cenozoic tectonic and climatic events in southern Iberian Peninsula: Implications for the evolutionary history of freshwater fish of the genus Squalius (Actinopterygii, Cyprinidae)

Southern Iberian freshwater ecosystems located at the border between the European and African plates represent a tectonically complex region spanning several geological ages, from the uplifting of the Betic Mountains in the Serravalian-Tortonian periods to the present. This area has also been subjected to the influence of changing climate conditions since the Middle-Upper Pliocene when seasonal weather patterns were established. Consequently, the ichthyofauna of southern Iberia is an interesting model system for analyzing the influence of Cenozoic tectonic and climatic events on its evolutionary history. The cyprinids Squalius malacitanus and Squalius pyrenaicus are allopatrically distributed in southern Iberia and their evolutionary history may have been defined by Cenozoic tectonic and climatic events. We analyzed MT-CYB (510 specimens) and RAG1 (140 specimens) genes of both species to reconstruct phylogenetic relationships and to estimate divergence times and ancestral distribution ranges of the species and their populations. We also assessed their levels of genetic structure and diversity as well as the amount of gene flow between populations. To investigate recent paleogeographical and climatic factors in southern Iberia, we modeled changes-through-time in sea level from the LGM to the present. Phylogenetic, geographic and population structure analyses revealed two well-supported species (S. malacitanus and S. pyrenaicus) in southern Iberia and two subclades (Atlantic and Mediterranean) within S. malacitanus. The origin of S. malacitanus and the separation of its Atlantic and Mediterranean populations occurred during the Serravalian-Tortonian and Miocene-Pliocene periods, respectively. These divergence events occurred in the Middle Pliocene and Pleistocene in S. pyrenaicus. In both species, Atlantic basins possessed populations with higher genetic diversity than Mediterranean, which may be explained by the Janda Lagoon. The isolation of S. malacitanus was earlier and related to the rising of the Betic Mountains. Divergence of its Atlantic and Mediterranean populations was associated with the creation of the freshwater systems of southern Iberia close to the Gibraltar Strait. The presence of S. pyrenaicus in southern Iberia may be the result of recent colonization associated with river capture, as demonstrated our biogeographic reconstruction.

Genetic Structure is Inconsistent with Subspecies Designations in the Western Massasauga Sistrurus tergeminus

The massasauga Sistrurus catenatus was historically divided into three subspecies, but this long-standing taxonomy has recently been called into question. Genetic research now recognizes a split of the species into the eastern massasauga S. catenatus and western massasauga S. tergeminus, with the latter split into two subspecies, the desert massasauga S. t. edwardsii and the prairie massasauga S. t. tergeminus. Although the distinction between geographically isolated populations of S. catenatus and S. tergeminus is well-supported genetically, the geographic relationships among populations of S. t. tergeminus and S. t. edwardsii remain unresolved because of incomplete sampling throughout the species’ range. This poses a difficult challenge for conservation and management of this species. Sistrurus t. tergeminus does not have state or federal conservation status, but S. t. edwardsii has been petitioned for listing under the U.S. Endangered Species Act. In this study, we used nuclear and mitochondrial DNA from 52 individuals from 7 states to explore the taxonomic and geographic relationships between S. t. tergeminus and S. t. edwardsii populations. Maximum likelihood and Bayesian inference frameworks for both nuclear and mtDNA genes indicated that S. t. tergeminus and S. t. edwardsii populations were genetically indistinguishable. However, at the species level, we did find eight well-supported mtDNA clades within S. tergeminus, including individuals from five peripheral populations in 1) Arizona and western New Mexico, 2) Colorado and Kansas, 3) Missouri, 4) Oklahoma, and 5) southern Texas. These peripherally isolated populations surrounded a larger population of individuals from north-central Texas, New Mexico, and Oklahoma that was contiguous with three additional genetically distinct populations. We conclude that the putative subspecies S. t. tergeminus and S. t. edwardsii, as currently defined, most likely represent polytypic phenotypes of S. tergeminus rather than discrete taxonomic entities. Instead, we suggest that S. tergeminus existed historically as a large, contiguous collection of populations that only recently became fragmented into several, as opposed to two, potentially discrete taxonomic entities.

Semi-permeable species boundaries in Iberian barbels (Barbus and Luciobarbus, Cyprinidae)

BACKGROUND

The evolution of species boundaries and the relative impact of selection and gene flow on genomic divergence are best studied in populations and species pairs exhibiting various levels of divergence along the speciation continuum. We studied species boundaries in Iberian barbels, Barbus and Luciobarbus, a system of populations and species spanning a wide degree of genetic relatedness, as well as geographic distribution and range overlap. We jointly analyze multiple types of molecular markers and morphological traits to gain a comprehensive perspective on the nature of species boundaries in these cyprinid fishes.

RESULTS

Intraspecific molecular and morphological differentiation is visible among many populations. Genomes of all sympatric species studied are porous to gene flow, even if they are not sister species. Compared to their allopatric counterparts, sympatric representatives of different species share alleles and show an increase in all measures of nucleotide polymorphism (S, Hd, K, π and θ). High molecular diversity is particularly striking in L. steindachneri from the Tejo and Guadiana rivers, which co-varies with other sympatric species. Interestingly, different nuclear markers introgress across species boundaries at various levels, with distinct impacts on population trees. As such, some loci exhibit limited introgression and population trees resemble the presumed species tree, while alleles at other loci introgress more freely and population trees reflect geographic affinities and interspecific gene flow. Additionally, extent of introgression decreases with increasing genetic divergence in hybridizing species pairs.

CONCLUSIONS

We show that reproductive isolation in Iberian Barbus and Luciobarbus is not complete and species boundaries are semi-permeable to (some) gene flow, as different species (including non-sister) are exchanging genes in areas of sympatry. Our results support a speciation-with-gene-flow scenario with heterogeneous barriers to gene flow across the genome, strengthening with genetic divergence. This is consistent with observations coming from other systems and supports the notion that speciation is not instantaneous but a gradual process, during which different species are still able to exchange some genes, while selection prevents gene flow at other loci. We also provide evidence for a hybrid origin of a barbel ecotype, L. steindachneri, suggesting that ecology plays a key role in species coexistence and hybridization in Iberian barbels. This ecotype with intermediate, yet variable, molecular, morphological, trophic and ecological characteristics is the local product of introgressive hybridization of L. comizo with up to three different species (with L. bocagei in the Tejo, with L. microcephalus and L. sclateri in the Guadiana). In spite of the homogenizing effects of ongoing gene flow, species can still be discriminated using a combination of morphological and molecular markers. Iberian barbels are thus an ideal system for the study of species boundaries, since they span a wide range of genetic divergences, with diverse ecologies and degrees of sympatry.

Metapopulations in temporary streams - the role of drought-flood cycles in promoting high genetic diversity in a critically endangered freshwater fish and its consequences for the future

Genetic factors have direct and indirect impacts in the viability of endangered species. Assessing their genetic diversity levels and population structure is thus fundamental for conservation and management. In this paper we use mitochondrial and nuclear markers to address phylogeographic and demographic data on the critically endangered Anaecypris hispanica, using a broad sampling set which covered its known distribution area in the Iberian Peninsula. Our results showed that the populations of A. hispanica are strongly differentiated (high and significant ФST and FST values, corroborated by the results from AMOVA and SAMOVA) and genetically diversified. We suggest that the restricted gene flow between populations may have been potentiated by ecological, hydrological and anthropogenic causes. Bayesian skyline plots revealed a signal for expansion for all populations (tMRCA between 68kya and 1.33Mya) and a genetic diversity latitudinal gradient was detected between the populations from the Upper (more diversified) and the Lower (less diversified) Guadiana river basin. We postulate a Pleistocenic westwards colonization route for A. hispanica in the Guadiana river basin, which is in agreement with the tempo and mode of paleoevolution of this drainage. The colonization of River Guadalquivir around 60kya with migrants from the Upper Guadiana, most likely by stream capture, is also suggested. This study highlights the view that critically endangered species facing range retreats (about 47% of its known populations have disappeared in the last 15years) are not necessarily small and genetically depleted. However, the extinction risk is not negligible since A. hispanica faces the combined effect of several deterministic and stochastic negative factors and, moreover, recolonization events after localized extinctions are very unlikely to occur due to the strong isolation of populations and to the patchily ecologically-conditioned distribution of fish. The inferred species distribution models highlight the significant contribution of temperature seasonality and isothermality to A. hispanica occurrence in Guadiana environments and emphasize the importance of stable climatic conditions for the preservation of this species. Given the strong population structure, high percentage of private haplotypes and virtual absence of inter-basin gene flow we suggest that each A. hispanica population should be considered as an independent Operational Conservation Unit and that ex-situ and in-situ actions should be conducted in parallel to allow for the long-term survival of the species and the preservation of the genetic integrity of its populations.