Intensive Commercialization of Endangered Sharks and Rays (Elasmobranchii) Along the Coastal Amazon as Revealed by DNA Barcode

Elasmobranchs represent a well-defined group, composed of about 1,150 species inhabiting diverse aquatic environments. Currently, several of these species have been classified as threaten due to overexploitation. Therefore, we used DNA barcode to identify traded species of sharks and stingrays in the municipality of Bragança (Amazon coastal region), a major fishery landing site in northern Brazil. We collected a total of 127 samples labeled into 24 commercial nomenclatures over 1-year period. Twenty species were discriminated and 13 of them are recognized under some threatening status. In relation to sharks, Carcharhinidae showed the highest number of species, half of them classified as endangered, followed by hammerhead sharks (Sphyrnidae), with four species also regarded as threatened with extinction. The Rhinopteridae and Dasyatidae rays were the most abundant groups, with trade records of the following threatened species: Rhinoptera brasiliensis, Rhinoptera bonasus, Hypanus berthalutzae, and Fontitrygon geijskesi. It is noteworthy that threatened species of elasmobranchs have been frequently and regularly exploited because of inefficient fishery management policies. Therefore, effective inspection practices need to be incorporated in fisheries, including the use of DNA barcode to enable a reliable method of species authentication and to assure the proper commercialization.


INTRODUCTION
Sharks and rays (Subclass Euselachii) are elasmobranchs whose origin dates back to the Devonian (∼400 million years ago) (Compagno, 1990). Currently, this group is represented by nearly 1,150 valid species, found mainly in marine habitats even though some of them are typical of freshwater environments (Nelson et al., 2016).
In spite of their wide distribution, sharks and rays have been severely threatened with extinction, mainly because of overfishing (Stevens et al., 2000;Pacoureau et al., 2021). The fishing pressure over this fish group has determined remarkable reductions of population size particularly because elasmobranchs present peculiar biological traits, such as longevity, late sex maturation and low fecundity rates (Stevens et al., 2000;Oliver et al., 2015), making this group even more vulnerable to overfishing. Global data provided by Food and Agriculture Organization of United Nations pointed out that nearly 30% of elasmobranchs are vulnerable, endangered or critically threatened with extinction in nature (IUCN, 2021). The population drops of sharks and rays by overfishing are directly associated with the commercialization of fins in Asian markets (Pough et al., 2008), bycatch and artisanal fishing (Robbins et al., 2006;Graham et al., 2010;Ferrette et al., 2019).
Both artisanal and industrial fishing practices to supply the demand of elasmobranch-derived products in regional markets take place regularly in several points along the northern coast of Brazil, where nearly 69 species of shark and rays are recognized, many of them classified as threatened Marceniuk et al., 2019;Rodrigues-Filho et al., 2020). As a matter of fact, many species of threatened sharks and rays are commercialized in coastal Amazon region . In the municipality of Bragança, northeastern Pará, these species are locally sold under the denomination of "cação" (sharks) and "arraias" (rays) in street markets, representing an important source of income (Freire et al., 2011;Martins et al., 2021).
During commercialization, the morphological identification of sharks and rays is usually jeopardized because body parts, such as fins, are removed right after catches in open sea (Rodrigues-Filho et al., 2009;Feitosa et al., 2018). In addition, the Brazilian laws lack a standardized nomenclature for the commercialization of these species, thus favoring illegal trade since the commercial common names often refer to a group of species (Ardura et al., 2010;Feitosa et al., 2018;Martins et al., 2021).
Inasmuch as useful meristic and morphometric traits are unavailable in the commercialized sharks and rays, alternative approaches should be implemented to assure the unambiguous identification of species. Accordingly, DNA barcode based on the sequencing of a portion of the mitochondrial gene Cytochrome C oxidase-subunit I (COI) has been widely used to discriminate and identify species of several groups, including fishes and their products that could not be properly validated by morphological traits (Ward et al., 2005;Steinke et al., 2017;Feitosa et al., 2018;Gomes et al., 2019;Carvalho et al., 2020).
Bragança represents a major fishing market in northeastern Pará in Amazon coastal region, northern Brazil (Freire et al., 2011;Santana et al., 2020;Martins et al., 2021). Fish landing occur daily in this municipality (Isaac et al., 2008) and most of products are sold in the street market of Bragança that encompasses indoor ("Mercado") and outdoor ("Feirinha") sites where a high diversity and number of species are commercialized (Freire et al., 2011;Martins et al., 2021). Over more than a decade, studies in this region have shown that Bragança lands and commercializes a large diversity of elasmobranchs with variation in the abundance and richness of sharks and rays species sold in the street market (Carmona et al., 2008;Rodrigues-Filho et al., 2009Guimarães-Costa et al., 2020;Feitosa et al., 2018;Marceniuk et al., 2019;Souza-Araujo et al., 2021). The elasmobranchs labeled as "cação" and "arraia" are among the most commercialized fish groups in this region but the lack of standard commercial common names and the mischaracterization of specimens constrain their morphological identification and justify the utilization of other methods based on molecular validation of species, like DNA barcode Rodrigues-Filho et al., 2020;Martins et al., 2021).
Therefore, in the present study, the dynamics of commercialization of sharks and rays was investigated in northeastern Pará, in order to identify the commercialized species based on DNA barcode and potential records of threatened species in local markets.

