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.

Comparative pattern of genetic structure in two Mediterranean killifishes Aphanius fasciatus and Aphanius iberus inferred from both mitochondrial and nuclear data.

In this study, genetic variation was assessed in Aphanius fasciatus and Aphanius iberus characterized by similar ecological traits but with very different distribution ranges in the Mediterranean area. Five populations of A. iberus and five of A. fasciatus were analysed using five polymorphic microsatellite loci and partial mitochondrial control region (D-loop) sequences. Congruent results were found with both nuclear and mitochondrial molecular markers. The results showed that similar levels of genetic divergence, based on mitochondrial control region sequences, are present among populations of A. iberus and among populations of A. fasciatus despite the very different geographic distance existing among the examined populations of the two species (low geographic distance in A. iberus and high in A. fasciatus). A possible explanation could be that the populations of A. iberus were isolated for a longer time than the populations of A. fasciatus supporting the hypothesis that the split in the lineage leading to A. iberus is older than the split in the lineage leading to A. fasciatus. The possibility that the wide circum-Mediterranean distribution of A. fasciatus ensures the high connectivity of its populations, preventing, in some cases, local differentiation, however, cannot be ruled out.

Comparative historical biogeography of three groups of Nearctic freshwater fishes across central Mexico.

Biogeographic patterns of the three main Nearctic groups of continental fishes inhabiting river drainages in central Mexico (livebearing goodeids, southern Mexican notropins and species of Algansea, the last two representing independent lineages of cyprinids) were obtained and compared by following two approaches: an estimate of divergence times and using a well-defined biogeographic method. Three concordant biogeographic events were identified among the three groups, showing some evidence of a partially congruent evolutionary history. The analysed groups show at least three independent colonization events into central Mexico: two western routes, followed by the Goodeinae and members of Algansea, and an early Plateau route followed by southern notropins. The most recent common ancestor (MRCA) of each of the three freshwater fish groups diversified in central Mexico in the Late Miocene. The lack of a strong congruence in their biogeographic patterns, and the differences in species richness among the three clades might be evidence for distinct patterns of diversification.

Different stocks of brook lamprey in Spain and their origin from Lampetra fluviatilis at two distinct times and places.

Using the mitochondrial non-coding region I, it was shown that the two Spanish Lampetra planeri populations (Cares-Deva and Olabidea-Ugarana) correspond to different genetic units. The Cares-Deva population is probably a recent offshoot of Lampetra fluviatilis, being the Olabidea-Ugarana population less diverse and of older origin.

Genetic diversity shaped by historical and recent factors in the live-bearing twoline skiffia Neotoca bilineata.

The endangered twoline skiffia Neotoca bilineata, a viviparous fish of the subfamily Goodeinae, endemic to central Mexico (inhabiting two basins, Cuitzeo and Lerma-Santiago) was evaluated using genetic and habitat information. The genetic variation of all remaining populations of the species was analysed using both mitochondrial and microsatellite markers and their habitat conditions were assessed using a water quality index (I(WQ)). An 80% local extinction was found across the distribution of N. bilineata. The species was found in three of the 16 historical localities plus one previously unreported site. Most areas inhabited by the remaining populations had I(WQ) scores unsuitable for the conservation of freshwater biodiversity. Populations showed low but significant genetic differentiation with both markers (mtDNA φ(ST) = 0.076, P < 0.001; microsatellite F(ST) = 0.314, P < 0.001). Borbollon, in the Cuitzeo Basin, showed the highest level of differentiation and was identified as a single genetic unit by Bayesian assignment methods. Rio Grande de Morelia and Salamanca populations showed the highest genetic diversity and also a high migration rate facilitated by an artificial channel that connected the two basins. Overall, high genetic diversity values were observed compared with other freshwater fishes (average N(a) = 16 alleles and loci and mean ±S.D. H(o) = 0.63 ± 0.10 and nucleotide diversity π = 0.006). This suggests that the observed genetic diversity has not diminished as rapidly as the species' habitat destruction. No evidence of correlation between habitat conditions and genetic diversity was found. The current pattern of genetic diversity may be the result of both historical factors and recent modifications of the hydrological system. The main threat to the species may be the rapid habitat deterioration and associated demographic stochasticity rather than genetic factors.

The evolutionary history of the allopolyploid Squalius alburnoides (Cyprinidae) complex in the northern Iberian Peninsula.

