InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Biochemistry, Genetics and Molecular Biology » "Microsatellite Markers", book edited by Ibrokhim Y. Abdurakhmonov, ISBN 978-953-51-2798-7, Print ISBN 978-953-51-2797-0, Published: November 30, 2016 under CC BY 3.0 license. © The Author(s).

Chapter 8

Microsatellites for the Amazonian Fish Hypophthalmus marginatus

By Emil J. Hernández‐Ruz, Evonnildo C. Gonçalves, Artur Silva, Rodolfo A. Salm, Isadora F. de França and Maria P.C. Schneider
DOI: 10.5772/65655

Article top

Microsatellites for the Amazonian Fish Hypophthalmus marginatus

Emil J. Hernández‐Ruz1, Evonnildo C. Gonçalves2, Artur Silva3, Rodolfo A. Salm1, Isadora F. de França1 and Maria P.C. Schneider3
Show details


We isolated 41 and characterized 17 microsatellite loci for evaluating the genetic structure of the Amazonian fish Hypophthalmus marginatus, from the Tocantins and Araguaia River in the Eastern Amazonia. Of the 17 selected microsatellite sequences, 15 were dinucleotide repeats, 9 of which were perfect (5–31 repetitions) and 6 were composite motifs. Among these 17 microsatellites, only two were polymorphic. The average number of alleles (Na) observed in the five examined populations ranged from 3.5 to 4.5, while the average observed heterozygosity (Ho) ranged from 0.3 to 0.6. The allelic frequency was less homogeneous at the locus Hm 5 than that for the Hm 13. Genetic diversity was measured in three upstream and two downstream populations under the influence of the Tucuruí Hydroelectric Dam. Our findings provide evidence for low levels of genetic diversity in H. marginatus of the Tocantis basin possibility related to the Dam construction. The Fst and Rst analysis fits well with migratory characteristics of H. marginatus, suggesting the existence of a gene flow mainly in the upstream or downstream directions. To test the hypothesis that the Dam was responsible for the detected reduction on this species genetic diversity, a large number of genetic markers are recommended, covering geographic distribution range of the fish species.

Keywords: hydroelectric dam influence, migratory fishes, population genetic structure

1. Introduction

Migratory freshwater fishes are vulnerable to a variety of anthropogenic impacts, including harvesting, pollution, and other types of habitat disturbance. The impact of dams is not only limited to the transformation of habitats from lentic to lotic but they also isolate the population from their places of spawning to feeding grounds [13]. Therefore, a survival and reproduction of migratory fish can be directly affected by changing thermal and hydrodynamic conditions in their habitat [4]. In the Brazilian Amazonia, over 10 million hectares (ha) of forests are expected to become permanently flooded after the construction of new‐planned dams [5]. Long‐term monitoring of fish populations is available only for Tucuruí Dam in the Tocantins River [610], where alterations following impoundment reduced the fish diversity on the reservoir resulting in the increase of predators such as Cichla spp. Schneider 1801, and Serrasalmus spp. Lacepède 1803 [11]. Furthermore, a drastic reduction in fish production has been observed downstream of the dam, probably due to the low oxygen content of water that runs through the turbines and the blocking of fish migration [6, 12]. Harvesting of freshwater shrimp downstream the dam has dropped from 179 tons in 1981, before the dam construction, to only 62 tons in 1988 three years after dam construction, while fish landings declined from 4726 to 831 tons (the dam was built between 1984 and 1985) [13]. Catches in the reservoir increased to pre‐flooding levels by the early 1990s [8], although nowadays migratory species such as Hypophthalmus marginatus are still not such abundant as before. Therefore, due to the great economic importance to local fisheries, the current genetic structure and species/biodiversity conservation of H. marginatus stocks in the low to medium Tocantins River is a matter of concern and was investigated here. Besides, helping to understand the impact of the construction of the Tucuruí Dam on the genetic variability, our results should contribute to the eventual development of population management strategies for the studied species and will hopefully arise concerns about the building of future dams in the Amazon.

2. Materials and methods

The Tocantins is largely a plateau river, flowing for most of its length within an enclosed valley, draining an area of 343,000 km2. Over the past three decades, its basin has suffered from huge anthropogenic pressures, including widespread deforestation, mining, and the construction of the Tucuruí between 1984 and 1985. This dam is one of the world's largest hydroelectric dams, which has a reservoir of 2840 km2, most of which was originally covered with primary terra firme forest [9, 10].

