Skip to main content

Development and characterization of microsatellite markers in Gaultheria pumila Lf. (Ericaceae)



Polymorphic microsatellite markers were developed for Gaultheria pumila (Ericaceae) to evaluate genetic diversity and population structure within its native range in Chile. This is a very important Ericaceae endemic to Chile with a large commercial potential. Its resistance to different abiotic conditions makes it a valuable target for genetic improvement.


Ten polymorphic simple sequence repeat (SSR) loci were isolated from Gaultheria pumila using new-generation 454 FLX Titanium pyrosequencing technology. The mean number of alleles per locus ranged from 2 to 4. Observed and expected heterozygosity ranged from 0.00 to 1.0 and 0.00 to 0.64, respectively.


From 10 SSR markers developed for G. pumila, 9 markers are promising candidates for analyzing genetic variation within or between natural populations of G. pumila and other species from the same genus.


The Gaultheria genus belongs to the Ericaceae family and is adapted to tropical and temperate conditions. In Chile, all of these species are commonly known as “chaura” or “mutilla del zorro” [1]. Gaultheria pumila is one of the most important species of this genus in Chile, regarding population size. In the wild, it is considered as a primary colonizing plant [2], capable of associating with symbiotic mycorrhizas [3]. The species shows a high degree of phenotypic diversity and even a botanic variety has been identified (Gaultheria pumila var. leucocarpa). It is an interesting species for domestication since it has shown a high level of phenotypic variability and high levels of polyphenols and anthocyanins [4, 5]. However, the lack of basic knowledge of agronomically important traits as well its genetic variability are serious limitations to its domestication process for its agronomic and commercial use. For this reason, the aim of this study was to develop a set of polymorphic microsatellite markers to be used as a tool to assess the diversity and genetic structure of G. pumila and other species of the Gaultheria genus.


In order to carry out complete sequencing, four samples of G. pumila were collected in two places. The Villarrica volcano in Región de la Araucanía, Chile, S: 39°21′570, W: 71′57′865 (El Playón sector), and near the Puyehue volcano, in Región de Los Lagos, Chile, S: 40°41′484, W: 72°32′719 (Orilla del Río sector). All the samples from the protected natural areas were collected under the authorization of the regulatory Chilean Agency, Corporación Nacional Forestal (CONAF).

The samples for sequencing were selected according to fruit color: red, pink and white. Total genomic DNA of G. pumila was extracted from young leaves with a modified CTAB (cetyltrimethylammonium bromide) method [6]. The quality and quantity of isolated DNA was determined using a Nanodrop and agarose gel electrophoresis. The library preparation and shotgun pyrosequencing of 5 μg DNA aliquot on a 454 GS-FLX instrument (Roche Applied Science, MACROGEN, Ltd. Seoul, South Korea) was prepared with the kit of Titanium Pico Titer Plate (Roche Diagnostics), following the manufacturer’s protocols.