Understanding the population structure, population dynamics and processes that give rise to polyploidy and helps to maintain it is central to our knowledge of the evolution of asexual vertebrates. Previous studies revealed high genetic diversity and several reproductive pathways in the southern populations of the Squalius alburnoides hybrid complex. In contrast, lower genetic variability and the associated limited chance of introducing new genetic combinations may threaten the survival of the northern Mondego populations. We analysed the genetic diversity and structure of nine populations of S. alburnoides in the Iberian Peninsula using microsatellite loci to provide further insights on the evolutionary history of this complex. Special attention was given to the less-studied northern populations (Mondego and Douro basins). Marked population structure, a high frequency of private alleles and a high diversity of some biotypes in the Douro basin indicate that some northern populations may not be at high risk of extinction, contrary to what was expected. The genetic diversity found in the northern Douro populations contradicts the general trend of remarkable genetic impoverishment northwards that occurs in other species and regions. The results indicate the possible existence of a glacial refugium in the Rabaçal River, corroborating findings in other species of this region. Historical events seem to have affected the geographical patterns of genetic variability found among and within the northern and southern populations of this complex and contributed to different patterns of genome composition. Therefore, historical events might have a major role in the long-term persistence of some polyploid hybrid taxa.

Exploring the effect of microsatellite size homoplasy on reconstruction of phylogenetic relationships of picote splitfin Zoogoneticus quitzeoensis.

This study explores the effects of microsatellite size homoplasies on the reconstruction of phylogenetic relationships and estimates of population parameters as the fixation index (F(ST) ) using as a case study a truncated microsatellite from the picote splitfin Zoogoneticus quitzeoensis. The results suggest that the use of imperfect microsatellites may have only a minor effect in phylogenetic and population studies.

On the occurrence of Anaecypris hispanica, an extremely endangered Iberian endemism, in the Guadalquivir River basin.

The jarabugo Anaecypris hispanica, considered endemic to the Guadiana River basin, has been found in the Guadalquivir River. First genetic data showed a high degree of similarity to those of the Guadiana River populations. The genetic study recovered five different groups of haplotypes, the Guadalquivir River specimens belong to the largest and most widely extended group.

Phylogenetic analysis of peri-Mediterranean blennies of the genus Salaria: molecular insights on the colonization of freshwaters.

In this paper, the phylogenetic relationships of the marine blenny Salaria pavo and the freshwater S. fluviatilis and S. economidisi were analyzed using four molecular markers: the mitochondrial 12S rRNA, 16S rRNA, and the control region and the nuclear first intron of the S7 ribosomal protein. The monophyly of Salaria is supported, as well as that of S. pavo and that of all the freshwater members of Salaria. Thus, the present results support a single origin for all freshwater Mediterranean blenniids. Our results reject the placement of the species of Salaria in the genus Lipophrys as proposed in previous studies. Using a molecular clock calibrated with trans-Isthmian geminate blenniid species, the split between the ancestor of the freshwater lineage and the ancestor of S. pavo is tentatively placed in the Middle Miocene (well before the Messinian). The marine S. pavo displays a very low level of intraspecific sequence divergence consistent with a Pleistocene bottleneck. S. fluviatilis is a paraphyletic entity with S. economidisi nested within it. A Moroccan population of S. fluviatilis is more divergent than S. economidisi, both in nuclear and mitochondrial genes. Fish from Israel together with some Turkish samples represent the second oldest split. It is argued that these populations may represent cryptic species. Thus, further studies on the taxonomy of these freshwater blennies are urgently needed.

Body shape evolution among ploidy levels of the Squalius alburnoides hybrid complex (Teleostei, Cyprinidae).

Hybridization, ploidy level and genomic constitution may be important to respond to different environments, by producing different phenotypes and thus reducing competitive interaction. Through geometric morphometrics, we examined variation in body size and shape among biotypes of the Squalius alburnoides hybrid complex and their sperm donor (Squalius carolitertii). Results showed that S. carolitertii is significantly larger in size than the biotypes of the complex. No significant relationship was observed between ploidy and body size among S. alburnoides biotypes. Significant variation in body shape was found between S. carolitertii and S. alburnoides, and between tetraploids and the other biotypes. These differences in biotypes may reduce resource competition, highlighting the potential importance of resource availability favouring one biotype over another. In S. alburnoides, the adaptation to different trophic niches through modification of trophic morphology, body shapes, and feeding behaviour, may result from an increase in ploidy and genomic constitution. This adaptation may account also for the formation and maintenance of this nonsexual complex.

Hybridogenetic reproduction and maternal ancestry of polyploid Iberian fish: the Tropidophoxinellus alburnoides complex.

Iberian minnows collectively known as the Tropidophoxinellus alburnoides Steindachner complex comprise diploid and polyploid forms with highly female biased sex ratios. Previous investigators suggested that all-female clonal reproduction and interspecific hybridization may occur in this complex. We examined nuclear (allozymes) and cytoplasmic genes (mtDNA) to assess the evolutionary origins, relationships, and reproductive modes of T. alburnoides from western Spain. The multi-locus allozyme data clearly revealed the hybrid nature of all polyploid forms of this fish and some diploid forms as well. Diagnostic markers identified fish from the genus Leuciscus as the paternal ancestor of hybrids in the Duero and Guadiana River Basins. Additionally, analysis of nuclear markers revealed that hybridogenetic reproduction occurs in the diploid and triploid hybrids. The hybrids fully express the paternal Leuciscus genome and then discard it during oogenesis. Hybridogenetic ova contain only maternal nuclear genes and mtDNA from a non-hybrid T. alburnoides ancestor. Apparently diploid and triploid hybrids of T. alburnoides persist as sperm parasites on males of a sexually reproducing Leuciscus host species.

