Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Molecular Phylogeny and Biogeography of the Amphidromous Fish Genus Dormitator Gill 1861 (Teleostei: Eleotridae)

  • Sesángari Galván-Quesada ,

    galygq@hotmail.com (SGQ); goodeido@yahoo.com.mx (ODD)

    Affiliations Programa Institucional de Doctorado en Ciencias Biológicas, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán, México, Laboratorio de Biología Acuática, Facultad de Biología, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán, México

  • Ignacio Doadrio,

    Affiliation Departamento de Biodiversidad y Biología Evolutiva, Museo Nacional de Ciencias Naturales, CSIC, Madrid, Spain

  • Fernando Alda,

    Current address: Museum of Natural Science, Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana, United States of America

    Affiliation Smithsonian Tropical Research Institute, Apartado 2072, Balboa, Republic of Panama

  • Anabel Perdices,

    Affiliations Departamento de Biodiversidad y Biología Evolutiva, Museo Nacional de Ciencias Naturales, CSIC, Madrid, Spain, Smithsonian Tropical Research Institute, Apartado 2072, Balboa, Republic of Panama

  • Ruth Gisela Reina,

    Affiliation Smithsonian Tropical Research Institute, Apartado 2072, Balboa, Republic of Panama

  • Martín García Varela,

    Affiliation Instituto de Biología, Universidad Nacional Autónoma de México, Distrito Federal, México

  • Natividad Hernández,

    Affiliation Instituto de Medicina Tropical Pedro Kourí, Apartado 601, La Habana, Cuba

  • Antonio Campos Mendoza,

    Affiliation Laboratorio de Biología Acuática, Facultad de Biología, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán, México

  • Eldredge Bermingham,

    Current address: Patricia and Phillip Frost Museum of Science, Miami, Florida, United States of America

    Affiliation Smithsonian Tropical Research Institute, Apartado 2072, Balboa, Republic of Panama

  • Omar Domínguez-Domínguez

    galygq@hotmail.com (SGQ); goodeido@yahoo.com.mx (ODD)

    Affiliation Laboratorio de Biología Acuática, Facultad de Biología, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Michoacán, México

Abstract

Species of the genus Dormitator, also known as sleepers, are representatives of the amphidromous freshwater fish fauna that inhabit the tropical and subtropical coastal environments of the Americas and Western Africa. Because of the distribution of this genus, it could be hypothesized that the evolutionary patterns in this genus, including a pair of geminate species across the Central American Isthmus, could be explained by vicariance following the break-up of Gondwana. However, the evolutionary history of this group has not been evaluated. We constructed a time-scaled molecular phylogeny of Dormitator using mitochondrial (Cytochrome b) and nuclear (Rhodopsin and β-actin) DNA sequence data to infer and date the cladogenetic events that drove the diversification of the genus and to relate them to the biogeographical history of Central America. Two divergent lineages of Dormitator were recovered: one that included all of the Pacific samples and another that included all of the eastern and western Atlantic samples. In contrast to the Pacific lineage, which showed no phylogeographic structure, the Atlantic lineage was geographically structured into four clades: Cameroon, Gulf of Mexico, West Cuba and Caribbean, showing evidence of potential cryptic species. The separation of the Pacific and Atlantic lineages was estimated to have occurred ~1 million years ago (Mya), whereas the four Atlantic clades showed mean times of divergence between 0.2 and 0.4 Mya. The splitting times of Dormitator between ocean basins are similar to those estimated for other geminate species pairs with shoreline estuarine preferences, which may indicate that the common evolutionary histories of the different clades are the result of isolation events associated with the closure of the Central American Isthmus and the subsequent climatic and oceanographic changes.

Introduction

Fish display a wide array of life-history strategies, many of which affect their dispersal potential and their subsequent geographic differentiation and intraspecific variability [13]. The implications and differences in ecology, evolutionary history and biogeography associated with different life-history strategies have been the focus of scientific debate from the time of Darwin [4] until the present [58]. For example, diadromy is a life-history behavior in which individuals spend predictable phases of their life cycle in freshwater rivers or in the ocean, typically for feeding or reproduction [9]. Because of their intermediate relationships with the continental and marine realms, fishes with this life-history strategy are usually characterized by high dispersal potential, widespread distributions and shallow genetic differentiation [5,10].

Brackish environments in the Neotropics host a substantial amount of diadromous fish fauna, including gobies, mullets, snooks and sleepers. However, few studies have been conducted on these groups of fish, which represent an important component of the biodiversity of these ecosystems [11]. Investigations of the evolutionary history of diadromous taxa are even rarer despite the value of studying their genetic variation for understanding the evolution of coastal environments and its influence on species geographic distribution [12]. For example, a recent study of Agonostomus monticola highlighted the importance of major vicariant events in generating cryptic diversity in diadromous species [10].

In the Neotropics, the closing of the Central American Seaway represented a dramatic event separating marine and coastal organisms and facilitating the emergence of geminate species pairs on each side of the Isthmus that followed independent evolutionary trajectories [13]. However, comparisons of divergence times among geminate species pairs do not support a single and simultaneous divergence time for all taxa, suggesting that species might have responded differently to the complex geological evolution of the Isthmus and its new habitat development [1416]. For instance, geminate pairs of gastropods inhabiting high intertidal mangroves show shallower genetic divergences compared to those inhabiting lower intertidal and subtidal environments [17]. Hence, ecological differences, such as habitat depth, may influence the timing of separation between geminate clades [17,18].

Fish species of the genus Dormitator Gill 1861 (Teleostei: Eleotridae) are amphidromous, which is a form of diadromy in which adults live and reproduce in freshwaters. After the eggs hatch, the larvae of this kind of species drift to the sea, where they spend a variable amount of time, potentially dispersing long distances, before returning to freshwater streams [19]. Their common name “sleepers” arises from the apparent lack of motility of these species as adults [20]. Dormitator inhabit freshwater and brackish environments along the tropical and subtropical coasts and estuaries of the eastern Pacific and Atlantic Oceans [21,22]. On the Pacific coast, the genus is distributed from the Gulf of California to Peru [22], including the Galapagos Islands [23]. On the western Atlantic coast, the distribution area ranges from North Carolina to Brazil [22], including the Antilles. On the eastern Atlantic coast, the genus ranges from Senegal to Angola [24]. Four species are recognized in the Atlantic [2427]: Dormitator maculatus (Bloch 1792), which is found in the western Atlantic from southern USA to Central America and southeastern Brazil; Dormitator cubanus Ginsburg 1953 in western Cuba (Pinar del Río); Dormitator lophocephalus Hoedeman 1951 in Suriname; and Dormitator lebretonis (Steindachner 1870) occurring from Senegal to Namibia. In contrast, only one species is found on the eastern Pacific coast: Dormitator latifrons (Richardson 1844). Within this region, Dormitator latifrons mexicanus Ginsburg 1953, restricted to the Pacific coast of Mexico, has been described as a distinct subspecies from the nominal Dormitator latifrons latifrons.

Probably because of its wide distribution range, Dormitator species have suffered from considerable taxonomic instability, including the use of at least 20 synonyms [2427]. However, Dormitator have been the subject of few systematic studies [2830], and their genetic diversity and phylogenetic relationships are unknown. Furthermore, because of the species distribution on both sides of the Atlantic and their allopatric distribution on the Pacific and the Atlantic slopes of the Americas, it could be hypothesized that the splitting of continentally disjunct lineages of Dormitator coincides with the break-up of Gondwana in the Early Cretaceous and that this genus may include at least one pair of geminate species: D. latifrons and D. maculatus [13,31]. The genus Dormitator is therefore a good model to 1) study the effects of geological events and life-history strategies (amphidromy) on species divergence and distribution patterns and 2) uncover hidden diversity of putative cryptic and geminate species pairs by using molecular methods. To address these objectives, we inferred a molecular (mitochondrial and nuclear) phylogenetic hypothesis of Dormitator species and used molecular clock analyses to investigate the patterns and timescale of lineage divergence across its distribution in America and Africa.

