Acta Zoologica (Stockholm) 91: 242–248 (April 2010) doi: 10.1111/j.1463-6395.2009.00404.x

A prime inference on genetic diversity (RAPDs) in the marine fish Atherinella brasiliensis (Teleostei, Atherinopsidae) from Southern Brazil Maria Cristina da Silva Cortinhas1, Chirlei Glienke1,AlbertoJose´ Prioli2,3, Rafael Bueno Noleto1, Daniele Aparecida Matoso1 and Marta Margarete Cestari1

Abstract 1Laborato´rio de Citogene´tica , DaSilvaCortinhas,M.C.,Glienke,C.,Prioli,A.J.,Noleto,R.B.,Matoso, Departamento de Gene´tica, Setor de D. A. and Cestari, M. M. 2010. A prime inference on genetic diversity Cieˆncias Biolo´gicas, Universidade Federal (RAPDs) in the marine fish Atherinella brasiliensis (Teleostei, Atherinopsidae) do Parana´, Centro Polite´cnico, Curitiba, from Southern Brazil. — Acta Zoologica (Stockholm) 91: 242–248 Parana´, Brazil, CEP 81531-990; 2Departamento de Biologia Celular e As a result of the importance of Atherinella brasiliensis in estuarine environ- Gene´tica, Universidade Estadual de ments, random amplified polymorphic DNA (RAPD) markers were used to Maringa´,Maringa´,Parana´, Brazil, CEP verify the genetic diversity in A. brasiliensis from two different places in Paran- 87020-900; 3Nu´cleo de Pesquisas em agua´ Bay (Parana´ State) and one from the Conceic¸a˜o Lagoon (Santa Catarina Limnologia, Ictiologia e Aqu¨icultura (NUPELIA), Universidade Estadual de State). Cytogenetic data have shown a high karyotypic diversity in some popu- Maringa´,Maringa´,Parana´, Brazil, lations, although in others this peculiarity demonstrates rearrangements such as 87020-900 heterochromatinization. In the present study, a low level of genetic structuring between the samples from Conceic¸a˜o Lagoon compared with the others was Keywords: observed through principal coordinate analysis (PCO), analysis of molecular Atherinella brasiliensis, silversides, genetic variance and Mantel test according to 79 RAPD markers. As this specie does variability, conservation genetics, RAPD not perform horizontal migration and the individuals of Conceic¸a˜o Lagoon are isolated, three hypotheses are proposed to explain the results: (i) similar envi- Accepted for publication: ronments may show homogeneous populations not depending on the geo- 17 February 2009 graphical distance, (ii) because vicariant events that formed the bays occurred in a recent period, the fragmentation effects over the structuring of the genetic diversity may still be low and not totally detectable by the RAPD technique and (iii) the isolation time or the number of generations may not be enough to promote a possible differentiation and genetic structuring between the speci- mens of these three places. The specimens of these places present a low level of differentiation and genetic structuring so we can consider them as a unique homogeneous population. MariaCristinadaSilvaCortinhas,Laborato´rio de Citogene´tica Animal, Departamento de Gene´tica, Setor de Cieˆncias Biolo´gicas, Universidade Federal do Parana´, Centro Polite´cnico, Jardim das Ame´ricas, Curitiba, Parana´, Brazil, CEP 81531-990. E-mail: [email protected]

In Brazil, the king-fish of the Atherinopsidae family are Introduction represented by the Odontesthes and Atherinella genera. Within The order is a monophyletic group known as the Atherinella , Atherinella brasiliensis (Quoy and silversides or king-fish, and it is composed of six families Gaimard 1824) can be found from Pernambuco State (among them Atherinopsidae) and 49 genera (Dyer 2006). (Franc¸a 2005) to Rio Grande do Sul State (Bemvenuti These fish are distributed in marine, fresh and estuarine 1987), mainly in coastal waters and especially in river waters and in tropical and temperate environments all around mouths and brackish water regions (Figueiredo and Menezes the world (Beheregaray and Sunnucks 2001). 1978).Thisfishhasbeenwidelyusedbothasasourceof