2.1. Samples

Eighty‐two samples (14–19 per site) obtained from muscle or liver tissue of H. marginatus were stored in absolute ethanol and frozen at -20°C for future genomic DNA extraction, which was performed using Sambrook standard protocol [14]. The specimens were preserved in 4% formaldehyde and deposited in the ichthyological collection of the Museu Paraense Emílio Goeldí (MPEG 13375, MPEG 17486, MPEG 17499 and MPEG 17578). Microsatellite loci were characterized in H. marginatus individuals from four different points of the Tocantins River: Itupiranga (05°06ʹ51.5ʺS, 49°21ʹ34.9ʺW), Tucuruí (04°18ʹ43.06ʺS, 49°19ʹ58.1ʺW), Cametá (02°03ʹ27.5ʺS, 49°20ʹ31.9ʺW), and Abaetetuba (01°40ʹ42.6ʺS, 49°00ʹ16.6ʺW). Additionally, samples from one point of Araguaia River (Conceição do Araguaia, 07°58ʹ10.8ʺ S, 49°11ʹ0.6ʺ W) were included in the analysis (Figure 1).


Figure 1.

Distribution of the five stocks of Hypophthalmus marginatus sampled on the Tocantins and Araguaia Rivers in eastern Amazonia.

2.2. Microsatellite development

Molecular tools such as microsatellites markers can give us a picture of the distribution of the genetic variability of the natural populations of Amazonian fish [15, 16]. To assess the genetic parameters of the Hypophthalmus in the Amazon basin, we developed a partial genomic library of these fishes enriched for microsatellites following the method of selective hybridization with biotinylated probe types (CT)8, conjugated to streptavidin‐coated magnetic beads [17]. After hybridization, the microsatellite‐enriched sequences were amplified by polymerase chain reaction (PCR), ligated into the pGEM‐T Easy Vector (Promega Corp., Madison, USA) and transformed into Escherichia coli TOP 10 electroporated‐competent bacteria. The transformed bacteria were plated on solid medium containing LB‐ampicillin (100 mg/ml) + X‐gal (2%), and after growth, the white colonies containing inserts were transferred and grown in 96‐well plates in a liquid Tartoff‐Hobbs Broth/ampicillin medium. A total of 96 positive clones were sequenced in both directions using the BigDye® Terminator v3.1 Cycle Sequencing Kit, according to the manufacturer's instructions (Applied Biosystems, Carlsbad, USA). The sequences were edited and aligned in the software BioEdit [18].

A total of 96 clones were sequenced, and 17 of these were selected for primer design using the software Primer3 [19]. In all primer pairs, the forward primers had an M13(‐21) tail added to its 5ʹ end [20]. An optimal annealing temperature was inferred using a gradient PCR with temperatures set between 52.3 and 70.5°C (Table 1).

LocusPrimer sequences (5ʹ–3ʹ)Repeat motif (5ʹ–3ʹ)T (°C)Size range

Table 1.

Characteristics of microsatellite loci isolated from Hypophthalmus marginatus.

[i] - *Polymorphic.

Genotyping reactions were carried out in the Biocycler Thermal Cycler MJ96+/MJ96G (Applied Biosystems), in a final volume of 10 μL. Each reaction contained 3.8 μL of MilliQ water, 1.2 μL of 50 mM MgCl2, 1.0 μL of 10 mM dNTPs, 1.0 μL of PCR buffer (100 mM Tris‐HCl, pH 8.5, 500 mM KCl), 0.4 μL of 2 μM tailed forward primer, 0.4 μL of 2 μM fluorescently labeled primer, 0.8 μL of 2 μM reverse primer, 0.3 μL of 2.5 U Taq DNA polymerase, and 1 μL of DNA (50–100 ng/μL). Reactions were submitted to the following cycling profile: hot start at 94°C for 60 s followed by 25 cycles of denaturing at 94°C for 30 s, annealing for 30 s at locus specific temperatures (Table 1), and extension at 68°C for 30 s; labeling step consisted of 20 cycles of denaturing at 94°C for 20 s, annealing at 52°C for 30 s, and extension at 72°C for 60 s; final extension was performed at 72°C for 30 min. The fragments present in 1 μL of PCR product were separated in 8% polyacrylamide denaturing gel in an ALFexpressTM II (Amersham Biosciences, Freiburg, Germany) automatic DNA sequencer. Amplicon size was estimated by the Allele Locator 1.03 software (Amersham Biosciences, India), based on the internal and external standards provided by the manufacturer.