Materials and Methods

Ethics statement

Field collections did not involve endangered or protected species. Field and laboratory protocols used in this study, including sampling procedures, were reviewed and approved by the Mexican Ministry of Environmental and Natural Resources (SEMARNAT), under collection permit number FAUT 0202. Further approval by an ethics committee was not necessary because this research did not include animal experimentation. Samples requested from other institutions, such as the Smithsonian Tropical Research Institute Neotropical Fish Collection (STRI, Panama) and Museo Nacional de Ciencias Naturales (CSIC, Spain), were also used. These samples were collected following sampling procedures reviewed and approved by STRI’s “Institutional Animal Care and Use Committee” and by CSIC’s “Ethics Committee”, respectively.

Specimen collection

Specimens were captured with seine, gill nets or electrofishing when possible. Immediately after capture, all individuals were anesthetized using Tricaine methanesulfonate (MS-222) to alleviate suffering. Tissue samples (~3 mm2 fin clips) were obtained, preserved in 95% ethanol or DMSO buffer and stored at 4°C. Once tissue samples were collected, fish were released at the same collecting site after confirming recovery of total motility (n = 100), or fish were humanely euthanized with an overdose of MS-222 (n = 158). Death was confirmed after no gill movement was observed for at least 10 minutes. Specimens were then individually tagged, fixed in formalin and then transferred to 70% ethanol for long-term storage. Voucher specimens are deposited in the Colección Ictiológica de la Facultad de Biología (CPUM-Universidad Michoacana de San Nicolás de Hidalgo, Mexico), Museo Nacional de Ciencias Naturales (CSIC Spain), Smithsonian Tropical Research Institute Neotropical Fish Collection (Panama) and Colección Nacional de Peces, Instituto de Biología (Universidad Nacional Autónoma de México, Mexico).

Dormitator specimens were collected from coastal environments along the Pacific and Atlantic slopes of the Americas and the eastern Atlantic coast of Africa. According to the described species distribution, our sampling, comprising 84 collecting locations, included specimens of D. latifrons (both putative subspecies), D. maculatus, D. cubanus and D. lebretonis (Fig 1, Table 1). Additionally, samples of Gobiomorus dormitor and Eleotris senegalensis were included in the analyses as outgroups because of their close evolutionary relationship with Dormitator [2830].

thumbnail
Fig 1. Distribution of Dormitator species and sampling locations.

Shaded areas represent the distribution range of Dormitator. Blue shading indicates the Pacific distribution (D. latifrons), pink shading indicates the western Atlantic distribution (including D. maculatus, D. cubanus and D. lophocephalus), and purple shading indicates the eastern Atlantic distribution areas (D. lebretonis). Dots represent sampling locations, and the colors represent the monophyletic groups recovered by our phylogenetic hypothesis as detailed below. Dot color code: blue = Pacific; red = Gulf of Mexico; purple = Cameroon; green = Caribbean; and yellow = West Cuba.

https://doi.org/10.1371/journal.pone.0153538.g001

thumbnail
Table 1. Summary of Dormitator samples and locations analyzed in this study.

https://doi.org/10.1371/journal.pone.0153538.t001

DNA sequencing

Whole genomic DNA was extracted from all of the tissue samples using a standard proteinase K and phenol-chloroform protocol [32] or Qiagen DNeasy Tissue Kits (Qiagen, Inc., Valencia, CA, USA). The mitochondrial Cytochrome b gene (Cytb) was amplified by polymerase chain reaction (PCR) using the primers GluDG [33] and H16460 [34]. Additionally, a subset of samples representative of the genetic variation found in the Cytb gene tree were sequenced using the primers RhF193 and RhR1039 [35] and the primers BactFor and BactRev [36] for the nuclear Rhodopsin (Rh) and Beta actin (β-actin) genes, respectively. PCRs were carried out in 25 μL volume reactions containing the following: 10X reaction buffer, 0.5 μM each primer, 0.2 mM dNTP, 2 mM MgCl2 and 1U of Taq DNA polymerase (Invitrogen). Thermocycling conditions consisted in an initial denaturation step at 94°C (2 min) followed by different time and temperature cycles depending on the gene: 35 cycles of denaturation at 94°C (45 s), annealing at 46°C (1 min) and extension at 72°C (90 s) for Cytb; 5 cycles of denaturation at 94°C (30 s), annealing at 50°C (45 s) and extension at 72°C (45 s), followed by 35 cycles of denaturation at 94°C (30 s), annealing at 54°C (45 s) and extension at 72°C (45 s) for Rh; 35 cycles of denaturation at 94°C (30 s), annealing at 55°C (40 s), extension at 72°C (90 s) for β-actin and a final extension at 72°C (5 min) in all cases. All gene fragments were sequenced in both directions using the same PCR primers. Sequencing was performed by MACROGEN Inc. (Korea) sequencing service, High-Throughput Genomics Unit sequencing service (USA) and the Smithsonian Tropical Research Institute sequencing facility (Panama). Chromatograms were visually examined and then edited and assembled using Bioedit 7.2.5 [37]. DNA sequences are available in the GenBank database under the following accession numbers: KU764787-KU765046 for Cytb, KU765047-KU765129 for Rh and KU958382-KU958464 for β-actin.

Phylogenetic inference

The evolutionary substitution models that best fit our data were determined for each gene using jModeltest2 [38] and the Akaike information criterion (AIC, Table 2). Once best-fit models were determined, they were used in all of the subsequent analyses. Phylogenetic hypotheses were independently inferred for each molecular marker, nuclear data set (Rh + β-actin) and the complete concatenated data set (Cytb + Rh + β-actin). Concatenated data sets only included individuals for which all genes were successfully sequenced. Maximum likelihood (ML) trees were generated using RAxMLBlackBox [39]. Genes were considered as individual partitions and treated independently with respect to evolutionary models and the optimization of branch lengths. Node support was assessed using 100 bootstrap ML replicates. Bayesian inference (BI) analyses were performed using MrBayes v.3.2.2 [40] via the CIPRES portal [41]. Each gene was considered as a distinct partition with unlinked maximum likelihood models. Two simultaneous Markov chain Monte Carlo (MCMC) searches were completed with four chains for 1 x 107 generations, and trees were sampled every 1000 generations with the first 25% of the trees discarded as burn-in. Convergence between runs was assessed by monitoring the standard deviation of split frequencies with MrBayes v.3.2.2 and by using the effective sampling size (ESS) criterion in Tracer v.1.6 [42].

thumbnail
Table 2. Gene information and selected evolutionary models.

https://doi.org/10.1371/journal.pone.0153538.t002

Mean uncorrected genetic p-distances (Dp) and their associated standard errors (S.E., 1000 bootstrap replicates) were calculated for the three genes independently and for the complete concatenated data set (Cytb + Rh + β-actin) within and between clades using MEGA v.6 [43].

Species tree and estimates of divergence time

The time to the most recent common ancestor (TMRCA) and confidence intervals (95% highest posterior density: HPD) were estimated for each clade in the Cytb data set using a relaxed molecular clock with an uncorrelated lognormal distribution of rates in BEAST v.1.8.0 [44]. A Yule speciation model was assumed as a tree prior (i.e., a constant rate of speciation per lineage [45]). Two independent analyses were performed by running the MCMC for 5 x 107 generations, with trees sampled every 5000 generations. The runs were examined for convergence and adequate ESS using Tracer v.1.6 and combined using LogCombiner v.1.8.0 with a burn-in fraction of 10%. The final consensus tree was produced using TreeAnnotator v.1.8.0.

To calibrate the molecular clock, we used previously published sequence and fossil record data of Gobiiformes (S1 Table). Because the Yule speciation model assumes that each tip of the tree represents one species, we selected one sequence from each of the main lineages of Dormitator here generated and included 37 sequences from 32 genera of Gobiiformes and two genera of Kurtiformes (S2 Table) to infer a molecular phylogeny anchored by four calibration points derived from six fossil species of Gobiiformes, with mean ages ranging from 12.5 to 52 million years ago (Mya) (S1 Table). To account for uncertainties in fossil dates or conflicts between fossils and molecules [46], fossil data points were included as soft calibration points using lognormal prior distributions in the stem nodes of interest.

Also, we estimated divergence times among Dormitator clades using the multispecies coalescent method *BEAST implemented in BEAST v.1.8.0. This method estimates a species tree while taking into account variation among gene trees [47]. We ran *BEAST using all sequences from all individuals and assigning them to five species, corresponding to the monophyletic clades obtained in the previous gene-tree phylogenetic hypothesis.

We used a relaxed lognormal molecular clock and a Birth-Death speciation model for each gene tree. We calibrated the molecular clock, incorporating as normal prior the substitution rate of 5.49 x 10−2 substitutions/site/million years (S.D. = 2.94 x 10−2) estimated for the Cytb gene in the fossil calibrated phylogeny of Gobiiformes described above, and estimated the substitution rate of Rh and β-actin genes relative to Cytb. We performed two independent MCMC runs, each for 5 x 107 generations, sampling every 5000 iterations. Each run was checked for convergence and adequate ESS sampling in Tracer v.1.8, combined using LogCombiner v.1.8.0 and summarized with TreeAnnotator v.1.8.0. All BEAST v.1.8.0 analyses were run in the CIPRES portal.

Results

Mitochondrial and nuclear gene trees

The mitochondrial Cytb gene was sequenced for 258 individuals of Dormitator (Table 1), which showed 306 variable positions among the 1041 base pairs (bp) sequenced (see Table 2 for additional polymorphism information and selected evolutionary models). The ML and BI Cytb gene trees produced identical topologies (Fig 2) that showed two largely divergent (Dp = 8.8%, S.E. = 0.8%) and highly supported lineages (bootstrap support, Bs = 100, and Bayesian posterior probability, Pp = 100). One lineage included all specimens identified as D. latifrons from the Pacific coast, and the other lineage included all specimens from the Atlantic coast, including D. maculatus, D. cubanus and D. lebretonis. The Pacific lineage showed a complete lack of phylogeographic structure from North Mexico to Ecuador. Conversely, four geographically structured clades were recovered in the Atlantic lineage with a within-group mean sequence divergence of Dp = 3% (S.E. = 0.3%). One clade in the Gulf of Mexico was formed by D. maculatus and represented the sister group to the other three clades (Dp = 4.7–6.7%, Table 3). The West Cuba clade included individuals from Isla de la Juventud and the province of Pinar del Río on the island of Cuba, the type locality of D. cubanus. The West Cuba clade was sister to the Caribbean clade (Dp = 2.8%, S.E. = 0.5, Table 3), which included all of the remaining Western Atlantic and Caribbean D. maculatus samples from the Dominican Republic, Cuba, Honduras, Nicaragua, Panama and Venezuela. The sample of D. lebretonis from Cameroon was recovered as the sister of the West Cuba and Caribbean clades; however, the support for this finding was low.

thumbnail
Fig 2. Phylogenetic hypothesis of Dormitator based on mitochondrial Cytb gene sequences.

The bullets on the nodes indicate posterior probabilities (Pp) rendered by the BI analysis (upper half) and bootstrap support (Bs) for the ML analysis (lower half). Black: high support (Pp ≥ 0.99, Bs ≥ 90%), gray: good support (Pp ≥ 0.90, Bs ≥ 65%), and white: low (Pp < 0.90, Bs < 65%). Numbers at the left of the nodes indicate the estimated mean TMRCA in Mya for each node, and between parenthesis the 95% HPD confidence intervals for each estimated date. The colors of the clades correspond to the geographic origin of the sample locations (colored dots in Fig 1).

https://doi.org/10.1371/journal.pone.0153538.g002

thumbnail
Table 3. Mean p-distances between and within Dormitator clades obtained in the present phylogenetic analysis.

https://doi.org/10.1371/journal.pone.0153538.t003

The Rh gene (831 bp) was sequenced in 81 individuals and included 38 variable sites, whereas the β-actin gene (972 bp) was sequenced in 81 individuals and included 52 variable sites (Tables 1 and 2). The selected evolutionary models for the nuclear markers are shown in Table 2. The ML and BI reconstructions revealed the same topology for each nuclear gene separately (data not shown). Here, we present the phylogenetic hypothesis generated by the Rh + β-actin concatenated data set (Fig 3).

thumbnail
Fig 3. Phylogenetic hypothesis of Dormitator based on nuclear genes (Rh + β-actin) sequences.

The bullets on the nodes indicate posterior probabilities (Pp) rendered by the BI analysis (upper half) and bootstrap support (Bs) for the ML analysis (lower half). Black: high support (Pp ≥ 0.99, Bs ≥ 90%) and gray: good support (Pp ≥ 0.90, Bs ≥ 65%). Colors of the clades: blue = Pacific, yellow = West Cuba, and brown = Gulf of Mexico + Cameroon + Caribbean.

https://doi.org/10.1371/journal.pone.0153538.g003

The nuclear (Rh + β-actin) phylogenetic hypothesis was consistent with the Cytb gene tree, revealing the two Pacific and Atlantic major lineages with high support (Pp ≥ 0.99, Bs ≥ 90%) and with genetic divergences ranging from Dp = 0.8% (S.E. = 0.3%) to Dp = 1.6% (S.E. = 0.4%) for Rh and β-actin, respectively. However, within the Atlantic lineage, only two clades were recovered by both nuclear genes: the West Cuba clade and the remaining Atlantic Dormitator samples, including D. maculatus from the Gulf of Mexico and the Caribbean and D. lebretonis from Cameroon. Average Dp values between these Atlantic clades ranged from 0.5% (S.E. = 0.2%) to 0.4% (S.E. = 0.2%) for Rh and β-actin, respectively.

The complete concatenated data set (Cytb + Rh + β-actin) analysis (Fig 4) produced a topology that was congruent with that of the previously described Cytb gene tree (Fig 2), with the highest support for the African D. lebretonis sister relationship occurring with the West Cuba and Caribbean clades. The Pacific and Atlantic lineages showed an overall Dp of 4.0% (S.E. = 0.4%).

thumbnail
Fig 4. Phylogenetic hypothesis of Dormitator based on the concatenated mitochondrial and nuclear (Cytb + Rh + β-actin) genes sequences.

The bullets on the nodes indicate posterior probabilities (Pp) rendered by the BI analysis (upper half) and bootstrap support (Bs) for the ML analysis (lower half). Black: high support (Pp ≥ 0.99, Bs ≥ 90%) and gray: good support (Pp ≥ 0.90, Bs ≥ 65%). The colors of the clades correspond to the geographic origin of the sample locations (colored dots in Fig 1).

https://doi.org/10.1371/journal.pone.0153538.g004

Species tree and molecular dating

The molecular clock analysis results indicated a mean TMRCA between the Pacific and Atlantic lineages of 1.5 Mya (95% HPD: 0.3–3.1 Mya), which corresponded to the Mid- to Late Pliocene period (1.75–5.3 Mya) (Fig 2, S1 Fig). The divergence times among the Atlantic lineage clades were all recent and showed overlapping confidence intervals, indicating that all of these cladogenetic events occurred during the Pleistocene (95% HPD ranging between 0.1–1.7 Mya). The TMRCA ranged from 0.8 Mya (95% HPD: 0.2–1.7 Mya) for the split between the Gulf of Mexico and remaining Atlantic clades (Cameroon, West Cuba and Caribbean) to 0.4 Mya (95% HPD: 0.1–0.9 Mya) for the West Cuba and Caribbean clades. The mean divergence of the Africa and Western Atlantic clades was estimated at 0.7 Mya (95% HPD: 0.1–1.2 Mya) (Fig 2, S1 Fig).

