Open Journal of Genetics
Vol.4 No.5(2014), Article ID:48816,7 pages
DOI:10.4236/ojgen.2014.45032
Isolation and Characterization of 26 Microsatellite Loci for the Shortfin Silverside Fish Chirostoma humboldtianum Valenciennes 1835 (Atherinopsidae: Menidiinae) Derived from Next Generation Sequencing and Their Cross-Amplification in Central Mexican Atherinopsids
Barriga-Sosa Irene de los Angeles1*, Rosa María García-Martínez1, Jhoana Díaz-Larrea2, Oscar Adrián Lozano Garza3, Francisco Javier García-De-León3
1Laboratorio de Genética y Biología Molecular, Planta Experimental de Producción Acuícola, Departamento de Hidrobiología, Universidad Autónoma Metropolitana Iztapalapa, México D.F., México
2Departamento de Hidrobiología, Universidad Autónoma Metropolitana Iztapalapa, México D.F., México
3Laboratorio de Genéticapara la Conservación, Centro de Investigaciones Biológicas del Noroeste, Instituto Politécnico Nacional, La Paz, México
Email: *ibs@xanum.uam.mx
Copyright © 2014 by authors and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
Received 10 June 2014; revised 9 July 2014; accepted 8 August 2014
ABSTRACT
The endemic silverside fish C. humboldtianum is of great ichtyologycal, economical and cultural relevance in central Mexico and it has been suggested that it is among a group of other “peces blancos”, the most ancestral species. Here we characterized a set of 26 microsatellite loci from the species in order to further assess population and phylogeographic issues that aid in evaluating their highly impacted populations. From 58 primer pairs tested on specimens from Villa Victoria dam (Rio Balsas Basin), 26 loci positively amplified on the species and cross-primed on specimens of the closely related and threatened Poblana alchichica, C. promelas and C. riojai. Loci resolved medium and high genetic variation (mean NA = 14.3, SD = 9.38; mean HO & HE = 0.47, SD = 0.32 and 0.58, SD = 0.32, respectively) and it is expected that these polymorphic loci are also useful in studing threatened atherinopsid species.
Keywords:Endemic, Polymorphism, Villa Victoria, Poblana alchichica, C. promelas, C. riojai
1. Introduction
Up until the first half of the 20th century, atherinopsids were the main fishing product of the inland waters of Mexico and counted among the top scale fish species. The representative endemic genera of Atherinopsidae in Central Mexico are Chirostoma and Poblana, locally recognized as “charales” and “peces blancos” (silversides), C. humboldtianum included. The species is of great economic importance, because it has been used for decades as food by humans inhabiting the shores of Pátzcuaro and Chapala lakes and reservoirs in the States of Michoacán and Mexico. Although it is true that atherinopsids are of great interest as a substantial screed in epicontinental fisheries capture, it is also true that the development and growth of artisanal fisheries in these regions lack of scientific basis and are rather based on the needs to supply high quality animal protein in short term. The species is composed by two morphotypes, clearly segregated into western and eastern populations and characterized by having differences in size, and in the mean values for median lateral and predorsal scales [1] [2] ; andrecognized as the hypothetical ancestral form of the “humboldtianum group” [1] [3] . During the last decades the populations of the species have been reduced or extirpated due to the introduction of non-native fish to the western basins [1] [2] [4] . Similarly, habitat loss, pollution and overfishing have played important roles in the decline and extinction of local populations of related species (C. riojai, C. promelas, Poblana alchichica) listed in the Official Mexican Standard NOM-059-Semarnat in 2010 and in the IUCN Red List.
A recent study on the population genetics of the species and based on sequences of the mitochondrial control region, described the species as highly differentiated and with intermediate levels of diversity [5] ; and although the authors proposed the conservations of all the populations studied, they also pointed out the possible evolutionary stochastic effects and the limitation of defining conservation units only on the bases of mtDNA. Therefore, eight poymorphic microsatellite loci were developed [6] based on an enrichment method [7] in order to further assess and confirm the diversity and structure of the species. Furthermore and with the aim of developing a larger set of genetic codominant markers available to study the silversides of México, a recently diverged group of fishes [8] [9] , here we characterized 26 additional polymorphic loci for the species and using Next Generation Sequencing and succesfully cross-amplified them in closely related species with the aim of their use in assessing further population genetics and phylogeographic issues that can aid in the conservation and management of atherinopsids.
2. Materials and Methods
2.1. Microsatellite Isolation
Genomic DNA was extracted from muscle tissue of one individual of C. humboldtinuam using a salt extraction protocol [10] and sent to the Georgia Genomics Facility (University of Georgia, Athens, USA) for sequencing. The specimen was deposited in the Colección Nacional de Peces Dulceacuícolas, Escuela Nacional de Ciencias Biológicas-IPN (ENCB-IPN-P6423). DNA was sheared using a Covaris S2, and Illumina adapters were ligated on using methods derived from [11] but using adapters equivalent to Illumina TruSeq with 10 nt indexes [12] . Libraries were pooled and run on the Illumina HiSeq 2000 (paired-end 100 reads). Resulting fastq files were demultiplexed, and reads were filtered and run through the PALfinder pipeline for microsatellite identification and primer design [13] . Over 10 million reads at 100 bp were obtained, covering approximately 1.312 bases and resolving 3× genome coverage. Detected reads with microsatellites were 108, 228, of which 0.16% presented primers. Further primer sorting was as follow: 1) all loci with primers that were found more than once were deleted, particularly most of the really long (>500 bp) and short dinucleotide repeats (<24 bp); 2) removal of the longest repeats. Suitable primer design was possible in 1153 pairs with repeats of mid-length. Fifty eight pairs of primers were selected on the bases of the following parameters: a) length (18 - 30 bp), b) G-C content (40% - 60%), c) melting temperature (55˚C - 65˚C), d) 3’-stability, e) avoid hairpin or dimer formation and f) no selfpriming, and tested for polymorphism.
2.2. Sampling
Mexican government kindly issued permit number DGOPA.07343.310810.4128 to conduct this research. Organisms from Villa Victoria were obtained from commercial artisanal catches. C. riojai and P. alchichica specimens were kindly donated by Gerardo Figueroa Lucero and Héctor Espinoza, respectively. Gill lamella tissue samples from C. promelas, C. sphyraena, C. grandocule and C. jordani were obtained from the Tissue collection of the Laboratorio de Genética y Biología Molecular, Planta Experimental de Producción Acuícola, Universidad Autónoma Metropolitana Iztapalapa.
2.3. DNA Extraction, PCR Amplification and Analysis
Total DNA was isolated as previously described from 30 C. humboldtianum from Villa Victoria Dam, State of Mexico (19˚26'N, 100˚00'W) and 22 specimens of Chirostoma and Poblanaalchichica (Table 1) to test crossamplification.
Table 1. Cross-amplification of 26 microsatellite loci of Chirostoma humboldtianum in five species of Chirostoma and one Poblana species. Sample size in parenthesis.
For C. promelas ++ = 2 genotypes recorded; + = 1 genotype; for C. jordani, C. sphyraena, C. riojai, C. grandoculeand P. alchichica ++ = 3 - 4 genotype products recorded; + = 1 - 2 genotypes; - = no amplification success when n = 4 or n = 2.
PCR were performed in 10 μl reaction, containing 40 ng DNA, 1× PCR buffer (20 mM Tris-HCl, pH 8.4, 50 mM KCl), 0.2 mM dNTPs, 0.2 μM each primer, 2.5 mM MgCl2, and 0.25 - 0.35 U Taq DNA polymerase (Invitrogen, Carlsbad, CA, USA). Cycling conditions included: 95˚C 5 min, 35 cycles of 30 s at 94˚C, 30 s at the locus-specific annealing temperature (Tm) (Table 2), 40 s at 72˚C, and a final 72˚C extension for 15 min and using an iCycler (BIORAD Laboratories, Hercules, CA, USA). PCR products were visualized by silver staining in 6% polyacrylamide gels (7.5 M urea). Allele sizes were determined using 10 bp ladder (Invitrogen). Cross-species PCR amplifications required a reduction on the stringency of reactions on some loci, carried out lowering 1 or 2 grades the Tm.
Table 2. Primer sequences and characteristics of 26 microsatellite amplified in 30 individuals of Chirostoma humboldtianum (Ch) from Villa Victoria dam.
*Significant deviation from HWE (p ˂ 0.05); AN = Accession number; ǂForward primers were 5’ labeled with the M13 sequence (5’TGTAAAACGA CGGCCAGT), for further fluorescence detection.
Potential presence of null alleles was assessed in Free NA [14] . Alleles per locus (NA) and observed (HO) and expected (HE) heterozygosities were calculated using deviations from Hardy-Weinberg equilibrium, linkage disequilibrium was estimated and tested using GENEPOP 4.0.10 [15] under default parameters. Multiple hypothesis tests p values were adjusted by the False Discovery Rate (FDR) method [16] . In order to gain information on the possible utility of these loci in the closely related species assayed, the percentage of polymorphism per species (%) was obtained by accounting only those loci that resolved more than one allele per locus and wihtin each species.
3. Results and Discussion
Among 58 primers tested, 26 microsatellites screened successfully, 25 resolved as polymorphic and Chum53 was monomorphic in specimens from Villa Victoria. Mean number of alleles in the polymorphic loci was 14.3, SD = 9.38, a mean number 2× higher than the resolved earlier for the species (see [6] ; 6.4, SD = 2.5), but similar to those characterized closely related species [17] [18] ; and HO and HE 0.47, SD = 0.32 and 0.58, SD = 0.32, respectively. These later values are similar to those previously reported (see [6] ; 0.61, SD = 0.19 and 0.63, SD = 0.12, for HO and HE, respectively). Linkage disequilibrium was detected for a pair of loci within specimens from the population (Chum15 and Chum39, p = 0.00). Null allele frequency estimates were negligible for half of the loci, moderate (>0.05 and <0.20) for 12 loci (Chum11, Chum14, Chum18, Chum30, Chum33, Chum39, Chum48, Chum59, Chum63, Chum68, Chum69 and Chum70) and high at one locus (Chum27 = 0.357). Significant deviation from HWE was detected in six loci (p ˂ 0.05, Table 2). This deviation might result from the heterozygote deficiency found in some loci (Chum11, 18, 27, 30 & 39), which in turn might be caused by endogamy as has been earlier suggested for the species [5] and in other species of the genus [19] [20] .
Cross-PCR amplifications of 22 individuals of Chirostoma and Poblana resulted in 16 polymorphic loci and the remaining 10 loci resolved genotypes in most of the species (see Table 1). For instance Chum 07, 18, 33, 39, 53, 57, 60, 62, 63 and 69 resolved genotypes in 3 or 4 individuals (++) of each of the six screened species, excepted for C. promelas, for which the number of specimens analyzed was two; Chum 12, 27, 44, 48, 56 and 65, rendered genotypes (++) in five of the species and genotypes (+) in one species per locus. Chum 11, 58 and 68, resolved genotypes (++) in four species and (+) in two species; Chum 14, 30, 55 and 59, also resolved genotypes (++) in four species and (+) in another species, however these loci did not amplified neither in C. promelas or C. sphyraena; the remaining loci amplified genotypes (++ or +) in one to three species and did not amplified in other two or three species. Although eight microsatellite loci have been recently developed for the species (see [6] ), these 26 new microsatellite loci are the first to be tested in cross-amplification between the Mesa Central Atherinopsids and will be helpful tools for genetic population and phylogeographic studies of C. humboldtianum, as well as in paternal assignments for culture purposes, studies all that can aid in developing a conservation plan for this important species, which although is among the group known as “pecesblancos” the species with the widest distribution (see [19] ), some of its populations have been reduced or extirpated due to habitat loss, pollution and overfishing (see [19] ). These markers will also aid in investigating the integrity of stocks/populations for management and conservation of related endangered species.
Acknowledgements
Funding for this research was provided by Consejo Nacional de Ciencia y Tecnología-CB-2009-01-130220 and UAM.147.07.03/147.09.01 granted to IDLABS. This work was completed by Irene D. L. A. Barriga-Sosa during her sabbatical leave in CCA-UAA.
References
- Barbour, C.D. (1973) A Biogeographical History of Chirostoma (Pisces: Atherinidae): A Species Flock from the Mexican Plateau. Copeia, 3, 533-556. http://dx.doi.org/10.2307/1443118
- Barbour, C.D. (1973) The Systematic and Evolution of the Genus Chirostoma Swaison (Pisces: Atherinidae). Tulane Studies in Zoology and Botany, 18, 97-141.
- Echelle, A.A. and Echelle, A.F. (1984) Evolutionary Genetics of a Species Flock: Atherinid Fishes on the Mesa Central of México. In: Echelle, A.A. and Kornfield, I., Eds., Evolution of Fish Species Flocks, University of Maine, Orono, 93-110.
- Alvarez, J. and Navarro, L. (1957) Los peces del valle de México. Secretaría de Marina, Instituto Nacional de Pesca, Dirección de Pesca e Industrias Conexas, México City.
- García-Martínez, R.M., Mejía, O., García de León, F.J. and Barriga-Sosa, I.D.L.A. Extreme Genetic Divergence in the Threatened Endemic Fish Chirostoma humboldtianum (Valenciennes, 1835): Implications for Its Conservation.
- García-Martínez, R.M., García De León, F.J., Mejía, O. and Barriga Sosa, I.D.L.A. (2014) Isolation and Characterization of Microsatellite Loci in the “Charal de Xochimilco” Chirostoma humboldtianum. Revista Mexicana de Biodiversidad.
- Glenn, T.C. and Schable, N.A. (2005) Isolating Microsatellite DNA Loci. Methods in Enzymology, 395, 202-222. http://dx.doi.org/10.1016/S0076-6879(05)95013-1
- Bloom, D.D., Weir, J.T., Piller, K.R. and Lovejoy, N.R. (2013) Do Freshwater Fishes Diversify Faster than Marine Fishes? A Test Usingstate-Dependent Diversification Analyses and Molecular Phylogenetics of New World Silversides (Atherinopsidae). Evolution, 67, 2040-2057. http://dx.doi.org/10.1111/evo.12074
- Bloom, D.D., Piller, K.R., Lyons, J., Mercado-Silva, N. and Medina-Nava, M. (2009) Systematics and Biogeography of the Silverside Tribe Menidiini (Teleostomi: Atherinopsidae) Based on the Mitochondrial ND2 Gene. Copeia, 2, 408-417. http://dx.doi.org/10.1643/CI-07-151
- Aljanabi, S.M. and Martinez, I. (1997) Universal and Rapid Salt-Extraction of High Quality Genomic DNA for PCR-Based Techniques. Nucleic Acids Research, 25, 4692-4693. http://dx.doi.org/10.1093/nar/25.22.4692
- Fisher, S. Barry, A., Abreu, J., Minie, B., Nolan, J., Delorey, T.M., Young, G., Fennell, T.J., Allen, A., Ambrogio, L., Berlin, A.M., Blumenstiel, B., Cibulskis, K., Friedrich, D., Johnson, R., Juhn, F., Reilly, B., Shammas, R., Stalker, J., Sykes, S.M., Thompson, J., Walsh, J., Zimmer, A., Zwirko, Z., Gabriel, S., Nicol, R. and Nusbaum, Ch. (2011) A Scalable, Fully Automated Process for Construction of Sequence-Ready Human Exome Targeted Capture Libraries. Genome Biology, 12, R1. http://dx.doi.org/10.1186/gb-2011-12-1-r1
- Faircloth, B.C. and Glenn, T.C. (2012) Not All Sequence Tags Are Created Equal: Designing and Validating Sequence Identification Tags Robust to Indels. PLoS One, 7, Article ID: e42543. http://dx.doi.org/10.1371/journal.pone.0042543
- Castoe, T.A., Poole, A.W., de Koning, A.P.J., Jones, K.L., Tomback, D.F., Oyler-McCance, S.J., Fike, J.A., Lance, S.L., Streicher, J.W., Smith, E.N. and Pollack, D.D. (2012) Rapid Microsatellite Identification from Illumina Paired-End Genomic Sequencing in Two Birds and a Snake. PLoS One, 7, Article ID: e30953. http://dx.doi.org/10.1371/journal.pone.0030953
- Chapuis, M.-P. and Estoup, A. (2007) Microsatellite Null Alleles and Estimation of Population Differentiation. Molecular Biology and Evolution, 24, 621-631. http://dx.doi.org/10.1093/molbev/msl191
- Rousset, F. (2008) Genepop’007: A Complete Re-Implementation of the Genepop Software for Windows and Linux. Molecular Ecology Resources, 8, 103-106. http://dx.doi.org/10.1111/j.1471-8286.2007.01931.x
- Benjamini, Y. and Hochberg, Y. (1995) Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B, 57, 289-300.
- Sbrocco, E.J. and Barber, P.H. (2011) Ten Polymorphic Microsatellite Loci for the Atlantic Silverside, Menidiamenidia. Conservation Genetic Resources, 3, 585-587.
- Beheregaray, L.B. and Sunnucks, P. (2000) Microsatellite Loci Isolated from Odontesthes argentinensis and the O. Perugiae Species Group and Their Use in Other South American Silverside Fish. Molecular Ecology, 9, 629-644. http://dx.doi.org/10.1046/j.1365-294x.2000.00882.x
- Barriga-Sosa, I.D.L.A., Ibañez, A.A.L. and Arredondo-Figueroa, J.L. (2002) Morphological and Genetic Variation on Seven Species of the Endangered Chirostoma “humboldtianum Species Group” (Atheriniformes: Atherinopsidae). International Journal of Tropical Biology and Conservation/Revista de Biología Tropical, 50, 199-216.
- Barriga-Sosa, I.D.L.A., Eguiarte, L.E. and Arredondo-Figueroa, J.L. (2004) Low but Significant Population Subdivision of Chirostomagrandocule from Lake Patzcuaro, Michoacan, Mexico. Biotropica, 36, 85-98.
NOTES
*Corresponding author.