Sampling
The collections were carried out by incursions in 15-day intervals to the street market of Bragança, state of Pará, Brazil over a year (February 2017 to January 2018) to cover the species that are traded throughout time space. The elasmobranchs traded in the Bragança street market are caught, mainly, as bycatch, in fisheries of acoupa weakfish Cynoscion acoupa (Lacèpede, 1801), pink-shrimp Penaeus sp., piramutaba catfish Brachyplatystoma vaillantii (Valenciennes, 1840) and several fish (Mourão et al., 2018;Marceniuk et al., 2019). These fisheries are conducted on a small and large scale by different fishing systems operating in different areas of the northern Brazilian coast (Amapá, Pará, and Maranhão), using a wide variety of fishing gear (Braga et al., 2006;Espírito Santo and Isaac, 2012).
The collection schedule was defined based on a previous report in the same location (Martins et al., 2021). A total of 127 samples of elasmobranchs were collected, 35 labeled as "ray" and 92 under the label "shark" (Supplementary Material 1). The samples were acquired in indoor and outdoor market ("feirinha") with available supply of rays and sharks on collection expeditions, avoiding duplication of samples from a single individual, especially in the case of large sharks that are usually sold by more than one fishmonger.
Specimens were previously identified using the commercial nomenclatures used by traders at the time of collection (Supplementary Material 1) to verify if there was a relationship between genetic identification and commercial nomenclatures for sharks and rays (Supplementary Material 1).
A large portion of the biological tissue samples collected were from uncharacterized specimens (sharks: headless or finless samples; rays: headless samples, filets, or only the pectoral fins) (Supplementary Figure 1). Only a small portion of samples encompassed whole individuals.
The whole specimens of rays (Hypanus guttatus, n = 2; Gymnura micrura, n = 3) and shark (Rhizoprionodon porosus, n = 1) were identified by specialized literature (Compagno, 1984;Mc Eachran and Carvalho, 2002) and used as vouchers to compose a reference sample database to assist in the identification process from public sequence banks (Genbank and BOLD Systems). These specimens were deposited in the ichthyological collection of the Laboratório de Genética Aplicada (LAGA), Instituto de Estudos Costeiros (IECOS), UFPA-Bragança.

Laboratory Procedures
Total DNA was extracted using the Wizard Genomic DNA Purification kit (Promega Corporation, Madison, WI, United States), according to the manufacturer's instructions. To amplify the barcode portion of the COI gene we used the primers FishF1 (5 TCAACCAACCACAAAGACATTGGCAC 3 ) and FishR1 (5 TAGACTTCTGGGTGGCCAAAGAATCA 3 ) (Ward et al., 2005).
Positive PCRs were purified with PEG (Polyethylene Glycol) according to Paithankar and Prasad (1991). The sequencing was performed by the dideoxy terminal method (Sanger et al., 1977) with the Big Dye kit (ABI Prism Dye Terminator Cycle Sequencing Reading Reaction-Thermo Fisher Scientific) according to the manufacturer's instructions, using FishF1 primer. After precipitation, the purified PCR products were submitted to capillary electrophoresis in ABI 3500 Genetic Analyzer automatic sequencer (Thermo Fisher Scientific).