2009 The Authors 242 Journal compilation 2009 The Royal Swedish Academy of Sciences Acta Zoologica (Stockholm) 91: 242–248 (April 2010) da Silva Cortinhas et al. • Genetic population structure Atherinella brasiliensis complementary wealth and even as a food (Hostim-Silva et al. used as biomarkers of loss or gain of genetic diversity in 1995) because it is a resident estuarine species (Pessanha organisms that live in contaminated environments (Wolf et al. and Arau´jo 2001) and is constant and abundant in these 2004), and are also useful in genotoxicity and carcinogenesis environments (Paiva Filho and Gianinni 1990; Pessanha et al. studies (Atienzar and Jha 2006). The threat to A. brasiliensis is 2000; Vendel et al. 2000). of great concern because the biological diversity involves not Da Silva Cortinhas et al. (2003) performed a cytogenetic only the variety of species, but also the intraspecific genetic study in A. brasiliensis from the Paranagua´ Bay (Parana´ State) variation. and Conceic¸a˜o Lagoon (Santa Catarina State) detecting Atherinella brasiliensis does not perform horizontal migration intra- and interpopulational chromosomal variation, probably (Hostim-Silva et al. 1995) living only in estuary environments. determined by pericentric inversions. Most of the analysed This feature represents a barrier to gene flow, increasing the individuals presented a diploid number of 48 chromosomes chances that genetic differences may accumulate with time. It constituted by four metacentric, 10 submetacentric, 14 sub- is relevant to point out that both contaminated environments telocentric and 20 acrocentric chromosomes. The pattern and non-horizontal migration may act as independent factors karyotype had a fundamental number equal to 76. However, that might promote differentiation and a possible genetic five different cytotypes were found in the specimens from structuring among fish from these areas. Paranagua´ Bay (Parana´ State) and only one different cytotype Considering the possible isolation of A. brasiliensis popula- in the specimens from the Conceic¸a˜o Lagoon (Santa Catarina tions because of the lack of horizontal migration, and also the State). contamination observed in the areas inhabited by this species, According to Sczepanski et al. (2007), variations in hetero- the diagnosis and interpretation of the genetic variation in this chromatinization of short arms may produce hetero- species is extremely important. This knowledge will support morphisms within some pairs, explaining the differences in protection and conservation programmes for the maintenance the fundamental number in this group, although for Sola of the diversity of this species. et al. (1988) this feature in Atherinopsidae species is partially related to the activity of nucleolar organizing Materials and Methods region (NORs) as well as to structural differences in hetero- chromatic regions. Sample collection For a better understanding of the genetic diversity of A. brasiliensis from Pontal do Parana´ and Laranjeiras Bay Atherinella brasiliensis was collected in two different places in (Parana´ State) and Conceic¸a˜o Lagoon (Santa Catarina Parana´ State; in Pontal do Parana´,attheentranceofParanag- State) the random amplified polymorphic DNAs (RAPDs) ua´ Bay (2534¢ S, 4821¢ W), 30 specimens were collected, were chosen because of their wide use in population genetics and 29 were obtained in Laranjeiras Bay (2535¢ S/ and the availability of documentation for many different 4819¢ W), which is connected to the Paranagua´ Bay. organisms, including fish (Dergam et al. 1998). This tech- Another 30 specimens were collected in Conceic¸a˜oLagoon nique allows the detection of a large number of polymor- (2735¢ S, 4815¢ W), city of Floriano´polis, Santa Catarina phismsinDNAsegmentsdistributed throughout the genome State (Fig. 1). of the studied individuals, mainly because of the quantity of A sample of approximately 100 mg of the dorsal muscle available primers used to detect variations in some regions was removed from each sample. The muscular tissue was of the genome (Bardakci and Skibinski 1994). This analysis stored in ethanol 96% and preserved in a freezer at –20 C. may be more sensitive than traditional methods because it involves non-coding neutral regions, which allow high muta- DNA extraction and purification tion levels (Cagigas et al. 1999). This marker has been very useful in the management of neotropical fish in places The total DNA was extracted by the phenol/chloroform affected by human activities, relating genetic variability data method according to Sambrook et al. (1989) with some modi- with environmental parameters (Hatanaka and Galetti fications. The sample of the muscle was macerated in liquid 2003; Matoso et al. 2004; Sofia et al. 2006) and identifying nitrogen and homogenized in 500 lLPSbuffer[Tris–HCl population structuring (Prioli et al. 2002; Almeida et al. 0.2 M pH 9.0, ethylenediaminetetraacetic acid (EDTA) 2003; Leuzzi et al. 2004; Torres and Ribeiro 2009). 30 mM pH 8.0, sodium dodecyl sulphate 2% and saccharose The contamination of some areas where A. brasiliensis is 5%], 500 lLTHbuffer(Tris–HCl10mM pH 7.5, NaCl found is harming the populations, and loss of genetic varia- 60 mM,EDTA10mM pH 8.0, saccharose 5%, spermine tion, as observed in Paranagua´ Bay (Parana´ State) near 0.15 mM and spermidine 0.15 mM)and6lL K proteinase Paranagua´ Seaport, is one of the problems. According to (20 mg/mL) and incubated in a water-bath for 90 min at Wolf et al. (2004), toxic compounds may interfere in the 50 C. Afterwards, the DNA was purified by extraction with genetic constitution of populations, either directly through phenol/chloroform (1 : 1, v : v) and then chloroform, after mutagenic activity, or indirectly via population-mediated which it was precipitated with saline solution (NaCl 5 M)and processes (i.e. selection, bottlenecks). The RAPDs can be 100% cold ethanol followed by an overnight incubation at

2009 The Authors Journal compilation 2009 The Royal Swedish Academy of Sciences 243 Genetic population structure Atherinella brasiliensis • da Silva Cortinhas et al. Acta Zoologica (Stockholm) 91: 242–248 (April 2010)

Fig. 1—Location where Atherinella brasiliensis was collected: 1, Conceic¸a˜o Lagoon (SC); 2, Pontal do Parana´ (PR);3,LaranjeirasBay(PR).

–20 C. The pellet was re-suspended in 60 lL1/10TEbuffer 1.4% agarose gel submerged in 1 · TBE buffer for 4 h at (Tris–HCl 10 mM pH 7.5, EDTA 1 mM pH 8.0) containing 60 V. Afterwards, the products were stained with ethidium RNAse 10 mg/mL and incubated in a water-bath at 37 Cfor bromide (0.05%), visualized under ultraviolet light and 90 min. The quantity of DNA in each sample was estimated photographed. Each band observed in the gel was considered through the comparison with the phage k DNA, for which the a Mendelian segregation locus and a binomial matrix (1/0) concentration was previously calculated using electrophoresis was generated based on presence/absence of bands, in 0.8% agarose gel for 90 min at 100 V. Afterwards, the respectively, with the aim of calculating the intra- and inter- DNA was stained with ethidium bromide (0.05%), visualized populational diversities. under ultraviolet light and photographed using an Olympus The similarity matrices were generated through the NTSYS camera. software version 2.1 (Rohlf 2000), adopting the Jaccard (1901) coefficient which plotted a principal coordinate analy- sis (PCO) to verify the dispersion of the collected individuals. RAPD and data analysis Using the same software, the cophenetic correlation was per- After the tests to accurately determine the concentration of formed to compare the two symmetric matrices, and grouping the main reagents of the polymerase chain reaction were per- analyses were obtained by the UPGMA method (Sneath and formed, a band pattern was obtained for each primer and the Sokal 1973). Bootstrap values for the UPGMA method were reactions were repeated several times to ensure the reliability obtained through the BOOD software version 3.03 (Coelho of the results. 2004). To confirm the results obtained by the UPGMA method, The amplification of the RAPD bands was performed grouping analyses through the neighbour-joining method according to the Bardakci and Skibinski (1994) conditions, were also performed using the software FREETREE (Pavlicek with a reaction mixture of 13 lL total volume containing: et al. 1999). The Mantel test (1967) was used through the 10 ng DNA, 0.46 lM primer, 1 U Taq polymerase, 2 mM MANTEL STRUCT 1.0 software (Miller 1999) to evaluate the

MgCl2,200lM dNTP, Tris–KCl buffer (Tris–HCl 20 mM genetic differentiation among the specimens of A. brasiliensis, pH 8.0 and KCl 50 mM) and ultra-pure water enough to through the comparison between the similarity and the model complete the volume. The Techne thermocycler was used matrix that expresses the origin of the individuals from pair- and programmed according to the conditions proposed by Al- wise comparisons between the groups. meida (1998), with an initial stage of 4 min at 92 C, followed A genetic distance Euclidean matrix was used to generate by 40 cycles of 1 min at 92 C, 90 s at 40 Cand2minat an analysis of molecular variance (AMOVA) according to 72 C. The final stage was 72 C for 5 min, followed by cool- Excoffier et al. (1992) through the ARLEQUIN software version ing at 20 Cfor20minand4Cfor30min. 2000 (Schneider et al. 2000) to estimate a possible genetic Twenty-fourprimersoftheOPA,OPE,OPK,OPXand structuring among the individuals from the places studied.

