1. Introduction
The Brazil has a coastal area with more than 8,000 kilometers, with a large volume of fresh water, is one of the largest hydrographic basins in the world, with a great marine biodiversity, due this scenario the Brazilian marine fish farming is in constant expansion and economic growth [
1,
2]. In 2012, the world production of fish was 158 million tons, moving US
$ 600 billion / year - US
$ 136 billion/year in exports. While one hectare of land can produce 0.12 tons/year of meat, the same area of water can produce 100 to 320 tons/year of fish. Making fish farming a more sustainable activity. With this potential, Brazil intends to increase aquaculture production [
2].
Among the various marine species with commercial interest and potential to be exploited, the species of the genus
Centropomus spp stand out.
Centropomus spp is the sole genus of the predominantly marine fish family Centropomidae (Teleostei) and have 12 species, commonly known as snooks (
http://researcharchive.calacademy.org/research/ichthyology/catalog/SpeciesByFamily.asp). Species of this genus can be found in both the eastern Pacific and western Atlantic [
3,
4], which was separated around 1.6 to 3 Mya after the Panama Isthmus closure [
5]. Among them, 6 of them inhabit the Brazilian coasts and we highlighted
Centropomus undecimalis (common snook, sergeant fish or robalo) and
Centropomus parallelus (fat snook, smallscale fat snook, little snook or chucumite) as the of highest interest for economic and sport fishing, which have wide distribution from southern Florida (Gulf of Mexico) to southern Brazil (Rio Grande do Sul).
These species are carnivores feeding mainly on fish, shrimp and crustaceans, however, there is a certain seasonal variation of feeding because they seem to be quite opportunistic predating depending on the availability of the environment in which they are [
6,
7], in view of the nobility of its meat, with excellent organoleptic quality, reaching high prices in the domestic and foreign markets, these attributes make these species to be exploited with sportive, artisanal and professional fishing and lack of data for catch reports generate a deficiency for their conservation [8-12] .
C. undecimalis is one of the largest members of the genus (140 cm and 22 kg) [
13] and has been included in several aquaculture programs in coastal areas of the Atlantic Ocean [4, 14, 15].
C. parallelus may reach up to 75 cm and weigh up to 4 kg with females being larger than males of the same age [
16]. Snook species are normally commercialized with size is from 0.5 to 1.0 kg [
16].
For the advancement of efficiency in the production and profitability of aquaculture, a diversified range of genetic resources is necessary. Research on genomics, genetics and breeding in aquaculture has the ultimate objective of improving production efficiency, sustainability, product quality, food safety, consumer protection and improving profitability, all with the goal of benefiting consumers [
17]. Thus, knowledge of genetic diversity of species is a tool to ensure the sustainability of future populations within their geographic range, despite habitat alterations, as well as captive breeding [18-20].
Population genetics are linked to allelic distribution and change in natural selection, mutation, migration, gene flow, geographic distribution, and genetic drift [
21]. Loss of genetic variability leads to problems in fish adaptation and survival, genetic knowledge of the population structure of common and fat snook leads to gains in fish production and conservation. The genetic diversity between the natural reserves and/or between the different fish farming systems is fundamental for the correct management of the species [
22]. There are several tools for genetic analysis today, mainly the analysis of the genome or regions in the genome and among them, microsatellites stand out. Microsatellite markers, also known as Short Tandem Repeats (STR), are 90–350 bp long non-transcriptional variable number of tandem repeat sequences at certain points of the genome. They are co-dominant and generally high polymorphic and are widely used for genetic diversity analysis and quantitative trait locus [23- 26].
Thus, in this study, we used microsatellite markers to investigate the genetic variability and the genetic characterization of C. undecimalis and C. parallelus populations, identifying if the citations species are geographically connected, using 10 loci of microsatellites, bringing strategies of management, sustainability and conservation for species citations, and providing genetic subsidies needed to be used in breeding programs for snooks.
3. Results and Discussion
The analyzed microsatellites showed high polymorphism of the ten loci used to verify the genetic variability of
C. udecimalis and
C. paralellus populations. Microsatellite loci have been evaluated in many species of fish, studies with these markers are being carried out, since the number of alleles per loco is quite variable. Microsatellites are abundant in the genomes of higher organisms [
34]. Thousands of such microsatellite markers have been identified in humans [34, 35]. For
C. undecimalis, microsatellites were described by Seyoum et al. [
28] and a study was performed by Hernández-Vidal et al. [
36]. No study with microsatellite markers has been performed in
C. parallelus so far. The alleles length used in this study from
C. undecimalis vary from 80 bp for locus Cun21B to 293 bp for locus Cun09. In the study by Seyoum et al. [
28] the variation was 100 bp for locus Cun21B to 243 bp for locus Cun14.
Numbers of alleles vary in the 3 populations of C. udecimalis. In in uSC, the locus Cun08 and Cun11 showed the highest number of all alleles (22), while loco Cun05B showed the lowest number of alleles (7). In uSP, the locus Cun12 showed the highest number of all alleles (12), while loco Cun05B showed the lowest number of alleles (3). In uCap, the locus Cun01 showed the highest number of all alleles (32), while loco Cun05B showed the lowest number of alleles (13). In C. paralellus in pSC, the locus Cun08 showed the highest number of all alleles (24), while loco Cun05B showed the lowest number of alleles (4). In pSP, the locus Cun14 showed the highest number of all alleles (15), while loco Cun16 showed the lowest number of alleles (6).
Although the loci with the highest number of alleles varied, the locus that showed the lowest number of alleles was Cun05B in all populations of C. paralellus and C. udecimalis ranging from 3 to 13.
A high number of different alleles per locus is a good indicator of high polymorphism. However, this is often strongly related to the sample size of each group studied [37, 38], for example, in the case of uCap, which showed the largest number of alleles at all loci. Thus, there are other parameters to be considered in a population genetics analysis. Ho and He are parameters that are quite reliable, and the locus Cun05B also showed the lowest Ho in general. In the uSP population, this same locus showed a Ho of 0.5, while all other loci showed a Ho of 1.
In genetics of populations the presence of null alleles has been commonly reported and may interfere with estimates of genetic variability and genetic distances [
39]. Due to the high He index observed for all loci, it was estimated there is no null allele in the locus analyzed in this study of the populations of
C. undecimalis and
C. paralellus. However, null allele was observed in
C. undecimalis in the locus Cun5B, Cun16 and Cun21B in uSP; in the locus Cun14, Cun21B and Cun22 in uSC; and in Cun10A, Cun11, Cun14 and Cun22 in uCap; while in
C. paralellus in the locus Cun10A, Cun11 and Cun16 in pSC; and in 14, 21B and 22 in pSP.
In the present study, genotypic and genotypic differentiation values, based on allele and genotype frequencies observed, indicated a possible genetic differentiation in the analyzed populations. Population differentiation studies in fish have been carried out to identify and manage different genetic pools. Population differentiation can occur in different forms (from genetically distinct populations to panmitic populations) and in this way, the basic conservation unit would be represented by a group genetically homogeneous individuals, such as the existence of different local populations [
40].
Regarding the populations of C. undecimalis, in fSP population the heterozygosity observed (HO) ranged from 1 to 0.5 with a mean of 0.95 whereas the heterozygosity expected (HE) ranged from 0.8 to 0.91 with an average of 0.82. In fSC population the HO ranged from 1 to 0.39 with an average of 0.89 whereas the HE ranged from 0.9 to 0.94 with an average of 0.89. In fCap population the HO ranged from 1 to 0.68 with an average of 0.93 whereas the HE ranged from 0.95 to 0.85 with an average of 0.92.
Regarding the populations of C. paralellus, in pSC population the HO ranged from 1 to 0.55 with an average of 0.86 while the HE ranged from 0.9 to 0.94 with an average of 0.89. In pSP population the HO ranged from 1 to 0.77 with an average of 0.95 whereas the HE ranged from 0, 9 to 0.85 with a mean of 0.84.
F
IS (coefficient of inbreeding) for those who did not present negative values indicating that there is no excess of heterozygotes in these populations. F
IS, varied among all locus of the populations and the two species (
Table 2 and
Table 3). The distance and the genetic identity among the
C. undecimalis populations was evidenced by the F
ST (fixation index) value of 0.0376 with an IC of (0.0250-0.0532), which according to Wright [
41] are indications of low genetic differentiation. PCA data showed that there is a genetic distinction between populations, especially between pSC and pSP populations of
C. paralellus, and in uSC and uSP of
C. undecimalis, however uCap is more like the uSP population than to uSC (
Figure).
The locus (Cun01, Cun05B, Cun08, Cun09, Cun10A, Cun11, Cun14, Cun16, Cun21B, and Cun22) did not exhibit Hardy-Weinberg equilibrium according to populations and species [
42]. There was no linkage disequilibrium between the locus, suggest that all are genetically independent.
Its capture is associated with migratory movements to freshwater ecosystems and, also, to spawning on the coast. The environmental differences found along the coast play an ecological and evolutionary role, favoring genetic differentiation. So far, there has been no study comparing the species and populations of C. undecimalis and C. parallelus between the coast of the Brazilians federated states of São Paulo and Santa Catarina [43- 45].
Figure 1.
and second (PC2) principal components based on the analysis of the individual C. undecimalis (A) and C. parallelus (B) using 10 microsatellites loci. The ellipses, shown in different colours, represent the analysed individuals in two different analyses.
Figure 1.
and second (PC2) principal components based on the analysis of the individual C. undecimalis (A) and C. parallelus (B) using 10 microsatellites loci. The ellipses, shown in different colours, represent the analysed individuals in two different analyses.
Although the scientific literature lacks sources about these two species, especially regarding their genetics, some considerations can be made. In summary, it was demonstrated by our data that these populations are found with high genetic variability, including the recently domesticated population of C. parallelus. The populations showed distinct allelic frequencies of some alleles, which demonstrates that they are consolidated populations. The Figure demonstrates that there are some points of convergence between populations of C. parallelus, which may indicate some reproductive contact. In addition, the same figure demonstrates that the domesticated population of C. undecimalis is more related to the population of the state of São Paulo than to Santa Catarina.