ORIGINAL RESEARCH article

Front. Genet., 12 October 2017

Sec. Evolutionary, Population, and Conservation Genetics

Volume 8 - 2017 | https://doi.org/10.3389/fgene.2017.00149

Revealing Hidden Diversity of the Underestimated Neotropical Ichthyofauna: DNA Barcoding in the Recently Described Genus Megaleporinus (Characiformes: Anostomidae)

  • 1. Laboratório de Biodiversidade Molecular e Conservação, Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Paulo, Brazil

  • 2. Departamento de Biologia Animal e Vegetal, Universidade Estadual de Londrina, Londrina, Brazil

  • 3. Laboratório de Genética da Conservação, Programa de Pós-Graduação em Biologia de Vertebrados, PUC Minas, Belo Horizonte, Brazil

  • 4. Departamento de Ciências Biológicas, Universidade Estadual do Sudoeste da Bahia, Jequié, Brazil

  • 5. Departamento de Biologia e Zoologia, Universidade Federal de Mato Grosso, Cuiabá, Brazil

  • 6. Departamento de Ictiología, Museo de Historia Natural, Universidad Nacional Mayor de San Marcos, Lima, Peru

  • 7. Facultad de Ciencias Exactas y Naturales, Universidad Surcolombiana, Huila, Colombia

  • 8. Grupo de Investigación en Genética y Reproducción Animal, Universidad de los Llanos, Villavicencio, Colombia

Article metrics

View details

67

Citations

6,7k

Views

1,9k

Downloads

Abstract

Molecular studies have improved our knowledge on the neotropical ichthyofauna. DNA barcoding has successfully been used in fish species identification and in detecting cryptic diversity. Megaleporinus (Anostomidae) is a recently described freshwater fish genus within which taxonomic uncertainties remain. Here we assessed all nominal species of this genus using a DNA barcode approach (Cytochrome Oxidase subunit I) with a broad sampling to generate a reference library, characterize new molecular lineages, and test the hypothesis that some of the nominal species represent species complexes. The analyses identified 16 (ABGD and BIN) to 18 (ABGD, GMYC, and PTP) different molecular operational taxonomic units (MOTUs) within the 10 studied nominal species, indicating cryptic biodiversity and potential candidate species. Only Megaleporinus brinco, Megaleporinus garmani, and Megaleporinus elongatus showed correspondence between nominal species and MOTUs. Within six nominal species, a subdivision in two MOTUs was found, while Megaleporinus obtusidens was divided in three MOTUs, suggesting that DNA barcode is a very useful approach to identify the molecular lineages of Megaleporinus, even in the case of recent divergence (< 0.5 Ma). Our results thus provided molecular findings that can be used along with morphological traits to better define each species, including candidate new species. This is the most complete analysis of DNA barcode in this recently described genus, and considering its economic value, a precise species identification is quite desirable and fundamental for conservation of the whole biodiversity of this fish.

Introduction

Neotropical freshwater fishes have a remarkable diversity, exceeding 8000 species (Reis et al., 2016), however, much taxonomic uncertainty exists leading to underestimated diversity (Pereira et al., 2013; Reis et al., 2016). Molecular studies have been crucial to improve our knowledge on the ichthyofauna, and DNA barcoding has successfully been used in fish species identification and in detecting species of taxonomic concerns or cryptic diversity (Pereira et al., 2013; Gomes et al., 2015; Ramirez and Galetti, 2015; Machado et al., 2016). Within the neotropical freshwater fishes, the order Characiformes represents more than 30% of the known species, and Anostomidae is one of the most species-rich families, occurring in all major hydrographic basins, with trans- and cis-Andean distribution in South America (Reis et al., 2003).

Comprising approximately 150 described species, distributed in 15 genera (Garavello and Britski, 2003; Sidlauskas and Vari, 2008; Ramirez et al., 2017), the known diversity of the Anostomidae has increased in recent years. For instance, 14 species and 1 genus were described only in the last 5 years (Birindelli et al., 2013; Burns et al., 2014). DNA barcoding has revealed taxonomic uncertainties within the genus Laemolyta (Ramirez and Galetti, 2015), and molecular phylogeny has helped to provide an understanding of the evolutionary history of the Anostomidae (Ramirez and Galetti, 2015; Ramirez et al., 2016, 2017).

Recently, the genus Megaleporinus (Ramirez et al., 2017) was described to include 16 lineages, corresponding to 10 nominal species, previously recognized in Leporinus or Hypomasticus (Ramirez et al., 2017). Megaleporinus is supported by cytogenetic, molecular, and morphological data. It is characterized by having a unique ZZ/ZW sex chromosome system (Galetti et al., 1995), while most cytogenetically known Leporinus species have no sex chromosomes (Galetti et al., 1981, 1991). Its monophyly is also well supported by mitochondrial and nuclear markers, which identified it as the sister group to Abramites (Ramirez et al., 2017). Concerning its morphology, Megaleporinus is characterized by being relatively large (adults usually reaching more than 35 cm standard length, including the largest species of the family), three teeth on each premaxillary and dentary bones, and a color pattern of one to three dark mid-lateral blotches (Ramirez et al., 2017). Because of its large size, Megaleporinus has an economic importance in subsistence fisheries and aquaculture (Garavello and Britski, 2003).

Recent studies indicate that there is a hidden biodiversity within Megaleporinus that needs to be better understood (Avelino et al., 2015; Ramirez et al., 2017). A study based on mitochondrial and nuclear markers, but using few individuals for each species, showed that several nominal species allocated to this genus comprise two or more molecular lineages allopatrically distributed in different basins (Ramirez et al., 2017).

