Skip to main content

Historical biogeography reveals new independent evolutionary lineages in the Pantosteus plebeius-nebuliferus species-group (Actinopterygii: Catostomidae)

  • The Correction to this article has been published in BMC Evolutionary Biology 2018 18:206



The Pantosteus plebeius-nebuliferus species-group is a group of freshwater fishes distributed in endo- and exorheic drainage basins in the Mexican Sierra Madre Occidental mountain range system and central North Mexico. The geological history of this region is considered an important factor in explaining the evolutionary history of low vagility animals like freshwaters fishes. The aim of this study was to examine the phylogenetic relationships and describe the evolutionary history of the species-group. We hypothesized that the genetic structure and distribution of the main clades of Pantosteus plebeius-nebuliferus are associated with the geological history of Northern Mexico. To this end, we obtained DNA sequences of mitochondrial and nuclear genes and performed phylogenetic and phylogeographic analyses. Divergence time estimation and ancestral area reconstruction were also carried out to propose a biogeographical hypothesis, and species boundaries within the species-group were also tested.


We identified four clades within the Pantosteus plebeius-nebuliferus species-group in both markers. Divergence ranged from 5.9% to 9.2% for cytb and 0.1% to 0.9% for GHI. We observed significant genetic structure and no shared haplotypes between clades. We estimated that the clades diverged during the last 5.1 Myr, with a biogeographic scenario suggesting eight vicariant and four dispersal events through the historic range of the species-group. We found that the best species-delimitation model is when four species are assumed, which correspond to the main clades. We identified nine evolutionary significance units (ESUs), pertinent to the conservation of the group, each representing populations present in distinct drainage basins.


The evolutionary history of the Pantosteus plebeius-nebuliferus species-group is characterized by vicariant post-dispersal processes, linked to geological changes in the Sierra Madre Occidental and central Northern Mexico since the Pliocene. This is congruent with biogeographic patterns described for other co-distributed fish species. We propose a new phylogenetic hypothesis for the species-group, clarifying the taxonomy of this evolutionarily complex group. Our results suggest that the species-group consists of at least four clades with independent evolutionary histories, two of which may represent new undescribed species. Our identification of ESUs provides a basis upon which conservation measures can be developed for the species-group.


The structure of current North American biodiversity is related to geological (tectonic and volcanic) and climatic changes that occurred mainly during the Neogene (ca. 33 Million years ago (Mya) and Quaternary (ca. past 2.5 Mya) [1, 2]. The Sierra Madre Occidental (SMOC) mountain range in Mexico is an area of high faunal and floristic endemism, and is considered an important biogeographic corridor [3]. The formation of the SMOC during the Oligocene (ca. 33-23 Mya) until the present [2, 4] is an important factor influencing the evolutionary history of several taxa, particularly organisms with low vagility such as freshwater fish species [3, 5,6,7,8,9,10,11,12].

Although the biodiversity of northwestern Mexico is mainly related to tecto-volcanic activity of the Tertiary, the SMOC is also considered a Pleistocenic refuge, and associated with the expansion and contraction of the distribution of numerous species [13]. The effects of climate changes on epicontinental waters has significantly affected the distribution of the fish fauna of this region [3]. The evolutionary history of the sucker fish family Catostomidae in Mexico may have been influenced by these processes, especially in the case of the Catostominae subfamily, the predominant fish group of the family in Mexico [14].

Within the Catostominae, two closely related species of the genus Pantosteus Cope & Yarrow, 1875; inhabit northern Mexico including the SMOC: Pantosteus plebeius (Baird & Girard, 1854) and Pantosteus nebuliferus (Garman, 1881) [14, 15]. Pantosteus nebuliferus has a restricted range, and is endemic to the endorheic Nazas and Aguanaval drainage basins. In contrast, Pantosteus plebeius is widespread and occurs in several drainage basins across the SMOC: the Mezquital, Piaxtla, Fuerte, and Yaqui basins; the Central Guzman hydrographic system (Santa Maria, Casas Grandes, and Del Carmen basins); the Upper Rio Grande basin (Conchos River in Mexico and Rio Grande River in USA); and the Mimbres basin in New Mexico, USA (Fig. 1) [16]. Pantosteus plebeius is also distributed in the Gila River of the Colorado basin, but the origin of this population is uncertain, some consider it to have been artificially introduced [17, 18], while others suggest this population is the result of stream capture from the Mimbres River [19].

Fig. 1

Distribution of Pantosteus plebeius-nebuliferus species-group and sampling localities. The colours correspond to the eight main basins and numbers to the localities: (1) Santa Gertrudis, (2) Arroyo Las Bayas, (3) Pino Suarez, (4) La Barranca, (5) Puente Mimbres, (6) El Cuarto, (7) Peñón Blanco, (8) Covadonga, (9) Arroyo Torreones, (10) El Peñasco, (11) El Olote, (12) Atotonilco, (13) Las Vegas, (14) Fuerte, (15) Oteros, (16) Bocoyna, (17) Conchos, (18) Ureyna, (19) South Fork Palomas Creek, (20) Escalariado, (21) Casas Grandes, (22) Ignacio Zaragoza, (23) Buenaventura and (24) Santa Clara

As with other Pantosteus species in Mexico, P. plebeius is believed to have colonized several drainage basins during the postglacial period [20]. Dispersal during anastomosis of rivers, followed by vicariance after river isolation, has been hypothesized as the main factor influencing genetic divergence among Pantosteus populations. It has been suggested that the populations of some drainage basins could represent one or more independent evolutionary lineages, such as those of the Mezquital, Piaxtla, Fuerte, and upper Rio Grande basin, each displaying significant genetic structure compared with Mimbres and Gila River populations [21, 22]. Given the complex evolutionary history of P. plebeius it has been suggested that this lineage may represent several distinct species, one or more of which may be synonymous with P. nebuliferus [23].

Campostoma ornatum (Girard, 1856) and Codoma ornata (Girard, 1856) are codistributed with the Pantosteus plebeius-nebuliferusc species-group throughout its range. Phylogeographic studies of these species have shown repeated isolation scenarios followed by merge or convergence in response to range expansion (reticulated biogeographic history) [3, 11]. This pattern was related to repeated events of dispersal and isolation associated with tecto-volcanic activities occurring in the SMOC since the early Pliocene, as well as to climatic fluctuations during glacial and interglacial periods. Population structure in Campostoma ornatum and Codoma ornata is reflected in high genetic differentiations among populations inhabiting the drainage basin in the northern SMOC, including the Guzman system, and the Río Grande, Yaqui, Mayo, Fuerte, and Nazas basins [3, 11].

The aim of this study was to 1) examine the phylogeography and reconstruct the evolutionary history of the Pantosteus plebeius-nebuliferus species-group in Mexico; 2) to explore the genetic diversity and structure of the species-group; and 3) determine population divergence times using nuclear and mitochondrial markers. We hypothesize that the genetic structure and distribution of the main clades of Pantosteus plebeius-nebuliferus are associated with the geological history of the SMOC and central Northern Mexico.


Sampling and dna isolation

We collected 87specimens with permission of relevant authorities in July 2013 from 20 localities in 8 drainage basins (Fig. 1 and Table 1) using electrofishing and hand nets. We included 6 samples from 4 localities from GenBank (GU937833, KJ441235 and KJ441237-KJ441240) for a total of 24 localities in 9 basins (localities 14, 18, 19 and 20 in Fig. 1).

Table 1 Sampled populations of the Pantosteus plebeius-nebuliferus species-group

Despite a considerable sampling effort, we were unable to obtain any specimens from the Aguanaval and Yaqui basins. We obtained fin clips and preserved them preserved in absolute ethanol and stored them at -75 °C. We fixed several specimens in formalin for identification and deposited in the Fish Collection of the Universidad Michoacana de San Nicolas de Hidalgo CPUM, Morelia, Mexico. We released the remaining captured fish alive.

We performed isolation of genomic DNA with BioSprint DNA Blood Kit QIAGEN according to the manufacturer’s instructions. We amplified the complete cytochrome b mitochondrial gene (cytb) using the primers GLuDG [24] and H16460 [25] and the nuclear 3rd intron of the growth hormone copy I (GHI), using the primers GHI3F and GHI3R [26].

