- Short Report
Characterization of polymorphic microsatellite loci in the Antarctic krill Euphausia superba
BMC Research Notesvolume 7, Article number: 73 (2014)
The Antarctic krill, Euphausia superba is a pelagic crustacean, abundant in high-density swarms (10 000 – 30 000 ind/m2) with a circumpolar distribution and a key role in the food web of the Southern Ocean. Only three EST derived microsatellite markers have been used in previous genetic studies, hence we developed additional highly polymorphic microsatellite markers to allow robust studies of the genetic variability and population differentiation within this species.
The microsatellite markers described here were obtained through an enriched genomic library, followed by 454 pyrosequencing. A total of 10 microsatellite markers were tested in 32 individuals from the Antarctic Peninsula. One of the tested loci was fixed for one allele while the other was variable. Of the remaining nine markers, seven showed no departure from Hardy-Weinberg equilibrium. The mean number of alleles was 14.9.
These markers open perspectives for population genetic studies of this species to unravel genetic structure, dispersal and population biology, vital information for future conservation.
The Southern Ocean, comprising around 20% of the oceanic surface of the Earth, is considered to play a crucial role in the regulation of our planet’s climate  and is extremely vulnerable to global climatic change. One of Southern Ocean’s most important species is the Antarctic krill (Euphausia superba, Dana 1852), a shrimp-like crustacean that plays a central role in this ecosystem, being both a grazer of phytoplankton, bacteria and micro-zooplankton, and prey for vertebrates such as fish, seabirds, seals, penguins and whales [2, 3]. It has a circumpolar distribution and is found in high-density swarms of 10,000 to 30,000 ind/m2. Krill is an important economic resource with estimated catches of 125,000 to 150,000 tons/year (CCAMLR, Commission for the Conservation of Antarctic Marine Living Resources). There is a serious concern about the combined effects of climate change, ocean acidification and fishing pressure on krill, and consequently on the vulnerability of the communities that depend on them .
For conservation and management purposes, knowledge of the genetic diversity and population differentiation of the target species is fundamental. To date some studies have addressed population differentiation in E. superba, using allozymes, mitochondrial DNA and EST-linked microsatellite markers [5–11]. However, the availability of microsatellite markers is limited and their development has proven difficult, as reported by other authors . This could possibly be linked to the low GC content of the krill genome (32%), and its susceptibility to damage via UV-B radiation . The main findings of the genetic studies performed to date suggest genetic homogeneity at large geographical scales. However several exceptions have been found, such as evidence of small-scale genetic heterogeneity found using mitochondrial SNP’s , and weak but significant genetic differentiation between populations . Some authors suggest that temporal variability might explain the genetic differentiation found .
Here we report the development and characterization of polymorphic microsatellite loci for Euphausia superba. These markers will be useful for population genetic studies of this species, with the potential to provide fundamental information for conservation studies and management.
Whole genomic DNA was isolated using the CTAB method  from three individuals of E. superba, from muscle tissue kept at -80°C. To isolate the microsatellite sequences, a combination of an SSR-enrichment protocol with 454 pyrosequencing was performed by a commercial company (Ecogenics GmbH, Zürich, Switzerland). One CT/GT enriched library was generated, the insert size of the libraries was 500-800 bp and the length of 454 reads was 233 bp. A total of 38396 sequences were obtained, of these 1286 had microsatellite repeats. A total of 210 primer pairs were designed using primer3 core  of which 48 primer pairs were delivered by the service provider following amplification tests on three individuals. We tested all primers for polymorphism using a panel of 7 individuals from the Antarctic Peninsula. An M13-tail (TGTAAAACGACGGCCAGT) was added at the 5’ end of all forward primers to enable fluorescent-dye labelling . Of all the primer pairs tested, we were able to optimize 10 polymorphic markers (Table 1).
Amplification reactions in 10 μL contained 10 ng of genomic DNA, 1x Qiagen HotStart Taq buffer, 200 μM of dNTP’s, 0.04 μM of forward primer, 0.16 μM of reverse primer and fluorescently-labeled M13 primer, and 0.5 U of HotStart Taq polymerase (Qiagen). PCRs were conducted in a Perkin-Elmer GeneAmp7200 (Waltham, MA, USA) with the following program: 15 min at 95°C; 30 cycles composed of 30 s denaturation at 95°C, 45 s at the annealing temperature (Table 1) and 45 s elongation at 72°C, followed by an additional 8 cycles composed of 30 s of denaturation at 95°C, 45 s at 53°C and 45 s elongation at 72°C. A final 30 min elongation step at 72°C was performed. PCR products were amplified with M13 primers end-labelled with different fluorescent dyes, FAM, ATT550 or HEX to allow multiplexing. Fragments were separated on an ABI3130XL automated sequencer (Applied Biosystems, Foster City, CA, USA) with Rox350 size standard. Alleles were scored using Peak Scanner 1.0 (Applied Biosystems).
The variability of the markers was tested on 32 individuals sampled around the Antarctic Peninsula (between 69.1489°S; 73.526°W and 61.740°S; 53.7692°W). The number of alleles per locus (A), the observed (HO) and expected (HE) heterozygosities (Table 1) and the heterozygote deficiency (FIS) were calculated using the software GENETIX 4.05 . The majority of the optimized markers (9) were highly polymorphic while in locus Esup5 one allele (147) is present in all individuals while there is variation in the other allele sizes, potentially limiting its usefulness for population genetic studies. The number of alleles found for the 9 loci ranged from 5 (Esup2) to 28 (Esup9); HO ranged from 0.36 to 0.90 and HE from 0.37 to 0.96. Significant heterozygote deficiency was observed for 2 markers (Esup6 and 9). Null alleles might occur at these loci, as confirmed by further analysis using the MICRO-CHECKER software . We tested for linkage disequilibrium (LD) between all pairs of loci using the software GENETIX 4.05 . The significance of the results was tested with 10000 permutations at the 5% level, and one pair showed significant LD (Esup3-Esup6, p = 0.09).
The high variability displayed by these microsatellite loci should be useful for assessing the genetic structure of E. superba at different spatial and temporal scales. In particular these novel markers should allow testing controversial hypotheses of spatial genetic structure versus panmixia in Antarctic krill, seasonal and inter-annual population variability, and effects of fishing pressure on population genetic diversity.
Availability of supporting data
The microsatellite sequences are available through the National Centre for Biotechnology Information (see http://www.ncbi.nlm.nih.gov/). The accession numbers on the repository are the following: GenBank accession no. KF648623 through KF648632.
Sarmiento JL, Hughes TMC, Stouffer RJ, Manabe S: Simulated response of the ocean carbon cycle to anthropogenic climate warming. Nature. 1998, 393: 245-249. 10.1038/30455.
Mauchline J: The biology of mysids and euphausiids. Adv Mar Biol. 1980, 18: 1-681.
de Santana CN, Rozenfeld AF, Marquet PA, Duarte CM: Topological properties of polar food webs. Mar Ecol Prog Ser. 2013, 474: 15-26.
Kawaguchi S, Ishida A, King R, Raymond B, Waller N, Constable A, Nicol S, Wakita M, Ishimatsu A: Risk maps for Antarctic krill under projected Southern Ocean acidification. Nature Clim Change. 2013, -doi:10.1038/nclimate1937
Fevolden S, Ayala F: Enzyme polymorphism in Antarctic krill (Euphausiacea); genetic variation between populations and species. Sarsia. 1981, 66: 167-181.
Kuhl S, Schneppenheim R: Electrophoretic investigation of genetic variation in two krill species Euphausia superba and E. crystallorophias (Euphausiidae). Polar Biol. 1985, 6: 17-23.
Fevolden SE, Schneppenheim R: Genetic homogeneity of krill (Euphausia superba Dana) in the Southern Ocean. Polar Biol. 1989, 9: 533-539. 10.1007/BF00261038.
Zane L, Ostellari L, Maccatrozzo L, Bargelloni L, Battaglia B, Patarnello T: Molecular evidence for genetic sudivision of Antarctic krill populations. Proc Royal Soc London B. 1998, 265: 2387-2391. 10.1098/rspb.1998.0588.
Goodall-Copestake WP, Perez-Espona S, Clark MS, Murphy EJ, Seear PJ, Tarling GA: Swarms of diversity at the gene cox1 in the Antarctic krill. Heredity. 2010, 104: 513-518. 10.1038/hdy.2009.188.
Bortolotto E, Bucklin A, Mezzavilla M, Zane L, Patarnello T: Gone with the currents: lack of genetic differentiation at the circum-continental scale in the Antarctic krill Euphausia superba. BMC Genet. 2011, 12: 32-
Jarman S, Elliott N, Nicol S, McMinn A, Newman S: The base composition of the krill genome and its potential susceptibility to damage by UV-B. Antarct Sci. 1999, 11: 23-26.
Batta-Lona PG, Bucklin A, Wiebe PH, Copley NJ, Patarnello T: Population genetic variation of the Southern Ocean krill, Euphausia superba, in the Western Antarctic Peninsula region based on mitochondrial single nucleotide polymorphisms (SNPs). Deep-Sea Res II. 2011, 58: 1652-1661. 10.1016/j.dsr2.2010.11.017.
Doyle JJ, Doyle JL: Isolation of plant DNA from fresh tissue. Focus. 1990, 12: 13-15.
Untergrasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG: Primer3- new capabilities and interfaces. Nucleic Acids Res. 2012, 40 (15): e115-10.1093/nar/gks596.
Schuelke M: An economic method for the fluorescent labelling of PCR fragments. Nature Biotech. 2000, 18: 233-234. 10.1038/72708.
Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F: Laboratoire Génome, Populations, Interactions, CNRS UMR 5000, Université de Montpellier II. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. 1996-2004, Montpellier (France)
Van Oosterhout C, Hutchinson W, Wills DPM, Shipley P: Micro-Checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004, 2: 377-379.
This is a contribution to the International Polar Year program. This study was supported by the Portuguese Science Foundation (FCT) through project PTCD/MAR/72630/2006 (GAP) and postdoctoral fellowship SFRH/BPD/39097/2007 (ST) and by the project ATOS funded by the Spanish Ministry of Economy and Competitivity (POL2006-00550/CTM). We thank the crew of R/V Hespérides, technicians from the UTM and cruise participants for help and support.
The authors declare they have no competing interests.
GAP was responsible for the design and implementation of the study, supervision of the work and processing interpretation of the results. RC and ST participated in data analysis, microsatellite marker validation and drafted the manuscript. CMD coordinated the field sampling and participated in data interpretation. All authors read and approved the final manuscript.