8urn:lsid:arphahub.com:pub:A71ED5FC-60ED-5DA3-AC8E-F6D2BB5B3573urn:lsid:zoobank.org:pub:C8FA3ADA-5585-4F26-9215-A520EE683979Comparative CytogeneticsCCG1993-07711993-078XPensoft Publishers10.3897/CompCytogen.v11i2.1157911579Research ArticleAnimaliaChordataMugilidaePiscesBiogeographyEvolutionary biologyPacificPacific OceanThe Mugilcurema species complex (Pisces, Mugilidae): a new karyotype for the Pacific white mullet mitochondrial lineageNirchioMauro12OliveiraClaudio3Siccha-RamirezZoila R.3SeneViviani F. de3SolaLuciana4MilanaValentina4valentina.milana@uniroma1.itRossiAnna Rita4Escuela de Ciencias Aplicadas del Mar, Núcleo de Nueva Esparta, Universidad de Oriente, Apartado 174, Porlamar, Isla de Margarita, VenezuelaUniversidad de OrienteIsla de MargaritaVenezuelaUniversidad Técnica de Machala, Av. Panamericana km 5½, Via Pasaje, Machala, El Oro, EcuadorUniversidad Técnica de MachalaMachalaEcuadorDepartamento de Morfologia, Instituto de Biociências Universidade Estadual Paulista, 18618-970 Botucatu, São Paulo, BrazilUniversidade Estadual PaulistaBotucatuBrazilDipartimento di Biologia e Biotecnologie “C. Darwin”, Sapienza - Università di Roma, Via Alfonso Borelli 50, 00161, Rome, Italy“La Sapienza” Rome UniversityRomeItaly
2017070420171122252377B5EFFBE-FFBE-440F-2B35-FF9BFFFEFFF56E4E2B0D-DB63-48E5-B9C4-2A53245953105756240301201710022017Mauro Nirchio, Claudio Oliveira, Zoila R. Siccha-Ramirez, Viviani F. de Sene, Luciana Sola, Valentina Milana, Anna Rita RossiThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.http://zoobank.org/6E4E2B0D-DB63-48E5-B9C4-2A5324595310
Recent molecular phylogenetic analyses have shown that the Mugilcurema Valenciennes, 1836 species complex includes M.incilis Hancock, 1830, M.thoburni (Jordan & Starks, 1896) and at least four “M.curema” mitochondrial lineages, considered as cryptic species. The cytogenetic data on some representatives of the species complex have shown a high cytogenetic diversity. This research reports the results of cytogenetic and molecular analyses of white mullet collected in Ecuador. The analyzed specimens were molecularly assigned to the Mugil sp. O, the putative cryptic species present in the Pacific Ocean and showed a 2n = 46 karyotype, which is composed of 2 metacentric and 44 subtelocentric/acrocentric chromosomes. This karyotype is different from the one described for M.incilis (2n = 48) and from those of the two western Atlantic lineages Mugilcurema (2n = 28), and Mugilmargaritae (2n = 24). Data suggest the need for a morphological analysis to assign a species name to this Pacific lineage.
Nirchio M, Oliveira C, Siccha-Ramirez ZR, de Sene VF, Sola L, Milana V, Rossi AR (2017) The Mugil curema species complex (Pisces, Mugilidae): a new karyotype for the Pacific white mullet mitochondrial lineage. Comparative Cytogenetics 11(2): 225–237. https://doi.org/10.3897/CompCytogen.v11i2.11579
Introduction
The family Mugilidae currently comprises 20 genera and 74 species (Eschmeyer and Fong 2016), which are widely distributed in various tropical, subtropical and temperate coastal regions of the world (Thomson 1997, Froese and Pauly 2016). These fishes show highly conserved morphological and anatomical characteristics, which are often associated with wide distribution ranges and, thus, the family has undergone many taxonomic revisions, both at the genus and species levels (Thomson 1997).
In the last decade, molecular phylogenetic and phylogeographic analyses have revealed that the morphological features commonly used to identify species seem to be insufficient, both to describe the great diversity of species within Mugilidae and to infer the phylogenetic relationships among the species (Durand et al. 2012, Durand and Borsa 2015). More specifically, Durand et al. (2012) showed that a proportion of the species with large distribution ranges, such as Mugilcephalus Linnaeus, 1758 and M.curema Valenciennes, 1836 consists of cryptic species. Specifically referring to M.curema, different mtDNA lineages had been previously identified along the American Atlantic coasts by Heras et al. (2006, 2009) and Fraga et al. (2007). Unfortunately these studies did not adopt a uniform nomenclature for the lineages (see Rossi et al. 2016 for a detailed review) and did not cover the entire species range, that includes both the Eastern and Western Atlantic coasts and the Eastern Pacific coast (Froese and Pauly 2016). Durand et al. (2012) and Durand and Borsa (2015) showed that these lineages are part of a Mugilcurema species complex which includes M.incilis Hancock, 1830 and M.thoburni (Jordan & Starks, 1896), and at least four “M.curema” mitochondrial lineages, considered as cryptic species. The first lineage is distributed along the Atlantic coast of the Americas and retains the name of M.curema, as the type locality of the original M.curema is Bahia, Brazil; the second lineage is present along the Atlantic African coasts and is indicated as Mugil sp. M. The third lineage, indicated as Mugil sp. N, is present in Venezuela and has recently been formally described as a new species, named M.margaritae Menezes, Nirchio, Oliveira & Siccha-Ramirez, 2015 (Menezes et al. 2015). The fourth lineage is distributed along the Pacific coast of the Americas, from the USA to Ecuador, and is indicated as Mugil sp. O.
Cytotaxonomy has been proven to be a powerful tool in revealing different lineages/species within Mugilidae. For example, the presence of different cytogenetic features (Nirchio et al. 2003) provided the basis for the identification of an undescribed species, M.rubrioculus Harrison, Nirchio, Oliveira, Ron & Gaviria, 2007 (Harrison et al. 2007), as well as the first hints about the existence of cryptic species among mullets, which, until then, had been reported under the name of M.curema (Nirchio et al. 2005). Although only two of the four mitochondrial lineages of “M.curema” have been cytogenetically investigated to date, they have been found to differ from each other in diploid number and chromosome formula, as well as differing from all the other mugilids investigated to date (see Rossi et al. 2016 for a review). M.curemasensu strictu shows a karyotype composed of 2n = 28 chromosomes (LeGrande and Fitzsimons 1976, Nirchio et al. 2005) and M.margaritae shows a karyotype composed of 2n = 24 chromosomes (Nirchio and Chequea 1998, Nirchio et al. 2005).
This paper reports the cytogenetic analysis of samples of the white mullet collected in Ecuador (Pacific Ocean); according to Durand and Borsa (2015), it is reasonable to assume that they belong to the Mugil sp. O. The aim of the study is to describe the karyotype and the cytogenetic features of this Mugil sp. and to detect specific cytotaxonomic markers that could be useful for its identification. To verify that all the analyzed specimens belong to the Mugil sp. O, sequences of the mitochondrial cytochrome oxidase I (COI) gene were also produced and compared to those previously obtained by Durand et al. (2012) and Durand and Borsa (2015).
Materials and methods
Seventeen juvenile specimens (undetermined sex), morphologically classified as white mullet (Mugilcurema) according to Harrison (1995), were caught by cast net at Puerto Hualtaco, at the border between Ecuador and Perú (3°26’S; 80°13’W), and transported alive to the laboratory. The fishes were sacrificed with an overdose of benzocaine (250 mg/l), following the guidelines of the AVMA (2013). Small pieces of muscle and cephalic kidneys were removed from all specimens, and nine individuals preserved in 70% ethanol were deposited as voucher specimens (Table 1).
Total genomic DNA was extracted from muscle according to Aljanabi and Martínez (1997).
GenBank accession number (A.N.), sampling areas and references of the Mugil sp. COI sequences used in phylogenetic analyses.
Individual (Voucher number)
A.N.
Sampling area
Reference
0102 (UTMACH0102)
KU504271
Ecuador
Present paper
0103 (UTMACH0103)
KU504271
Ecuador
Present paper
104
KU504271
Ecuador
Present paper
105
KU504271
Ecuador
Present paper
119
KU504272
Ecuador
Present paper
120
KU504271
Ecuador
Present paper
121
KU504271
Ecuador
Present paper
122
KU504271
Ecuador
Present paper
123
KU504271
Ecuador
Present paper
124
KU504271
Ecuador
Present paper
76104 (LBP 76104)
KU504271
Ecuador
Present paper
76105 (LBP 76105)
KU504271
Ecuador
Present paper
76107 (LBP 76107)
KU504271
Ecuador
Present paper
76129 (LBP 76129)
KU504271
Ecuador
Present paper
76130 (LBP 76129)
KU504271
Ecuador
Present paper
76131(LBP 76131)
KU504271
Ecuador
Present paper
76132 (LBP 76132)
KU504271
Ecuador
Present paper
415
JQ060604
El Salvador
Durand et al. 2012
426
JQ060600
El Salvador
Durand et al. 2012
429
JQ060601
El Salvador
Durand et al. 2012
430
JQ060602
El Salvador
Durand et al. 2012
432
JQ060603
El Salvador
Durand et al. 2012
293
JQ060573
Western Panama
Durand et al. 2012
294
JQ060574
Western Panama
Durand et al. 2012
413
JQ060592
Perù
Durand et al. 2012
420
JQ060595
Ecuador
Durand et al. 2012
423
JQ060597
Western Mexico
Durand et al. 2012
425
JQ060599
Western Mexico
Durand et al. 2012
406
JQ060588
Western Mexico
Durand et al. 2012
422
JQ060596
Western Mexico
Durand et al. 2012
396
JQ060580
Togo
Durand et al. 2012
397
JQ060581
Togo
Durand et al. 2012
390
JQ060575
Senegal
Durand et al. 2012
391
JQ060576
Senegal
Durand et al. 2012
392
JQ060577
Senegal
Durand et al. 2012
393
JQ060578
Benin
Durand et al. 2012
394
JQ060579
Benin
Durand et al. 2012
399
JQ060582
Venezuela
Durand et al. 2012
400
JQ060583
Venezuela
Durand et al. 2012
401
JQ060584
Venezuela
Durand et al. 2012
403
JQ060585
Venezuela
Durand et al. 2012
414
JQ060593
Venezuela
Durand et al. 2012
408
JQ060590
Brazil
Durand et al. 2012
411
JQ060591
Guadeloupe
Durand et al. 2012
419
JQ060594
Belize
Durand et al. 2012
404
JQ060586
Eastern USA
Durand et al. 2012
407
JQ060589
Eastern Mexico
Durand et al. 2012
417
JQ060605
Uruguay
Durand et al. 2012
418
JQ060606
Uruguay
Durand et al. 2012
405
JQ060587
Honduras
Durand et al. 2012
6435
JX559534
Galapagos Is.
Durand et al. 2012
6445
JX559535
Galapagos Is.
Durand et al. 2012
299
JQ060609
French Guyana
Durand et al. 2012
302
JQ060608
French Guyana
Durand et al. 2012
780
HQ285928
Venezuela
Hett et al. 2011
782
HQ285929
Venezuela
Hett et al. 2011
785
HQ285930
Venezuela
Hett et al. 2011
786
HQ285931
Venezuela
Hett et al. 2011
788
HQ285927
Venezuela
Hett et al. 2011
LBP: Laboratório de Biologia e Genética de Peixes, UNESP, Botucatu (São Paulo State, Brazil); UTMACH: Laboratorio de Acuicultura, Universidad Técnica de Machala, Ecuador.
A 546 base-pair (bp) fragment of the mitochondrial cytochrome oxidase subunit I gene (COI) was amplified by PCR using primers FishF1 and FishR2 (Ward et al. 2005) and the procedures reported in Milana et al. (2011). The obtained sequences were aligned using the program MEGA5 (Tamura et al. 2011) and submitted to the GenBank database (https://www.ncbi.nlm.nih.gov/Genbank) under accession numbers KU504271–KU504272 (see Table 1 for details). BLAST (Basic Local Alignment Search Tool) software was used for similarity searching of the COI sequences in GenBank.
Tree reconstructions were conducted using neighbor-joining (NJ), maximum-likelihood (ML) and Bayesian inference (BI) analyses. The NJ and ML analyses (1000 bootstrap replicates) were performed using MEGA5 and PhyML v3.0 (Guindonet al. 2010), respectively. The Bayesian analyses were carried out as implemented in MrBayes v3.1.2 (Huelsenbeck and Ronquist 2001); two independent runs of four Markov chains each for 106 generations were performed. Modeltest v3.7 (Posada and Crandall 1998) and MrModeltest v2.3 (Nylander 2008) were used to select the evolutionary models for the ML and the BI analyses, respectively, according to the Akaike information criterion. All 37 COI sequences of Mugil sp. obtained by Durand et al. (2012) and Durand and Borsa (2015), and five COI sequences of M.incilis, previously obtained from Venezuelan specimens by our research group (Hett et al. 2011), were included in the phylogenetic analyses (Table 1).
Cell suspensions were obtained from the cephalic kidney, following the procedure reported by Nirchio and Oliveira (2006). Nucleolus organizer regions (NORs) were identified by silver (Ag) nitrate staining (Howell and Black 1980), and C-banding patterns were obtained following the protocol described by Sumner (1972).
Fluorescence in situ hybridization (FISH) was accomplished according to Pinkel et al. (1986). (TTAGGG)n, major (18S rDNA) and minor (5S rDNA) ribosomal probes were amplified by a polymerase chain reaction (PCR) from the genomic DNA of Eigenmannia sp. 2, using primers available from the literature (Ijdo et al. 1991, White et al. 1990, Pendás et al. 1994, respectively). The 18S rDNA sequences were labelled during PCR with Digoxigenin-11-dUTP; the 5S rDNA and (TTAGGG)n probes were labelled with biotin-16-dUTP. The detection of hybridization signals was performed using conjugated avidin-fluorescein (FITC) for the 18S rDNA probe and anti-digoxigenin-rhodamine for the 5S rDNA and (TTAGGG)n probes. The chromosomes were counterstained with 4',6-diamidino-2-phenylindole (DAPI) and propidium iodide.
The mitotic figures were photographed using an Olympus BX61 photomicroscope equipped with the appropriate selective filters for FISH and with a DP70 digital camera. The images were digitally edited with Adobe Photoshop CS6 Extended.
Results
Similarity searching of the obtained COI sequences in the GenBank database, using the BLAST function, provided 99.6–100% similarity with those obtained by Durand et al. (2012) and Durand and Borsa (2015) for the Pacific white mullet, i.e., the Mugil sp. O. These data were confirmed by the phylogenetic tree topology, obtained by NJ, ML and BI analyses (Fig. 1): all the sequences of white mullets from Ecuador collected in this study clustered within the Pacific Mugil sp. O.
0DBAEC91-DDC2-552F-A086-F79BD01BBA1E
Neighbor-joining tree based on COI sequences. At each node, bootstrap values > 70% (NJ and ML) and posterior probabilities > 0.9 (BI) are shown. Stars indicate sequences obtained in this study; the remaining sequences are from Durand et al. (2012), Durand and Borsa (2015) and from Hett et al. (2011) (see Table 1). For each lineage, the karyotype (2n), the fundamental number (FN) and the chromosome formula are indicated. m: metacentric chromosomes; sm: submetacentric chromosomes; st/a: subtelocentric/acrocentric chromosomes.
https://binary.pensoft.net/fig/130645
In all the individuals, the karyotype is composed of 46 chromosomes, 2 metacentric and 44 subtelocentric/acrocentric, with a fundamental number (FN) of 48 (Fig. 2). The metacentric chromosome pair number 1 was clearly identifiable, whereas the homologues belonging to the subtelocentric and acrocentric series could not be unequivocally identified, due to their uniformly decreasing size. The only exception is the acrocentric chromosome pair classified as number 15 because its homologues show a more or less pronounced terminal achromatic region that is positively stained with AgNO3 (Fig. 2, inset).
57E24DBC-86C6-5676-9BB7-C172F63C6DEF
Conventional Giemsa-stained karyotype of the Pacific white mullet. In the inset, the acrocentric chromosome pair n. 15 sequentially Ag-stained; m: metacentric chromosomes; st/a: subtelocentric/acrocentric chromosomes. Scale bar: 10 μm.
https://binary.pensoft.net/fig/130646
C-banding (Fig. 3) revealed the presence of constitutive heterochromatin at the centromeres of most chromosomes and at the telomeres of eight of them. The metacentric chromosome pair number 1 shows C-positive blocks both in the centromeric and in the terminal location; the acrocentric chromosome pair number 15 shows conspicuous heterochromatic blocks in the terminal region.
Dual FISH (Fig. 4a, b) revealed that the 18S rDNA probe yielded two hybridization signals on the same location detected by silver staining on chromosome pair number 15, whereas the 5S rDNA probes hybridized on one smaller medium-sized subtelo/acrocentric chromosome pair (likely number 20) proximal to the centromere.
Mapping of the (TTAGGG)n telomeric repeats showed the presence of positive signals on both telomeres of all chromosomes. No additional, interstitial or centromeric (TTAGGG)n positive signals were detected (Fig. 4c), even on metacentric chromosome pair number 1.
F34E5E7A-CC68-5076-BDD6-95CF323237C9
Somatic C-banded metaphases of the Pacific white mullets. Arrowheads indicate chromosome pair number one. Arrows indicate terminal heterochromatic blocks on chromosome pair 15. Scale bar: 10 μm.
Somatic metaphases of the Pacific white mullet showing positive sites after FISH (a) with 18S rDNA (arrows) and 5S rDNA (asterisks) probes and (c) with telomeric repeats. Arrowheads indicate chromosome pair number one. In (b) enlargement of selected samples of chromosome pairs 15 and 20, after DAPI staining and FISH with rDNA probes, showing 18S (above) and 5S (below) positive sites, respectively. Scale bar: 10 μm.
https://binary.pensoft.net/fig/130648Discussion
Most of the approximately 20 species of Mugilidae cytogenetically investigated so far (see Rossi et al. 2016 for a review) show a conservative 48 uniarmed (subtelo-, acrocentric chromosomes) karyotype, as well as a conserved FN = 48. Even among the 15 cryptic species identified within the Mugilcephalus species complex (Durand and Borsa 2015), the six cytogenetically investigated lineages share not only the chromosome formula but also the major cytogenetic features (see Rossi et al. 2016for references). The only exceptions to this picture of 2n = 48 are represented by Liza (currently Planiliza) abu (Heckel, 1843) (Değer et al. 2013) and the two “M.curema” lineages so far investigated (LeGrande and Fitzsimons 1976, Nirchio and Cequea 1998, Nirchio et al. 2005). Planilizaabu (Değer et al. 2013), which has a limited Asian distribution, shows a karyotype characterized by 2 large metacentric and 46 subtelo-, acrocentric chromosomes, a diploid number of 48 and FN = 50; a pericentric inversion is invoked to interpret the origin of the metacentric chromosome pair. The two cytogenetically known cytotypes of “M.curema”, i.e., M.curemasensu strictu and M.margaritae, instead, are both characterized by FN = 48 and by a massive presence of biarmed chromosomes, likely derived from extensive Robertsonian centric fusions of subtelo- and acrocentric chromosomes. The M.curema karyotype is composed of 20 metacentric, 4 subtelo- and 4 acrocentric chromosomes, while the M.margaritae karyotype is composed of 22 metacentric and 2 submetacentric chromosomes.
The specimens analyzed in this study, molecularly assigned to the Pacific Mugil sp. O (Durand and Borsa 2015), show a still-undescribed karyotype in the family, i.e., a diploid number of 46 chromosomes, two of which are metacentrics and 44 of which are subtelo-acrocentrics. Most of the chromosomes of this karyotype are uniarmed, as in the other species belonging to the genus Mugil, as well as in different genera of Mugilidae (see Rossi et al. 2016). Nonetheless, the presence of biarmed chromosomes, of a reduced diploid number and of a conserved FN, which are shared with M.curemasensu strictu and M.margaritae, suggests that this karyotype originated by a limited number of centric fusions, i.e., only two uniarmed chromosome pairs were involved. These data confirm that a diploid number different from 48 characterizes all, and so far exclusively, the three investigated “M.curema” mitochondrial lineages, with a number of fusions covering the entire range of possibilities (Fig. 1). In M.margaritae, all uniarmed chromosomes underwent fusion, in M.curema most uniarmed chromosomes underwent fusion, and in the Mugil sp. lineage O, only two pairs. At present, it is not possible to discriminate whether the acquisition of the cytogenetic features and/or heterochromatin sequences that promote centric fusions occurred when the ancestor of the M.curema species complex split from the other Mugil species. Thus, it is not possible to ascertain whether these features were lost in some lineages, or, alternatively, were not acquired at that stage so that not all the molecular lineages were involved. In any case, as the lack of additional telomeric sequences is usually interpreted as a stabilizing factor for fusions (Slijepcevic 1998), the absence of telomeric sequences in a pericentromeric or interstitial position in all the cytogenetically studied “M.curema” lineages suggests that Robertsonian fusions are irreversible.
The Mugil sp. O described in this study shows the presence of NORs on a single chromosome pair, as well as minor ribosomal genes carried by a single chromosome pair. These features are common to most of the mugilids, including all the Mugil species (Rossi et al. 2016). Nevertheless, their location appears to be variable in different species/lineages of the M.curema species complex and does not allow any inference on the direction of chromosomal changes within the species complex.
Further analyses are required to draw a comprehensive picture of the chromosomal evolution within the M.curema species complex. Data on the karyotype of M.thoburniand of the white mullet Mugil sp. M from the East Atlantic (Durand and Borsa 2015) are still missing, as well as data on the molecular analysis of the satellite DNA of the whole complex. It is worth noting that in the phylogenetic trees, the node separating “M.curema” lineages from Mugilincilis was not resolved (Fig. 4, Durand et al. 2012, Durand and Borsa 2015), and the latter species shows a karyotype (Hett et al. 2011) that is the closest to the “typical” all uniarmed mullet karyotype from which, presumably, the “M.curema” Robertsonian karyotypes derived. In a very recent paper (Xia et al. 2016) based both on molecular and diagnostic morphological characters, M.incilis appears to be the sister species to “Mugilcurema” lineages, and M.thoburni is external to them. Unfortunately, only two of the “Mugilcurema” lineages (Durand and Borsa 2015) were included in the analysis.
Data, although preliminary, strongly suggest that each of the “Mugilcurema” lineages within the species complex has its own karyotype. This evidence, and the absence of intermediate karyotypes in the geographic area where different lineages/cytotypes are in sympatry, supports Durand and Borsa’s hypothesis (2015) that chromosomal differences probably prevent interbreeding and indicate the actual reproductive isolation of cryptic species. In this context, a morphological analysis is now needed to assign a species name to the here-examined Pacific Mugil sp. O and possibly to the remaining allopatric East Atlantic Mugil sp. M.
Finally, it needs to be verified whether the karyotype observed in the specimens from Ecuador is also shared by specimens belonging to the Mugil sp. O. from other sampling sites along the American Pacific coast. In particular, a karyotypic analysis is needed for the western Mexican coast, because in the phylogenetic trees two individuals from this region are grouped in a subcluster that is highly divergent from the one that includes the remaining Pacific specimens.
Acknowledgments
This work was supported by Centro de Investigación, Universidad Técnica de Machala -UTMACH-, Ecuador (to MN); Sapienza University, Rome, Italy – Progetto di Università 2016 (to LS, VM, AR); Fundação de Amparo à Pesquisa do Estado de São Paulo -FAPESP-, Brazil and Conselho Nacional de DesenvolvimentoCientífico e Tecnológico -CNPq-, Brazil (to CO).
The authors wish to thank Dr Omar Rogerio Sánchez Romero (UTMACH) and Dr Ricardo Britzke (UTMACH) for their assistance in sampling.
ReferencesAljanabiSMMartinezI (1997) Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques.25(22): 4692–4693. https://doi.org/10.1093/nar/25.22.4692AVMA (2013) Guidelines for the euthanasia of animals. https://www.avma.org/KB/Policies/Documents/euthanasia.pdfDeğerDÜnlüEGaffaroğluM (2013) Karyotype of mullet Lizaabu Heckel, 1846 (Pisces: Mugilidae) from the Tigris River, Turkey.29(1): 234–236. https://doi.org/10.1111/j.1439-0426.2012.02028.xDurandJDBorsaP (2015) Mitochondrial phylogeny of grey mullets (Acanthopterygii: Mugilidae) suggests high proportion of cryptic species.338(4): 266–277. https://doi.org/10.1016/j.crvi.2015.01.007DurandJDShenKNChenWJJamandreBWBlelHDiopKNirchioMGarciade León FJWhitfieldAKChangCWBorsaiP (2012) Systematics of the Mugilidae (Teleostei: Mugiliformes): molecular phylogenetic evidence challenges two centuries of morphology-based taxonomy.64(1): 73–92. https://doi.org/10.1016/j.ympev.2012.03.006EschmeyerWNFongJD (2016) Species by family/subfamily. http://researcharchive.calacademy.org/research/ichthyology/catalog/SpeciesByFamily.asp [Accessed 6 March 2016]FragaESchneiderHNirchioMSanta-BrigidaERodrigues-FilhoLFSampaioI (2007) Molecular phylogenetic analyses of mullets (Mugilidae, Mugiliformes) based on two mitochondrial genes.23(5): 598–604. https://doi.org/10.1111/j.1439-0426.2007.00911.xFroeseRPaulyD (2016)GuindonSDufayardJFLefortVAnisimovaMHordijkWGascuelO (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0.59(3): 307–21. https://doi.org/10.1093/sysbio/syq010HarrisonIJ (1995) Mugilidae. In: FischerWKruppFSchneiderWSommerCCarpenterKENiemVH (Eds) Guia FAO para la identificacion de species para los fines de la pesca., 1293–1298. http://www.fao.org/docrep/010/v6250s/v6250s00.htmHarrisonIJNirchioMOliveiraCRonEGaviriaJ (2007) A new species of mullet (Teleostei: Mugilidae) from Venezuela, with a discussion on the taxonomy of Mugilgaimardianus. Journal of Fish Biology 71 (suppl. A): 76–97. https://doi.org/10.1111/j.1095-8649.2007.01520.xHerasSGonzalezCastro MRoldanMI (2006) Mugilcurema in Argentinean waters: combined morphological and molecular approach.261(2): 473–478. https://doi.org/10.1016/j.aquaculture.2006.07.003HerasSRoldánMIGonzalezCastro M (2009) Molecular phylogeny of Mugilidae fishes revised.19(2): 217–231. https://doi.org/10.1007/s11160-008-9100-3HettAKNirchioMOliveiraCSicchaZRRossiARSolaL (2011) Karyotype characterization of Mugilincilis Hancock, 1830 (Mugiliformes: Mugilidae), including a description of an unusual co-localization of major and minor ribosomal genes in the family.9(1): 107–112. https://doi.org/10.1590/S1679-62252011005000005HowellWMBlackDA (1980) Controlled silver staining of nucleolus organizer regions with a protective colloidal developer: a 1-step method.36(8): 1014–1015. https://doi.org/10.1007/BF01953855HuelsenbeckJPRonquistF (2001) MRBAYES: Bayesian inference of phylogenetic trees.17(8): 754–755. https://doi.org/10.1093/bioinformatics/17.8.754IjdoJWWellsRABaldiniAReedersST (1991) Improved telomere detection using a telomere repeat probe (TTAGGG)n generated by PCR. Nucleic Acids Research 19(17): 4780. https://doi.org/10.1093/nar/19.17.4780LeGrandeWHFitzsimonsJM (1976) Karyology of the mullets Mugilcurema and Mugilcephalus (Perciformes: Mugilidae) from Louisiana.1976(2): 388–391. https://doi.org/10.2307/1443973MenezesNNirchioMOliveiraCSiccha-RamirezR (2015) Taxonomic review of the species of Mugil (Teleostei: Perciformes: Mugilidae) from the Atlantic South Caribbean and South America, with integration of the morphological, cytogenetic and molecular data.3918(1): 1–38. https://doi.org/10.11646/zootaxa.3918.1.1MilanaVFusariARossiARSolaL (2011) Molecular and morphological identification of an uncommon centrolophid fish.6(3): 440–445. https://doi.org/10.2478/s11535-011-0016-xNirchioMCequeaH (1998) Karyology of Mugilliza and M.curema from Venezuela.27(1): 45–50.NirchioMOliveiraC (2006) Universidad de Oriente, Nueva Esparta, Porlamar,Venezuela, 203 pp.NirchioMCervigónFPortoJIRPerezJEGomezJAVillalazJ (2003) Karyotype supporting Mugilcurema Valenciennes, 1836 and Mugilgaimardianus Desmarest, 1831 (Mugilidae: Teleostei) as two valid nominal species.67(1): 113–115. https://doi.org/10.3989/scimar.2003.67n1113NirchioMCiprianoRRCestariMMFenocchioAS (2005) Cytogenetical and morphological features reveal significant differences among Venezuelan and Brazilian samples of Mugilcurema.3(1): 99–102. https://doi.org/10.1590/S1679-62252005000100006NylanderJAA (2008) MrModeltest v2.3. Program distributed by the author. Evolutionary Biology Centre, Uppsala University, Sweden.PendásAMMóranPFreijeJPGarcia-VásquezE (1994) Chromosomal mapping and nucleotide sequence of two tandem repeats of Atlantic salmon 5S rDNA.67(1): 31–36. https://doi.org/10.1159/000133792PinkelDStraumeTGrayJW (1986) Cytogenetic analysis using quantitative, high-sensitivity, fluorescence hybridization.83(9): 2934–2938. https://doi.org/10.1073/pnas.83.9.2934PosadaDCrandallKA (1998) MODELTEST: testing the model of DNA substitution.14(9): 817–818. https://doi.org/10.1093/bioinformatics/14.9.817RossiARCrosettiDLiviS (2016) Genetics of Mugilidae. In: CrosettiDBlaberSJM (Eds) Biology, Ecology and Culture of Grey Mullet (Mugilidae)., 349–397. https://doi.org/10.1201/b19927-16SlijepcevicP (1998) Telomeres and mechanisms of Robertsonian fusion.107(2): 136–140. https://doi.org/10.1007/s004120050289SumnerAT (1972) A simple technique for demonstrating centromeric heterocromatin.75(1): 304–306. https://doi.org/10.1016/0014-4827(72)90558-7TamuraKPetersonDPetersonNStecherGNeiMKumarS (2011) MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.28(10): 2731–2739. https://doi.org/10.1093/molbev/msr121ThomsonJM (1997) The Mugilidae of the world.41: 457–562.WardRDZemlakTSInnesBHLastPRHerbertPND (2005) DNA barcoding Australia’s fish species.360(1462): 1847–1857. https://doi.org/10.1098/rstb.2005.1716WhiteTJBrunsTLeeSTaylorJW (1990) Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: InnisMAGelfandDHSninskyJJWhiteTJ (Eds) PCR Protocols: A Guide to Methods and Applications., 315–322. https://doi.org/10.1016/b978-0-12-372180-8.50042-1XiaRDurandJDFuC (2016) Multilocus resolution of Mugilidae phylogeny (Teleostei: Mugiliformes): implications for the family’s taxonomy.96: 161–177. https://doi.org/10.1016/j.ympev.2015.12.010