OPW (Operon Technologies Inc., Alameda, CA, USA) kits This analysis provides the hST, which is analogous to the were tested for the oligonucleotide selection used in the pres- Wright (1978) FST, which represents one measure of the ent report, but only 12 were chosen based on repeatability, genetic differentiation degree or interpopulational sub- band clear pattern and presence of polymorphic bands. The division (Slatkin 1993). According to Wright (1978), to obtain amplification products were submitted to electrophoresis in a significant estimative different from zero, results between 0

2009 The Authors 244 Journal compilation 2009 The Royal Swedish Academy of Sciences Acta Zoologica (Stockholm) 91: 242–248 (April 2010) da Silva Cortinhas et al. • Genetic population structure Atherinella brasiliensis

Table 1 Random amplified polymorphic DNA (RAPD) primers, Both the phenogram obtained through the neighbour-join- corresponding to sequences and number of bands associated with ing method and the one generated through the UPGMA each primer method (not shown) presented low bootstrap values. Besides, the phenogram obtained through the UPGMA method Number of bands presented a low cophenetic correlation (0.69), suggesting a short correspondence between the generated matrices Primer Sequence (5¢–3¢) n n t p through the phenogram and the matrix obtained through the OPW03 GTCCGGAGTG 5 4 observations. We can therefore infer that the generated phen- OPW04 CAGAAGCGGA 8 6 ograms were not the best representation of the genetic diver- OPW06 AGGCCCGATG 7 3 sity pattern among the specimens of the different places OPW09 GTGACCGAGT 7 5 studied. OPW17 GTCCTGGGTT 8 5 ThePCO,basedontheJaccardsimilaritymatrix,wasthe OPX04 CCGCTACCGA 13 11 best indicator of genetic diversity among the studied individu- OPX07 GAGCGAGGCT 6 4 OPX17 GACAGGGACC 10 9 als (Fig. 3). Through the PCO it was possible to observe that OPX20 CCCAGCTAGA 3 2 the specimens from Pontal do Parana´ and Laranjeiras (Parana´ OPA06 GGTCCCTGAC 5 3 State) comprised a homogeneous population. This hypothesis OPE09 CTTCACCCGA 3 2 wasconfirmedbytheAMOVA results and by the Mantel test OPK14 CCCGCTACAC 4 1 (Tables 2 and 3), with significant P-values. Because these two Total 79 55 regions are very close, gene flow between the individuals of both places is likely. The total number of scorable bands and the number of polymorphic fragments are represented by nt and np, respectively. Comparing the specimens from Conceic¸a˜oLagoon(Santa Catarina State) and from Parana´ State, there was an increase in the genetic diversity (r ¼ 0.1627 and r ¼ 0.1530, respec- tively) observed through the Mantel test (Table 3), and also and 0.05 indicate low genetic structuring; between 0.05 and in the genetic structuring (hST ¼ 0.07519 and hST ¼ 0.25 indicate moderate genetic structuring; and over 0.25 0.07432, respectively) observed through AMOVA (Table 2), indicate a high genetic structuring. with highly significant P. According to Wright (1978) it is pos- sible to deduce the occurrence of incipient genetic structuring among the individuals of this area compared with others previ- Results ously analysed. However, the results are still too few to geneti- Seventy-nine bands with variable polymorphism between cally differentiate the specimens of the three places, so they 2000 and 300 base pairs (bp) (Table 1) were obtained, can be considered as a unique population. which indicates a high intrapopulational genetic variation compared with the interpopulational relations (Table 2). Discussion Figure 2 shows the electrophoretic profile obtained with the OPX07 primer and the polymorphic band variation is indi- The interpopulational genetic variation levels of marine fish cated by the arrow (~600 bp). This figure corroborates the are, in theory, lower than those of freshwater fish. This finding results observed for other primers, which showed that the seems to be associated with the population structures of the variability occurred predominantly within the areas studied marine species, which present a less restricted gene flow; and not between them. also, the formation of geographical barriers is much less

Table 2 Analysis of molecular variance (AMOVA –Excoffieret al., 1992) with P values obtained from 1023 permutations

Variation

Locations Source Percentage hST P-value

Pontal do Parana´ (PR) · Laranjeiras (PR) Among loc. 2.02 0.02016 0.01173 Within loc. 97.98 Pontal do Parana´ (PR) · Conceic¸a˜o Lagoon (SC) Among loc. 7.52 0.07519 0.00000 Within loc. 92.48 Laranjeiras (PR) · Conceic¸a˜o Lagoon (SC) Among loc 7.43 0.07432 0.00000 Within loc. 92.57

Loc., localities.