The species-tree hypothesis agreed with the mitochondrial and concatenated gene trees in recovering the Pacific and Atlantic Dormitator as reciprocally monophyletic groups but differed regarding the relationships among Atlantic Dormitator taxa (Fig 5). Dormitator lebretonis from Cameroon was recovered in a basal position from the remaining Atlantic species (Pp = 0.75). Also, D. maculatus from the Gulf of Mexico and the Caribbean were recovered as sister taxa (Pp = 0.96), and D. cubanus was sister to them (Pp = 0.75). Divergence time estimates for these cladogenetic events were younger, although similar and with overlapping 95% HPD, than the estimates obtained for the fossil calibrated Cytb gene tree (Fig 2, S1 Fig). The main event that split the Atlantic and Pacific lineages was dated ~0.98 Mya (95% HPD: 0.45–1.34), and the divergence between D. lebretonis and the western Atlantic taxa was dated at 0.43 Mya (95% HPD: 0.28–0.63). Within the western Atlantic, mean divergence times ranged between 0.19 and 0.35 Mya.

thumbnail
Fig 5. Time-calibrated Dormitator species tree hypothesis, based on the multispecies coalescent analyses of Cytb, Rh and β-actin genes.

Numbers at the left of the nodes indicate the estimated mean TMRCA in Mya for each node, and horizontal grey bars at nodes represent the 95% HPD intervals for each estimated date. Node circles indicate posterior probabilities values (Pp) for the *BEAST analysis. The scale bar below the tree shows time in Mya. The colors of the clades correspond to the geographic origin of the sample locations (colored dots in Fig 1). Image of D. latifrons was obtained from the CPUM-Universidad Michoacana de San Nicolás de Hidalgo, México; voucher specimen No. 24632.

https://doi.org/10.1371/journal.pone.0153538.g005

Discussion

Species show different responses to common geological and environmental processes. As a consequence, their resulting evolutionary patterns may also differ because they depend on the relative roles of geological and ecological factors and their interaction with species-specific life-history traits [48,49]. These differences have been widely explored and compared across a broad range of taxa in Central America (see Bacon et al. [16] for a recent review), including fish (e.g. [5054]), and reconciled multiple pulses of dispersal and vicariance. The timing of these events was most likely associated not only with the variable availability of land and marine and freshwater corridors but also with the establishment of suitable climates and environments during the formation of the Isthmus [16]. However, how these factors affect the evolutionary patterns of amphidromous fish like Dormitator species had never been addressed before.

Species with an amphidromous life-history strategy would be expected to show different evolutionary patterns compared to strictly marine or freshwater species [10,55,56]. Contrary to this prediction, and despite the supposed capability of its larvae for long-distance dispersal through marine environments, Dormitator showed many evolutionary similarities with other freshwater fish, such as an absence of phylogeographic structure on the Pacific coast and a high phylogeographic structure in the Atlantic. Furthermore, the estimated time of vicariance between ocean basins was among the most recent reported to date, suggesting that gene flow existed between oceans until the very last stages of formation of the Central American Isthmus in the late Pliocene-Pleistocene period.

Pacific-Atlantic split

The closure of the Central American Isthmus is the major vicariant event that shaped the phylogenetic structure of Dormitator. The finding of a Pacific-Atlantic divergence rather than the expected Africa-America divergence rejects the hypothesis that the evolutionary patterns of Dormitator can be explained only by vicariance following the break-up of Gondwana. In contrast, this pattern suggests a widespread ancestral distribution of the genus or an American origin with dispersal to Africa following the closure of the Isthmus. Early transisthmian differentiation led to the formation of two main lineages of Dormitator: one lineage composed of the species from the Pacific basin (D. latifrons) and the other lineage composed of the species sampled in the Atlantic basin (D. maculatus, D. cubanus and D. lebretonis). This relationship was independently recovered by both mitochondrial and nuclear genes, even though the latter have been considered to evolve at a rate not fast enough to resolve transisthmian divergences [15]. Here, despite their slower substitution rate and higher effective population size, nuclear loci provided a clear pattern of divergence that indicated a complete interruption of gene flow between oceans.

The divergence between the Pacific and Atlantic Cytb lineages of Dormitator was estimated at between 0.3 and 3.1 Mya, which is consistent with a final closure of the Central American Seaway during the Late Pliocene (ca. 2.8 Mya, sensu Coates et al. [57]). The species tree analysis, on the other hand, estimated a younger date for this event (0.45-1-34 Mya). This pattern is expected according to the multispecies coalescent, which predicts that gene divergence will predate species divergence [58] and furthermore suggests that gene flow could have occurred across the Isthmus of Panama until very recent times during the Pleistocene. Divergence times between geminate species pairs, however, remain a subject of debate. According to Lessios [15], 30% of 115 species pairs of geminate clades (including echinoids, crustaceans, mollusks and fish) were likely to have diverged ca. 2.8 Mya, approximately 63% were separated at some point earlier during the long period of geological upheavals associated with the rising Isthmus, and 7% may have maintained genetic connections after the Isthmus closure. More recently, Bacon et al. [16] provided evidence that transisthmian divergence was a complex process that proceeded over at least two stages at approximately 24 Mya and 8 Mya. Furthermore, a significant decrease in the migration rates of marine organisms was detected at 2 Mya (1.03–4.35 Mya) [16], which is highly consistent with the divergence time estimated for Dormitator. Apparently, the large variation in divergence times is not predicted by intrinsic species biological factors. Rather, it seems that these differences are better explained by extrinsic features controlling habitat formation and availability during the geological development of the region [16]. This differentiation has been tested in the gastropod genera Cerithium and Cerithidea, in which geminate species inhabiting high intertidal mangrove habitats exhibit less evolutionary divergence than those that inhabit lower intertidal and subtidal habitats [17]. It is assumed that mangrove or estuarine habitats were the last to disappear during the Central American Seaway closure, thus allowing species that inhabited these environments, such as Dormitator, to maintain gene flow until the final closure of the seaway. Therefore, the divergence times of these species pairs most accurately correspond to the final completion of the Isthmus of Panama because they represent the milestone for the last connectivity events between oceans.

Pacific lineage

Dormitator in the Pacific showed a total absence of phylogeographic structure from Northern Mexico to Ecuador, indicating high gene flow across the distribution range of D. latifrons. Therefore, molecular data do not support the morphological differentiation of the subspecies of D. latifrons mexicanus in the eastern Pacific [59]. The genetic homogeneity across the Pacific distribution of Dormitator could be a product of their amphidromous life-history strategy, which potentially allows for long-distance dispersal during the marine larval stages, and their tolerance to salinity, which could also facilitate marine dispersal during adult stages [60,61]. Reduced genetic structure across the lower Central American Pacific coast has also been reported among primary freshwater fish species, which, contrary to amphidromous fish, do not tolerate salinity [51,54,62]. In such cases, it has been proposed that the lower sea levels during the last glacial maxima may have favored population dispersal and expansion events via river anastomosis [62]. If accepted, this hypothesis would indicate that the common pattern of low phylogeographical structure across the Pacific coast is not explained by biotic factors, such as amphidromy, but rather by abiotic factors, such as geology and climatic changes.

Atlantic lineage differentiation

The Atlantic lineage of Dormitator contained three species: D. maculatus, D. cubanus and D. lebretonis, all of which were reciprocally monophyletic, although their relationships showed conflicting results between the gene trees and the species tree analyses. For example, in the Cytb and concatenated gene trees, D. maculatus from the Gulf of Mexico and the Caribbean represented two non-sister clades, while the species tree recovered these clades as sisters with high support. Additionally, the large genetic divergence between the two clades of D. maculatus (Cytb Dp = 4.7 ± 0.6) might warrant a taxonomic revision because they might constitute cryptic species. Although cryptic evolutionary lineages are expected in widely distributed species, including amphidromous species [12], further analyses that incorporate missing species of the genus, such as D. lophocephalus from Suriname, could help to determine the validity of these species as well as their distribution ranges. The four Atlantic clades showed distinct and geographically delimited distributions, suggesting that phylogenetic structure may be the result of long-term vicariance and barriers to dispersal among basins within the Caribbean. This pattern is expected and consistent with the higher historical isolation of the Caribbean compared to the Pacific drainages [50,54,62].

