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

Multigene phylogeny supports diversification of four-eyed fishes and one-sided livebearers (Cyprinodontiformes: Anablepidae) related to major South American geological events

  • Pedro F. Amorim ,

    Contributed equally to this work with: Pedro F. Amorim, Wilson J. E. M. Costa

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    pedro_f_a@hotmail.com

    Affiliation Laboratory of Systematics and Evolution of Teleost Fishes, Institute of Biology, Federal University of Rio de Janeiro, Rio de Janeiro, Rio de Janeiro, Brazil

  • Wilson J. E. M. Costa

    Contributed equally to this work with: Pedro F. Amorim, Wilson J. E. M. Costa

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Laboratory of Systematics and Evolution of Teleost Fishes, Institute of Biology, Federal University of Rio de Janeiro, Rio de Janeiro, Rio de Janeiro, Brazil

Abstract

The high diversity of Neotropical fishes has been attributed to major South American palaeogeographic events, such as Andean uplift, rise of the Isthmus of Panama and marine transgressions. However, the unavailability of temporal information about evolution and diversification of some fish groups prevents the establishment of robust hypotheses about correlations between species diversification and proposed palaeogeographical events. One example is the Anablepidae, a family of teleost fishes found mostly in coastal habitats of Central and South America, but also in some inner river basins of South America. Historical aspects of the distribution patterns of the Anablepidae were never analysed and no accurate estimation of time of its origin and diversification is presently available. A multi-gene analysis was performed to estimate Anablepidae phylogenetic position, age and biogeography, comprising seven nuclear genes. The suborder Cyprinodontoidei was recovered in three major clades, one comprising all the Old World Cyprinodontoidei and two comprising New World lineages. Anablepidae was recovered as the sister group of the New World Poeciliidae, with the Amazonian genus Fluviphylax as their sister group. The ages found for the origin and diversification of Cyprinodontiformes were congruent with the pattern recorded for other vertebrate groups, with an origin anterior to the Cretaceous-Paleogene (K-Pg) transition and diversification during the Paleogene. The age estimated for the split between the Atlantic and Pacific lineages of Anableps was congruent with the rise of Panamanian Isthmus. The results suggest Miocene marine transgressions as determinant to the current distribution of Jenynsia.

Introduction

The family Anablepidae is a group of cyprinodontiform fishes comprising three recent genera, Anableps Scopoli, 1777, Jenynsia Günther, 1866 and Oxyzygonectes Fowler, 1916, commonly known as four-eyed fishes, one-sided livebearers, and white-eye, respectively [1]. Two species of Jenynsia and all those included in Anableps and Oxyzygonectes live in coastal salt-water or brackish-water habitats, whereas the remaining 12 species of Jenynsia inhabit freshwaters [2, 3, 4, 5, 6].

Regarding the fossil record, a monotypic extinct genus, Carrionellus White, 1927 from the Early Miocene of Loja, Ecuador, was first described as a member of the Cyprinodontidae [7], and later transferred to Anablepidae species [2]. However, after a detailed analysis of the type series of Carrionellus diumortuus White, 1927, the genus was then realocated in the Cyprinodontidae family [8]. In a recent study two new monotypic genera from Middle Miocene of Northern Argentina were described in the family Anablepidae, Tucumanableps Sferco, Herbst, Aguilera & Mirande, 2017, recovered as the sister group of Anableps, and Sachajenynsia Sferco, Herbst, Aguilera & Mirande, 2017, the sister of all other genera of Anablepidae except Oxyzygonectes [9].

Species of Anableps exhibit one of the oddest adaptations among Cyprinodontiformes, comprising prominent eyes that rise above their heads, with horizontally divided pupils with an hourglass-like shape. This unique eye morphology allows them to simultaneously see aerial and aquatic environments when swimming with the middle line of the eye at the water surface [10]. These unusual morphological and behavioural traits attracted attention since the first European travellers arrived in America, so records of Anableps date from the early 17th century [11], more than 150 years prior to the formal description of the first species of the genus [12]. The complexity of their eye makes these animals important organisms for studies of embryogeny and gene expression (e.g. [13, 14, 15]). Furthermore, due to the presence of a modified anal fin of males in a tubular gonopodium and a placenta-like tissue in females, which together allow internal fertilization and viviparity [16, 17], Anableps and Jenynsia are also models of studies involving reproduction and ontogeny (e.g. [18, 19, 20]).

The geographical distribution of Anablepidae is restricted to Neotropical areas (Fig 1): Oxyzygonectes is a monotypic genus found in the Pacific coast of Central America; Anableps comprises three species inhabiting both the Pacific coast of Central America and the Atlantic coast of northern South America; and Jenynsia, the most species diverse genus of Anablepidae, comprising 14 species inhabiting the Paraná-Paraguay river system and coastal drainages of south-eastern South America [2, 4, 21]. However, biogeographical historical events responsible for the present distribution of the Anablepidae were never deeply explored [2]. Previous authors discussed the distribution of Anableps which could be related to the rise of the Panama Isthmus or ancient events, but no conclusion was proposed [1, 2, 5]. Jenynsia presents a huge disjunct distribution when compared to other anablepid genera (Fig 1), but there is no available explanation for this distribution pattern based on analytical tools for biogeographical inference. Here we first provide a time-calibrated multi-gene phylogeny of the Anablepidae, discussing two possible past historical scenarios that have shaped the present distribution of the main Anablepidae components: (1) invasion of the most recent common ancestor into the Amazon basin through Miocene marine transgressions followed by dispersion to south through ancient connections between the Amazon and Paraná river basins or (2) a dispersion along the eastern South America coast with an invasion to the freshwater habitats in the mouth of the La Plata river, during Miocene, through marine transgression forming the Paranean Sea.

thumbnail
Fig 1. Geographical distribution of the genera of Anablepidae.

Oxyzygonectes in red, Anableps in yellow and Jenynsia in green.

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

Age estimates for the Anablepidae and other Neotropical cyprinodontoid groups are still uncertain. Previous studies proposed the origin of the clade comprising Anablepidae and Poeciliidae (sensu [1]) dating back to the Lower Cretaceous (130–110 Ma), before separation between Africa and South America [22, 23, 24], but a recent time-calibrated analysis [25] estimated the origin of Poecilinae and Anablepidae clade between the Upper Cretaceous and the Early Palaeocene (62–67 Ma). In the current study, we analysed the largest dataset of genetic information for Anablepidae to test its phylogenetic position among the Cyprinodontoidei and to estimate time of origin and diversification of the family and their genera using three calibrating points based on well-known cyprinodontoid fossils.

Materials and methods

Taxonomic sampling, DNA extraction, PCR, and sequencing

A total of 82 terminal taxa were used, including one species of Beloniformes, the sister group to Cyprinodontiformes (e.g. [26, 27]), in which the analyses were rooted, and 81 species of Cyprinodontiformes, among which are represented all species of Anableps and Oxyzygonectes, and eight of the 14 species of Jenynsia. Molecular data were obtained from specimens euthanized with buffered solution of ethyl-3-amino-benzoat-methansulfonat (MS-222) at a concentration of 250 mg/l, for a period of 10 minutes or more, until completely ceasing opercular movements, following the methods for euthanasia approved by CEUA-CCS-UFRJ (Ethics Committee for Animal Use of Federal University of Rio de Janeiro; permit number: 01200.001568/2013-87). These individuals were fixed and preserved in absolute ethanol. The dataset was complemented with sequences obtained from GenBank, a complete list of taxa used in the analysis is presented in S1 Table. Seven nuclear genes were sampled: ectodermal-neural cortex (ENC1), glycosyltransferase (GLYT), cardiac muscle myosin heavy chain 6 alpha (MYH6), recombination activating gene 1 (RAG1), rhodopsin (RHO), SH3 and PX domain-containing 3-like protein (SH3PX3) and tyrosine kinase exons 8–10 and introns 8–9 (X-SRC). The introns of X-SRC were excluded due to the high degree of variability found, which produce highly ambiguous alignments.

Total genomic DNA was extracted from muscle tissue of the right side of the caudal peduncle using the DNeasy Blood & Tissue Kit (Qiagen), according to the manufacturer instructions. To amplify the fragments of the DNA, the primers provided in the S2 Table were used. Polymerase chain reactions (PCR) were performed in 30μl reaction mixtures containing 5x Green GoTaq Reaction Buffer (Promega), 3.2 mM MgCl2, 1 μM of each primer, 75 ng of total genomic DNA, 0.2 mM of each dNTP and 1U of Taq polymerase. The thermocycling profile was: (1) 1 cycle of 4 minutes at 94°C; (2) 35 cycles of 1 minute at 92°C, 1 minute at 50–63°C and 1 minute at 72°C; and (3) 1 cycle of 4 minutes at 72°C. In all PCR reactions, negative controls without DNA were used to check contaminations. Amplified PCR products were purified using the Wizard SV Gel and PCR Clean-Up System (Promega). Sequencing reactions were made by Macrogen Inc. (South Korea) using the BigDye Terminator Cycle Sequencing Mix (Applied Biosystems). Cycle sequencing reactions were performed in 10 μl reaction volumes containing 1 μl BigDye 2.5, 1.55 μl 5x sequencing buffer (Applied Biosystems), 2 μl of the amplified products (10–40ng), and 2 μl primer. The thermocycling profile was: (1) 35 cycles of 10 seconds at 96°C, 5 seconds at 54°C and 4 minutes at 60°C. The sequencing reactions were purified and denatured and the samples were run on an ABI 3130 Genetic Analyzer. Sequences were edited using MEGA 6 [28] and aligned using ClustalW [29]. The DNA sequences were translated into amino acids residues to test for the absence of premature stop codons or indels using the program MEGA 6.0.

Phylogenetic analysis

The dataset was analysed in the program PartitionFinder [30] to determine the best partitioning scheme and nucleotide substitution models. The optimal partition strategy is shown in S3 Table. The phylogenetic analyses were performed in the programs Garli 2.0 [31], for Maximum Likelyhood (ML) and MrBayes 3.2.6 [32], for Bayesian inference (BI). The support values of the ML analysis were calculated by 1000 bootstrap replications [33]. For BI analysis, four independent Markov Chain Monte Carlo (MCMC) were performed with 30 million generations each, sampling one of every 1,000 trees. The support values of the BI analysis are given by posterior probability. The quality of the MCMC chains was evaluated in Tracer 1.6 [34], and a 25% burn-in was removed.

Divergence-time estimation

The divergence date analyses were performed in Beast v.1.8 [35, 36], using the concatenated dataset and the same partitions as described above, and a normal uncorrelated relaxed clock model, which emphasizes the minimum age and has been considered appropriate for fossil calibration points [37]. Bayesian inference was performed with 30 million generations, a sampling frequency of 1000. The value of parameters of the analysis, convergence of the MCMC chains, sample size and the stationary phase of the chains were evaluated using Tracer v. 1.6 [34]. A Yule speciation process for the tree prior [38] was used, establishing three calibration points. The first one was placed at the stem comprising the genera Aphanius and Valencia. It corresponds to the origin of the crown European cyprinodontoid clade [39], which was estimated to have occurred at least 33 Ma on the basis of the oldest identifiable clade member, the fossil Prolebias stenoura Sauvage, 1874 from the Lower Stampien (Lower Oligocene) of Puy-de-Dôme, France [40] (prior setting: lognormal distribution, mean = 33 and standard deviation = 0.5). The second point was the node where Aphanius and Valencia diverge, corresponding to the most ancient record with recognizable synapomorphies of Aphanius, the fossil Aphanius chebianus (Obrhelová, 1985) of the Ottnangian (Lower Miocene) of the Cheb basin, Czech Republic ([41]; WJEMC, pers. observ.), with 17 Ma (prior setting: lognormal distribution, mean = 17 and standard deviation = 0.5). The third point was placed in the divergence between Anableps and Jenynsia, this position corresponds to the Tucumanableps cionei Sferco, Herbst, Aguilera & Mirande, 2017. This species is recognised as the sister goup of Anabelps from the Meddle Miocene of the Río Sali formation, Argentina [9], with 12 Ma (prior settings: lognormal distribution, mean = 12 and standard deviation = 0.5). The position of each calibration points is also shown in Fig 2.

thumbnail
Fig 2. Time-calibrated analysis performed in Beast 1.8.3.

Black circles indicate the position of the calibrations used. Bars represent the 95% highest posterior credibility intervals of divergence times.

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

Reconstructing ancestral habitat

An analysis of ancestral states was used to infer possible past habitat scenarios without a priori assumptions about states relationships [42]. Reconstruction of ancestral states were conducted in the program RASP 3.02 [43] using the analytical approach Bayesian Binary MCMC (BBM) [44]. The analysis was directed to historical pattern of habitats, based on: (A) freshwater and (B) brackish water.

Results

Phylogenetic relationships

The analyses of BI and ML (Fig 3) presented identical topologies and high supports for most nodes. Both suborders of Cyprinodontiformes, Cyprinodontoidei, and Aplocheiloidei, were recovered with high support values. Three main clades were found in Cyprinodontoidei, here referred as clades A, B and C, clades A and C comprising New World taxa and clade B restricted to Old World taxa (Fig 3). Both analyses did not recover Poeciliidae (sensu [1]) as monophyletic, with their taxa appearing in three unrelated lineages: one, including African procatopodines as part of the Old-World clade B, and two, comprising the subfamily Poeciliinae and the Amazonian genus Fluviphylax as unrelated parts of the New World clade A (Fig 3). The analyses highly supported Anablepidae as sister to the internal fertilizing family Poeciliidae (sensu [45]), as well as a clade containing Anablepidae, Poeciliidae, and Fluviphylax (clade A). In Anablepidae the genus Oxyzygonectes was recovered as the sister group of Anableps and Jenynsia, and both subgenera of Jenynsia were also corroborated [2].

thumbnail
Fig 3. Maximum likelihood (ML) tree.

ML tree with the same topology of the Bayesian Inference (BI) tree. Values on the nodes are relative to the Bootstrap of ML and posterior probability of BI analyses respectively. Asterisks mean maximum values and dashes mean values under 60%.

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

Time-calibrated analysis

The Beast 1.8.3 time-calibrated analysis is shown in Fig 2. Divergence times and their estimated 95% highest posterior density (HPD) intervals place the origin of Cyprinodontiformes in the Upper Cretaceous (77.6 Ma). Subsequent diversification of Cyprinodontoidei in the clades A, B and C occurred in Eocene (between 41.4–37.2 Ma). The currently recognised families originated between Eocene and Miocene (38.1–19.6 Ma).