2009 The Authors Journal compilation 2009 The Royal Swedish Academy of Sciences 245 Genetic population structure Atherinella brasiliensis • da Silva Cortinhas et al. Acta Zoologica (Stockholm) 91: 242–248 (April 2010)

Fig. 2—Profile of random amplified polymorphic DNA electrophoresis of Atherinella brasiliensis from Pontal do Parana´ (Parana´ State), Laranjeiras (Parana´ State) and Conceic¸a˜o Lagoon (Santa Catarina State) showing the amplified fragments by primer OPX 07. The first and the last columns contain the molecular weight marker fragments (a ladder of 100 base pairs, PB-L Productos Bio-Lo´gicos), presenting a first additional band of 1434 base pairs. Lanes 1 to 10, the specimens from Pontal do Parana´ (Parana´ State); lanes 11 to 19, the specimens from Laranjeiras (Parana´ State); and lanes 20 to 29, the specimens from Conceic¸a˜o Lagoon (Santa Catarina State).

common than in continental waters (Ward et al. 1994). How- ever, with the increasing number of studies using molecular markers in marine fish, this observation may not be a general rule and the genetic structure of these populations may depend on several other factors such as ecology, demography, behaviour, genetics, oceanography, climate and tectonism (Affonso 2004). In the samples of the Pontal do Parana´ (Parana´ State) area, variations in the intensity, loss or gain of bands were not detected among the 79 analysed RAPD markers compared with fish of other areas. This finding suggests that ecological factors (i.e. polluted environment) did not interfere in the alteration of the genetic variability in the fish from this area. The 55 polymorphic loci occurred in the individuals from the three areas showing high intrapopulational variability and a low interpopulational variability. The results of the present report are in accordance with other cytogenetic data available in the literature (Sola et al. 1988; Sczepanski et al. 2007), in Fig. 3—Principal coordinate analysis based on Jaccard genetic which the karyotypic diversity in Atheriniformes is related to similarity (1901). The specimens from Pontal do Parana´ (Parana´ heteromorphisms within some pairs explaining the differences State) are represented in white; the specimens from Laranjeiras in the fundamental number. (Parana´ State) are in grey and the specimens from Conceic¸a˜oLagoon Fish, in common with most , often show a (Santa Catarina State) are in black. subdivision of each species into local populations which are isolated to varying degrees. Among freshwaters, populations inhabiting isolated lakes or rivers with physical barriers Table 3 Mantel test (1967) between are reasonably easy to define, whereas populations of marine fish are usually less clearly delineated. If localized Populations rPR2 populations inhabit similar environments or remain interconnected through migration and gene flow, they may Pontal do Sul (PR) x Laranjeiras (PR) 0.0394 0.0188 0.0016 display largely homogeneous arrays of phenotypic or genetic Pontal do Sul (PR) x Lagoa da Conceic¸a˜o (SC) 0.1627 0.0001 0.0265 traits (Carvalho 1993). This fact can be observed in the Laranjeiras (PR) x Lagoa da Conceic¸a˜o (SC) 0.1530 0.0001 0.0234 samples of Pontal do Sul and Laranjeiras Bay (Parana´ State) Correlations (r), respective Type I (P) and determination coefficients (R2), where the gene flow is intense and, because of the proximity obtained in with 10.000 random permutations. of these areas, they tend to become homogeneous and, as a

