Original Research ARTICLE
Genetic Variation of the Endangered Neotropical Catfish Steindachneridion scriptum (Siluriformes: Pimelodidae)
- 1Laboratório de Biologia e Cultivo de Peixes de Água Doce, Departamento de Aquicultura, Universidade Federal de Santa Catarina, Florianópolis, Brazil
- 2Programa de Pós-Graduação em Aquicultura, Centro de Ciências Agrárias, Universidade Federal de Santa Catarina, Florianópolis, Brazil
Steindachneridion scriptum is an important species as a resource for fisheries and aquaculture; it is currently threatened and has a reduced occurrence in South America. The damming of rivers, overfishing, and contamination of freshwater environments are the main impacts on the maintenance of this species. We accessed the genetic diversity and structure of S. scriptum using the DNA barcode and control region (D-loop) sequences of 43 individuals from the Upper Uruguay River Basin (UUR) and 10 sequences from the Upper Paraná River Basin (UPR), which were obtained from GenBank. S. scriptum from the UUR and the UPR were assigned in two distinct molecular operational taxonomic units (MOTUs) with higher inter-specific K2P distance than the optimum threshold (OT = 0.0079). The COI Intra-MOTU distances of S. scriptum specimens from the UUR ranged from 0.0000 to 0.0100. The control region indicated a high number of haplotypes and low nucleotide diversity, compatible with a new population in recent expansion process. Genetic structure was observed, with high differentiation between UUR and UPR basins, identified by BAPS, haplotype network, AMOVA (FST = 0.78, p < 0.05) and Mantel test. S. scriptum from the UUR showed a slight differentiation (FST = 0.068, p < 0.05), but not isolation-by-distance. Negative values of Tajima’s D and Fu’s Fs suggest recent demographic oscillations. The Bayesian skyline plot analysis indicated possible population expansion from beginning 2,500 years ago and a recent reduction in the population size. Low nucleotide diversity, spatial population structure, and the reduction of effective population size should be considered for the planning of strategies aimed at the conservation and rehabilitation of this important fisheries resource.
Freshwater ecosystems are among the most endangered ecosystems (Dudgeon et al., 2006). Habitat degradation, hydrologic alterations, habitat fragmentation, sediment deposition, and overfishing are the principal causes of declines and extinctions of freshwater fishes (Dudgeon et al., 2006; Agostinho et al., 2008; Helfman, 2008; Hoeinghaus et al., 2009). In addition, species with geographically restricted distribution are more susceptible to erosion of the genetic diversity due to habitat fragmentation (Vrijenhoek et al., 1985). Fishery resources are an integral part of most societies and make important contributions to economic and social health and well-being in many countries and areas (FAO, 2002). The understanding of the genetic diversity and structure of wild populations of fish species are important to the regulation of fisheries and conservation management strategies (Carvalho and Hauser, 1994; Iervolino et al., 2010).
Steindachneridion scriptum (Miranda Ribeiro, 1918) is a large catfish belonging to the family Pimelodidae. This potamodromous fish species presents restricted distribution in the Upper Uruguay River (UUR) and Upper Paraná River (UPR) basins (Marques et al., 2002). S. scriptum is an important fishing resource to the riverine fishermen (Schork et al., 2013); nonetheless, it was recently classified as an endangered species according to the Chico Mendes Institute for Biodiversity Conservation (ICMBIO, 2014). Human activities (e.g., damming of rivers, illegal fishing, and industrial waste) are the main threats to S. scriptum in the UUR (Fontana et al., 2003) and are the principle reason for population reduction in this basin (Beux and Zaniboni-Filho, 2008).
The correct management of fish stocks depends on the precise identification of the target species, which may present very similar morphological characteristics to other species. The latest revision of the genera Steindachneridion recognizes six valid species (Garavello, 2005). Using museum specimen, Garavello (2005) suggests the presence of S. punctatum in the UUR and UPR with few characteristics that differ from S. scriptum. Despite the description of the two species of the genus Steindachneridion, ichthyofauna studies never report the presence of S. punctatum in the UUR (Schork et al., 2012, 2013). In addition to the correct ichthyofauna and management of fish stocks, uncertainties of taxonomic identification may be a problem for stock management (e.g., formation of in vivo and in vitro banks, restocking programs) and adequate fisheries control.
Taxonomic uncertainties are common in fishes (Pereira et al., 2013) and can be investigated using DNA barcode methodology, which permits the unambiguous identification of the majority of fish species (Ward, 2009). Particularly for endangered species, prior knowledge of the distribution of genetic variability within and among natural populations as well as the implementation of an efficient management plan based on genetic features are important measures for its maintenance and recovery (Cross, 2000). The distribution of the genetic variation within and between populations can be assessed using the mitochondrial control region (Sivasundar et al., 2001; Iervolino et al., 2010; Ochoa et al., 2015). Given these findings, we tested the null hypothesis that S. scriptum from UUR represent a single molecular operational taxonomic unit (MOTU). Posteriorly, we investigate the genetic diversity and population structure of S. scriptum from the UUR and UPR basins using the mitochondrial control region.
Materials and Methods
Study Area and Sampling
The Uruguay River (Uru) originates in Brazilian territory (in the Serra Geral Mountains) in Southern Brazil, together with the Paraná and Paraguay Rivers form the La Plata Basin (Zaniboni-Filho and Schulz, 2003). Samples of S. scriptum were collected by scientific fishing and local fishermen between 2006 and 2015, with authorization of the Brazilian Institute of the Environment and Renewable Natural Resources (IBAMA; protocol number: 02026.005762/2004-71). A total of 19 individuals from the Uru River and 24 from the Canoas River (Can), UUR basin, were sampled (Figure 1). Tissues were preserved in 95% ethanol until extraction. This research was conducted under Animal Care Protocol PP00788 of the Federal University of Santa Catarina (UFSC).
FIGURE 1. Map showing the locations of sampling sites of the Steindachneridion scriptum in the Upper Uruguay River Basin. Black dots indicate sample points.
DNA Extraction and Amplification of the Mitochondrial Fragments
Total DNA was obtained with fin clips following a salt extraction method (Aljanabi and Martinez, 1997). For the DNA barcode analyses, a fragment of 652 bp of Cytochrome Oxidase subunit I (COI) was amplified through polymerase chain reaction (PCR), using primers FishF1/FishR1 (Ward et al., 2005) and following Bellafronte et al. (2013). The partial amplification of the mitochondrial control region (D-loop) was performed using primers FTTP-L and DLR1-H according to Huergo et al. (2011). The PCR products were checked for amplification using gel electrophoresis with 1% agarose gels purified using PEG 20% (Lis, 1980). Sequencing reactions were performed using BigDye TM Terminator v 3.1 (Applied Biosystems), and the PCR products were sequenced for both strands in ABI 3500XL (Applied Biosystems).
DNA sequences from each individual of both genes were edited using Geneious 5.4.4 (Wu and Drummond, 2011) to generate a consensus sequence. For the DNA barcode analysis, we combined the COI reference sequences of S. scriptum (access FUPR686-09, PDCAP027-14, and PDCAP028-14) from the UPR basin, S. parahybae (access FPSR293-10–FPSR297-10) from the Paraíba do Sul Basin and one specimen of Pseudoplatystoma corruscans and one specimen of Zungaro jahu to root our phylogenetic analyses. All sequences available in Barcode of Life Data System (BOLD). Intra- and inter-specific genetic distances based on the Kimura 2-parameter (K2P) evolution model were calculated using Mega 6 (Tamura et al., 2013).
We used the phylogenetic General Mixed Yule Coalescent (GMYC) approach based on single-locus data that is a relative robust tool for species delimitation (Pons et al., 2006; Fujisawa and Barraclough, 2013). The ultrametric tree was generated in BEAST v.2.2.1 (Bouckaert et al., 2014), with the substitution model calculated in the JModelTest 2.1.4 (HKY+G; Darriba et al., 2012), using relaxed molecular clock with a lognormal distribution and birth–death model. Three independent runs were carried out with 20 million generations each. Posteriorly, the runs were combined using the LogCombiner v.1.8.3 (Drummond et al., 2012), with a burn-in of 25. Data mixing and effective sample size (ESS) were verified in Tracer v1.5. GMYC was carried out in Species Limits by Threshold Statistics (SPLITs; Monaghan et al., 2009) with RStudio1, using the unique threshold method to detect the transition point between intra- and inter-specific relationships.
In addition to the standard threshold adopted to Neotropical fishes (Pereira et al., 2013), we calculated an optimum threshold (OT; Collins et al., 2012) directly from all dataset, using the local minima function in the R package SPIDER (SPecies IDentity and Evolution in R; Brown et al., 2012). The OT value was used to define the MOTUs using the software jMOTU (Jones et al., 2011). The graphical representation of the MOTUs was performed by a neighbor-joining analysis (NJ) using the K2P model with Mega 6.6 (Tamura et al., 2013). The support of the clades was tested by the bootstrap method with 10,000 pseudo-replicates.
DNA D-loop sequences were combined with 16 S. scriptum sequences downloaded from GenBank, of specimens collected between 1995 and 2002 that corresponded to the Uru River (access EU930029.1–EU930038.1) and specimens from Tibaji River (UPR) (access EU930039.1–EU930044.1). The overall genetic diversity was estimated using the following DnaSP software (Rozas et al., 2003) parameters: nucleotide diversity (π) (Nei, 1987), haplotype diversity (Hd) (Nei and Tajima, 1981), and number of polymorphic sites (S). Genetic diversity within and between sample sites was hierarchically tested by Analysis of Molecular Variance (AMOVA) (Excoffier et al., 1992) with 10,000 permutations to test the pairwise population comparison (FST) using Arlequin 3.1 (Excoffier and Lischer, 2010). Spatial genetic structure was inferred using Bayesian Analysis of Population Structure 6.0 (BAPS) software (Corander et al., 2013). First, BAPS was run with 10 replicates for every level of k (1–6) without origin information (“clustering of individuals”) and subsequently using “clustering of groups of individuals.” We used Mantel tests as implemented in the Alleles In Space (AIS) 1.0 (Miller, 2005) to test for a correlation of geographic stream distance and genetic distance (isolation-by-distance; IBD) expressed as FST, with 10,000 permutations to assess significance. Tajima’s (1989) D, Fu’s (1997) Fs, and mismatch distributions were estimated with DnaSP (Rozas et al., 2003). A median-joining haplotype network was generated through PopART (Bandelt et al., 1999). Demographic history was investigated using Bayesian Skyline Plot (BSP) in BEAST 2.1.3 (Drummond et al., 2012) with the evolutionary model obtained in the Jmodeltest program. The graphic was generated in Tracer 1.66 (Rambaut et al., 2014).
The COI sequences were deposited in BOLD systems (accession UUR001-17–UUR052-17), and DNA D-loop sequences were uploaded in GenBank (accession MF045370–MF045412). The voucher corresponding to S. scriptum from the UUR basin was deposited in the Zoology Museum of the Universidade Estadual de Londrina MZUEL 15569 (Type locality: Itaqui, RS – Uru Basin).
Mitochondrial DNA Barcoding
The consensus alignment of 43 COI sequences were obtained from samples identified morphologically as S. scriptum, resulting in a total length of 611 bp, with 7 polymorphic sites and 6 haplotypes defined. The inter-specific nucleotide frequencies were 26.02% of Cytosine, 28.97% of Thymine, 27.19% of Adenine, and 17.82% of Guanine. No stop codons, insertions, or deletions were observed in the COI sequences, indicating that they represent fragments of functional mitochondrial genes and not nuclear mitochondrial pseudo-genes (Numts). Considering P. corruscans and Z. jahu as the outgroup, the maximum likelihood for the GMYC model was significantly superior (L = 628.3242) to the likelihood of the null model (Lo = 614.7156, p < 0.0001). The single-threshold GMYC model suggested the presence two clusters (confidence interval 2–13) of four ML entities (GMYC ‘species,’ named ‘MOTUs’ herein) with a confidence interval of 4–17 (S. scriptum, S. parahybae, and two outgroups).
The optimal threshold calculated for all S. scriptum sequences used in this study was OT = 0.0079 (0.79%) of divergence. From the set of the genus Steindachneridion sequences available on BOLD, three MOTUs were identified using the OT value and software jMOTU: S. scriptum from the UUR, S. scriptum from the UPR, and S. parahybae. The COI inter-MOTU between S. scriptum from the UUR and S. scriptum from the UPR showed values larger than the OT (mean 0.012, minimum 0.010), while intra-MOTU values for fishes from both basins were lower than the average OT (UUR = 0.000 to 0.010; UPR = 0.000) (Table 1). The minimum inter-specific distances between S. scriptum MOTUs (UUR and UPR) and S. parahybae were 0.100 (10%) and 0.090 (9%), respectively.
TABLE 1. Intra-MOTU distances (in bold) and inter-MOTU genetic distances using the COI gene and K2p model.
The NJ-K2P (Supplementary Figure SI 1) and Bayesian Inference topologies (Figure 2) were clustered in two principal clusters, corresponding to S. scriptum and S. parahybae species. The S. scriptum clade appeared divided into two well-supported sub-clades formed by UUR and UPR individuals, in both NJ and BI phylogenetic trees. Based on these results, all specimens from Uru were composed of a single MOTU, named here as S. scriptum.
FIGURE 2. Phylogenetic tree of the genus Steindachneridion based on the mitochondrial gene COI inferred with BEAST. Values before nodes are Bayesian posterior probabilities. Codes FUPR, PDCAP, and FPSR were downloaded from the BOLD system. Pseudoplatystoma corruscans and Z. jahu are outgroups. Vertical lines in Blue = jMOTU MOTUs, and vertical lines in Green indicate GMYC MOTUs.
Mitochondrial DNA D-Loop
The final alignment size of the 59 D-loop consensus sequences of S. scriptum specimens from the UUR and UPR (43 newly sequenced, 16 downloaded from GenBank) was 865 bp. A total of 56 variable sites were found in the region defining a total of 36 haplotypes. Overall, 30 haplotypes identified in the analyses were unique and exclusive. Of these, 20 were from the Uru, 4 from Can, and 6 from UPR basin (Tibagi River). The most common haplotype was H8, which was recorded 12 times and was shared by samples from the Uru and Can rivers, and haplotypes H4 and H27, both with three records, were exclusive to the Uru and Can rivers, respectively. The haplotype network (Figure 3) revealed a high degree of similarity between the specimens from the Uru and Can rivers, even though with slight differences in the haplotype frequencies and exclusive haplotypes from each river. These specimens were also differentiated from the samples of UPR by 12 mutations.
FIGURE 3. Median-joining network of Steindachneridion scriptum, based on haplotypes of mtDNA control region. The colors indicate locality according to the legend, size of the circles illustrates the number of identical haplotypes, and small black circles, hypothetical ancestors or unsampled haplotypes. Hatch marks represent the number of mutations by which haplotypes differ.
The average nucleotide frequencies were 34.38% of Adenine, 33.17% of Thymine, 12.19% of Guanine, and 20.26% of Cytosine. Genetic variability, expressed as Hd) and nucleotide diversity (π), was higher in S. scriptum from the Uru (Hd = 0.959/π = 0.007) in comparison with the Can River (Hd = 0.837/π = 0.004) (Table 2). In addition, the samples from UPR have high diversity indices, in comparison with samples from UUR (N = 6; Nh = 6; Hd = 1.0000, π = 0.00698; D = -0.88901, p = 0.24002; Fs = -1.81313, p = 0.07768).
TABLE 2. Genetic diversity of Steindachneridion scriptum from the Uruguay and Canoas rivers (Upper Uruguay River Basin) estimated using D-loop control region.
The patterns of genetic variability found within and between populations in the AMOVA were based on the two principal clusters: Uruguay vs. Paraná Basins and Uru vs. Can rivers (Table 3). When the populations were considered as two basin groups, the AMOVA among groups was 78.17% and FST value was highly significant (FST = 0.781; p = 0.000). Genetic divergence between individuals from the Uruguay Basin (Uru vs. Can rivers) was low but significant (FST = 0.0682, p = 0.00475). The population groupings generated by the BAPS, without origin information of the samples, revealed the existence of three clusters (K = 3, Supplementary Figure SI 2), with slight differentiation between Can and Uru individuals. On the other hand, analysis with the individuals identified by sample group indicated two clusters (K = 2, Figure 4), corresponding to Uruguay and Paraná Basins. The IBD analysis showed a significant positive correlation (r2 = 0.85, p < 0.001) between the geographical distance and corresponding FST for S. scriptum from UUR and UPR (Supplementary Figure SI 3). On the other hand, FST values plotted over distance no reveal patter of isolation by distance for S. scriptum from UUR (r = -0.052, p = 0.682).
TABLE 3. Hierarchical AMOVA analysis and FST values for Steindachneridion scriptum according to their geographical location, estimated using D-loop control region.
FIGURE 4. Estimate of the probable groups of populations produced by the BAPS between Steindachneridion scriptum from Upper Uruguay and Upper Parana Basins, assigned to two clusters (K = 2).
Tajima’s (1989) D-neutrality tests, applied to detect evidence of strong selective pressures, and Fu’s (1997) Fs-tests, used specifically to detect population expansion, revealed significant negative values for all individuals from the UUR (D = -1.907, p < 0.05; Fs = -19.246, p < 0.01; Table 3). Non-significant negative values were estimated for D and Fs indexes for specimens from the Can River and UPR, whereas fish sampled in the Uru showed significant negative values estimated for D and Fs indexes. The BSP analysis (Supplementary Figure SI 4), used to explore previous demographic signals of S. scriptum, indicated early demographic expansion approximately 2,500 years ago as well as a fairly recent population reduction.
The DNA barcode confirmed the identification of all the individuals of S. scriptum from UUR as a single MOTU. The different methods, jMOTU and GMYC, were congruent in delimiting S. scriptum and S. parahybae. However, within S. scriptum, jMOTU methodology identified two distinct MOTUs between UUR basin and UPR basin, while GMYC only one, despite strongly supported clade. Inter-MOTU divergence between S. scriptum from the UUR and the UPR was higher than the OT and the mean intra-specific divergence found for freshwater fish (0.3%) (Ward et al., 2005; Lara et al., 2010; Pereira et al., 2013). The MOTUs do not necessarily represent species (Blaxter et al., 2005) but can indicate molecular entities (Casiraghi et al., 2010). The two MOTUs estimated between S. scriptum from the UUR and S. scriptum from the UPR could be explained by geographic isolation between watersheds that occurred during the Miocene epoch (between 5 and 24 million years ago) when the UUR and the UPR became isolated (Albert and Reis, 2011). Based on the results, S. scriptum from these two hydrographic systems are most likely in the process of incipient allopatric speciation since the genetic structuring in fish is in fact often evidenced and influenced principally by geological, ecological, and behavioral factors (Allan and Flecker, 1993). Congeners S. scriptum and S. parahybae species showed mean inter-MOTUs 10 times greater than the OT, indicating the existence of the barcode gap (Hebert et al., 2003, 2004) that allows us to assign an unknown Steindachneridion specimen to its species using a genetic distance criterion with an insignificant error rate.
Although they belong to the same watershed, individuals from the Can River and the Uru River showed a slight genetic differentiation probably due to the topography of the region and the interaction between the species’ biology and environmental characteristics (e.g., the Can River is located at a higher altitude with a lower water temperature than the Uru River). The haplotype network indicated a greater genetic similarity between the specimens from the Can and Uru, whereas is possible to observe differences in the haplotype frequencies, and exclusive haplotypes for each river. Recent studies with potamodromous fish species reported genetic structure in hydrographic systems without apparent physical barriers, resulting from behaviors related to IBD (Hardy and Vekemans, 1999; Primmer et al., 2006; Han et al., 2010), homing (Windle and Rose, 2005; Batista and Alves-Gomes, 2006; Neville et al., 2006), and isolation-by-time (IBT) (Hendry and Day, 2005; Braga-Silva and Galetti, 2016; Ribolli et al., 2017).
Nucleotide diversity of S. scriptum from the UUR was low in comparison with the values found in neotropical freshwater fish (π = 1.5%) (Batista and Alves-Gomes, 2006; Iervolino et al., 2010; Ashikaga et al., 2015). The high haplotype diversity and low nucleotide diversity seem quite compatible with a new population in recent expansion process, similar to what is shown in BSP analyses. This pattern may be a signature of such expansion that is long enough to examine a change in the haplotypes resulted from the mutation, but is not long enough to accumulate large differences between sequences (Avise, 2000).
Individuals from the Uru River were genetically more diverse than fishes from the Can River, indicating that the main channel of the Uru River allows the meeting of individuals of different areas (or tributaries), favoring the maintenance of the highest level of genetic diversity. Although it is a relevant fishing resource, useful molecular markers such as microsatellites, extensively employed in fish genetic studies, are still not developed for any Steindachneridion species, and the knowledge about genetic characteristics of the genus is incipient. RAPD markers indicated low genetic diversity for S. scriptum from the UUR, as reported by Ramella et al. (2006), as well as for S. melanodermatum from the Iguaçu River Basin (Matoso et al., 2011). Low genetic diversity may indicate recent or historic reduction of this diversity; however, some endangered populations may have a historic maintenance of small effective population sizes (Matocq and Villablanca, 2001). In this way, the low diversity detected in S. scriptum from the UUR can be attributed to the following: (1) population history of this species: a Bayesian skyline plot analysis revealed a subtle increase in effective population size over time and demographic swings in the recent past, with a notable increase in effective population size between 1,000 and 2,500 years ago and a subsequent reduction in the effective size of females; and (2) evolutionary history of the species: according to Garavello (2005) and Swarça et al. (2005), S. scriptum notably maintain conserved morphological and cytogenetic patterns.
In general, concerns and actions of conservation are more related to the perception of the disappearance of a given species than of genetic diversity reduction (Frankham, 2005). Therefore, given the low genetic diversity associated with the current scenario of fragmentation of the UUR basin and the population reduction of S. scriptum in some stretches of the Itá and Machadinho reservoirs (Schork et al., 2012, 2013), this study highlights the necessity of mitigation measures and more intense monitoring of illegal fishing to avoid the collapse of this important fishing resource. In addition, our results were congruent, identifying great differentiation between individuals from the UPR and UUR Basins. Further studies with a larger number of samples morphological analyze to better define the taxonomic status of endangered S. scriptum.
RP made substantial contributions to the design of the work and acquisition, analysis, and interpretation of data for the work, and drafted it until the approval for publication of the content. JR made substantial contributions to the conception and design of the work, analysis and interpretation of the data, and drafted and revised it critically for important intellectual content. EZ-F made substantial contributions to the conception of the work, revised it critically for important intellectual content, and provided approval for publication of the content. The authors agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
This work was supported by the Engie Tractebel Energia, Consórcio Machadinho e Consórcio Itá. RP thanks Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES); JR acknowledges the Programa Nacional de Pós-Doutorado (PNPD/CAPES); EZ-F thanks Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, Grant No. 302860/2014-2).
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.
We are grateful to Laboratório de Biologia e Cultivo de Peixes de Água Doce (LAPAD) of Universidade Federal de Santa Catarina (UFSC), Pedro Iaczinki, and local fishermen for help with fish collections. We thank Dra. Carolina Machado, Leonardo Porto Ferreira, and Ueslei Lopes for their contributions. We also thank the two reviewers for their valuable contributions on improving the manuscript and for the suggestions. Research was conducted under Animal Care Protocol PP00788.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00048/full#supplementary-material
Agostinho, A. A., Pelicice, F. M., and Gomes, L. C. (2008). Dams and the fish fauna of the Neotropical region: impacts and management related to diversity and fisheries. Braz. J. Biol. 68, 1119–1132. doi: 10.1590/S1519-69842008000500019
Ashikaga, F. Y., Orsi, M. L., Oliveira, C., and Senhorini, J. A. (2015). The endangered species Brycon orbignyanus: genetic analysis and definition of priority areas for conservation. Environ. Biol. Fishes 98, 1845–1855. doi: 10.1007/s10641-015-0402-8
Batista, J. S., and Alves-Gomes, J. A. (2006). Phylogeography of Brachyplatystoma rousseauxii (Siluriformes - Pimelodidae) in the Amazon Basin offers preliminary evidence for the first case of “homing” for an Amazonian migratory catfish. Genet. Mol. Res. 5, 723–740.
Bellafronte, E., Mariguela, T. C., Henrique, L., Pereira, L. H. G., Oliveira, C., and Moreira-Filho, O. (2013). DNA barcode of Parodontidae species from the La Plata river basin - applying new data to clarify taxonomic problems. Neotrop. Ichthyol. 11, 497–506. doi: 10.1590/S1679-62252013000300003
Beux, L., and Zaniboni-Filho, E. (2008). “Produção pesqueira do reservatório de Itá. Alto Rio Uruguai,” in Reservatório de Itá: Estudos Ambientais. Desenvolvimento de Tecnologias de Cultivo e Conservação da Ictiofauna, ed. E. Z. Filho (Florianópolis: UFSC), 87–108.
Blaxter, M., Mann, J., Chapman, T., Thomas, F., Whitton, C., Floyd, R., et al. (2005). Defining operational taxonomic units using DNA barcode data. Philos. Trans. R. Soc. Lond. B Biol. Sci. 360, 1935–1943. doi: 10.1098/rstb.2005.1725
Bouckaert, R., Heled, J., Kühnert, D., Vaughan, T., Wu, C. H., Xie, D., et al. (2014). BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 10:e1003537. doi: 10.1371/journal.pcbi.1003537
Braga-Silva, A., and Galetti, P. M. (2016). Evidence of isolation by time in freshwater migratory fish Prochilodus costatus (Characiformes, Prochilodontidae). Hydrobiologia 765, 159–167. doi: 10.1007/s10750-015-2409-8
Brown, S. D. J., Collins, R. A., Boyer, S., Lefort, M. C., Malumbres-Olarte, J., Vink, C. J., et al. (2012). Spider: an R package for the analysis of species identity and evolution, with particular reference to DNA barcoding. Mol. Ecol. Resour. 12, 562–565. doi: 10.1111/j.1755-0998.2011.03108.x
Casiraghi, M., Labra, M., Ferri, E., Galimberti, A., and de Mattia, F. (2010). DNA barcoding: a six-question tour to improve users’ awareness about the method. Brief. Bioinform. 11, 440–453. doi: 10.1093/bib/bbq003
Collins, R. A., Boykin, L. M., Cruickshank, R. H., and Armstrong, K. F. (2012). Barcoding’s next top model: an evaluation of nucleotide substitution models for specimen identification. Methods Ecol. Evol. 3, 457–465. doi: 10.1111/j.2041-210X.2011.00176.x
Dudgeon, D., Arthington, A. H., Gessner, M. O., Kawabata, Z. I., Knowler, D. J., Lévêque, C., et al. (2006). Freshwater biodiversity: importance, threats, status and conservation challenges. Biol. Rev. Camb. Philos. Soc. 81, 163–182. doi: 10.1017/S1464793105006950
Excoffier, L., and Lischer, H. E. L. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847.x
Excoffier, L., Smouse, P. E., and Quattro, J. M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131, 479–491.
Fujisawa, T., and Barraclough, T. G. (2013). Delimiting species using single-locus data and the Generalized Mixed Yule Coalescent approach: a revised method and evaluation on simulated data sets. Syst. Biol. 62, 707–724. doi: 10.1093/sysbio/syt033
Han, Y., Hung, C., Liao, Y., and Tzeng, W. (2010). Population genetic structure of the Japanese eel Anguilla japonica: panmixia at spatial and temporal scales. Mar. Ecol. Prog. Ser. 401, 221–232. doi: 10.3354/meps08422
Hardy, O. J., and Vekemans, X. (1999). Isolation by distance in a continuous population: reconciliation between spatial autocorrelation analysis and population genetics models. Heredity 83, 145–154. doi: 10.1046/j.1365-2540.1999.00558.x
Hebert, P. D. N., Ratnasingham, S., and Jeremy, R. (2003). Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc. Biol. Sci. 270, 96–99. doi: 10.1098/rsbl.2003.0025
Hoeinghaus, D. J., Agostinho, A. A., Gomes, L. C., Pelicice, F. M., Okada, E. K., Latini, J. D., et al. (2009). Effects of river impoundment on ecosystem services of large tropical rivers: embodied energy and market value of artisanal fisheries. Conserv. Biol. 23, 1222–1231. doi: 10.1111/j.1523-1739.2009.01248.x
Huergo, G. M., Filgueiras-Souza, R. J., Batista, J. D. S., Formiga-Aquino, K., and Alves-Gomes, J. A. (2011). Molecular genetics as a tool for fisheries management in the Brazilian Amazon: Piraíba (Brachyplatystoma filamentosum and Brachyplatystoma capapretum) (Siluriformes: Pimelodidae) in white-water rivers. Pan Am. J. Aquat. Sci. 6, 280–289.
Iervolino, F., Kawakami, E., Resende, D., Wagner, A., and Hilsdorf, S. (2010). The lack of genetic differentiation of pacu (Piaractus mesopotamicus) populations in the Upper-Paraguay Basin revealed by the mitochondrial DNA D-loop region: implications for fishery management. Fish. Res. 101, 27–31. doi: 10.1016/j.fishres.2009.09.003
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
Marques, A. A. B., Fontana, C. S., Vélez, E., Bencke, G. A., Schneider, M., and Reis, R. E. (2002). Lista das Espécies da Fauna Ameaçadas de Extinção no Rio Grande do Sul. Porto Alegre: Fundação Zoobotânica do Rio Grande do Sul.
Matoso, D. A., da Silva, M., Cortinhas, M. C. S., Cestari, M. M., de Almeida, M. C., Vicari, M. R., et al. (2011). Two genetic stocks of Steindachneridion melanodermatum living in sympatry in nature and genetic variability of wild parents and F1 generation. Genet. Mol. Res. 10, 2606–2612. doi: 10.4238/2011.October.25.8
Monaghan, M. T., Wild, R., Elliot, M., Fujisawa, T., Balke, M., Inward, D. J. G., et al. (2009). Accelerated species inventory on Madagascar using coalescent-based models of species delineation. Syst. Biol. 58, 298–311. doi: 10.1093/sysbio/syp027
Neville, H. M., Isaak, D. J., Dunham, J. B., Thurow, R. F., and Rieman, B. E. (2006). Fine-scale natal homing and localized movement as shaped by sex and spawning habitat in Chinook salmon: insights from spatial autocorrelation analysis of individual genotypes. Mol. Ecol. 15, 4589–4602. doi: 10.1111/j.1365-294X.2006.03082.x
Ochoa, L. E., Pereira, L. H. G., Costa-Silva, G. J., Roxo, F. F., Batista, J. S., Formiga, K., et al. (2015). Genetic structure and historical diversification of catfish Brachyplatystoma platynemum (Siluriformes: Pimelodidae) in the Amazon basin with implications for its conservation. Ecol. Evol. 5, 2005–2020. doi: 10.1002/ece3.1486
Pereira, L. H. G., Hanner, R., Foresti, F., and Oliveira, C. (2013). Can DNA barcoding accurately discriminate megadiverse Neotropical freshwater fish fauna? BMC Genet. 14:20. doi: 10.1186/1471-2156-14-20
Pons, J., Barraclough, T. G., Gomez-Zurita, J., Cardoso, A., Duran, D. P., Hazell, S., et al. (2006). Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst. Biol. 55, 595–609. doi: 10.1080/10635150600852011
Primmer, C. R., Veselov, A. J., Zubchenko, A., Poututkin, A., Bakhmet, I., and Koskinen, M. T. (2006). Isolation by distance within a river system: genetic population structuring of Atlantic salmon, Salmo salar, in tributaries of the Varzuga River in northwest Russia. Mol. Ecol. 15, 653–666. doi: 10.1111/j.1365-294X.2005.02844.x
Rambaut, A., Suchard, M., Xie, D., and Drummond, A. (2014). Tracer v1.6 | BEAST. Available at: http://beast.bio.ed.ac.uk/tracer [accessed March 27, 2017].
Ramella, M. S., Kroth, M. A., Meurer, S., De, A. P., Nuñer, O., Filho, E. Z., et al. (2006). Genetic variability in four fish species (Pimelodus maculatus, Prochilodus lineatus, Salminus brasiliensis and Steindachneridion scripta) from Uruguay River Basin. Braz. Arch. Biol. Technol. 49, 589–598. doi: 10.1590/S1516-89132006000500008
Ribolli, J., Hoeinghaus, D. J., Johnson, J. A., Zaniboni-Filho, E., de Freitas, P. D., and Galetti, P. M. (2017). Isolation-by-time population structure in potamodromous Dourado Salminus brasiliensis in southern Brazil. Conserv. Genet. 18, 67–76. doi: 10.1007/s10592-016-0882-x
Rozas, J., Sánchez-DelBarrio, J. C., Messeguer, X., and Rozas, R. (2003). DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19, 2496–2497. doi: 10.1093/bioinformatics/btg359
Schork, G., Hermes-Silva, S., Beux, L. F., Zaniboni-Filho, E., Nuñer, A. P., and de, O. (2012). Diagnóstico da pesca artesanal na usina hidrelétrica de machadinho, alto Rio Uruguai - Brasil. Bol. Inst. Pesca 38, 97–108.
Schork, G., Hermes-Silva, S., and Zaniboni-Filho, E. (2013). Analysis of fishing activity in the Itá reservoir, upper uruguay river, in the period 2004-2009. Braz. J. Biol. 73, 559–571. doi: 10.1590/S1519-69842013000300014
Sivasundar, A., Bermingham, E., and Ortí, G. (2001). Population structure and biogeography of migratory freshwater fishes (Prochilodus: Characiformes) in major South American rivers. Mol. Ecol. 10, 407–417. doi: 10.1046/j.1365-294X.2001.01194.x
Swarça, A. C., Fenocchio, A. S., Cestari, M. M., and Dias, A. L. (2005). First chromosome data on Steindachneridion scripta (Pisces, Siluriformes, Pimelodidae) from Brazilian rivers: Giemsa, CBG, G-, and RE banding. Genet. Mol. Res. 4, 734–741.
Ward, R. D., Zemlak, T. S., Innes, B. H., Last, P. R., and Hebert, P. D. N. (2005). DNA barcoding Australia’s fish species. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 360, 1847–1857. doi: 10.1098/rstb.2005.1716
Wu, C. H., and Drummond, A. J. (2011). Joint inference of microsatellite mutation models, population history and genealogies using transdimensional Markov Chain Monte Carlo. Genetics 188, 151–164. doi: 10.1534/genetics.110.125260
Zaniboni-Filho, E., and Schulz, U. H. (2003). “Migratory fishes of the Uruguay River,” in Migratory Fishes of South America. Biology, Fisheries and Conservation Status, ed. J. Carolsfeld (Washington, DC: World Bank), 157–194.
Keywords: control region mitochondrial DNA, conservation of natural resources, DNA barcode, endangered species, freshwater fishes
Citation: Paixão RV, Ribolli J and Zaniboni-Filho E (2018) Genetic Variation of the Endangered Neotropical Catfish Steindachneridion scriptum (Siluriformes: Pimelodidae). Front. Genet. 9:48. doi: 10.3389/fgene.2018.00048
Received: 17 October 2017; Accepted: 31 January 2018;
Published: 19 February 2018.
Edited by:Rodrigo A. Torres, Universidade Federal de Pernambuco, Brazil
Reviewed by:Jorge Luis Ramirez, Federal University of São Carlos, Brazil
Jorge Abdala Dergam, Universidade Federal de Viçosa, Brazil
Copyright © 2018 Paixão, Ribolli and Zaniboni-Filho. 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: Rômulo V. Paixão, firstname.lastname@example.org