| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Universitetet i Bergen, Botanisk Institutt, Allégaten 41, N-5007 Bergen, Norway
| ABSTRACT |
|---|
|
|
|---|
A 753771 bp long intronic sequence from the mitochondrial cox 1 gene of Cladonia subcervicornis (Cladoniaceae, Lecanorales, Ascomycota) was amplified with newly designed PCR primers. The cox 1 intron sequence, which apparently has not been used for phylogenetic or population genetic research in fungi, displays high infraspecific variation. Sequences were obtained from 124 specimens from four neighboring localities in coastal Hordaland, western Norway. An exact test of population differentiation and population pairwise fixation indices FST show significantly reduced gene flow between the northernmost locality and the other three populations. Although Cladonia subcervicornis frequently produces apothecia, we conclude that dispersal by ascospores over long distances is rather ineffective in this species.
Key words: Cladonia, dispersal, lichenized Ascomycota, maturase, PCR primers
| INTRODUCTION |
|---|
|
|
|---|
The development of PCR-based methods is beginning to change this picture because it finally grants lichenologists access to molecular markers, but other problems still impede population genetic studies. Lichens are symbiotic organisms, and DNA extractions of all but the largest thalli almost inevitably contain both fungal and algal or cyanobacterial DNA. Typical datasets in population genetics often comprise hundreds of individuals, but cost-effective methods, such as RAPD-PCR or ISSR-PCR, cannot discriminate among DNA of the symbionts. Unless algal-free mycobiont cultures are used, a costly and laborious approach, the results potentially are flawed. To our knowledge, fungal specific microsatellite primers have not been published for lichens. As a result, nucleotide sequences obtained by fungal-specific PCR-primers are the only reasonable alternative, in most cases. The high cost associated with this approach has not promoted lichen-population genetics, and the number of available marker sequences consequently is low (Zoller et al 1999
, Printzen 2002
, Printzen and Ekman 2002
). Other studies addressing questions of infraspecific variation in lichens with the help of nucleotide sequences include DePriest (1993
, 1994
) and Beard and DePriest (1996)
, who employed PCR-RFLP to demonstrate variability of number, size and position of Group 1 introns in the nuclear SSU rDNA of Cladonia chlorophaea and C. subtenuis within populations. However, most nonribosomal gene loci used in lichen phylogeny and population genetics do not contain introns, which makes it unlikely that sufficient variability can be uncovered by PCR-RFLP.
The first objective of our study was to extend the array of molecular markers for population studies in lichens. Fungal mitochondrial genomes are evolutionarily derived and are characterized by accelerated sequence divergence (Gray et al 1999
). Considerable length variation has been observed among different species, largely due to the frequent occurrence of Group 1 and Group 2 introns. Both types of introns are self-splicing, mobile elements. Many of them contain ORFs encoding for homing endonucleases, maturases that assist in the splicing of introns, or proteins that function as both (Gimble 2000
). An estimated 30% of Group 1 introns contain such ORFs, while they seem to be less common in Group 2 introns (Chevalier and Stoddard 2001
). Group 1 introns have been shown to spread rapidly within populations and among species. This mobility causes coding as well as noncoding parts of the intron sequence to be relatively unaffected by selection against deleterious mutations. Instead, they undergo regular cycles of transmission, fixation, degeneration and loss (Goddard and Burt 1999
). Mitochondrial sequences thus can be assumed to contain the amount of variability necessary for studies at the population level. The few available studies indeed showed appreciable amounts of infraspecific variation (Zoller et al 1999
, Printzen 2002
).
Although lichens produce an astonishingly varied array of diaspores, very little is known about their dispersal. Sexual reproduction in lichens occurs via ascospores or basidiospores, while asexual dispersal can occur either through fungal propagules, such as conidia or structures containing both fungal and algal cells. Büdel and Scheidegger (1996)
give a comprehensive overview of the modes of reproduction in lichens. Because fungal ascospores are orders of magnitude smaller than seeds of vascular plants, their dispersal is unlikely to be limited physically. This indeed was supported by findings of lichen propagules in pollen traps on the Atlantic Ocean (Harmata and Olech 1991
). These findings, however, did not prove effective dispersal, i.e., dispersal of propagules and establishment in a new environment (Cain et al 2000
). Indirect evidence for effective spore dispersal was provided by Tibell (1994)
and Wedin (1995)
, who found statistically significant correlations between small spores and wide geographical distribution in species of the Caliciales and Sphaerophoraceae. On the other hand, investigations of isozyme patterns of Umbilicaria mammulata (Hageman and Fahselt 1992
) showed more limited genetic exchange among subpopulations than within, even within a rather limited geographical area. But because Umbilicaria mammulata is very rarely fertile and did not produce vegetative propagules in the study area, the results shed little light on the ability of diaspore dispersal in lichens. In any case, the authors challenged the widespread opinion that lichens disperse easily and effectively and, for the first time, demonstrated that lichen dispersal could be assessed indirectly by measuring genetic differentiation between populations.
To contribute data on a lichen species that frequently produces ascomata and to test the suitability of our new marker for population studies, we compared populations of the macrolichen Cladonia subcervicornis in a windswept, open environment on a local scale. Lichens of the genus Cladonia belong to the most conspicuous and most frequently collected lichens. Cladonia subcervicornis occurs in oceanic parts of northern and western Europe from the British Isles to northern Norway (Purvis and James 1992
, Santesson 1993
). In Norway it occurs predominantly in coastal heaths, where it forms compact cushions over humus in rock crevices. Cladonia subcervicornis is common in these habitats and may dominate the vegetation of exposed rocks. Strong winds are recorded throughout the year near the coast of southwestern Norway. Calm conditions represent only 6% of all observations at Flesland weather station, and predominant winds change from southeasterly in autumn and winter to northwesterly in summer (Bjørbæk 1993
). Ascospores thus should disperse widely in open coastal habitats.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Fungal specific primers were designed following the method proposed by Carbone and Kohn (1999)
. Primers were checked for possible hairpin structures or dimer formation with Primer Premier MINI version 4.04 (Premier Biosoft, URL: http://www.PremierBiosoft.com). Out of eight forward and six backward primers originally designed, a combination of these two primers amplified part of the cytochrome oxidase subunit 1 sequence in various Cladonia species: (5959F5') 5'-TCT TAA CGT TGC TGT ATG CTG-3'; (6711R3') 5'-GAA CCG AAA CTA GTA GAA CCA TA-3'. Primer numbers refer to positions in the cox 1 sequence of Neurospora crassa (Genbank accession X14669).
A total of 10 µL of DNA extractions were used in 50 µL PCR reactions together with 5 µL HerculaseTM Reaction Buffer (Stratagene 60026052), 1 mM MgCl2, 2.5 mM each dNTP, 2.5 U polymerase (HerculaseTM, Stratagene 600264), and 0.8 µM each of the 5'- and 3'- primers. PCR cycling conditions were: 95 C (5'), six cycles of 95 C (45''), 5246 C (touchdown, 45''), 72 C (1'45''), 34 cycles of 95 C (30''), 46 C (30''), 72 C (1'45''), and a final extension of 72 C (10'). PCR products were purified with the QIAquickTM PCR Purification Kit (QIAGEN 28106) or the QIAquickTM Gel Extraction Kit (QIAGEN 28706). Purified DNA was labelled with the BigDye TerminatorTM Kit (Applied Biosystems) and cycle sequenced at 94 C (30''), and 29 cycles at 95 C (15''), 50 C (15''), 60 C (4'). Sequences were determined on an ABI PRISM® 3700 DNA Analyzer (Applied Biosystems), assembled with SeqManTM II, version 4.05 (DNASTAR) and manually aligned in BioEdit, version 4.8.8 (Hall 1999
). Sequences of each observed haplotype are deposited in GenBank (accession numbers AY148463AY148469).
Sequence identification
To verify the identity of the nucleotide sequences amplified by us, a BLAST search in GenBank was performed. Cox 1 sequences from lecanoralean ascomycetes have not been published before. To rule out the possibility that we were sequencing populations of parasitic fungi, we also sequenced Cladonia polydactyla (Flörke) Sprengel, collected in different habitats, for comparison. Cladonia sequences were aligned with ClustalW as implemented in Bioedit version 4.8.8 (Hall 1999
). One ORF within our sequences was detected and translated with the same program. The protein sequence of the ORF was subsequently subjected to a protein BLAST search to find similar protein sequences.
Data analysis
To compare the variability of cox 1 intron sequences with that of other sequences that have been used for infraspecific studies in lichens, average pairwise differences were calculated for the entire dataset using SITES (Hey and Wakeley 1997
) and compared with unpublished datasets of IGS and ITS sequences from Cavernularia hultenii.
The recovered sequences were collapsed into haplotypes and the occurrence of the haplotypes noted for each locality. Genealogy information was not used in the subsequent data analysis because only seven different haplotypes were observed (see below). The dataset was tested for population subdivision by two different approaches using the program Arlequin version 2.001 (Schneider et al 2001
). Population pairwise fixation indices FST were calculated and their significance tested with a nonparametric permutation approach with 1000 permutations of haplotypes among localities. An exact test of population differentiation (Raymond and Rousset 1995
) was performed and significance tested using Markov chain Monte Carlo with 1 000 000 steps, 50 000 of which were dememorization steps.
| RESULTS |
|---|
|
|
|---|
|
= 0.02173, Var (
) = 1.18 x 10-4). Sequences of the internal transcribed spacer (ITS) and intergenic spacer (IGS) regions of the nuclear ribosomal DNA studied in 150 (ITS), viz. 135 (IGS) collections of Cavernularia hultenii from Newfoundland, Norway and the Pacific Northwest of North America were considerably more conserved (ITS:
= 0.00346, Var (
) = 5.06 x 10-6; IGS:
= 0.00669, Var (
) = 1.67 x 10-5; Printzen and Ekman unpubl. data). The cox1 intron sequences presented here apparently detect more infraspecific variability than other known fungal DNA sequences. Despite the high variability, only seven different haplotypes were recognized in our samples. The genetic relationships of these haplotypes and their distribution among the four sampled subpopulations is shown in Fig. 3. The results of the test for population differentiation are shown in Table I. Population pairwise fixation indices, as well as the exact test, suggest significant population differentiation between Alvøyni and the other three populations. Neither test indicates significant differentiation between the other three localities.
|
|
| DISCUSSION |
|---|
|
|
|---|
Three of the four populations studied by us display little genetic differentiation, with pairwise FST values well below 0.05. This result is not surprising in the case of a fertile lichen in an open environment. The northernmost population on Alvøyni, however, is genetically very different from the other three populations (Table I). Such pronounced population differentiation in a fertile lichen species on a small geographic scale is unexpected. The question remains whether this differentiation is really the result of reduced gene flow among populations. Inferring levels of gene flow from FST values relies on demographic and population genetic models, the applicability of which is difficult to judge (Ouborg et al 1999
). These difficulties are especially pronounced in lichens, where information on reproductive biology, ecology and population dynamics either is lacking or rudimentary. McCauley et al (1995)
demonstrated that, in plant metapopulations with high extinction-recolonization dynamics, the genetic differentation between populations indicates much lower dispersal than is actually present.
The few reports on growth rates in Cladonia (e.g., Jahns and Ott 1982
) indicate that all specimens collected for this study were several years old. Consequently, it is unlikely that the population at Alvøyni was established recently. The fact that the more-or-less barren rock surfaces on which C. subcervicornis grows are not easily colonized by vascular plants further supports the fact that populations of C. subcervicornis are not subject to high extinction-recolonization dynamics. Regardless of the underlying mathematical model and the fact that we cannot determine Nmthe number of migrants per generation between populationsgenetic differentiation hence is likely to result from limited dispersal.
Figure 3 shows that the haplotypes sampled by us cluster into three groups of two closely related alleles. These groups are separated by 621 steps. A seventh, single haplotype is separated from the rest of the network by 14 steps. Cryptic speciation recently has been detected in the Letharia vulpina-columbiana species complex (Kroken and Taylor 2001
). Thus it is conceivable that the individuals sampled for this study in fact belong to different phylogenetic species. Instead of observing genetic differentiation among populations, in this case we would observe differences in the distribution of four species. Cladonia subcervicornis occurs over a range of 3000 km along the European coast. In the absence of any ecological and morphological variability among localities, the probability of encountering four different cryptic species by chance when sampling over a latitudinal range of just 45 km appears minute, unless the number of cryptic species is very high or the cryptic species are all sympatric. It cannot be concluded from the single example of Letharia that cryptic speciation is common among lichens. In fact, the "cryptic" species detected within that genus by Kroken and Taylor (2001)
already had been described as infraspecific taxa (under L. vulpina) based on morphological differentiation (Schade 1955
). Without further evidence, we consider it more likely that we are dealing with populations of a single species.
Our results corroborate a study by Hageman and Fahselt (1992
), who reported among-stand variation of isozyme patterns on a similar geographical scale. Those authors used Umbilicaria mammulata, a species that rarely is found fertile and that lacks specialized vegetative propagules, as their study object. Our results indicate that ascospores do not help to disperse C. subcervicornis over longer distances, although it frequently is found with apothecia. Two recent studies by Sillett et al (2000
) and Hilmo and Såstad (2001
) add to our results. In transplantation and growth experiments, these authors showed that some lichens are restricted to old-growth forests, not because of ecological specialization but because of dispersal limitations. Perhaps effective dispersal is much rarer among many lichens than would be predicted by the abundant production of sexual and asexual diaspores. Using molecular markers in studies on the population dynamics of lichens could yield new insights into their dispersal biology. With only few nucleotide markers for infraspecific lichen studies available, cox 1 intron sequences could prove a valuable addition to this small arsenal of tools.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Accepted for publication November 25, 2002.
| LITERATURE CITED |
|---|
|
|
|---|
Bjørbæk G., 1993 Nasjonalatlas for Norge, kartblad 3.1.8. Vind 1:3,5 mill. Statens kartverk
Brown D, Cureton-Brown MM, Kershaw KA., 1989 Electrophoretic variation in a population of the lichen Peltigera rufescens (Weis.) Mudd. New Phytologist 112:561-567
Büdel B, Scheidegger C., 1996 Thallus morphology and anatomy. In: Nash III TH, ed. Lichen biology. Cambridge: Cambridge University Press. p 3764
Cain ML, Milligan BG, Strand AE., 2000 Long-distance dispersal in plant populations. Am J Bot 87:1217-1227
Carbone I, Kohn LM., 1999 A method for designing primer stes for speciation studies in filamentous ascomycetes. Mycologia 91:553-556
Chevalier BS, Stoddard BL., 2001 Homing endonucleases: structural and functional insight into the catalysts of intron/intein mobility. Nucleic Acids Res 29:3757-3774
Cummings DJ, McNally KL, Domenico JM, Matsuura ET., 1990 The complete DNA sequence of the mitochondrial genome of Podospora anserina. Curr Gen 17:375-402[Medline]
DePriest PT., 1993 Small subunit rDNA variation in a population of lichen fungi due to optional group-I introns. Gene 134:67-74[Medline]
. 1994 Variation in the Cladonia chlorophaea complex II: ribosomal DNA variation in a Southern Appalachian population. Bryologist 97:117-126
Fahselt D, Hageman C., 1983 Isozyme banding patterns in two stands of Cetraria arenaria Kärnef. Bryologist 86:129-134
Gimble FS., 2000 Invasion of a multitude of genetic niches by mobile endonuclease genes. FEMS Microbiol Lett 185:99-107[Medline]
Goddard MR, Burt A., 1999 Recurrent invasion and extinction of a selfish gene. Proc Natl Acad Sci USA 96:13880-13885
Gonzalez P, Barroso G, Labarère J., 1998 Molecular analysis of the split cox 1 gene from the Basidiomycota Agrocybe aegerita: relationship of its introns with homologous Ascomycota introns and divergence levels from common ancestral species. Gene 220:45-53[Medline]
Gray MW, Burger G, Lang BF., 1999 Mitochondrial evolution. Science 283:1476-1481
Hageman C, Fahselt D., 1984 Intraspecific variability of isozymes of the lichen Umbilicaria mammulata. Can J Bot 62:617-623
. 1990 Enzyme electromorph variation in the lichen family Umbilicariaceae: within-stand polymorphism in umbilicate lichens of Eastern Canada. Can J Bot 68:2636-2643
. 1992 Geographical distribution of enzyme polymorphisms in the lichen Umbilicaria mammulata. Bryologist 95:316-323
Hall TA., 1999 BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 41:95-98
Harmata K, Olech M., 1991 Transect for aerobiological studies from Antarctica to Poland. Grana 30:458-463
Hey J, Wakeley J., 1997 A coalescent estimator of the population recombination rate. Genetics 145:833-846[Abstract]
Hilmo O, Såstad SM., 2001 Colonization of old-forest lichens in a young and an old boreal Picea abies forest: an experimental approach. Biol Cons 102:251-259
Jahns HM, Ott S., 1982 Flechtenentwicklung an dicht benachbarten Standorten. Herzogia 6:201-241
Kroken S, Taylor JW., 2001 A gene genealogical approach to recognize phylogenetic species boundaries in the lichenized fungus Letharia. Mycologia 93:38-53
Mattsson JE., 1994 Lichen proteins, secondary products and morphology: a review of protein studies in lichens with special emphasis on taxonomy. J Hattori Bot Lab 76:235-248
McCauley DE, Raveill J, Antonovics J., 1995 Local founding events as determinants of genetic structure in a plant metapopulation. Heredity 75:630-636
Ouborg NJ, Piquot Y, van Groenendael JM., 1999 Population genetics, molecular markers and the study of dispersal in plants. J Ecol 87:551-568
Paquin B, Laforest M-J, Forget L, Roewer I, Zhang W, Longcore J, Lang BF., 1997 The fungal mitochondrial genome project: evolution of fungal mitochondrial genomes and their gene expression. Genetics 31:380-395
Printzen C., 2002 Fungal specific primers for the amplification of mitochondrial LSU in lichens. Mol Ecol Notes (2002) 2:130-132
, Ekman S., 2002 Genetic variability and its geographical distribution in the widely disjunct Cavernularia hultenii. Lichenologist 34:101-111
Purvis OW, James PW., 1992 Cladonia. In: Moore DM, ed. The Lichen Flora of Great Britain and Ireland. London: Natural History Museum Publications. p 188210
Raymond M, Rousset F., 1995 An exact test for population differentiation. Evolution 49:1280-1283
Santesson R., 1993 The lichens and lichenicolous fungi of Sweden and Norway. Lund: SBT-förlaget. 240 p
Schade A., 1955 Letharia vulpina (L.) Vain.II. Ihr Vorkommen in der Neuen Welt und ihr Verhältnis zu Letharia californica (Lev.) Hue em. DR. Feddes Repert 58:179-197
Schneider S, Roessli D, Excoffier L., 2001 Arlequin: A software for population genetics data analysis. Ver. 2.001. Genetics and Biometry Lab, Dept. of Anthropology, University of Geneva
Sillett SC, McCune B, Peck JE, Rambo TE, Ruchty A., 2000 Dispersal limitations of epiphytic lichens result in species dependent on old-growth forests. Ecol Appl 10:789-799
Tibell L., 1994 Distribution patterns and dispersal strategies of Caliciales. Bot J Linn Soc 116:159-202
Wedin M., 1995 The lichen family Sphaerophoraceae (Caliciales, Ascomycotina) in temperate areas of the southern hemisphere. Symb Bot Upsal 31:1-102
Zoller S, Lutzoni F, Scheidegger C., 1999 Genetic variation within and among populations of the threatened lichen Lobaria pulmonaria in Switzerland and implications for its conservation. Mol Ecol 8:2049-2059[Medline]
This article has been cited by other articles:
![]() |
S. Werth and V. L. Sork Local genetic structure in a North American epiphytic lichen, Ramalina menziesii (Ramalinaceae) Am. J. Botany, May 1, 2008; 95(5): 568 - 576. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. R. Naesborg Taxonomic revision of the Lecania cyrtella group based on molecular and morphological evidence. Mycologia, May 1, 2008; 100(3): 397 - 416. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |