SciELO - Scientific Electronic Library Online

vol.17 issue2Historical biogeography of fishes from coastal basins of Maranhão State, northeastern BrazilThe bony fishes (Teleostei) caught by industrial trawlers off the Brazilian North coast, with insights into its conservation author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand




Related links


Neotropical Ichthyology

Print version ISSN 1679-6225On-line version ISSN 1982-0224

Neotrop. ichthyol. vol.17 no.2 Maringá  2019  Epub June 19, 2019 

Original article

Species or population? Systematic status of Vieja coatlicue (Teleostei: Cichlidae)

1Field Museum of Natural History, 1400 S. Lake Shore Dr., 60605 Chicago, IL, USA.,, (corresponding author)

2Universidad de Ciencias y Artes de Chiapas, Museo de Zoología, Colección de Ictiología, Libramiento Norte Poniente 1150, Col. Lajas Maciel, C.P. 29039 Tuxtla-Gutiérrez, CH, Mexico.,

3LSU Museum of Natural Science, Department of Biological Sciences, Louisiana State University, 119 Foster Hall, 70803 Baton Rouge, LA, USA.,

4Department of Biological Sciences, Southeastern Louisiana University, SLU 10736, 70402 Hammond, LA, USA., 0000-0003-1289-9351


Pacific and Atlantic-slope rivers of the Isthmus of Tehuantepec in Mexico contain two cichlid species of the genus Vieja relevant to the current study, V. zonata and V. coatlicue, respectively. The Atlantic-slope species was initially considered a population of uncertain taxonomic status and recently described as a distinct species based on three putatively diagnostic morphological characters. The objective of this study was to combine new and existing genetic data, along with reanalysis of morphological characters, to test the distinctiveness of V. coatlicue. Little genetic variability and no geographic structuring were recovered. Additionally, putatively diagnostic morphological character states were present across both species, failing to separate the forms. The synthesis of molecular and morphological data supports the recognition of V. coatlicue as a junior synonym of V. zonata.

Keywords: Cichlid; Freshwater; Mexico; River; Taxonomy


Dos especies del género Vieja se encuentran distribuidas a lo largo Istmo de Tehuantepec en México; V. zonata se distribuye en los ríos de la vertiente del Pacífico y V. coatlicue en los ríos de la vertiente del Atlántico. La especie que se distribuye en la vertiente del Atlántico era considerada como una población de estado taxonómico incierto, pero recientemente fue descrita como nueva especie en base a tres caracteres morfológicos diagnósticos. El objetivo de este estudio fue utilizar datos moleculares nuevos con datos moleculares previamente publicados en combinación con una reanálisis de los caracteres morfológicos para testear la validez taxonómica de la especie V. coatlicue. Nuestros análisis moleculares no recobran estructura geográfica y además muestran baja variabilidad genética. Además, los tres caracteres morfológicos diagnósticos se encontraron presentes en individuos de ambas especies, nuestro trabajo sugiere que estos caracteres morfológicos diagnósticos no proven información que permita separar a las formas que se distribuyen en la vertiente del Pacífico y la vertiente del Atlántico. Nuestra síntesis de datos moleculares y morfológicos provee evidencia para reconocer a V. coatlicue como una sinónima mas reciente de V. zonata.

Palabras clave: Agua dulce; Cíclido; México; Río; Taxonomía


The genus Vieja Fernández-Yépez, 1969 is a clade of predominantly Northern Middle American cichlids, with highest diversity and endemism occurring in southern Mexico and Guatemala. One species in the genus, Vieja zonata (Meek, 1905), is endemic to Mexico and found in rivers along the Isthmus of Tehuantepec from the Río Tequisistlán in Oaxaca west of Tehuantepec to Río Tapantepec near the Chiapas border likely near the town of San Pedro Tapantepec (Kullander, 2003; Miller et al., 2005; Matamoros et al., 2015). In addition, an Atlantic slope population of this species exists in the Río Coatzacoalcos drainage (Conkel, 1993; Miller et al., 2005) that Miller et al. (2005) claimed did not represent V. zonata but a “population of undetermined taxonomic status”. Del Moral-Flores et al. (2018) recognized the Río Coatzacoalcos population as a distinct species, describing it as V. coatlicue based on three diagnostic morphological characters. Compared to V. zonata, the species was diagnosed based on an elongate body, dorsal-fin origin equal or posterior to opercular border, and a separation between the ends of non-filamentous rays of the pelvic fins and the anus.

In a molecular systematic study of Vieja, McMahan et al. (2010) recovered low genetic differentiation among individuals of sister species V. zonata and V. guttulata (Günther, 1864) in their phylogeny of the genus. Species-level analyses were beyond the scope of that study; however, Atlantic slope (Río Coatzacoalcos) sequences of V. zonata (currently V. coatlicue) were included in their molecular phylogeny and all populations were recovered in the same clade. Given the lack of genetic distinctiveness of the population in Río Coatzacoalcos, the objective of this study was to synthesize genetic data and reassess the morphological characters provided by del Moral-Flores et al. (2018) to test morphological distinctiveness of V. coatlicue in the Río Coatzacoalcos drainage.

Material and Methods

Molecular data. Sequence data for V. coatlicue and several populations of V. zonata have been published in previous studies (e.g. McMahan et al., 2010; Říčan et al., 2016). We sequenced three additional individuals to provide data from other localities for both species (Fig. 1; Tab. 1). Whole genomic DNA was extracted from tissue samples (muscle or fin clip) using the DNeasy kit (Qiagen, Inc.). The mitochondrial cytochrome b (cyt b) gene was amplified and sequenced for samples following protocols described in Perdices et al. (2002). Cytochrome b is an appropriate marker to address the main objective of this study as it has been quite useful to provide resolution at intra- and interspecific levels within Middle American cichlids (López-Fernández et al., 2010; McMahan et al., 2010, 2017a; Harrison et al., 2014; Barrientos-Villalobos et al., 2018). Additionally, various nuclear markers have failed to delineate taxa at this level (López-Fernández et al., 2010; McMahan et al., 2010; Říčan et al., 2016). Amplification products were visualized on a 1% agarose gel, followed by purification with ExoSap-IT, and then sequenced on an ABI 3730 in the Pritzker Laboratory for Molecular Systematics and Evolution at the Field Museum of Natural History. An individual of V. guttulata was sequenced as an additional ingroup taxon to test monophyly of V. zonata and V. coatlicue. Vieja fenestrata (Günther, 1860) was used as an outgroup following phylogenetic relationships of McMahan et al. (2010) and Říčan et al. (2016). All new sequences were deposited in GenBank (Tab. 1).

Fig. 1 Map of Mexico displaying localities for specimens (triangles) examined and tissue samples (circles) for Vieja coatlicue (green) and Vieja zonata (black) used in this study. Star indicates locality for Vieja guttulata sample. 

Tab. 1 Specimens/tissues examined. No. refers to number by tip in phylogeny (Fig. 2). Acronyms follow Sabaj (2016). 

Species No. Tissue No. Voucher No. GenBank Drainage/Locality Latitude Longitude
Vieja coatlicue 2 - - KU854739 Río Coatzacoalcos 18.152194 -94.790194
Vieja coatlicue 5 - - KU854741 Río Coatzacoalcos 16.767778 -95.021806
Vieja coatlicue 7 SLU-TC 347 SLU 5010 FJ668647 Río Coatzacoalcos 16.77148 -95.01855
Vieja coatlicue 10 SLU-TC 341 SLU 5010 FJ668639 Río Coatzacoalcos 16.77148 -95.01855
Vieja coatlicue 11 SLU-TC 348 SLU 5010 FJ668638 Río Coatzacoalcos 16.77148 -95.01855
Vieja coatlicue 13 - TNHC 29065 AY324023 Río Coatzacoalcos; Río Dos Caños 18.168106 -94.916667
Vieja coatlicue 15 - TNHC 29054 AY324024 Río Coatzacoalcos; Río Sarabia 17.102197 -95.040522
Vieja coatlicue 16 - - KU854740 Río Coatzacoalcos 16.767778 -95.021806
Vieja coatlicue 17 - MNCN 147114 DQ990728 Río Coatzacoalcos; Río Sarabia 17.102197 -95.040522
Vieja coatlicue - UNICACH 6702 - Río Coatzacoalcos; Río Jaltepec 17.387582 -95.057214
Vieja coatlicue - UMMZ 184759 - Río Coatzacoalcos; Río Almoloya 16.8333333 -95
Vieja coatlicue - UMMZ 178548 - Río Coatzacoalcos; Donaji 17.386577 -95.05784
Vieja fenestrata - - DQ990726 Río Coatzacoalcos - -
Vieja guttulata SLU-TC 2272 SLU 8117 MK837107 Río Coatán; Tapachula 14.896442 -92.323668
Vieja zonata 1 - TNHC 29059 AY324026 Río los Perros 16.557617 -95.091117
Vieja zonata 3 - - KU854738 Río Tehuantepec 16.339361 -95.241917
Vieja zonata 4 - - MK837015 Río los Perros 16.442715 -95.028022
Vieja zonata 6 SLU-TC 343 SLU 5014 FJ668637 Río Tehuantepec 16.377621 -95.241244
Vieja zonata 8 SLU-TC 344 SLU 5014 FJ668646 Río Tehuantepec 16.377621 -95.241244
Vieja zonata 9 SLU-TC 342 SLU 5014 FJ668642 Río Tehuantepec 16.377621 -95.241244
Vieja zonata 12 - - MK837016 Río Ostuta 16.499821 -94.432838
Vieja zonata 14 - TNHC 29058 AY324025 Río los Perros 16.557617 -95.091117
Vieja zonata - UMMZ 178573 - Río Tehuantepec 16.415523 -95.597819
Vieja zonata - UMMZ 184746 - Río Novillero, Tapantepec 16.390642 -94.133723
Vieja zonata - UNICACH 6356 - Río Ostuta 16.43171 -94.24725
Vieja zonata - FMNH 3776 - Niltepec - -
Vieja zonata - FMNH 3777 - Niltepec - -