Phylogenetic relationships of Iberian cyprinids: systematic and biogeographical implications.

The phylogenetic relationships among all Iberian endemic cyprinids were inferred using the complete nucleotide sequence of the cytochrome b gene. The inferred molecular phylogeny included representatives from Central European, Asian and North African species, and is highly congruent with previous phylogenies based on osteological characters. Iberian cyprinids were grouped into only five, very speciose lineages (with the exception of the monotypic Anaecypris): Barbus, Luciobarbus, Chondrostoma, Leuciscus and Anaecypris. The existence of such a relatively small number of Iberian cyprinid lineages can be explained by the historical isolation of the Iberian Peninsula. North African and Asian barbels are the sister group of Iberian Luciobarbus, supporting a south-eastern route of colonization of the Iberian Peninsula for this subgenus. With leuciscins, Anaecypris hispanica was considered a relict species as it could not be related to any other Iberian cyprinid. The phylogenetic relationships among the main lineages of Iberian cyprinids based on cytochrome b sequence data supported the traditional division of the Cyprinidae into two subfamilies: Cyprininae and Leuciscinae.

Evolutionary history and speciation modes in the cyprinid genus Barbus.

Phylogenetic relationships and evolutionary patterns in the genus Barbus were examined through the analysis of the complete sequences of three mitochondrial genes: ATPases 8 and 6, which overlap slightly, and cytochrome b. This complex genus includes diploid as well as tetraploid and hexaploid species that are distributed throughout the Palaearctic, Ethiopian and Asiatic biogeographical regions. Given that genome duplication is an important evolutionary mechanism in eukaryotes, in the present report we attempt to describe its role in the evolution of the genus Barbus, as well as drawing systematic and phylogenetic conclusions. The phylogenetic results indicated the splitting of the current Barbus genus into five main mitochondrial lineages corresponding to (i) the genus Barbus sensu stricto (tetraploid, which is subdivided into the subgenera Barbus and Luciobarbus), (ii) the hexaploid species, (iii) the Ethiopian tetraploid species, (iv) the African diploid species, and (v) the Asian diploid species. The branching of 'foreign' genera as sister groups of some of these monophyletic assemblages (such as Aulopyge is to Barbus sensu stricto or Varicorhinus is to the hexaploid barbels) demonstrates the polyphyly of the group. Moreover, the relationships between the proposed lineages also show that genome duplication may be considered as a homoplasic character since it must have occurred over at least three independent periods and/or in three independent areas. In relation to the possible saltational evolutionary model for the polyploid species examined here, it was found that, although feasible at the nuclear level, the mitochondrial markers looked at do not appear to have undergone this type of evolution. Rather, they seem to have experienced more or less constant change through time.

Origin, radiation, dispersion and allopatric hybridization in the chub Leuciscus cephalus.

The phylogenetic relationships of 492 chub (Leuciscus cephalus) belonging to 89 populations across the species' range were assessed using 600 base pairs of cytochrome b. Furthermore, nine species belonging to the L. cephalus complex were also analysed (over the whole cytochrome b) in order to test potential allopatric hybridization with L. cephalus sensu stricto (i.e. the chub). Our results show that the chub includes four highly divergent lineages descending from a quick radiation that took place three million years ago. The geographical distribution of these lineages and results of the nested clade analysis indicated that the chub may have originated from Mesopotamia. Chub radiation probably occurred during an important vicariant event such as the isolation of numerous Turkish river systems, a consequence of the uplift of the Anatolian Plateau (formerly covered by a broad inland lake). Dispersion of these lineages arose from the changes in the European hydrographic network and, thus, the chub and endemic species of the L. cephalus complex met by secondary contacts. Our results show several patterns of introgression, from Leuciscus lepidus fully introgressed by chub mitochondrial DNA to Leuciscus borysthenicus where no introgression at all was detected. We assume that these hybridization events might constitute an important evolutionary process for the settlement of the chub in new environments in the Mediterranean area.

Southern Mexican minnows of the genus Notropis (Actinopterygii, Cyprinidae): genetic variation, phylogenetic relationships and biogeographical implications.