Database and Analysis
The sequences were edited in the software BioEdit v. 7.1.3.0 (Hall, 1999), by visual inspection of chromatogram files in ABI format, using a threshold of ≥30 for PHRED quality score of base calling for trimming. DNA sequences were aligned using CLUSTAL W (Thompson et al., 1994;Larkin et al., 2007), available in the software BioEdit v. 7.1.3.0 (Hall, 1999). Presence of stop codons was evaluated in the BOLD Systems platform (Ratnasingham and Hebert, 2007) To help identify the taxa, the haplotypes were defined in the software DNAsp (Librado and Rozas, 2009). This identification was performed based on the similarity of our sequences with those available in GenBank (NCBI-National Center for Biotechnology Information) 2 and BOLD (Barcode of Life Data Systems-reference libraries), 3 using the highest similarity values (99-100%) (Supplementary Material 1) to discriminate species (Hebert et al., 2003). Nineteen reference sequences from GenBank (Supplementary Material 2) were used. The maximum and average barcode gap distances were evaluated in the BOLD Systems platform using the Barcode Gap Analysis tool (Ratnasingham and Hebert, 2019).
Bayesian inference (BI) was performed in the software BEAST v. 1.8.4 (Drummond et al., 2012) assuming the following parameters: strict clock, Yule speciation process, HKY + I + G evolutionary model as selected by jModelTest2, 50 million generations and 10% of burn-in. The log files were checked in Tracer v. 1.7  to assess the burn-in length and the convergence based on ESS (Effective Sample Size) parameters above 200. The generated trees were summarized in TreeAnnotator v. 1.8.4 to obtain the consensus tree, which was visualized in FigTree v. 1.4.4 (Rambaut, 2018). The p distance was calculated in MEGA X (Kumar et al., 2018) to compare intra and interspecific distances and to assess the genetic divergence of closely related species (Kumar et al., 2018).
After the genetic identification of traded species, the conservation status was verified in the IUCN Red List of Threatened Species (IUCN, 2021) and Brazil Red Book of Threatened Species of Fauna (ICMBio, 2018a). The graph of threatened categories and the graph relating the traded categories and the identified species were built using the online tool RAWGraphs 4 .

Molecular Identification and Genetic Similarity
A dataset of 127 samples (92 of sharks and 35 of rays) comprising COI fragments of 599 base pairs were obtained, referring to 24 commercial nomenclatures and totaling 31 haplotypes (Supplementary Material 1). All generated sequences were stored in GenBank under accession numbers MZ303812-MZ303842.
A total of 20 species was discriminated by DNA barcode. Nineteen of them were identified based on their genetic similarity with sequences from public datasets, most of them exhibiting values equal to 100% (Table 1).

Identification Based on Phylogenetic Inference
The Bayesian tree topology showed 20 clades, being 9 of them related to rays (order Myliobatiformes) and 11 shark taxa from orders Orectolobiformes and Carcharhiniformes. The representatives of Myliobatiformes encompassed two species of Rhinopteridae, three species of Dasyatidae, one species of Potamotrygonidae, one species of Gymnuridae and one species of Aetobatidae, as well as specimen identified as Myliobatiformes 1. The order Orectolobiformes was represented only by species of family Ginglymostomatidae (one species) while the Carcharhiniformes were represented by six species of Carcharhinidae and four species of Sphyrnidae.
In the case of rays, the most representative families were: Rhinopteridae (N = 11) with two species-Rhinoptera brasiliensis   (Figure 1B).
The specimen designated here as Myliobatiformes 1 was not suitable for molecular identification at the species level because it showed genetic similarity with sequences from public banks of up to 88% with species of the families Rhinopteridae and Mobulidae. The IB tree showed the grouping of this taxon with the species Mobula birostris (Walbaum, 1792) and Mobula hypostoma (Bancroft, 1831) and R. brasiliensis, R. bonasus (Supplementary Material 2). However, the distance p recorded for Myliobatiformes 1 to M. birostris and M. hypostoma was 13.11 and 12.57%, respectively, and from Myliobatiformes 1 to R. brasiliensis and R. bonasus was 13, 98%, and 12.62%, respectively.

Commercial Labels and the Supply of Threatened Species
Fifteen commercial nomenclatures were recorded for sharks, representing 11 species, while nine commercial labels were recorded for rays (Figure 2) with discrimination of nine taxa (eight identified at species level and one of them at the order level). Thirteen out of the 20 listed species are considered endangered and five of them are regarded as nearly threatened (ICMBio, 2018a;IUCN, 2021). Amongst the identified species, five are classified as data deficient (DD), one as not evaluated (NE), one as endangered (EN), one as vulnerable (VU) according to Brazilian lists (ICMBio, 2018a), while according to international red lists, the samples comprised one least concern (LC), five near threatened (NT), and five critically endangered (CR) species (IUCN, 2021; Figure 3).
The collected of sharks and rays over a 1-year period revealed regular supplies for some species, such as S. mokarran and S. lewini, sampled nearly every month. Other species were collected occasionally, including R. brasiliensis, C. plumbeus, C. porosus, S. tudes, S. tiburo, and C. falciformis.

Diversity of Elasmobranchs in the Amazon Coast
A total of 20 taxa of elasmobranchs were discriminated via DNA barcode among the samples collected over a year in the street market of Bragança. Nineteen of them could be identified at species-level revealing the high efficiency of this method to assess the diversity of commercialized sharks and rays. Previous molecular reports based on molecular markers carried out in northern Brazilian coast identified similar numbers of elasmobranch species in fish landings and local markets (N = 26), comprising 17 taxa of sharks and nine species of rays Rodrigues-Filho et al., 2020).
Ten out of the 17 shark species reported in northern Brazilian coast  were found in the street Market of Bragança. In addition, we recorded for the first time the presence of Carcharhinus plumbeus in local markets from this region (Rodrigues-Filho et al., 2009;Feitosa et al., 2018).
A total of 11 species of sharks are commercialized in the street market of Bragança, as similarly reported by Rodrigues-Filho et al. (2009) using 12S and 16S mitochondrial DNA markers. Nonetheless, the species composition between both reports differed inasmuch as we identified some species that have never been reported before such as C. plumbeus, G. cirratum, S.  As for rays, nine species have been commercialized in Bragança, including two unrecorded species in previous studies (Carmona et al., 2008): H. berthalutzae and R. brasilenses.
However, in both reports some taxa of rays could not be precisely identified. A potential obstacle to a proper identification of the species was the lack of whole specimens to be used as a voucher for morphological identification as well as the absence of sequences in public datasets for comparative analyses (Rodrigues-Filho et al., 2020). In recent years, studies based on molecular identification of stingray species have allowed both the identification of valid species and the registration of new occurrences and new species (Guimarães-Costa et al., 2020;Pentean et al., 2020). In a recent study, Rodrigues-Filho et al. (2020) identified nine ray species on the Amazon coast using DNA barcodes. However, the unidentified stingray taxon (Myliobatiformes 1) discovered here has no genetic match with any of the species presented by Rodrigues-Filho et al. (2020), nor with any other species previously reported in this region with COI sequences in public banks.
Even though 69 species of elasmobranchs are recognized for the Amazon coast (Marceniuk et al., 2019), the diversity of sharks and rays from this region remains overlooked. As a matter of fact, as long as the number of studies in this region has increased, new discoveries have been reported. Recently, morphological analyses reported the following species for the first time in Amazon coast: Narcine bancroftii, Pseudobatos horkelii, Urobatis jamaicensis, and Narcine sp. (Marceniuk et al., 2019). Similarly, molecular data identified Narcine brasiliensis and Mustelus canis (Guimarães-Costa et al., 2020;Rodrigues-Filho et al., 2020) while integrative analyses based on morphological, genetic and ecological studies recorded H. berthalutzae (Pentean et al., 2020) for this region. These data highlight the importance of focusing on the diversity of elasmobranchs in the north coast of Brazil, as recent studies have identified distinct results on the species composition in this group of fishes.

Commercialization of Threatened Species
The present results showed that 13 out of the 20 taxa of elasmobranchs are threatened with extinction and five of them are nearly threatened according to the Red List of Threatened Brazilian Fauna (ICMBio, 2018a) and the International Red List of Threatened Species (IUCN, 2021). When both national and international lists are compared, we noticed species classified under different threatening statuses. For instance, R. bonasus is regarded as DD by ICMBio (2018a) but as VU according to IUCN (2021), while C. falciformis is classified as NT and VU by ICMBio (2018a) and IUCN (2021), respectively.
Elasmobranchs tend to be sold uncharacterized and often using commercial names, usually referring to a general fish category. The lack of standardization in nomenclature can also make it easier for species with different threat status to be sold together. A case, frequently observed in the present study, was the utilization of the local term "cação panã" to refer to the following species: R. porosus, C. porosus, S. tiburo, S. tudes, S. mokarran, and S. lewini. Similarly, the name long-nose stingray or "arraia bicuda" was indiscriminately used to refer to H. berthalutzae, H. guttatus and F. geijskesi. These data highlight that commercial labels are not a reliable source of information for the species identification, as previously reported by other authors (Clarke et al., 2006;Rodrigues-Filho et al., 2009;Domingues et al., 2013).