In this study, we used a DNA barcoding approach to generate a reference library for Megaleporinus, assessing all nominal species and lineages previously described. We included a broad sampling for most of the species. Our hypothesis is that DNA barcoding support the observation that some of the nominal species represent species complexes with most molecular operational taxonomic units (MOTUs) allopatrically distributed in different basins, as proposed by Ramirez et al. (2017). Identifying such hidden biodiversity within this genus, this paper will contribute to a more complete understanding of its diversity and to the conservation of this important fish group.

Materials and Methods

Sampling

Animals were collected on public land, handled and killed under permission (ICMBIO/MMA N° 32215) provided by the Environment Ministry (MMA). This study did not involve endangered or protected species. Fish were collected by fishing rods and gillnets. No ethics committee approval is required for these organisms in Brazil. Fish were killed in the field using cold water and immediately transferred onto ice. Tissue samples were collected after fish death was confirmed through lack of operculum movement.

Specimens from several populations of all Megaleporinus species were used in this study, totaling 79 samples of the 10 nominal species, and comprising the 16 molecular lineages described by Ramirez et al., 2017 (Figures 1, 2 and Table 1). Voucher numbers are provided for the specimens (Table 1). Additionally, previous DNA barcode sequences of specimens from the São Francisco (Carvalho et al., 2011), Paraná (Pereira et al., 2013), Paranapanema (Frantine-Silva et al., 2015), and lower Paraná basins (Díaz et al., 2016) were included in our data set giving a total of 116 sequences (Figures 1, 2 and Table 1).

FIGURE 1

FIGURE 2

Table 1