Reconstructing ancestral habitat

The reconstruction of ancestral habitats recovered the ancestor of Anablepidae as a brackish water species (Fig 4). In the genus Jenynsia, the BBM analysis recovered the most recent common ancestor as a freshwater species, a condition kept in derived lineages, except in the clade comprising J. lineata and J. darwini, wherein the analysis recovered the ancestor inhabiting brackish waters.

thumbnail
Fig 4. Analysis of evolution of ancestral habitat.

BBM analysis based on the Bayesian Inference (BI) tree. Blue indicates freshwater (A), and yellow indicates brackish water (B). Only Anablepidae is shown.

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

Discussion

Phylogenetic relationships of Anablepidae

The present study is congruent with most recent analyses of Cyprinodontiformes (e.g. [25, 39, 45]), where Poeciliidae (sensu [1],) is not corroborated. Similarly, the phylogenetic analyses recovered three unrelated clades, presently receiving family status, Fluviphylacidae, Poeciliidae and Procatopodidae [45] (Fig 3). Once again, the African clade Procatopodidae was recovered as the sister group of a clade comprising Valenciidae and Aphaniidae [39, 45, 46], composing the Clade B, and the New World lineages Poeciliidae and Fluviphylacidae, recovered in the Clade A along with the Anablepidae (Fig 3). The proposed sister group relationships between Anablepidae and Poeciliinae have been reported since the first molecular phylogenetic analyses of Cyprinodontiformes [25, 47, 48], but the position of Fluviphylax as the sister group of this clade was only suggested in recent studies [25, 45].

In addition, this study supports relationships among genera of Anablepidae as proposed in analyses exclusively based on morphological characters [1, 2], with Oxyzygonectes as the sister group of the clade containing Jenynsia and Anableps.

Timing of origin and diversification of the Anablepidae

In contrast to previous phylogenetic studies where the distribution of the Cyprinodontoidei in both the New and Old World was assumed to be result of continental drift between Africa and South America in the Cretaceous [22, 24], the current time-calibrated analysis (see below) indicated that the whole order Cyprinodontiformes had its origin in the Upper Cretaceous (77.6 Ma, 95% HPD 106.6–57.0 Ma) and the diversification of the Cyprinodontoidei occurred in the Eocene (41.4 Ma, 95% HPD 47.0–37.3 Ma), therefore much more recently than the complete separation of Africa and South America (about 100 Ma). The origin of the Cyprinodontiformes clade just before the Cretaceous-Paleogene transition (K–Pg), followed by a great diversification during the Paleogene, is a pattern also found in other animal groups, such as placental mammals, crocodyliforms, birds and some other groups of teleost fishes (e.g. [49, 50, 51, 52, 53]). The K–Pg event has been associated with a mass extinction event that affected great part of Earth's biodiversity, as well as considered as an important event for the subsequent diversification of vertebrates [49, 54, 55, 56]. The results of the present analyses and the biogeographic information previously published for Poeciliidae [25] indicates that the origin of Anablepidae, and the clade composed by both families, have probably occurred in South America. The origin of the Anablepidae is here estimated to have occurred during the Oligocene (29.6 Ma, 95% HPD 34.7–24.5 Ma), followed by a split between the oviparous brackish and marine genus Oxyzygonectes and the clade of viviparous anablepids comprising the brackish and marine genus Anableps and the predominantly freshwater genus Jenynsia (27.9 Ma, 95% HPD 15.0–12.3 Ma), which diverged from each other in the Miocene (12.5 Ma, 95% HPD 28.7–16.2 Ma). The last time estimate is thus temporally congruent with the Miocene marine transgressions in the Amazon basin [57, 58, 59]. The transition of South American marine fish lineages to freshwater environments have been associated with Miocene marine transgressions, followed by species diversification in freshwater habitats, as the Pebas Mega-wetland system in the Amazon river basin and the Paranean Sea in the Paraná river basin [57, 58]. Those systems were complex ecosystems with different kinds of habitats, promoting diversification of vertebrate and invertebrate lineages and adaptation to freshwater environments [60, 61, 62]. More recently an analysis concluded that marine to freshwater transition in South American drums and pufferfishes was temporally coincident with Miocene marine transgressions in the Amazon river basin [63]. Considering the result of the BBM analysis and that species of Anableps and Oxyzygonectes are restricted to coastal marine habitats of Central America and northern South America (Fig 4), we postulate that Jenynsia is a marine-derived lineage, the origin of the genus is associated with a transition of the ancestor from marine to freshwater environments favoured by Miocene marine transgressions.

