Translator Disclaimer
3 November 2014 Microsatellite Loci in Two Epiphytic Lichens with Contrasting Dispersal Modes: Nephroma laevigatum and N. parile (Nephromataceae)
Rocío Belinchón, Christopher J. Ellis, Rebecca Yahr
Author Affiliations +

Dispersal among populations can explain gene flow patterns, can help in understanding the processes driving genetic structure and diversity, and can explain population persistence under different ecological and environmental scenarios (Lowe et al., 2004). Dispersal traits (i.e., reproductive strategy and propagule size) reflect a trade-off between dispersal ability and establishment success, and they have important consequences for population genetic structure (Hartl and Clark, 1997). In lichens, the contrast between sexual vs. asexual reproductive modes has provided a convenient system for evaluating consequent patterns of genetic structure and diversity (Werth, 2010); this is especially the case for lichen epiphytes that occur on trees mimicking discrete habitats (islands) of measurable size and age, and that are linked through dispersal dynamics to a network of surrounding trees, and at a larger scale, to adjacent woodland stands. Asexually reproducing lichen species disperse both their symbionts together (fungus and photosynthetic alga or cyanobacteria), but their propagules are larger compared to sexually reproducing species, which have a higher likelihood of long-distance dispersal (Seaward, 2008). However, this broad generalization in dispersal ability and its consequences for gene flow remain contentious owing to contradictory results (Werth et al., 2014). In particular, the majority of epiphyte studies have focused on a single model species—Lobaria pulmonaria (L.) Hoffm. (e.g., Scheidegger and Werth, 2009)—that can reproduce both sexually and asexually, making it therefore problematic to unequivocally partition the ecological consequences of contrasting dispersal modes.

The focus of this paper is to describe microsatellite markers developed for two codistributed epiphytic lichens with different reproductive strategies and dispersal modes: i.e., Nephroma laevigatum Ach. and N. parile (Ach.) Ach. Both species are morphologically very similar and have a Nostoc photobiont, but they have developed contrasting reproductive traits. Nephroma laevigatum typically reproduces via sexual ascospores (spores = 18.5 × 6 µm), while the predominantly asexually reproducing N. parile normally produces marginal soredia (soredia = 150 µm) and only very infrequent apothecia (Smith et al., 2009). It has been shown using multigene phytogenies that these ecologically similar foliose lichens are each monophyletic and closely related to each other, although they are not sister species (Lohtander et al., 2002; Sérusiaux et al., 2011). Both species have oceanic to boreal-montane distributions, are common in ancient woodlands in Europe, and are sensitive to SO2 pollution (Smith et al., 2009).


For primer design, fungal genomic DNA was extracted from eight individuals per lichen species collected from four different populations across a climatic gradient in Scotland (Appendix 1). DNA was isolated from 20 mg of dried thallus material using the DNeasy Plant Mini Kit (QIAGEN, Hilden, Germany), following the manufacturer's protocol modified by increasing initial lysis incubation to 1 h at 70°C. Microsatellite sequences were isolated and primers were designed by Ecogenics GmbH (Zurich, Switzerland) using magnetic streptavidin beads and biotin-labeled CT and GT repeat oligonucleotides. The genomic library was enriched for simple sequence repeats (SSRs) prior to 454 sequencing. The SSR content is species dependent, but usually less than 2% of the sequence read pool contains a microsatellite insert without enrichment (Ecogenics GmbH). The SSR-enriched library was analyzed on a Roche 454 platform using GS FLX Titanium reagents (454 Life Sciences, a Roche Company, Branford, Connecticut, USA). Based on fragment length, repeat motif, and no poly N present in each amplicon, the most promising markers were chosen for screening in N. laevigatum and N. parile using a Primer3 core code (available from; Rozen and Skaletsky, 2000). The remaining sequences were not suitable for primer design because flanking regions were absent or too short. In the end, the most promising markers were proved (1) to amplify, and (2) to be repeatable and polymorphic (i.e., bearing at least two alleles per population); the remaining were monomorphic, nonamplifying in some individuals, or too unspecific. For N. laevigatum, a total of 19,431 reads had an average length of 404 base pairs. Of these, 630 contained a microsatellite insert with a tetra- or a trinucleotide of at least six repeat units or a dinucleotide of at least 10 repeat units. Suitable primer design was possible in 265 reads (i.e., read sequence is not starting or ending with the repeat motif). In the case of N. parile, a total of 21,136 reads had an average length of 396 base pairs. Of these, 1386 contained a microsatellite insert with a tetra- or a trinucleotide of at least six repeat units or a dinucleotide of at least 10 repeat units. Suitable primer design was possible in 488 reads. Primers were designed for 12 microsatellite inserts for N. laevigatum and 14 microsatellite inserts for N. parile and tested for polymorphism. Although few simple sequence repeats are expected in bacterial genomes, an additional BLAST search was included for all primer candidates to exclude any that matched the available Nostoc reference genome before testing for polymorphism.