MOTURiver (Locality)BasinBINGenBankMuseum ID
Megaleporinus brincoContas12ContasADB0463KU134850MZUSP – 118670
M. brincoContas12ContasADB0463KX925449MZUSP – 118670
M. brincoContas12ContasADB0463KX925450MZUSP – 118670
M. cf. conirostrisDoce (Governador Valadares)18DoceACL4264KF568977MCNI-PUCMG-0186
M. cf. conirostrisDoce (Governador Valadares)18DoceACL4264KX925451
M. cf. conirostrisDoce (Governador Valadares)18DoceACL4264KX925452
M. cf. conirostrisDoce (Baixo Guandú)19DoceACL4264KX925453LISDEBE 6971
M. cf. muyscorumMeta (Puerto Lopez)2OrinocoADB0512KU134851
M. cf. obtusidens ParaguayCuiaba (Santo Antonio de Leverger)8ParaguaiACL3942KU134861MZUSP – 118668
M. cf. obtusidens São FranciscoSão Francisco (Tres Marias)20São FranciscoABZ0928HM405029
M. cf. obtusidens São FranciscoPandeiros14São FranciscoABZ0928HM405142
M. cf. obtusidens São FranciscoUrucuia (Urucuia)13São FranciscoABZ0928HM906022
M. cf. obtusidens São FranciscoUrucuia (Urucuia)13São FranciscoABZ0928HM906023
M. cf. obtusidens São FranciscoSão Francisco (Tres Marias)16São FranciscoABZ0928HM405028
M. cf. obtusidens São FranciscoSão Francisco (Pirapora)27São FranciscoABZ0928KX925498LISDEBE 6973
M. cf. obtusidens São FranciscoSão Francisco (Pirapora)27São FranciscoABZ0928KX925499
M. cf. obtusidens São FranciscoSão Francisco (Pirapora)27São FranciscoABZ0928KX925500
M. cf. obtusidens São FranciscoSão Francisco (Pirapora)27São FranciscoABZ0928KX925501
M. cf. obtusidens São FranciscoSão Francisco (Pirapora)27São FranciscoABZ0928KX925502
M. cf. obtusidens São FranciscoPandeiros (Pandeiros)14São FranciscoABZ0928KU134862MCP – 44805
M. cf. obtusidens São FranciscoUrucuia (Urucuia)13São FranciscoABZ0928KX925503
M. cf. obtusidens São FranciscoUrucuia (Urucuia)13São FranciscoABZ0928KX925504MCP – 44076
M. cf. reinhardtiItapicurú11Itapicurú-mirimAAD1729KU134849
M. cf. reinhardtiItapicurú11Itapicurú-mirimAAD1729KX925454
M. cf. reinhardtiItapicurú11Itapicurú-mirimAAD1729KX925455
M. cf. reinhardtiItapicurú11Itapicurú-mirimAAD1729KX925456UESB-8206
M. cf. trifasciatusAraguaia (Ouro fino)10TocantinsACL3074KX925457GEPEMA – 4975
M. cf. trifasciatusAraguaia (Barra do Garças)10TocantinsACL3074KF568998GEPEMA – 5095
M. cf. trifasciatusAraguaia (Barra do Garças)10TocantinsACL3074KX925458GEPEMA – 5594
M. conirostrisParaibuna20Paraiba do SulACL3731KU134852
M. conirostrisParaibuna20Paraiba do SulACL3731KX925459
M. elongatusItacambiruçu (Grão Mogol)15JequitinhonhaABY2894KX925463
M. elongatusItacambiruçu (Grão Mogol)15JequitinhonhaABY2894KU134853MCNI-PUCMG-0375
M. elongatusJequitinhonhaABY2894KX925464MCNI-PUCMG-0221
M. elongatusJequitinhonha (UHE Irapé)29JequitinhonhaABY2894KU134854
M. elongatusJequitinhonha (UHE Irapé)29JequitinhonhaABY2894KX925465MCNI-PUCMG-0299
M. elongatusJequitinhonha (UHE Irapé)29JequitinhonhaABY2894KX925466MCNI-PUCMG-0300
M. elongatusRio Pardo (Águas Vermelhas)28PardoABY2894KX925460MCNI-PUCMG-4451
M. elongatusRio Pardo (Águas vermelhas)28PardoABY2894KX925461MCNI-PUCMG-5175
M. elongatusRio Pardo (Águas vermelhas)28PardoABY2894KX925462MCNI-PUCMG-5176
M. garmaniItacambiruçu (Grão Mogol)15JequitinhonhaACL3227KU134855MCNI-PUCMG-0021
M. garmaniItacambiruçu (Grão Mogol)15JequitinhonhaACL3227KX925467MCNI-PUCMG-0020
M. garmaniItacambiruçu (Grão Mogol)15JequitinhonhaACL3227KX925468MCNI-PUCMG-0021
M. garmaniItacambiruçu (Grão Mogol)15JequitinhonhaACL3227KX925469MCNI-PUCMG-0021
M. garmaniItacambiruçu (Grão Mogol)15JequitinhonhaACL3227KX925470MCNI-PUCMG-0374
M. macrocephalusManhuaçu (São José do Mantimento)DoceAAE5328KX925475MCNI-PUCMG-0460
M. macrocephalusCuiaba (Santo Antonio de Leverger)8ParaguaiAAE5328KU134856MZUSP – 118667
M. macrocephalusCuiaba (Santo Antonio de Leverger)8ParaguaiAAE5328KX925471LISDEBE 6972
M. macrocephalusCuiaba (Cuiaba)7ParaguaiAAE5328KX925474
M. macrocephalusCuiaba (Barão de Melgaço)9ParaguaiAAE5328KX925476LISDEBE 6974
M. macrocephalusCuiaba (Barão de Melgaço)9ParaguaiAAE5328KX925477LISDEBE 6974
M. macrocephalusTiete (Barra Bonita)ParanáAAE5328KX925473
M. macrocephalusPandeirosSão FranciscoAAE5328HM906021
M. macrocephalusAraguaia (Ouro fino)TocantinsAAE5328KX925472GEPEMA – 4974
M. macrocephalus ParanáCinzas (Bandeirantes)25ParanáACO1303KM897611
M. macrocephalus ParanáCinzas (Bandeirantes)25ParanáACO1303KM897537
M. macrocephalus ParanáCinzas (Bandeirantes)25ParanáACO1303KM897575
M. macrocephalus ParanáCinzas (Bandeirantes)25ParanáACO1303KM897296
M. macrocephalus ParanáPiracicaba (Tamanduá)23ParanáACO1303JN988999LBPV-19469
M. muyscorumMagdalena (Neiva)1MagdalenaADB0701KX925478ICN-19072
M. muyscorumMagdalena (Neiva)1MagdalenaADB0701KX925479ICN-19073
M. muyscorumMagdalena (Neiva)1MagdalenaADB0701KU134857ICN-19074
M. obtusidensJacuí (Jacuizinho Foz)26JacuíAAB8578KU134859MCP-25476
M. obtusidensPiracicaba (Tamanduá)23ParanáAAB8578JN988985LBPV-19849
M. obtusidensPiracicaba (Tamanduá)23ParanáAAB8578JN988984LBPV-19850
M. obtusidensPiracicaba (Tamanduá)23ParanáAAB8578JN988983LBPV-19852
M. obtusidensParanapanema (Canoas)22ParanáAAB8578KM897227
M. obtusidensParanapanema (Canoas)22ParanáAAB8578KM897138
M. obtusidensParanapanema (Canoas)22ParanáAAB8578KM897434
M. obtusidensCinzas (Bandeirantes)25ParanáAAB8578KM897597
M. obtusidensTurvo (Icem)21ParanáAAB8578KX925480LISDEBE 6969
M. obtusidensTurvo (Icem)21ParanáAAB8578KX925481
M. obtusidensTurvo (Icem)21ParanáAAB8578KX925482
M. obtusidensTurvo (Icem)21ParanáAAB8578KX925483
M. obtusidensTurvo (Icem)21ParanáAAB8578KU134858MZUSP – 113982
M. obtusidensTurvo (Icem)21ParanáAAB8578KF568987
M. obtusidensParaná (Porto Camargo)24ParanáAAB8578KX925484
M. obtusidensParaná (Porto Camargo)24ParanáAAB8578KX925485
M. obtusidensIbicui (BR 472)25UruguayAAB8578KU134860MCP-28917
M. piavussuPiracicaba (Tamanduá)23ParanáAAB8569JN989005LBPV-15587
M. piavussuPiracicaba (Tamanduá)23ParanáAAB8569JN989004LBPV-19851
M. piavussuPiracicaba (Tamanduá)23ParanáAAB8569JN989003LBPV-19854
M. piavussuParanapanema (Canoas)22ParanáAAB8569KM897529
M. piavussuParanapanema (Canoas)22ParanáAAB8569KM897489
M. piavussuParanapanema (Canoas)22ParanáAAB8569KM897621
M. piavussuCinzas (Bandeirantes)25ParanáAAB8569KM897419
M. piavussuParanapanema (Canoas II)31ParanáAAB8569KM897506
M. piavussuParanapanema (Canoas II)31ParanáAAB8569KM897347
M. piavussuParanapanema (Canoas II)31ParanáAAB8569KM897192
M. piavussuTurvo (Icem)21ParanáAAB8569KF568991MZUSP – 113981
M. piavussuTurvo (Icem)21ParanáAAB8569KX925486LISDEBE 6968
M. piavussuTurvo (Icem)21ParanáAAB8569KX925487LISDEBE 6970
M. piavussuTurvo (Icem)21ParanáAAB8569KX925488LISDEBE 6970
M. piavussuTurvo (Icem)21ParanáAAB8569KX925489
M. piavussuTurvo (Icem)21ParanáAAB8569KX925490
M. piavussuTurvo (Icem)21ParanáAAB8569KX925491
M. piavussuParaná (Porto Camargo)24ParanáAAB8569KX925492
M. piavussuParaná (Pauliceia)30ParanáAAB8569KX925493
M. cf. piavussu lower ParanáParaná (Rosario)32ParanáAAB8569KU288864
M. cf. piavussu lower ParanáParaná (Rosario)32ParanáAAB8569KU288865
M. cf. piavussu lower ParanáParaná (Rosario)32ParanáAAB8569KU288866
M. cf. piavussu lower ParanáParaná (Rosario)32ParanáAAB8569KU289030
M. reinhardtiPandeiros14São FranciscoAAD1729HM906025
M. reinhardtiCurimataí17São FranciscoAAD1729HM405147MCP – 44776
M. reinhardtiCurimataí17São FranciscoAAD1729HM906026
M. reinhardtiUrucuia (Urucuia)13São FranciscoAAD1729HM906027
M. reinhardtiUrucuia (Urucuia)13São FranciscoAAD1729HM906028
M. reinhardtiSão Francisco (Três Marias)16São FranciscoAAD1729KX925494
M. reinhardtiUrucuia (Urucuia)13São FranciscoAAD1729KX925495
M. reinhardtiUrucuia (Urucuia)13São FranciscoAAD1729KX925496
M. reinhardtiCurimataí17São FranciscoAAD1729KX925497MCP - 44770
M. trifasciatusMadeira6AmazonasACL3073KU134864UFRO-I 4902
M. trifasciatusUcayali (Pucallpa)5AmazonasACL3073KU134865MUSM - 47351
M. trifasciatusUcayali (Pucallpa)5AmazonasACL3073KX925505MUSM - 47351
M. trifasciatusAmazonas (Belen)3AmazonasACL3073KX925506MUSM - 47364
M. trifasciatusLago Catalão (Manaus)4AmazonasACL3073KX925507INPA 11641
M. reinhardtiUrucuia (Urucuia)13São FranciscoAAD1729HM906028
M. reinhardtiSão Francisco (Três Marias)16São FranciscoAAD1729KX925494
M. reinhardtiUrucuia (Urucuia)13São FranciscoAAD1729KX925495
M. reinhardtiUrucuia (Urucuia)13São FranciscoAAD1729KX925496
M. reinhardtiCurimataí17São FranciscoAAD1729KX925497MCP - 44770
M. trifasciatusMadeira6AmazonasACL3073KU134864UFRO-I 4902
M. trifasciatusUcayali (Pucallpa)5AmazonasACL3073KU134865MUSM - 47351
M. trifasciatusUcayali (Pucallpa)5AmazonasACL3073KX925505MUSM - 47351
M. trifasciatusAmazonas (Belen)3AmazonasACL3073KX925506MUSM - 47364
M. trifasciatusLago Catalão (Manaus)4AmazonasACL3073KX925507INPA 11641