Using the program Micro‐Checker v2.2.3 [21], we verified that null alleles were not present in the data set. Genetic parameters were obtained with the program Arlequin 3.5 [22].

2.3. Data analysis

To measure the genetic variability within each population, the number of alleles per locus (A), effective number of alleles per locus (Ae), the observed heterozygosity (Ho), and expected (He) under Hardy‐Weinberg equilibrium (HWE) for each locus and their averages were calculated using Popgene 32 software [23].

The Genepop 3.1b program [24] was used to test whether the data obtained are significant deviations from Hardy‐Weinberg equilibrium and the occurrence of connection imbalance between the loci analyzed. This program uses a method of Markov chain to get an unbiased estimate of Fisher's exact test to detect a significant deficiency or excess of heterozygotes [25].

The differentiation between the populations was evaluated by comparing peer‐to‐peer two ways: 1. Estimates FST [26], using the Arlequin 3.5 program [22] and 2. Estimates RST [27], by program RSTCalc [28].

RST is an analog of Wright FST, which is based on the stepwise mutation model (SMM) for microsatellite loci and has the particularity of not displaying associated bias to differences in the size of the alleles in the samples of populations and / or differences the variance between loci. Some authors argue that the estimates of FST show lower when compared to RST estimates in same analysis [29].

The existence of linkage disequilibrium between the loci was evaluated by Genepop 3.4 program [24].

3. Results

Among a total of 98 sequenced recombinant clones, 61 clones presented microsatellites. After sequence analysis, it was found that 41 clones showed more than four repeats of microsatellites and we have designed and purchased primers for 17 of SSRs.

Among the 17 selected microsatellite sequences were a hexanucleotide, a tetranucleotide, and 15 dinucleotide repeats. Of 15 dinucleotide repeats, nine were perfect (5–31 repetitions) and six were composite dinucleotide repeat type (Table 1). Of the 17 loci, only two (Hm 5 and Hm 13) were polymorphic among studied samples.

The average number of alleles (A) observed in five populations examined ranged from 3.5 to 4.5, while the average observed heterozygosity (Ho) ranged from 0.3 to 0.6 (Table 2). The allelic frequency was less homogeneous to the Hm 5 locus than for the Hm 13; its most frequent allele was the same for all populations (Table 3).

Hm 5Hm 13
Size (bps)Relative frequenciesSize (bps)Relative frequencies

Table 2.

Average allelic frequencies for two microsatellite loci in five populations of Hypophthalmus marginatus.

Population (N)LocusAAeHoHeP – EHW
Abaetetuba (14)Hm
Mean (SD)4.0 (0.0)4.0 (0.0)0.3 (0.1)0.5 (0.2)0.2 (0.1)
Cametá (19)Hm
Mean (SD)4.5 (2.1)4.2 (1.8)0.5 (0.1)0.6 (0.1)0.2 (0.0)
Itupiranga (15)Hm
Mean (SD)3.5 (2.1)4.0 (1.5)0.6 (0.2)0.6 (0.2)0.5 (0.6)
Tucuruí (19)Hm
Mean (SD)4.0 (1.4)3.9 (1.3)0.5 (0.1)0.6 (0.1)0.3 (0.1)
Conceição (15)Hm
Mean (SD)4.5 (0.7)4.4 (0.6)0.5 (0.0)0.6 (0.1)0.6 (0.3)

Table 3.

Variability intrapopulation in five populations of Hypophthalmus marginatus.

[i] - N = sample size; A = number of alleles; Ae = allelic richness; Ho = observed heterozygosity; He = expected heterozygosity second Nei (1973); SD = Standard deviation.

Significant values of Fst were observed in all comparisons including Abaetetuba and the comparison of Tucuruí × Conceição (Table 4) gave negative Rst values representing interactions where the variance within a population exceeds the variance between populations.

Conceição****0.28 (0.000)0.05 (0.05)0.04 (0.09)0.10 (0.001)
Abaetetuba0.71 (0.0000)****0.24 (0.00)0.28 (0.00)0.27 (0.000)
Itupiranga0.35 (0.0006)0.26 (0.0061)****0.005 (0.41)0.006 (0.350)
Cametá0.20 (0.0100)0.47 (0.0002)0.04 (0.19)****0.008 (0.350)
Tucuruí0.41 (0.0000)0.42 (0.0000)-0.02 (0.63)0.02 (0.223)****

Table 4.

Comparisons of peer to peer between the populations of H. marginatus analyzed, with FST values (above the diagonal) calculated [21]according to Weir and Cockerham (1984), with their respective P‐values (in quotes) and values of Rst (below the diagonal) calculated according to Michalakis and Excoffier (1996) with their respective P‐values (in quotes).

4. Discussions

Although tested just by two polymorphic SSRs, genetic variability of H. marginatus population is under the influence of the Tucuruí Dam on the Tocantis River and as measured by heterozygosity test it is relatively low, regarding the number of alleles (absolute and effective). Few studies use less than six microsatellite markers in population studies [30]. Although we would like to stress that the results presented here are preliminary, they are consistent with results obtained with other molecular markers. The importance of our work is that it described a library that can be used to test for the polymorphism kind of economic importance to the Amazon.

Some populations of H. marginatus, such as those of Abaetetuba and Conceição, had remarkable low values. We could indicate that the heterozygosity was similar to that shown by other Siluriformes [31].

A comparison made elsewhere of H. marginatus cytochrome b gene with the same genes in other freshwater fishes [32, 33] indicates very low genetic diversity levels in H. marginatus, possibly reflecting a low mutation rate in this species [34] or a characteristic of Pimelodidae, as found in other studies [35].

The heterozygosity values provided by two microsatellite loci in H. marginatus in five populations analyzed did not differ between populations upstream and downstream of Tucuruí Dam. In these populations, the observed heterozygosity was lower than expected heterozygosity for all populations, thus indicating that there is no evidence of population bottleneck. On the other hand, the distribution of heterozygosity did not vary much in comparison with other species of commercial fish such as arowana (Osteoglossum bicirrhosum (Vandelli 1829)) [15] or pirarucu (Arapaima gigas) [16].

There is no much information on populations of H. marginatus, or other Hypophthalmus, as this knowledge would have been instrumental in understanding the effects of the Tucuruí Dam on the populations of H. marginatus Tocantins and Araguaia. The values of Fst and Rst showed little differentiation among populations of the same side of the current course of the Tocantins River, for example, populations of Conceição and Itupiranga upstream of Tucuruí Dam or Abaetetuba and Cametá downstream of the same dam gave values of Fst low, as between populations separated by the dam as Abaetetuba and Conceição, accented Fst values indicating greater differentiation among populations, probably generated by low levels of gene flow.

Laroche and Durand [36] studied the genetic structure of the Percidae Zinger asper populations, an endangered endemic species, affected by the construction of a dam built on the River Rhone in France, and found significant genetic differences between upstream and downstream populations.

5. Preliminary conclusions

If low differentiation is prevalent, the geographic distribution range of large samples would be recommended for genetic analyses. It would also be useful to assess the levels of genetic variability within and among populations from different basins for a better understanding of population dynamics of these species. Although our conclusions were made by using only two microsatellite loci analyses, results are consistent with data from mitochondrial markers.


This study was supported by CAPES, and CNPq. We are also grateful to IBAMA, for providing the authorization 013‐2007 to capture, collect, and transport the biological material to the Centrais Elétricas do Norte do Brasil S.A. (ELETRONORTE S.A.) as well as for logistic support at Tucuruí. We thank Cassio Andrade (EMATER/PA) and Comunidade Jaraqüera Grande for providing logistic support at Cametá; Aparecida, Jhonis, and Leo for providing logistic support at Conceição de Araguaia. Special thanks to Soraya Andrade, Silvanira Ribeiro Barbosa for helping in laboratory. EJHR thanks Mauricio Papa de Arruda for collaboration with ALFexpress management.


