Comparisons of the genetic structure of Squalius populations (Teleostei, Cyprinidae) from rivers with contrasting histories, drainage areas and climatic conditions based on two molecular markers

Phylogeography and ecological differentiation of strictly Mediterranean taxa: the case of the Iberian endemic Odontites recordonii

PREMISE

Ecological drivers for genetic differentiation in Mediterranean climates are still underexplored. We have used the strictly Mediterranean endemic Odontites recordonii as a model species to address this question. This species is one of the three Iberian representatives of the O. vernus group, which are morphologically similar. Thus, it was additionally necessary to clarify their phylogenetic relationships.

METHODS

We used amplified fragment length polymorphisms to reveal phylogenetic relationships within O. vernus group, and to reconstruct the phylogeographic patterns within O. recordonii. Additionally, ecological niche models were generated to detect refugia along the Quaternary climatic oscillations. And finally, alleles under natural selection were identified, and correlations between allele presences and environmental variables were calculated in order to shed light on the ecological drivers promoting differentiation.

RESULTS

The three species from the O. vernus group were recovered as distinct species. Three genetic groups were found within O. recordonii and a putative refugium was detected for each one. Eighty-one alleles could be under diversifying selection, and 58 alleles showed significant correlations with environmental variables, especially with temperature and precipitation seasonality and summer drought.

CONCLUSIONS

The three Iberian species of the O. vernus group are reciprocal monophyletic taxa. The three genetic groups of O. recordonii could have been restricted to narrow refugia during the Quaternary and displayed present distributions in accordance with bioclimatic conditions. Temperature and precipitation seasonality and the intensity of summer drought are definitory climatic parameters of Mediterranean-type climates, and they could have acted as drivers of genetic differentiation on O. recordonii.

Genomic data and multi-species demographic modelling uncover past hybridization between currently allopatric freshwater species

Evidence for ancient interspecific gene flow through hybridization has been reported in many animal and plant taxa based on genetic markers. The study of genomic patterns of closely related species with allopatric distributions allows the assessment of the relative importance of vicariant isolating events and past gene flow. Here, we investigated the role of gene flow in the evolutionary history of four closely related freshwater fish species with currently allopatric distributions in western Iberian rivers-Squalius carolitertii, S. pyrenaicus, S. torgalensis and S. aradensis-using a population genomics dataset of 23,562 SNPs from 48 individuals, obtained through genotyping by sequencing (GBS). We uncovered a species tree with two well-differentiated clades: (i) S. carolitertii and S. pyrenaicus; and (ii) S. torgalensis and S. aradensis. By using D-statistics and demographic modelling based on the site frequency spectrum, comparing alternative demographic scenarios of hybrid origin, secondary contact and isolation, we found that the S. pyrenaicus North lineage is likely the result of an ancient hybridization event between S. carolitertii (contributing ~84%) and S. pyrenaicus South lineage (contributing ~16%), consistent with a hybrid speciation scenario. Furthermore, in the hybrid lineage, we identify outlier loci potentially affected by selection favouring genes from each parental lineage at different genomic regions. Our results suggest that ancient hybridization can affect speciation and that freshwater fish species currently in allopatry are useful to study these processes.

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.

Genetic legacy of tertiary climatic change: a case study of two freshwater loaches, Schistura fasciolata and Pseudogastromyzon myersi, in Hong Kong

The high biodiversity and strong population structure of freshwater fauna has often been attributed to historical geological and climatic alterations. The impact of these historical changes on obligate freshwater species on a small geographical scale has not been well understood due to the lack of fine-scale comparative phylogeographic studies. Strong population structure has been reported in a goby and a caridean shrimp in Hong Kong, a small but highly developed city in South China, but the common drivers of population differentiation in freshwater fauna in this region remain unclear. This study examined the fine-scale phylogeographic patterns of two freshwater loaches, Schistura fasciolata and Pseudogastromyzon myersi in Hong Kong, using sequence data of mitochondrial control region and two nuclear markers (interphotoreceptor retinoid binding protein gene 2 and ribosomal protein S13 gene). Results show that they exhibit pronounced population structure as supported by high and significant Φ_ST. Phylogenetic analyses based on the control region reveal six and three distinct lineages in S. fasciolata and P. myersi, respectively. Phylogeographic structure of both species generally follows the paleodrainage pattern, though P. myersi shows a shallower structure on the Mainland, perhaps due to their higher mobility. Most of these lineages diverged during the Pliocene and Late Pleistocene, a period with marked sea-level fluctuations. In a broader context, this suggests that sea-level fluctuation played an important role in shaping even the fine-scale population structure of freshwater fish in South China, implying that the genetic diversity of this fauna may be higher than expected.