Sampling information and GenBank accession for specimens included in the analysis.

Obtained from BOLD.

DNA Extraction, Amplification, and Sequencing

Total DNA was extracted from tissues (fins, muscle, or liver) by the standard phenol–chloroform method (Sambrook et al., 1989). A fragment of Cytochrome Oxidase subunit I (COI; 698 bp) was amplified via polymerase chain reaction (PCR) using primers AnosCOIF and AnosCOIR (Ramirez and Galetti, 2015). PCR products were sequenced for both strands using an ABI 3730xl (Applied Biosystems, Waltham, MA, United States) automatic sequencer. Contigs were assembled and edited using BioEdit (Hall, 1999). All sequences were evaluated manually, deleting regions of low quality. All sequences were verified to represent the COI gene and were checked for indels and stop codons. GenBank (Benson et al., 2017) accession numbers are given in Table 1. All information about specimen, sequences, and electropherograms were deposited in a data set of The Barcode of Life Database platform (BOLD) with code DS-MGLEP.

DNA Barcode Analysis

The general mixed Yule coalescent (GMYC) model (Pons et al., 2006) with a single threshold, implemented in the splits packages in the R 3.3.3 statistical software (R Core Team, 2017), was used to infer MOTUs. For the GMYC input, an ultrametric tree was generated using Beast 2.4.3 (Bouckaert et al., 2014), with a lognormal relaxed clock, a birth and death model, and a GTR+G substitution model, chosen using jModeltest 2 (Darriba et al., 2012), using 50 million MCMC generations and a burn-in of 10%. Poisson tree processes (PTP) model (Zhang et al., 2013) was used for MOTUs delimitation through the bPTP server1, using default values. The bPTP server includes a Bayesian implementation of the PTP model and the original maximum likelihood PTP. For the PTP input, a tree was generated using Beast 2.4.6 (Bouckaert et al., 2014), with a strict clock, a birth and death model, and the GTR+G substitution model, using 50 million MCMC generations and a burn-in of 10%.

Additionally, two cluster algorithms were used, the Barcode Index Number System (BIN) (Ratnasingham and Hebert, 2013) and Automatic Barcode Gap Discovery (ABGD) (Puillandre et al., 2012). The BIN was automatically determined in the BOLD Workbench, while the ABGD was performed using Kimura-2-parameter (K2P) distance and default values through the web interface2.

COI intraspecific and interspecific genetic distances were estimated using the K2P model implemented in Mega 6.0 (Tamura et al., 2013). These values were used to calculate the mean, minimum, and maximum values for intra- and inter-MOTU distances, and intra- and interspecific distances (nominal species). A genetic distance neighbor-joining (NJ) tree analysis was performed based on the K2P substitution model in Mega 6.0 (Tamura et al., 2013).