Sequence data were visually inspected and aligned using the software Geneious v.9.1 (Kearse et al., 2012). A multiple sequence alignment was generated using the Muscle algorithm (Edgar, 2004) implemented in Geneious using default parameters. The final alignment consisted of 1,050 base pairs. Partition scheme and model selection were performed using the software PartionFinder (Lanfear et al., 2012) using a greedy search. The best partition scheme and substitution model were selected using the Akaike Information Criterion (AIC).

Phylogenetic hypotheses were inferred under a Bayesian framework using the software Mr. Bayes v.3.2.6 (Hulsenbeck et al., 2001) implemented in the CIPRES portal (Miller et al., 2010). Three independent analyses of 12,000,000 generations of two replicated runs each were executed, with tree sampling every 6,000 generations. Sampling stationarity was assessed via visual inspection in Tracer v.1.5 (Rambaut et al., 2018) and observed average standard deviations of split frequencies (<0.01). The initial 25% of sampled trees were discarded as burn-in. Bayesian posterior probabilities were calculated using the post burn-in trees, and the three independent runs were compared to evaluate topological congruence.

Due to missing sequence data in one sample (KU854740), estimation of genetic diversity and reconstruction of a haplotype network relied on a slightly reduced dataset (1,011 base pairs). Nucleotide (π; Nei, Li, 1979), genetic (θw; Watterson, 1975), and haplotype (Hd) diversity indices were calculated in the software DnaSP v.5.10 (Librado, Rozas, 2009). Matrilineal haplotypes were reconstructed using statistical parsimony analysis in the software TCS (Templeton et al., 1992) with a 95% confidence limit for connected haplotypes. The recovered network with TCS was used to infer gene genealogies using a median-joining network (Bandelt et al., 1999) implemented in the software POPART (Leigh, Bryant, 2015).

Morphological data. Material examined covered four localities for V. zonata (N= 29) along Pacific drainages in the Isthmus of Tehuantepec, including the holotype and paratype, and three populations of V. coatlicue (N= 22) in the Río Coatzacoalcos (Tab. 1; Fig. 1). Our goal was to reassess the putative morphological distinctiveness of V. coatlicue relative to V. zonata; thus, characters examined herein focused on the three described by del Moral-Flores et al. (2018) to diagnose V. coatlicue as a distinct species from V. zonata. Standard length (SL) and maximum body depth (BD) were measured for each specimen. The ratio of SL:BD, also known as the fineness ratio, was calculated and used as a measure of body elongation, with greater values indicating more elongate bodies (Aguirre et al., 2016). Additionally, data on anterior insertion of the dorsal fin relative to posterior of the operculum (anterior/posterior to or in-line with) and presence/absence of a space between non-filamentous portions of the pelvic rays and the anus were recorded. These data were assessed to determine if the proposed diagnostic characters could separate the two species and if morphological differences were congruent with any observed patterns of genetic differentiation.