Allozyme analyses were performed to determine patterns of genetic variation and phylogenetic relationships within the genus Notropis in southern Mexico. Products of 28 gene loci were resolved in 24 geographic samples belonging to four putative species. These species represent the southern limit of the Cyprinidae in North America. Five loci were found monomorphic and 11 were diagnostic among species when the outgroup, Phoxinus phoxinus, was not taken into account. Four groups were identified, but these do not correspond perfectly to the four nominal species. Notropis sallaei was identified by allozyme analyses being the most basal among southern Mexican Notropis and the most genetically divergent. Notropis imeldae, was also diagnosed on the basis of allozymic variation, however a population sample representing the Balsas drainage was genetically divergent and is considered as an undescribed new species (N. n. sp.). No genetic differences were found between the samples of N. boucardi and N. moralesi. Therefore, we recommend that N. moralesi should be considered as a junior synonym of N. boucardi. A broad geographic sampling strategy was employed across all the distribution range of N. boucardi representing rivers from three different drainages, which cover both Mexican slopes. The slight divergence found among the headwater populations of N. boucardi permitted us to predict a model of paleohydrographic relationships of these three drainages. Two alternative hypotheses are postulated to explain the current distribution pattern of populations of N. boucardi in southern Mexico.

Spatial heterogeneity in the Mediterranean Biodiversity Hotspot affects barcoding accuracy of its freshwater fishes.

Incomplete knowledge of biodiversity remains a stumbling block for conservation planning and even occurs within globally important Biodiversity Hotspots (BH). Although technical advances have boosted the power of molecular biodiversity assessments, the link between DNA sequences and species and the analytics to discriminate entities remain crucial. Here, we present an analysis of the first DNA barcode library for the freshwater fish fauna of the Mediterranean BH (526 spp.), with virtually complete species coverage (498 spp., 98% extant species). In order to build an identification system supporting conservation, we compared species determination by taxonomists to multiple clustering analyses of DNA barcodes for 3165 specimens. The congruence of barcode clusters with morphological determination was strongly dependent on the method of cluster delineation, but was highest with the general mixed Yule-coalescent (GMYC) model-based approach (83% of all species recovered as GMYC entity). Overall, genetic morphological discontinuities suggest the existence of up to 64 previously unrecognized candidate species. We found reduced identification accuracy when using the entire DNA-barcode database, compared with analyses on databases for individual river catchments. This scale effect has important implications for barcoding assessments and suggests that fairly simple identification pipelines provide sufficient resolution in local applications. We calculated Evolutionarily Distinct and Globally Endangered scores in order to identify candidate species for conservation priority and argue that the evolutionary content of barcode data can be used to detect priority species for future IUCN assessments. We show that large-scale barcoding inventories of complex biotas are feasible and contribute directly to the evaluation of conservation priorities.

Permanent genetic resources added to Molecular Ecology Resources Database 1 December 2010-31 January 2011.

This article documents the addition of 238 microsatellite marker loci to the Molecular Ecology Resources Database. Loci were developed for the following species: Alytes dickhilleni, Arapaima gigas, Austropotamobius italicus, Blumeria graminis f. sp. tritici, Cobitis lutheri, Dendroctonus ponderosae, Glossina morsitans morsitans, Haplophilus subterraneus, Kirengeshoma palmata, Lysimachia japonica, Macrolophus pygmaeus, Microtus cabrerae, Mytilus galloprovincialis, Pallisentis (Neosentis) celatus, Pulmonaria officinalis, Salminus franciscanus, Thais chocolata and Zootoca vivipara. These loci were cross-tested on the following species: Acanthina monodon, Alytes cisternasii, Alytes maurus, Alytes muletensis, Alytes obstetricans almogavarii, Alytes obstetricans boscai, Alytes obstetricans obstetricans, Alytes obstetricans pertinax, Cambarellus montezumae, Cambarellus zempoalensis, Chorus giganteus, Cobitis tetralineata, Glossina fuscipes fuscipes, Glossina pallidipes, Lysimachia japonica var. japonica, Lysimachia japonica var. minutissima, Orconectes virilis, Pacifastacus leniusculus, Procambarus clarkii, Salminus brasiliensis and Salminus hilarii.

Insights on speciation patterns in the genus Iberochondrostoma (Cyprinidae): evidence from mitochondrial and nuclear data.

In this paper, the patterns of cladogenesis in the cyprinid fish genus Iberochondrostoma were analysed using a mitochondrial (cytochrome b) and a nuclear (beta-actin) gene fragment. The two genes yielded discordant results. While the cytochrome b gene yielded a fully dichotomous tree, where all species of the genus are monophyletic, the much slower beta-actin gene yielded star-like relationships. However, when information from both genes was considered together, the data suggested the persistence of a very large central unit from which at least two peripheral clades arose at different times. This pattern which is akin to peripatric speciation was shown to be compatible with the paleogeographical information available. It is suggested that combining the techniques of phylogeny and phylogeography and the use of multiple markers varying in their rate of evolution may enrich our understanding of speciation and evolution of clades beyond species level.