Results

The alignment of COI sequences resulted in 600 characters with 158 parsimony informative sites (included in the Supplementary Material). The GMYC analysis resulted in 18 MOTUs (Confidence interval: 16–18) (Table 2). The GMYC model was preferred over the null model (likelihood ratio = 73.49, P < 0.0001), indicating that GMYC results are reliable. The PTP analyses (maximum likelihood and Bayesian implementation) resulted in the same 18 MOTUs obtained in GMYC. The ABGD analysis found six partitions with 27 (P = 0.001) to 16 groups (P = 0.01), including a partition with the same 18 MOTUs (P = 0.005) obtained in the GMYC and PTP analyses. The BOLD system determined 16 BINs (Table 2), showing discordance with our MOTUs in only two BINs, AAB8569 [M. piavussu (Britski et al., 2012) and M. cf. piavussu lower Paraná] and AAD1729 [M. reinhardti (Lütken, 1875) and M. cf. reinhardti]. The clustering of the MOTUs obtained by the analyses is shown in Figure 3.

Table 2

Mean intra-Maximum intra-NNDistance to NN
MOTUs
Megaleporinus brinco00M. obtusidens6.78
M. conirostris00M. cf. conirostris3.99
M. cf. conirostris00M. conirostris3.99
M. elongatus0.040.17M. cf. obtusidens São Francisco2.74
M. garmani00M. obtusidens7.68
M. macrocephalus00M. macrocephalus Paraná1.86
M. macrocephalus Paraná00M. macrocephalus1.86
M. muyscorum00M. reinhardti11.6
M. cf. muyscorumM. trifasciatus7.48
M. obtusidens0.140.5M. cf. obtusidens São Francisco2.84
M. cf. obtusidens ParaguayM. piavussu Lower Paranáa2.9
M. cf. obtusidens São Francisco00M. elongatus2.74
M. piavussu0.060.17M. piavussu Lower Paraná0.67
M. piavussu Lower Paraná0.080.17M. piavussu0.67
M. reinhardti00M. cf. reinhardti0.67
M. cf. reinhardti00M. reinhardti0.67
M. trifasciatus00M. macrocephalus4.52
M. cf. trifasciatus00M. trifasciatus6.33
Nominal
Megaleporinus brinco00M. obtusidens6.78
M. conirostris2.133.99M. obtusidens5.6
M. elongatus0.040.17M. obtusidens2.74
M. garmani00M. obtusidens7.68
M. macrocephalus0.861.86M. trifasciatus4.52
M. muyscorum7.6615.31M. trifasciatus7.48
M. obtusidens1.946.72M. elongatus2.74
M. piavussu0.261.01M. obtusidens2.9
M. reinhardti0.310.7M. conirostris6.14
M. trifasciatus3.396.33M. macrocephalus4.52

Genetic K2P distances of Megaleporinus species.

The mean and the maximum of intra-group distances, the nearest neighbor (NN), and the minimum distance to the NN for MOTUs (ABGD, GMYC, and PTP) and Nominal species.

FIGURE 3

Only Megaleporinus brinco (Birindelli and Britski, 2013), Megaleporinus garmani (Borodin, 1929), and Megaleporinus elongatus (Valenciennes, 1850) showed correspondence between nominal species and MOTUs. Within six nominal species, a subdivision in two MOTUs was found, while Megaleporinus obtusidens (Valenciennes, 1837) was divided in three MOTUs (Table 2).

The mean of intra-MOTU and maximum intra-MOTU distances, the nearest neighbor (NN), and the minimum distance to the NN are shown in Table 2, for both GMYC MOTUs and nominal species.

The overall mean of intra-MOTU distances was 0.03%, the maximum intra-MOTU distance was 0.5% (M. obtusidens), and the mean of inter-MOTU distances was 9.19%. The lowest and highest values of inter-MOTU distances were 0.67 and 15.31%, respectively. Considering these values, there is a barcode gap that allowed identifying successfully all MOTUs using COI distance. In contrast, when only the nominal species were considered, the maximum intraspecific distance increased to 15.31% [M. muyscorum (Steindachner, 1900)], and, in addition, no barcode gap was found.

Discussion

Our hypothesis that some of the nominal species represent species complexes separated in different basins could not be rejected by DNA barcoding analysis, revealing taxonomic uncertainties and a hidden diversity within this recently described genus. The DNA barcode analyses identified 16 (ABGD and BIN) to 18 (ABGD, GMYC, and PTP) different MOTUs (Figure 3), with two new MOTUs (M. macrocephalus Paraná and M. cf. piavussu lower Paraná) not analyzed by Ramirez et al. (2017). This high number of MOTUs contrasts with the 10 nominal species recognized in the genus thus far, showing several potential target for cryptic species to be described, reinforcing the general idea that there is still a lot of undocumented diversity within the neotropical ichthyofauna (Reis et al., 2016). The difference between the number of MOTUs detected is due to the lower genetic distance value (0.67%) between two pairs of MOTUs: M. reinhardti and M. cf. reinhardti, separating the genetic lineages from São Francisco and Itapicuru, respectively, and between M. piavussu and M. cf. piavussu lower Paraná. These lower genetic distance values are likely due to a recent divergence between these MOTUs [<0.5 Ma for M. reinhardti and M. cf. reinhardti according to Ramirez et al. (2017)]. Of note, besides presenting an allopatric distribution, these MOTUs were also recovered by the monophyly criterion (Figure 3). MOTUs with recent origin have less time to accumulate genetic differences than species with ancient origin, hindering their identification. Despite this low genetic distance, the species delimitation methods could delimit these MOTUS, especially those based on phylogenetic trees (GMYC and PTP).