The present analysis found that the current geographical range of species of Jenynsia, (Fig 1) [2, 21], may be a consequence of an invasion to the freshwater environments and subsequent dispersal through the Pebas Mega-wetland system, combined with an ancient connection between the Amazon and Paraná river basins [59, 64, 65]. Despite some divergences [61, 62, 63, 64, 65, 66], most authors agree that this connection lasted until the Late Miocene (11–10 Ma) [59, 64, 65]. According to the time-calibrated analysis (Fig 2), species diversification of Jenynsia occurred in the Early Pliocene (4.3 Ma, 95% HPD 5.9–3.0 Ma), after separation of the Amazon and Paraná river basins. The most recent ancestor of Jenynsia must have dispersed along this ancient connection and subsequently invading the Paranean Sea, then reaching the current distribution of the genus, which is partially coincident with the geographical range of this marine transgression in the Paraná river basin from Middle to Late Miocene [21, 67]. The diversification of Jenynsia in Early Pliocene (4.6 Ma, 95% HPD 6.5–3.1 Ma) is congruent with the regression of the Paranean Sea [59, 68], so that the current range of the species of Jenynsia would be the result of the decrease of this marine transgression, as previously proposed for some of the species of genus [67]. The present distribution pattern of Jenynsia is also similar to that found for other aquatic organisms which the distribution is attributed to the Paranean Sea, including the freshwater Crustacean families Palaemonidae, Sergestidae and Trichodactylidae [69]. Both known fossil genera of Anablepiadae were found in an area previously occupied by the Paranean Sea [9]. The presence of Tucumanableps, the sister group of Anableps, and Sachajenynsia, the sister group of all genera of Anablepidae except Oxyzygonectes, in this area also corroborates the ancient connections between Amazon and Paraná river basins. This also indicates that the ancestors of Jenynsia and other Anablepidea fossil species inhabited Pebas Mega-wetland and Paranean Sea.The time-calibrated analysis inferred the divergence between A. dowei, endemic to the Pacific coast of Central America, and the clade comprising A. anableps and A. microlepis, two species endemic to northern South American coast, occurred in the Early Pliocene (3.5 Ma, 95% HPD 5.5–2.0 Ma). This age is coincident with the final stages of the rise of the Panamanian Isthmus [70], as previously suggested [2, 5]. However, due to the complex geological history of Central America and the absence of time information, previously authors argued that this disjunctive distribution pattern could also be caused by older events. The present time-calibrated analysis is the first to highly supports that disjunction as a result of the rise of the Panamanian Isthmus.

The only known fossil record of a recent anablepid genus consists of a premaxilla belonging to an unidentified species of Jenynsia from the Late Pleistocene (230–125 Kyr) of Centinela del Mar, Argentina [71]. The only extant anablepid species found in this area is J. lineata, for which an age of 1.4 Ma (95% HPD 2.4–0.6 Ma; Pleistocene) is herein estimated. Considering the premaxilla morphology, the site localization, and the fossil age, we conclude that this fossil probably represents a specimen of J. lineata or an extinct closely related species.

Supporting information

S1 Table. List of analysed species with respective GenBank accession numbers for each analysed gene.

Bold accession numbers are sequences first used in this study.

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

(DOC)

S3 Table. Best evolutive model for each partition found by PartitionFinder.

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

(DOC)

Acknowledgments

We are thankful to L.S. Macedo, F.S. Rangel-Pereira, M.A. Barbosa, J.L.O. Mattos, A.M. Katz, C. Feltrin, and O.C. Simões for help during collecting trips; to J.L.O. Mattos and B.M. Camisão for assistance in analyses; and to M. Loureiro for the loan of specimens. This study was funded by FAPERJ (Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro) and CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico). Collections were made with permits provided by IBAMA (Instituto Brasileiro de Meio Ambiente e dos Recursos Naturais Renováveis) and ICMBio (Instituto Chico Mendes de Conservação da Biodiversidade).

