Evolutionary history of the fish genus Astyanax Baird & Girard (1854) (Actinopterygii, Characidae) in Mesoamerica reveals multiple morphological homoplasies
© Ornelas-García et al; licensee BioMed Central Ltd. 2008
Received: 23 May 2008
Accepted: 22 December 2008
Published: 22 December 2008
Mesoamerica is one of the world's most complex biogeographical regions, mostly due to its complex geological history. This complexity has led to interesting biogeographical processes that have resulted in the current diversity and distribution of fauna in the region. The fish genus Astyanax represents a useful model to assess biogeographical hypotheses due to it being one of the most diverse and widely distributed freshwater fish species in the New World. We used mitochondrial and nuclear DNA to evaluate phylogenetic relationships within the genus in Mesoamerica, and to develop historical biogeographical hypotheses to explain its current distribution.
Analysis of the entire mitochondrial cytochrome b (Cytb) gene in 208 individuals from 147 localities and of a subset of individuals for three mitochondrial genes (Cytb, 16 S, and COI) and a single nuclear gene (RAG1) yielded similar topologies, recovering six major groups with significant phylogeographic structure. Populations from North America and Upper Central America formed a monophyletic group, while Middle Central America showed evidence of rapid radiation with incompletely resolved relationships. Lower Central America lineages showed a fragmented structure, with geographically restricted taxa showing high levels of molecular divergence. All Bramocharax samples grouped with their sympatric Astyanax lineages (in some cases even with allopatric Astyanax populations), with less than 1% divergence between them. These results suggest a homoplasic nature to the trophic specializations associated with Bramocharax ecomorphs, which seem to have arisen independently in different Astyanax lineages. We observed higher taxonomic diversity compared to previous phylogenetic studies of the Astyanax genus. Colonization of Mesoamerica by Astyanax before the final closure of the Isthmus of Panama (3.3 Mya) explains the deep level of divergence detected in Lower Central America. The colonization of Upper Mesoamerica apparently occurred by two independent routes, with lineage turnover over a large part of the region.
Our results support multiple, independent origins of morphological traits in Astyanax, whereby the morphotype associated with Bramocharax represents a recurrent trophic adaptation. Molecular clock estimates indicate that Astyanax was present in Mesoamerica during the Miocene (~8 Mya), which implies the existence of an incipient land-bridge connecting South America and Central America before the final closure of the Isthmus of Panama (~3.3 Mya).
Mesoamerica is one of the most complex biogeographical areas in the world [1–5]. This complexity reflects the confluence of Neotropical and Nearctic biotas and a long history of geological activity, stretching from the Miocene to the present, during which movements of the Cocos, North American, Pacific and Caribbean Plates [6, 7] created barriers and land-bridges that have affected the distribution of freshwater fishes [8–13]. For example, the Pliocene (~3.3 Mya) closure of the Panama Strait has been postulated to be one of the most important causes of faunal interchange between Neartic and Neotropical regions . Climatic changes have also been invoked to explain the distribution of Mesoamerican fish fauna . Distinguishing between climatic and geological effects requires information on phylogeny and species boundaries in a diversity of taxa [8, 9, 16–19].
Special attention has been devoted to understanding the number and timing of colonizations of Mesoamerica by freshwater fishes from South America: a topic which remains somewhat controversial [8, 13, 20–22]. The most widely accepted theories support two waves of colonization: 1) an ancient episode (70 – 80 Mya) through a proto-Antillean arc and 2) a more recent, Cenozoic episode via the Antillean islands and/or a continental corridor [14, 23, 24]. Molecular data suggest colonization of Mesoamerica by primary freshwater fishes about 4–7 Mya [8, 12, 13, 22]. This is incongruent with the geological data, which does not support the existence of a continental land-bridge before the closure of the Panama Strait (3.3 Mya). Older colonization events have been postulated for secondary freshwater fishes: for example, Early to Mid Miocene (12.7–23 Mya) colonization for Synbranchidae , 14–24 Mya for heroinid cichlids (10 Mya for Mesoamerican lineages) and 18.4–20 Mya for rivulids .
The absence of primary freshwater genera (e.g., Hypopomus, Pimelodella, Rhamdia and Roeboides) from Mesoamerica and the Antillean islands argues against an ancient colonization route through a proto-Antillean arc. Instead, it supports a colonization route through an incipient land-bridge formed during the gradual uplifting of the Panama Isthmus [8, 12], over a time span of 3–20 Mya, combined with changes in sea level . This is supported by molecular studies of arthropods , amphibians , and marine geminate species pairs on either side of the Panamanian Isthmus .
Support for a more recent colonization of Mesoamerica by primary freshwater fish through the Panama Strait comes from phylogeographical studies of Characids (e.g., Brycon, Bryconamericus, Eretmobrycon, and Cyphocharax). These studies indicate multiple waves of rapid expansion from South America during the Pliocene ~3.3 Mya. .
The genus Astyanax provides an ideal model to investigate the relative importance of vicariance and dispersal on biogeographical patterns. This is partly because it is widely distributed across the region , and because its dispersal is confined to freshwater routes and dependent, therefore, on the formation of land-bridges.
Characiforms are generally assumed to have a Gondwanan (South American) origin [30–32], as supported by the fossil record , so the presence of Characidae in Northern America is viewed as a consequence of dispersal.
Astyanax comprises more than 107 recognized species and is, together with Hyphesobrycon (105 species), the largest and most diverse characiform genus [34, 35]. Moreover, Astyanax has the widest distribution of American characids, being found from the Nearctic (Colorado River in Texas and New Mexico) to the Neotropics (Negro River in Patagonia) .
Previous phylogenetic studies [36, 37] of the biogeography of Astyanax used a small number of samples from Mexico, Belize and Guatemala, and did not find geographical congruence for some of the groups recovered (i.e. Yucatan and Belizean populations were not the most closely related despite their geographical proximity). Furthermore, conspecific cave and nearest surface populations formed two separate lineages, in agreement with an earlier study of the genus . This was attributed to at least two separate colonizations of Mesoamerica from South America during the Pleistocene. Estimated colonization times based on the cytochrome b gene were 1.8 and 4.5 Mya (3.1 Mya), with an estimated divergence rate of 1.5% per pairwise comparison per million years, which coincides with the closure of the Panama Strait (3.3 Mya) . However, incomplete sampling (only few samples were included from upper Central America and Mexico) could lead to erroneous interpretations. This study provides a phylogeographical analysis based on a comprehensive distribution-wide sampling regime and more extensive sampling, and thus should provide new insights into the evolutionary history of the genus.
The genus Astyanax is characterised by high phenotypic plasticity and a capacity to adapt to diverse habitats [36, 38–41]. There is clear evidence of extremely rapid adaptations of fish to new habitats and environments, with ecological specialization and morphological differentiation, generally in accordance with genetic divergence [15, 42–44]. Considerable attention has been given to the evolution of developmental mechanisms and adaptation to cave environments [36, 38, 39], but less attention has been given to other habitat associated morphological plasticity. In this regard, Bramocharax, which is sympatric with Astyanax, is characterized by conspicuous trophic specializations, including differences in the number of premaxillary teeth, the presence of diastemas on the maxillary teeth, as well as differences in the shape and number of cuspids on the premaxillary, maxillary and dentary teeth, with some species (B. caballeroi and B.baileyi) having intermediate states between the morphotypes of Astyanax and Bramocharax species [45–47].
In this study we used mitochondrial and nuclear DNA sequences to develop a robust phylogenetic hypothesis for Astyanax and Bramocharax. This allows us to test biogeographical hypotheses for the Mesoamerican fish fauna, including the relative importance of historical geology and climatic factors.
Three mitochondrial (Cytb, COI and 16 S) genes and one nuclear gene (RAG-1) were sequenced, giving a total of 3862 characters (2350 mitochondrial and 1512 nuclear).
Primers and PCR conditions.
RAG1 (a, s)
RAG9 (a, s)
Cytochrome oxidase I (COI)
FISHF1 (a, s)
Cytochrome b (Cytb)
Glu- F (a, s)
Thr- R (a, s)
CbCHR (a, s)
16SAR (a, s)
All analyses supported the polyphyly of Bramocharax, with species of Bramocharax being sister groups to different clades of Astyanax (Figure 1), making Astyanax paraphyletic.
Phylogenetic Performance of each gene.
Geographic structuring within the major phylogenetic groupings
GROUP I (Mexico and Upper Central America)
Four main clades were recovered from Group I. Clade I comprised most of the Mexican and Upper Central American (Guatemala and Belize) populations. Clades II-IV represented fewer populations with a patchy distribution over the range of Clade I.
Within Clade I we found geographical structure corresponding to the following four lineages: Lineage Id, from the Chiapas region of Mexico, was sister to a clade comprising Lineage Ie from the Candelaria region and a pair of sister lineages (1a and 1b) that occupy a wider region from the Yucatan Peninsula to the Bravo-Conchos basin.
Lineage Ib contains Astyanax and Bramocharax from southern TMVB and Belize. It is subdivided into four sublineages. It ranges from the Media Luna Lagoon (Panuco basin, Mexico) to the Mopan Basin (Belize) on the Atlantic slope, and from the Armeria – Coahuayana Basin to the Balsas Basin (both in Mexico) on the Pacific slope. This group is a good example of morphological plasticity with low levels of genetic differentiation (see further details in Additional file 2).
The first sublineage included Astyanax fasciatus from Puente Nacional to Grijalva – Usumacinta basins (including Papaloapan and Coatzacoalcos basins) on the Atlantic slope, and populations from the type localities of A. armandoi and A. altior (Palenque in the Grijalva – Usumacinta basin and Cenote of Noc-Ac in the Yucatan Peninsula, respectively). Both of these are junior synonyms of A. fasciatus [35, 48]. We also found shared haplotypes or low divergence (see Additional file 2) between Bramocharax caballeroi and sympatric A. fasciatus from Lake Catemaco (see region C, Los Tuxtlas, Figure 4). Additionally these distances were lower than those observed within Astyanax populations.
Lineage Ic, (region D, Figure 3) comprised populations from the Atlantic slope: Astyanax fasciatus (a synonym of A. mexicanus sensu Lima et al. ) from the type locality of A. petenensis (Peten – Itza Lake, Yucatan Peninsula) and A. fasciatus from the Candelaria karst region in Guatemala. This lineage included Bramocharax dorioni from its type locality (Semococh river). Similar to other instances of sympatry between Bramocharax and Astyanax, these two morphotypes showed low levels of genetic differentiation (less 0.5%, see Additional file 3).
Lineage Id (Chiapas) grouped Pacific slope populations of A. fasciatus (Figure 1) from the Pichoacan basin (Oaxaca, Mexico) to the El Jococal Lagoon (El Salvador), including the Guatemalan coast, with very low genetic divergence within the clade.
Clade II comprised two lineages. (Ie and If; mean divergence of = 2.13% ± 1.21). Lineage Ie (referred to as "Sabinos-Aguanaval-Mezquital") grouped troglobitic morphotypes (A. jordani) from the Piedras, Tinaja, La Curva, and Sabinos caves (type locality of A. hubbsi, synonym of A. jordani sensu Lima et al. ) with the surface-dwelling populations of A. mexicanus from the Mezquital and Nazas – Aguanaval basins, with a low level of differentiation between surface and cave populations (mean of = 1.4% ± 0.9).
Lineage If included populations of A. mexicanus from the Rascon valley and the Panuco basin in Tamposa. This lineage was highly divergent with regard to the rest of the lineages in this group, and those populations closest geographically (see Additional file 2)
Clade III grouped Atlantic slope populations of A. fasciatus from the La Palma and Maquinas basins (see Los Tuxtlas region, Figure 3) with those from the upper Polochic basin (Cahabon river) and the Grijalva-Usumacinta basin.
Clade IV contained populations from the Montebello Lagoons (Montebello) in south-eastern Mexico. This group clustered with the Maquinas and La Palma populations in the analysis of Cytb alone, but this clustering was not supported by the combined nuclear and mitochondrial analysis (Figures 2 and 4).
GROUP II (Middle Central America)
Our analyses did not resolve relationships among the five main clades recovered within Group II. This group occurs widely over Middle Central America, ranging from Belize to Nicaragua and Costa Rica (E, Figure 3).
Clade V included A. fasciatus from the Macal Basin (Belize). It is highly divergent with respect to the other clades from Group II (D K81uf = 3.6% ± 1.39).
Clade VII comprised two lineages (IIa and IIb). Lineage IIa included Astyanax nicaraguensis and Bramocharax bransfordii. As in other clades containing both morphotypes there were low levels of genetic differentiation (D K81uf = 0.9% ± 0.61). This lineage occurs from the Jigüina Basin to the Sarapiqui basin (Nicaragua) on the Atlantic slope, and included populations from the great Lakes of Nicaragua (Managua and Nicaragua, region E, figures 1 and 4) and the Senacapa and Brito basins (Nicaragua) on the Pacific slope. Lineage IIb included A. fasciatus from Lake Nicaragua and the San Juan basin (Nicaragua), as well as the Barbilla and Reventaza basins on the Atlantic slope of Costa Rica.
Clade VIII included A. fasciatus from the Pacific tributaries of Nicaragua and A. nasutus from the Managua Lake Basin.
Clade IX included A. fasciatus from a very restricted area comprising the Ciruelas and Tempisque basins on the Pacific slope of Costa Rica ("Ciruelas-Tempisque", Figure 5).
Lower Central America
Four major groups (groups III-VI) in Lower Central America were characterized by more restricted ranges relative to Groups I and II. Group III comprised A. fasciatus populations from the Ciruelas and Chires basins on the Pacific Slope of Costa Rica (figures 1 and 2). Two highly divergent (D K81uf = 2.5%) haplotypes (corresponding to Groups II and III) were found in sympatry in the Ciruelas basin (figures 1 and 2).
Group IV included a single and well-differentiated population of A. orthodus from the Sixaola basin on the Atlantic slope of Costa Rica. Group V contained Pacific slope A. fasciatus populations from the Puntarenas basin (Costa Rica) to the Lagarto basin on the Panama-Costa Rica border. Group VI included A. fasciatus from the Chagres region on the Atlantic slope of Panama.
Discussion and conclusion
Systematics of the genera Astyanax and Bramocharax
Our analyses do not support the previously proposed monophyly of Bramocharax based on morphological analyses [45–47]. Moreover, Bramocharax specimens were present in two of the seven major Astyanax clades, with low levels of genetic differentiation when both morphotypes were found in sympatry (such as in lake Catemaco where it was possible to find haplotypes shared between individuals from both genera). Differentiation was equally reduced between allopatric populations of Bramocharax and Astyanax.
The genus Astyanax has been considered to be monophyletic in Mesoamerica , but polyphyletic in South America  on the basis of molecular analyses. Our results support the monophyly of Mesoamerican Astyanax only if we consider Bramocharax species to be morphotypes of Astyanax within the range of its phenotypic plasticity. This hypothesis is supported by the low genetic divergence between specimens of Bramocharax and Astyanax, and the evidence of recurrent evolution of the Bramocharax morphotype within Astyanax (Figure 3). This morphotype is associated with lacustrine habitats, suggesting that its recurrent evolution is a result of morphological convergence to similar ecological factors; similar patterns have been shown in other freshwater fishes [37, 38, 41, 45, 49]. If the "recurrent convergence" hypothesis is considered to be correct, then the taxonomy of Bramocharax needs to be revised and the evolutionary mechanisms giving rise to these morphological homoplasies need further investigation. Our analyses question the taxonomic utility of trophic characters (e.g., teeth shape or jaw modification), as previously done by Rosen  on the basis of intermediate morphological states between Astyanax fasciatus and Bramocharax baileyi.
Further incidences of morphological convergence were found in troglobitic morphotypes of Astyanax jordani (this has been noted by previous authors [37, 38, 50]), providing further evidence of independent (at least two different times, see Figure 6) adaptation to troglobitic habitats. The presence of recurrent morphological convergence in Astyanax [50, 51] makes the delimitation of species and genera difficult. Thus the absence of congruence between phylogenetic relationships uncovered in this study and previous taxonomic classifications for Astyanax and Bramocharax from Mesoamerica [35, 37, 52] is not surprising. In addition, our results are not in agreement with the idea that Astyanax (including samples from Mexico and Upper Central America) is a single species (i.e., A. fasciatus) as has previously proposed .
Although not a main goal of this study, we propose a provisional taxonomic nomenclature for Astyanax populations from Mesoamerica (see Additional file 3). The nomenclature proposed is based on well-defined monophyletic groups, high genetic divergences with Cytb (>2% K81uf), and in agreement with geographical distributions. In ascribing species names we gave priority to previous species descriptions and diagnostic morphological traits. Where monophyletic lineages could not be assigned to a valid species name, they were assigned to their own monophyletic group as Astyanax sp.
Genetic and Time Divergences
Genetic distances in percentage among major groups (below diagonal, uncorrected p sequence divergences; above diagonal, ML distances K81uf)
We found a pattern of north-south phylogeographical structuring. The major phylogenetic groups were mostly non-overlapping, with the exception of Groups I and II, which overlap in the upper part of the Polochic basin of Guatemala, and Groups II and III, which overlap in the Ciruelas basin of Costa Rica. This north-south pattern is similar to that reported for other freshwater fishes [12, 13, 22]. We explain the observation of sympatric lineages of Astyanax in terms of niche overlap and lineage turnover, similar to that proposed in biogeographical models for other characins in Mesoamerica .
Summary of the timings of the main geological events in Mesoamerica based on freshwater fauna.
Roeboides, Hypopomus and Pimelodella* 
4–7 Mya Mio-Pliocene
3 Mya Pliocene
2.5–2.9 Mya guatemalensis clade
1.5% (Cytb) HKY85
6.5–5.6 Mya laticauda clade
Synbranchus and Ophisternon** 
1.3% (ATP6&8) K2P
12.7 – 23 Mya Miocene
1.5% (Cytb) TrN+G
Colonization of the Mesoamerica from South America
10 Mya Miocene
18–20 Mya Miocene
Brycon, Bryconamericus, Eretmobrycon and Cyphocharax * 
< 3.1 Mya Plio-Pleistocene
MtDNA and nuclear genes
Mateos (2002) calibration
Before to Panama closure Cretaceous-Miocene
8–16 Mya to 2.8–6.4 Mya Mio-Pliocene
Dispersal hypothesis on the origin of genus Astyanaxin Mesoamerica
We accept the widely held hypothesis of a South American origin for Astyanax and other Central American characids [15, 30, 33]. This is supported by the observation that Lower Central America lineages were most closely related to South American samples from Brazil and Argentina. Considering the widely used Cytb calibration rate for fish (1.09%/my HKY distances) [22, 38, 54] and our mean rate of 0.8%/my K81uf distances from R8s (Figure 6 and Table 4), levels of divergence for populations of South and Central America imply a period of Mesoamerican colonization/expansion of Astyanax from South America about 7.8–8.1 Mya, before the final uplift of the Isthmus of Panama ~3.3 Mya [7, 20, 56]. The inclusion of more Astyanax samples from both sides of the Sierra de Perija and Merida Andes in further studies could improve this dating scenario.
The colonization of Central America prior to Late Cenozoic closure of the Panama Strait is incongruent with the geological data, and with other studies of characid genera (Brycon, Bryconamericus, Eretmobrycon and Cyphocharax) , including a previous study of Astyanax , all of which propose that closure of the strait ~3.3 Mya provided the first opportunity for colonization of Central America from South America.
An earlier colonization of Mesoamerica has been proposed for other freshwater fishes (Table 3) [[8, 12], , ]. For example, ancient colonization events have been proposed for the family Poeciliidae (Cretaceous – Rosen model )  and the Cichlidae, Rivulidae and Synbranchidae families (Miocene – GAARlandia model as proposed by Iturralde-Vinent and MacPhee ) [8, 12, 13, 22]. However in contrast with Astyanax, these families also occur in the Caribbean islands, and as secondary fishes could have crossed through a shallow passage between South America and Central America during the Miocene [13, 20, 22].
Divergence times similar to those found in this study have been reported in molecular studies of primary freshwater fauna [8, 12]. For example, the Bermingham and Martin model  proposes a colonization of Mesoamerica between 4 and 7 Mya, prior to the final closure of the Panama Isthmus (~3.3 Mya), based on comparative phylogeography of three genera [Roeboides (Characidae), Pimelodella (Pimelodidae) and Hypopomus (Hypopomidae)] Moreover, Rhamdia seems to have colonized Mesoamerica in two different waves, one prior to the closure of the Panama Strait (6.5 to 5.6 Mya, R. laticauda group) . These estimates are in agreement with our study, and coincident with other primary and secondary freshwater fish fauna [12, 13, 22], as well as divergence times for invertebrates (9 Mya in pseudoscorpions)  and benthic foraminifera fossils (8 Mya) .
We found evidence of a unique biogeographical pattern involving multiple waves of expansion of Group I (Clades II – V) Astyanax in the upper part of the Polochic-Motagua fault (Mexico and Chiapas region). These clades have a restricted distribution overlapping that of Clade I, and in general occupy relatively stable ecological environments (springs or lakes), which can be less affected by climate change. Niche overlap and lineage turnover could explain this pattern, except in stable habitats where two lineages are found in sympatry (lineages Ia of Clade I and the lineage of Clade II in the Huasteca region, and in the Mezquital and Nazas-Aguanaval basin).
Main vicariant events in Astyanaxpopulations from Mesoamerica
The vicariance events involving Astyanax in Mesoamerica occurred during the Plio-Miocene (4–8 Mya), occurring earlier in Lower Central America (Panama and mainly Costa Rica) than in Central and Upper Mesoamerica.
A pattern of restricted geographic ranges in Lower Central America (Groups IV-VI) supports pronounced geographical fragmentation as a consequence of tectonics movements [7, 20], which eventually resulted in closure of the Panama Strait ~3.3 Mya. Bermingham and Martin  have implicated multiple range fragmentation during the Miocene in patterns of diversity in other taxa of the primary freshwater fauna. With our data, five main vicariant events were identified for Lower Central America (Figure 6). These are related to changes in eustatic sea level (5–8 Mya)  and to the formation of inter-oceanic biogeographical barriers  during the Middle-Late Miocene (8 Mya).
In Upper Central America, the main volcanic activity was produced by the Trans-Mexican Volcanic Belt (TMVB). This region was affected by periods of intense geological activity between 3 and 12 Mya, with some volcanic activity still occurring today [53, 59, 60]. The geographic structuring evident in Clade I of Astyanax indicates that the TMVB formed an effective geographic barrier during its development during the late Miocene 4 – 6 Mya (Figure 6). This date is in agreement with the geology of the region and previous studies of several groups of vertebrates [1, 5, 10, 11, 61, 62].
Other biogeographical patterns
Other biogeographical patterns were obtained in Lower and Middle Central America. These cannot be explained by geological barriers, but are in accordance with the main biogeographical regions proposed for other freshwater fishes [15, 22, 63].
While our results contrast somewhat with the Mesoamerican faunal regions recognized by Bussing , they do support his Isthmian Region, but with a fragmented pattern similar to that reported for other primary freshwater fishes [8, 12, 15]. Furthermore, we found that some Belize and Guatemala Astyanax clades (Clades VII and VIII in Figure 5) were joined to the Central America group (group II), a pattern shared with other Mesoamerican Cichlids . Our Chiapas-Nicaraguan lineages (Lineage Id in Figure 4) did not reach the Pacific cost of Costa Rica, but instead have their southern limits in El Salvador. The distribution of lineages from San Juan Region differs from Bussing's proposal , occurring from the Barbilla basin (Costa Rica) to Belize (Atlantic slope) and on the Pacific slope from the Rio Grande Basin to the Ciruelas basin in the Nicoya Gulf (Costa Rica). This pattern has been previously reported by other authors [22, 63].
Finally, we dated the separation of Groups I and II to about 6 and 7.8 Mya (Table 4). These coincide geographically with Polochic-Motagua Fault , reported as a transition region for other freshwater fish groups [22, 13]. In addition, we observed the presence of two well differentiated Lineages of Astyanax (Clade IV in Group I and Clade V in Group II) in sympatry in the Polochic basin, a finding that has not been reported for other characids . We explain this pattern in terms of river capture whereby the Cahabon tributary was diverted to the Polochic river as a consequence of tectonic activity (Sierra de Chiapas ), while separating the Cahabon river from the Grijalva-Usumacinta (Clade IV in Group I). This has been proposed for Rhamdia .
Tissue collection and DNA extractions
A total of 208 specimens of the Astyanax and Bramocharax from 141 localities from Panama to the Mexican-USA border (Additional file 3; Figure 3) were analyzed, corresponding to 10 species of Astyanax: A. aeneus, A. altior, A. armandoi, A. fasciatus, A. jordani, A. mexicanus, A. nasutus, A. nicaraguensis, A. orthodus and A. petenensis; and three species of Bramocharax: B. caballeroi, B. dorioni and B. bransfordii. Samples of Roeboides bouchellei (from El Salvador), Astyanax bimaculatus (from Argentina) and Astyanax fasciatus (GenBank sequence from Brazil) were used as outgroups. We collected individuals from the type localities of most of the species (10 nominal species) described in Mesoamerica and considered valid by Lima et al . Specimens were sampled by electro-fishing and netting, individually tagged, and preserved in DMSO/EDTA buffer  or 95% ethanol. DNA voucher specimens and their associated lots were subsequently preserved in 10% buffered formalin and deposited in the Museo Nacional de Ciencias Naturales of Madrid, Spain (MNCN), and the Universidad Michoacana de San Nicolas de Hidalgo, Michoacan, Mexico (UMSNH, see Additional file 3).
DNA extraction and Sequencing
Estimated dates derived from the r8s molecular dating analyses, along with standard deviations of the node ages derived from Penalized Likelihood bootstrap analyses.
Strict Molecular Clock
Relaxed Molecular Clock NPRS with r8s Mya
95% of confidence
Merida-Perija Mountains 8–12 Mya
24.6 ML (18 – 23)
32.3 ML (23.7–42.8)
Colossoma macropomum 15 Mya
11 p (9.8 – 12.4)
14.6 p (12.8–16.3)
7.7 ML (6.36–9.32)
10.05 ML (8.35–12.23)
7.9 SD 0.165 (7.8–8.1)
5.5 p (4.9–6.4)
7.3 p (6.4–8.5)
4.6 ML (2.9–7.9)
6 ML (3.8–10)
7.1* SD 0.34 (6.0–7.8)
3.7 p (2.5–5.7)
4.9 p (3.3–7.5)
Radiation in North Mesoamerica Lineages
3.8 ML (2–6.2)
4.83 ML (2.7–8.2)
6.7* SD 0.37 (5.7–7.5)
3.1 p (1.9–4.8)
4.1 p (2.5–6.3)
Central America Radiation
3.7 ML (2–6.2)
6.9* SD 0.39 (5.9–7.7)
3.1 p (1.9–4.8)
4.2 p (0.2–5.5)
TMVB 3–6 Mya
3.6 ML (2.3–5.4)
4.7 ML (3.1–7)
5.25* SD 0.37 (4.4–6.2)
3.07 p (2.2–4.3)
4.03 p (2.8–5.7)
Chromatograms and alignments were visually checked and verified. Saturation for transition and transversion substitutions was checked by plotting the absolute number of changes at each codon position against patristic distances for coding genes only.
Phylogenetic reconstruction was performed for Bayesian Inference (BI) using MrBayes version 3.1.2 . We used Modeltest 3.07  to find the best-fit model of evolution for each gene fragment using the Bayesian Information Criterion (BIC)  (see Additional file 4). BI was performed on two data sets as follows: (1) Cytb gene only with separate best-fit models for each codon position and (2) using the three mtDNA and RAG1 (separate best-fit models for each codon position were used) genes with a separate best-fit model of evolution for each gene fragment (partition). Analyses of best-fit models of evolution and BI were performed for a subset of the data.
Bayesian analyses were performed using two independent runs of four Metropolis-coupled chains of 10 million generations each to estimate the posterior probability distribution. The first 10,000 trees were discarded as burn-in. The program Tracer v1.4  was used to assess run convergence and determine burn-in.
Sequence data were also analysed using maximum parsimony (MP) as implemented in PAUP* 4.0 b10 , NONA version 2.0  and WINCLADA version 1.00.08 . MP analyses in PAUP* and NONA/WINCLADA were done under the same heuristic search strategy. Statistical support for recovered clades was assessed using bootstrap (1000 pseudo-replications). We applied different weights for transversions and transitions according to the empiric criterion obtained in PAUP* 4.0 b10 . The two datasets run for BI were also run for MP.
Analysis 1 (for which most species and populations of the Astyanax and Bramocharax genera from Mesoamerica were represented) was used to infer phylogenetic relationships among populations. Analysis 2 (for which only a subset of the species/populations were available) was used to infer relationships among the main lineages identified in analysis 1.
Molecular clock and divergence time
Since we have a more complete data matrix, and in order to make our data comparable with previous studies in the region, we calibrated the molecular clock based on Cytb data alone rather than the combined data matrix (which included two more mitochondrial genes and one nuclear DNA gene).
Rate heterogeneity within the dataset was assessed using the likelihood-ratio test (LRT) [74, 75]. LRTs were determined by comparing the log likelihood of the optimal topology recovered by maximum likelihood analysis (using appropriate models identified by Modeltest), while enforcing the molecular clock to the log likelihood of the optimal topology recovered by one that did not. The likelihood ratio statistic is twice the difference between the two log likelihoods. This statistic is compared to a χ2 distribution with degrees of freedom equal to the number of terminals minus two following .
Because the molecular clock hypothesis was rejected, we conducted a non-parametric rate smoothing approach (NPRS) of divergence time estimation with the r8s package  to estimate divergence between the taxa. NPRS relaxes the molecular clock assumption by applying a least squares smoothing of estimates of substitutions rates.
Standard errors of divergence dates were estimated using the boot strapping procedure outlined in, and implemented by, Perl scripts in the r8s bootkit provided by Torsten Eriksson . The first 100 bootstrapped datasets were created from the original Cytb dataset with the program Mesquite v. 2.01 . Branch lengths were then re-estimated for each bootstrapped dataset in PAUP* using the original ML parameters. The resulting trees, with branch lengths, were then imported into r8s. The TN (Truncated Newton) algorithm was implemented.
We constrained the molecular clock considering two main points in the tree topology (see figure 6). First, our ingroup (node 1, Figure 6) was calibrated with the isolation of the Maracaibo basin as a result of the rise of Sierra de Perija and Merida Andes 8–12 Mya: [7, 79, 80]. Additionally for this node, we also used the oldest fossil record for a Characid (Colossoma macropomum, 15 Mya) in the Magdalena river system . This allowed us to determine the minimum and maximum values for this node (8 and 15 Mya, respectively). We also constrained the molecular clock using the final closure of the Trans-Mexican Volcanic Belt (TMVB) 3–6 Mya (node 6 in the Figure 6) [53, 60].
We thank Jerry Johnson and Petter Unmack for providing material. Thanks also to Lourdes Alcaraz for her assistance with the laboratory work, and Carlos Pedraza, Paul Bloor, Alejandro Zaldivar and Guy Reeves for helpful suggestions on an early version of the manuscript. We also thank an anonymous referee for his helpful and meticulous work in revising of this article. CPOG was supported by CONACyT – Fundación Carolina grant 196833.
- Contreras-Balderas SHO, Lozano-Vilano ML: Punta de Morro, an interesting barrier for distributional patterns of Continental fishes in North and Central Veracruz, Mexico. Publ Biol Fac Cienc Biol Univ Autóm Nuevo Leon México. 1996, 16: 37-42.Google Scholar
- Domínguez-Domínguez O, Doadrio I, Pérez-Ponce de León G: Historical biogeography of some river basins in Central Mexico evidenced by their goodeine freshwater fishes: A preliminary hypothesis using secondary Brooks Parsimony Analysis (BPA). J Biogeogr. 2006, 1437-1447. 10.1111/j.1365-2699.2006.01526.x. 33
- Huidobro L, Morrone JJ, Villalobos JL, Alvarez F: Distributional patterns of freshwater taxa (fishes, crustaceans and plants) from the Mexican Transition Zone. J Biogeogr. 2006, 33 (4): 731-741. 10.1111/j.1365-2699.2005.01400.x.View ArticleGoogle Scholar
- Morrone JJ: Biogeographical regions under track and cladistic scrutiny. J Biogeogr. 2002, 29 (2): 149-152. 10.1046/j.1365-2699.2002.00662.x.View ArticleGoogle Scholar
- Zaldivar-Riveron A, Leon-Regagnon V, de Oca ANM: Phylogeny of the Mexican coastal leopard frogs of the Rana berlandieri group based on mtDNA sequences. Mol Phylogenet Evol. 2004, 30 (1): 38-49. 10.1016/S1055-7903(03)00141-6.View ArticlePubMedGoogle Scholar
- Guzman-Speziale M, Valdes-Gonzalez C, Molina E, Gomez JM: Seismic activity along the Central America Volcanic Arc: Is it related to subduction of the Cocos plate?. Tectonophysics. 2005, 400 (1–4): 241-254. 10.1016/j.tecto.2005.03.006.View ArticleGoogle Scholar
- Iturralde-Vinent MA, MacPhee RDE: Paleogeography of the Caribbean region: Implications for cenozoic biogeography. B Am Mus Nat Hist. 1999, 1-95. 238
- Bermingham E, Martin AP: Comparative mtDNA phylogeography of neotropical freshwater fishes: testing shared history to infer the evolutionary landscape of lower Central America. Mol Ecol. 1998, 7 (4): 499-517. 10.1046/j.1365-294x.1998.00358.x.View ArticlePubMedGoogle Scholar
- Martin AP, Bermingham E: Systematics and evolution of lower Central American cichlids inferred from analysis of cytochrome b gene sequences. Mol Phylogenet Evol. 1998, 9 (2): 192-203. 10.1006/mpev.1997.0461.View ArticlePubMedGoogle Scholar
- Mateos M: Comparative phylogeography of livebearing fishes in the genera Poeciliopsis and Poecilia (Poeciliidae : Cyprinodontiformes) in central Mexico. J Biogeogr. 2005, 32 (5): 775-780. 10.1111/j.1365-2699.2005.01236.x.View ArticleGoogle Scholar
- Mateos M, Sanjur OI, Vrijenhoeck RC: Historical Biogeography of the Livebearing Fish genus Poeciliopsis (Poecilidae: Cyprinodontiformes). Evolution. 2002, 56: 972-984.View ArticlePubMedGoogle Scholar
- Perdices A, Bermingham EAM, Doadrio I: Evolutionary history of the genus Rhamdia (Teleostei: Pimelodidae) in Central America. Mol Phylogenet Evol. 2002, 25: 172-189. 10.1016/S1055-7903(02)00224-5.View ArticlePubMedGoogle Scholar
- Perdices A, Doadrio I, Bermingham E: Evolutionary history of the synbranchid eels (Teleostei: Synbranchidae) in Central America and the Caribbean islands inferred from their molecular phylogeny. Mol Phylogenet Evol. 2005, 37 (2): 460-473. 10.1016/j.ympev.2005.01.020.View ArticlePubMedGoogle Scholar
- Bussing WA: Patterns of distribution of the Central American ichthyofauna. The Great American Biotic Interchange. Edited by: Stehli FG, Webb SD. 1985, New York: Plenum Press, New York, 453-473.View ArticleGoogle Scholar
- Reeves RG, Bermingham E: Colonization, population expansion, and lineage turnover: phylogeography of Mesoamerican characiform fish. Biol J Linn Soc. 2006, 88 (2): 235-255. 10.1111/j.1095-8312.2006.00619.x.View ArticleGoogle Scholar
- Avise JC: Molecular Markers, Natural History and Evolution. New York, NY. 1994Google Scholar
- Joseph L, Moritz C, Hugall A: Molecular Support for Vicariance as a Source of Diversity in Rain-Forest. Proc R Soc Lond [Biol]. 1995, 260 (1358): 177-182. 10.1098/rspb.1995.0077.View ArticleGoogle Scholar
- Patton JL, Dasilva MNF, Malcolm JR: Gene Genealogy and Differentiation among Arboreal Spiny Rats (Rodentia, Echimyidae) of the Amazon Basin – a Test of the Riverine Barrier Hypothesis. Evolution. 1994, 48 (4): 1314-1323. 10.2307/2410388.View ArticleGoogle Scholar
- Templeton AR, Routman E, Phillips CA: Separating Population-Structure from Population History – a Cladistic-Analysis of the Geographical-Distribution of Mitochondrial-DNA Haplotypes in the Tiger Salamander, Ambystoma-Tigrinum. Genetics. 1995, 140 (2): 767-782.PubMed CentralPubMedGoogle Scholar
- Coates A, Oblando JA: The geologic evolution of the Central America Isthmus. Evolution and Environmental in Tropical America. Edited by: Jackson JBC, Budd AF, Coates AG. 1996, Chicago: Chicago University Press, Chicago, 21-56.Google Scholar
- Myers GS: Derivation of Freshwater Fish Fauna of Central America. Copeia. 1966, 766-10.2307/1441405. 4
- Concheiro Perez GA, Rican O, Orti G, Bermingham E, Doadrio I, Zardoya R: Phylogeny and biogeography of 91 species of heroine cichlids (Teleostei: Cichlidae) based on sequences of the cytochrome b gene. Mol Phylogenet Evol. 2007, 43 (1): 91-110. 10.1016/j.ympev.2006.08.012.View ArticleGoogle Scholar
- Rosen DE: Vicariance Model of Caribbean Biogeography. Syst Zool. 1975, 24 (4): 431-464. 10.2307/2412905.View ArticleGoogle Scholar
- Rosen DE: Vicariant patterns and historical explanation in biogeography. Syst Zool. 1978, 27: 159-188. 10.2307/2412970.View ArticleGoogle Scholar
- Haq BU, Hardenbol J, Vail PR: CHRONOLOGY OF FLUCTUATING SEA LEVELS SINCE THE TRIASSIC. Science. 1987, 235 (4793): 1156-1167. 10.1126/science.235.4793.1156.View ArticlePubMedGoogle Scholar
- Zeh JA, Zeh DW, Bonilla MM: Phylogeography of the harlequin beetle-riding pseudoscorpion and the rise of the Isthmus of Panama. Mol Ecol. 2003, 12 (10): 2759-2769. 10.1046/j.1365-294X.2003.01914.x.View ArticlePubMedGoogle Scholar
- Crawford AJ, Smith EN: Cenozoic biogeography and evolution in direct-developing frogs of Central America (Leptodactylidae: Eleutherodactylus) as inferred from a phylogenetic analysis of nuclear and mitochondrial genes. Mol Phylogenet Evol. 2005, 35 (3): 536-555. 10.1016/j.ympev.2005.03.006.View ArticlePubMedGoogle Scholar
- Marko PB: Fossil calibration of molecular clocks and the divergence times of geminate species pairs separated by the Isthmus of Panama. Mol Biol Evol. 2002, 19 (11): 2005-2021.View ArticlePubMedGoogle Scholar
- Baird SF, Girard CF: Descriptions of new species of fishes collected in Texas, New Mexico and Sonora, by Mr. John H. Clark, on the U. S. and Mexican Boundary Survey, and in Texas by Capt. Stewart Van Vliet, U.S.A. Proc Nat Acad Sci. 1854: 24-29.
- Calcagnotto D, Schaefer SA, DeSalle R: Relationships among characiform fishes inferred from analysis of nuclear and mitochondrial gene sequences. Mol Phylogenet Evol. 2005, 36 (1): 135-153. 10.1016/j.ympev.2005.01.004.View ArticlePubMedGoogle Scholar
- Ortí G, Meyer A: The radiation of characiform fishes and the limits of resolution of mitochondrial ribosomal DNA sequences. Syst Biol. 1997, 46 (1): 75-100. 10.2307/2413637.View ArticlePubMedGoogle Scholar
- Otero O, Gayet M: Palaeoichthyofaunas from the Lower Oligocene and Miocene of the Arabian Plate: palaeoecological and palaeobiogeographical implications. Palaeogeography Palaeoclimatology Palaeoecology. 2001, 165 (1–2): 141-169. 10.1016/S0031-0182(00)00158-9.View ArticleGoogle Scholar
- Gayet M, Marshall LG, Sempere T, Meunier FJ, Cappetta H, Rage JC: Middle Maastrichtian vertebrates (fishes, amphibians, dinosaurs and other reptiles, mammals) from Pajcha Pata (Bolivia). Biostratigraphic, palaeoecologic and palaeobiogeographic implications. Palaeogeography Palaeoclimatology Palaeoecology. 2001, 169 (1–2): 39-68. 10.1016/S0031-0182(01)00214-0.View ArticleGoogle Scholar
- Eschmeyer WN: Catalog of fishes. Updated database version of June 2007. FishBase. 2007Google Scholar
- Lima FCT, Malabarba LR, Buckup PA, Pezzi Da Silva JF, Vari RP, Harold A, Benine R, Oyakawa OT, Pavanelli CS, Menezes NA, et al: Genera Incertae Sedis in Characidae. Checklist of the Freshwater Fishes of South and Central America. Edited by: Reis RE, Kullander SO, Ferraris CJ Jr. 2003, Porto Alegre Brasil: EDIPUCRS, 106-168.Google Scholar
- Strecker U, Bernatchez L, Wilkens H: Genetic divergence between cave and surface populations of Astyanax in Mexico (Characidae, Teleostei). Mol Ecol. 2003, 12 (3): 699-710. 10.1046/j.1365-294X.2003.01753.x.View ArticlePubMedGoogle Scholar
- Strecker U, Faundez VH, Wilkens H: Phylogeography of surface and cave Astyanax (Teleostei) from Central and North America based on cytochrome b sequence data. Mol Phylogenet Evol. 2004, 33 (2): 469-481. 10.1016/j.ympev.2004.07.001.View ArticlePubMedGoogle Scholar
- Dowling TE, Martasian DP, Jeffery WR: Evidence for Multiple Genetic Forms with Similar Eyeless Phenotypes in the Blind Cavefish, Astyanax mexicanus. Mol Biol Evol. 2002, 19 (4): 446-455.View ArticlePubMedGoogle Scholar
- Jeffery WR: Cave fish as a model system in evolutionary developmental biology. Developmental Biology. 2001, 231: 1-12. 10.1006/dbio.2000.0121.View ArticlePubMedGoogle Scholar
- Lozano-Vilano ML, Contreras-Balderas S: Astyanax armandoi, n. sp. from Chiapas, Mexico (Pisces, Ostariophysi: Characidae) with a Comparison to the nominal species A. aeneus and A. mexicanus. Universidad y Ciencia. 1990, 7: 95-107.Google Scholar
- Paulo-Maya J: Análisis morfométrico del género Astyanax (Pisces: Characidae) en México. 1994, México, D. F. : Instituto Politécnico NacionalGoogle Scholar
- Bernatchez L, Chouinard A, Lu GQ: Integrating molecular genetics and ecology in studies of adaptive radiation: whitefish, Coregonus sp., as a case study. Biol J Linn Soc. 1999, 68 (1–2): 173-194. 10.1111/j.1095-8312.1999.tb01165.x.View ArticleGoogle Scholar
- Brunner PC, Douglas MR, Osinov A, Wilson CC, Bernatchez L: Holarctic phylogeography of Arctic charr (Salvelinus alpinus L.) inferred from mitochondrial DNA sequences. Evolution. 2001, 55 (3): 573-586. 10.1554/0014-3820(2001)055[0573:HPOACS]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Danley PD, Kocher TD: Speciation in rapidly diverging systems: lessons from Lake Malawi. Mol Ecol. 2001, 10 (5): 1075-1086. 10.1046/j.1365-294X.2001.01283.x.View ArticlePubMedGoogle Scholar
- Valdez-Moreno ME: A checklist of the freshwater ichthyofauna from El Peten and Alta Verapaz, Guatemala, with notes for its conservation and management. Zootaxa. 2005, 43-60. 1072
- Rosen DE: A New Tetragonopterine Characid Fish From Guatemala. Am Mus Novit. 1970, 1-17. 2435
- Rosen DE: Origin of the Characid Fish Genus Bramocharax and Description of a Second, More Primitive, Species in Guatemala. Am Mus Novit. 1972, 1-21. 2500
- Schmitter-Soto JJ, Valdez-Moreno ME, Rodiles-Hernandez R, González-Díaz AA: Astyanax armandoi, a Junior Synonym of Astyanax aeneus (Teleostei: Characidae). Copeia. 2008, 409-413. 10.1643/CI-07-012. 2
- Contreras-Balderas S, Lozano-Vilano ML: Problemas nomenclaturales de las formas mexicanas del género Astyanax (Pisces: Characidae). Zoología Informa. 1988, 38: 1-13.Google Scholar
- Protas ME, Hersey C, Kochanek D, Zhou Y, Wilkens H, Jeffery WR, Zon LI, Borowsky R, Tabin CJ: Genetic analysis of cavefish reveals molecular convergence in the evolution of albinism. Nat Genet. 2006, 38 (1): 107-111. 10.1038/ng1700.View ArticlePubMedGoogle Scholar
- Wilkens H, Strecker U: Convergent evolution of the cavefish Astyanax (Characidae, Teleostei): genetic evidence from reduced eye-size and pigmentation. Biol J Linn Soc. 2003, 80 (4): 545-554. 10.1111/j.1095-8312.2003.00230.x.View ArticleGoogle Scholar
- Miller R: Freshwater Fishes of México. 2005, Chicago: The University of Chicago Press, 1:Google Scholar
- Ferrari L, López-Martínez M, Aguirre-Díaz G, Carrasco-Núñez G: Space Time patterns of Cenozoic arc volcanism in central México: from the Sierra Madre Occidental to the Mexican Volcanic Belt. Geology. 1999, 27: 303-306. 10.1130/0091-7613(1999)027<0303:STPOCA>2.3.CO;2.View ArticleGoogle Scholar
- Doadrio I, Perdices A: Phylogenetic relationships among the Ibero-African cobitids (Cobitis, cobitidae) based on cytochrome b sequence data. Mol Phylogenet Evol. 2005, 37 (2): 484-493. 10.1016/j.ympev.2005.07.009.View ArticlePubMedGoogle Scholar
- Murphy WJ, Collier GE: Phylogenetic relationships within the aplocheiloid fish genus Rivulus (Cyprinodontiformes, Rivulidae): Implications for Caribbean and Central American biogeography. Mol Biol Evol. 1996, 13 (5): 642-649.View ArticlePubMedGoogle Scholar
- Bartoli G, Sarnthein M, Weinelt M, Erlenkeuser H, Garbe-Schonberg D, Lea DW: Final closure of Panama and the onset of northern hemisphere glaciation. Earth Planet Sci Lett. 2005, 237 (1–2): 33-44. 10.1016/j.epsl.2005.06.020.View ArticleGoogle Scholar
- Hrbek T, Seckinger J, Meyer A: A phylogenetic and biogeographic perspective on the evolution of poeciliid fishes. Mol Phylogenet Evol. 2007, 43 (3): 986-998. 10.1016/j.ympev.2006.06.009.View ArticlePubMedGoogle Scholar
- Collins L, Coates A, Berggren W, Aubry M, Zhang J: The late Miocene Panama isthmian strait. Geology. 1996, 24 (8): 687-690. 10.1130/0091-7613(1996)024<0687:TLMPIS>2.3.CO;2.View ArticleGoogle Scholar
- Ferrari L, Conticelli S, Potrone CM, Manetti P: Late Miocene volcanism and intra-arc tectonics during the early development of the Trans-Mexican Volcanic Belt. Tectonophysics. 2000, 318: 161-185. 10.1016/S0040-1951(99)00310-8.View ArticleGoogle Scholar
- Ferrari L, Tagami T, Eguchi M, Orozco-Esquivel MT, Petrone CM, Jacobo-Albarran J, Lopez-Martinez M: Geology, geochronology and tectonic setting of late Cenozoic volcanism along the southwestern Gulf of Mexico: The Eastern Alkaline Province revisited. J Volcanol Geotherm Res. 2005, 146 (4): 284-306. 10.1016/j.jvolgeores.2005.02.004.View ArticleGoogle Scholar
- Mulcahy DG, Mendelson JR: Phylogeography and Speciation of the Morphologically Variable, Widespread Species Bufo valliceps, Based on Molecular Evidence from mtDNA. Mol Phylogenet Evol. 2000, 17 (2): 173-189. 10.1006/mpev.2000.0827.View ArticlePubMedGoogle Scholar
- Mulcahy DG, Morrill BH, Mendelson JR: Historical biogeography of lowland species of toads (Bufo) across the Trans-Mexican Neovolcanic Belt and the Isthmus of Tehuantepec. J Biogeogr. 2006, 33 (11): 1889-1904. 10.1111/j.1365-2699.2006.01546.x.View ArticleGoogle Scholar
- Smith SA, Bermingham E: The biogeography of lower Mesoamerican freshwater fishes. J Biogeogr. 2005, 32 (10): 1835-1854. 10.1111/j.1365-2699.2005.01317.x.View ArticleGoogle Scholar
- Guzman-Speziale M: Active seismic deformation in the grabens of northern Central America and its relationship to the relative motion of the North America-Caribbean plate boundary. Tectonophysics. 2001, 337 (1–2): 39-51. 10.1016/S0040-1951(01)00110-X.View ArticleGoogle Scholar
- Seutin G, White BN, Boag PT: Preservation of Avian Blood and Tissue Samples for DNA Analyses. Can J Zool. 1991, 69 (1): 82-90. 10.1139/z91-013.View ArticleGoogle Scholar
- Sambrook J, Fritsch E, Maniatis T: Molecular cloning: A laboratory manual. 1989, New York: Cold Spring LaboratoryGoogle Scholar
- Quenouille B, Bermingham E, Planes S: Molecular systematics of the damselfishes (Teleostei : Pomacentridae): Bayesian phylogenetic analyses of mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol. 2004, 31 (1): 66-88. 10.1016/S1055-7903(03)00278-1.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MrBayes: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14 (9): 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Rambaut A, Drummond A: Tracer [computer program]. 2007, [http://tree.bio.ed.ac.uk/software/tracer]1.4Google Scholar
- Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods). 1998, Sunderland, Massachussetts: Sinauer AssociatesGoogle Scholar
- Goloboff PA: NONA. Noname (a bastard son of Pee-Wee). 2.0 (32 bit version) edn. New York Program and documentation. 1993, Computer program distributed by J.M. Carpenter, Department of Entomology, American Museum of Natural History New York, 1993Google Scholar
- Nixon KC: Winclada (BETA). 1999, New York: Published by the author, ITHACA, NY, 0.9.9Google Scholar
- Goldman N: Simple Diagnostic Statistical Tests of Models for DNA Substitution. Journal of Molecular Evolution. 1993, 37 (6): 650-661.PubMedGoogle Scholar
- Sanderson MJ: A nonparametric approach to estimating divergence times in the absence of rate constancy. Mol Biol Evol. 1997, 14 (12): 1218-1231.View ArticleGoogle Scholar
- Huelsenbeck JP, Rannala B: Phylogenetic methods come of age: Testing hypotheses in an evolutionary context. Science. 1997, 276 (5310): 227-232. 10.1126/science.276.5310.227.View ArticlePubMedGoogle Scholar
- r8s bootkit 2. [http://www.bergianska.se/index_forskning_soft.html]
- Mesquite: a modular system for evolutionary analysis V. 2.01.
- Díaz de Gamero ML: The changing course of the Orinoco River during the Neogene: A review. Palaeogeography Palaeoclimatology Palaeoecology. 1996, 123 (1–4): 385-402. 10.1016/0031-0182(96)00115-0.View ArticleGoogle Scholar
- Lundberg JG: The temporal context for the diversification of Neotropical Fishes. Phylogeny and Classification of Neotropical Fishes. Edited by: Malabarba LR, Reis RE, Vari RP, Lucena ZMS, Lucena CAS. 1998, Porto Alegre, Brasil EDIPUCRS, 49-68.Google Scholar
- Lundberg JG, Machadoallison A, Kay RF: Miocene Characid Fishes from Colombia – Evolutionary Stasis and Extirpation. Science. 1986, 234 (4773): 208-209. 10.1126/science.234.4773.208.View ArticlePubMedGoogle Scholar
- Ferrari L, Lopez-Martinez M, Aguirre-Diaz G, Carrasco-Nunez G: Space-time patterns of Cenozoic arc volcanism in central Mexico: From the Sierra Madre Occidental to the Mexican Volcanic Belt. Geology. 1999, 27 (4): 303-306. 10.1130/0091-7613(1999)027<0303:STPOCA>2.3.CO;2.View ArticleGoogle Scholar
- Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN: DNA barcoding Australia's fish species. Philos Trans R Soc Lond, B. 2005, 360 (1462): 1847-1857. 10.1098/rstb.2005.1716.View ArticleGoogle Scholar
- Zardoya R, Doadrio I: Phylogenetic relationships of Iberian cyprinids: systematic and biogeographical implications. Proc R Soc Lond [Biol]. 1998, 265 (1403): 1365-1372. 10.1098/rspb.1998.0443.View ArticleGoogle Scholar
- Palumbi S, Martin AP, Romano S, McMillan WO, Stice L, Grabowski G: The Simple Fool's Guide to PCR. 1991, Special publication, Honolulu: University Hawaii Press, 2.0Google Scholar
- Murphy WJ, Thomerson JE, Collier GE: Phylogeny of the neotropical killifish family Rivulidae (Cyprinodontiformes, Aplocheiloidei) inferred from mitochondrial DNA sequences. Mol Phylogenet Evol. 1999, 13 (2): 289-301. 10.1006/mpev.1999.0656.View ArticlePubMedGoogle Scholar