Skip to main content

Development and evaluation of 16 new microsatellite loci for the rock ptarmigan (Lagopus muta) and cross-species amplification for the willow grouse (L. lagopus)



The genetic markers designed for this study can facilitate future genetic studies on the rock ptarmigan (Lagopus muta). To our knowledge no microsatellite markers have ever been developed specifically for this species before. These new microsatellite markers will be useful for population genetics studies and for future conservation projects.


Using Next Generation Sequencing 6252 potential microsatellite sequences were found. Sixteen nonpalindromic tetranucleotide microsatellites and their respective primers were selected. The markers were tested on both the rock ptarmigan and the willow grouse (L. lagopus). The number of alleles varied between 2 and 18 for the rock ptarmigan, and between 3 and 13 for the willow grouse. Expected heterozygosity was in the range 0.1244–0.8692 and 0.1358–0.8722 for the rock ptarmigan and the willow grouse, respectively.


The rock ptarmigan (Lagopus muta) and the willow grouse (L. lagopus) are two closely related bird species that inhabit the alpine zone throughout the northern hemisphere. In Scandinavia they are important game species. Both have experienced a drastic decline in recent years [1] and there is an urgent need for tools to investigate the underlying causes of the decline, including markers for population genetic studies. To our knowledge only two studies on population genetics have been carried out on the rock ptarmigan, but none of the markers used were species-specific [2, 3].

Microsatellite markers have many applications [4], and with only a limited number of markers needed, they are among the most popular type of genetic markers used in ecological studies [5]. Although single-nucleotide polymorphisms (SNPs) are becoming increasingly popular for population genetic studies, due to the ease with which they can be read, microsatellites still have some advantages [5]. Microsatellite markers are more polymorphic than SNPs [6], so a smaller number of loci is needed for the same level of precision [7], decreasing the cost of the analysis. Moreover some studies suggest that microsatellite markers are better than SNPs to detect recent population structure events [8, 9]. However, in order to prevent high scoring error rates, only long motifs (tetranucleotides), which are less subject to slippage, should be selected [10]. This also aids inter-laboratory comparison.

To avoid the problem of null alleles, polymerase chain reaction (PCR) based markers require primer sequences that bind to highly conserved regions with few or no inter-individual differences. To minimize errors due to primer site mutations it is best to use species-specific markers [11].

The number of markers required is inversely correlated with the degree of genetic differentiation across populations. Thus, for large populations or species with a high migration rate (as is the case for the rock ptarmigan), a larger number of microsatellites is required [4]. In this study, we have designed microsatellite markers specifically for the rock ptarmigan and also tested their potential with the closely-related willow grouse.

Main text

Materials and methods

For the library construction approximatively 1 mg of ethanol-preserved liver sample from a rock ptarmigan originating from the western part of the Hardangervidda plateau, Norway, was used. This was sent to Cornell University (Ithaca, New York, USA), where the following analyses were conducted as a commercial service. At the Evolutionary Genetics Core Facility a library was prepared; this was then sent to the Sequencing and Genotyping Facility (Cornell Life Sciences Core Laboratory Center) for Titanium 454 sequencing. The sequences obtained were assembled under stringent condition with a match size of 120 base pair (bp), a match percentage of 94% and a minimum sequence length of 150 bp. From this dataset all di-, tri-, tetra-, penta- and hexanucleotides repeats were screened with the software msatcommander [12, 13] and sequences less than 100 bp and/or with a small number of repeats were excluded. Primers were designed with Primerselect software (part of the Lasergene package [14]). Every microsatellite sequence chosen was run through the online version of the Basic Local Alignment Search Tool for nucleotide (BLASTn) [15] in order to check for similarity with already published sequences.

To examine whether the selected microsatellites exhibited enough polymorphism, additional rock ptarmigan samples were collected from two distinct areas 324 km apart (Table 1). DNA from tested populations came from feathers (from hunted birds), and from fecal pellets found on the ground [16].

Table 1 Summary of the genetic samples collected in southern Norway and in Sweden during 2015/16