Molecular data. A total of 10 individuals of V. coatlicue (four localities) and seven individuals of V. zonata (five localities) were analyzed based on new and previously published sequences. The best partition scheme was by codon with independent substitution models (1st codon: HKY, 2nd codon: F81, 3rd codon: GTR+gamma). A single clade containing V. zonata and V. coatlicue was recovered based on phylogenetic analysis of cyt b data, with an average 0.2% sequence divergence observed across all individuals of both species. We recovered no geographic structuring within this clade (Fig. 2). Vieja guttulata was recovered as sister to the clade containing V. zonata and V. coatlicue, with an average sequence divergence of 2.2%. This result is consistent with McMahan et al. (2015) and Říčan et al. (2016). Sequence divergence between the outgroup (V. fenestrata) and ingroup was 2.7%.

Fig. 2 Bayesian phylogeny based on cytochrome b sequence data. Number above clades is posterior probability. 

We recovered four haplotypes in the network generated by statistical parsimony (Fig. 3). One haplotype was common and widely distributed across the Coatzacoalcos and Tehuantepec river basins and shared between individuals of V. coatlicue and V. zonata (Fig. 3). Finally, all samples from the clade containing V. coatlicue and V. zonata showed low levels of nucleotide (π=0.493, standard deviation 0.131), genetic (θw =0.0016, standard deviation 0.0078), and haplotype (Hd=0.0031, standard deviation 0.00139) diversity across populations and river basins.

Fig. 3 Haplotype network recovered from analysis of the cyt b dataset for populations of Vieja coatlicue (Atlantic) and Vieja zonata (Pacific). Dashes lines represent mutational steps. 

Morphological data. The average fineness ratio for V. coatlicue was 2.11 (N=22, standard deviation 0.11, minimum 1.87, maximum 2.31), compared to 2.28 (N=29, standard deviation 0.1, minimum 2.04, maximum 2.45) for V. zonata (Tab. 2). This reveals a slightly more streamlined or elongate overall body shape in V. zonata, inconsistent with the purported elongate body diagnostic of V. coatlicue. Across both species we found that most individuals possessed a dorsal-fin origin marginally anterior to- or in-line with the posterior margin of the operculum, although a moderately higher number of individuals of V. zonata exhibited a dorsal-fin origin slightly posterior to the operculum (Tab. 3; Fig. 4). Dorsal-fin origin was only vaguely anterior to the posterior end of the operculum in individuals with this character state and could easily be considered in-line with posterior of the operculum. Additionally, we observed these character states to be skewed in bent specimens, depending on if gathering data from the right or left sides. Absence of a space between non-filamentous rays of the pelvic fin and the anus was observed across the majority of individuals of both V. coatlicue and V. zonata (Tab. 3; Fig. 5).

Fig. 4 Dorsal-fin origin relative to operculum in Vieja coatlicue (top; UMMZ 178548) and Vieja zonata (bottom; FMNH 3776 (holotype), UMMZ 184746), illustrating variability within and between species. Scale bar equals 1cm. 

Fig. 5 Pelvic fin relative to anus in Vieja coatlicue (left of vertical line; UMMZ 178548) and Vieja zonata (right of vertical line; FMNH 3776 (holotype), UMMZ 184746), illustrating variability within and between species. Asterisk indicates holotype of V. zonata

Tab. 2 Fineness ratios (SL:BD) for Vieja coatlicue and Vieja zonata. SD = standard deviation. 

N Mean SD Minimum Maximum
V. coatlicue 22 2.11 0.11 1.87 2.31
V. zonata 29 2.28 0.10 2.04 2.45

Tab. 3 Number of individuals of Vieja coatlicue and Vieja zonata with examined character states dorsal-fin origin relative to operculum and presence/absence of space between non-filamentous ends of the pelvic fins and anus. 

Dorsal-fin origin relative to operculum
anterior equal posterior
V. coatlicue 10 8 4
V. zonata 9 9 11
Separation between pelvic fin and anus
absent present
V. coatlicue 18 4
V. zonata 25 4


We found a lack of monophyly for V. zonata and V. coatlicue, as well as absence of any reasonable level of genetic divergence across species or populations. Additionally, reanalysis of putatively diagnostic morphological characters fails to unambiguously separate the two species, with all character states being observed across the range of the two species. Fineness ratios show the opposite trend from description of V. coatlicue, with slightly more elongate bodies in populations of V. zonata. The origin of the dorsal fin relative to posterior of the operculum exhibits variability between the two species; however, in nearly all specimens the origin is very close to in-line with the posterior edge of the operculum. Few individuals, across both species, exhibit a clearly anterior or posterior condition. We found that specimens of either species rarely possess a clearly defined space or gap between tips of the (non-filamentous) pelvic fin and the anus. High levels of variability were consistently observed for this character across all material examined. Indeed, the holotype of V. zonata possesses the most evident space between the non-filamentous pelvic rays and the anus (Fig. 5). The variability in these morphological characters is high and consistent across populations of both species. Thus, the combination of morphological and genetic data supports the recognition of V. coatlicue as a junior synonym of V. zonata.

Many studies have demonstrated high levels of intraspecific morphological variability within cichlid fishes (McMahan et al., 2017b; Barrientos-Villalobos et al., 2018; Gómez-González et al., 2018). While genetic data are certainly not a required component of detection and description of new species, these independent sources of data can offer additional ways to test hypotheses regarding species delineation. While McMahan et al. (2010) provided sequence data from populations of V. zonata from the Río Tehuantepec and Río Coatzacoalcos, admittedly the authors did not discuss this in text given it was outside the scope of that study. However, locality data and sequences for these specimens were published and deposited in GenBank. Říčan et al. (2016) later found the same lack of genetic distinctiveness between V. zonata from both slopes of the Isthmus of Tehuantepec. Anzueto-Calvo et al. (2016) additionally report the occurrence of V. zonata in the Río Grijalva, which should be studied in future analyses of this and related species.

The Isthmus of Tehuantepec represents the narrowest geographic point between the Pacific and Atlantic basins in the northern portion of Middle America and has been demonstrated to be an important region for organismal dispersal across Mexico, permitting taxa to expand between rivers along both slopes (Quiroz-Martínez et al., 2014). This hypothesis is additionally supported by occurrence of a population of V. zonata in the Atlantic Río Coatzacoalcos. Thus, V. zonata in the Atlantic represent a recently dispersed population and not a distinct lineage.


We are thankful to Susan Mochel (FMNH) and Kevin Swagel (FMNH) for invaluable help with logistics and specimen loans. Doug Nelson (UMMZ) and Hernán López-Fernández (UMMZ) kindly provided loans of specimens supporting this work. Many thanks to Adán Gómez González for valuable discussions and fieldwork assistance. This work was supported, in part, with a grant from the Negaunee Foundation to CDM. Molecular work was carried out in the Pritzker Laboratory for Molecular Systematics and Evolution at the Field Museum of Natural History, with support from the Pritzker Foundation.


Aguirre WE, Navarrete R, Malato G, Calle P, Loh MK, Vital WF, Valadez G, Vu V, Shervette VR, Granda JC. Body shape variation and population genetic structure of Rhoadsia altipinna (Characidae: Rhoadsiinae) in Southwestern Ecuador. Copeia. 2016; 104(2):554-69. [ Links ]

Anzueto-Calvo MJ, Velázquez-Velázquez E, Gómez-González AE. Peces de la Reserva de la Biosfera Selva El Ocote y presa Nezahualcóyotl (Malpaso) Chiapas, México. Rev Mex Biodivers. 2016; 87(3):972-79. [ Links ]

Bandelt HJ, Foster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999; 16(1):37-48. [ Links ]

Barrientos-Villalobos J, Schmitter-Soto JJ, Espinosa de los Monteros AJ. Several subspecies or phenotypic plasticity? A geometric morphometric and molecular analysis of variability of the Mayan cichlid Mayaheros urophthalmus in the Yucatan. Copeia . 2018; 106(2):268-78. [ Links ]

Conkel D. Cichlids of North and Central America. New Jersey: T.F.H. Publications; 1993. [ Links ]

Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32(5):1792-97. [ Links ]

Gómez-González AE, Alvarez F, Matamoros WA, Velázquez-Velázquez E, Schmitter-Soto JJ, Gonzalez-Diaz AA, McMahan CD. Redescription of Vieja hartwegi (Taylor & Miller 1980) (Teleostei: Cichlidae) from the Grijalva River basin, Mexico and Guatemala, with description of a rheophilic morph. Zootaxa. 2018; 4375(3):371-91. [ Links ]

Harrison E, Trexler JC, Collins TM, Vasquez-Domínguez E, Razo-Mendivil U, Matamoros WA, Barrientos C. Genetic evidence for multiple sources of the non-native fish Cichlasoma urophthalmus (Günther; Mayan cichlids) in Southern Florida. PLoS ONE. 2014; 9(9):e104173. [ Links ]

