The genus Algansea is one of the most representative freshwater fish groups in central Mexico due to its wide geographic distribution and unusual level of endemicity. Despite the small number of species, this genus has had an unsettled taxonomic history due to high levels of intraspecific morphological variation. Moreover, several phylogenetic hypotheses among congeners have been proposed but have had the following shortcomings: the use of homoplasious morphological characters, the use of character codification and polarisation methods that lacked objectivity, and incomplete taxonomic sampling. In this study, a phylogenetic analysis among species of Algansea is presented. This analysis is based upon two molecular markers, the mitochondrial gene cytochrome b and the first intron of the ribosomal protein S7 gene.
Bayesian analysis based on a combined matrix (cytochrome b and first intron S7) showed that Algansea is a monophyletic group and that Agosia chrysogaster is the sister group. Divergence times dated the origin of the genus around 16.6 MYA, with subsequent cladogenetic events occurring between 6.4 and 2.8 MYA. When mapped onto the molecular phylogenetic hypothesis, the character states of three morphological characters did not support previous hypotheses on the evolution of morphological traits in the genus Algansea, whereas the character states of the remaining six characters partially corroborated those hypotheses.
Monophyly of the genus Algansea was corroborated in this study. Tree topology shows the genus consists of three main lineages: Central-Eastern, Western, and Southern clades. However, the relationships among these clades remained unresolved. Congruence found between the available geological and climatic history and the divergence times made it possible to infer the biogeographical history of Algansea, which suggested that vicariance events were responsible for the evolutionary history of the genus. Interestingly, this pattern was shared with other members of the freshwater fish fauna of central Mexico. In addition, molecular data also show that some morphological traits alleged to represent synapomorphies in previous studies were actually homoplasies. Others traits were corroborated as synapomorphies, particularly in those species of a subgroup corresponding with the Central-Eastern clade within Algansea; this corroboration is interpreted as a result of evolutionary adaptations.
Mexico lies between the Nearctic and Neotropical biogeographical zones and is considered to be a transitional zone. Because of this, it is possible to find fauna with different evolutionary origins . One of the world's great tropical-subtropical highlands is the massive uplift known as the Mesa Central of Mexico (MCM) and its southern limit, the Trans Mexican Volcanic Belt (TMVB). Since the Miocene, this zone has experienced an active geological history, which has promoted a complex surface configuration, including a wide variety of ecosystems. The freshwater fish fauna of the MCM [2, 3] is unique, with around 78 species, and is represented by an unusual level of endemicity (70%) [4, 5]. Most of the endemic species are represented by monophyletic groups that have undergone a diversification process within central Mexico [6, 7] such as the entire subfamily Goodeinae (41 species), the Atherinopsid genera Chirostoma (19 species), and three endemic genera belonging to the family Cyprinidae: Algansea (7 species), Evarra (3 species), and Yuriria (3 species) [6, 8–12].
Knowledge about the diversification processes of freshwater fish in central Mexico, including their origin as well as their evolutionary and biogeographical history, is still incomplete. Several hypotheses regarding the biogeography of freshwater fish in the region have been discussed in some detail in several studies [2, 3, 13–17]. These authors described general patterns using occurrence data and morphological comparisons. More recently, studies that incorporated molecular approaches in a phylogenetic context have been conducted to elucidate the biogeographical and evolutionary history of fish in central Mexico among groups such as poecilids , goodeids [7, 9, 19], and cyprinids [10, 12]. All of these studies agree that the historical biogeography of central Mexico and its freshwater fish fauna is linked to the intense geological activity since the early Miocene. This activity has generated a complex hydrological system characterised by high dynamism and the formation and destruction of drainages. This dynamism promoted vicariance, taxon-pulse, and species-pulse events . However, the complexity of these biogeographic patterns and the few fish groups studied thus far make it necessary to study other co-distributed fish groups, such as members of the genus Algansea, to formulate a more robust biogeographical scenario of the area.
The first attempts to uncover phylogenetic patterns of freshwater fish in central Mexico were based on just a few morphological traits, resulting in a non-robust hypothesis regarding the evolutionary history of the groups [2, 3, 21, 22]. More recently, phylogenetic studies on different freshwater fish families occurring in the region that were based on various molecular markers revealed results contradictory to those results obtained with morphological characters [9, 11, 12, 23]. Such molecular approaches demonstrated that fish diversity was in fact underestimated because several new species were described; in addition, some morphological characters commonly used in the classification of these groups were homoplasies. The genus Algansea is one of these groups that exemplify the particular problem of using morphological characters as the only source of information to establish a phylogenetic hypothesis. Although this genus possesses a few representative species, Algansea aphanea, Algansea avia, Algansea barbata, Algansea lacustris, Algansea monticola, Algansea popoche, Algansea tincella, , and Algansea amecae , its taxonomic history has been unstable, particularly due to high levels of intraspecific morphological variability found in the most widely distributed species within the genus, A. tincella .
The first significant attempt to elucidate the systematic relationships among species of Algansea was made by Barbour and Miller . These authors studied only six valid species and two subspecies and found two groups: (1) species with maxillary barbels (A. monticola - including A. m. monticola and A. m. avia -, A. barbata and A. aphanea) and (2) species without maxillary barbels (A. lacustris, A. popoche and A. tincella). In addition, based on the presence of a plate-like dermosphenotic bone, the genus Gila was determined to be the closest relative to Algansea . However, current molecular studies of Mexican cyprinids do not support this relationship, and these analyses established the genus Agosia as the sister group of Algansea [10, 12]. Subsequent to Barbour and Miller's hypothesis, the subspecies A. m. avia was validated as an independent species, A. avia. This species, together with A. monticola, A. barbata, and A. aphanea, are the barbeled species. The last taxonomic change implied that the highly differentiated population of A. tincella from the Ameca River should be recognised as A. amecae . All the aforementioned suggest that the morphological characters on which those first analyses were based have evolved independently several times throughout the evolutionary history of cyprinids , resulting in controversial classification schemes .
Molecular markers have proven to be very useful tools when elucidating the phylogenetic relationships of vertebrate groups. In the last few years, the mitochondrial gene cytochrome b (cyt b) and the first intron of the S7 ribosomal protein gene have been used in different phylogenetic analyses to provide insights into the fish evolutionary history at various taxonomic levels. Both genes have proven their utility not only for inferring the phylogenetic relationships among closely related species but also for investigating intraspecific variation and even for establishing species boundaries [27–31]. Accordingly, in this paper we used both cytb and the first intron of S7 to infer the phylogenetic relationships among species and populations of Algansea throughout their entire distribution range in central Mexico. Additionally, the results of these analyses were used to: 1) compare with the previous phylogenetic hypotheses based on morphological data; 2) characterise the evolution of the main morphological characters used for identifying natural groups; and 3) establish possible biogeographic scenarios in which species and populations may have evolved in the context of the geological and climatic history and contrast this history with other co-distributed taxa.
Sequence data patterns
In the cytb sequences (1,140 bp), a total of 425 sites (37%) were variable, and 279 (24%) were parsimony informative. As expected for a protein-coding gene, third codon positions were the most variable (320), followed by the first (74) and the second (31) positions. For the S7 intron, the size was 965 bp including gaps, and there were 279 variable sites, 103 of which were parsimony informative. Comparison of the individual genes' phylogenetic performance under taxon sampling of the combined analysis is found in Table 1. Estimated parameters using Modeltest  are shown in Table 1. As in other North American cyprinids' cytb, a low Guanine frequency in the third codon positions was found [12, 23, 33]. S7 sequences were AT-rich, as has been previously reported for other families of freshwater fishes [34, 35] and other North American cyprinids . Despite such apparent bias, the χ2 test for base homogeneity indicated that the base frequency distribution was always homogeneous among taxa.
Parameters and statistics summary of phylogenetic analyses
cyt b + S7
Number of sites
Number of variables sites
Most parsimonious trees
A = 24%
A = 20%
A = 35%
A = 26%
A = 31%
C = 25%
C = 26%
C = 34%
C = 29%
C = 15%
G = 26%
G = 13%
G = 11%
G = 17%
G = 20%
T = 25%
T = 40%
T = 20%
T = 28%
T = 34%
χ2 test of base frequencies
χ2 = 4.62
χ2 = 1.28
χ2 = 1.28
χ2 = 21.68
χ2 = 3.71
gl = 153
gl = 153
gl = 153
gl = 153
gl = 54
p = 1.00
p = 1.00
p = 1.00
p = 1.00
p = 1.00
BIC = Bayesian information
Details of the Maximum Parsimony (MP) and Bayesian Inference (BI) analyses are summarised in Table 1. All the analyses (mitochondrial, nuclear, and combined data) showed almost identical and well-supported topologies; however, the MP tree was less resolved than the Bayesian consensus tree. Thus, we focus our discussion on the more resolved Bayesian tree and only summarise the results of the MP analysis.
All the cladograms retrieved Algansea as monophyletic with high support values (Figures. 1, 2 and 3). Also, all the species formed monophyletic assemblages, except A. tincella and A. amecae in the nuclear gene cladogram (Figure 2). In addition, all the trees showed Agosia chrysogaster as the sister group to Algansea (Figures 1, 2 and 3).
Sister group relationships among species of Algansea were generally consistent with the cytb or combined data sets. Three main clades corresponding to the Central-Eastern, Western, and Southern regions of the whole distribution range of the genus Algansea were recovered with high support values. The Central-Eastern clade is represented by the species A. barbata, A. amecae, A. lacustris, and A. tincella (Figures 1, 2 and 3). The Western clade included the two species A. monticola and A. avia from the northwestern headwaters of the Santiago River and the lower Santiago, respectively (Figures 1 and 3). The last clade is represented only by A. aphanea (Figures 1, 2 and 3), a species that inhabits the Southern Pacific river basins of Ayuquila-Armería and Coahuayana-Tamazula (see Additional File 1 and Figure 4). The basal relationships among the three main clades were not resolved.
When intron S7 was used, there were no differences between both phylogenetic methods (MP and BI); inclusion or exclusion of indels in the MP analyses gave the same topology. However, with respect to the cytb and the combined data sets, a lower resolution tree was formed with a resulting basal polytomy (Figure 2).
The main difference between the three data sets was the position of the species A. barbata with respect to all the congeners. In the cytb and combined data analyses, A. barbata is included in the Central-Eastern clade (Figures 1 and 3). However, in the S7 analysis, A. barbata is associated with a member of the Western clade (Figure 2). The most resolved and best-supported hypothesis was obtained with combined Bayesian analyses (Figure 3).
For cytb, the genetic divergences found among species ranged from = 0.020 ± 0.004 found between A. tincella and A. lacustris to 0.090 ± 0.008 found between A. aphanea and A. avia. Genetic divergences among species of the Central-Eastern clade showed relatively low values, with the highest value (Dp = 0.038 ± 0.006) observed between A. lacustris and A. amecae and the lowest value (Dp = 0.020 ± 0.004) observed between A. lacustris and A. tincella. Divergences between the species of the Western clade accounted for the higher values, with = 0.078 ± 0.008. Among the three main clades, the observed divergence values were similar: = 0.083 ± 0.008 between the Southern and Central-Eastern clades, followed by = 0.086 ± 0.007 between the Western and Central-Eastern clades, and = 0.088 ± 0.007 between the Southern and Western clades. Within species, the lowest divergence value was = 0.002 ± 0.001 between both analysed populations of A. aphanea. The highest divergences within species were found in populations of A. tincella, particularly between the population from the Quitupan River (Balsas Basin), which had a divergence value of = 0.014-0.018 ± 0.004 with respect to the remaining A. tincella populations. For the S7 intron, the lowest divergences between species were found between A. tincella and A. lacustris, with values ranging from = 0.005 ± 0.002 to = 0.008 ± 0.004. The highest divergences were found between A. aphanea and A. monticola when compared to the remaining species of Algansea, which had values ranging from Dp = 0.020 ± 0.005 to = 0.036 ± 0.006.
Molecular clock estimates among the genera Algansea and Agosia placed the time to the most recent common ancestor (TMRCA) at 14.7 MYA (16.3-10.6 MYA). The estimated age of the diversification of the Central-Eastern, Western, and Southern clades is 6.2 MYA (7.8-4.7 MYA). The split between A. avia and A. monticola occurred 4.5 MYA (6.2-3.2 MYA). Within the Central-Eastern clade, the oldest splitting event occurred about 2.9 MYA (4.1-2.1 MYA) between A. barbata and the clade containing the ancestor of the A. amecae, A. lacustris, A. tincella. Following the cladogenetic events, A. amecae and the clade A. lacustris, A. tincella were dated to 2.4 MYA (3.1-1.9 MYA), and A. lacustris and A. tincella were dated to 1.9 MYA (2.6-1.3 MYA).
Mapping some morphological traits commonly used in the classification of Algansea onto the molecular phylogenies clearly shows that the orientation of the mouth, either terminal or upturned, represents a homoplasious character (Figure 5A). Except for a subgroup of the Central-Eastern clade formed by A. amecae, A. lacustris, and A. tincella, the supraethmoid orientation and the gut flexure were homoplasies as well (Figure 5D, F). The lack of maxillary barbels, a higher number of gill rakers (Figure 5B), an irregular supraethmoid margin, and well-developed epiotic bones (Figure 5C) were synapomorphies for the aforementioned subgroup. In addition, the remaining two characters mapped, i.e., standard length (as was coded in the present study) and the neurocranium dome, were synapomorphies for the whole genus, supporting the monophyly of the Central-Eastern, Western and Southern clades (Figure 5E, G).
Barbour and Miller's  and Jensen and Barbour's  analyses demonstrated that the standard length was not useful for establishing sister-group relationships among species of Algansea. In our study, this character was re-evaluated and was found to be informative with the proper codification, i.e., smaller (< 100 mm) or larger size (> 100 mm). Standard length character was found to be a synapomorphy for a species group formed by Southern and Western clades (smaller group) and the Central-Eastern clade (larger group) (Figure 5E).
The monophyly of the genus Algansea was corroborated using molecular markers. The position of the monotypic genus Agosia as a sister group of Algansea was consistent in all analyses. The relationship of the genus Gila with Algansea as previously suggested by Barbour and Miller , was rejected. Our results support previous studies in which a close phylogenetic relationship was found between Algansea and Agosia [10, 12, 37].
Sequence variation for the mitochondrial cytb and nuclear S7, independently or combined, did not provide enough information to resolve the relationships among the three well-supported groups: the Central-Eastern, Western, and Southern clades. A lower supported relationship, Western clade + Southern clade, was found when the cyt b and the combined matrix were used (Figures 1 and 3), whereas a basal polytomy was found using the S7 intron (Figure 2). Phylogenetic reconstruction shows no resolution in basal members of Algansea, which might be the result of distinct evolution rates of the markers used in this study. The evolution rate of cytb was almost twice as high as the evolution rate of nuclear S7 intron. This fact had been previously reported for other fish [30, 31, 38]. However, Schönhuth et al.  recently showed that, although the evolution rates of cytb and S7 differ markedly, both markers had a similar performance in reconstructing the phylogenetic relationships among species and genera of North American cyprinids.
Therefore, this lack of resolution, especially at the base of the cladogram of the genus, as shown by separate and combined analyses, may be explained as a result of inadequate sampling of data, sampling of taxa, or a compositional base bias. However, such phylogenetic pattern could be caused by a process implicit to the species' evolutionary history as well, such as a rapid or simultaneous speciation .
The combined data analysis provided the best resolution for establishing the phylogenetic relationships among members of Algansea, and that analysis also shows the highest support values. Particularly for the Central-Eastern clade, the combined analysis represented the best estimate for the relationships among the included species (Figure 3).
Otherwise, the inconsistent phylogenetic position of A. barbata, in particular in the mitochondrial and nuclear trees, is the main difference found between topologies obtained by using mitochondrial and nuclear genes. The lack of congruence among molecular markers could be explained by gene duplication, hybridisation (ancient or recent), or incomplete lineage sorting [40, 41].
Despite the existence of pseudogenes in the S7 ribosomal protein gene in mice , the amplification of single PCR products for the first and second introns among distant fish species suggests that the process of duplication of the S7 gene is considerably lower in fish than in mammals . In addition, no duplication cases have been found in previous studies on fish where the S7 intron was used [12, 28, 30, 31, 35, 38, 42–44]; therefore, a duplication event probably did not occur within Algansea.
Distinguishing between hybridisation and deep coalescence is complicated because both processes generate similar phylogenetic patterns [40, 41, 45, 46]. Although introgressive hybridisation have played an important role in the evolutionary history of some North American cyprinids [47–50], Barbour and Miller  recognised hybrids in a zone where A.tincella and A. popoche are sympatric. There are two reasons for discarding a recent and/or historical hybridisation as the cause of incongruence between mitochondrial and nuclear genes-based hypotheses (due to the unstable position of A. barbata). On the one hand, disjunctive historical and current distribution patterns of A. barbata and A. monticola (see Biogeographic Implications and Figure 4) suggest that no contact between the two species occurred. Second, as in other North American cyprinids [48–50] and other fish groups [51, 52], when a hybridisation event is detected by phylogenetic analysis, it needs further corroboration using data from different sources. However, data gathered thus far on the evolutionary history of Algansea show contradicting results with respect to the relationships among A. barbata and other congeners. For instance, morphological characters support the relationships between A. barbata and A. aphanea , while topology obtained through the cyt b gene placed A. barbata within the Central-Eastern clade as the sister species of the not-barbeled species (Figure 1). The topology obtained through the first S7 intron shows that A. barbata and A. monticola are sister taxa (Figure 2).
Based on the premise that the S7 intron shows a slower evolution rate, it is possible that inconsistencies may reflect insufficient time to complete the lineage sorting. An incomplete lineage sorting is found in relatively recent diversification events, and in particular when nuclear genes are employed a longer time is required to reach reciprocal monophyly . Given the rapid speciation of the genus Algansea in central Mexico, relationship between A. barbata and A. monticola obtained in the S7 intron topology is attributed to the retention of an ancestral polymorphism. In fact, incongruence was observed between the three kinds of evidence, morphological, and the mitochondrial and nuclear genes, reflecting the random nature of lineage sorting [40, 46].
Based on the sister group relationships between Algansea and Gila as found by Barbour and Miller , it was postulated that the genus Algansea derived from a widespread ancestor in the highlands of western Mexico during the late Tertiary. Our findings indicate that the sister group of Algansea is the genus Agosia. This genus is distributed from the Western Sierra Madre Occidental in Mexico to the Rocky Mountains in the United States (Figure 4), supporting a western origin and a colonisation route to central Mexico by the ancestor of the genus Algansea.
According to the estimated divergence times, the cladogenetic event involving the ancestor of Agosia and Algansea was dated about 14.7 MYA (16.3-10.6 MYA), during the middle Miocene (Figure 6). This event is associated with the uplift of the western part of the MCM and the southern part of the Sierra Madre Occidental. This uplift was promoted by tecto-volcanic activity in the region and was related to the opening of the Proto Gulf of California during the lower and middle Miocene [54, 55]. This biogeographical pattern is in agreement with findings made in other freshwater fishes, such as goodeids (the separation of the Characodontini tribe from the remaining tribes of the subfamily Goodeinae), and with the separation of the two divergent groups of the poecilid Poeciliopsis spp., which occurred about 15.5 MYA and between 8 and 16 MYA, respectively [18, 20].
The separation of the three main clades of Algansea is dated at some point between the Miocene and Pliocene, about 6.2 MYA (7.8-4.7 MYA), a period recognised for a series of profound global palaeoclimatic changes and high geological activity in central Mexico . In particular, during the formation of the western part of the TMVB, several events were responsible for the isolation of major groups within Algansea and within other groups of freshwater fish (Figure 6). Including the genus Poeciliopsis , the formation of the two main lineages within the cyprinid Notropis of central Mexico , and the diversification of three tribes among the Goodeinae , all of which coincided with these events.
Within the Western clade, the cladogenesis of A. avia and A. monticola, both inhabitants of the Santiago River, was dated to the Pliocene, 4.5 MYA (6.2-3.2 MYA). Currently, these species are separated by a deep canyon in the main channel of the Santiago River; therefore, the formation of this geological structure may be responsible for the separation of these species . The origin of the Santiago River Canyon is related to the Santa Rosa faulting, which is dated to the early Pliocene  (Figure 6). The presence of A. avia in the Compostela River in Nayarit can be attributed to a river capture event frequently occurring in the upper part of rivers in central México due to headwaters erosion. This event is dated to less than 1 MYA  and coincides with the distribution pattern of the goodeid Xenotoca eiseni.
The event that originated the separation between A. barbata and the ancestor of the clade formed by A. tincella, A. lacustris, and A. amecae is calculated to have occurred 2.9 MYA (4.1-2.1 MYA), reinforcing the idea that climatic and geological changes in the MCM during that time led to the diversification of Algansea. During this period, intense tecto-volcanic activity combined with a dry period produced a strong compartmentalisation of palaeolakes in central Mexico. This activity could have been responsible for population isolation and subsequent speciation events [56, 59]. Likewise, this distributional pattern is in agreement with the separation of two lineages of the goodeid tribe Girardinichthyini that occurred in the palaeolakes of the region and are dated between the upper Miocene and lower Pleistocene .
The second event within the Central-Eastern clade was the isolation of A. amecae in the Ameca River from the ancestor of the remaining members of the Central-Eastern clade, A. tincella and A. lacustris. This cladogenetic event also took place in the upper Pliocene, approximately 2.4 MYA (3.1-1.9 MYA). This pattern clearly indicates an ancient connection between the upper Ameca River and the Lerma-Chapala River system, which was further disrupted by geological activity along the Ameca and San Marcos faults, which occurred throughout the Pliocene (Figure 6) . It is also related to the dry period reported in the area between 3.5 to 1.8 MYA . Ferrari et al.  mentioned the presence of volcano-lacustrine deposits along these faults, which probably indicates that the formation of these depressions occurred during the early Pliocene. These results are in agreement with findings in other groups of freshwater fishes that currently occur in both basins, such as the cyprinids Yuriria amatlana, Zoogoneticus tequila, and Ameca splendens (all endemic to the Ameca River basin) and their sister groups, inhabitants of the Lerma-Chapala system [7, 11].
The cladogenetic event of the species pair A. lacustris, endemic to Patzcuaro Lake, and A. tincella, a widespread species in the Lerma, Panuco, Balsas, and Verde de Santiago rivers basins and Cuitzeo Lake, is dated to the upper Pliocene, 1.9 MYA (2.6-1.3 MYA). Although the formation and evolution of Patzcuaro Lake remains unclear and controversy surrounds the age of the lake's origin, volcanic activity in the Tarasco Corridor has been previously proposed as one of the causes of diversification of other fish species in the region in the last 3 MYA .
The historical biogeography of the genus Algansea seems to be closely linked to the intense geological and climatic activity of central Mexico, which is one of the most complex and active geological areas in the world [61–63]. Based on both historical and geological data, most of the cladogenetic events within the genus Algansea were the result of vicariant events (Figure 6). This intense activity has been proposed as the main cause of diversification in other freshwater species of central Mexico, such as the Goodeids [7, 20, 64], Poeciliids [18, 65], helminth parasites [66–70], snakes , frogs [72, 73], and ambystomatid salamanders .
Most of the aforementioned vicariant events are highly consistent with the topology obtained from the combined analyses. Still, more groups must be studied in the same context for a better understanding of the processes that led to the diversification of the complex biota that occurs in this transitional zone.
Morphological evolution and taxonomic implications
North American cyprinids show a diverse and complex morphological evolution; this complexity is probably the main reason to consider theses fishes as one of the most taxonomically difficult groups on the continent .
According to our results, the barbeled species group (A. monticola, A. aphanea and A. barbata) that was based solely on morphology from previous phylogenetic hypotheses of the genus Algansea [21, 36] was not recovered as monophyletic; this result corroborates the conflictive role that the maxillary barbels have for the classification of the family Cyprinidae [25, 26]. However, even though maxillary barbels seem to have evolved independently several times , the presence of such structures in A. barbata is regarded as a plesiomorphic character for the clade to which this species belong. Instead, the lack of barbels in the subgroup of the Central-Eastern clade, together with a higher number of gill rakers, six other character states (except for body size smaller than 100 mm and the rounded neurocranium dome), are explained as derived characters (Figure 5A-G).
In particular, possessing a larger body is an adaptation in North American minnows that inhabit open-water habitats . This suggests that the morphological traits found in the members of Algansea that constitute the Central-Eastern clade are clearly associated with lacustrine environments. This adaptation might have been promoted when species existed in the paleolakes in central Mexico (see Biogeographic implications).
Barbour and Miller  described some characteristics of the ecological conditions of the habitats where species of Algansea occur. These data, and our own observations in the sampling localities, allow us to conclude that the species forming the Western and Southern clades are restricted to the upper parts of rivers, which are characterised by swift currents with rocky and gravel bottoms and clear and shallow water (less than one meter deep). On the other hand, within the member of the Central-Eastern clade, A. lacustris, is restricted to lakes; A. amecae and A. barbata live in streams with slow to moderate water flow, sandy silt bottoms, and water depths reaching more than one meter. Furthermore, A. tincella occurs in a variety of habitats ranging from small streams to lakes. These ecological aspects are closely tied with larger body size and, therefore, are found mainly in the members of the Central-Eastern clade.
Although our study did not include the species A. popoche (endemic to Chapala Lake), we could assume, based on the distribution of this species, that it would likely be included within the Central-Eastern clade. Barbour and Miller  proposed A. popoche as the sister species of A. lacustris. However, we were unable to collect specimens of A. popoche because this species may have gone extinct, and we could not further corroborate this idea. However, most of the recognised species within the genus Algansea (A. aphanea, A. avia, A. barbata, A. lacustris, A. monticola, A. tincella, , and A. amecae ) were corroborated by our phylogenetic analyses. Two species regarded as subspecies of A. monticola by Barbour and Miller , A. monticola avia and A. monticola monticola, exhibited genetic divergence levels in our study that allowed us to further test the hypothesis that A. avia is an independent species. These tests granted species rank to A. avia and corroborated the morphological differentiation found by Jensen and Barbour  and Barbour and Miller .
The present study is based on two independent sources of evidence (mitochondrial and nuclear molecular markers) that corroborate the monophyly of the genus Algansea. Despite the low resolution found in the phylogenies within this genus, three well-supported clades were recovered: the Southern clade, including A. aphanea; the Western clade, including A. avia and A. monticola; and the Central-Eastern clade, including A. barbata, A. amecae, A. lacustris, and A. tincella. In particular, the monophyly of the Central-Eastern clade was found, and our results did not support the groups proposed by Barbour and Miller  based on the presence of barbels. Historical biogeography patterns of this clade are shared with other groups of co-distributed fishes in central Mexico. Additionally, a larger body size in members of the Central-Eastern clade suggests an evolutionary adaptation that arose when these species colonised the lacustrine environments in which they occur in central Mexico.
We corroborate the western origin of the genus Algansea; however, the sister species is the monotypic Agosia. Most of the cladogenetic events are associated with vicariance occurring in the region as a result of tectonic activity and climatic changes. Such events seem to have also shaped the evolutionary history of other freshwater fish groups.
Species and Sampling
Fin clips were obtained from 48 specimens in 20 localities (Additional File 1 and Figure 4) belonging to seven putative species of the genus Algansea (A. aphanea, A. avia, A. barbata, A. lacustris, A. monticola, A. tincella, and A. amecae). Fishes were collected by electrofishing and seine nets. The only species not included in the analysis was A. popoche, an endemic species from Chapala Lake . This species has not been found in the locality despite a great sampling effort made by different research groups, raising the possibility that the species is now extirpated or extinct . Additionally, fin clips were collected from representative species of the cyprinid genera Agosia, Campostoma, Dionda, and Gila (Table 1). To test the monophyly of the genus Algansea, and due to the relationships found between Algansea and other North American cyprinids [10, 12], the first three genera were included as ingroups, leaving Gila as an outgroup. Voucher specimens of all species were deposited in the Museo Nacional de Ciencias Naturales (MNCN), Madrid, Spain, and at the Colección de Peces de la Universidad Michoacana de San Nicolas de Hidalgo (CPUM), Michoacán, México.
DNA extraction, PCR and sequencing
Total genomic DNA was extracted from ethanol-preserved fin clips according to standard CTAB and phenol-chloroform extraction procedures . Both the complete sequence of cytb and the first intron of ribosomal protein gene S7 were amplified via polymerase chain reaction (PCR). For both genes, reaction amplifications were carried out in 25 μl reactions containing: 2.5 μl 10× buffer with MgCl2 (Biotools), 0.5 μl of each dNTP (10 μM), 0.3 μl of each primer, 1-2 μl genomic DNA (50 ng/Ml), 1 unit of Taq DNA polymerase (Biotools), and distilled water to bring the final reaction volume to 25 μl.
For the cytb gene, the primers used were those primers in Machordom and Doadrio . For the first intron of S7, we used the primers of Chow and Hazama . The amplification was performed using the following conditions. For cytb, 35 cycles were used: denaturation at 94°C for 45 sec, annealing at 46°C for 1 min, and extension at 72°C for 1:30 min. A final extension at 72°C for 5 min was performed to completely extend the amplified product. For the S7 first intron, 30 cycles were performed using a denaturation at 94°C for 1 min, annealing at 54°C for 1:30 min, and an extension at 72°C for 2 min; then, a final extension was performed at 72°C for 7 min.
PCR products were purified with the QIAquick (QIAGEN) kit, checked on 1% agarose gels, and sequenced using the Big Dye Deoxy Terminator cycle-sequencing kit (Applied Biosystems Inc.) in an ABI PRISM 3700 DNA analyser. For cytb, chromatograms and alignments were visually checked and the S7 sequences were aligned using Clustal × 1.83 with default parameters . All sequences were revised and verified using MEGA 3.1 . GeneBank accession numbers are presented in Additional File 1. We obtained the complete mitochondrial cytb (total of 1140 bp) of 47 specimens from 25 localities (Additional File 1). Based on at least one representative specimen from each mitochondrial clade, we sequenced the first intron of the ribosomal protein S7 (965 bp, including gaps) for a subset of 19 specimens from 16 localities (Additional File 1). The cyt b matrix was increased with previously published sequences for Algansea species (Schönhuth et al. : DQ324088-DQ324092).
Homologous regions were aligned manually against previously published cytb sequences of Algansea . The nucleotide composition and base frequencies were examined, and the homogeneity test of base frequencies (for each single codon position and all the positions) was carried out for all taxa and both genes with the program PAUP 4.1 . Furthermore, the saturation of transition and transversion changes was checked for each gene by plotting the absolute number of changes of each codon position against their patristic distances (p). There was no evidence of saturation for any of the sequence data sets (data not shown).
Phylogenetic analyses were carried out for both genes separately and as a combined data set using Maximum Parsimony (MP) and Bayesian Inference (BI) analyses. The MP analyses were performed in PAUP* version 4.1 using heuristic searches (TBR Branch swapping; MULPARS option in effect) with 10 random stepwise additions of taxa. Analyses of the S7 gene and combined matrix were done including or excluding the indels.
For the BI analyses, the best-fit model for the different genes and each codon position (for cyt b in particular) were selected using Modeltest 3.7  based on the Bayesian Information Criterion (BIC). Bayesian analyses were performed in MrBayes 3.1.2 , using two independent runs of four Metropolis-Coupled Markov Chain Monte Carlo (MCMC) of 1,000,000 generations each to estimate the posterior probability distribution. The combined sequence matrices were partitioned per gene fragment, and independent model parameters were estimated for each partition. Topologies were sampled every 100 generations. Once the average standard deviation of split frequencies was less than 0.01, as suggested by MrBayes 3.1.2 , convergence between runs was checked. This was accomplished by comparing the 50% majority rule consensus tree for each run (after discarding the first 20,000, 25,000 and 10,000 generations for the cytb, S7, and combined gene data sets, respectively), and no incongruence between the runs were found.
The robustness of the clades was assessed using bootstrapping (1,000 pseudoreplicates) for the MP analyses and Bayesian posterior probabilities for the BI analyses.
The ancestral character-state reconstruction of the evolution of some of the morphological characters consisted of a subset of nine characters from a data matrix of twenty-six characters used in Barbour and Miller's  and Jensen and Barbour's  phylogenies (Additional File 2). These characters were selected because they were synapomorphies supporting particular clades in those studies. The software MacClade version 3.06  was used to reconstruct character evolution.
Molecular clock and divergences
The uncorrected p distances and absolute changes were calculated using all the specimens analysed. The average p distances between the main clades were also calculated for both markers using the program MEGA v.3.1 .
A relaxed molecular clock Bayesian approach using cytb gene (1140 bp) was performed to infer the time to the most recent common ancestor (TMRCA) in the different species and groups within Algansea and with respect to Agosia. Because of the lack of reliable fossil records for the genus Algansea, a molecular clock of 1.05% per million year, which has been established for cytb in North American Phoxinini by Dowling et al.  and widely applied to Euro Asiatic cyprinids [85–89], was used. The molecular clock was used to i) estimate the divergence times of the main cladogenetic events implicated in the origin and diversification of Algansea and ii) test if these estimates are in agreement with the historical geologic events in the central Mexico region that may be responsible for the origin and diversification of other co-distributed groups of fishes. Divergence times and their confidence intervals were estimated using a relaxed clock model in BEAST v1.4.6 , with a strategy that included branch rates drawn from an uncorrelated log-normal distribution . The estimates of divergence times was done with a Yule tree prior, the branch length substitution rate sampled from a prior normal distribution (with a mean value of 0.010 and a standard deviation of 0.001 ) and included just one representative of each species. The TrN+G substitution model was used, and three MCMC for 90 × 106 generations were run. The remaining parameters were set as default options and changed as recommended by the BEAST output file. We checked for the effective sample size (ESS) convergence and the stationary of the different analyses in Tracer 1.4  and combined the results in the BEAST module LogCombiner 1.4.4. After removing 10% of the generations from each analysis as "burn-in", the ESS exceeded 210 for all parameters.
The authors wish to thank the Tiacaque fish farm workers for providing the A. barbata specimen. We also thank the fisherman from the Uranden natural reserve and Patzcuaro Lake for providing specimens of A. lacustris. We thank Rogelio Rosas-Valdez and Jaquelina Bravo-Arteaga for their help during field work. We thank H. Gante for your revisions and comments to this work. This study was partially funded by grants from CGL2006-12325/BOS, the program PAPIIT-UNAM-IN209608, and CONACYT (grant number 83043) to GPPDL. RPR thank the Consejo Nacional de Ciencia y Tecnología for the scholarship.
Posgrado en Ciencias Biológicas, Instituto de Biología, Universidad Nacional Autónoma de México
Laboratorio de Biología Acuática, Facultad de Biología, Universidad Michoacana de San Nicolás de Hidalgo
Departamento de Zoología. Ap. Postal 70-153, Instituto de Biología, Universidad Nacional Autónoma de México
Departamento de Biodiversidad y Biología Evolutiva, Museo Nacional de Ciencias Naturales, CSIC
Morrone JJ: Hacia una síntesis biogeográfica de México.Revista Mexicana de Biodiversidad 2005, 76: 207–252.
Barbour CD: The systematics and evolution of the genus Chirostoma Swainson (Pisces, Atherinidae).Tulane Studies in Zoology and Botanical 1973, 18: 97–141.
Barbour CD: A biogeographical history of Chirostoma (Pisces: Atherinidae): A species flock from the Mexican Plateau.Copeia 1973, 1973: 533–556.View Article
Guzmán–Arroyo F: Osteología y variación no geográfica de la suspensión de la aleta anal de Goodea luitpoldi, (Osteichthyes:Goodeidae).Universidad, Ciencia y Tecnología 1994, 3: 33–41.
Miller RR, Minkley WL, Norris SM: Freshwater fishes of Mexico.1 Edition Chicago: University of Chicago Press 2005.
Miller RR, Smith ML: Origin and geography of the fishes of central Mexico.The Zoogeography of North American Freshwater Fishes(Edited by: Hocutt CH, Wiley EO). New York: Wiley–Interscience Publications 1986, 487–519.
Domínguez–Domíguez O, Doadrio I, Pérez–PoncedeLeon G: Historical biogeography of some river basins in central Mexico evidenced by their goodeine freshwater fishes: a preliminary hypothesis using secondary Brooks parsimony analysis.Journal of Biogeography 2006, 33: 1437–1447.View Article
Espinosa H, Gaspar M, Fuentes P: Listados faunísticos de México. III. Los peces dulceacuícolas mexicanos México DF: Instituto de biología, Universidad Nacional Autónoma de México 1993.
Doadrio I, Domínguez O: Phylogenetic relationships within the fish family Goodeidae based on cytochrome b sequence data.Molecular Phylogenetic and Evolution 2004, 31: 416–430.View Article
Schönhuth S, Doadrio I, Mayden R: A biogeographic perspective on the phylogeny of mexican cyprinids (Actinopterygii: Cyprinidae).Studies of north American desert fishes: in honor of E P (Phil) Pister, conservationist1 Edition(Edited by: Lozano–Vilano ML, Contreras–Balderas S). Monterrey, NL: Universidad Autónoma de Nuevo León 2006, 102–124.
Domínguez–Domínguez A, Pompa–Domínguez A, Doadrio I: A new species of the Genus Yuriria Jordan & Evermann, 1896 (actinopterygii, cyprinidae) from the Ameca Basin if the Central Mexican Plateau.Graellsia 2007, 63: 259–271.View Article
Schönhuth S, Doadrio I, Domínguez–Domínguez O, Hillis DM, Mayden R: Molecular evolution of southern North American Cyprinidae (Actinopterygii), with the description of the new genus Tampichthys from central Mexico.Molecular Phylogenetics and Evolution 2008, 47: 729–756.View ArticlePubMed
De Buen F: Los Lagos Michoacanos I. Caracteres generales. El Lago de Zirahuén.Revista de la Sociedad Mexicana de Historia Natural 1943, 4: 211–232.
Álvarez del Villar J: Ictiología Michoacana V. Origen y distribución de la ictiofauna dulceacuícola Michoacana.Anales de la Escuela Nacional de Ciencias Biológicas 1972, 19: 155–161.
Parenti L: A phylogenetic and biogeographic analysis of cyprinidontiform fishes (Teleostei, Atherinomorpha).Bulletin of the American Museum of Natural History 1981, 168: 335–557.
Echelle AA, Echelle AF: Evolutionary genetics of a "species flock": atherinid fishes on the Mesa Central of Mexico.Evolution of fish species flocks(Edited by: Echelle AA, Kornfield I). Orono: University of Maine Press 1984, 93–110.
Moncayo–Estrada R, Israde–Alcantara I, Garduño–Monroy VH: La cherehuita Hubbsina turneri De Buen (1941) (Pisces, Goodeidae): origen, distribución y su uso en la regionalización de la cuenca del lerma.Hidrobiológica 2001, 11: 1–13.
Mateos M, Sanjur OI, Vrijenhoek C: Historical biogeography of the livebearing fish genus Poeciliopsis (Poecilidae: Cyprinodontiformes).Evolution 2002, 56: 972–984.PubMed
Domínguez–Domínguez O, Alda F, Pérez–PoncedeLeón G, García–Garitagoitia JL, Doadrio I: Evolutionary history of the endangered fish Zoogoneticus quitzeoensis (Bean, 1898) (Cyprinodontiformes: Goodeidae) using a sequential approach to phylogeography based on mitochondrial and nuclear DNA data.BMC Evololutionary Biology 2008, 8: 161.View Article
Domínguez–Domínguez O, Pedraza–Lara CP, Gurrola–Sánchez N, Pérez–Rodríguez R, Alcaraz L, Perea S, Ornelas CP, Israde–Alcántara I, Garduño–Monroy VH, Doadrio I, Pérez–PoncedeLeón G, Brooks DR: Historical biogeography of the Goodeinae (Cyprinodontiforms).Viviparous fishes II(Edited by: Uribe–Aranzabal MC, Grier H). Florida: New Life Publications, in press.
Barbour CD, Miller RR: A revision of the genus Algansea.Miscellaneous Publications Museum of Zoology, University of Michigan 1978, 155: 1–72.
Lyons J, Polaco OJ, Cochran PA: Morphological variations among the Mexican lampreys (Petromyzontidae: Lampetra: subgenus Tetrapleurodon).Southwestern Naturalists 1996, 41: 365–374.
Schönhuth MS, Doadrio I: Phylogenetic relationships of Mexican minnows of the genus Notropis (Actinopterygii: Cyprinidae).Biological Journal of the Linnean Society 2003, 80: 323–337.View Article
Pérez–Rodríguez R, Pérez–PoncedeLeón G, Domínguez–Domínguez O, Doadrio I: A new species of the genus Algansea Girard, 1856 (Actinopterygii: Cyprinidae) from the Ameca River basin, in Central Mexico.Revista Mexicana de Biodiversidad8: 340.
Howes GJ: Systematics and biogeography: an overview.Cyprinids fishes: Systematics, biology and exploitation(Edited by: Winfield IJ, Nelson JS). London: Chapman & Hall 1991, 1–33.
Simons AM, Berendzen PB, Mayden RL: Molecular systematics of North American phoxinin genera (Actinopterygii: Cyprinidae) inferred from mitochondrial 12S and 16S ribosomal RNA sequences.Zoological Journal of the Linnean Society 2003, 139: 63–80.View Article
Lydeard C, Roe K: The phylogenetic utility of the mitochondrial cytochrome b gene for inferring relationships among Actinopterigian fishes.Molecular systematics of fishes(Edited by: Kocher TD, Stepien CA). New York Academic Press 1997, 285–302.
Chow S, Hazama K: Universal PCR primers for S7 ribosomal protein gene introns in fish.Molecular Ecology 1998, 7: 1247–1263.View Article
Farias IP, Ortí G, Sampaio I, Schneider H, Meyer A: The cytochrome b gene as a phylogenetic marker: the limits of resolution for analyzing relationships among cichlid fishes.Journal of Molecular Evolution 2001, 53: 89–103.PubMed
Bernardi G, Bucciarelli G, Costagliola D, Robertson DR, Heiser JB: Evolution of coral reef fish Thalassoma spp. (Labridae). 1. Molecular phylogeny and biogeography.Marine Biology 2004, 144: 369–375.View Article
Johnson JB, Dowling TE, Belk MC: Neglected taxonomy of rare desert fishes: congruent evidence for two species of leatherside chub.Systematic Zoology 2004, 53: 841–855.
Posada D, Crandall KA: Modeltest: testing the model of DNA substitution.Bioinformatics 1998, 14: 817–818.View ArticlePubMed
Bielawski JP, Gold JR: Phylogenetic Relationships of Cyprinid Fishes in Subgenus Notropis Inferred from Nucleotide Sequences of the Mitochondrially Encoded Cytochrome b Gene.Copeia 2001, 656–667.
Orti G, Petry P, Porto JIR, Jégu M, Meyer A: Patterns of nucleotide change in mitochondrial ribosomal RNA genes and the phylogeny of piranhas.Journal of Molecular Evolution 1996, 42: 169–182.View ArticlePubMed
Lavoué S, Sullivan JP, Hopkins CD: Phylogenetic utility of the first two introns of the S7 ribosomal protein gene in african electric fishes (Mormyroidea: Teleostei) and congruence with other molecular markers.Biological Journal of the Linnean Society 2003, 78: 273–292.View Article
Jensen RJ, Barbour CD: A phylogenetic reconstruction of the mexican cyprinid fish genus Algansea.Systematic Zoology 1981, 30: 41–57.View Article
Minckley WL, Hendrickson DA, Bond CE: Geography of western North American freshwater fishes: description and relationships to intracontinental tectonism.Zoogeography of North American Freshwater fishes(Edited by: Hocutt CH, Wiley EO). New York: Wiley–Interscience Publications 1986, 519–613.
Sullivan JP, Lavoué S, Hopkins CD: Discovery and phylogenetic analysis of a riverine species flock of african electric fishes (Mormyridae: Teleostei).Evolution 2002, 56: 597–616.PubMed
Goncalves H, Martínez–Solano I, Ferrand N, García–París M: Conflicting phylogenetic signal of nuclear vs mitocondrial DNA markers in midwife toads (Anura, Discoglossidae, Alytes): Depp coalescence or ancestral hybridization?Molecular Phylogenetics and Evolution 2007, 44: 494–500.View ArticlePubMed
Holland BR, Benthin S, Lockhart PJ, Moulton V, Huber KT: Using supernetworks to distinguish hybridization from lineage–sorting.BMC Evolutionary Biology 2008, 8: 202.View ArticlePubMed
Annilo T, Jelina J, Pata I, Metspalu : Isolation and characterization of the Mouse ribosomal proteína S7 gene.Biochemistry and Molecular Biology International 1998, 46: 287–295.PubMed
Chow S, Scholey VP, Nakazawa A, Margulies JB, Wexler JB, Olson RJ, Hazama K: Direct Evidence for Mendelian Inheritance of the Variations in the Ribosomal Protein Gene Introns inYellowfin Tuna (Thunnus albacares).Marine Biotechnology 2001, 3: 22–26.View ArticlePubMed
He S, Mayden RL, Wang X, Wang W, LTang K, Chen WJ, Chen Y: Molecular phylogenetics of the family Cyprinidae (Actinopterygii: Cypriniformes) as evidenced by sequence variation in the first intron of S7 ribosomal protein–coding gene: Further evidence from a nuclear gene of the systematic chaos in the family.Molecular Phylogenetics and Evolution 2008, 46: 818–829.View ArticlePubMed
Holder MT, Anderson JA, Holloway AK: Difficulties in detecting hybridization.Systematic Biology 2001, 50: 978–982.View ArticlePubMed
Buckley TR, Cordeiro M, Marshall DC, Simon C: Differentiating between hypotheses of lineage sorting and introgression in New Zealand alpine cicadas (Maoricicada Dugdale).Systematic Biology 2006, 55: 411–425.View ArticlePubMed
Miller DL, Behnke RJ: Two new intergeneric cyprinid hybrids from the Bonneville Basin, Utah.Copeia 1985, 1985: 509–515.View Article
Dowling TE, Smith GR, Brown WM: Reproductive isolation and introgression between Notropis cornutus and Notropis chrysocephalus (family Cyprinidae): comparson of morphology allozymes, and mitochondrial DNA.Evolution 1989, 43: 620–634.View Article
DeMarias BD, Dowling TE, Douglas ME, Minckley WL, Marsh PC: Origin of Gila seminuda (Teleostei:Cyprinidae) through introgressive hybridization: Implications for evolution and conservation.Proceedings of the National Academy of Sciences of the United States of America 1992, 89: 2747–2751.View Article
Dowling TE, DeMarias BD: Evolution significance of introgressive hybridization in cyprinid fishes.Nature 1993, 362: 444–446.View Article
Bostrom MA, Collete BB, Luckhurst BE, Reece KS, Graves JE: Hybridization between two serranids, the coney (Cephalopholis fulva) and the creole–fish (Paranthias furcifer), at Bermuda.Fishery Bulletin 2002, 100: 651–661.
Rosenthal GG, Rosa–Reyna XF, Kazians S, Stephens MJ, Morizot DC, Ryan MJ, León JGd: Dissolution of sexual signal complexes in a hybrid zone between the swordtails Xiphophorus birchmanni and Xiphophorus malinche (Poeciliidae).Copeia 2003, 2003: 299–307.View Article
Moore WS: Inferring phylogenies from mtDNA variation: mitochondrial–gene trees versus nuclear–gene trees.Evolution 1995, 49: 718–726.View Article
Henry CD, Aranda–Gómez JJ: Plate interactions control middlelate Miocene, proto–Gulf and Basin and Range extension in the southern Basin and Range.Tectonophysics 2000, 318: 1–26.View Article
Ferrari L, López–Martínez M, Rosas–Elguera J: Ingnimbrite flare–up and deformation in the southern Sierra Madre Occidental, western México. Implications for the late subduction history of Farallon plate.Tectonics 2002, 21: 1–24.View Article
Israde–Alcántara I, Wade M, Garduño–Monroy VH, Barron J: Estratigrafía y encuadramiento geodinámico de las cuencas lacustres del centro de México.Unión Mexicana de estudios del Cuaternario Fondo de Cultura Económica, in press.
Barbour CD, Miller RR: Diversification in the mexican cyprinid fish Algansea monticola (Pisces: Cyprinidae), with description of a new subspecies.Copeia 1994, 662–676.
Ferrari L, Rosas–Elguera J: Late Miocene to Quaternary extension an the northern boundary of the Jalisco block, western Mexico: The Tepic–Zacoalco rift revised.Geological Society of America Special Paper 1999, 334: 1–23.
Israde–Alcántara I, Garduño–Monroy VH: Lacustrine record in a volcanic intra–arc setting: the evolution of Late Neogene Cuitzeo basin system (central–western Mexico, Michoacan).Palaeogeography, Palaeoclimatology and Palaeoecology 1999, 151: 209–227.View Article
Ferrari L, López–Martínez M, Aguirre–Díaz G, Carrasco–Núñez G: Space–time patterns of Cenozoic arc volcanism in central Mexico: from the sierra madre occidental to the Mexican Volcanic Belt.Geology 1999, 27: 303–306.View Article
Tamayo LT, West RC: The hydrology of Middle America.Handbook of Middle America indians(Edited by: West RC). Texas: University of Texas Press 1964, 84–121.
Ferrusquia–Villafranca I: Geología de México: una synopsis.Diversisdad Biológica de México: orígenes y distribución(Edited by: Rmamoorthy TP, Bye R, Lot A, Fa J). México DF: Instituto de biología, Universidad Nacional Autónoma de México 1998, 1–107.
Aranda–Gómez JJ, Henry CD, Luhr JF: Evolución tectomagmática post–paleocénica de la Sierra Madre Occidental y de la porción meriodonal de la provincia tectónica de Cuencas y Sierras, México.Boletín de la Sociedad Geológica Mexicana 2000, 53: 57–59.
Gesundheit P, Macias–García C: Biogeografía cladísta de la familia goodeidae (Cyprinodontiformes).Una prespectiva latinoamericana de la biogeografía(Edited by: Llorente J). México DF: Comisión Nacional para el Conocimiento y Uso de la Biodiverisidad 2006, 319–337.
Mateos M: Comparative phylogeography of livebearig fishes in the genera Poeciliopsis and Poecilia (Poeciliidae: Cyprinodontiformes) in central Mexico.Journal of Biogeography 2005, 32: 775–780.View Article
Pérez–PoncedeLeón G: Biodiversity and biogeographic patterns in the Mesa Central of México: insights from host–parasities systems.Journal of Parasitology 2003, 89: 126–133.
Aguilar–Aguilar R, Contreras–Medina R, Salgado–Maldonado G: Parsimony analysis of endemicity (PAE) of Mexican hydrological basins based on helminth parasites of freshwater fishes.Journal of Biogeography 2003, 30: 1861–1872.View Article
Pérez–PoncedeLeón G: processes AC BohpoffiMtsfpa: Biogeography of helminth parasites of freshwater fishes in Mexico: the search for patterns and processes.Journal of Biogeography 2005, 32: 645–659.View Article
Mejia–Madrid HM, Vázquez–Domínguez E, Pérez–PoncedeLeón G: Phylogeography and freshwater basins in central Mexico: recent history as revealed by the fish parasite Rhabdochona lichtenfelsi (Nematoda).Journal of Biogeography 2007, 34: 787–801.View Article
Rosas–Valdéz R, Domíngue–Domínguez O, Choudhury A, Pérez–PoncedeLeón G: Helminth parasites of the Balsas catfish Ictalurus balsanus in several localities of the Balsas River Drainage, Mexico: Species composition and biogeographical affinities.Comparative parasitology 2007, 74: 204–210.View Article
Contant R: Observations on Garter Snakes of the Thamnophis eques Complex in the Lakes of Mexico's Transvolcanic Belt, with Descriptions of New Taxa.American Museum Novitates, American Museum of Natural Histroy, NY 2003, 3406: 1–64.
Mulcahy DG, Mendelson JR III: Phylogeography and speciation of the morphologically variable, widespread species Bufo valliceps, based on molecular evidence from mtDNA.Molecular Phylogenetics and Evolution 2000, 17: 173–189.View ArticlePubMed
Zaldívar–Riverón A, León–Regagnon V, Nieto–MontesdeOca A: Phylogeny of the Mexican coastal leopard frogs of the Rana berlandieri group based on mtDNA sequences.Molecular Phylogenetics and Evolution 2004, 30: 38–49.View ArticlePubMed
Weisrock DW, Shaffer HB, Storz BL, Storz SR, Voss SR: Multiple nuclear gene sequences identify pfulogenetic species boundaries in the rapidly radiating clade of Mexican ambystomatid salamanders.Molecular Ecology 2006, 15: 2489–2503.View ArticlePubMed
Chan MD: Fish ecomorphology: predicting habitat preferences of stream fishes from their body shape.PhD thesis Faculty of the Virginia Polytechnic Institute and the State University 2001.
Lyons J, González–Hernández G, Soto–Galera E, Guzmán–Arroyo M: Decline of freshwater fishes and fisheries in selected drainages of west–central Mexico.Fisheries 1998, 23: 10–18.View Article
Sambrook J, Fritsch EF, Maniatis T: Molecular cloning: a laboratory manual New York: Cold Spring Harbor Laboratory Press 1989.
Machordom A, Doadrio I: Evidence of a Cenozoic Betic–Kabilian connection based on freshwater fish phylogeography (Luciobarbus, Cyprinidae).Molecular Phylogenetics and Evolution 2001, 18: 252–263.View ArticlePubMed
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.Nucleic Acids Researh 1997, 24: 4876–4882.View Article
Kumar S, Tamura K, Nei M: MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment.Briefings in Bioinformatics 2004, 5: 150–163.View ArticlePubMed
Swofford DL: Phylogenetic analysis using parsimony (* and other methods) V 4.0b 10 Sunderland MA, Sinauer Associates 2004.
Maddison DR, WPMeC: McClade. Analysis of phylogeny and character evolution v.4.05 OS X Sunderland MA, Sinauer Associates 2002.
Dowling TE, Tibbets CA, Minckley WL, Smith GR: Evolutionary relatioships of the plagopterins (Teleostei:Cyprinidae) from cytochrome b sequences.Copeia 2002, 2002: 665–678.View Article
Doadrio I, Carmona JA: Testing freshwater Lago Mare dispersal theory on the phylogeny relationships of Iberiancyprinid genera, Chondrostoma and Squalius.Graellsia 2003, 59: 457–473.
Durand JD, Bianco PG, Laroche J, Gilles A: Insight into the origin of endemic Mediterranean ichthyofauna – Phylogeography of Chondrostoma genus (Teleostean, Cyprinidae).Journal of Heredity 2003, 94: 315–328.View ArticlePubMed
Doadrio I, Carmona JA: Phylogenetic relationships of the genus Chondrostoma using cytochrome b sequences.Molecular Phylogenetics and Evolution 2004, 13: 2807–2817.
Robalo JI, Sousa–Santos CS, Carvalho–Almada V, Doadrio I: Paleobiogeography of Two Iberian Endemic Cyprinid Fishes (Chondrostoma arcasii–Chondrostoma macrolepidotus)Inferred from Mitochondrial DNA Sequence Data.Journal of Heredity 2006, 97: 143–149.View ArticlePubMed
Robalo JI, Carvalho–Almada V, Levya A, Doadrio I: Re–examination and phylogeny of the genus Chondrostoma based on mitochondrial and nuclear data and the definition of 5 new genera.Molecular Phylogenetics and Evolution 2006, 42: 362–372.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.