Isolation and characterization of 24 polymorphic microsatellite loci for the study of genetic population structure of the sheepshead Archosargus probatocephalus (Actinopterygii, Perciformes, Sparidae)
© Seyoum et al. 2016
Received: 29 January 2016
Accepted: 21 April 2016
Published: 29 April 2016
The sheepshead (Archosargus probatocephalus) is found in nearshore waters from Nova Scotia, Canada, to Rio Grande do Sul, Brazil. In the southeastern United States two subspecies are recognized based on a number of meristic characters, primarily counts of melanistic pigment bars. The only previous study based on mtDNA control-region sequence found limited divergence between those subspecies and isolation by distance among 15 locations from Florida (Atlantic Ocean) to Texas (Gulf of Mexico). In the same study, using six sparid microsatellite markers, Bayesian analysis showed that the Gulf and Atlantic sheepshead form a single population. To reinvestigate the fine-scale genetic population structure and examine genetic support for the morphologically classified subspecies, a set of species-specific microsatellite markers was needed.
Here we report on 24 polymorphic microsatellite markers isolated from sheepshead and screened in 57 specimens from the Indian River, Florida. The average number of alleles per locus was 13.1; mean observed and expected heterozygosities were 0.68 and 0.73, respectively. Nine sparid markers screened for the same specimens showed an average of 8.6 alleles per locus; mean observed and expected heterozygosities were 0.46 and 0.55, respectively.
The polymorphic markers reported here can be used to search for genetic evidence for the morphologically defined subspecies, to elucidate the fine-scale genetic population structure of this broadly distributed coastal species, and to provide an opportunity to directly compare results of population delineation between nonspecific and species-specific markers.
The sheepshead (Archosargus probatocephalus) is an economically important, estuarine–marine teleost fish that is widely distributed from Nova Scotia, Canada, to Rio Grande do Sul, Brazil [1–3]. They migrate offshore to spawn in late winter and return to estuaries in early spring . Limited movement of sheepshead along the coast could lead to discrete populations among spawning groups, but assuming pelagic eggs disperse freely, this may facilitate gene flow and thwart the formation of population structure .
Geographic variation in sheepshead bar counts and growth rates is evident [2, 5]. Variation in melanistic bar patterns initially led to the designation of two subspecies in North America , one in the western and northern Gulf of Mexico and another in the eastern Gulf of Mexico and along the Eastern Seaboard. A recent statistical re-evaluation of morphometric data confirmed that these putative subspecies exhibited significantly different numbers not only of melanistic bars but also of several other meristic characters (i.e., scales, gill rakers, and fin rays [2, 7]). This new analyses also indicated the presence of a hybrid zone in the northeastern Gulf .
The only study of the genetic population structure of the sheepshead based on six microsatellite loci developed for other sparid showed that the Gulf and Atlantic populations belong to a single panmictic population . In the same study based on mtDNA control region sequence, no specific genetic boundaries were evident that corresponded with the two morphologically defined subspecies mentioned above. In samples from northeast Florida to Texas, mtDNA and microsatellite differentiation was attributed to isolation by distance rather than independent genetic stocks . The mtDNA, however, as a single locus may not reflect enough genealogical histories among populations and because of this has a low-resolution power and often fails to reveal fine-scale population structure. Nonspecific markers may also fall in the category of low resolution molecular techniques, particularly if few are used. Species-specific microsatellite markers, however, may represent genealogical record from the source organism with which to observe population structure in fishes that reside in open, coastal habitats. We postulate that in direct competition species-specific markers should be superior to nonspecific markers and can be helpful in re-examining genetic evidence for the validity of the subspecies and the genetic stock structure of the sheepshead.
Microsatellite loci were isolated following the PIMA (PCR-based isolation of microsatellite arrays) method of Lunt et al. , modified by Seyoum et al. . Nuclear DNA (nDNA) was first purified from liver tissue from a single sheepshead via density-gradient ultracentrifugation  to minimize competition with mitochondrial DNA during random amplified polymorphic DNA (RAPD) PCRs. RAPD PCRs were conducted in 50-µl reactions containing 15–25 ng of the purified DNA, 50 μM of dNTP mix, 0.25 μl of 0.1–mg/ml BSA; two or three primers randomly chosen from a set of 120, 10-mer RAPD primers (Qiagen Operon Inc.); 5 μl of Taq polymerase buffer (10×), 2.5 mM MgCl2 (Promega, final concentration); and 1.25 units of Go Taq DNA polymerase (Promega Corporation). The reaction profile was 94 °C for 2 min, 30 × (94 °C for 40 s, 35 °C for 40 s, 72 °C for 45 s), and final extension at 72 °C for 30 min. The purified PCR products (Agilent Technologies) were T-A cloned  into plasmid vectors Bluescript PBC KS-Agilent Technologies) that had been tailed with homemade dTTP . About 50 recombinant colonies from each of the PCR products were screened by performing PCR (12.5 μl total reaction volume) containing T3 and T7 vector primers and four repeat-specific primers (5′-[AC]10-3′, 5′-[AG]10-3′, 5′-[AGC]5-3′, 5′-[ACT]12-3′). Here, the reaction profile was 94 °C for 2 min, 35 × (94 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s), and final extension 72 °C for 7 min. The PCR products of the clones were run through a 1.5 % low-EEO agarose gel. Colonies that showed two or more bands were further amplified using only the vector primers, the products gel-purified and then cycle-sequenced from both directions using BigDye (version 3.1; Applied Biosystems), and the sequencing products visualized on an Applied Biosystems Prism™ 3130-Avant Genetic Analyzer. Primers were designed for candidate loci using OligoPerfect (Thermo Fisher Scientific); annealing temperature was adjusted to 58–60 °C and fragment size to three categories, 95–115, 125–165, and 185–250, to facilitate multiplex PCR and minimize overlapping of fragment sizes during visualization. Multiplex PCR amplifications for each specimen were carried out in an Eppendorf thermal cycler containing 50–100 ng of total DNA and three optimally selected primers, each forward primer labeled with a unique fluorescent dye. The multiplex PCR reaction consisted of a step-down profile and was as follows: 94 °C for 2 min, 5 × (94 °C for 45 s, 61 °C for 45 s, 72 °C for 45 s); 8 × (94 °C for 40 s, 59 °C for 40 s, 72 °C for 40 s); 10 × (94 °C for 35 s, 57 °C for 35 s, 72 °C for 35 s) 12 × (94 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s) and a final extension at 72 °C for 15 min. The fragments were visualized on an ABI 3130 XL genetic analyzer and genotyped using GeneMapper (version 4.0, Applied Biosystems Inc.). For fragment assays, we used a Gene Scan-500 ROX-labeled size standard.
Characteristics of 24 polymorphic microsatellite loci in 57 specimens of sheepshead (Archosargus probatocephalus) from the Indian River, Florida
Primer sequence (5′→3′) forward/reverse
Allele size range
GenBank accession no.
Microsatellite-marker GENEPOP data was generated using the Microsatellite Marker Toolkit Excel add-on. Genotypic disequilibrium among loci was examined using GENEPOP, version 3.4 . Hardy–Weinberg equilibrium (HWE) expectation, observed and unbiased expected heterozygosity estimates with Bonferroni correction, number of alleles, and polymorphic information content, a value that is indicative of a measure of the informativeness of a genetic marker for linkage studies were estimated using the program CERVUS, version 3.0.7; .
There were significant departures from HWE after Bonferroni correction at two loci and no linkage disequilibrium at any pair of loci. Analyses using the program CERVUS suggested that the observed nonconformance to HWE may have resulted from the presence of null alleles at those loci. The average number of alleles per locus was 13.1 (range 2–27); the mean observed and expected heterozygosities were, respectively, 0.68 (range 0.09–0.97) and 0.73 (range 0.29–0.95). The relative informativeness of each marker ranged from 0.24 to 0.93 (mean 0.70), with 11 of the 24 loci between 0.82 and 0.94.
Characteristics of 9 adopted sparid microsatellite DNA loci in 53 specimens of sheepshead (Archosargus probatocephalus) from the Indian River, Florida
Primer sequence (5′→3′) forward/reverse
Allele size range
GenBank accession no.
Costs of developing species-specific markers remain a concern. Compared with the enrichment protocol of developing microsatellite markers, the PIMA method requires less expertise, less time, and, so, less expense. The method’s drawback is that it results mostly in dinucleotides, which are more plentiful because their mutation rate is at least six times that of other short tandem repeats (STR). Although tri- and tetranucleotides may not be more variable than dinucleotide, they’re definitely easier to score and lead to a lower genotyping error. The various kinds of tri- tetra- penta- and hexanucleotides can be obtained with the PIMA method, but with specific designs and much more laborious search.
Nonspecific markers are used to circumvent the expertise and the expense required to develop specific markers, but they may not have adequate resolution power to reveal fine-scale population structure. In a direct comparative study, the use of 11 highly polymorphic red drum (nonspecific) markers failed to delineate the eastern and western Gulf spotted seatrout (Cynoscion nebulosus) samples as belonging to different clusters compared to only three spotted seatrout markers that accomplished the task (Seyoum et al. in preparation). While nonspecific microsatellite loci are widely used [15, 16] and could reveal strong genetic breaks between populations, there has been no study to directly compare specific vs. nonspecific microsatellite loci on equal terms of allelic variability in revealing fine genetic breaks between populations. The inadequate (and possibly misleading) results of the previous work on the study of the genetic population structure of the sheepshead that used non-specific markers may be offset by these new, species-specific markers developed in this study. We believe that these suite of species-specific markers will shed light on the stock structure of the sheepshead and on the validity of the morphologically classified subspecies. The stock structure of this fish ought to be determined now, especially because as attention increases toward regulation of other Gulf species, the sheepshead is now the target of greater demand and pressure of in the Gulf of Mexico . Updated knowledge of sheepshead stock structure will help direct relevant management actions to protect this species at the appropriate spatial scale across its range.
SS and CP did the study and carried out the computational work for locus identification and the statistical analysis. All authors read and approved the final manuscript.
We are thankful to the Tequesta Field Laboratory of the Fish and Wildlife Research Institute for collecting fish. This work was supported by the US Department of the Interior, US Fish and Wildlife Service, under the Federal Aid in Sport Fish Restoration Program, Grant F–69. The statements, findings, and conclusions are those of the authors and do not necessarily reflect the views or policies of the Department of the Interior. Mention of trade names or commercial products does not constitute their endorsement by the US Government.
Sample collection and availability of raw data
Permission to sample the fish was given by the Fish and Wildlife Research Institute and actual collection was carried out by the Fisheries Independent Monitoring program of the Institute. Actual DNA sequences of the microsatellite markers were all deposited in GenBank and the raw data is available upon request from the contact author.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Robins C, Richard C, Ray GC, Douglas J, Freud R. A field guide to Atlantic coast fishes of North America. Peterson field guide series no 32. Boston: Houghton Mifflin Co.; 1986.Google Scholar
- Caldwell DK. Systematics and variation in the sparid fish Archosargus probatocephalus. Bull South Calif Acad Sci. 1965;64:89–100.Google Scholar
- Tremain DM, Harnden CW, Adams DH. Multidirectional movements of sportfish species between an estuarine no-take zone and surrounding waters of the Indian River Lagoon, Florida. Fish Bull. 2004;102:533–44.Google Scholar
- McBride R. Managing a marine stock portfolio: stock identification, structure, and management of 25 fishery species along the Atlantic coast of the United States. N Am J Fish Manag. 2014;34:710–34.View ArticleGoogle Scholar
- Dutka-Gianelli J, Murie DJ. Age and growth of sheepshead, Archosargus probatocephalus (Pisces: Sparidae), from the northwest coast of Florida. Bull Mar Sci. 2001;68:69–83.Google Scholar
- Ginsburg I. Eight new fishes from the Gulf Coast of the United States, with two new genera and notes on geographic distribution. J Wash Acad Sci. 1952;42:84–101.Google Scholar
- Anderson JD, Karel JW, Anderson KA, Roper-Foo PA. Genetic assessment of sheepshead stock structure in the northern Gulf of Mexico: morphological divergence in the face of gene flow. N Am J Fish Manag. 2008;28:592–606.View ArticleGoogle Scholar
- Lunt D, Hutchinson W, Carvalho G. An efficient method of PCR-based isolation of microsatellite arrays (PIMA). Mol Ecol. 1999;8:891–4.View ArticleGoogle Scholar
- Seyoum S, Tringali MD, Sullivan JG. Isolation and characterization of 27 polymorphic microsatellite loci for the common snook, Centropomus undecimalis. Mol Ecol Notes. 2005;5:192–4.View ArticleGoogle Scholar
- Lansman RA, Shade RO, Shapira JF, Avise JC. The use of restriction endonucleases to measure mitochondrial DNA sequence relatedness in natural populations. III. Techniques and potential applications. J Mol Evol. 1981;17:214–26.View ArticlePubMedGoogle Scholar
- Zhou MY, Gomez-Sanchez CE. Universal TA cloning. Curr Issu Mol Biol. 2000;2:1–7.Google Scholar
- Marchuk D, Drumm M, Saulino A, Collins FS. Construction of T–vectors, a rapid and general system for direct cloning of unmodified PCR products. Nucl Acid Res. 1990;19(5):1154.View ArticleGoogle Scholar
- Rousset F. Genepop’007: a complete reimplementation of the GenePop software for Windows and Linux. Mol Ecol Res. 2008;8:103–6.View ArticleGoogle Scholar
- Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16:109–1106. doi:10.1111/j.1365-294X.2007.03089.x.View ArticleGoogle Scholar
- Barbará T, Palma-Silva C, Paggi GM, Bered F, Fay MF, Lexer C. Cross-species transfer of nuclear microsatellite markers: potential and limitations. Mol Ecol. 2007;16:3759–67.View ArticlePubMedGoogle Scholar
- Maduna SN, Rossouw C, Roodt-Wilding R, Bester-van der Merwe AE. Microsatellite cross-species amplification and utility in southern African elasmobranchs: a valuable resource for fisheries management and conservation. BMC Res Notes. 2014;7:1–12.View ArticleGoogle Scholar
- VanderKooy SJ, editor. The sheepshead fishery of the Gulf of Mexico, United States: a fisheries profile, United State. Ocean Springs: Gulf States Marine Fisheries Commission, Publication Number 143; 2006.Google Scholar
- Chen SL, Liu YG, Xu MY, Li J. Isolation and characterization of polymorphic microsatellite loci from an EST-library of red sea bream (Chrysophrys major) and cross-species amplification. Mol Ecol Notes. 2005;5:215–7.View ArticleGoogle Scholar
- Yap ES, Spencer PB, Chaplin JA, Potter IC. The estuarine teleost, Acanthopagrus butcheri (Sparidae), shows low levels of polymorphism at five microsatellite loci. Mol Ecol. 2000;9:2225–6.View ArticlePubMedGoogle Scholar
- Roques S, Galarza JA, Macpherson E, Turner GF, Rico C. Isolation and characterization of nine polymorphic microsatellite markers in the two-banded sea bream (Diplodus vulgaris) and cross-species amplification in the white sea bream (Diplodus sargus) and the saddled bream (Oblada melanura). Mol Ecol Notes. 2007;7:661–3.View ArticleGoogle Scholar
- Adcock GJ, Ramirez JHB, Hauser L, Smith P, Carvalho GR. Screening of DNA polymorphisms in samples of archived scales from New Zealand snapper. Fish Biol. 2000;56:1283–7.View ArticleGoogle Scholar
- Takagi M, Taniguchi N, Cook D, Doyle RW. Isolation and characterization of microsatellite loci from red sea bream Pagrus major and detection in closely related species. Fish Sci. 1997;1997(63):199–204.Google Scholar