Gulf of Mexico clade.

The youngest cladogenetic event within the Dormitator Atlantic lineage, as inferred by the multispecies coalescence analysis, corresponds to the divergence of the Gulf of Mexico and Caribbean clades approximately 0.19 Mya (95% HPD: 0.11–0.28 Mya). This divergence could have resulted from a combination of factors and a series of Pleistocene oceanographic changes that followed the closure of the Panama Isthmus. Oceanographic changes included variations of the sea level and changes in salinity, temperature and current patterns, such as the intensification of the Loop Current, the Florida Current, the Equatorial Undercurrent and the Gulf Stream [63]. All of these factors had major effects on marine paleobiogeography and paleoceanography [64] and may have also caused the isolation of populations in the Gulf of Mexico. In particular, the Loop Current has been proposed to play an important role as a barrier between the Gulf of Mexico and the Caribbean [10]. This current moves north between the Yucatán Peninsula and Cuba and recirculates to travel between Florida and Cuba [65]. Additionally, because currents along the coast of the Gulf of Mexico are not very strong, it is unlikely that Dormitator could be dispersed over long distances during their marine larval stage. In addition, the Yucatán Peninsula has been largely recognized as an area with a meager fish fauna because of its lack of surface river systems [22]. Consequently, this peninsula could have constituted a barrier to dispersal between the Caribbean and the Gulf of Mexico basins [10].

The large divergence observed in Dormitator for the Gulf of Mexico clade is also consistent with its delimitation as a separate biogeographic region from the Greater Caribbean, which is based on shore-fishes distribution data [66] and the genetic discontinuity observed in multiple co-distributed marine taxa [6770]. Recent climatic events are hypothesized to have strongly influenced the distributions and genetic signatures of taxa between the Atlantic and Gulf of Mexico ocean basins [68,71]. For example, the reduced sea levels associated with the glacial maxima might have led to the isolation of the western part of the Gulf of Mexico, which represents a refuge area for populations of estuarine species [67]. However, this expected phylogenetic separation between the Gulf of Mexico and the Atlantic Ocean basins is not observed in several near and offshore marine species and may not be, due to these species’ greater dispersal capabilities [67] (also see Soltis et al. [70]). Thus, it appears that vicariant scenarios alone cannot explain this phylogenetic break and the current phylogeographic pattern is likely the result of the interplay among the shared paleoclimatic history, ecological factors and species-specific life-history traits [67,68,70].

Cameroon.

The divergence of D. lebretonis from the western Atlantic species occurred between 0.28 and 0.63 Mya. The short elapsed time since the split, as well as the retained ancestral polymorphisms or recent trans-Atlantic gene flow, could explain why this divergence was only recovered by the mitochondrial and combined data sets and not by any combination of the nuclear markers analyzed [12,15]. Similarly, differences in effective population size and gene flow rates between populations on both sides of the Atlantic might explain the discordant results of the gene and species trees [58]. Unfortunately, these hypotheses remain untested because only one sample of D. lebretonis could be included in the study. Therefore, these results must be interpreted cautiously.

The sister relationship of the African D. lebretonis and the rest of the American Atlantic clades both rejects the vicariant hypothesis of a Gondwanan origin of Dormitator and postdates by almost 20 million years Rosen’s [72] Eastern Atlantic/Western Atlantic generalized track hypothesis, proposing that sister species distributed across the Atlantic regions were the result of Cenozoic (65–20 million years ago) spreading across the Atlantic basin. Several groups of fish share this amphi-Atlantic distribution [7376], and despite genetic differentiation between Atlantic coasts, certain authors consider the mid-Atlantic Barrier (> 3500 km of deep ocean) to act as a soft barrier for marine species [76]. In fact, patterns showing a greater resemblance between taxa on both sides of the Atlantic than between the western Atlantic and eastern Pacific are not uncommon for estuarine taxa. For example, for the mangroves Avicennia germinans and Rhizophora mangle, which also show current trans-Atlantic gene flow [77], a strong role of contemporary currents has been proposed in shaping their genetic landscape. The main currents are the North Equatorial Current and South Equatorial Current, which travel from east to west, and the North Equatorial Counter Current, which travels from west to east. All of them were established following the breakup of West Gondwana and have maintained the same direction connecting the Americas with West Africa since the final closure of the Central American Isthmus [7881]. Because of the ecological and biogeographical resemblances among Rhizophora, Avicennia and Dormitator, it is likely that these species, as well as other amphidromous or estuarine fish (e.g., Eleotris), are affected by similar processes that allow long-distance dispersal through marine current drift. However, the degree of isolation and contemporary gene flow across the Atlantic should be addressed by examining a wider sample of African populations.

Caribbean and West Cuba clades.

The island of Cuba harbored two clades of Dormitator, with non-overlapping distributions. The West Cuba clade was restricted to Juventud Island and Pinar del Río Province locations in the western region of Cuba, which fits with the distribution of D. cubanus according to Ginsburg [59]. The Caribbean clade included the eastern region of Cuba as well as Dominican Republic, Honduras, Nicaragua, Panama and Venezuela, and corresponded to the distribution range of D. maculatus.

The gene trees showed a closer relationship between the West Cuba clade and the Caribbean clade, while the species tree showed West Cuba as the sister taxon of D. maculatus from the Caribbean and the Gulf of Mexico clades. Additionally, the West Cuba clade was the only clade recovered by the nuclear gene trees as distinct within the Atlantic basin. This discrepancy between molecular markers and methods might suggest a selective or demographic process through which nuclear alleles evolved or were sorted faster than mitochondrial alleles. One possibility is that this process could have been favored by drift in a population with a small effective size after the isolation of western Cuba as a result of climatic changes, current shifts and sea level variations [82].

Considering the large and homogeneous distributions of most species or clades of Dormitator, the small geographic scale of differentiation and the highly restricted distribution of D. cubanus is remarkable. Yet the pattern and timing of cladogenesis in Cuban Dormitator match previous biogeographic studies on fish and amphibians that support the paleogeographic hypothesis that Juventud Island has been isolated from Cuba at least since the Oligocene and more recently during the Pleistocene period [72,8386]. Within the endemic Cuban poeciliid tribe Girardinini, the isolation of Juventud Island and western Cuba (Guanahacabibes Peninsula at Pinar del Río Province) has been proposed as a vicariant event promoting speciation in this group of fish [85]. For example, the presence of Quintana, the sister group of all Girardinini, in Juventud Island could explain an ancient isolation event between the two islands. Additionally, a more recent isolation during the late Pliocene-Pleistocene would be supported by the well-differentiated populations of Glaridichthys falcatus in Juventud Island and central Cuba, and by Girardinus rivasi, which is endemic to Juventud Island [85]. These two species and populations of Juventud Island present Cytb divergences from their closest relatives in Cuba (Dp = 1.7 to 2.4% between G. falcatus populations and Dp = 1.1 to 1.8% between G. rivasi and G. microdactylus) that were similar to those observed between D. cubanus and D. maculatus (Dp = 2.8±0.4%), suggesting coeval cladogenetic events. A concordant pattern was observed in the Peltophryne toad radiation, with reciprocally monophyletic populations in Juventud Island and Cuba (Peltophryne empusa and Peltophryne peltocephala) and species with distributions limited to western Cuba (Peltophryne fustiger). In this case, genetic divergences were shallow but significant (16S Dp = 0.8–0.4%, COI Dp = 2.1–0.6%) and suggest a recent Pleistocene-Holocene separation between the islands of Juventud and Cuba [86]. A similar speciation pattern might also exist in marine taxa from Juventud Island. For example, the coral reef fish Gramma dejongi was described as an endemic species from this region, but contrary to Dormitator, it shows no mitochondrial genetic differentiation relative to co-occurring congeneric species. Therefore, this difference has been proposed as a case of incipient speciation or even a local color variant of the sympatric and widespread Gramma loreto [87].