Commercialization of Shark Species
Six species of the family Carcharhinidae (R. porosus, C. porosus, C. leucas C. falciformis, C. plumbeus, and C. acronotus) were commercialized, presenting the highest species richness but with lower frequency in relation to those representatives of the family Sphyrnidae.
Three out of the six Carcharhinidae species herein identified are classified as threatened, like C. porosus whose status is defined as "critically endangered." This species was rarely reported in the present study, thus differing from a previous report (Rodrigues-Filho et al., 2009) that placed C. porosus as the most commercialized shark species in fish markets of Bragança. The population decline of C. porosus, formerly abundant in catches along the northern Brazilian coast, has increased from 1990 to 2004 (Rodrigues-Filho et al., 2009;ICMBio, 2018b), associated with cases of local extinction (ICMBio, 2018b;Dulvy et al., 2021).
The species R. porosus was observed being commercialized, similarly other reports have detected the increased commercialization of this species Martins et al., 2018). Even though R. porosus is not considered a threatened species, population approaches are encouraged since most of samples collected in the present study referred to newborns removed from pregnant females.
On the other hand, C. falciformis, C. plumbeus, and C. leucas, usually exploited in fin markets (Clarke et al., 2006), encompassed a few samples in this study, differing from other reports in southeastern Brazilian coast, where C. falciformis was the most abundant shark species (51%) in fish landings (Domingues et al., 2013). These data show that the supply of these species is sporadic, possibly because they derive from bycatches or these species may be being marketed in other locations (Lessa et al., 1999).
In relation to the family Ginglymostomatidae, the species G. cirratum is regarded as vulnerable according to Brazilian official reports (ICMBio, 2018a). Other reports have indicated that stocks of G. cirratum along the Brazilian coast have decreased over the last years (SBEEL, 2005).
Based on these results, we observed a change in the pattern of shark commercialization with the highest abundance for Sphyrnidae while previous reports indicated representatives of Carcharhinidae as the most commercialized species (Rodrigues-Filho et al., 2009;Feitosa et al., 2018). Within the family Sphyrnidae, S. mokarran was the most frequently exploited species as also observed in another study carried out along the northern Brazilian coast that placed this shark as this shark is the fourth most found and sampled species . These data are worrying since all species of the genus Sphyrna identified in the present study are threatened with extinction.
For instance, during the development of this study the status of S. mokarran has changed from endangered to critically endangered (IUCN, 2021). These results are worrisome, since overexploited species tend to have their populations reduced and may lead to extinction (Dulvy et al., 2014).

Commercialization of Ray Species
We identified two species belonging to the family Rhinopteridae, being Rhinoptera brasiliensis considered threatened for being sensitive species to fishing pressure, as the females of this species have only one offspring per gestation (Vooren et al., 2005;IUCN, 2021). The occurrence of the species on the market was not regular, suggesting an occasional trade. This scenario corroborates with previous studies that reported a low number of samples found in Bragança (Rodrigues-Filho et al., 2020). These data are in agreement with previous studies showing that R. brasiliensis, which was previously considered restricted to the southern coast of America (Lessa et al., 1999), actually presents a wider range of occurrence, reaching the northern coast of Brazil (Marceniuk et al., 2019;Rodrigues-Filho et al., 2020).
Another representative of the family Rhinopteridae, Rhinoptera bonasus, was also represented by a few samples in the present work. This species was also reported in another study along the northern coast (Rodrigues-Filho et al., 2020).
The species G. micrura in the present study showed a genetic divergence of 10.10% in relation to the sequences of the G. micrura populations from the Mexican coast (Sarmiento-Camacho and Valdez-Moreno, 2018), suggesting a restriction of the gene flow between populations from the north and south of the Atlantic. On the other hand, the genetic similarity was equal to 100% when the present sequences of this species were compared to other samples from northern Brazil (Rodrigues-Filho et al., 2020). This result might indicate that this taxon encompasses two distinct groups (or even species). In fact, Yokota and De Carvalho (2017) suggested a taxonomic revision of G. micrura, recommending the denomination Gymnura lessae for the populations distributed along United States coast and Gulf of Mexico, while G. micrura would be restricted to central and southern regions of Western Atlantic (Brazil, northern Venezuela, Trinidad, and Tobago).
Dasyatidae family was represented by the commercialization of H. berthalutzae, H. guttatus, and F. geijskesi, species that have been previously considered as data deficient but after a recent reevaluation (IUCN, 2021) are presently regarded as vulnerable, nearly endangered and critically endangered, respectively. According to ICMBio (2018a), H. guttatus is a species of least concern by presenting no significant threats, while H. berthalutzae remains not evaluated (NE) because of its recent description (Pentean et al., 2020) and F. geijskesi was classified as DD (ICMBio, 2018a).
In the present study, we recorded the commercialization of H. berthalutzae (formerly recognized as H. americanus) for the first time. The geographic distribution of this species ranges from the mouth of Amazonas River to the coast of São Paulo, southeastern Brazil, probably representing an endemic species to the Brazilian coast (Pentean et al., 2020). We also identified the trade of the endangered species (IUCN, 2021) S. schmardae belonging to the Potamotrygonidae family.
In the sample labeled as ray, we also found a specie that could not be attributed to any taxon Myliobatiformes 1 (Arr63F136). After evaluating and comparing with all those stingray species described for the North coast with COI sequences available in public banks, Myliobatiformes 1 cannot be precisely identified (Supplementary Figure 2).