DNA extraction, PCR and genotyping of the tested populations was done at the genetic laboratory of the University College of Southeast Norway. Two different DNA extraction kits were used depending on the type of sample; the “DNeasy Blood & Tissue kit” (Qiagen, Cat. No. 69506) for feathers, and the “QIAamp DNA Stool Mini Kit” (Qiagen, Cat. No. 51504) for fecal pellets. In order to increase the quantity of DNA extracted during the “DNeasy Blood & Tissue” procedure, 5 µL of DTT 1 M was added to the ATL buffer during the first step of the protocol. In the “QIAamp DNA Stool Mini Kit” protocol, the quantity of buffer ASL was increased from 1.4 to 1.6 mL. In order to assess the efficiency of DNA extraction, the concentration of DNA was measured with a Qubit™ dsDNA HS Assay kit (Thermo Scientific™ Cat. No. Q32854) on a Qubit 3.0 Fluorometer (Thermo Scientific™) for 30 samples. Each sample was amplified with species-specific mitochondrial DNA primers Lagsp3F, Lag3R and Mut3R [see 17] in order to distinguish L. muta and L. lagopus, whose fecal pellets are not distinguishable morphologically, and whose habitats overlap. PCR conditions follow Nyström et al. 2006 [17], and the results were read via gel electrophoresis [16]. Rock ptarmigans were identified by an amplicon size of 212 bp, and willow grouse by an amplicon size of 154 bp.

For microsatellite analysis, PCR amplification was performed with forward primers labeled with FAM NED, PET or VIC fluorescent dyes. The reaction was made with 6µL of “Qiagen Taq PCR mastermix” (Qiagen, Cat. No. 201445) with a final concentration of 1×, 0.25 µL of a 20 µM stock of each primers, 2 µL of DNA template and ultrapure water to a total volume of 12 µL. PCR was done in an Eppendorf Mastercycler® gradient thermal cycler: 15 min at 95 °C, 40 cycles consisting of 30 s at 94 °C, 90 s at 60 °C, 60 s at 72 °C, and finally 10 min at 72 °C. Microsatellite markers were multiplexed in five different PCR reactions. For every reaction a negative control was included and analysed in parallel with the samples. If any positive signal was found in the negative control the results were rejected and the entire run was repeated. In order to check for consistency of the data, 28 rock ptarmigan samples were tested twice.

For the genotyping 1.5 µL of PCR product was added in a mix containing 9.7 µL of formamide (Thermo Scientific™, Cat. No. 17899) and 0.3 µL of GeneScan 500LIZ dye Size Standard (Applied Biosystems™, Cat. No. 4322682). Genotyping was realized on a 3130xl Genetic analyzer (Applied Biosystems). Allele scoring was done using GeneMapper software V5.0 and visually controlled (Applied Biosystems). To avoid analyzing the same individual twice, an individual identity analysis was carried out with the software Cervus V3.0.7 [18]. The minimum number of matching loci was set to 7, and 1 fuzzy matching was allowed. Microsatellite markers were then tested for the presence of null alleles with the software microchecker [19]. Linkage disequilibrium and Hardy–Weinberg equilibrium were tested with the software GENEPOP [20]. In order to control for multiple testing, the p value threshold was adjusted with Bonferroni correction. Number of alleles Na, observed heterozygosity Ho and expected heterozygosity He were calculated with the software Genetix V4.05 [21, see Table 2]. The same procedures were also performed on samples from the willow grouse, which were collected from two different areas 124 km apart in Norway and in Sweden (Table 1).


After the Titanium 454 sequencing 31037 sequences were obtained, which resulted in 6252 potential microsatellite sequences after screening with msatcommander. Sixteen microsatellite nonpalindromic tetranucleotide repeats with lengths between 85 and 226 bp, and with number of repeats between 5 and 15 were selected. BLASTn analysis of the sequences gave two positive matches with red jungelfowl (Gallus gallus) sequences. The first sequence named CH261-100D18 (Accession Number AC147597) had an 87% match with Mut18, but the repeat motif was not present on the Gallus galllus sequence. The second sequence named 8F6 (Accession Number X78623.1 [22]), had a 90% match with Mut01 and the repeat motif was identical. Because we couldn’t find any publications using this microsatellite on grouse, we kept it for analysis.

In total 63 rock ptarmigan individuals were tested (Table 1). Measured DNA quantity for the 30 samples tested ranged from 0.06 to 44 ng/µL (average: 7 ng/µL). Identity analysis with Cervus detected three redundant individuals that were removed prior to analysis. Replicate data was available for 28 rock ptarmigan samples, representing a total of 371 pairs of allele identification (failed reaction were not counted). There was in total nine mismatches between replicates, one mismatch for Mut18, Mut4 and Mut22, two mismatches for Mut12 and three mismatches for Mut20. The software microchecker detected one microsatellite marker (Mut12) that presented a high frequency of null alleles. No significant linkage disequilibrium was found after Bonferroni-correction. A significant deviation (p < 0.0031) from Hardy–Weinberg (heterozygosity deficit) was found after Bonferroni correction on markers Mut01 and Mut12. The total number of alleles per microsatellite ranged from 2 (Mut24) to 18 (Mut17). Expected heterozygosity was between 0.1244 (Mut24) and 0.8332 (Mut08), and observed heterozygosity ranged between 0.1333 (Mut24) and 0.8571 (Mut04) (Table 2).