The deep phylogenetic differentiation observed in the Atlantic lineage supports the current species of Dormitator and provides evidence for possible cryptic species. Even the lowest sequence divergence among Atlantic clades, which was observed between the West Cuba and Caribbean clades (Cytb Dp = 2.8%), largely exceeds the 2% sequence divergence generally accepted as a cut-off value between sister species of vertebrates [88]. Our molecular data support the validity of D. cubanus, the distribution of which is consistent with the original description of the species [59] from Pinar del Río and extends its range to Juventud Island. The parapatric occurrence of D. cubanus with D. maculatus in Cuba was also noted, although more recent works had exclusively indicated the presence of D. maculatus in all the Cuban territory [89,90].

Conclusion

This work is the first study to address the evolutionary history of the amphidromous fish genus Dormitator. Despite the expected long-distance dispersal capability of these species based on their amphidromous life-style, vicariant geological processes rather than life-history traits appear to be the main drivers of speciation in Dormitator. First, the presence of geminate lineages of Dormitator across the Pacific and Atlantic slopes confirms the importance of the Central American Isthmus in the evolutionary history of the genus. Second, the four geographically delimited Atlantic clades appear to be shaped by isolation processes that are largely congruent with geologic and oceanographic events related to the closure of the Central America Isthmus. The high genetic structure and isolation of the Caribbean basin is consistent with patterns observed in both primary freshwater and marine shore fish and further contrasts with the genetic homogeneity along the Pacific coast. This pattern is shared among other species with much lower dispersal capabilities and is likely caused by enhanced river connectivity during periods with lower sea levels in the last glacial maximum. Comparative studies of other fish sympatric to Dormitator could provide further information on the relative roles of abiotic (i.e., geological and climatic) and biotic processes (i.e., species ecology and life-history traits) in the evolutionary history and biogeographic patterns of amphidromous neotropical fishes.

Supporting Information

S1 Fig. Time-calibrated phylogenetic reconstruction of Dormitator and 37 Gobiiformes taxa based on mitochondrial Cytb sequences.

Within Dormitator (colored area), horizontal grey bars represent the 95% highest posterior density of the estimated time to the most recent common ancestor (TMRCA) in Mya for each node. Numbers at the left of the grey bars indicate the estimated mean TMRCA. The colors of the clades correspond to the geographic origin of the sample locations (as in Fig 1). Bullets represent posterior probability support for each node. Black bullets represent Pp = 1, and white bullets represent Pp = 0.6. Red circles represent the four fossil calibrations used to time-constrain each node in the molecular clock analysis (S2 Table). The scale bar below the tree shows time in Mya.

https://doi.org/10.1371/journal.pone.0153538.s001

(EPS)

S1 Table. Fossil species of Gobiiformes used to calibrate the Dormitator molecular clock analysis.

https://doi.org/10.1371/journal.pone.0153538.s002

(DOCX)

S2 Table. Species, GenBank accession number and fossil calibration points for each individual included in the Dormitator molecular clock analysis based on cytochrome b (Cytb) sequences.

https://doi.org/10.1371/journal.pone.0153538.s003

(DOCX)

Acknowledgments

We would like to acknowledge all of the people who helped with field trips and/or laboratory work: Paloma Garzón, Adolfo de Sostoa, Silvia Perea, Lourdes Alcaraz, Diushi Keri Corona, Rodolfo Pérez, Salvador Romero, Georgina Palacios, Carmen Pedraza, Eloísa Torres, Edgar Sandoval, Fernando Mar, Adrián García, Alfrancis Arredondo, Carlos Pinacho, Arturo Angulo, Enrique Barraza, Oscar Lasso, Eyda Gómez, Heidi Banford, Claudia De Jesús and Xavier Madrigal. We thank A.J. Turner for his help revising the English. We also thank the Agencia Nacional del Ambiente de Panamá and all the authorities and institutions that issued the permits necessary to collect specimens for this research. The authors appreciate the helpful comments and discussion provided by the academic editor and two anonymous reviewers that significantly improved this paper.

Author Contributions

Conceived and designed the experiments: SGQ ID NH ODD. Performed the experiments: SGQ FA AP RGR. Analyzed the data: SGQ ID FA ODD. Contributed reagents/materials/analysis tools: ID MGV ACM EB ODD NH. Wrote the paper: SGQ ID FA ODD. Coordinated the study and supervised the SGQ PhD thesis: ODD. Provided molecular laboratory space: ID EB ODD. Provided technical support: ID RGR MGV NH ACM EB ODD. Conducted field work: SGQ ID AP MGV NH ACM EB ODD. Discussed the results: SGQ ID FA AP RGR MGV ACM EB ODD.

References

  1. 1. Shulman MJ, Bermingham E. Early life histories, ocean currents, and the population genetics of Caribbean reef fishes. Evolution. 1995;49: 897–910.
  2. 2. Chubb AL, Zink RM, Fitzsimons JM. Patterns of mtDNA variation in Hawaiian freshwater fishes: the phylogeographic consequences of amphidromy. J Hered. 1998;89: 8–16. pmid:9487675
  3. 3. McDowall RM. Diadromy, diversity and divergence: implications for speciation processes in fishes. Fish Fish. 2001;2: 278–285.
  4. 4. McDowall RM. Diadromy in fishes: migrations between freshwater and marine environments. London, UK: Croom Helm; 1988.
  5. 5. McDowall RM. Diadromy and genetic diversity in Nearctic and Palearctic fishes. Mol Ecol. 1999;8: 527–528.
  6. 6. Leathwick JR, Elith J, Chadderton WL, Rowe D, Hastie T. Dispersal, disturbance and the contrasting biogeographies of New Zealand’s diadromous and non-diadromous fish species. J Biogeogr. 2008;35: 1481–1497.
  7. 7. Darwin C. On the origin of species by means of natural selection, or the preservation of favoured races in the struggle for life. 6th ed. London: Murray; 1872.
  8. 8. Béguer M, Beaulaton L, Rochard E. Distribution and richness of diadromous fish assemblages in Western Europe: large-scale explanatory factors. Ecol Freshwat Fish. 2007;16: 221–237.
  9. 9. Cook BD, Bernays S, Pringle CM, Hughes JM. Marine dispersal determines the genetic population structure of migratory stream fauna of Puerto Rico: evidence for island-scale population recovery processes. J North American Benthol Soc. 2009;28: 709–718.
  10. 10. McMahan CD, Davis MP, Domínguez-Domínguez O, García de León FJ, Doadrio I, Piller KR. From the mountains to the sea: phylogeography and cryptic diversity within the mountain mullet, Agonostomus monticola (Teleostei: Mugilidae). J Biogeogr. 2013;40: 894–904.
  11. 11. Violante-González J, Rojas-Herrera A, Aguirre-Macedo ML. Seasonal patterns in metazoan parasite community of the “fat sleeper” Dormitator latifrons (Pisces: Eleotridae) from Tres Palos lagoon, Guerrero, México. Rev Biol Trop. 2008;56: 1419–1427. pmid:19419054
  12. 12. Taillebois L, Castelin M, Ovenden JR, Bonillo C, Keith P. Contrasting genetic structure among populations of two amphidromous fish species (Sicydiinae) in the Central West Pacific. PLOS ONE. 2013;8: e75465. pmid:24130714
  13. 13. Jordan DS. The law of geminate species. Am Nat. 1908;42: 73–80.
  14. 14. Bermingham E, McCafferty SS, Martin AP. Fish biogeography and molecular clocks: perspectives from the Panamian isthmus. In: Kocher TD, Stepien CA, editors. Molecular systematics of fishes. San Diego, California, USA: Academic Press; 1997. pp. 113–128.
  15. 15. Lessios HA. The great American schism: divergence of marine organisms after the rise of the Central American isthmus. Annu Rev Ecol Evol Syst. 2008;39: 63–91.
  16. 16. Bacon CD, Silvestro D, Jaramillo C, Smith BT, Chakrabarty P, Antonelli A. Biological evidence supports an early and complex emergence of the isthmus of Panama. Proc Natl Acad Sci U S A. 2015;112: 6110–6115. pmid:25918375
  17. 17. Miura O, Torchin ME, Bermingham E. Molecular phylogenetics reveals differential divergence of coastal snails separated by the Isthmus of Panama. Mol Phylogenet Evol. 2010;56: 40–48. pmid:20399869
  18. 18. Knowlton N, Weigt LA. New dates and new rates for divergence across the isthmus of Panama. Proc R Soc London Ser B. 1998;265: 2257–2263.
  19. 19. Myers GS. Usage of anadromous, catadromous, and allied terms for migratory fishes. Copeia. 1949;2: 89–97.
  20. 20. Nordlie FG, Haney DC. Euryhaline adaptations in the fat sleeper Dormitator maculatus. J Fish Biol. 1993;43: 433–439.
  21. 21. Nelson DJ. Fishes of the world. 2nd ed. New York: John Wiley and Sons; 1994.
  22. 22. Miller RR, Minckley WL, Norris SM. Freshwater fishes of México. Chicago: The University of Chicago Press; 2005.
  23. 23. Massay S, Mosquera R. Presence of chame Dormitator latifrom (Richardson, 1844) (Pisces: Eleotridae), in the Galapagos Islands, Ecuador. J Fish Biol. 1992;40: 815–816.
  24. 24. Harrison J, Miller P. Eleotridae. In: Lévêque C, Paugy D, Teugels GG, editors. Faune des poisons d’eaux douces et saumâtres d’Afrique de l’Ouest. Tomo 2. Paris, France: ORSTOM Editions; 1992.
  25. 25. Kullander SO. Gobiidae (Gobies). In: Reis RE, Kullander SO, Ferraris CJ Jr, editors. Checklist of the freshwater fishes of South and Central America. Porto Alegre, Brasil: EDIPUCRS; 2003. pp. 657–665.
  26. 26. Eschmeyer WN, editor. Catalog of fishes. California Academy of Sciences; 2015. Available: http://researcharchive.calacademy.org/research/ichthyology/catalog/fishcatmain.asp. Accessed 11 June 2015.
  27. 27. Froese R, Pauly D, editors. FishBase. 2015. Available: www.fishbase.org. Accessed 11 June 2015.
  28. 28. Akihito Iwata AA, Kobayashi T, Ikeo K, Imanishi T, Ono H, et al. Evolutionary aspects of gobioid fishes based upon a phylogenetic analysis of mitochondrial cytochrome b genes. Gene. 2000;259: 5–15. pmid:11163956
  29. 29. Thacker CE, Hardman MA. Molecular phylogeny of basal gobioid fishes: Rhyacichthyidae, Odontobutidae, Xenisthmidae, Eleotridae (Teleostei: Perciformes: Gobioidei). Mol Phylogenet Evol. 2005;37: 858–871. pmid:15975831
  30. 30. Agorreta A, San Mauro D, Schliewen U, Van Tassell JL, Kovačić M, Zardoya R, et al. Molecular phylogenetics of Gobioidei and phylogenetic placement of European gobies. Mol Phylogenet Evol. 2013;69: 619–633. pmid:23911892
  31. 31. Follett WI. The fresh-water fishes-their origins and affinities. Syst Zool. 1960;9: 212–232.
  32. 32. Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual. Cold Spring Harbor, NY: Cold Spring Harbor Laboratory; 1989.
  33. 33. Palumbi SR. Nucleic acids II: The polymerase chain reaction. In: Hillis DM, Moritz C, Mable BK, editors. Molecular systematics. 2nd ed. Sunderland MA: Sinauer Associates; 1996. pp. 205–247.
  34. 34. Perdices A, Bermingham E, Montilla A, Doadrio I. Evolutionary history of the genus Rhamdia (Teleostei: Pimelodidae) in Central America. Mol Phylogenet Evol. 2002;25: 172–189. pmid:12383759
  35. 35. Chen WJ, Bonillo C, Lecointre G. Repeatability of clades as a criterion of reliability: a case study for molecular phylogeny of Acanthomorpha (Teleostei) with larger number of taxa. Mol Phylogenet Evol. 2003;26: 262–288. pmid:12565036
  36. 36. Robalo JI, Sousa Santos C, Levy A, Almada VC. Molecular insights on the taxonomic position of the paternal ancestor of the Squalius alburnoideshybridogenetic complex. Mol Phylogenet Evol. 2006;39: 276–281. pmid:16213165
  37. 37. Hall TA. BioEdit: a user-friendly biological sequence alignment Editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41: 95–98.
  38. 38. Darriba D, Taboada GL, Doallo R, Posada D. JModelTest 2: More models, new heuristics and parallel computing. Nat Methods. 2012;9: 772.
  39. 39. Stamatakis A, Hoover P, Rougemont J. A rapid bootstrap algorithm for the RAxML web servers. Syst Biol. 2008;57: 758–771. pmid:18853362
  40. 40. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61: 539–542. pmid:22357727
  41. 41. Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES science gateway for inference of large phylogenetic trees. In: Proceedings of the Gateway Computing Environments Workshop (GCE), 14 Nov. 2010, New Orleans, LA. Piscataway, NJ: IEEE; 2010. pp. 1–8.
  42. 42. Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v 1.6. 2014. Available: http://beast.bio.ed.ac.uk/Tracer.
  43. 43. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30: 2725–2729. pmid:24132122
  44. 44. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUTi and the BEAST 1.7. Mol Biol Evol. 2012;29: 1969–1973. pmid:22367748
  45. 45. Drummond AJ, Ho SY, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLOS Biol. 2006;4: e88. pmid:16683862
  46. 46. Yang Z, Rannala B. Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds. Mol Biol Evol. 2006;23: 212–226. pmid:16177230
  47. 47. Heled J, Drummond AJ. Bayesian inference of species trees from multilocus data. Mol Biol Evol. 2010;27: 570–580. pmid:19906793
  48. 48. Bloom DD, Lovejoy NR. The evolutionary origins of diadromy inferred from a time-calibrated phylogeny for Clupeiformes (herring and allies). Proc Biol Sci. 2014;281: 20132081. pmid:24430843
  49. 49. Paz A, Ibáñez R, Lips KR, Crawford AJ. Testing the role of ecology and life history in structuring genetic variation across a landscape: a trait-based phylogeographic approach. Mol Ecol. 2015;24: 3723–3737. pmid:26080899
  50. 50. Bermingham E, Martín AP. Comparative mtDNA phylogeography of neotropical freshwater fishes: testing shared history to infer the evolutionary landscape of lower Central America. Mol Ecol. 1998;7: 499–517. pmid:9628002
  51. 51. Reeves RG, Bermingham E. Colonization, population expansion, and lineage turnover: phylogeography of Mesoamerican characiform fish. Biol J Linn Soc. 2006;88: 235–255.
  52. 52. 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: 460–473. pmid:16223677
  53. 53. Ornelas-García CP, Domínguez-Domínguez O, Doadrio I. Evolutionary history of the fish genus Astyanax Baird & Girard (1854) (Actinopterygii, Characidae) in Mesoamerica reveals multiple morphological homoplasies. BMC Evol Biol. 2008;8: 340. pmid:19102731
  54. 54. Picq S, Alda F, Krahe R, Bermingham E. Miocene and Pliocene colonization of the Central American isthmus by the weakly electric fish Brachyhypopomus occidentalis (Hypopomidae, Gymnotiformes). J Biogeogr. 2014;41: 1520–1532.
  55. 55. Tringali MD, Bert TM, Seyoum S, Bermingham E, Bartolacci D. Molecular phylogenetics and ecological diversification of the transisthmian fish genus Centropomus (Perciformes: Centropomidae). Mol Phylogenet Evol. 1999;13: 193–207. pmid:10508552
  56. 56. Craig MT, Hastings PA, Pondella DJ II. Speciation in the Central American seaway: the importance of taxon sampling in the identification of trans-isthmian geminate pairs. J Biogeogr. 2004;31: 1085–1091.
  57. 57. Coates AG, McNeill DF, Aubry MP, Berggren WA, Collins LS. An introduction to the geology of the Bocas del Toro archipelago, Panama. Carib J Sci. 2005;41: 374–391.
  58. 58. Edwards SV, Beerli P. Perspective: gene divergence, population divergence, and the variance in coalescence time in phylogeographic studies. Evolution. 2000;54: 1839–1854. pmid:11209764
  59. 59. Ginsburg I. Ten new American gobioid fishes in the United States National Museum, including additions to a revision of Gobionellus. Jour Wash Acad Sci. 1953;43: 18–26.
  60. 60. Ancieta DF, Landa A. Reseña taxonómica y biológica de los peces cultivados en el área andina incluyendo la Costa del Perú. FAO Inf Pesca. 1977;2: 106–113.
  61. 61. Castro-Aguirre JL, Schmitter-Soto JJ, Espinoza Pérez H. Ictiofauna estuarino-lagunar y vicaria de México. CICIMAR-IPN. México, D. F.: Limusa Noriega; 1999.
  62. 62. Smith SA, Bermingham E. The biogeography of lower Mesoamerican freshwater fishes. J Biogeogr. 2005;32: 1835–1854.
  63. 63. Schmidt DN. The closure history of the Central American seaway: evidence from isotopes and fossils to models and molecules. In: William M, Haywood AM, Gregory FJ, Schmidt DN, editors. Deep-time perspectives on climatic change: marrying the signal from computer models and biological proxies. London, UK: Geological Society Publishing House; 2007. pp. 427–442.
  64. 64. Cronin TM, Dowsett HJ. Biotic and oceanographic response to the Pliocene closing of the Central American isthmus. In: Jackson JBC, Budd AF, Coates AG, editors. Evolution and environment in tropical America. Chicago, IL: Chicago University Press; 1996. pp. 76–104.
  65. 65. Coats DA. Chapter 6. The loop current. In: Milliman JD, Imamura E, editors. The physical oceanography of the U.S. Atlantic and Eastern Gulf of Mexico. Herndon, VA, USA: U.S. Department of the Interior, Mineral Management Service, Atlantic OCS Region; 1992.
  66. 66. Robertson DR, Cramer KL. Defining and dividing the greater Caribbean: insights from the biogeography of shorefishes. PLOS ONE. 2014;9: e102918. pmid:25054225
  67. 67. Bowen BW, Avise JC. Genetic structure of Atlantic and Gulf of Mexico populations of sea bass, menhaden, and sturgeon: influence of zoogeographic factors and life-history patterns. Mar Biol. 1990;107: 371–381.
  68. 68. Avise JC. Phylogeography, the history and formation of species. Cambridge, MA: Harvard University Press; 2000.
  69. 69. Palumbi SR. Genetic divergence, reproductive isolation, and marine speciation. Annu Rev Ecol Syst. 1994;25: 547–572.
  70. 70. Soltis DE, Morris AB, McLachlan JS, Manos PS, Soltis PS. Comparative phylogeography of unglaciated eastern North America. Mol Ecol. 2006;15: 4261–4293. pmid:17107465
  71. 71. Avise JC. Population structure and biogeographic history of a regional fauna: a case history with lessons for conservation biology. Okies. 1992;63: 62–76.
  72. 72. Rosen DE. Geological hierarchies and biogeographic congruence in the Caribbean. Ann Mo Bot Gard. 1985;72: 636–659.
  73. 73. Haug GH, Tiedemann R. Effect of the formation of the Isthmus of Panama on Atlantic Ocean thermohaline circulation. Nature. 1998;393: 673–676.
  74. 74. Muss A, Robertson DR, Stepien CA, Wirtz P, Bowen BW. Phylogeography of Ophioblennius: the role of ocean currents and geography in reef fish evolution. Evolution. 2001;55: 561–572. pmid:11327163
  75. 75. Carlin JL, Robertson DR, Bowen BW. Ancient divergences and recent connections in two tropical Atlantic reef fishes Epinephelus adscensionis and Rypticus saponaceous (Percoidei: Serranidae). Mar Biol. 2003;143: 1057–1069.
  76. 76. Floeter SR, Rocha LA, Robertson DR, Joyeux JC, Smith-Vaniz WF, Wirtz P, et al. Atlantic reef fish biogeography and evolution. J Biogeogr. 2008;35: 22–47.
  77. 77. Cerón-Souza I, González EG, Schwarzbach AE, Salas-Leiva DE, Rivera-Ocasio E, Toro-Perea N, et al. Contrasting demographic history and gene flow patterns of two mangrove species on either side of the Central American isthmus. Ecol Evol. 2015;5: 3486–3499. pmid:26380680
  78. 78. Richardson PL, Walsh D. Mapping climatological seasonal variations of surface currents in the tropical Atlantic using ship drifts. J Geophys Res. 1986;91: 10537–10550.
  79. 79. Parrish JT. The paleogeography of the opening south Atlantic. In: Georgeand W, Lavocat R, editors. The Africa-South America connection. Oxford: Clarendon Press; 1993. pp. 8–41.
  80. 80. Stramma L, Schott F. The mean flow field of the tropical Atlantic Ocean. Deep Sea Res Part II Topical Stud Oceanogr. 1999;46: 279–303.
  81. 81. Renner SS. Plant dispersal across the tropical Atlantic by wind and sea currents. Int J Plant Sci. 2004;16S: S23–S33.
  82. 82. Taylor MS, Hellberg ME. Genetic evidence for local retention of pelagic larvae in a Caribbean reef fish. Science. 2003;299: 107–109. pmid:12511651
  83. 83. Rauchenberger M. Historical biogeography of Poeciliid fishes in the Caribbean. Syst Biol. 1988;37: 356–365.
  84. 84. Iturralde-Vinent MA, MacPhee RDE. Paleogeography of the Caribbean region: implications for Cenozoic biogeography. Bull Am Mus Nat Hist. 1999;238: 1–95.
  85. 85. Doadrio I, Perea S, Alcaraz L, Hernández N. Molecular phylogeny and biogeography of the Cuban genus Girardinus Poey, 1854 and relationships within the tribe Girardinini (Actinopterygii, Poeciliidae). Mol Phylogenet Evol. 2009;50: 16–30. pmid:18854217
  86. 86. Alonso R, Crawford AJ, Bermingham E. Molecular phylogeny of an endemic radiation of Cuban toads (Bufonidae:Peltophryne) based on mitochondrial and nuclear genes. J Biogeogr. 2012;39: 434–451.
  87. 87. Victor BC, Randall JE. Gramma dejongi, a New Basslet (Perciformes: Grammatidae) from Cuba, a sympatric sibling species of G. loreto. Zool. Stud. 2010;49: 865–871.
  88. 88. Johns GC, Avise JC. A comparative summary of genetic distances in the vertebrates from the mitochondrial cytochrome b gene. Mol Biol Evol. 1998;15: 1481–1490. pmid:12572611
  89. 89. Vergara RR. Principales características de la ictiofauna dulceacuícola cubana. La Habana, Cuba: Editorial Academia; 1992.
  90. 90. Lara A, Ponce de León JL, Rodríguez R, Casane D, Côté G, Bernatchez L, et al. DNA barcoding of Cuban freshwater fishes: evidence for cryptic species and taxonomic conflicts. Mol Ecol Resour. 2010;10: 421–430. pmid:21565041