Multilocus approach reveals cryptic lineages in the goby Rhinogobius duospilus in Hong Kong streams: Role of paleodrainage systems in shaping marked population differentiation in a city

Drainage history is a well-demonstrated factor that influences the population structure of freshwater inhabitants over a broad geographic scale. However, there has been little research undertaken on such a relationship with freshwater fish on a small geographical scale, especially in Asia. In this study, we investigated the role of local, small drainage systems in affecting the population genetic structure of a freshwater goby, Rhinogobius duospilus, in Hong Kong streams using a multilocus approach. Analyses on nine genetic markers (2 mitochondrial and 7 nuclear markers, including 5 microsatellite markers) reveal prominent and intensive genetic structuring (2.1-5.4% mtDNA sequence divergence) in R. duospilus in Hong Kong. The lineages and clusters recovered from mtDNA data and assignment analysis of nuclear markers coincide with the paleodrainage networks. Furthermore, marked population subdivision between streams located on different side branches (<20km apart) within the same paleodrainage area is observed and gene flow occurs only between closely situated streams that share common paleodrainage confluences. In an extreme case, gene flow is limited between streams that are less than 5km apart. Apparently, such an intensive population structure is attributed to the regional paleodrainage pattern, together with the highly sedentary life style of R. duospilus, which reduces contemporary gene flow and dispersal between populations in neighbouring streams.

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.

Phylogeography, historical demography and habitat suitability modelling of freshwater fishes inhabiting seasonally fluctuating Mediterranean river systems: a case study using the Iberian cyprinid Squalius valentinus

The Mediterranean freshwater fish fauna has evolved under constraints imposed by the seasonal weather/hydrological patterns that define the Mediterranean climate. These conditions have influenced the genetic and demographic structure of aquatic communities since their origins in the Mid-Pliocene. Freshwater species in Mediterranean-type climates will likely constitute genetically well-differentiated populations, to varying extents depending on basin size, as a consequence of fragmentation resulting from drought/flood cycles. We developed an integrative framework to study the spatial patterns in genetic diversity, demographic trends, habitat suitability modelling and landscape genetics, to evaluate the evolutionary response of Mediterranean-type freshwater fish to seasonal fluctuations in weather. To test this evolutionary response, the model species used was Squalius valentinus, an endemic cyprinid of the Spanish Levantine area, where seasonal weather fluctuations are extreme, although our findings may be extrapolated to other Mediterranean-type species. Our results underscore the significant role of the Mediterranean climate, along with Pleistocene glaciations, in diversification of S. valentinus. We found higher nuclear diversity in larger drainage basins, but higher mitochondrial diversity correlated to habitat suitability rather than basin size. We also found strong correlation between genetic structure and climatic factors associated with Mediterranean seasonality. Demographic and migration analyses suggested population expansion during glacial periods that also contributed to the current genetic structure of S. valentinus populations. The inferred models support the significant contribution of precipitation and temperature to S. valentinus habitat suitability and allow recognizing areas of habitat stability. We highlight the importance of stable habitat conditions, fostered by typical karstic springs found on the Mediterranean littoral coasts, for the preservation of freshwater species inhabiting seasonally fluctuating river systems.

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.

Comparative cytogenetics of two endangered leuciscine fish, Squalius aradensis and S. torgalensis (Teleostei, Cyprinidae), from the Iberian Peninsula

In this study, the description of the karyotypes of the endangered chubs Squalius aradensis (Coelho, Bogutskaya, Rodrigues and Collares-Pereira, 1998) and Squalius torgalensis (Coelho, Bogutskaya, Rodrigues and Collares-Pereira, 1998) is presented by means of conventional (Giemsa-staining, Chromomycin A3 (CMA3)-fluorescence, Silver-impregnation (Ag-NORs)) and molecular (fluorescence in situ hybridization (FISH) with 18S rDNA probe) protocols. These endemic sister-species have an allopatric but adjacent distribution in the most southwestern part of the Iberian Peninsula. Diploid chromosome number was invariably 2n = 50 and karyotypes of both species were grossly similar, composed of metacentric and submetacentric elements with a reduced number of acrocentric pairs. Sequential staining using FISH with an 18S rDNA probe, CMA3 and Ag-NORs treatments revealed consistent positive signals located at the end of the short arms of a submetacentric chromosome pair, likely homologous in both species. While providing useful cytogenetic comparative data against other members of the genus Squalius Bonaparte, 1837, the work aimed to draw attention towards the conservation of two narrow-range and highly confined fish species.