Table 2 Characterization of 16 polymorphic microsatellites for the rock ptarmigan (regular font) and the willow grouse (italics)

For the willow grouse the microsatellite markers Mut03 and Mut24 were hard to read and therefore removed from analysis, but the other microsatellites produced amplicon of size and pattern equivalent to the rock ptarmigan. Two microsatellite markers, Mut06 and Mut09, presented a high frequency of null alleles. No significant linkage disequilibrium was found after Bonferroni-correction. A significant deviation (p < 0.0033) from Hardy–Weinberg (heterozygosity deficit) was found after Bonferroni correction for markers Mut02, Mut06 and Mut09. The total number of alleles per microsatellite for the willow grouse ranged from 3 (Mut18) to 13 (Mut20), expected heterozygosity from 0.1358 (Mut14) to 0.8722 (Mut08), and observed heterozygosity from 0.0952 (Mut14) to 0.8421 (Mut04).


Most of the tested microsatellite markers produced good results, except for Mut01 and Mut12, which both presented a heterozygosity deficit. The heterozygosity deficit of Mut01 could be explained by the Wahlund-effect [23], as we pooled two different populations together for the Hardy–Weinberg equilibrium calculation. When calculating heterozygosity deficit separately for each of the two populations, only Mut12 presented a significant heterozygosity deficit after Bonferroni correction. The significant heterozygosity deficit for Mut12 is probably due to the presence of null alleles.

The need for high quality and species-specific microsatellite markers resulted in the development of a new library for the rock ptarmigan. The microsatellites have primers optimised for the target species and are designed to have a low amount of stutter, by selecting tetranucleotide repeat motifs, which are significantly less prone to slippage and a repeat number of five in order to achieve high degrees of polymorphism and heterozygosity, [24]. The high level of heterozygosity and the important number of alleles makes them suitable for studies on population structure and/or dispersal behaviour.


  • The number of published microsatellite markers could be too low for population structure analysis at a local scale.

  • While tetranucleotide motif are less prone to slippage, they are also less polymorphic than di- or trinucleotide motif.

  • These microsatellite markers have only been tested on two species of grouse.



polymerase chain reaction


base pair


Basic Local Alignment Search Tool for nucleotide


number of alleles

Ho :

observed heterozygosity

He :

