Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 05 February 2018
Sec. Evolutionary and Population Genetics
This article is part of the Research Topic Genetics, Evolution, and Conservation of Neotropical Fishes View all 34 articles

Identification of a New Mullet Species Complex Based on an Integrative Molecular and Cytogenetic Investigation of Mugil hospes (Mugilidae: Mugiliformes)

  • 1Facultad de Ciencias Agropecuarias, Universidad Técnica de Machala, Machala, Ecuador
  • 2Departamento de Morfologia, Instituto de Biociências, Universidade Estadual Paulista “Júlio de Mesquita Filho”, São Paulo, Brazil
  • 3Dipartimento di Biologia e Biotecnologie “C. Darwin”, Sapienza Università di Roma, Rome, Italy

Mullets are very common fishes included in the family Mugilidae, (Mugiliformes), which are characterized by both a remarkably uniform external morphology and internal anatomy. Recently, within this family, different species complexes were molecularly identified within Mugil, a genus which is characterized by lineages that sometimes show very different karyotypes. Here we report the results of cytogenetic and molecular analyses conducted on Mugil hospes, commonly known as the hospe mullet, from Ecuador. The study aims to verify whether the original described species from the Pacific Ocean corresponds to that identified in the Atlantic Ocean, and to identify species-specific chromosome markers that can add new comparative data about Mugilidae karyotype evolution. The karyotype of M. hospes from Ecuador is composed of 48 acrocentric chromosomes and shows two active nucleolar organizer regions (NORs). In situ hybridization, using different types of repetitive sequences (rDNAs, U1 snDNA, telomeric repeats) as probes, identified species-specific chromosome markers that have been compared with those of other species of the genus Mugil. Cytochrome c oxidase subunit I (COI) sequence analysis shows only 92–93% similarity with sequences previously deposited under this species name in GenBank, all of which were from the Atlantic Ocean. Phylogenetic reconstructions indicate the presence of three well-supported hospe mullet lineages whose molecular divergence is compatible with the presence of distinct species. Indeed, the first lineage includes samples from Ecuador, whereas the other two lineages include the Atlantic samples and correspond to M. brevirostris from Brazil and Mugil sp. R from Belize/Venezuela. Results here provided reiterate the pivotal importance of an integrative molecular and cytogenetic approach in the reconstruction of the relationships within Mugilidae.

Introduction

Mullets is the popular name of fishes included in the Mugilidae, a species rich family that is the only representative of the order Mugiliformes. These fishes are distributed in several coastal aquatic habitats in tropical, subtropical and temperate regions of the world, where they are ecologically, recreationally and commercially important (Thomson, 1966). According to different authors (see González-Castro and Ghasemzadeh, 2016 and references herein), the family has approximately 26 genera, but Eschmeyer and Fong (2017) ascribe to Mugilidae 20 genera and 75 valid species.

In Mugilidae, most of the classical morphological characters used in species identification and/or systematics have poor diagnostic power and morphometric variability is limited (Schultz, 1946; Thomson, 1997; González-Castro, 2007; González-Castro and Ghasemzadeh, 2016). These characteristics are associated with the wide distribution of most of the species, which raises questions about their actual taxonomic status. Cytogenetic and molecular studies have provided important data for understanding the systematic relationships and evolutionary pathways among mullet species (Harrison et al., 2007; Sola et al., 2007; Durand et al., 2012; Durand and Borsa, 2015). These studies have also shown that it is necessary to use integrative approaches to study mugilids. Indeed, the use of repetitive sequences such as ribosomal genes (18S rDNA and 5S rDNA) as probes in FISH mapping has been shown to be a very informative cytotaxonomic tool in revealing different lineages/species within Mugilidae (Nirchio et al., 2007, 2017; Sola et al., 2007). On the other hand, the utility of molecular markers in this family to identify the species, better define the genera, and reconstruct their phylogenetic relationships, is well-represented by the huge amount of literature on this topic published in the last 15 years (see Rossi et al., 2016 for a review and Durand et al., 2017). In addition, molecular phylogenetic analyses have been used successfully in the investigation of chromosome evolution in some fish groups as those of the genus Characidium (Pansonato-Alves et al., 2014) and Triportheus (Yano et al., 2014), and in Geophagus brasiliensis (Alves-Silva and Dergam, 2015).

Mugil, which presently includes 16 valid species (Eschmeyer and Fong, 2017), is the most cytogenetically studied genus among the Mugilidae. Nine species have been investigated to date (see section “Discussion” and Figure 6). Nonetheless, the number of species is probably underestimated currently, as recent molecular data have indicated that there are different species complexes within this genus. For example, the cosmopolitan M. cephalus was found to be composed of 15 well supported mitochondrial lineages (Durand and Borsa, 2015), including the one sampled in the type-locality (Mediterranean Sea); six of these lineages have already been cytogenetically analyzed (Rossi et al., 1996, 2016). However, these lineages lack formal descriptions and species name attribution.

Very recently, Durand et al. (2017) reported the presence of two well-supported mitochondrial lineages in the hospe mullet “Mugil hospes,” a species that, according to Barletta and Dantas (2016), is distributed in the western Atlantic from Belize to Brazil and in the eastern Pacific from Mexico to Ecuador. The first molecular lineage includes sequences from Brazil and corresponds to the resurrected species Mugil brevirostris, which is distributed from the northern Brazilian coast (Amapá) to the southern Brazilian coast (Rio Grande do Sul) (Menezes et al., 2015); the second lineage is represented by haplotypes collected in the Gulf of Mexico (Belize/Venezuela) and was named Mugil sp. R. Samples from the eastern Pacific were not included in these analyses or in any other molecular study. The karyotype of the species remains undescribed.

In this research, specimens of M. hospes from Ecuador have been collected and their morphological characters accurately analyzed to make sure of the correct species identification. Cytogenetic and mitochondrial cytochrome c oxidase subunit I (COI) sequence analyses were performed aiming to (a) verify whether the original described M. hospes from the Pacific Ocean corresponds to one of the two lineages identified in the Atlantic Ocean or represents a third lineage, (b) estimate if the divergence among lineages is sufficient to attribute them to different species, (c) identify species-specific chromosome markers and add new comparative data that allow cytotaxonomic inferences on Mugilidae karyotype evolution.

Materials and Methods

Fourteen specimens of Mugil hospes (four males, four females, six immature), were collected with a cast net from a reservoir that provides water to a shrimp pool located at Barbones, El Oro Province, Ecuador (3°09′14.0″ S 79°53′53.1″ W). Fishes were transported to the laboratory in sealed plastic bags (32′) containing two gallons of water, and the air in the bags was replaced with pure oxygen. All 14 individuals were used to prepare cell suspensions. A subsample of eight individuals was used for molecular and morphological analyses. Voucher specimens were deposited in the fish collection of the Laboratório de Biologia e Genética de Peixes (LBP), UNESP, Botucatu (São Paulo State, Brazil) (collection numbers LBP 23325) and Universidad Técnica de Machala (UTMACH-174-UTMACH-182; UTMACH-187; UTMACH-191-UTMACH-194). All experiments were conducted according to the Ethical Committee of Instituto de Biociências/UNESP/Botucatu, under protocol number 1057.

Morphological Analysis

Each fish was measured. Measurements and counts were taken as described by Menezes et al. (2010). Mouth width and mouth depth were measured as described by Thomson (1997). Twenty morphometric characters (Supplementary Table 1) and nine meristic characters (Supplementary Table 2) were recorded for each fish.

Molecular Analysis

Genomic DNA was extracted from muscle tissue that was preserved in 95% ethanol. DNA samples were obtained for eight specimens (one male, two female, five immature), according to procedures described by Aljanabi and Martínez (1997). A 655 bp fragment of the mitochondrial COI was amplified by PCR and sequenced using primers and protocols reported by Nirchio et al. (2017). DNA sequences were aligned using the software Clustal X (Thompson et al., 1997) and deposited in GenBank (Accession numbers: KY964500-KY964504). The basic local alignment search tool (BLAST1) was used to search for similar sequences to confirm species assignment.

For phylogenetic tree reconstruction, a subset of the COI sequences of Mugil, previously analyzed by Durand et al. (2017), was considered. Those sequences that showed greater than 90% similarity (i.e., the six sequences of M. brevirostris and the seven sequences of Mugil sp. R) were also included; Agonostomus monticola (Bancroft, 1834) (JQ060401) was used as an outgroup.

Three types of phylogenetic reconstructions were conducted: neighbour-joining (NJ), maximum-likelihood (ML) and Bayesian inference (BI) analyses. NJ and ML analyses (1000 bootstrap pseudoreplicates) were performed using MEGA7 (Kumar et al., 2016) and PhyML 3.0 (Guindon et al., 2010), respectively. Bayesian analyses were carried out as implemented in MrBayes 3.1.2 (Huelsenbeck and Ronquist, 2001). Two independent runs of four Markov chains, each for 1,000,000 generations were performed. ModelTest 3.7 (Posada and Crandall, 1998) and MrModelTest 2.3 (Nylander, 2008) were used to select, according to the Akaike information criterion, the evolutionary models that best fit the data set for the ML (GTR + I + G, with nst = 6, gamma shape = 4.682, and proportion of invariant sites = 0.637) and the BI (GTR + I + G) analyses, respectively. Genetic distances were calculated with MEGA7 using the Kimura-2-parameters substitution model (Kimura, 1980).

Cytogenetic Analysis

Each fish received an intra-abdominal injection of 0.0125% colchicine (1.0 ml/100 g body weight) 50 min before being sacrificed by administering a numbing overdose of benzocaine (250 mg/L) as recommended by the Guidelines for the Euthanasia of Animals of the American Veterinary Medical Association (AVMA, 2013). Kidney cells were suspended, and chromosomes were prepared by following the conventional air-drying method, as described by Nirchio and Oliveira (2006). Classical staining techniques (Giemsa, Ag-staining, C-banding) and fluorescence in situ hybridization (FISH) were used to map ribosomal gene clusters (5S rDNA and 18S rDNA) and U1 snRNA gene clusters (U1 snRNA is a non-coding RNA that forms part of the spliceosome) (Nilsen, 2003). Telomeric probes were also applied. For the conventional karyotype, slides were stained for 20 min with 10% Giemsa in phosphate buffer at pH 6.88. Active nucleolus organizer regions (NORs) were revealed by silver (Ag) staining as described by Howell and Black (1980); this was performed after Giemsa staining (Rábová et al., 2015). C-banding was performed following the method of Sumner (1972).

The 5S rDNA, 18S rDNA, U1 snRNA genes and telomeric repeats were mapped onto chromosomes by FISH using the method described by Pinkel et al. (1986). Sequences of 5S rDNA, 18S rDNA, U1 snDNA and telomeric repeats were obtained by polymerase chain reaction (PCR) from the genome of Hypsolebias flagellatus and used as probes. The primers used for amplification were 5SA (5′-TCA ACC AAC CAC AAA GAC ATT GGC AC-3′) and 5SB (5′-TAG ACT TCT GGG TGG CCA AAG GAA TCA-3′) (Pendás et al., 1995), 18S6F (5′-CTC TTT CGA GGC CCT GTA AT-3′) and 18S6R (5′-CAG CTT TGC AAC CAT ACT CC-3′) (Utsunomia et al., 2016), U1F (5′-GCA GTC GAG ATT CCC ACA TT-3′) and U1R (5′-CTT ACC TGG CAG GGG AGA TA-3′) (Silva et al., 2015) and (TTAGGG)5 and (CCCTAA)5 (Ijdo et al., 1991). The 5S rDNA and telomeric probes were labeled with biotin-16-dUTP (2′-deoxyuridine 5′-triphosphate), and the 18S rDNA and U1 snRNA gene probes were labeled by including digoxigenin-11-dUTP in the PCR. Hybridization was detected with fluorescein-conjugated avidin (FITC, Sigma–Aldrich2) and anti-digoxigenin-rhodamine conjugate (Roche Applied Science3), respectively. Chromosomes were counterstained with 4,6-diamidino-2-phenylindole (DAPI), which was included in the Vectashield mounting medium (Vector Laboratories4).

Conventionally stained metaphase cells were photographed using a Motic B400, equipped with a Moticam 5000C digital camera using Motic Images Plus 2.0 ML software. FISH images were captured with an Olympus BX61 photomicroscope equipped with a DP70 digital camera using Image-Pro plus 6.0 software (Media Cybernetics). Images were merged and edited for optimization of brightness and contrast using Photoshop (Adobe Systems, Inc.) Version 2015.0.0.

Results

Meristic and Morphometric Characters

The fresh specimens were gray on the dorsal side and white/silver on the ventral side. The pelvic fins had a yellowish tone, and the base of each pectoral fin had a visible dark spot. The dorsal fins and caudal fins were dusky. The distal tips of the anterior rays of the second dorsal fin were slightly darker. The pelvic and anal fins were pale. The body was elongated, with a slightly pointed snout (see Figure 1). The origin of the first dorsal fin was midway between the tip of the snout and the base of the caudal fin. The second dorsal fin and anal fin were profusely covered with scales. One row of small teeth was visible on the upper and lower lips (viewed under the microscope). There were adipose eyelids and widely separated spiny-rayed dorsal fins with four spines in the first dorsal fin and one spine plus eight soft rays in the second dorsal fin (small specimen with nine soft rays). Pelvic fins were sub-abdominal with one spine and 5–6 branched soft rays (commonly I+5). Pectoral fins were long, reaching the level of the origin of the first dorsal fin or extending just beyond, with two spines (the first spine very small) and 11–13 soft rays (commonly 12 rays). The anal fin had three spines and nine soft rays (first spine very short, and hidden by overlying scales). There was a large pectoral axillary scale, with 37–38 scales in longitudinal series (commonly 38), 11–14 scales in an oblique row extending to the origin of the pelvic fin (commonly 13) and 13 scales in a transversal series, as well as 17–22 scales in a circum-peduncular series (commonly 19) (Supplementary Table 2).

FIGURE 1
www.frontiersin.org

FIGURE 1. Specimen of Mugil hospes from Ecuador.

Molecular Analysis

BLAST was used to show that COI nucleotide sequences from GenBank have 92–93% similarity with specimens originally identified as Mugil hospes, and with M. trichodon and Mugil sp., which were all collected in the Atlantic Ocean (from Brazil and Belize). Similarity values with other Mugil species were all below 90%.

The phylogenetic tree obtained by NJ, ML (lnL = -3276.75382), and BI (lnL = -3986.503923) analyses (Figure 2) shows three well-supported lineages of “M. hospes.” The first two correspond to the M. brevirostris (Brazil) and Mugil sp. R (Belize and Venezuela) lineages identified by Durand et al. (2017), whereas the third, referred to hereafter as Mugil hospes (see Discussion), includes all the sequences from Ecuador obtained in this study (Figure 2). The genetic distance is 0.077 between M. hospes/M. brevirostris and between M. hospes/Mugil sp. R, and 0.073 between M. brevirostris/Mugil sp. R.

FIGURE 2
www.frontiersin.org

FIGURE 2. Neighbor-joining tree based on COI sequences. At each node bootstrap values (NJ and ML) and posterior probabilities (BI) are shown. Black dots at the nodes indicate bootstrap values > 70% (NJ and ML) and posterior probabilities > 0.9 (BI). Stars indicate sequences obtained in this study; the remaining sequences are from Durand et al. (2017).

Cytogenetic Analysis

All individuals showed a diploid number of 2n = 48 and a karyotype composed entirely of uniformly decreasing acrocentric chromosomes. Thus, the fundamental arm number (NF) was 48. Only two pairs of homologous chromosomes can be identified with certainty: pair 5, due to a clear interstitial secondary constriction, and pair 24, which is distinctly small (Figure 3A). Sequential Giemsa-silver (Ag) nitrate staining enabled the identification of two actively transcribing NORs, interstitially located on the secondary constriction of chromosome pair 5 (Figures 3B,C). C-banding showed that constitutive heterochromatin is restricted to the centromeric regions of all chromosomes, and there is a pericentromeric heterochromatin block on the secondary constriction of chromosome pair 5 (Figure 3D).

FIGURE 3
www.frontiersin.org

FIGURE 3. (A) Conventional Giemsa-stained karyotype of Mugil hospes. Sequential Giemsa (B) and AgNO3-staining (C) of the metaphase plate (inset shows NOR-bearing chromosomes); arrows indicate the NOR-bearing chromosome. (D) C-banded chromosomes; arrows indicate positive heterochromatic blocks on the interstitial secondary constriction of chromosome pair 5.

Double FISH experiments using 5S and 18S rDNA as probes revealed two positive sites detected for each probe (Figure 4A), located on different chromosome pairs. The 18S rDNA positive sites correspond to the AgNO3 sites on the secondary constriction of chromosome pair 5. The 5S rDNA probes hybridized interstitially on a pair of medium-sized chromosomes. Double FISH using the U1 snDNA and 5S rDNA probes revealed positive U1 snDNA signals on the telomeric region of a pair of medium-sized chromosomes distinct from the 5S-rDNA-bearing chromosome pair (Figure 4B). When chromosomes treated by double FISH were sorted by decreasing size, it was possible to assign the 5S rDNA sites to pair 15 and the U1 snDNA signals to pair 10 (Figure 5). Telomeric repeats were located at both ends of each chromosome, although signal intensities varied between chromosomes (Figure 4C).

FIGURE 4
www.frontiersin.org

FIGURE 4. Somatic metaphase chromosomes of Mugil hospes assayed by FISH and counterstained with DAPI: (A) 5S rDNA (arrows) and 18S rDNA (arrowheads); (B) U1 snRNA (asterisks), 5S rDNA (arrows), chromosome pair 5 (circle), and (C) telomeric repeats. Enlargement of selected samples of chromosome pairs after DAPI staining (left) and FISH (right) are shown in the insets: (A) chromosome pairs 5 and 15, with probes showing 18S rDNA (above) and 5S rDNA (below) positive sites; (B) chromosome pair 5 (above), chromosome pair 10 (center) showing U1 snDNA positive sites and chromosome pair 15 (below) showing 5S rDNA.

FIGURE 5
www.frontiersin.org

FIGURE 5. FISH karyotype. Interstitial secondary constriction corresponding to NOR (chromosome pair 5); 5S rDNA (chromosome pair 15) and U1 snDNA (chromosome pair 10) positive sites are evident.

Discussion

Meristic and morphometric data of samples from Ecuador agree with the original description of M. hospes (Jordan and Culver 1895 in Jordan, 1895): this species has pectoral fins whose tips reach and extend slightly past the vertical line passing through the origin of the first dorsal fin (with four spines). This morphological character is shared with M. brevirostris, which inhabits the opposite side of the Americas (i.e., the Atlantic coast).

Sequence analysis showed that samples of the hospe mullet from Ecuador are genetically very different from those collected in the Atlantic Ocean, all of which were originally identified with the same name M. hospes. Thus, in addition to the two lineages identified by Durand et al. (2017) in the Atlantic Ocean, M. brevirostris and Mugil sp. R, a third lineage is present in the Pacific Ocean. The genetic distances between the Pacific and the two Atlantic lineages are higher than the COI 2% threshold value that discriminates different species (Ward, 2009), and in the phylogenetic reconstruction, the three species form a monophyletic and well-supported clade. Cryptic species are defined as distinct evolutionary lineages not detectable with traditional taxonomic approaches, due to the absence of morphological differences (Avise, 2000; Mallet, 2010). In the last decade, barcoding methods based on COI sequences have made possible their identification in several marine and freshwater fish species (Ward et al., 2008; Lara et al., 2010; Puckridge et al., 2013; Mateussi et al., 2017; Okamoto et al., 2017; Ramirez et al., 2017; Shimabukuro-Dias et al., 2017). In Mugilidae evidences of cryptic species were inferred from mitochondrial tree topology, independent data from nuclear markers, and on the base of the geographic distribution of sister lineages (Durand and Borsa, 2015). Our results indicate that the genetic distances between the different hospe mullet are comparable to those reported among species within both the M. cephalus and M. curema species complexes (Durand et al., 2017), and the three lineages inhabits different geographic areas. Thus, we hypothesize that besides the Mugil cephalus and Mugil curema species complexes, there is an additional putative one, which should be identified as the M. hospes species complex; the name M. hospes should be kept by the Pacific samples, being Mazatlán (in the eastern Pacific) the species type-locality.

Cytogenetic analysis shows that the 48 acrocentric chromosome karyotype detected in M. hospes is consistent with the generally available data on diploid chromosome number and karyotype structure in Mugilidae (Sola et al., 2007; Rossi et al., 2016). This confirms that the only exception is represented by the mullets belonging to the Mugil curema species complex (Nirchio et al., 2017).

Apart from the number of chromosomes, many microstructural changes are evident looking at the variability in the locations of ribosomal genes in Mugil. For example, 5S rDNA cistrons are always localized to an interstitial position, although they are on different chromosomes in different species (Figure 6). The 18S rDNA cistrons seems to be more variable and can be found in the telomeric or interstitial regions of a long chromosome arm, or even on the short arms of different chromosomes. The variability in the localization of the major ribosomal genes could be attributable to their association with heterochromatinas that is observed in Mugil cephalus, M. margaritae (formerly M. curema), M. rubrioculus, M. curema, M. liza, M. trichodon, M. incilis, Mugil sp. O (Rossi et al., 1996, 2005; Nirchio et al., 2005a,b, 2007, 2017; Hett et al., 2011), and M. hospes (present study). Heterochromatin is known to evolve rapidly, and its composition, that includes highly repetitive simple sequences like satellite DNA and transposable elements, is often different even between closely related species. This characteristic might promote rearrangements of the associated genes and might play an important role in reproductive isolation between sister species (Hughes and Hawley, 2009).

FIGURE 6
www.frontiersin.org

FIGURE 6. Idiograms of the karyotypes observed in Mugil species. Re-drawn and updated from Rossi et al. (2016). Black circles: major ribosomal gene locations. Yellow circles: minor ribosomal gene locations.

Cytogenetic mapping of U1 snDNA probes in M. hospes showed the presence of a single U1 gene cluster, located in the terminal position of a chromosomal pair different from the 18S rDNA and 5S bearing chromosomes (Figures 4B, 5). There are no data available on the chromosome mapping of these sequences in other Mugilidae; thus, it is not possible, at this stage, to compare our results with those of other species in the family. However, the mapping of these sequences, combined with other repetitive sequences in other Mugil, might allow the identification of other chromosome re-arrangement. The analysis of chromosome localization of these sequences is restricted to a few other fish species. In Merluccius merluccius (Merlucciidae), multiple interstitial U1 sites are present (García-Souto et al., 2015). In 19 species of cichlids (Cabral-de-Mello et al., 2012), these sites could be either interstitial or terminal on a single st/a chromosome pair, and represent good chromosomal markers that allow the detection of many microstructural chromosomal rearrangements. On the contrary, in five species of Astyanax (Characidae), there is a conserved pattern in the number of U1 sites per genome, and these sequences are frequently associated with 5S rDNA sequences (Silva et al., 2015).

Telomeric DNA repeat sequences were found at the very ends of chromosomes, as observed in 15 different orders of teleosts (Ocalewicz, 2013). In mugilids, telomeric repeats have been mapped in nine species (Gornung et al., 2004; Rossi et al., 2005; Nirchio et al., 2017) and were found also to be interspersed in NORs. Signal intensity variability between chromosomes, as observed in M. hospes, has been previously reported in other fishes (Rocco et al., 2002; Ocalewicz and Dobosz, 2009; Pomianowski et al., 2012), including Mugil species such as M. cephalus (Gornung et al., 2004), M. liza and M. margaritae (Rossi et al., 2005), and Mugil sp. O (Nirchio et al., 2017). This variability is probably due to differing copy numbers of these sequences in the different sites (Lansdorp et al., 1996).

Conclusion

The data presented here confirm that a complex dynamic has played in the karyotype evolution of Mugil, and they reiterate the usefulness of cytogenetic and molecular data in the reconstruction of relationships among taxa within Mugilidae. Species of this family are usually characterized by morphological features that are “insufficient to describe its actual species diversity” (Durand and Borsa, 2015). The combined use of morphological, molecular and cytogenetic analysis is necessary in these fishes to avoid species misidentification and to reconcile the confused picture obtained by morphology-based taxonomy with molecular-based taxonomy. In the case of the M. hospes species complex, Mugil sp. R, which is distributed in the Caribbean Sea, still deserves a formal morphological description and specific name attribution. This species, along with the Brazilian species M. brevirostris, also requires a karyotype description. Thus, at this stage, it is not possible to determine whether this complex is characterized by karyotypes that differ in the total number and morphology of chromosomes, like the M. curema species complex (Nirchio et al., 2017), or whether it is characterized by karyotype homogeneity, like M. cephalus species complex (Rossi et al., 1996).

Author Contributions

MN, CO, and AR designed the study. FP and VM conducted the lab work. MN, AR, and VM designed and conducted the analyses. All authors analyzed the results and wrote the manuscript.

Funding

Funding was provided by: Centro de Investigación of Universidad Técnica de Machala, Ecuador; Fundação de Amparo à Pesquisa do Estado de São Paulo -FAPESP-, Brazil; Conselho Nacional de Desenvolvimento Científico e Tecnológico -CNPq-, Brazil; Sapienza University, Rome (Progetto Ricerca Università 2016).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors are grateful to Dr. Naercio Menezes and Dr. Ricardo Britzke, who helped with the identification of M. hospes. All experiments comply with the current laws of Ecuador and Italy.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00017/full#supplementary-material

Footnotes

  1. ^https://blast.ncbi.nlm.nih.gov/Blast.cgi
  2. ^www.sigma-aldrich.com
  3. ^https://lifescience.roche.com/
  4. ^https://vectorlabs.com/

References

Aljanabi, S. M., and Martínez, L. (1997). Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Res. 25, 4692–4693. doi: 10.1093/nar/25.22.4692

PubMed Abstract | CrossRef Full Text | Google Scholar

Alves-Silva, A. P., and Dergam, J. A. (2015). Cryptic speciation within the Neotropical Cichlid Geophagus brasiliensis (Quoy & Gaimard, 1824) (Teleostei Cichlidae): a new paradigm in karyotypical and molecular evolution. Zebrafish 12, 91–101. doi: 10.1089/zeb.2014.0999

PubMed Abstract | CrossRef Full Text | Google Scholar

Avise, J. C. (2000). Phylogeography. Boston, MA: Harvard University Press.

Google Scholar

AVMA (2013). Guidelines for the Euthanasia of Animals. Version 2013.0.1. Available at: https://www.avma.org/KB/Policies/Documents/euthanasia.pdf (accessed May 27, 2016).

Bancroft, E. N. (1834). “The class Pisces, arranged by the Baron Cuvier, with supplementary additions,” in The Animal Kingdom, eds E. Griffith and C. H. Smith (London: Whittaker & Co.), 680.

Google Scholar

Barletta, M., and Dantas, D. V. (2016). “Biogeography and distribution of Mugilidae in the Americas,” in Biology, Ecology and Culture of Grey Mullet (Mugilidae), eds D. Crosetti and S. J. M. Blaber (Boca Raton, FL: CRC Press), 42–62. doi: 10.1201/b19927-4

CrossRef Full Text

Cabral-de-Mello, D. C., Valente, G. T., Nakajima, R. T., and Martins, C. (2012). Genomic organization and comparative chromosome mapping of the U1 snRNA gene in cichlid fish, with an emphasis in Oreochromis niloticus. Chromosome Res. 20, 279–292. doi: 10.1007/s10577-011-9271-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Durand, J. D., and Borsa, P. (2015). Mitochondrial phylogeny of grey mullets (Acanthopterygii: Mugilidae) suggests high proportion of cryptic species. C. R. Biol. 338, 266–277. doi: 10.1016/j.crvi.2015.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Durand, J. D., Chen, W. J., Shen, K. N., Jamandre, B. W., Blel, H., Diop, K., et al. (2012). Systematics of the grey mullets (Teleostei: Mugiliformes: Mugilidae): molecular phylogenetic evidence challenges two centuries of morphology-based taxonomy. Mol. Phylogenet. Evol. 64, 73–92. doi: 10.1016/j.ympev.2012.03.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Durand, J. D., Hubert, N., Shen, K. N., and Borsa, P. (2017). DNA barcoding grey mullets. Rev. Fish Biol. Fish. 27, 233–243. doi: 10.1007/s11160-016-9457-7

CrossRef Full Text | Google Scholar

Eschmeyer, W. N., and Fong, J. D. (2017). Catalog of Fishes, Species by Family/Subfamily. Electronic Version. Available at: http://researcharchive.calacademy.org/research/ichthyology/catalog/SpeciesByFamily.asp (accessed May 30, 2017).

Google Scholar

García-Souto, D., Troncoso, T., Pérez, M., and Pasantes, J. J. (2015). Molecular cytogenetic analysis of the European hake Merluccius merluccius (Merlucciidae, Gadiformes): U1 and U2 snRNA Gene Clusters Map to the Same Location. PLOS ONE 10:e0146150. doi: 10.1371/journal.pone.0146150

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Castro, M. (2007). Los Peces Representantes de la Familia Mugilidae en Argentina. Ph.D. thesis, Universidad Nacional de Mar del Plata, Argentina.

Google Scholar

González-Castro, M., and Ghasemzadeh, J. (2016). “Morphology and morphometry based taxonomy of Mugilidae,” in Biology, Ecology and Culture of Grey Mullet (Mugilidae), eds D. Crosetti and S. J. M. Blaber (Boca Raton, FL: CRC Press), 1–20.

Google Scholar

Gornung, E., Mannarelli, M. E., Rossi, A. R., and Sola, L. (2004). Chromosomal evolution in Mugilidae (Pisces, Mugiliformes): FISH mapping of the (TTAGGG)n telomeric repeat in the six Mediterranean mullets. Hereditas 140, 158–159. doi: 10.1111/j.1601-5223.2004.01784.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Guindon, S., Dufayard, J. F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321. doi: 10.1093/sysbio/syq010

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrison, I. J., Nirchio, M., Oliveira, C., Ron, E., and Gaviria, J. (2007). A new species of mullet (Teleostei: Mugilidae) from Venezuela, with a discussion on the taxonomy of Mugil gaimardianus. J. Fish Biol. 71, 76–97. doi: 10.1111/j.1095-8649.2007.01520.x

CrossRef Full Text | Google Scholar

Hett, A. K., Nirchio, M., Oliveira, C., Siccha, Z. R., Rossi, A. R., and Sola, L. (2011). Karyotype characterization of Mugil incilis Hancock, 1830 (Mugiliformes: Mugilidae), including a description of an unusual co-localization of major and minor ribosomal genes in the family. Neotrop. Ichthyol. 9, 107–112. doi: 10.1590/S1679-62252011005000005

CrossRef Full Text | Google Scholar

Howell, W. M., and Black, D. A. (1980). Controlled silver staining of nucleolus organizer regions with a protective colloidal developer: a 1-step method. Experientia 36, 1014–1015. doi: 10.1007/BF01953855

PubMed Abstract | CrossRef Full Text | Google Scholar

Huelsenbeck, J. P., and Ronquist, F. (2001). MRBAYES: bayesian inference of phylogenetic trees. Bioinformatics 17, 754–755. doi: 10.1093/bioinformatics/17.8.754

CrossRef Full Text | Google Scholar

Hughes, S. E., and Hawley, R. S. (2009). Heterochromatin: a rapidly evolving species barrier. PLOS Biol. 7:e1000233. doi: 10.1371/journal.pbio.1000233

PubMed Abstract | CrossRef Full Text | Google Scholar

Ijdo, J. W., Wells, R. A., Baldini, A., and Reeders, S. T. (1991). Improved telomere detection using a telomere repeat probe (TTAGGG)n generated by PCR. Nucleic Acids Res. 19:4780. doi: 10.1093/nar/19.17.4780

PubMed Abstract | CrossRef Full Text | Google Scholar

Jordan, D. S. (1895). The fishes of Sinaloa. Proc. California Acad. Sci. 5, 377–514.

Google Scholar

Kimura, M. (1980). A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Biol. 16, 111–120. doi: 10.1007/BF01731581

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 33, k1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Lansdorp, P. P. M., Verwoerd, N. P., Van den Rijke, F. M., Dragowska, V., Little, M. T., Dirks, R. W., et al. (1996). Heterogeneity in telomere length of human chromosomes. Hum. Mol. Genet. 5, 685–691. doi: 10.1093/hmg/5.5.685

CrossRef Full Text | Google Scholar

Lara, A., Ponce de León, J. L., Rodríguez, R., Casane, D., Côté, G., Bernatchez, L., et al. (2010). DNA barcoding of Cuban freshwater fishes: evidence for cryptic species and taxonomic conflicts. Mol. Ecol. Resour. 10, 421–430. doi: 10.1111/j.1755-0998.2009.02785.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mallet, J. (2010). Why was Darwin’s view of species rejected by twentieth century biologists? Biol. Philos. 25, 497–527. doi: 10.1007/s10539-010-9213-7

CrossRef Full Text | Google Scholar

Mateussi, N. T. B., Pavanelli, C. S., and Oliveira, C. (2017). Molecular identification of cryptic diversity in species of cis-Andean Mylossoma (Characiformes: Serrasalmidae). Mitochondrial DNA A DNA Mapp. Seq. Anal. 28, 778–780. doi: 10.1080/24701394.2016.1180515

PubMed Abstract | CrossRef Full Text | Google Scholar

Menezes, N. A., Nirchio, M., Oliveira, C., and Siccha Ramirez, R. (2015). Taxonomic review of the species of Mugil (Teleostei: Perciformes: Mugilidae) from the Atlantic South Caribbean and South America, with integration of morphological, cytogenetic and molecular data. Zootaxa 3918, 1–38. doi: 10.11646/zootaxa.3918.1.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Menezes, N. A., Oliveira, C., and Nirchio, M. (2010). An old taxonomic dilemma: the identity of the western south Atlantic lebranche mullet (Teleostei: Perciformes: Mugilidae). Zootaxa 2519, 59–68. doi: 10.5281/zenodo.196200

CrossRef Full Text | Google Scholar

Nilsen, T. W. (2003). The spliceosome: the most complex macromolecular machine in the cell? Bioessays 25, 1147–1149.

Google Scholar

Nirchio, M., Cipriano, R. R., Cestari, M. M., and Fenocchio, A. S. (2005a). Cytogenetical and morphological features reveal significant differences among Venezuelan and Brazilian samples of Mugil curema. Neotrop. Ichthyol. 3, k99–102. doi: 10.1590/S1679-62252005000100006

CrossRef Full Text | Google Scholar

Nirchio, M., Ron, E., and Rossi, A. R. (2005b). Karyological characterization of Mugil trichodon Poey, 1876 (Pisces: Mugilidae). Sci. Mar. 69, 525–530. doi: 10.3989/scimar.2005.69n4525

CrossRef Full Text | Google Scholar

Nirchio, M., and Oliveira, C. (2006). Citogenética de Peces. Cumaná: Rectorado de la Universidad de Oriente.

Nirchio, M., Oliveira, C., Ferreira, I. A., Pérez, J. E., Gaviria, J. I., Harrison, I., et al. (2007). Comparative cytogenetic and allozyme analysis of Mugil rubrioculus and M. curema (Teleostei: Mugilidae) from Venezuela. Interciencia 32, k757–762.

Google Scholar

Nirchio, M., Oliveira, C., Siccha Ramirez, Z. R., de Sene, V. F., Sola, L., Milana, V., et al. (2017). The Mugil curema species complex: a new karyotype for the Pacific white mullet mitochondrial lineage. Comp. Cytogenet. 11, 225–237. doi: 10.3897/CompCytogen.v11i2.11579

PubMed Abstract | CrossRef Full Text | Google Scholar

Nylander, J. A. A. (2008). MrModeltest v2.3. Program Distributed by the Author. Uppsala: Uppsala University.

Google Scholar

Ocalewicz, K. (2013). Telomeres in Fishes. Cytogenet. Genome Res. 141, 114–125. doi: 10.1159/000354278

PubMed Abstract | CrossRef Full Text | Google Scholar

Ocalewicz, K., and Dobosz, S. (2009). Karyotype variation in the albino rainbow trout (Oncorhynchus mykiss Walbaum). Genome 52, 347–352. doi: 10.1139/g09-009

PubMed Abstract | CrossRef Full Text | Google Scholar

Okamoto, N., Maung-Saw-Htoo-Thaw, Sein-Thaung, Koike, K., Shibata, J., and Tomiyama, T. (2017). Species composition of hairtails (Trichiuridae) in Myanmar. Reg. Stud. Mar. Sci. 17, 73–77. doi: 10.1016/j.rsma.2017.11.007

CrossRef Full Text | Google Scholar

Pansonato-Alves, J. C., Serrano,É. A., Utsunomia, R., Camacho, J. P. M., Costa Silva, G. J. D., Vicari, M. R., et al. (2014). Single origin of sex chromosomes and multiple origins of B chromosomes in fish genus Characidium. PLOS ONE 9:e107169. doi: 10.1371/journal.pone.0107169

PubMed Abstract | CrossRef Full Text | Google Scholar

Pendás, A. M., Morán, P., Martínez, J. L., and Garcia-Vásquez, E. (1995). Applications of 5S rDNA in Atlantic salmon, brown trout, and in Atlantic salmon x brown trout hybrid identification. Mol. Ecol. 4, 275–276. doi: 10.1111/j.1365-294X.1995.tb00220.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pinkel, D., Straume, T., and Gray, J. W. (1986). Cytogenetic analysis using quantitative, high-sensitivity, fluorescence hybridization. Proc. Natl. Acad. Sci. U.S.A. 83, 2934–2938. doi: 10.1073/pnas.83.9.2934

CrossRef Full Text | Google Scholar

Pomianowski, L., Jankun, M., and Ocalewicz, K. (2012). Detection of interstitial telomeric sequences in the Arctic charr (Salvelinus alpinus, Linnaeus 1758) (Teleostei, Salmonidae). Genome 55, 26–32. doi: 10.1139/G11-073

PubMed Abstract | CrossRef Full Text | Google Scholar

Posada, D., and Crandall, K. A. (1998). MODELTEST: testing the model of DNA substitution. Bioinformatics 14, 817–818. doi: 10.1093/bioinformatics/14.9.817

PubMed Abstract | CrossRef Full Text | Google Scholar

Puckridge, M., Andreakis, N., Appleyard, S. A., and Ward, R. D. (2013). Cryptic diversity in flathead fishes (Scorpaeniformes: Platycephalidae) across the Indo-West Pacific uncovered by DNA barcoding. Mol. Ecol. Resour. 13, 32–42. doi: 10.1111/1755-0998.12022

PubMed Abstract | CrossRef Full Text | Google Scholar

Rábová, M., Völker, M., Pelikánová, S., and Ráb, P. (2015). “Sequential chromosome banding in fishes,” in Fish Cytogenetic Techniques – Ray-fin fishes and Chondrichthyans, eds C. Ozouf-Costaz, E. Pisano, F. Foresti, and L. F. Almeida-Toledo (Boca Raton, FL: CRC Press), 92–102. doi: 10.1201/b18534-15

CrossRef Full Text | Google Scholar

Ramirez, J. L., Birindelli, J. L., Carvalho, D. C., Affonso, P. R. A. M., Venere, P. C., Ortega, H., et al. (2017). Revealing hidden diversity of the underestimated neotropical ichthyofauna: DNA barcoding in the recently described genus Megaleporinus (Characiformes: Anostomidae). Front. Genet. 8:149. doi: 10.3389/fgene.2017.00149

PubMed Abstract | CrossRef Full Text | Google Scholar

Rocco, L., Morescalchi, M. A., Costagliola, D., and Stingo, V. (2002). Karyotype and genome characterization in four cartilaginous fishes. Gene 295, 289–298. doi: 10.1016/S0378-1119(02)00730-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossi, A. R., Crosetti, D., Gornung, E., and Sola, L. (1996). Cytogenetic analysis of global populations of Mugil cephalus (striped mullet) by different staining techniques and fluorescent in situ hybridization. Heredity 76, 77–82. doi: 10.1038/hdy.1996.10

CrossRef Full Text | Google Scholar

Rossi, A. R., Crosetti, D., and Livi, S. (2016). “Genetics of Mugilidae,” in Biology, Ecology and Culture of Grey Mullet (Mugilidae), eds D. Crosetti and S. J. M. Blaber (Boca Raton, FL: CRC Press), 349–397. doi: 10.1201/b19927-16

CrossRef Full Text | Google Scholar

Rossi, A. R., Gornung, E., Sola, L., and Nirchio, M. (2005). Comparative molecular cytogenetic analysis of two congeneric species, Mugil curema and M. liza, characterized by significant karyotype diversity. Genetica 125, 27–32. doi: 10.1007/s10709-005-2668-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Schultz, L. P. (1946). A revision of the genera of mullets, fishes of the genera Mugilidae, with description of three new genera. Proc. U. S. Natl. Mus. 96, 377–395. doi: 10.5479/si.00963801.96-3204.377

CrossRef Full Text | Google Scholar

Shimabukuro-Dias, C. K., da Costa Silva, G. J., Ashikaga, F. Y., Foresti, F., and Oliveira, C. (2017). Molecular identification of the fish fauna from the pantanal flood plain area in Brazil. Mitochondrial DNA A DNA Mapp. Seq. Anal. 28, 588–592. doi: 10.3109/24701394.2016.1149826

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva, D. M., Utsunomia, R., Pansonato-Alves, J. C., Oliveira, C., and Foresti, F. (2015). Chromosomal mapping of repetitive DNA sequences in five species of Astyanax (Characiformes, Characidae) reveals independent location of U1 and U2 snRNA sites and association of U1 snRNA and 5S rDNA. Cytogenet. Genome Res. 146, 144–152. doi: 10.1159/000438813

PubMed Abstract | CrossRef Full Text | Google Scholar

Sola, L., Gornung, E., Mannarelli, M. E., and Rossi, A. R. (2007). “Chromosomal evolution in Mugilidae, Mugilomorpha: an overview,” in Fish Cytogenetics, eds E. Pisano, C. Ozouf-Costaz, F. Foresti, and B. G. Kapoor (Enfield, NH: Science Publishers), 165–194.

Google Scholar

Sumner, A. T. (1972). A simple technique for demonstrating centromeric heterochromatin. Exp. Cell Res. 75, 304–306. doi: 10.1016/0014-4827(72)90558-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, J. D., Gibson, T. J., Plewniak, F., Jeanmougin, F., and Higgins, D. G. (1997). The CLUSTAL X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25, 4876–4882. doi: 10.1093/nar/25.24.4876

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, J. M. (1966). The grey mullets. Oceanogr. Mar. Biol. 4, 301–335.

Google Scholar

Thomson, J. M. (1997). The Mugilidae of the world. Mem. Queensl. Mus. 43, 457–562.

Google Scholar

Utsunomia, R., Silva, D. M. Z. A., Ruiz-Ruano, F. J., Araya-Jaime, C., Pansonato-Alves, J. C., Scacchetti, P. C., et al. (2016). Uncovering the ancestry of B chromosomes in Moenkhausia sanctaefilomenae (Teleostei, Characidae). PLOS ONE 11:e0150573. doi: 10.1371/journal.pone.0150573

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, R. D. (2009). DNA barcode divergence among species and genera of birds and fishes. Mol. Ecol. Res. 9, 1077–1085. doi: 10.1111/j.1755-0998.2009.02541.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, R. D., Costa, F. O., Holmes, B. H., and Steinke, D. (2008). DNA barcoding of shared fish species from the North Atlantic and Australasia: minimal divergence for most taxa, but Zeus faber and Lepidopus caudatus each probably constitute two species. Aquat. Biol. 3, 71–78. doi: 10.3354/ab00068

CrossRef Full Text | Google Scholar

Yano, C. F., Poltronieri, J., Bertollo, L. A. C., Artoni, R. F., Liehr, T., and Cioffi, M. B. (2014). Chromosomal mapping of repetitive DNAs in Triportheus trifurcatus (Characidae, Characiformes): insights into the differentiation of the Z and W chromosomes. PLOS ONE 9:e90946. doi: 10.1371/journal.pone.0090946

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: fish cytogenetics, fish molecular phylogeny, COI, chromosomal evolution, FISH, Mugilidae

Citation: Nirchio M, Paim FG, Milana V, Rossi AR and Oliveira C (2018) Identification of a New Mullet Species Complex Based on an Integrative Molecular and Cytogenetic Investigation of Mugil hospes (Mugilidae: Mugiliformes). Front. Genet. 9:17. doi: 10.3389/fgene.2018.00017

Received: 08 July 2017; Accepted: 15 January 2018;
Published: 05 February 2018.

Edited by:

Roberto Ferreira Artoni, Ponta Grossa State University, Brazil

Reviewed by:

Marcelo R. S. Briones, Federal University of São Paulo, Brazil
Marcelo De Bello Cioffi, Federal University of São Carlos, Brazil

Copyright © 2018 Nirchio, Paim, Milana, Rossi and Oliveira. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Claudio Oliveira, claudio@ibb.unesp.br

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.