Impacts Caused by the Exploitation of Elasmobranchs
When we compared our results with previous studies, the reported species diversity changed. For example, Sphyrna zygaena, Rhizoprionodon terraenovae, and Negaprion brevirostris were found to be traded in the open market in Bragança (Rodrigues-Filho et al., 2009), while the landing of Pristis pristis in Bragança was reported by Palmeira et al. (2013). These referenced species were not recorded in this study throughout the sampling period. The absence of these species can be explained by their commercialization in other places or even that they are showing a possible decline in their natural populations, since these taxa are classified as critically threatened by overfishing (IUCN, 2021).
The evaluation of exploitation of sharks and rays, whether at a regional level or on a global scale, is difficult, once catch volume estimates and landing statistical data in fisheries are often incomplete and superficial (Dent and Clarke, 2015). Fishing practices are particularly harmful to these fishes since most elasmobranchs are k-strategist animals that are more susceptible to exploitation and have slow recoveries after stock depletions (Stevens et al., 2000;Oliver et al., 2015;Pacoureau et al., 2021).
In order to avoid overexploitation of these fisheries resources, deep changes are necessary, such as strict monitoring of what species have been caught and landed, implementation of onboard and landing inspectors, as well as the use of more selective fishing tools to target species thus reducing bycatches (Barreto et al., 2017). Popular names also need to be re-evaluated in future fisheries statistics as they vary greatly among regions and fail in providing reliable information about which species are being exploited (Ardura et al., 2010;Feitosa et al., 2018).
Moreover, surveillance strategies need to be expanded since a large number of elasmobranchs are widely distributed, including nearly 150 species of shark that migrate across national borders, while a quarter of threatened species can migrate across at least 18 countries (Dulvy et al., 2014). Therefore, global public policies are required to the proper conservation of these organisms.
In this context, DNA-based methodologies, such as the use of DNA barcode, as well as the enrichment of public banks with reference DNA sequences for several taxa are particularly useful to obtain accurate data. This information can be further used to support reliable statistical data for the development of efficient management plans in fisheries (Dulvy et al., 2014;Oliver et al., 2015;Barreto et al., 2016).

CONCLUSION
Twenty-four commercial nomenclatures for elasmobranchs were recorded in the street market of Bragança, referring to 11 species of sharks and 9 taxa of rays, as discriminated by DNA barcode. Sphyrnidae and Carcharhinidae were the most sampled families of sharks, while Rhinopteridae and Dasyatidae were the most represented groups of rays. We also identified the regular trade of threatened species, particularly in the family Sphyrnidae (e.g., S. mokarran and S. lewini). This scenario is largely related to the lack of standardized commercial nomenclatures and mischaracterization of traded specimens thus revealing the importance of using molecular tools, especially DNA barcode in species validation.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.

ETHICS STATEMENT
Ethical review and approval was not required for the animal study because the animals were bought at the fish market. They were already dead.

AUTHOR CONTRIBUTIONS
TM and GE-G participated in the study design and structure. TM, PS, and GE-G obtained and organized the database. IS and GE-G generated the computational data. TM, AG-C, and RS performed the statistical analysis. TM, PS, ÍL, AG-C, and GE-G wrote the draft manuscript. TM, PS, ÍL, and GE-G reviewed the manuscript. TM and ÍL performed the final formatting of the manuscript. All authors contributed to the article and approved the submitted version.

FUNDING
This work was supported by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Financial Code 001, by granting a master's scholarship to TM, to CNPq for financing this study (Universal project #439113/2018-0).