We amplified DNA samples by Polymerase Chain Reaction (PCR), using 12 μl volume reactions with final concentrations of 0.2 μM of each primer, 0.25 mM of each dNTP, 1.5 M of MgCL2, and 1 U of Taq DNA Polymerase. The PCR procedure consisted of 2 min at 95 °C followed by 35 cycles of 45 s at 94 °C for denaturation, 1.5 min at 46 °C and 52.2 ° C for primer alignment for cytb and GHI respectively, 2 min at 72 °C for synthesis, and a final extension of 5 min at 72 °C. We quantified the PCR products by electrophoresis on 1.5% agarose gel and submitted to Macrogen Inc. (Korea) and htSEQ Inc (High-Throughput Sequencing, University of Washington, USA) for sequencing.

Genetic variation

We estimated the genetic diversity within drainage basins for the two gene markers based on nucleotide diversity (π), haplotype diversity (h), and the proportion of segregating sites (ΘS), using the software DNAsp 5.0 [27]. We obtained the genetic p-distances (DP) between populations using Mega 5.2 [28]. A bootstrapping process was implemented with 1000 repetitions.

We estimated the genetic structure at the interspecific and intraspecific level of the Pantosteus plebeius-nebuliferus species-group with fixation indices (ΦST) for both molecular markers. We used spatial analysis of molecular variance of geographically homogeneous K groups (SAMOVA) with SAMOVA 2.0 [29] and AMOVA (for K=1 to identify genetic structure without a priori information and the genetic variance explained within population) with Arlequin [30] to define groups of populations (basins) that were maximally differentiated without constraints of their geographic distribution. For both SAMOVA and AMOVA analysis we used 10 000 iterations from each of 100 random initial conditions. We tested for each K value from 2 to 8 in SAMOVA. We applied a Bonferroni correction [31] to each p-value obtained in the paired test of genetic differentiation.

Phylogenetic analysis and haplotype network reconstruction

We manually aligned DNA sequences in Mega 5.2 and examined them using chromatograms. We phased GHI sequences with point mutations using DNAsp and applied a test of recombination using a coalescent-based Bayesian method (10 000 replicates) in the same software.