Hulsenbeck JP, Ronquist F, Nielsen R, Bollback JP. Bayesian inference of phylogeny and its impact on evolutionary biology. Science. 2001; 294(5550):2310-14. [ Links ]

Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Meintjes P, Drummond A. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012; 28(12):1647-49. [ Links ]

Kullander SO. Family Cichlidae (Cichlids). In: Reis RE, Kullander SO, Ferraris CJ Jr., organizers. Check list of the freshwater fishes of South and Central America. Porto Alegre: Edipucrs; 2003. p.605-54. [ Links ]

Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol . 2012; 29(6):1695-701. [ Links ]

Leigh JW, Bryant D. PopART: Full-feature software for haplotype network construction. Methods Ecol Evol. 2015; 6(9):1110-16. [ Links ]

Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics . 2009; 25(11):1451-52. [ Links ]

López-Fernández H, Winemiller KO, Honeycutt RL. Multilocus phylogeny and rapid radiations in Neotropical cichlid fishes (Perciformes: Cichlidae: Cichlinae). Mol Phylogenet Evol. 2010; 55(3):1070-86. [ Links ]

Matamoros WA, McMahan CD, Chakrabarty P, Albert JS, Schaefer JF. Derivation of the freshwater fish fauna of Central America revisited: Myers’s hypothesis in the twenty-first century. Cladistics. 2015; 31(2):177-88. [ Links ]

McMahan CD, Geheber AD, Piller KR. Molecular systematics of the enigmatic Middle American genus Vieja (Teleostei: Cichlidae). Mol Phylogenet Evol . 2010; 57(3):1293-300. [ Links ]

McMahan CD, Matamoros WA, Piller KR, Chakrabarty P. Taxonomy and systematics of the herichthyins (Cichlidae: Tribe Heroini), with the description of eight new Middle American Genera. Zootaxa . 2015; 3999(2):211-34. [ Links ]

McMahan CD, Ginger L, Cage M, David KT, Chakrabarty P, Johnston M, Matamoros WA. Pleistocene to Holocene expansion of the black-belt cichlid in Central America, Vieja maculicauda (Teleostei: Cichlidae). PLoS ONE . 2017a; 12(5):e0178439. [ Links ]

McMahan CD, Kutz J, Murray CM, Chakrabarty P, Geheber AD, Elias D. Objectively measuring subjectively described traits: geographic variation in body shape and caudal coloration pattern within Vieja melanura (Teleostei: Cichlidae). Rev Biol Trop. 2017b; 65(2):623-31. [ Links ]

Miller RR, Minckley WL, Norris SM. Freshwater fishes of México. Chicago (IL): University of Chicago Press; 2005. [ Links ]

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), 2010; New Orleans, LA. 2010. p. 1-8. [ Links ]

del Moral-Flores LF, López-Segovia E, Hernández-Arellano T. Vieja coatlicue sp. nov., una nueva especie de cíclido (Actinopterygii: Cichlidae) de la Cuenca del Río Coatzacoalcos, México. Rev Zool, Universidad Nacional Autónoma de México. 2018; 29:1-12. [ Links ]

Nei M, Li WH. Mathematical model for studying genetic variation in terms of restriction endonucleases. PNAS. 1979; 76(10):5269-73. [ Links ]

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

Quiroz-Martínez B, Álvarez F, Espinosa H, Salgado-Maldonado G. Concordant biogeographic patterns among multiple taxonomic groups in the Mexican freshwater biota. PLoS ONE . 2014; 9(8):e105510. [ Links ]

Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior summarization in Bayesian phylogenetics using Tracer 1.7. Syst Biol. 2018; 67(5):901-04. [ Links ]

Říčan O, Piálek L, Dragová K, Novák J. Diversity and evolution of the Middle American cichlid fishes (Teleostei: Cichlidae) with revised classification. Vertebr Zool. 2016; 66(1):1-102. [ Links ]

Sabaj MH, editor. Standard symbolic codes for institutional resource collections in herpetology and ichthyology: an online reference. Version 6.5 [Internet]. American Society of Ichthyologists and Herpetologists; 2018 [updated 16 Aug 2016]. Available from: ]

Templeton AR, Crandall KA, Sing CF. A cladistics analysis of phenotypic association with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics. 1992; 132(2):619-33. [ Links ]

Watterson GA. On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975; 7(2):256-76. [ Links ]

Received: January 22, 2019; Accepted: April 15, 2019

Edited by

Hernán López-Fernández

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License