expected heterozygosity


  1. Henriksen S, Hilmo O. Norsk rødliste for arter 2015. Norge: Artsdatabanken; 2015. p. 193.

    Google Scholar 

  2. Caizergues A, Bernard-Laurent A, Brenot J-F, Ellison L, Rasplus JY. Population genetic structure of rock ptarmigan Lagopus mutus in Northern and Western Europe. Mol Ecol. 2003;12:2267–74.

    Article  CAS  PubMed  Google Scholar 

  3. Bech N, Boissier J, Drovetski S, Novoa C. Population genetic structure of rock ptarmigan in the ‘sky islands’ of French Pyrenees: implications for conservation. Anim Conserv. 2009;12:138–46.

    Article  Google Scholar 

  4. Selkoe KA, Toonen RJ. Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers. Ecol Lett. 2006;9:615–29.

    Article  PubMed  Google Scholar 

  5. Guichoux E, Lagache L, Wagner S, Chaumeil P, Léger P, Lepais O, et al. Current trends in microsatellite genotyping. Mol Ecol Resour. 2011;11:591–611.

    Article  CAS  PubMed  Google Scholar 

  6. Fabbri E, Caniglia R, Mucci N, Thomsen HP, Krag K, Pertoldi C, et al. Comparison of single nucleotide polymorphisms and microsatellites in non-invasive genetic monitoring of a wolf population. Arch of Biol Sci. 2012;64:321–35.

    Article  Google Scholar 

  7. Gärke C, Ytournel F, Bed’hom B, Gut I, Lathrop M, Weigend S, et al. Comparison of SNPs and microsatellites for assessing the genetic structure of chicken populations. Anim Genet. 2012;43:419–28.

    Article  PubMed  Google Scholar 

  8. DeFaveri J, Viitaniemi H, Leder E, Merila J. Characterizing genic and nongenic molecular markers: comparison of microsatellites and SNPs. Mol Ecol Resour. 2013;13:377–92.

    Article  CAS  PubMed  Google Scholar 

  9. Granevitze Z, David L, Twito T, Weigend S, Feldman M, Hillel J. Phylogenetic resolution power of microsatellites and various single-nucleotide polymorphism types assessed in 10 divergent chicken populations. Anim Genet. 2014;45:87–95.

    Article  CAS  PubMed  Google Scholar 

  10. Flores-Renteria L, Whipple AV. A new approach to improve the scoring of mononucleotide microsatellite loci. Am J Bot. 2011;98:e51–3.

    Article  PubMed  Google Scholar 

  11. Glenn TC, Schable NA. Isolating microsatelite DNA Loci. Methods Enzymol. 2005;395:202–22.

    Article  CAS  PubMed  Google Scholar 

  12. Faircloth BC. Msatcommander: detection of microsatellite repeat arrays and automated, locus-specific primer design. Mol Ecol Resour. 2008;8:92–4.

    Article  CAS  PubMed  Google Scholar 

  13. Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. In: Misener S, Krawetz SA, editors. Methods in molecular biology. Totowa, NJ: Humana Press; 2000. p. 365–386.

    Google Scholar 

  14. DNASTAR. Lasergene package. Accessed 25 Sept 2017.

  15. NCBI. Basic Local Alignment Search Tool. Accessed 20 Sept 2017.

  16. Bergan F, Sæbø M, Parker H. Evaluation of in-field stability of mitochondrial and nuclear DNA in snow sampled fecal pellets from Rock ptarmigan (Lagopus muta). Oecol Mont. 2016;25:1–14.

    Google Scholar 

  17. Nyström J, Dalén L, Hellström P, Ekenstedt J, Angleby H, Angerbjörn A. Effect of local prey availability on gyrfalcon diet: DNA analysis on ptarmigan remains at nest sites. J Zool. 2006;269:57–64.

    Article  Google Scholar 

  18. Kalinowski ST, Taper ML, Marshall TC. Revising how the computer program cervus accommodates genotyping error increases success in paternity assignment. Mol Ecol. 2007;16:1099–106.

    Article  PubMed  Google Scholar 

  19. Van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.

    Article  Google Scholar 

  20. Raymond M, Rousset F. GENEPOP (Version 1.2): population genetics software for exact tests and ecumenicism. J Hered. 1995;86:248–9.

    Article  Google Scholar 

  21. Belkhir K, Borsa P, Chikhi L, Raufaste F, Bonhomme N. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Montpellier: Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II; 2004.

    Google Scholar 

  22. Gibbs M, Dawson DA, McCamley C, Wardle AF, Armour JA, Burke T. Chicken microsatellite markers isolated from libraries enriched for simple tandem repeats. Anim Genet. 1997;28:401–17.

    Article  CAS  PubMed  Google Scholar 

  23. Allendorf FW, Luikart G. Conservation and the Genetics of Populations. Hoboken: Blackwell Publishing; 2007. p. 664.

    Google Scholar 

  24. Edwards A, Civitello A, Hammond HA, Caskey CT. DNA typing and genetic mapping with trimeric and tetrameric tandem repeats. Am J Hum Genet. 1991;49:746–56.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Author’s contributions

ØS originated with the project idea, supervised the project and secured the funding. MS developed the study design and the methods. FB, collected and prepared the sample for the library construction and sequencing. FB and JMC did the sample collection. FB selected the microsatellite loci and developed the primers. AJ participated in the study design, selection of the microsatellite loci and the primers developments. FB and JMC did the laboratory work for the genotyping. JMC did the data analysis and the writing of the manuscript. All authors, participated in the draft. All authors read and approved the final manuscript.


We would like to thank Pauline Casters, Sylvain Clot, Richard Hančinský, Inger Marie Hernes, Zuzana Krajova, Iver Lodden Boine, Damien Picard, Tereza Prikrylova, Christian Robstad, Johanna Skrutvold, Juraj Sovcik and Michal Torma for assistance in the field, and Martin Hagen Ring for helping out in the field and in the laboratory. Finally, we would like to thank all the hunters that sent us feather samples.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The 6252 potential microsatellite sequences that were screened, their respective primers, and the 31037 raw sequences obtained from the Titanium 454 sequencing are available on demand to the corresponding author. Genotype data for every individual and microsatellite loci are also available on demand to the corresponding author.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


The study was funded by the University College of Southeast Norway.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Jean-Marc Costanzi.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Costanzi, JM., Bergan, F., Sæbø, M. et al. Development and evaluation of 16 new microsatellite loci for the rock ptarmigan (Lagopus muta) and cross-species amplification for the willow grouse (L. lagopus). BMC Res Notes 11, 147 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Birds
  • Lagopus muta
  • Lagopus lagopus
  • Microsatellites
  • Population genetics