We used two algorithms for phylogenetic reconstructions, using Catostomus catostomus as outgroup (GenBank Accession: AF454871 for cytb and GU937824 for GHI). We conducted independent Maximum Likelihood (ML) analyses with RAxMLGUI 1.3.1 [32, 33] for both genes, performing 10 000 bootstrap repetitions and using the evolutionary substitution model estimated with PartitionFinder [34]. These models were General Time-Reversible [35] + gamma (GTR+G) for the complete cytb and Tamura-Nei [36] + gamma (TrN+G) for GHI (see Additional file 1). We conducted Bayesian phylogenetic (BI) reconstruction with the software MrBayes 3.2.6 [37], using the above selected evolutionary substitution model and implementing two runs of four Markov Chain Monte Carlo (MCMC) processes with 7 million generations for cytb and 15 million generations for GHI, sampling every 100 generations in both cases. The difference in chain length was because GHI required more generations to obtain sufficient values of effective sample size (> 200). We evaluated the convergence of the log-likelihood (-InL) values of the two runs, with 10% of reconstructions discarded as burn-in, to construct the consensus tree (σ < 0.005). We used the posterior probabilities obtained based on a confidence limit of 95% (highest posterior density-HPD) to evaluate the support of nodes. We used a incongruence length difference test (Partition homogeneity test) [38] in PAUP* 4.0b10 [39] to evaluate the significance of conflict among data sets. We used 1000 resampling characters, and performed concatenated analyses of nuclear and mitochondrial data for ML and BI with the parameters previously mentioned. We constructed haplotype networks for both genes using PopART (available at htt:// and applying the Median-Joining method [40].

Divergence times estimation and ancestral area reconstruction

We estimated divergence times by a coalescent-based method using the software BEAST 1.8.0 [41], based on the described evolutionary substitution model for both genes (see Additional file 1). We used a relaxed molecular clock, using *Beast (including nuclear an mitochondrial data) with Yule speciation tree model. The clock was calibrated using the fossil record (in hard minimum bound for Lognormal distribution), which was located in the root of the Pantosteus plebeius-nebuliferus species-group. This was represented by fossils of Pantosteus asitus, dated from the Miocene-Pliocene (7.5-2.5 Mya) [42]. We conducted another molecular clock analysis using the mutation rate range for cytb estimated for teleosts of 0.76-2.2%/Myr with uniform distribution in BEAST [43,44,45,46].

We implemented three independent *Beast analyses, each with 70 000 000 generations, sampling every 1000 generations. After evaluating the posterior parameter values based on effective sample size, and according to the convergence data, using Tracer 1.6 [47], we discarded 10% of runs and combined the three analyses to construct a maximum clade credibility tree using BEAST modules (LogCombiner and TreeAnnotator, respectively).

We performed an ancestral area reconstruction with RASP 3.2, 20160719 [48], applying Statistical-Dispersal Vicariance Analysis (S-DIVA) [49, 50] and the Statistical Dispersal-Extintion-Cladogenesis model (S-DEC) [51]. In both cases we used the chronograms obtained in BEAST to resample every 3000 samples from the total set to obtain a subset of 70 000 trees. We considered nine areas in the analysis corresponding mainly to the drainage basins sampled The areas included were the basins: A, Upper North Rio Grande (Palomas); B, Casas Grandes; C, Santa Maria; D, Del Carmen; E, Fuerte; F, Upper South Rio Grande (Conchos); G, Piaxtla; H, Nazas; and I, Mezquital. We applied an unconstrained model allowing any combination of geographic range in adjacency matrix with a maximum of five adjacent hydrographic basins, configuring the dispersal rate with the same probability for each population (basin) of the Pantosteus plebeius-nebuliferus species-group.

Species delimitation analysis

We performed species delimitation tests using Bayesian Phylogenetics and Phylogeography software BPP 3.4 [52] to generate posterior probabilities for two species-delimitation models in the Pantosteus plebeius-nebuliferus species-group,. We used the reversible-jump Markov chain Monte Carlo (rjMCMC) methodology [53] to implement an species delimitation analysis A10 in BPP using a fixed guide tree (in this case a species tree) for a species-delimitation model assuming four species according with the main clades observed, which was obtained in the divergence time estimation as mentioned above. We provided an additional guide tree for a species-delimitation model assuming two species, where we assigned individuals to lineages defined by recognized species P. plebeius and P. nebuliferus according with the distribution described for these species [16]. We performed an analysis A11 for unguided species delimitation analysis to determine the best species-delimitation model compared to a priori models mentioned above under A10 analyses.

Including both molecular markers, we analyze three combinations of parameters for population size (Θ) and species divergence time (τ): In large population size with deep divergences (Θ α=1, β=10; τ α=1, β=10); small population size with low divergences (Θ α=2, β=2000; τ α=2, β=2000); and large population size with low divergences (Θ α=1, β=10; τ α=2, β=2000). We ran the analyses for 500 000 generations, sampling each five generations, with a final burn-in of the first 25 000.


We obtained a total of 92 DNA sequences of the cytb gene (1140bp) (Genbank accession numbers: MG203619-MG203705) and 83 of the 3rd intron of GHI (654bp) (Genbank accession numbers: MG203706-MG203869) for Pantosteus plebeius (n = 63 of cytb and n = 49 of GHI) and P. nebuliferus (n = 29 of cytb and n = 33 of GHI) from 22 localities in 8 drainage basins, representing most of the species distribution range [16]. We sequenced 87and 82 sequences respectively for cytb and GHI and obtained5 of cytb and 1 for GHI from GenBank (Fig. 1 and Table 1). We did not detect any significant recombination in the nuclear GHI sequences (p > 0.05).

Genetic diversity

Although the haplotype diversity in all basins was relative high (h = 0.829-0.981) for cytb and GHI, the nucleotide diversity and the proportion of segregating sites were low for both markers (π = 0.001-0.010 and ΘS = 0.0005-0.010) (Table 2). The populations with the lowest nucleotide diversity (π = 0.003) were those from Fuerte basin. The Casas Grandes and Nazas basin populations showed the highest number of haplotypes. However, the number of haplotypes was different in each marker by population as result of the difference in variation that each molecular marker has.

Table 2 Genetic diversity for cytb|GHI markers in Pantosteus plebeius-nebuliferus populations

Genetic distances and structure

The highest absolute genetic distances (AGD) for cytb found between clades were between the Mezquital (Clade IV) and Piaxtla/Nazas basins (Clade III) (DP = 9.2%) (Table 3), and higher AGD than 5% between all clades was observed. Based on GHI the AGD highest distances (DP = 0.9%) were between the Mezquital and Fuerte and Del Carmen populations, and between the Piaxtla and Del Carmen and Fuerte populations. Between all clades AGD distances were lower than one percent. The lowest AGD were observed between populations of Casas Grandes and Santa Maria populations for cytb (DP = 1.3%) and between the Palomas and Casas Grandes and Santa Maria for GHI (DP = 0.1%), both comparisons within the Clade I (Guzman hydrographic system).

Table 3 Absolute pairwise un-corrected p-distances DP for cytb (under the diagonal) and GHI (above the diagonal using pairwise-deletion) of populations of Pantosteus plebeius-nebuliferus species-group

AMOVA revealed significant genetic structures among populations (ΦST = 0.9028 for cytb; ΦST = 0.6904 for GHI; p < 0.05) (Table 4). The highest genetic structure with SAMOVA was found when K = 8 (ΦCT = 0.8784 for cytb; ΦCT = 0.7460 for GHI, p < 0.05), and each basin was clustered in an independent group, with the exception of one population of the current Upper South Rio Grande, the Conchos population, which was clustered with the Fuerte population. In the case of GHI, no significant genetic structure was observed when groups of populations (basins) recognized for P. plebeius and P. nebuliferus were compared (ΦCT = 0.1707, p > 0.05).

Table 4 Spatial Analysis of Molecular Variance (SAMOVA) of the Pantosteus plebeius-nebuliferus species-group baes on cytb|GHI. Only the results for K = 2, 4, and 8 are shown, as they are meaningful in phylogenetic, systematic, and biogeographic terms

High ΦST and significant pairwise values were observed among basins in both molecular markers, with ΦST ranging form 0.461 to 0.968 for cytb and form 0.050 to 0.907 for GHI (p < 0.5) (Table 5). The highest value was observed between the populations of Piaxtla and Fuerte for cytb, while between Piaxtla and Del Carmen populations for GHI. The lowest but significant ΦST was between Piaxtla and Nazas for cytb and between Del Carmen and Santa Maria for GHI. In contrast to cytb ΦST values, no significant ΦST between populations of Santa Maria and Casas Grandes were observed for GHI.

Table 5 Genetic differentiation using pairwise ΦST for cytb (under the diagonal) and GHI (above the diagonal) among basins

Phylogenetic relationships and haplotype network analysis

The incongruence length difference test did not show significant differences (p < 0.05), indicating that both molecular markers presented the same phylogenetic signal. Bayesian and Maximum Likelihood phylogenetic analyses for the cytb and GHI sequences (see Additional files 2 and 3) and concatenated analysis (Fig. 2) recovered the same topology, with four highly supported clades (main clades). Clade I included samples from the Guzman hydrographic system (Casas Grandes, Santa Maria, and Del Carmen basins) and the Upper North Rio Grande system (Palomas). Clade II comprised the populations of the Fuerte basin and Conchos River (Upper South Rio Grande basin). Clade III, corresponding to the Nazas and Piaxtla basin populations and Clade IV consisted of the population from the Mezquital basin.

Fig. 2

Phylogenetic reconstruction Pantosteus plebeius-nebuliferus species-group based on the concatenated data of cytb gene and GHI region. The numerical values in nodes represent the posterior probabilities and the bootstrap values for Bayesian Inference and Maximum Likelihood, respectively

The same reciprocally monophyletic haplogroups were recovered in the haplotype networks reconstruction for both molecular markers, corresponding to the four main clades obtained by phylogenetic analysis, but differ in the relationships. In cytb, the Mezquital population (Clade IV) is related (72 mutations steps or MS) to the Guzman System/Palomas populations (Clade I), whereas the Nazas/Piaxtla populations (Clade III) are more closely related to the Conchos/Fuerte populations (72 MS) (Clade II) than to the Mezquital population (Fig. 3a and Fig. 3b). The number of MS separating the population of Fuerte from the Conchos River; and the population of Nazas from Piaxtla River, was 3 MS in both cases. We did not include as MS the median vectors or hypothetical haplotypes implied by the haplotype network reconstruction.

Fig. 3

Haplotype network for Pantosteus plebeius-nebuliferus species-group based on cytb (a) and GHI (b). The numbers inside the circles represent the number of mutations steps between haplotypes. The grey circles represent median vector or hypothesised sequences which is required to connect existing sequences within the network

In case of GHI, the range of MS between clades was 1-2 MS, geographic correspondence was observed, and shared haplotypes were identified between populations that are closely related phylogenetically (Casas Grandes/Del Carmen/Santa Maria/Palomas, Fuerte/Conchos and Piaxtla/Nazas) (Fig. 3b).

Divergence times and ancestral area reconstruction

The global likelihood scores with the highest probabilities were obtained using the S-DEC approach (Fig. 4). However total congruence in ancestral areas predicted was observed between both analyses. The scenario showed by the two models (S-DIVA and S-DEC) inferred an ancestral area (ABCDF) of the most recent common ancestor (MRCA) of the Pantosteus plebeius-nebuliferus species-group formed by the Guzman system, and the Upper South and North Rio Grande (Conchos and Palomas, respectively). From this first vicariant event dated at ca. 5.1-3.9 Mya (HPD ≥ 95%) to the upper Pliocene, a succession of 7 vicariant and 4 dispersal events were identified during the historical distribution of the species (Fig. 4).

Fig. 4

Divergence times estimation (based on species tree), S-DIVA and S-DEC analyses of the Pantosteus plebeius-nebuliferus species-group. Red labels represent the divergence times in millions years ago. The bold names represent the ancestral area estimated. Dispersal events and vicariant events are represented with - and | symbols, respectively. The * represent the value for posterior probability of 100% for each clade tested in the most sensitive species-delimitation model (four species) with values of Θ α=1, β=10; τ α=1, β=10 in the species delimitation analysis A10 and A11

Species limits

The species-delimitation model assuming four species was strongly sensitive to low parameter of Θ and τ, showing high posterior probability for the nodes in A10 and was the species-delimitation model in A11 analysis with the highest probability (P=1 S=4) (Fig. 4). The four species supported correspond to the main clades identified in this work for Pantosteus plebeius-nebuliferus (Fig. 2). The model of species-delimitation reducing the number of species to two (P. plebeius and P. nebuliferus) showed low posterior probabilities for the three combinations of parameters of Θ and τ (results not showed).


The main forces shaping cladogenetic events in the Pantosteus plebeius-nebuliferus species-group are consistent with geological processes occurring since the Pliocene associated with the formation and evolution of the SMOC, central North Mexico and the tectonic activity of the Rio Grande Rift. Our results revealed that the species-group comprises more than two independent evolutionary lineages with significant genetic structures between them in nuclear and mitochondrial markers, indicating a long history of isolation. Within clades mitochondrial information showed high divergence and several mutational steps between clades, unlike the nuclear marker that showed shared haplotypes between some basins; and low number of mutational steps between clades, reflecting the combined effects of low variation and incomplete lineage sorting on the nuclear loci [54]. As has occurred with the genus Moxotoma [12] and Campostoma ornatum [3], relevant vicariant events and probably river captures may have played an important role in range expansion shaping the spatial distribution of genetic variation in some drainage basins. Genetic variation within drainage basins seems to have been influenced by ecological characteristics of the species and climatic fluctuations since the Pleistocene. Further studies should seek to identify the roles of historical, intrinsic, and anthropogenic influences on genetic differentiation within basins [55].

Historical biogeography of pantosteus plebeius-nebuliferus

Pliocene events

We dated the most recent common ancestor of all clades of the Pantosteus plebeius-nebuliferus species-group to the Pliocene (5.1-3.9 Mya, HPD ≥ 95%) (Fig. 5), in an ancestral area comprised by the proto-Guzman system, proto-Upper South Rio Grande (Conchos) and proto-Upper North Rio Grande (Palomas). We identified a range expansion of the species-group to the Fuerte basin, which probably resulted from a basin catchment, allowing faunal interchange as proposed for Gila pulchra [10], and Codoma ornata [11]. These events occurred during the tecto-volcanic episodes in SMOC evolution, including repeated volcanic activity called “alkaline basalt events” [4, 56, 57]. Another dispersal event from the ancestral areas to Nazas was proposed. These events could have been promoted by alkaline basalts events during the Mapimi Bolson formation (Fig. 4 and Fig. 5) [12, 58, 59], which is partially congruent with patterns established for Codoma ornata [11] and Campostoma ornatum [7].

Fig. 5

Biogeographical scenario for the Pantosteus plebeius-nebuliferus species-group in the North Mexico

These dispersal events were followed by a vicariant event that isolated the ancestor of the Proto Guzman system/Palomas (Clade I) ca. 4.5 Mya, associated with the volcanism in the Babicora-Bustillos sector, disrupting the connection between the proto-Fuerte/Conchos and the proto-Guzman systems/Palomas by the periodic accumulation of sediments in the Conchos Valley [12]. This ancient isolation event is reflected in the high genetic divergence between clades I and II in cytb (DP = 6.3-7.0%, ΦST = 0.911-0.938 p < 0.05; Table 3) and in the well-structured haplogroups for both markers (Fig. 3). The split between clade I with respect to the clades II, III, and IV of the P. plebeius-nebuliferus species-group, is geographically, genetically and in some cases temporarily consistent with other co-distributed species complexes of fishes such as Cyprinella spp. [60], Campostoma ornatum [3, 61], Cyprinodon spp. [6, 62], Gila spp. [10], Codoma ornata [11] and Moxostoma cf. austrinum from Conchos [12]. This finding also corroborates the divergence time of the Rio Grande Pantosteus spp. population with respect to the Nazas population in a previous work [63].

The Mapimi Bolson desertification and the alkaline basalt activity at the central part of the SMOC, was likely responsible for the disconnection that promoted the vicariant event of the ancestor of the proto-Nazas/Piaxtla (Clade III)/Mezquital basin populations (Clade IV) from the proto-Fuerte/Conchos populations in the Pliocene, an event that is consistent with results found for Codoma ornata results [11]. Two range expansion events of the species-group from the Nazas to Piaxtla and Mezquital basins following the vicariant event, were found in the same period, which could be the result of river catchments associated with alkaline basalt events in the SMOC [4, 12, 57], an event that is consistent with results found in other aquatic fauna, such as Campostoma ornatum [3] and semi-aquatic snake species of the genus Thamnophis [64].

Due to the physiographic complexity and the occurrence of several overlapping geologic events in the region of Nazas, Piaxtla and Mezquital, is difficult to identify specific geologic events involved in the cladogenesis of clades III and IV occurred ca. 2.6 Mya. The isolation of the proto-Nazas/Piaxtla from the Mezquital populations possibly had its source in the alkaline volcanic activity occurring in the region during the Pliocene, especially south of Laguna Santiaguillo in the Nazas basin, in the southern region of the SMOC [4] (Fig. 4) separating the tributary El Tunal (Mezquital Basin) from the Nazas basin occurring in the Late Miocene and Pliocene [2, 65, 66].

Pleistocene isolation

The origin and evolution of the three basins within the Guzman system date from the Late Cretaceous and Cenozoic to the Late Pleistocene, and its complexity could explain the species diversity and genetic divergence of fishes in the region. During the Early Pleistocene, the activity of the Rio Grande Rift and the glacial stage (Kansan, Middle Pleistocene) [67] were likely involved in the fragmentation of the Cabeza de Vaca paleolake and the formation of the Palomas pluvial paleolake, a relict of the Cabeza de Vaca that was fragmented into several endorheic rivers, lagoons, and springs. Among these, the Casas Grandes and Upper Rio Grande basins were firstly isolated from Del Carmen and Santa Maria (ca. 1.8 Mya) by the formation of the Los Muertos Bolson [60, 68,69,70] (Fig. 5). The isolation of the Upper North Rio Grande, including its flow to the Gulf of Mexico and capture of the Conchos River, was dated to the Early Pleistocene [71,72,73,74], but the tectonic subsidence in the Rio Grande Rift and desertification of Los Muertos Bolson has continued during the Middle Pleistocene and the present [75,76,77] allowing the fragmentation of the Guzman system. This paleohydrological pattern is reflected by our data, since we date the isolation of Palomas from Casas Grandes as ca. 1.1 Mya and 1.2 Mya for the isolation of Del Carmen from Santa Maria (Fig. 4). This is also consistent with previous studies that found low genetic divergences between populations of Pantosteus plebeius from the Mimbres River (type locality of the species) and P. plebeius from the Palomas River (Upper North Rio Grande) [22].

The three basins within the Guzman system were isolated and reconnected several times during the Pre-Illinoian, Illinoian, and Pre-Holocene periods [78,79,80]. Remixing gene pools and low genetic structures also were observed in other co-ocurring fishes along the Guzman basin, including Gila spp. [10, 81] and Cyprinodon spp. [6, 62]. However, the high interbasin genetic distances and the lack of shared haplotypes for cytb is evidence of the ancestral isolation that occurred in Del Carmen and Casas Grandes populations of Campostoma ornatum (DP = 2.9%; ΦST = 0.935) [3], and the significant genetic differentiation (FST = 0.380) in Cyprinella formosa [60, 69, 82,83,84] which represents an isolation event that also occurred in the Pleistocene. In contrast, the low genetic differentiation and haplotypes shared among basins for GHI are evidence for the low variation and of the incomplete lineage sorting or a later genetic mixing between these basins. More studies are necessary to explain the configuration of the ichthyofauna of the Guzman system.

Finally, the vicariant events involved in the isolation of the Nazas from the Piaxtla and the Conchos from the Fuerte populations, were estimated to have occurred 1.4 Mya and 0.6 Mya in the Pleistocene (Fig. 4), respectively, and are associated with the most recent volcanic eruptions and tectonic movements in the SMOC [85]. This supports the results obtained for Campostoma ornatum [3], Gila pulchra (Conchos River & Fuerte basin) [10], and Codoma ornata (upper Conchos River & Fuerte basin) [11] that show low genetic differentiation as well as shared haplotypes among basins as we observed in the nuclear marker (Fig. 3).

Within-basin differentiation

A high number of mutational steps and genetic intrabasin distances where estimated for cytb in the Pantosteus plebeius-nebuliferus species-group. Surprisingly, moderate between-haplotype MS values (MS = 4-7) were found within the Mezquital, Nazas, Fuerte, Casas Grandes, and Santa Maria basin populations than between the Nazas and Piaxtla (MS = 3) and Conchos and Fuerte populations (MS = 3), reflecting considerable genetic differentiation in cytb (Fig. 3 and Tables 345). The intrapopulation genetic structure in those drainage basins could be explained in part by life history traits. Migration episodes associated with environmental fluctuations on spawning sites and food availability have been extensively described in other species of Catostomidae [86,87,88,89,90,91,92,93]. Accordingly, the historic climate fluctuations in the area [94,95,96] could promote similar fragmentation episodes in P. plebeius-nebuliferus populations, which could be responsible for the isolation of migrants in separate regions within the basins.

Taxonomic implications and conservation considerations

The Pantosteus plebeius-nebuliferus species-group comprises a set of distinct genetic groups, suggesting that it is composed of at least four clades with unique and independent evolutionary histories. We found cladogenesis to be associated with the high tecto-volcanic activity and climate fluctuations from the Pliocene to the present in North Mexico. This is supported by the phylogenetic relationships, species tree, species-delimitation test including both molecular markers. Also is supported by the genetic divergences observed between clades for cytb (see Additional file 2), which are higher than those reported among sibling vertebrate groups for cytb (~2%) [97,98,99,100] and between sister species within the family [98]. In the case of GHI, the genetic distances estimated between species of Catostomidae are 3.3% [101] taking into account different genera of the family, and which were higher than estimated in this work (Table 3, DP<1%). However, although the nuclear GHI showed a low number of mutational steps and genetic distances (Fig. 3 and Table 3), we detected significant structure (Table 4) and recovered the four clades in the phylogenetic inference (see Additional file 3) supporting the mitochondrial results. The genetic structure and divergences in the species-group were similar to those observed in other fish taxa that partially or completely occur in basins of northwest Mexico and the southwestern USA, such as Cyprinodon spp. [6, 62], Gila spp. [10, 81], Cyprinella spp. [60, 69, 84], or even among populations of other considered species complexes as Campostoma spp. [3, 61], and Moxotoma spp. [12]. Thus, we recommend an integrative taxonomic revision of the Pantosteus plebeius-nebuliferus to recognize the four independent evolutionary lineages as different species.

The population of the Piaxtla basin was recognized as Pantosteus plebeius [16] and suggested as an independent evolutionary lineage [20], our results show that this population is closely related to P. nebuliferus, as was also demonstrated by [63]. Low genetic divergences between populations of P. plebeius from the type locality Mimbres River and from the Upper North Rio Grande have been observed by other authors [22, 102]. Our results show that Paloma River is closely related with Guzman system populations, thus we suggest that Clade I must be considered populations of the recognized species P. plebeius. The Clade II (Fuerte basin and Conchos River populations) and the Clade IV (Mezquital River population) do not represent populations of any species recognized of the species-group studied, thus we suggest that the two clades could represent two independent evolutionary lineages that must be described and recognized as new species of the genus Pantosteus.

The identification of the Evolutionary Significant Units (ESU) is of particular interest especially in species considered threatened [103]. Such is the case of Pantosteus plebeius-nebuliferus species group, which is listed as in danger of extinction in the Norm-059-SEMARNAT-2010 [104]. We identified the group of populations whitin independent evolutionary lineages, which warrant separate management or priority for conservation. Based on the high genetic divergences; number of mutation steps; phylogenetic relationships; shared haplotypes; reproductive isolation and the subsequent absent exchangeability of populations, we suggest nine ESUs corresponding to each basin where the P. plebeius-nebuliferus species-group inhabit: (1) Upper North Grande River, (2) Casas Grandes, (3) Santa Maria, (4) Del Carmen, (5) Fuerte, (6) Upper South Grande River basin (Conchos River), (7) Piaxtla, (8) Nazas and (9) Mezquital. The definition of these nine ESUs should aid in establishing conservation measures of the species-group and of the evolutionary lineages of which they are configured.


Our study supports at least four independent evolutionary lineages in the Pantosteus plebeius-nebuliferus species-group. Their biogeographic history is linked to geological events occurring since the Pliocene, associated with the formation and evolution of the SMOC, central North Mexico and the tectonic activity of the Rio Grande Rift. Other freshwater fishes already studied and occurring partially or completely in the same range of northwest Mexico and the southwestern USA have similar genetic structures and divergences. We recommend an integrative taxonomic revision of the Pantosteus plebeius-nebuliferus species-group and propose nine ESUs for conservation purposes which are restricted to individual basins where the species-group is distributed. As an endangered species, we hope our survey contributes to the management and preservation of these fishes.

Change history

  • 27 December 2018

    Following publication of the original article



Absolute Genetic Distances


Analysis of Molecular Variance


Bayesian Inference

ca :



Fish Collection of the Universidad Michoacana de San Nicolás de Hidalgo

cytb :

Cytochrome b

DP :

Genetic p-distances


Evolutionary Significance Unit


3rd Intron of Growth Hormone copy I


General Time-Reversible + Gamma


Haplotype diversity


Number of haplotype


Highest Posterior Density




Markov Chain Monte Carlo


Maximum Likelihood


Most Common Ancestor


Mutation Steps


Million years ago


Million years


Polymerase Chain Reaction


Spatial Analysis of Molecular Variance


Statistical Dispersal-Extintion-Cladogenesis


Statistical-Dispersal Vicariance Analysis


Sierra Madre Occidental


Tamura-Nei + Gamma


Upper North River


United States of America


Upper South River

Θ :

Population size

ΘS :

Proportion of segregating sites


Nucleotide diversity


Species divergence time


Pairwise fixation indices


  1. 1.

    Riddle BR, Hafner DJ. A step-wise approach to integrating phyloeographic and phylogenetic biogeographic perspectives on the history of a core North American warm deserts biota. J Arid Environ. 2006;66:435–61.

    Google Scholar 

  2. 2.

    Aguirre-Díaz GJ, Labarthe-Hernández G, Tristán-González M, Nieto-Obregón J, Gutiérrez-Palomares I. The Ignimbrite Flare-Up and Graben Calderas of the Sierra Madre Occidental, Mexico. In: Martí JGJ, editor. Caldera Volcanism Analysis, Modelling and Response, vol. 10. Oxford: Dev Volcano, Elsevier; 2008. p. 143–80.

    Google Scholar 

  3. 3.

    Domínguez-Domínguez O, Vila M, Pérez-Rodríguez R, Remón N, Doadrio I. Complex evolutionary history of the Mexican stoneroller Campostoma ornatum Girard, 1856 (Actinopterygii: Cyprinidae). BMC Evol Biol. 2011;11:153.

    PubMed  PubMed Central  Google Scholar 

  4. 4.

    Ferrari L, Valencia-Moreno M, Bryan S. Magmatism and tectonics of the Sierra Madre Occidental and its relation with the evolution of the western margin of North America. Geol S Am S. 2007;422:1–39.

    Google Scholar 

  5. 5.

    Rinne JN, Hughes RM, Calamusso B. Historical changes in large river fish assemblages of the Americas. Bethesda: Am Fish S S; 2005. p. 45.

    Google Scholar 

  6. 6.

    Echelle AA, Carson EW, Echelle AF, Van Den Bussche RA, Dowling TE, Meyer A. Historical biogeography of the New-World Pupfish genus Cyprinodon (Teleostei: Cyprinodontidae). Copeia. 2005;2:320–39.

    Google Scholar 

  7. 7.

    Domínguez-Domínguez O, Alda F, Pérez-Ponce de León G, García-Garitagoitia JL, Doadrio I. Evolutionary history of the endangered fish Zoogoneticus quitzoensis (Bean, 1898) (Cyprinodontiformes: Goodeidae) using a sequential approach to phylogeography based on mitochondrial and nuclear DNA data. BMC Evol Biol. 2008;8:161.

    PubMed  PubMed Central  Google Scholar 

  8. 8.

    Domínguez-Domínguez O, Pérez-Ponce de León G. La mesa central de México es una provincia biogeográfica? Análisis descriptivo basado en componentes bióticos dulceacuícolas. Rev Mex Biodivers. 2009;80:835–52.

    Google Scholar 

  9. 9.

    Schönhuth S, Hills DM, Neely DA, Lozano-Vilano L, Perdices A, Mayden RL. Phylogeny, diversity, and species delimitation of the North American Round-Nosed Minnows (Teleostei: Dionda), as inferred from mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol. 2012;62:427–46.

    PubMed  Google Scholar 

  10. 10.

    Schönhuth S, Perdices A, Lozano-Vilano L, García-de-León FJ, Espinosa H, Mayden RL. Phylonetic relationship of North American western chubs of the genus Gila (Cyprinidae, Teleostei), with emphasis on southern species. Mol Phylogenet Evol. 2014;70:210–30.

    PubMed  Google Scholar 

  11. 11.

    Schönhuth S, Lozano-Vilano L, Perdices A, Espinosa H, Mayden RL. Phylogeny, genetic diversity and phylogeography of the genus Codoma (Teleostei, Cyprinidae). Zool Scr. 2015;44:11–28.

    Google Scholar 

  12. 12.

    Pérez-Rodríguez R, Domínguez-Domínguez O, Mar-Silva AF, Doadrio I, Pérez-Ponce de León G. The historical biogeography of the southern group of the sucker genus Moxostoma (Teleostei: Catostomidae) and the colonization of central Mexico. Zool J Linn Soc. 2016;177:633–47.

    Google Scholar 

  13. 13.

    Ceballos G, Arroyo-Cabrales J, Ponce E. Effects of Pleistocene environmental changes on the distribution and community structure of mammalian fauna of Mexico. Quaternary Res. 2010;73:464–73.

    Google Scholar 

  14. 14.

    Wei-Jen C, Mayden RL. Phylogeny of suckers (Teleostei: Cypriniformes: Catostomidae): further evidence of relationships provided by the single-copy nuclear gene IRBP2. Zootaxa. 2012;3586:195–210.

    Google Scholar 

  15. 15.

    Doosey MH, Bart HL Jr, Saitoh K, Miya M. Phylogenetic relationships of Catostomid fishes (Actinopterygii: Cypriniformes) based on mitochondrial ND4/ND5 gene sequences. Mol Phylogenet Evol. 2010;54:1028–34.

    PubMed  Google Scholar 

  16. 16.

    Miller RR, Minckley WL, Norris SM. Freshwater Fishes of México. Chicago: University Of Chicago Press; 2005.

    Google Scholar 

  17. 17.

    Buth DG, Crabtree CB. A comparison of the genetic characteristics of Zuni sucker (Catostomus discobolus yarrow), the bluehead sucker (Catostomus discobolus), and the Rio Grande sucker (Catostomus plebeius). Santa Fe: Report of the New Mexico Department of Game and fish; 1995. p. 34.

    Google Scholar 

  18. 18.

    Page LM, Burr BM. Peterson field guide to freshwater fishes of North America north of Mexico. 2nd ed. Boston: Houghton Mifflin Harcourt; 2011.

    Google Scholar 

  19. 19.

    Sublette JE, Hatch MD, Subblete M. The fishes of New Mexico. Albuquerque: University of New Mexico Press; 1990.

    Google Scholar 

  20. 20.

    Ferris SD, Buth DG, Whitt GS. Substantial genetic differentiation among populations of Catostomus plebeius. Copeia. 1982;2:444–9.

    Google Scholar 

  21. 21.

    Crabtree CB, Buth DG. Biochemical systematics of the catostomid genus Catostomus: assessment of C. clarki, C. plebeius, and C. discobolus including the Zuni sucker, C. d. yarrow. Copeia. 1987;4:843–54.

    Google Scholar 

  22. 22.

    McPhee MV, Osborne MJ, Turner TF. Genetic diversity, population structure, and demographic history of the Rio Grande Sucker, Catostomus (Pantosteus) plebeius, in New Mexico. Copeia. 2008;1:191–9.

    Google Scholar 

  23. 23.

    Smith GR. Phylogeny and biogeography of the Catostomidae, freshwater fishes of North America and Asia. In: Mayden RL, editor. Systematics, historical ecology and North American freshwater fishes. Stanford: Stanford University Press; 1992. p. 778–826.

    Google Scholar 

  24. 24.

    Palumbi SR, Martin A, Romano S, McMillan MO, Stice L, Grabowski G. Simple fool’s guide to PCR. Honolulu: University of Hawaii; 1991.

    Google Scholar 

  25. 25.

    Perdices A, Bermingham E, Montilla A, Doadrio I. Evolutionary history of the genus Rhamdia (Teleostei: Pimelodidae) in Central America. Mol Phylogenet Evol. 2002;25:172–89.

    CAS  PubMed  Google Scholar 

  26. 26.

    Clements MD. Molecular systematics and biogeography of the Moxostomatini (Cypriniformes:Catostomidae), with special reference to taxa in Texas and Mexico. Ph.D. Dissertation. New Orleans: Tulane University; 2008.

    Google Scholar 

  27. 27.

    Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.

    CAS  Google Scholar 

  28. 28.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Dupanloup I, Schneider S, Excoffier L. A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002;11:2571–81.

    CAS  PubMed  Google Scholar 

  30. 30.

    Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.

    Google Scholar 

  31. 31.

    Rice WR. Analyzing tables of statistical test. Evolution. 1989;43:223–5.

    PubMed  Google Scholar 

  32. 32.

    Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.

    CAS  PubMed  Google Scholar 

  33. 33.

    Silvestro D, Michalak I. RAxMLGUI: a graphical front-end for RAxML. Org Divers Evol. 2012;12:335–7.

    Google Scholar 

  34. 34.

    Lanfear R, Calcott B, Ho SY, Guindon S. PartitionFinder: Combined selection of partitioning schemes and substitutions models for phylogenetic analyses. Mol Biol Evol. 2012;29:1695–701.

    CAS  PubMed  Google Scholar 

  35. 35.

    Tavaré S. Some probabilistic and statistical problems in the analysis of DNA sequences. In Lectures on mathematics in the life sciences.1986;17:57-86.

  36. 36.

    Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10:512–26.

    CAS  PubMed  Google Scholar 

  37. 37.

    Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.

    CAS  PubMed  Google Scholar 

  38. 38.

    Farris JD, Källersjö M, Kluge AG, Bult C. Constructing a significance test for incongruence. Syst Biol. 1995;44:570–2.

    Google Scholar 

  39. 39.

    Swofford DL. PAUP*: Phylogenetic analysis using parsimony (*and other methods). Version 4. Sunderland: Sinauer Associates; 2003.

    Google Scholar 

  40. 40.

    Bandelt HJ, Forster P, Röhl A. Median-Joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 2000;16:37–48.

    Google Scholar 

  41. 41.

    Drummond AJ, Ho SYW, Rawlence N, Rambaut A. A rough guide to BEAST 1.4. United Kingdom: University of Edinburgh; 2007.

    Google Scholar 

  42. 42.

    Smith GR, Stewart JD, Carpenter NE. Fossil and recent mountain suckers, Pantosteus, and significance of introgression in catostomin fishes of Western United States. Occas Pap Mus Zool Univ Mich. 2013;743:1–59.

    Google Scholar 

  43. 43.

    Zardoya R, Doadrio I. Molecular evidence on the evolutionary and biogeographical patterns of European Cyprinids. J Mol Evol. 1999;49:227–37.

    CAS  PubMed  Google Scholar 

  44. 44.

    Machordom A, Doadrio I. Evidence of a Cenozoic Betic-Kabilian connection based on freshwater fish phylogeography (Luciobarbus, Cyprinidae). Mol Phylogenet Evol. 2001;18:252–63.

    CAS  PubMed  Google Scholar 

  45. 45.

    Perdices AA, Carmona JA, Fernández-Delgado C, Doadrio I. Nuclear and mitochondrial data reveal high genetic divergence among Atlantic and Mediterranean populations of the Iberian killifish Aphanius iberus (Teleostei: Cyprinodontidae). Heredity. 2001;87:314–24.

    CAS  PubMed  Google Scholar 

  46. 46.

    Near TJ, Benard MF. Rapid allopatric speciation in logperch darters (Percidae: Percina). Evolution. 2004;58:2798–808.

    PubMed  Google Scholar 

  47. 47.

    Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6. 2014. Available at Accessed Feb 2017.

  48. 48.

    Yu Y, Harrid AJ, Blair C, He XJ. RASP (Reconstruct Ancestral State in Phylogenies): A tool for historical biogeography. Mol Phylogenet Evol. 2015;87:46–9.

    PubMed  Google Scholar 

  49. 49.

    Yu Y, Harris AJ, He XJ. S-DIVA (Statistical Disperal-Vicariance Analysis): A tool for inferring biogeographic histories. Mol Phylogenet Evol. 2010;56:848–50.

    PubMed  Google Scholar 

  50. 50.

    Nylander JA, Olsson U, Alström P, Sanmartín I. Accounting for phylogenetic uncertainty in biogeography: a Bayesian approach to dispersal-vicariance analysis of the thrushes (Aves: Turdus). Syst Biol. 2008;57:257–68.

    PubMed  Google Scholar 

  51. 51.

    Ree RH, Smith SA. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst Biol. 2008;57:4–14.

    PubMed  Google Scholar 

  52. 52.

    Yang ZH. The BPP program for species tree estimation and species delimitation. Curr Zool. 2015;61:854–65.

    Google Scholar 

  53. 53.

    Yang ZH, Rannala B. Bayesian species delimitation using multilocus sequence data. Proc Natl Acad Sci USA. 2010;107:9254–69.

    Google Scholar 

  54. 54.

    Egge JJD, Nicholson PW, Stark AW. (Siluriformes: Ictaluridae), a Mississipi Embayment endemic: evidence for a cryptic lineage in the Hatchie River. J Fish Biol. 2015;86:493–526.

    CAS  PubMed  Google Scholar 

  55. 55.

    Husemann M, Ray JW, King RS, Hooser EA, Danley PD. Comparative biogeography reveals differences in population genetic structure of five species of stream fishes. Biol J Linn Soc. 2012;107:867–86.

    Google Scholar 

  56. 56.

    Henry CD, Aranda-Gómez JJ. Plate interactions control middle-late Miocene proto-Gulf and Basin and Tange extension in the southern Basin and Range. Tectonophysics. 2000;318:1–16.

    Google Scholar 

  57. 57.

    Aranda-Gómez JJ, Luhr JF, Housh TB, Valdez-Moreno G, Chávez-Cabello G. El vulcanismo tipo intraplaca del Cenozoico tardío en el centro y norte de México: una revisión. Bol Soc Geol Mex. 2005;57(3):187–225.

    Google Scholar 

  58. 58.

    Smith ML, Miller RR. The evolution of the Rio Grande basin as inferred from its fish fauna. In: Hocutt CH, Wiley EO, editors. The zoogeography of North American freshwater fishes. New York: Wiley; 1986. p. 457–85.

    Google Scholar 

  59. 59.

    Mayden RL, Matson RH, Hillis DM. Speciation in the North American genus Dionda (Teleostei: Cypriniformes). In: Mayden RL, editor. Systematics, historical ecology and North American freshwater fishes. Stanford: Stanford University Press; 1992. p. 710–46.

    Google Scholar 

  60. 60.

    Schönhuth S, Mayden RL. Phylogenetic relationships in the genus Cyprinella (Actinopterygii: Cyprinidae) based on mitochondrial and nuclear gene sequences. Mol Phylogenet Evol. 2010;55:77–98.

    PubMed  Google Scholar 

  61. 61.

    Schönhuth S, Blum MJ, Lozano-Vilano L, Neely DA, Varela-Romero A, Espinosa H, Perdices A, Mayden RL. Inter-basin exchange and repeated headwater capture across the Sierra Madre Occidental inferred from the phylogeography of Mexican stonerollers. J Biogeogr. 2011;38:1406–21.

    Google Scholar 

  62. 62.

    Echelle AA. The western North American pupfish clade (Cyprinodontidae: Cyprinodon): mitochondrial DNA divergence times and drainage history. Geol S Am S. 2008;439:27–38.

    Google Scholar 

  63. 63.

    Unmack PJ, Dowling TE, Laitinen NJ, Secor CL, Mayden RL, Sciozawa DK, Smith GR. Influence of introgression and geological processes on phylogenetic relationships of western North American Mountain Suckers (Pantosteus, Catostomidae). Plos One. 2014.

  64. 64.

    Conant R. Semiaquatic snakes of the genus Thamnophis from the isolated drainage system of the Rio Nazas and adjacent areas in Mexico. Copeia. 1963(3):473–99.

  65. 65.

    Damon PE, Nieto-Obregón J. Un plegamiento Neogénico en Nayarit y Jalisco y evolución geomorfológica del Rio Grande de Santiago. Asociación Ingenieros de Minas, Metalurgistas y Geólogos de México, Memoria, XII Convención Nacional. 1979;156-191.

  66. 66.

    Luhr JF, Henry CD, Housh TB, Aranda-Gómez JJ, McIntosh WC. Early extension and associated mafic alkali volcanism from the southern Basin and Range Province: geology and petrology of the Rodeo and Nazas volvanic fields, Durango (Mexico). Geol Soc Am Bull. 2001;113:760–73.

    CAS  Google Scholar 

  67. 67.

    Strain WS. Blancan mammalian fauna and Pleistocene formations. Hudspeth County, Austin: Texas Memorial Museum; 1966. p. 10.

    Google Scholar 

  68. 68.

    Reeves CC. Pluvial Lake Palomas, northwestern Chihuahua, Mexico. New Mexico Geol Soc Guidebook. 1969;20:143–55.

    Google Scholar 

  69. 69.

    Wood RM, Mayden RL. Speciation and anagenesis in the genus Cyprinella of Mexico (Teleostei: Cyprinidae): a case study of Model III allopatric speciation. Rev Fish Biol Fisher. 2002;12:253–71.

    Google Scholar 

  70. 70.

    Domínguez-Acosta M, Gill TE. PIXE based geochemical characterization of the pluvial lake Palomas system, Chihuahua, Mexico, Procceedings of the 11th International Congerence on PIXE and its Analytical Applications, Puebla, Mexico; 2007. p. 2–33.

    Google Scholar 

  71. 71.

    Gustavson TC. Arid Basin Depositional Systems and Paleosols: Fort Hancock and Camp Rice Formations (Pliocene-Pleistocene), Hueco Bolson, West Texas and Adjacent Mexico: The University of Texas at Austin, Bureau of Economic Geology, Report of Investigations; 1991. No. 198. p. 49.

  72. 72.

    Galloway WE, Ganey-Curry PE, Li X, Buffler RT. Cenozoic depositional history of the Gulf of Mexico Basin. Am Assoc Petr Geol B. 2000;84:1743–74.

    Google Scholar 

  73. 73.

    Rosenthal J, Forstner MRJ. Effects of Plio-Pleistocene barrier on Chihuahua Desert herpetofauna. In: Hoyt CA, Karges J, editors. Proceedings of the sixth symposium on the natural resources of the Chihuahua Desert region; 2004. p. 269–82.

    Google Scholar 

  74. 74.

    Leeder MR, Mack GH. Basin-fill incision, Rio Grande and Gulf of Corinth rifts: convergent response to climatic and tectonic drivers. In: Nichols G, Williams E, Paola C, Wiler Sons J, editors. Sedimentary processes, environments and basins: A tribute to Peter Friend (Special Publication 38 of the IAS); 2009. p. 9–27.

    Google Scholar 

  75. 75.

    Mack GH, James WC.. Control of basin symmetry on fluvial lithofacies. Camp Rice and Palomas Formations (Plio-Pleistocene) in the southern Rio Grande Rift, USA. In: Marzo M, Pulgdefabregas C, editors. Alluvial sedimentation. Special Publication of the International Association of Sedimentologists; 1993;17. p. 439-449.

  76. 76.

    Berggren WA, Hilgen FJ, Langereis CG, Kent DV, Obradovich JD, Raffi I, Raymo ME, Shackleton NJ. Late Neogene chronology: New pespectives in high-resolution stratigraphy. Geol Soc Am Bull. 1995;107:1272–87.

    Google Scholar 

  77. 77.

    Mack GH, Love DW, Seager WR. Spillover model for axial rivers in regions of continental extension: the Rio Mimbres and Rio Grande in the southern Rio Grande rift. USA. Sedimentology. 1997;44:637–52.

    Google Scholar 

  78. 78.

    Ortega-Ramírez JR, Valiente-Banuet A, Urrutia-Fugauchi J, Mortera-Gutiérrez CA, Alvarado-Valdéz G. Paleoclimatic changes during the Late Pleistocene-Holocene in Laguna Babícora, near the Chihuahua Desert, Mexico. Can J Earth Sci. 1998;35:1168–79.

    Google Scholar 

  79. 79.

    Castiglia PJ, Fawcett PJ. Large Holocene lakes and climatic change in the Chihuahua Desert. Geology. 2006;34:113–6.

    Google Scholar 

  80. 80.

    Roy PD, Jonathan MP, Pérez-Cruz LL, Sánchez-Córdova MM, Quiroz-Jiménez JD, Romero FM. A millennial-scale Late Pleistocene-Holocene palaeoclimatic record from the western Chihuahua Desert, Mexico. Boreas. 2012;41:707–17.

    Google Scholar 

  81. 81.

    Schönhuth S, Shiozawa DK, Dowling TE, Mayden RL. Molecular systematics of western North American cyprinids (Cyprinodontiformes: Cyprinidae). Zootaxa. 2012;3586:281–303.

    Google Scholar 

  82. 82.

    Girard CF. Researches upon the cyprinoid fishes inhabiting the fresh waters of the United States, west of the Mississippi Valley, from specimens in the museum of the Smithsonian Institution. Proc Nat Sci Philadelphia. 1856;8:165–213.

    Google Scholar 

  83. 83.

    Mayden RL, Hillis DM. Natural history and systematics of the largemouth shiner, Cyprinella bocagrande (Teleostei: Cypriniformes), with comments on conservation status. Copeia. 1990:004–1011.

  84. 84.

    Broughton RE, Gold JR. Phylogenetic relationships in the North American cyprinid genus Cyprinella (Actinopterygii: Cyprinidae) based on mitochondrial ND2 and ND4L gene sequences. Copeia. 2000:1–10.

  85. 85.

    Burr BM. A review of the Mexican Stoneroller, Campostoma ornatum Girard (Pisces: Cyprinidae). Trans San Diego Soc Nat Hist. 1976;18:127–43.

    Google Scholar 

  86. 86.

    Raney EC, Webster DA. The spring migration of the common White Sucker, Catostomus c. commersonnii (Lacépède), in Skaneateles Lake Inlet, New York. Copeia. 1942:139–48.

  87. 87.

    Geen GH. Reproduction of three species of suckers (Catostomidae) in British Columbia. Thesis Dissertation. Canada: University of British Columbia; 1958.

    Google Scholar 

  88. 88.

    Hohler DB. A dwarfed population of Catostomus rimiculus (Catostomidae: Pisces) in Jenny Creek, Jackson County, Oregon. Thesis Dissertation. USA: Oregon State University; 1981.

    Google Scholar 

  89. 89.

    Tyus HM, Karp CA. Spawning and movements of Razorback Sucker, Xyrauchen texanus, in the Green River basin of Colorado and Utah. Southwest Nat. 1990;35:427–33.

    Google Scholar 

  90. 90.

    Rinne JN. Reproductive biology of the Rio Grande sucker, Catostomus plebeius (Cypriniformes), in a montane stream, New Mexico. Southwest Nat. 1995;40:237–41.

    Google Scholar 

  91. 91.

    Bunt CM, Cooke SJ. Post-spawn movements and habitat use by greater redhorse, Moxostoma valenciennesi. Ecol Freshw Fish. 2001;10:57–60.

    Google Scholar 

  92. 92.

    Schmetterling DA, McFee JA. Migrations by fluvial largescale suckers (Catostomus macrocheilus) after transport upstream of Milltown Dam, Montana. Northwest Sci. 2006;80:18–25.

    Google Scholar 

  93. 93.

    Cole DD, Mock KE, Cardall BL, Crowl TA. Morphological and genetic structuring in the Utah Lake sucker complex. Mol Ecol. 2008;17:5189–204.

    CAS  PubMed  Google Scholar 

  94. 94.

    Metcalfe SE, O’Hara SL, Caballero M, Davies SJ. Records of Late Pleistocene-Holocene climatic change in Mexico - a review. Quaternary Sci Rev. 2000;19:699–721.

    Google Scholar 

  95. 95.

    Metcalfe SE. Late Quaternary environments of the Northern deserts and central Transvolcanic Belt of Mexico. Ann Mo Bot Gard. 2006;93:258–73.

    Google Scholar 

  96. 96.

    Mastretta-Yanes A, Moreno-Letelier A, Piñero D, Jorgensen TH, Emerson BC. Biodiversity in the Mexican highlands and the interaction of geology, geography and climate within the Trans-Mexican Volcanic Belt. J Biogeogr. 2015;42:1586–600.

    Google Scholar 

  97. 97.

    Brown WM, George MJR, Wilson AC. Rapid evolution of animal mitochondrial DNA. P Natl Acad Sci Usa. 1979;76:1967–71.

    CAS  Google Scholar 

  98. 98.

    Avise JC. The history and purview of phylogeography: A personal reflection. Mol Ecol. 1998;7:371–9.

    Google Scholar 

  99. 99.

    Bradley RD, Baker RJ. A test of the genetic species concept: cytochrome-b sequences and mammals. J Mammal. 2001;82:960–73.

    Google Scholar 

  100. 100.

    Berendzen PB, Simons AM, Wood RM. Phylogeography of the northern hogsucker, Hypentelium nigricans (Teleostei: Cypriniformes): Genetic evidence for the existence of the ancient Teays River. J Biogeogr. 2003;30:1139–52.

    Google Scholar 

  101. 101.

    Bart HL, Reneau PC, Doosey MH, Bell CD. Evolutionary divergence of duplicate copies of the growth hormone gene in suckers (Actinopterygii: Catostomidae). Int J Mol Sci. 2010;11:1090–102.

    CAS  PubMed  PubMed Central  Google Scholar 

  102. 102.

    Baird SF, Girard C. Descriptions of nre species of fishes collected in Texas, New Mexico and Sonora, by Mr. John H. Clark, on the US and Mexican Boundary Survey, and in Texas by Capt. Stuart Vliet. Proc Acad Nat Sci Phila. 1985;7:24–9.

    Google Scholar 

  103. 103.

    Allendorf FW, Luikart G. Conservation and the genetics of populations. Oxford: Blackwell Publishing Ltd.; 2007.

    Google Scholar 

  104. 104.

    SEMARNAT. Norma Oficial Mexicana 059: Protección ambiental-especies nativas de México de flora y fauna silvestres-categorías de riesgo y especificaciones para su inclusión, exclusión o cambio-lista de especies en riesgo. México: Diario Oficial de la Federación; 2010.

    Google Scholar 

  105. 105.

    Clements MD, Bart HL Jr, Hurley DL. A different perspective on the phylogenetic relationships of the Moxostomatini (Cyprinoformes: Catostomidae) based on cytochrome-b and growth hormone intron sequences. Mol Phylogenet Evol. 2012;63:159–67.

    PubMed  Google Scholar 

Download references


We thank Rogelio Rosas Valdez and Hugo Mejia Madrid for field support in obtaining fish samples. Thanks to Javier Pérez Tris and Silvia Perea for valuable comments, and to the reviewer Tom Jameson from Chester Zoo and reviewers from The Lucidus Consultancy, Alan Pike and Kathleen Hills, for English corrections and suggestions.


This work was partial supported by the Coordinación de la Investigación Científica (CIC-UMSNH 2014), the Ministerio de Economía y Competitividad of Spain (grant reference CGL2013-41375-P) and FEDER Spain (grant reference CGL2016-75262-P). Both resources were used to finance fieldwork, laboratory (DNA isolation and amplification of molecular markers) and sequencing service, as well as for the publication of this work. The Consejo Nacional de Ciencia y Tecnología CONACYT with sabbatical grant and the scholarship to ODD (grant number: 260596) and DKCS (number grant: 254270/411621) respectively. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

Availability of data and materials

The datasets generated and/or analyzed during the current study are included in the research article (and its additional files). Also, all DNA alignments and phylogenetic trees are available via Treebase ( All DNA sequences obtained in this work are available on GenBank under the Accession Numbers MG203619-MG203705 for cytb gene, and MG203706-MG203869 for GHI region.

Author information




All authors contributed intellectually to the work and approved the submission. ODD, RRP, LTM, YHD, ID and DKCS participated in fieldwork. DKCS, LTM and JRPB carried out the work lab. DKCS, LTM, JRPB, ODD and ID analyzed the data and wrote the manuscript. All authors have read and approved the manuscript for submission.

Corresponding author

Correspondence to Omar Domínguez-Domínguez.

Ethics declarations

Ethics approval and consent to participate

This research has been performed within an appropiate ethical frameworks. Field and laboratory protocols used in this study, including sampling procedures, were reviewed and approved by the Mexican Ministry of Environmental and Natural Resources (SEMARNAT), under collection permit number PPF/DGOPA-262/17. Approval by an ethics committee was not necessary because this research did not include animal experimentation, and the samples were obtained using tricaine mezylate (MS-222) in fishes as anesthetic, according with the Official Mexican Norm NOM-051-ZOO-1995 and NOM-033-SAG/ZOO-2014 for humanitarian treatment in the mobilization of animals.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional information

The original version of this article was revised: one of the authors’ names was incorrectly spelled. Originally the author name has been published as José Ramón Pardons-Blas. The correct author name is José Ramón Pardos-Blas

Additional files

Additional file 1:

Evolutionary substitution model and estimated parameters for cytb and GHI by Akaike Information Criterion. (DOC 35 kb)

Additional file 2:

Phylogenetic tree (Bayesian Inference/Maximum Likelihood) based on mitochondrial gene cytb for the Pantosteus plebeius-nebuliferus species-group. (EPS 1046 kb)

Additional file 3:

Phylogenetic tree (Bayesian Inference/Maximum Likelihood) based on nuclear gene GHI for the Pantosteus plebeius-nebuliferus species-group. (EPS 1151 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Corona-Santiago, D.K., Domínguez-Domínguez, O., Tovar-Mora, L. et al. Historical biogeography reveals new independent evolutionary lineages in the Pantosteus plebeius-nebuliferus species-group (Actinopterygii: Catostomidae). BMC Evol Biol 18, 173 (2018).

Download citation


  • cytb
  • GHI
  • Evolutionary history
  • Pantosteus plebeius-nebuliferus species-group
  • Sierra madre occidental
  • Biogeographic patters
  • Independent evolutionary history