1 - Barthem R, Lambert M, Petrere M. Life strategies of some long‐distance migratory catfish in relation to hydroelectric/dams in the Amazon Basin. Biological Conservation. 1991;55:339–345. DOI:10.1016/0006‐3207(91)90037‐A.
2 - Wei Q, Ke F, Zhang J, Zhuang P, Luo J, Zhou R, Yang W. Biology, fisheries, and conservation of sturgeons and paddlefish in China. Environmental Biology of Fishes. 1997;48:241–255. DOI: 10.1007/0‐306‐46854‐9_14.
3 - Ruban GI. Species structure, contemporary distribution and status of the Siberian sturgeon, Acipenser baerii. Environmental Biology of Fishes. 1997;48:221–230. DOI: 10.1007/0‐306‐46854‐9_12.
4 - Agostinho AA, Júlio HF Jr, Borghetti JR. Considerations on the impacts of dams on fish populations and measures for mitigation - a case study: Itaipu reservoir. Revista UNIMAR. 1992; 14:89-107
5 - Fearnside PM. Dams in the Amazon: Belo Monte and Brazil's hydroelectric development of the Xingu River Basin. Environment Management. 2006;38(1):16–27. DOI:10.1007/s00267‐005‐0113‐6.
6 - Carvalho JL, M豯na B. Studies on two migratory fish from lower Tocantins River, before closure of Tucuru럤am. Amazoniana. 1986;9:595-607.
7 - Mérona B, Carvalho JL, Bittencourt MM. The immediate effects of the closure dam Tucuruí (Brazil) on the downstream ichthyofauna. Revue d'Hydrobiologie Tropicale. 1987;20:73?84.
8 - Ribeiro MCLB, Petrere M, Juras AA. Ecological integrity and fisheries ecology of the Araguaia-Tocatins River Basin, Brazil. Rivers Research and Applications. 1995;11:325?350. DOI: 10.1002/rrr.3450110308.
9 - Cetra M, Petrere M. Small-scale fisheries in the middle River Tocantins, Imperatriz (MA), Brazil. Fisheries Management and Ecology. 2001;8:153-162.
10 - Santos GM, Mérona B, Juras AA, Jégu M. Fish and fishing in the Lower Tocantins River twenty years after Tucuruí Hydroelectric Plant . ELECTRONORTE. Bras쫩a. 2004. 216 p.
11 - Leite RAN, Bittencourt MM. Hydropower impact on the Amazon fish fauna: The sample Tucuruí. In: Val AL, Figiuolo R and Feldberg E, editors. Scientific Basis for Preservation Strategies and Development of the Amazon: Facts and Perspectives Vol 1. Instituto Nacional de Pesquisas da Amazônia (INPA). Manaus; 1991. p. 85-100.
12 - Odinetz-Collart O. The shrimps Macrobrachium amazonicum (Palaemonidae) in the Lower Tocantins, after closing the dam of Tucuru램Brazil).Revue d'Hydrobiologie Tropicale. 1987;20(2):131-144.
13 - Odinetz-Collart O. Ecology and fishing potential of shrimp cinnamon, Macrobrachium amazonicum, in the Amazon Basin In: Ferreira EJG, GM Santos, Leão ELM, Oliveira LA (eds) Scientific Basis for Preservation and Development of the Amazon Strategies, Vol 2 Instituto Nacional de Pesquisas da Amaz?nia (INPA), Manaus; 1993. p. 147-166.
14 - Sambrook J, EF Fritsch, Maniatis T. Molecular Cloning: A Laboratory Manual. 2nd edn. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, USA. 1989.
15 - Silva T J, Hrbek T, Farias I P. Microsatellite markers for the silver arowana (Osteoglossidae, Osteoglossiformes). Molecular Ecology Notes 2009;9:1019–1022. DOI: 10.1111/j.1755‐0998.2009.02556.x.
16 - Farias IP, Hrbek T, Brinkmann H, Sampaio I, Meyer A. Characterization and isolation of DNA microsatellite primers for Arapaima gigas, an economically important but severely over‐exploited fish species of the Amazon basin. Molecular Ecology Notes. 2003;3:128–130. DOI:10.1046/j.1471‐8286.2003.00375.x.
17 - Refseth UH, Fangan BM, Jakobsen KS. Hybridization capture of microsatellites directly from genomic DNA. Electrophoresis. 1997;18:1519–1523. DOI: 10.1002/elps.1150180905.
18 - Hall TA. Bioedit v709: Biological sequence alignment editor analysis program for Windows 95/98/Nt. Nucleic Acids Symposium Series 2007;41:95–98.
19 - Rozen S, Skaletsky HJ. Primer 3 on the www for general users and for biologist programmers. In: Krawetz S and Misener S, editors. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Humana Press, Totowa; NJ. 2000. p. 365–386.
20 - Schuelke M. An economic method for the fluorescent labeling of PCR fragments: A poor man's approach to genotyping for research and high‐throughput diagnostics. Nature Biotechnology 2000;18:233–234. DOI:10.1038/72708.
21 - Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICROCHECKER: software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes. 2004;4:535–538. DOI: 10.1111/j.1471‐8286.2004.00684.x.
22 - Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources. 2010;10:564–567. DOI: 10.1111/j.1755‐0998.2010.02847.x.
23 - Yeh FC, Rong‐cai Y, T Boyle. POPGENE VERSION 1.31: Microsoft Window‐based Freeware for Population Genetic Analysis. University of Albert, Centre for International Forestry Research. 1999.
24 - Raymond M, Rousset F. Genepop Version 1.2.: Population genetics software for exact tests and ecumenicism. Journal of Heredity. 1995;86:248–249.
25 - Eggert LS, Mundy NI, Woodruff DS. Population structure of loggerhead shrikes in the California Channel Islands. Molecular Ecology. 2004:13:2121–2133. DOI: 10.1111/j.1365‐294X.2004.02218.x.
26 - Wright S. Evolution & genetics for populations, vol.4. Variability Within & Among Natural Populations. University of Chicago Press, Chicago. 1978.
27 - Slatkin MA. Measure of population subdivision based on microsatellite allele frequencies. Genetics. 1995;139:457–462.
28 - Goodman SJ. R‐ST Calc: a collection of computer programs for calculating estimates of genetic differentiation from microsatellite data & determining their significance. Molecular Ecology. 1997;6:881–885. DOI: 10.1111/j.1365‐294X.1997.tb00143.x.
29 - Jarne P, Lagoda PJL. Microsatellites, from molecules to populations & back. Trends in Ecology & Evolution. 1996;11:424–429.
30 - Ji YJ, Zhang DX, Hewitt GM, Kang L, Li DM. Polymorphic microsatellite loci for the cotton bollworm Helicoverpa armigera (Lepidoptera: Noctuidae) and some remarks on their isolation. Molecular Ecology Notes. 2003;3:102–104. DOI: 10.1046/j.1471‐8286.2003.00366.x.
31 - Revaldaves E, Pereira LHG, Foresti F, Oliveira C. Isolation and characterization of microsatellite loci in Pseudoplatystoma corruscans (Siluriformes: Pimelodidae) and cross‐species amplification. Molecular Ecology Notes. 2005;5:463–465. DOI: 10.1111/j.1471‐8286.2005.00883.x.
32 - Sivasundar A, Bermingham E, Ortí G. Population structure and biogeography of migratory freshwater fishes (Prochilodus: Characiformes) in major South American rivers. Molecular Ecology. 2001;10:407–418.
33 - Turner TF, Mcphee MV, Campbell P, Winemiller KO. Phylogeography and intraspecific genetic variation of prochilodontid fishes endemic to rivers of northern South America. Journal of Fish Biology. 2004;64:186–201. DOI: 10.1111/j.1095‐8649.2004.00299.x.
34 - Hernández‐Ruz EJ, Goncalves EC, Silva A LC, Schneider MPC. Low genetic diversity of Hypophthalmus marginatus from the Tocantins River based on cytochrome b sequence data. International Journal of Genetic and Molecular Biology. 2013;5(6):71–77. DOI: 10.5897/IJGMB2013.0071.
35 - Du M, Liu YH, Niu BZ. Isolation and characterization of polymorphic microsatellite markers in Bagarius yarrelli using RNA‐Seq. Genetics and Molecular Research. 2015;14(4):16308–16311. DOI: 10.4238/2015.
36 - Laroche J, Durand JD. Genetic structure of fragmented populations of a threatened endemic percid of the Rhône River: Zingel asper. Heredity. 2004;92:329–334. DOI: 10.1038/sj.hdy.6800424.