A total of 164,000 sequences were read, with an average length of 417 bp, which were generated by three independent runs. The sequences were assembled using the software GS De Novo Assembler v2.9 (

The MSATCOMANDER ( software for locating microsatellites was used. MSATCOMMANDER is a program written to locate microsatellite repeats within fasta-formatted sequences or consensus files. MSATCOMMANDER will search for all di-, tri-, tetra-, penta, and hexa-nucleotide repeats. The search parameters were: mononucleotide repeat length: 10, dinucleotide repeat length: 6, trinucleotide repeat length: 4, tetranucleotide repeat length: 4, pentanucleotide repeat length: 4, and hexanucleotide repeat length: 4.

The primers were designed using the software PRIMER 3 [7]. The PCR fragment amplification and validation of selected SSRs were performed following the method developed by [8] which uses three primers: a forward primer with an M13 (-21) tail at its 5′ end, a normal reverse primer and the universal M13 (-21) primer labeled with either 6-FAM, VIC, PET or NED fluorochromes. The microsatellite information and GenBank accession numbers are listed in Table 1. PCR reactions were performed in a 20 μL reaction mixture with 10 ng templates DNA, 0.15 mM of each dNTP; 1× Taq polymerase reaction buffer; 1.5 mM MgCl2; 0.025 μM forward primer; 0.1 μM reverse primer; 0.1 μM M13 primer and 0.35 U Taq DNA polymerase.

Table 1 Characteristics of 10 microsatellite loci and primer pair’s development for G. pumila

PCR amplifications were performed in an Applied Biosystems Veriti (Life Technologies), under the following conditions: initial denaturation at 94 °C for 5 min; 30 cycles of 30 s (s) at 94 °C, annealing temperature specific to each primer pair for 45 s, extension at 72 °C for 45 s, followed by 8 cycles of 30 s at 94 °C, annealing at 53 °C for 45 s, extension at 72 °C for 45 s and a final extension at 72 °C for 30 min. The annealing temperature (°C) requirements of primers are specified in Table 1.

PCR products were resolved as following: 2 μL of PCR products were mixed with 10 μL HiDi formamide (Applied Biosystems) and 0.2 μL GeneScan 500LIZ size standard (Life Technologies, Foster City, CA) and separated by capillary electrophoresis on an ABI 3130xl Prism Genetic Analyzer with POP-7 polymer (Life Technologies, Foster City, CA) in the Genetic Resources Unit, La Platina-INIA, Chile. Allele sizes were automatically calculated with GeneMapper software v4.0 and manually checked.

Three populations of G. pumila were selected to evaluate the variability in the isolated loci: Región de la Araucanía, (Villarrica volcano, n = 10); Región de los Lagos, (Puyehue volcano, n = 10); Región de Magallanes y la Antártica chilena (Punta Arenas, n = 10). Genetic parameters such as observed number of alleles (Na) and observed and expected heterozygosity (HoHe) were estimated with PopGene [9]. The inbreeding coefficient (FIS) [10] was determined using Genetix [11]. Deviation from the Hardy–Weinberg equilibrium was determined with GENEPOP v 4.2 [12] (Table 2). The presence of null alleles was checked using MICRO-CHECKER version 2.2.3 [13].

Table 2 Estimated genetic parameters and Hardy–Weinberg test results in 10 microsatellites loci in three wild populations of G. pumila

In addition, cross-amplification was tested in two other species, Gaultheria mucronata and Gaultheria caespitosa. PCR reactions and electrophoresis were performed according to the conditions described above.


Primers were considered successful when one clear distinct band was detected on 2% agarose gel. Primers were designed with flanking repeat lengths of at least 10, 6 and 4 for mononucleotides, dinucleotides, and higher core repeats, respectively. Mononuclotide-containing repeats were the most frequent (59%), followed by dinucleotides (23%), trinucleotides (12%), tetranucleotides (3%) and pentanucleotides (2%).

Seventeen SSRs were tested and they included primers for four mononucleotides, five dinucleotides, four trinucleotides, three tetranucleotides, and one pentanucleotide. The first test consisted of selecting the best SSRs that showed amplification in a 2% agarose gel after electrophoresis. Ten SSRs were selected for further development and analysis. From these microsatellite loci, nine were variable and polymorphic between populations.

The microsatellite GP.7 was inconsistent in different tested condition, and it was not included further. The mean number of alleles per locus ranged from 1 to 4, with an average of 1.96. Observed and expected heterozygosity ranged from 0.00 to 1.0 and 0.00 to 0.67, respectively. In addition, three microsatellites loci GP.12, GP.13, GP.17 exhibited significant deviation from Hardy–Weinberg equilibrium, for the Puyehue population. In the Villarrica population GP.14 was monomorphic, while the loci GP.9, GP.10, GP. 12, GP.13 and GP.18 had no significant deviations. In Punta Arenas, four loci (GP.9, GP.10, GP.12 and GP.15) did not show significant deviations from HWE. Such deviations can be due to the small sample size, selfing, or to the substructure of the populations. Additionally, after visualization of the PCR products for each loci by gel electrophoresis, it was found that out of 10 SSR loci, 5 showed cross amplification in G. mucronata (GP.14, GP.15, GP.16, GP.17 and GP.18) and 6 G. caespitosa (GP.9, GP.13, GP.14, GP.15, GP.16 and GP.18) (Table 2). This suggests that the markers could be useful in carrying out studies of genetic diversity in other Gaultheria species.


Nine microsatellites primer pairs were developed and proved to be polymorphic in three populations of G. pumila. These SSRs will be useful for future studies of population structure and assessment of diversity in this species. Also, these findings provide a basis for starting a domestication program for this species, as well as other related species in the Ericaceae family. Further research will be needed to find new SSR microsatellite markers and test them among G. pumila and other species from this genus.



Corporación Nacional Forestal


cetyltrimethylammonium bromide


simple sequence repeat


polymerase chain reaction


  1. 1.

    Teillier S, Escobar F. Revisión del género Gaultheria L. (Ericaceae) en Chile. Gayana Botánica. 2013;70(1):136–53.

    Article  Google Scholar 

  2. 2.

    Luteyn JL. Diversity, adaptation, and endemism in neotropical Ericaceae: biogeographical patterns in the Vaccinieae. Bot Rev. 2002;68(1):55–87.

    Article  Google Scholar 

  3. 3.

    Medina J, Davinson F, Carrillo R, Rodriguez M. Simbiosis micorrícica presente en Nothofagus dombeyi (Mirb.) Oerst y Gaultheria pumila (L. fil) Middleton. Especies Colonizadoras de Escoriales Volcanicos en la Region de la Araucanía, Chile. In: I Congr. Nac. Conserv. la Biodivers. a través su uso, 23 al 25 Noviembre 2009, Sucre Boliv; 2009. p. 1–10.

  4. 4.

    Middleton DJ. A chemotaxonomic survey of flavonoids and simple phenols in the leaves of Gaultheria L. and related genera (Ericaceae). Bot J Linn Soc. 1992;110(4):313–24.

    Article  Google Scholar 

  5. 5.

    Villagra E, et al. Morphometric and phytochemical characterization of chaura fruits (Gaultheria pumila): a native Chilean berry with commercial potential. Biol Res. 2014;47(1):1–8.

    Article  Google Scholar 

  6. 6.

    Murray MG, Thompson WF. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 1980;8(19):4321–6.

    CAS  Article  Google Scholar 

  7. 7.

    Untergasser A, et al. Primer 3-new capabilities and interfaces. Nucleic Acids Res. 2012;40(15):1–12.

    Article  Google Scholar 

  8. 8.

    Schuelke M. An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000;18(2):233–4.

    CAS  Article  Google Scholar 

  9. 9.

    Yeh F, Rong-cai Y, Boyle T, Freeware MW. POPGENE, the user-friendly shareware for population genetic analysis. Mol Biol Biotechnol. Center, Univ. Alberta, Edmonton, Alberta, Canada, no. August; 1997. p. 1–29.

  10. 10.

    Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution (N. Y.). 1984;38(6):1358–70.

    CAS  Google Scholar 

  11. 11.

    Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F. Genetix 4.05 Logiciel sous Windows™ pour la genetique des populations. Montpellier: Université de Montpellier; 2003.

    Google Scholar 

  12. 12.

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

    Article  Google Scholar 

  13. 13.

    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(3):535–8.

    Article  Google Scholar 

Download references

Authors’ contributions

RGG led this investigation; JP and RGG, planned and analyzed the raw data as well administered the funds of this research; MS, HP, BCh and EG, performed all the sampling and field work; KQ, BC and PDSC contributed with the statistical and genetical analysis; JP, PC, BCh and RGG performed the laboratory work and sampling processes. All authors read and approved the final manuscript.


The authors would like to thank the Ecuadorian Government (Senescyt) and the Facultad de Agronomía of the Universidad Técnica de Manabí for the scholarship to José Pico Mendoza. The authors thank the Corporación Nacional Forestal for the support and authorization to access the protected natural areas where the study was developed.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The authors declare that all the raw data and materials can be accessed from the uploaded files accompanying this submission.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


We are grateful to Gobierno Regional del Maule for supporting this research through the project Laboratorio Regional de Biotecnologías Aplicadas (VITOTRECH II-FIC-BIC 30.136.372-0). This research was also partially supported by institutional funds of the Company Sociedad de Investigación y Servicios BioTECNOS Ltda., from Chile.

Publisher’s Note

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

Author information



Corresponding authors

Correspondence to Rolando Garcia-Gonzales or José Pico-Mendoza or Karla Quiroz.

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

Garcia-Gonzales, R., Pico-Mendoza, J., Quiroz, K. et al. Development and characterization of microsatellite markers in Gaultheria pumila Lf. (Ericaceae). Biol Res 51, 42 (2018).

Download citation


  • Microsatellite
  • Pyrosequencing
  • Chilean berry
  • Gaultheria pumila