Table 1.

Characteristics of 26 microsatellite loci developed in Nephroma laevigatum and N. parile.


Table 2.

Results of microsatellite screening in four populations of Nephroma laevigatum.a,b


Polymorphisms within the 26 microsatellite loci were determined by Ecogenics GmbH using a test set of 15 samples for each species collected from four Scottish populations (Appendix 1), and using M13-tailed forward primers and universal M13 fluorescent-labeled primer following Schuelke (2000). For these samples, genomic DNA was extracted from dried thallus material with the QIAGEN DNeasy Plant Mini Kit, and loci were PCR amplified in the Veriti 96-Well Thermal Cycler (Applied Biosystems, Carlsbad, California, USA). Primers that showed clear amplification profiles and reliable amplification and that were polymorphic across samples are described in Table 1.

Of the set of 26 polymorphic primers obtained from the polymorphism test, a larger screen of marker variation was completed using four populations for each species (Tables 2, 3), totaling more than 60 samples per species. All PCR products obtained were multiplexed (Table 1). PCR reactions were performed in a total volume of 10 µL containing 1 µL of genomic DNA, 2 µL of water, 5 µL of HotStarTaq (QIAGEN), and 2 µL of primer mix (1 µM each of forward and reverse primer concentration). Cycling conditions were: 15 min at 95°C, 35× (30 s at 94°C, 90 s at 56°C, 60 s at 72°C), and a final extension of 30 min at 72°C. For the fragment analyses, 1 µL of 1:20 diluted multiplex products was used. PCR products were sized on an ABI3730 Genetic Analyzer (Applied Biosystems), and genotypes were assigned with a GeneScan 500 LIZ Size Standard (Life Technologies, Carlsbad, California, USA). Of the initial 26 microsatellite primer pairs, four were not included in the results of this screen due to their low amplification success on the evaluated samples (>8% nonamplifications; NLae01, NLae06, NPar04, NPar14). The number of alleles and unbiased genetic diversity were generated using GenAlEx 6.501 (Peakall and Smouse, 2006; Tables 2, 3). For the sample of 66 specimens for N. laevigatum, all 10 loci evaluated were polymorphic with between three and 13 alleles per locus, a total of 68 alleles across loci, and altogether comprising 46 multilocus genotypes detected over all populations. From the sample of 61 specimens of N. parile, all 12 loci evaluated were polymorphic with two to six alleles per locus, a total number of 48 alleles across loci, and 16 multilocus genotypes detected over all populations. Average gene diversities varied from 0.38 to 0.71 and 0.01 to 0.50 for N. laevigatum and N. parile, respectively (Tables 2, 3). A small sample of specimens (four per species) was repeated to check for reproducibility, and all yielded identical allele lengths.

Table 3.

Results of microsatellite screening in four populations of Nephroma parile.a,b



Twelve and 14 microsatellite primers were characterized to investigate population structure and gene flow in N. laevigatum and N. parile, respectively. The markers demonstrate high-resolution variability at even a relatively small geographic sampling scale, with no shared genotypes among populations for the sexual species, in contrast to the predominantly asexual N. parile. We conclude that the markers can potentially provide insights contributing to an improved understanding of population genetic processes, and they are currently being used to analyze population genetic structure in N. laevigatum and N. parile for variable habitats across a steep climatic gradient.



D. L. Hartl , and A. G. Clark . 1997. Principles of population genetics. Sinauer, Sunderland, Massachusetts, USA. Google Scholar


K. Lohtander , I. Oksanen , and J. Rikkinen . 2002. A phylogenetic study of Nephroma (lichen-forming Ascomycota). Mycological Research 106: 777–787. Google Scholar


A. Lowe , S. Harris , and P. Ashton . 2004. Gene flow and mating system. In A. Lowe , S. Harris , and P. Ashton [eds.], Ecological genetics: Design, analysis, and application, 106–149. Blackwell Publishing, Oxford, United Kingdom. Google Scholar


R. Peakall , and P. E. Smouse . 2006. GenAlEx 6: Genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes 6: 288–295. Google Scholar


S. Rozen , and H. Skaletsky . 2000. Primer3 on the WWW for general users and for biologist programmers. In S. Misener and S. A. Krawetz [eds.]. Methods in molecular biology, vol. 132: Bioinformatics methods and protocols, 365–386. Humana Press, Totowa, New Jersey, USA. Google Scholar


C. Scheidegger , and S. Werth . 2009. Conservation strategies for lichens: Insights from population biology. Fungal Biology Reviews 23: 55–66. Google Scholar


M. Schuelke 2000. An economic method for the fluorescent labelling of PCR fragments. Nature Biotechnology 18: 233–234. Google Scholar


M. R. D. Seaward 2008. Environmental role of lichens. In T. H. Nash III [ed.]. Lichen biology, 274–298. Cambridge University Press, Cambridge, United Kingdom. Google Scholar


E. Sérusiaux , J. C. Villareal , T. Wheeler , and B. Goffinet . 2011. Recent origin, active speciation and dispersal for the lichen genus Nephroma (Peltigerales) in Macronesia. Journal of Biogeography 38: 1138–1151. Google Scholar


C. W. Smith , A. Aptroot , B. J. Coppins , A. Fletcher , O. L. Gilbert , P. W. James , and P. A. Wolseley . 2009. The lichens of Great Britain and Ireland. British Lichen Society, United Kingdom. Google Scholar


S. Werth 2010. Population genetics of lichen-forming fungi: A review. Lichenologist (London, England) 42: 499–519. Google Scholar


S. Werth , S. Cheenacharoen , and C. Scheidegger . 2014. Propagule size is not a good predictor for regional population subdivision or fine-scale spatial structure in lichenized fungi. Fungal Biology 118: 126–138. Google Scholar


Appendix 1.

Location information, number of individuals sampled, and voucher specimens of Nephroma species used in this study. All specimens were sampled as epiphytes from broadleaved woodlands.a



[1] The authors thank Markus Ruhsam, Laura Forrest, Michelle Hollingsworth, Alfredo García-Femández, María Prieto, and The GenePool (University of Edinburgh) for technical assistance. We also thank the Scottish Natural Heritage (SNH) and the landowners for access permission. This study was funded by a Marie Curie Intra-European Fellowship within the 7th European Community Framework Program to R.B. (PIEF-GA-2011-299330).

Rocío Belinchón, Christopher J. Ellis, and Rebecca Yahr "Microsatellite Loci in Two Epiphytic Lichens with Contrasting Dispersal Modes: Nephroma laevigatum and N. parile (Nephromataceae)," Applications in Plant Sciences 2(11), (3 November 2014).
Received: 21 August 2014; Accepted: 21 September 2014; Published: 3 November 2014

Back to Top