References

  1. 1. Parenti LR. A Phylogenetic and Biogeographic Analysis of Cyprinodontiform fishes (Teleostei, Atherinomorpha). Bull Am Mus Nat Hist.1981; 168: 335–557.
  2. 2. Ghedotti MJ. Phylogeny and classification of the Anablepidae (Teleostei: Cyprinodontiformes. In Malabarba LR, Reis RE, Vari RP, Lecena ZMS, Lucena CAS editors. Phylogeny and classification of Neotropical Fishes. Porto Alegre: Edipucrs: 1998. pp 561–582.
  3. 3. Calviño P, Alonso F. First record of the genus Jenynsia from marine water on the coast of Punta del Este, Maldonado, Uruguay (Cyprinodontiformes: Anablepidae). J Fish Biol. 2016; 88:1236–1240. pmid:26817617
  4. 4. Amorim PF. Jenynsia lineata species complex, revision and new species description (Cyprinodontiformes: Anablepidae). J Fish Biol. 2018; 92: 1312–1332. pmid:29516517
  5. 5. Miller RR. Ecology, habits and relationships of the middle american cuatro ojos, anableps dowi (pisces: Anablepidae). Copeia. 1979; 1979: 82–91.
  6. 6. Murase A, Angulo A, Miyazaki Y, Bussing W, López M. Marine and estuarine fish diversity in the inner Gulf of Nicoya, Pacific coast of Costa Rica, Central America. Check List. 2014; 10: 1401–1413.
  7. 7. White EI. On a fossil cyprinodont from Ecuador. J Nat Hist.1927; 20: 519–522.
  8. 8. Costa WJEM. Redescription and phylogenetic position of the fossil killifish Carrionellus diumortuus White from the Lower Miocene of Ecuador (Teleostei: Cyprinodontiformes). Cybium. 2011; 35: 181–187.
  9. 9. Sferco E, Herbst R, Aguilera G, Mirande JM. The rise of internal fertilization in the Anablepidae (Teleostei, Cyprinodontiformes): two new genera and species from the Miocene of Tucumán, Argentina. Pap Palaeontol.
  10. 10. Garman S. 1895. The cyprinodonts. Mem Mus Comp Zool. 1895; 19: 1–170.
  11. 11. Baughman JL. An early mention of Anableps. Copeia. 1947; 1947: 200.
  12. 12. Linnaeus C. Systema naturae per regna tria naturae, secundum classes, ordines, genera, species, cum characteribus, differentiis, synonymis, locis. Tomus I. 10th ed. Stockholm: Laurentius Salvius; 1758.
  13. 13. Swamynathan SK, Crawford MA, Robison WG, Kanungo J, Piatigorsky J. Adaptive differences in the structure and macromolecular compositions of the air and water corneas of the “four-eyed” fish (Anableps anableps). FASEB J. 2003; 17: 1996–2005. pmid:14597669
  14. 14. Owens GL, Windsor DJ, Mui J, Taylor JS. A fish eye out of water: ten visual opsins in the four-eyed fish, Anableps anableps. PloS one. 2009; pmid:19551143
  15. 15. Owens GL, Rennison DJ, Allison WT, Taylor JS. In the four-eyed fish (Anableps anableps), the regions of the retina exposed to aquatic and aerial light do not express the same set of opsin genes. Biol Lett. 2012; 8: 86–89. pmid:21775314
  16. 16. Turner CL. Adaptations of viviparity in jenynsiid fishes. J Morphol 1940; 67: 291–297.
  17. 17. Turner CL. Follicular pseudoplacenta and gut modifications in anablepid fishes. J Morphol 1940; 67: 91–104.
  18. 18. Bisazza A, Manfredi S, Pilastro A. Sexual competition, coercive mating and mate assessment in the one sided livebearer, Jenynsia multidentata: Are They predictive of sexual dimorphism?. Ethology. 2000; 106: 961–978.
  19. 19. Sandrini JZ, Rola RC, Lopes FM, Buffon HF, Freitas MM, Martins CDMG, et al. Effects of glyphosate on cholinesterase activity of the mussel Perna perna and the fish Danio rerio and Jenynsia multidentata: In vitro studies. Aquat Toxicol. 2013; 130: 171–173. pmid:23411353
  20. 20. Oliveira RDS, Macieira RM, Giarrizzo T. Ontogenetic shifts in fishes between vegetated and unvegetated tidepools: assessing the effect of physical structure on fish habitat selection. J Fish Biol. 2016; 89: 959–976. pmid:27271815
  21. 21. Aguilera G, Mirande JM, Calviño PA, Lobo LF. Jenynsia luxata, a new species from Northwestern Argentina, with additional observations of J. maculata Regan and phylogeny of the genus (Cyprinodontiformes: Anablepidae). Neotrop Ichthyol. 2013; 11: 565–572.
  22. 22. Ghedotti MJ. Phylogenetic analysis and taxonomy of the poecilioid fishes (Teleostei: Cyprinodontiformes). Zool J Linn Soc. 2000; 130: 1–53.
  23. 23. Rosen DE. A vicariance model of Caribbean biogeography. Syst Biol. 1975; 24: 431–464.
  24. 24. Hrbek T, Seckinger J, Meyer A. A phylogenetic and biogeographic perspective on the evolution of poeciliid fishes. Mol Biol Evol. 2007; 43: 986–998. pmid:17185005
  25. 25. Reznick DN, Furness AI, Meredith RW, Springer MS. The origin and biogeographic diversification of fishes in the family Poeciliidae. PloS one. 2017; pmid:28278162
  26. 26. Wainwright PC, Smith WL, Price SA, Tang KL, Sparks JS, Ferry LA, Kuhn KL, Eytan RI, Near TJ. 2012. The evolution of pharyngognathy: a phylogenetic and functional appraisal of the pharyngeal jaw key innovation in labroid fishes and beyond. Syst Biol. 2012; 61: 1001–1027. pmid:22744773
  27. 27. Betancur-R R, Broughton RE, Wiley EO, Carpenter K, López JA, Holcroft NI, et al. The tree of life and a new classification of bony fishes. PLOS Currents Tree of Life. 2013; pmid:23653398
  28. 28. 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
  29. 29. Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, et al. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003; 31: 3497–500. pmid:12824352
  30. 30. Lanfear R, Calcott B, Ho SY, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012; 29: 1695–1701. pmid:22319168
  31. 31. Zwickl, DJ. Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. Ph.D. Thesis, University of Texas. 2006.Avaliable from: https://repositories.lib.utexas.edu/handle/2152/2666
  32. 32. Ronquist F, Teslenko M, Van Der Mark P, Ayres DL, Darling A, Höhna S, Larget B, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012; 61: 539–42. pmid:22357727
  33. 33. Felsenstein J.Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985; 39: 783–791. pmid:28561359
  34. 34. Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6. 2014. Available from: http://beast.bio.ed.ac.uk/Tracer
  35. 35. 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
  36. 36. Heath TA. A hierarchical Bayesian model for calibrating estimates of species divergence times. Syst Biol. 2012; 61: 793–809. pmid:22334343
  37. 37. Parham JF, Donoghue PCJ, Bell CJ, Calway TD, Head JJ, Holroyd PA, et al. Best Practices for Justifying Fossil Calibrations. Syst Biol. 2012; 61: 346–359. pmid:22105867
  38. 38. Gernhard T, Hartmann K, Steel M. Stochastic properties of generalised Yule models, with biodiversity applications. Journal of mathematical biology. 2008; 57: 713–735. pmid:18509650
  39. 39. Pohl M, Milvertz FC, Meyer A, Vences M. Multigene phylogeny of cyprinodontiform fishes suggests continental radiations and a rogue taxon position of Pantanodon. Vertebr Zool. 2015; 65: 37–44.
  40. 40. Costa WJEM. Oligocene killifishes (Teleostei: Cyprinodontiformes) from southern France: relationships, taxonomic position, and evidence of internal fertilization. Vertebr Zool. 2012; 62: 371–386.
  41. 41. Gaudant J. Note complémentaire sur l'ichthyofaune oligocène de Seifhennersdorf (Saxe, Allemagne) et de Varnsdorf, Kundratice, Lbín, Skalice, Knížecí, etc.(Bohême, République tchèque). Ann Naturhist Mus Wien Ser A Mineral Petrogr Geol Palaontol Anthropol Prahis. 2009; 111A: 281–312.
  42. 42. Ronquist F. Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography. Syst Biol. 1997; 46: 195–203.
  43. 43. Yu Y, Harris AJ, Blair C, He X. RASP (Reconstruct Ancestral State in Phylogenies): a tool for historical biogeography. Mol Phylogenet Evol. 2015; 87, 46–49. pmid:25819445
  44. 44. Ronquist F, Huelsenbeck JP. MrBayes3: Bayesian phylogenetic inference undermixed models. Bioinformatics. 2003; 19:1572–1574. pmid:12912839
  45. 45. Bragança PHN, Amorim PF, Costa WJEM. Pantanodontidae (Teleostei, Cyprinodontiformes), the sister group to all other cyprinodontoid killifishes as inferred by molecular data. Zoosystematics Evol. 2018; 94, 137–145. doi.org/10.3897/zse.94.22173.
  46. 46. Freyhof J, Özuluğ M., & Saç G. Neotype designation of Aphanius iconii, first reviser action to stabilise the usage of A. fontinalis and A. meridionalis and comments on the family group names of fishes placed in Cyprinodontidae (Teleostei: Cyprinodontiformes). Zootaxa. 2017; 4294:573–585.
  47. 47. Meyer A, Lydeard C. The evolution of copulatory organs, internal fertilization, placentae, and viviparity in killifishes (Cyprinodontiformes), as inferred from a DNA phylogeny of the tyrosine kinase gene X-src Proc R Soc Lond B Biol Sci; 1993. 254: 153–162. pmid:8290610
  48. 48. Parker A. Combining molecular and morphological data in fishes systematic: Exemples from Cyprinodontiformes. In: Kocher TD, Stepin CA editors. Molecular systematic of fishes. San Diego: Academic Press; 1997. pp 163–188.
  49. 49. Ericson PG, Anderson CL, Britton T, Elzanowski A, Johansson US, Källersjö M, et al. Diversification of Neoaves: integration of molecular sequence data and fossils. Biol Lett. 2006; 2: 543–547. pmid:17148284
  50. 50. Friedman M. 2010. Explosive morphological diversification of spiny-finned teleost fishes in the aftermath of the end-Cretaceous extinction. Proc R Soc Lond B Biol Sci. 2010; pmid:20133356
  51. 51. Near TJ, Dornburg A, Eytan RI, Keck BP, Smith WL, Kuhn KL, et al. Phylogeny and tempo of diversification in the superradiation of spiny-rayed fishes. Proc Natl Acad Sci U S A. 2013; 110: 12738–12743. pmid:23858462
  52. 52. O’leary MA, Bloch JI, Flynn JJ, Gaudin TJ, Giallombardo A, Giannini NP, et al. The placental mammal ancestor and the post-K-Pg radiation of placentals. Science. 2013; 339: 662–667. pmid:23393258
  53. 53. Kellner AW, Pinheiro AE, Campos DA. A new sebecid from the Paleogene of Brazil and the crocodyliform radiation after the K-Pg boundary. PLoS One. 2014; pmid:24454686
  54. 54. Wible JR, Rougier GW, Novacek MJ, Asher RJ. Cretaceous eutherians and Laurasian origin for placental mammals near the K/T boundary. Nature. 2007; 447: 1003–1006. pmid:17581585
  55. 55. Longrich NR, Tokaryk T, Field DJ. Mass extinction of birds at the Cretaceous-Paleogene (K-Pg) boundary. Proc Natl Acad Sci U S A. 2011; 108: 15253–15257. pmid:21914849
  56. 56. Longrich NR, Bhullar BAS, Gauthier JA. Mass extinction of lizards and snakes at the Cretaceous-Paleogene boundary. Proc Natl Acad Sci U S A. 2012; 109: 21396–21401. pmid:23236177
  57. 57. Lovejoy NR, Bermingham E, Martin AP. Marine incursion into South America. Nature. 1998; 396: 421–422.
  58. 58. Lovejoy NR, Albert JS, Crampton WG. 2006. Miocene marine incursions and marine/freshwater transitions: Evidence from Neotropical fishes. J South Am Earth Sci. 2006; 21: 5–13.
  59. 59. Lundberg JG, Marshall LG, Guerrero J, Horton B, Malabarba MCSL, Wesselingh F, The stage for Neotropical fish diversification: a history of tropical South American rivers. In Malabarba LR, Reis RE, Vari RP, Lecena ZMS, Lucena CAS editors. Phylogeny and classification of Neotropical Fishes. Porto Alegre: Edipucrs: 1998. pp 13–48.
  60. 60. Wesselingh FP, Guerrero J, Räsänen ME, Romero-Pittman L, Vonhof HB. 2006. Landscape evolution and depositional processes in the Miocene Pebas lake/wetland system: Evidence from exploratory boreholes in northeastern Peru. Scr Geol. 2006; 133: 323–361.
  61. 61. Wesselingh FP, Räsänen ME, Irion G, Vonhof HB, Kaandorp R, Renema W, et al. Lake Pebas: A palaeoecological reconstruction of a Miocene, long-lived lake complex in western Amazonia. Cainozoic Res. 2002; 1: 35–81.
  62. 62. Bloom DD, Lovejoy NR. 2011. The biogeography of marine incursions in South America. In: Albert JS, Reis RE, editors. Historical biogeography of neotropical freshwater fishes. Berkely and Los Angeles: UC Press: 2011. pp 137–144.
  63. 63. Bloom DD, Lovejoy NR. On the origins of marine derived freshwater fishes in South America. J Biogeogr. 2017; 44: 1927–1938.
  64. 64. Montoya-Burgos JI. Historical biogeography of the catfish genus Hypostomus (Siluriformes: Loricariidae), with implications on the diversification of Neotropical ichthyofauna. Mol Ecol. 2003; 12: 1855–1867. pmid:12803637
  65. 65. Latrubesse EM, da Silva SA, Cozzuol M, Absy ML. Late Miocene continental sedimentation in southwestern Amazonia and its regional significance: Biotic and geological evidence. J South Am Earth Sci. 2007; 23: 61–80.
  66. 66. Sivasundar A, Bermingham E, Orti G. Population structure and biogeography of migratory freshwater fishes Prochilodus (Characiformes) in major South American rivers. Mol Ecol. 2001; 10: 407–417. pmid:11298955
  67. 67. Aguilera G, Mirande JM. A new species of Jenynsia (Cyprinodontiformes: Anablepidae) from northwestern Argentina and its phylogenetic relationships. Zootaxa. 2005; 1096: 29–39.
  68. 68. Marshall LG, Lundberg JG. Technical comment: Miocene deposits in the Amazon foreland basin. Science. 1996; 237: 123–124.
  69. 69. Collins PA, Giri F, Williner V. Biogeography of the freshwater decapods in the La Plata basin, South America. J Crustacean Biol. 2011; 31: 179–191.
  70. 70. Iturralde-Vinent MA, MacPhee RD. Paleogeography of the Caribbean region: implications for Cenozoic biogeography. Bull Am Mus Nat Hist. 1999; 238: 1–95.
  71. 71. Bogan S, de los Reyes ML, Cenizo MM. Primer registro del género Jenynsia Günther, 1866 (Teleostei: Cyprinodontiformes) en el Pleistoceno Medio tardío de la provincia de Buenos Aires (Argentina). Pap Avulsos Zool. 2009; 49: 81–86.