2009 The Authors 246 Journal compilation 2009 The Royal Swedish Academy of Sciences Acta Zoologica (Stockholm) 91: 242–248 (April 2010) da Silva Cortinhas et al. • Genetic population structure Atherinella brasiliensis result, they follow similar evolutionary paths. Regarding the References samples of Conceic¸a˜o Lagoon, where the gene flow is Affonso, P. R. A. M. 2004. Marcadores moleculares na ana´lise de espe´cies restricted, incipient structuring could be explained probably ecomposic¸a˜o populacional de peixes marinhos de recifes de corais da by the environmental similarity that selects individuals with Famı´lia Pomacanthidae (Perciformes). Tese (Doutorado em Gene´tica similar characteristics. eEvoluc¸a˜o). Centro de Cieˆncias Biolo´gicas e da Sau´de, Universid- Our data concerning RAPD markers indicate three possible ade Federal de Sa˜o Carlos, 158 pp. hypotheses: Almeida, F. S. 1998. Ana´lise da variabilidade gene´tica em Pimelodidae e (i) According to Chapman et al. (1999) there is no relation- Rhamdiidae (Pisces – Siluriformes) da bacia do rio Tibagi. Dissertac¸a˜o de Mestrado. Universidade Estadual de Londrina, 102 pp. ship between genetic differences and geographical distances Almeida, F. S., Sodre´,L.M.K.andContel,E.P.B.2003.Popula- among sites. What might be happening is that similar environ- tion structure analysis of Pimelodus maculatus (Pisces, Siluriformes) ments, even geographically distant, may show homogeneous from the Tieteˆ and Paranapanema Rivers (Brazil). – Genetics and populations. Molecular Biology 26: 301–305. (ii) The fluctuations that occurred in the relative sea level Angulo,R.J.,Lessa,G.C.andSouza,M.C.2006.Acriticalreview about 5150 years ago (Suguio et al. 2005), associated with cli- of mid- to late-Holocene sea-level fluctuations on the eastern Bra- zilian coastline. – Quaternary Science Reviews 25: 486–506. matic changes during the Quaternary (Angulo et al. 2006), Atienzar, F. A. and Jha, A. N. 2006. The random amplified polymor- played an important role in the formation of the sedimentary phic DNA (RAPD) assay and related techniques applied to geno- plains of the Brazilian coast. The transformations that led to toxicity and carcinogenesis studies: a critical review. – Mutation the formation of bays occurred in a recent period, so the frag- Research 613: 76–102. mentation effects on genetic diversity structuring may still be Bardakci, F. and Skibinski, D. O. F. 1994. Application of the RAPD low and not easily detectable by the RAPD technique. technique in Tila´pia fish species and subespecies identification. – (iii) The isolation time and/or the number of generations Heredity 73: 117–123. Beheregaray, L. B. and Sunnucks, P. 2001. Fine-scale genetic struc- may not be sufficient to promote a possible differentiation and ture, estuarine colonization and incipient speciation in the marine genetic structuring between the specimens of these three silverside fish Odontesthes argentinensis.–Molecular Ecology 10: places. 2849–2866. The results of the present report are very important because Bemvenuti, M. A. 1987. Abundaˆncia, distribuic¸a˜oereproduc¸a˜ode they are a first approach to the genetic diversity of A. brasilien- peixes-rei (Atherinidae) na regia˜o estuarina da Lagoa dos Patos, sis at the molecular level. Knowledge of the genetic variability RS, Brasil. – Atlaˆntica 9:5–32. is crucial because ignoring its distribution in and between the Cagigas, M. E., Vazquez, E., Blanco, G. and Sanchez, J. A. 1999. Combined assessment of genetic variability in populations of brown specimens of a species, render conservation and preservation trout (Salmo trutta L.) based on allozymes, microsatellites, and actions almost impossible. RAPD markers. – Marine Biotechnology 1: 286–296. The specimens from these locations present a low level Carvalho, G. R. 1993. Evolutionary aspects of fish distribution: of differentiation and genetic structuring so we can consider genetic variability and adaptation. – Journal of Fish Biology 43: 53– them as a unique homogeneous population. In a future 73. conservation and management study, therefore, we could Chapman, L. J., Chapman, C. A., Brazeau, D. A., McLaughlin, B. and Jordan, M. 1999. Papyrus swamps, hypoxia, and faunal diversi- keep protected specimens protected in an unpolluted fication: variation among populations of Barbus neumayeri. area by reintroducing them in these three places when – Journal of Fish Biology 54: 310–327. necessary. Coelho,A.S.G.2004.Bood – Software para avaliac¸a˜o de dendrogramas baseados em estimativas de distaˆncias/similaridades gene´ticas atrave´s do procedimento de bootstrap.Laborato´rio de Gene´tica e Genoˆmica de Conclusion Plantas – Departamento de Biologia Geral, Instituto de Cieˆncias Analysis with the RAPD marker inidcated that A. brasiliensis Biolo´gicas, Universidade Federal de Goia´s, Goiaˆnia/Goia´s/Brazil. Da Silva Cortinhas, M. C., Cestari, M. M., Swarc¸a,A.C.andFenoc- from Pontal do Parana´ and Laranjeiras (Parana´ State) and chio, A. S. 2003. First chromosome data about the silverside Atheri- from Conceic¸a˜o Lagoon (Santa Catarina State) were almost nella brasiliensis (Atheriniformes, Pisces) from the south coast of homogeneous. New studies are necessary to confirm the data Brazil. Conventional, C-NOR and CMA3 bandings and FISH obtained through this technique. studies. – Caryologia 56: 187–191. Dergam,J.A.,Suzuki,H.I.,Shibatta,O.A.,Duboc,L.F.,Ju´lio, H. F. Jr, Giuliano-Caetano, L. and Black, W. C. IV 1998. Molecular Acknowledgements biogeography of the neotropical fish Hoplias malabaricus (Erythrini- dae: Characiformes) in the Iguac¸u, Tibagi and Parana´ Rivers. – The authors would like to thank Professor Dr Sonia M. A. P. Genetics and Molecular Biology 21: 493–496. Prioli for the laboratory analyses performed at NUPELIA Dyer, B. S. H. 2006. Systematic revision of South American silver- (Universidade Estadual de Maringa´) and also for technical sides (Teleostei, Atheriniformes). – Biocell 30: 69–88. assistance; Professor Dr Alexandre Coelho (Universidade Excoffier,I.,Smouse,P.E.andQuattro,J.M.1992.Analysisof Federal de Goia´s) who critically reviewed this manuscript; molecular variance inferred from metric distances among DNA and CAPES (Coordenac¸a˜odeAperfeic¸oamento de Pessoal haplotypes application to human mitochondrial DNA restriction data. – Genetics 131: 479–491. de Nı´vel Superior) for the financial support.

2009 The Authors Journal compilation 2009 The Royal Swedish Academy of Sciences 247 Genetic population structure Atherinella brasiliensis • da Silva Cortinhas et al. Acta Zoologica (Stockholm) 91: 242–248 (April 2010)

Figueiredo, J. L. and Menezes, N. A. 1978. Manual de Peixes Marin- Identification of Astyanax altiparanae (Teleostei, Characidae) in the hos do Sudeste do Brasil. II. Teleostei (1). Museu de Zoologia da Uni- Iguac¸u River, Brazil, based on mitochondrial DNA and RAPD versidade de Sa˜oPaulo,Sa˜oPaulo. markers. – Genetics and Molecular Biology 25: 421–430. Franc¸a, E. J. 2005. Abundaˆncia distribuic¸a˜o e caracterizac¸a˜o morfolo´gica Rohlf, F. J. 2000. Numerical and Multivariate Analysis Sys- de larvas de Atherinella brasiliensis (Atherinopsidae, Atheriniformes) no tem. NTSYS-pc. Manual. Version 2.1. Applied Biostatistics Setau- estua´rio do rio Jaguaribe, Pernambuco, Brasil. Dissertac¸a˜o apresent- ket – Department of Ecology and Evolution, State University of ada ao Mestrado de Recursos Pesqueiros e Aqu¨icultura, Universid- New York, New York. ade Federal de Pernambuco, 42 pp. Sambrook, K. J., Fritsch, E. F. and Maniatis, T. 1989. Molecular Hatanaka,T.andGaletti,P.M.Jr.2003.RAPDmarkersindicatethe Cloning: a Laboratory Manual, 2nd edn. Cold Spring Harbor Labo- occurrence of structured populations in a migratory freshwater fish ratory Press, Cold Spring Harbor. species. – Genetics and Molecular Biology 26:19–25. Schneider, S., Roessli, D. and Excoffier, L. 2000. ARLEQUIN ver Hostim-Silva,M.,Clezar,L.,Ribeiro,G.C.andMachado,C.1995. 2000: A Software for Population Genetics Data Analysis. Genetics and Estrutura populacional de Xenomelaniris brasiliensis (Quoy & Biometry Laboratory Department of Anthropology and Ecology, Gaimard, 1824) (Osteichthyes-Atherinidae) na Lagoa da Con- University of Geneva, Geneva. ceic¸a˜o, SC, Brasil. – Arquivos de Biologia E Tecnologia 38: 949–960. Sczepanski,T.S.,Noleto,R.B.,Kantek,D.L.Z.,Cortinhas,M.C. Jaccard, P. 1901. Etude comparative de la distribution florale dans S. and Cestari, M. M. 2007. Classical and molecular cytogenetics une portion des Alpes et des Jura. – Bulletin del la Socie´te´ Vaudoise of Atherinella brasiliensis (Teleostei, Atherinoformes) from South des Sciences Naturelles 37: 547–579. coast of Brazil. – Journal of Fish Biology 71: 453–460. Leuzzi,M.S.P.,Almeida,F.S.andOrsi,M.L.2004.Analysisby Slatkin, M. 1993. Isolation by distance in equilibrium and nonequilib- RAPD of the genetic structure of Astyanax altiparanae (Pisces, rium populations. – Evolution 47: 264–279. Characiformes) in reservoirs of the Paranapanema River, Brazil. – Sneath, P. H. A. and Sokal, R. R. 1973. Numerical Taxonomy. Free- Genetics and Molecular Biology 27: 335–362. man, San Francisco. Matoso,D.A.,Artoni,R.F.andGaletti,P.M.Jr.2004.Genetic Sofia,S.H.,Silva,C.R.M.,Galindo,B.A.,Almeida,F.S.,Sodre, diversity of the small characid fish Astyanax sp., and its significance L. M. K. and Martinez, C. B. R. 2006. Population genetic structure for conservation. – Hydrobiologia 527: 223–225. of Astyanax scabripinnis (Teleostei, Characidae) from an urban Miller,M.P.1999.Mantel Struct: a Program for the Detection of Popula- stream. – Hydrobiologia 553: 245–254. tion Structure Through the Use of Mantel Tests. Department of Biolog- Sola, L., Natili, G. L. and Cataudella, S. 1988. Cytogenetical charac- ical Sciences – Box 5640, Northern Arizona University, Flagstaff terization of Odontesthes bonariensis (Pisces, Atherinidae) an Argenti- AZ. nae species introduced in Italy. – Genetica 77: 217–224. Paiva-Filho, A. M. and Giannini, R. 1990. Contribuic¸a˜oaoconheci- Suguio,K.,Angulo,R.J.,Carvalho,A.M.,Correˆa, I. C. S., Tomaz- mentodabiologiadopeixe-rei,Xenomelaniris brasiliensis (Quoy e elli, L. J., Willwock, J. A. and Vital, H. 2005. Paleonı´veis do mar e Gaimard, 1824) (Atherinidae) no complexo Baı´a Estuarino de San- paleolinhasdacosta.InSouza,C.R.G.,Suguio,K.,Oliveira,A. tos e Sa˜o Vicente. Brasil. – Boletim Do Instituto Oceanogra´fico 38: M.S.andOliveira,P.E.(Eds):Quaterna´rio do Brasil.HolosEdito- 1–9. ra, Ribeira˜oPreto. Pavlicek, A., Hrda, S. and Flegr, J. 1999. FreeTree – freeware pro- Torres, R. A. and Ribeiro, J. 2009. The remarkable species complex gram for construction of phylogenetic trees on the basis of distance Mimagoniates microlepsis (Characiformes: Glaucodinae) from the data and bootstrap/jackknife analysisofthetreerobustness.Appli- southern Atlantic rain forest (Brazil) as revealed by molecular sys- cation in the RAPD analysis of genus Frenkelia.–Folia Biologica tematic and population genetic analyses. – Hydrobiologia 617: 157– (Prague) 45: 97–99. 170. Pessanha, L. F. and Arau´jo, F. G. 2001. Recrutamento do peixe-rei Vendel,A.L.,Santos,C.,Nakayama,P.andSpach,H.L.2000.O Atherinella brasiliensis (Quoy e Gainard) (Atheriniformes, Atherino- uso de re´plica no estudo da ictiofauna de uma planı´ciedemare´.– pridal) na margem continental da Baı´a de Sepetiba, Rio de Janeiro, Acta Biolo´gica Paranaense 29: 177–186. Brasil. – Revista Brasileira de Zoologia 18: 1265–1274. Ward, R. D., Woordwark, M. and Skibinski, D. O. F. 1994. A com- Pessanha, A. L. M., Arau´jo, F. G., Azevedo, M. C. C. and Gomes, I. parison of genetic diversity levels in marine, freshwater, and anadro- D. 2000. Variac¸o˜es temporais e espaciais na composic¸a˜o e estrutura mus fishes. – Journal of Fish Biology 44: 213–232. da comunidade de peixes jovens da Baı´a de Sepetiba, Rio de Wolf, H. D., Blust, R. and Backeljau, T. 2004. The use of RAPD in Janeiro. – Revista Brasileira de Zoologia 17: 251–261. ecotoxicology. – Mutation Research 566: 249–262. Prioli, S. M. A. P., Prioli, A. J., Ju´lio, H. F. Jr, Pavanelli, C. S., Olive- Wright, S. 1978. Evolution and Genetics of Population. Vol. 2: The Theory ira, A. V., Carrer, H., Carraro, D. M. and Prioli, L. M. 2002. of Gene Frequencies. University of Chicago Press, Chicago, 511 pp.

2009 The Authors 248 Journal compilation 2009 The Royal Swedish Academy of Sciences