A key aspect implicit in the DNA barcoding analysis is the genetic distance threshold used to define MOTUs. COI distances of 1% (Hubert et al., 2008) to 2% (Pereira et al., 2013) have been claimed as threshold to fish DNA barcode analysis. However, such values were derived from comparative analyses among phylogenetically diverse groups. For instance, 2% was used to characterize DNA barcoding of a fish community of a given river (Pereira et al., 2013). However, when the DNA barcoding analyses have focused within a group of species closely related (e.g., a genus), lower threshold values have been reported (Carvalho et al., 2011; Pereira et al., 2011, 2013; Ramirez and Galetti, 2015). Particularly in Anostomidae, a lower threshold of 0.92% was reported to distinguish MOTUs within the genus Laemolyta (Ramirez and Galetti, 2015). Although most of the values obtained herein were above 2% (13 out of 18 MOTUs, Table 2), a maximum threshold of 0.67% for Megaleporinus was detected between the MOTUs obtained. It reinforces that lower genetic distance values might be obtained when intra-genus MOTUs are analyzed, mainly between recent divergent lineages.

Five nominal species, M. conirostris (Steindachner, 1875), M. macrocephalus (Garavello and Britski, 1988), M. muyscorum, M. obtusidens, and M. trifasciatus (Steindachner, 1876), showed high COI distance values (> 1.8%, Table 2) between individuals from different basins, indicating a scenario of potential allopatric speciation within these species.

In contrast to previous results (Avelino et al., 2015), evidence of local differentiation was not found here and all cryptic diversity correspond to inter-basin differentiation. Analyzing only two samples of M. reinhardti from the Três Marias (MG, Brazil) region (São Francisco basin), Avelino et al. (2015) reported an intraspecific distance of 3.8% between them, suggesting a local differentiation. Here we analyzed nine individuals, representing four different localities, including Três Marias region, and we found no genetic distance (0%) among them. Mitochondrial pseudogenes, sequencing errors, or misidentification could explain such discrepancies, and it would be more cautious to consider M. reinhardti from São Francisco as a single MOTU, as recovered here.

Similar discordance is observed for M. piavussu (upper Paraná). Avelino et al. (2015) included four samples from a single locality and reported a mean intraspecific distance of 2.8%. Our present data set for this species included 18 individuals obtained from six localities and showed a lower maximum intraspecific distance of 0.17%. It is strongly suggested that M. piavussu is also a single MOTU.

Incongruences were also observed within the nominal M. obtusidens. While four groups (A–D), showing 0.7–4.1% mean intraspecific distances, were previously reported (Avelino et al., 2015), we found three MOTUs showing 0–0.5% COI distances. The group D mentioned as part of M. obtusidens by Avelino et al. (2015), which included individuals caught downstream the Itaipú dam (Paraná basin), was recovered here as a sister group of M. piavussu, and was named M. cf. piavussu lower Paraná (Figure 3).

One particular aspect was highlighted in our results. Several individuals clustered in the M. macrocephalus clade were caught in different hydrographic basins, as Doce, São Francisco, Tocantins, and Paraná, outside of its original distribution in the Paraguay basin likely due to aquaculture releasing. Similar findings had already been described in the São Francisco basin (Carvalho et al., 2011). This species is a commercial important fish being extensively farmed throughout the Brazilian territory, and accidental or intentional releasing can occur (e.g., Langeani et al., 2007; Vieira, 2010). In such case, the use of DNA barcoding provides a rapid and accurate identification of this species and can be used in management and monitoring potential ecosystem disturbance caused by an invasive species.

In summary, the use of DNA barcoding points at the need for a taxonomic revision of this genus. A search for morphological traits able to support a taxonomic delimitation could be facilitated whether the MOTUs identified here are considered. A morphological trait showing a range of variation when searched within a given nominal species perhaps could be more informative if studied in each MOTU separately. In such case, our results would give an important contribution for the taxonomy of Megaleporinus facilitating the search for decisive taxonomic characters. This is the most complete analysis of DNA barcode in this recently described genus, and considering the economic value of this group, a precise species identification is quite desirable and fundamental for conservation of the whole biodiversity of this genus.

Statements

Author contributions

JR and PG designed the research. JR, DC, PA, PV, HO, MC-A, and JR-P collected data. JR performed the analyses. All authors contributed to the writing of the manuscript.

Funding

The authors thank Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for financial support (Universal 473474/2011-5 and 405309/2016-3 to PG, Universal 420255/2016-8 to JB and Rede BrBOL 564953/2010-5). JR received a fellowship grant from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP 2011/21836-4). Authors received productivity research grants from CNPq (304440/2009-4 to PG) and Fundação Araucária (641/2014 to JB).

Acknowledgments

We are grateful to C. Cramer, C. Doria, C. Nolorbe, D. Motta, H. Sanchez, J.C. Riofrio, and W. Troy for help to obtain part of the tissue samples and MMA/ICMBIO for collection authorization (32215-1). The authors thank the three reviewers for suggestion and comments which improved the manuscript.

Conflict of interest

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.

Supplementary material

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

References

  • 1

    AvelinoG. S.BritskiH. A.ForestiF.OliveiraC. (2015). Molecular identification of Leporinus from the south portion of South America.DNA Barcodes398109. 10.1515/dna-2015-0013

  • 2

    BensonD. A.CavanaughM.ClarkK.Karsch-MizrachiI.LipmanD. J.OstellJ.et al (2017). GenBank.Nucleic Acids Res.45D37D42. 10.1093/nar/gkw1070

  • 3

    BirindelliJ. L. O.BritskiH. A. (2013). Two new species of Leporinus (Characiformes: Anostomidae) from the Brazilian Amazon, and redescription of Leporinus striatus Kner 1858.J. Fish Biol.8311281160. 10.1111/jfb.12206

  • 4

    BirindelliJ. L. O.BritskiH. A.GaravelloJ. C. (2013). Two new species of Leporinus Agassiz (Characiformes?: Anostomidae) from eastern basins of Brazil, and redescription of L. melanopleura Günther.Neotrop. Ichthyol.11923. 10.1590/S1679-62252013000100002

  • 5

    BorodinN. A. (1929). Notes on some species and subespecies of the genus Leporinus Spix.Mem. Museum Comp. Zool.50269290.

  • 6

    BouckaertR.HeledJ.KühnertD.VaughanT.WuC.-H.XieD.et al (2014). BEAST 2: a software platform for bayesian evolutionary analysis.PLOS Comput. Biol.10:e1003537. 10.1371/journal.pcbi.1003537

  • 7

    BritskiH. A.BirindelliJ. L. O.GaravelloJ. C. (2012). A new species of Leporinus Agassiz, 1829 from the Upper Rio Paraná basin (Characiformes, Anostomidae) with redescription of L. elongatusValenciennes, 1850 and L. obtusidens (Valenciennes, 1837).Pap. Avulsos Zool.52441475.

  • 8

    BurnsM. D.FrableB. W.SidlauskasB. L. (2014). A new species of Leporinus (Characiformes: Anostomidae), from the Orinoco Basin, Venezuela.Copeia2014206214. 10.1643/CI-13-071

  • 9

    CarvalhoD. C.OliveiraD. A. A.PompeuP. S.LealC. G.OliveiraC.HannerR. (2011). Deep barcode divergence in Brazilian freshwater fishes: the case of the São Francisco River basin.Mitochondrial DNA228086. 10.3109/19401736.2011.588214

  • 10

    DarribaD.TaboadaG. L.DoalloR.PosadaD. (2012). jModelTest 2: more models, new heuristics and parallel computing.Nat. Methods9772. 10.1038/nmeth.2109

  • 11

    DíazJ.VillanovaG. V.BrancoliniF.del PazoF.PosnerV. M.GrimbergA.et al (2016). First DNA barcode reference library for the identification of South American freshwater fish from the Lower Paraná river.PLOS ONE11:e0157419. 10.1371/journal.pone.0157419

  • 12

    Frantine-SilvaW.SofiaS. H.OrsiM. L.AlmeidaF. S. (2015). DNA barcoding of freshwater ichthyoplankton in the Neotropics as a tool for ecological monitoring.Mol. Ecol. Resour.1512261237. 10.1111/1755-0998.12385

  • 13

    GalettiP. M.Jr.CesarA. C. G.VenereP. C. (1991). Heterochromatin and NORs variability in Leporinus fish (Anostomidae, Characiformes).Caryologia44287292. 10.1080/00087114.1991.10797193

  • 14

    GalettiP. M.Jr.LimaN. R. W.VenereP. C. (1995). A monophyletic ZW sex chromosome system in Leporinus (Anostomidae, Characiformes).Cytologia (Tokyo).60375382. 10.1508/cytologia.60.375

  • 15

    GalettiPMJrForestiF.BertolloL. A.MoreiraFilho O (1981). Heteromorphic sex chromosomes in three species of the genus Leporinus (Pisces, Anostomidae).Cytogenet. Genome Res.29138142. 10.1159/000131562

  • 16

    GaravelloJ. C.BritskiH. A. (1988). Leporinus macrocephalus sp. n. da bacia do rio Paraguai (Ostariophysi, Anostomidae).Naturalia136774.

  • 17

    GaravelloJ. C.BritskiH. A. (2003). “Family Anostomidae,” inCheck List of the Freshwater Fishes of South and Central AmericaedsReisR. E.KullanderS. O.JrFerrarisC. J.. (Porto Alegre: EDIPUCRS) 7184.

  • 18

    GomesL. C.PessaliT. C.SalesN. G.PompeuP. S.CarvalhoD. C. (2015). Integrative taxonomy detects cryptic and overlooked fish species in a neotropical river basin.Genetica143581588. 10.1007/s10709-015-9856-z

  • 19

    HallT. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.Nucleic Acids Symp. Ser.419598.

  • 20

    HubertN.HannerR.HolmE.MandrakN. E.TaylorE.BurridgeM.et al (2008). Identifying Canadian freshwater fishes through DNA barcodes.PLOS ONE3:e2490. 10.1371/journal.pone.0002490

  • 21

    LangeaniF.Corrêa e CastroR. M.OyakawaO. T.ShibattaO. A.PavanelliC. S.CasattiL. (2007). Diversidade da ictiofauna do Alto Rio Paraná: composição atual e perspectivas futuras.Biota Neotrop.7181197. 10.1590/S1676-06032007000300020

  • 22

    LütkenC. R. (1875). Velhas Flodens Fiske. Et bidrag til Brasiliens Ichthyologi. Elfter Prof. J. Reinhardt Indsamlinger og Optegnelser.Danske Vidensk. Selsk. Skr. Kjøbenhavn12122254.

  • 23

    MachadoC. D. B.IshizukaT. K.FreitasP. D.De ValiatiV. H.GalettiP. M. (2016). DNA barcoding reveals taxonomic uncertainty in Salminus (Characiformes).Syst. Biodivers.15372382. 10.1080/14772000.2016.1254390

  • 24

    PereiraL. H. G.HannerR.ForestiF.OliveiraC. (2013). Can DNA barcoding accurately discriminate megadiverse Neotropical freshwater fish fauna?BMC Genet.14:20. 10.1186/1471-2156-14-20

  • 25

    PereiraL. H. G.MaiaG. M. G.HannerR.ForestiF.OliveiraC. (2011). DNA barcodes discriminate freshwater fishes from the Paraíba do Sul River Basin, São Paulo, Brazil.Mitochondrial DNA22(Suppl. 1)7179. 10.3109/19401736.2010.532213

  • 26

    PonsJ.BarracloughT.Gomez-ZuritaJ.CardosoA.DuranD.HazellS.et al (2006). Sequence-based species delimitation for the DNA taxonomy of undescribed insects.Syst. Biol.55595609. 10.1080/10635150600852011

  • 27

    PuillandreN.LambertA.BrouilletS.AchazG. (2012). ABGD, automatic barcode gap discovery for primary species delimitation.Mol. Ecol.2118641877. 10.1111/j.1365-294X.2011.05239.x

  • 28

    R Core Team (2017). R: A Language and Environment for Statistical Computing.Vienna: R Foundation for Statistical Computing.

  • 29

    RamirezJ. L.BirindelliJ. L. O.GalettiP. M. (2017). A new genus of Anostomidae (Ostariophysi: Characiformes): diversity, phylogeny and biogeography based on cytogenetic, molecular and morphological data.Mol. Phylogenet. Evol.107308323. 10.1016/j.ympev.2016.11.012

  • 30

    RamirezJ. L.Carvalho-CostaL. F.VenereP. C.CarvalhoD. C.TroyW. P.GalettiP. M. (2016). Testing monophyly of the freshwater fish Leporinus (Characiformes, Anostomidae) through molecular analysis.J. Fish Biol.8812041214. 10.1111/jfb.12906

  • 31

    RamirezJ. L.GalettiP. M.Jr. (2015). DNA barcode and evolutionary relationship within Laemolyta Cope 1872 (Characiformes: Anostomidae) through molecular analyses.Mol. Phylogenet. Evol.937782. 10.1016/j.ympev.2015.07.021

  • 32

    RatnasinghamS.HebertP. D. N. (2013). A DNA-based registry for all animal species: the Barcode Index Number (BIN) System.PLOS ONE8:e66213. 10.1371/journal.pone.0066213

  • 33

    ReisR. E.AlbertJ. S.Di DarioF.MincaroneM. M. M.PetryP. L.RochaL. R. (2016). Fish biodiversity and conservation in South America.J. Fish Biol.891247. 10.1111/jfb.13016

  • 34

    ReisR. E.KullanderS. O.FerrarisC. J.Jr. (2003). Check List of the Freshwater Fishes of South and Central America.Porto Alegre: EDIPUCRS.

  • 35

    SambrookJ.FritishE. F.ManiatisT. (1989). Molecular Cloning: A Laboratory Manual.Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press.

  • 36

    SidlauskasB. L.VariR. P. (2008). Phylogenetic relationships within the South American fish family Anostomidae (Teleostei, Ostariophysi,Characiformes).Zool. J. Linn. Soc.15470210. 10.1111/j.1096-3642.2008.00407.x

  • 37

    SteindachnerF. (1875). Die Süsswasserfische des südöstlichen Brasilien (II).Sitzungsber. Akad. Wiss. Wien71211245.

  • 38

    SteindachnerF. (1876). Ichthyologische Beiträge (V).Sitzungsber. Akad. Wiss. Wien7449240.

  • 39

    SteindachnerF. (1900). Erstattungen eines vorlaüfigen Berichtes über einige von Ihrer königlichen Hoheit Frau Prinzessin Therese von Bayeren während einer Reise nach Südamerika 1898 gesammelte neue Fischarten.Anz. Akad. Wiss. Wien373.

  • 40

    TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: Molecular Evolutionary Genetics Analysis version 6.0.Mol. Biol. Evol.3027252729. 10.1093/molbev/mst197

  • 41

    ValenciennesA. (1837). “Poissons [plates],” inVoyage Dans l’Amérique Méridionaleed.d’OrbignyA. (Paris: Pitois-Levrault) 1834-42.

  • 42

    ValenciennesM. A. (1850). “Suite du livre vingt-deuxième. Suite de la famille des Salmonoïdes,” inHistoire Naturelle des Poissons. Tome Vingt-DeuxièmeedsCuvierM.ValenciennesM. A. (Malden, MA: Blackwell Science) 191.

  • 43

    VieiraF. (2010). Distribuição, impactos ambientais e conservação da fauna de peixes da bacia do rio Doce.MG Biota2522.

  • 44

    ZhangJ.KapliP.PavlidisP.StamatakisA. (2013). A general species delimitation method with applications to phylogenetic placements.Bioinformatics2928692876. 10.1093/bioinformatics/btt499

Summary

Keywords

cryptic species, freshwater fishes, allopatric speciation, South American basins, cytochrome oxidase subunit I

Citation

Ramirez JL, Birindelli JL, Carvalho DC, Affonso PRAM, Venere PC, Ortega H, Carrillo-Avila M, Rodríguez-Pulido JA and Galetti Jr. PM (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

Received

09 August 2017

Accepted

27 September 2017

Published

12 October 2017

Volume

8 - 2017

Edited by

Rodrigo A. Torres, Federal University of Pernambuco, Brazil

Reviewed by

Fábio Fernandes Roxo, Universidade Estadual Paulista Júlio de Mesquita Filho (UNESP), Brazil; Henrik R. Nilsson, University of Gothenburg, Sweden; Uedson Pereira Jacobina, Federal University of Alagoas, Brazil

Updates

Copyright

*Correspondence: Jorge L. Ramirez,

This article was submitted to Evolutionary and Population Genetics, a section of the journal Frontiers in Genetics

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics