Mycologia
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

DOI: 10.3852/mycologia.98.1.57
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zhou, Q.-M.
Right arrow Articles by Wei, J.-C.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Zhou, Q.-M.
Right arrow Articles by Wei, J.-C.
Agricola
Right arrow Articles by Zhou, Q.-M.
Right arrow Articles by Wei, J.-C.
Mycologia, 98(1), 2006, pp. 57-67.
© 2006 by The Mycological Society of America

A study of the genetic variability of Rhizoplaca chrysoleuca using DNA sequences and secondary metabolic substances


Qi-Ming Zhou
Shou-Yu Guo 1
Man-Rong Huang
Jiang-Chun Wei

     Systematic Mycology & Lichenology Laboratory, Institute of Microbiology, Academia Sinica, Bei-yi-tiao No. 13, Zhong-guan-cun, Beijing, 100080, P.R. China

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 

Specimens of Rhizoplaca chrysoleuca from Mount Wuling can be divided into two distinct groups based on obvious differences in morphological characters. Here we investigated 26 specimens of R. chrysoleuca from Mount Wuling, 10 specimens of this species from other areas and seven specimens of other Rhizoplaca species by analyzing morphology, chemistry and genetics. Nine chemotypes were detected among the specimens of R. chrysoleuca from Mount Wuling, and five of them were reported for the first time. Based on the ITS phylogenetic analysis, the chemotypes and the insertion distribution patterns in SSU rDNA, the samples of R. chrysoleuca from Mount Wuling were grouped in two distinct clades corresponding to two phenotypic groups and no gene flow was detected between these two groups. Our results establish all individuals of Rhizoplaca chrysoleuca are conspecific although some populations have been isolated on Mount Wuling, indicating that they are in the process of speciation. Our study also reveals that the relationships between genotypes and chemotypes are complicated and should be avoided, and we instead recommend using single individuals or few individuals from the same site to represent the population or whole species in systematics study. The results also indicate that Rhizoplaca chrysoleuca might provide a good model for studying the speciation of saxicolous lichenized fungi.

Key words: chemotype, group I intron, insertion distribution pattern, ITS, speciation, SSU rDNA


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Lichens are well known for their secondary metabolic substances that rarely occur elsewhere in nature (Culberson 1969Go, Culberson et al 1988Go). For more than 100 y lichenologists have used such unique organic compounds as an aid in lichen identification and classification (e.g. Narui et al 1996Go). Recently there has been conclusive evidence that the major secondary metabolism observed in lichens is due to mycobionts (Culberson and Armaleo 1992Go) and is inherited as a genetic trait (Culberson et al 1988Go). Some studies also have suggested that the recognition of chemically defined species could be supported by their distinct geographical distributions and ecological preferences (Hawksworth 1976Go). One of the classical examples was provided by Cladonia chlorophaea complex (DePriest 1992Go, 1994Go). However chemical similarity or difference is not an infallible indicator of systematic relationships because chemistry is known to be variable in some species (Culberson 1986Go, Culberson et al 1988Go) and sporadically shared by divergent groups. This phenomenon has been reported in Rhizoplaca with several chemotypes detected for R. chrysoleuca, R. melanophalma, R. peltata, and R. subdiscrepans (Leuckert et al 1977Go, Wei 1984Go, Ryan and Nash 1997Go, Arup and Grube 2000Go).

Molecular evidence, particularly variations in nucleotide sequences, offers a large number of potentially informative characters. The ITS regions were used generally in the comparisons of closely related species and genera, and their sequences information has been one of the primary criteria in ascertaining relationships at the specific level (e.g. Groner and LaGreca 1997Go, Goward and Goffinet 2000Go). Recently the ITS region also was used in population analysis (e.g. Presa et al 2002Go, Koch et al 2003Go). Other genes also have been applied in fungal systematics including chitin synthase genes (e.g. Mehmann et al 1994Go, Kroken and Taylor 2001Go), the ß-tubulin gene (e.g. Koufopanou et al 2001Go, Myllys et al 2001Go, Cruse et al 2002Go), the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene (e.g. Myllys et al 2002Go) and the microsatellite markers (e.g. Carter et al 2001Go, Walser et al 2005Go). These genes were used mostly in studying population structure or cryptic species. Studies on population differentiation and isolation in lichens have been carried out recently, especially for the corticolous and epiphytic lichens (e.g. Graphis scripta [Murtagh et al 2000Go], Letharia vulpina [Högberg et al 2002Go] and Lobaria pulmonaria [Walser et al 2005Go]).

Rhizoplaca is a small genus of saxicolous lichens with nine species usually containing usnic acid (Ryan and Nash 1997Go, Arup and Grube 2000Go, Kirk et al 2001Go). Species of Rhizoplaca are more or less umbilicate and fixed to the substrate or more rarely are fruticose or unattached, vagrant lichens. This genus formerly was recognized by Choisy (1929)Go, who described it under the synonym Omphalodina M. Choisy with R. chrysoleuca (Smith) Zopf (= Omphalodina chrysoleuca) as the type species. R. chrysoleuca is famous for its nearly worldwide distribution and its ability to produce ice nucleation active proteins (Kieft 1988Go, Kieft and Ruscetti 1990Go).

Sometimes populations of R. chrysoleuca that differ in morphology and secondary products can be found growing intermixed with typical material. For example on Mount Wuling (Hebei Province, north China, ca. 110 km east of Beijing), populations comprising various morphotypes and several distinct chemotypes are found. The presence of a particular nrDNA sequence pattern being unique to particular chemotypes would support chemotypic divergence. In contrast the random association of nrDNA sequence patterns and chemotypes would support chemotypic polymorphism within the species (DePriest, 1994Go). The differences in nrDNA sequences combined with the chemical variability were used to identify genetically distinct individuals and genetic diversity in the present study. We investigated the morphotypic and chemotypic polymorphisms of R. chrysoleuca from a molecular-genetic view and studied the population structure of R. chrysoleuca on Mount Wuling.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Samples.— – Lichen samples selected for this study are listed (TABLE IGo). Forty-three Rhizoplaca specimens were selected from five localities in Hebei Province, Jilin Province, Xinjiang Uygur Autonomous Region and Nei Mongol Autonomous Region, China (FIG. 1AGo). Among them 26 R. chrysoleuca samples were taken from nine lichen mats (each ca. 5–10 cm diam) in two sites at different elevation (1700 m and 2100 m, respectively) on Mount Wuling, Hebei Province, north China (FIG. 1B,CGo). All samples were carefully brushed to remove possible epiphytic contaminants. These samples are deposited in the Lichen Section of the Mycological Herbarium, Academia Sinica (HMAS-L), as voucher specimens for the ITS and SSU rDNA sequences.


View this table:
[in this window]
[in a new window]
 
TABLE I. Information for 43 samples included in this study

 

Figure 1
View larger version (10K):
[in this window]
[in a new window]
 
FIG. 1. Sketch maps of collecting sites. A, five localities in China (• Mount Wuling, Hebei; {blacktriangledown} Mount Arxan, Motianling and Mount Jiguan, Inner Mongol; {blacktriangleup} Mount Tuding, Jilin; {diamondsuit} Jimsar, Xinjinag; {blacksquare} Hanas and Guanyuting, Xinjiang). B and C are the loci at 2100 m and 1700 m elevations respectively. Each dot represents a single mat. Samples were taken from every mat. Samples’ serial numbers indicate their positions. For example, R. chrysoleuca111 was collected from 1-1 mat in B locus (2100 m), and R. chrysoleuca246 was collected from 2–4 mat in C locus (1700 m).

 
Morphology.— – Morphology of the thalli of R. chrysoleuca was examined first. Observations of external features were made under a dissecting microscope (Zeiss SV11). Two main morphotypes (TABLE IIGo) were established based on the morphological data that were obtained from ca. 200 specimens (including the samples used in the molecular study) collected in different regions (most of them from China and 17 specimens from Europe and North America).


View this table:
[in this window]
[in a new window]
 
TABLE II. Morphological variations in the Rhizoplaca chrysoleuca

 
Chemistry.— – Individual thalli of 36 R. chrysoleuca (TABLE IGo) were examined chemically (chemotypes given in TABLE IIIGo) with standard thin-layer chromatographic (TLC) methods (Culberson 1972Go). The compounds in each thallus were extracted with acetone, and the extracts were spotted on 10 x 20 cm silica gel plates (supported by glass, Qingdao Ocean Chemical Factory) and developed in solvent system C. Spots were visualized under UV light at 254 nm before spraying with 10% sulfuric acid and heating at 110 C for 10 min. Lethariella cladonioides which contains norstictic acid and atranorin was selected as a control sample.


View this table:
[in this window]
[in a new window]
 
TABLE III. The chemotypes of the Rhizoplaca chrysoleuca in Chinaa

 
PCR amplification and sequencing.— – Total DNA was extracted from all 43 individual Rhizoplaca thalli (TABLE IGo) using the modified CTAB method (Rogers and Bendich 1988Go). DNA extracts were used for PCR amplification of the ITS regions including the 5.8S gene of the nuclear rDNA with ITS5 and ITS4 (White et al 1990Go) as primers. PCR reactions were performed in a DNA Thermal Cycler (Biometra) as follows: initial denaturation at 95 C for 3 min, followed by 35 cycles of 30 s denaturation at 94 C, 45 s annealing at 58 C, 1 min extension at 72 C and completed with a final 8 min extension at 72 C. Products were purified with Gel Extraction Mini Kit (SABC). Sequencing reactions were carried out by Shanghai Genecore Corp. with an ABI 3700 Sequencer. Both complementary strands of each sample were sequenced.

Partial SSU rDNA sequences were obtained from all R. chrysoleuca samples after amplifying with primer pair NS21UCB and NS24UCB (Gargas and Taylor 1992Go) except for 5A and 17A, but the length of PCR products differed significantly among individuals (TABLE IGo). The fragments obtained from samples 111, 211, 221, 223, G8, 231, 234, 7A, 8A, 9A, 10A and 15A were chosen to represent the different length products and sequenced by Shanghai Bioasian Corp. using an ABI 377 Sequencer. However only the 5' and 3' flank sequences of the PCR product for 15A were obtained because it failed in walking reactions. The near complete SSU rDNA sequence of R. chrysoleuca 144 had been obtained in a previous study (AY530888 [GenBank] ).

Data analysis.— – ITS sequences (including ITS1, 5.8S rDNA and ITS2) of Rhizoplaca acquired in this study and those from GenBank (Parmelia sulcata, AF410840 [GenBank] ; R. cerebriformis U312, AF159942 [GenBank] ; R. chrysoleuca U192, AF159924 [GenBank] ; R. chrysoleuca U302, AF159940 [GenBank] ; R. cylindrica U305, AF159941 [GenBank] ; R. haydenii U285, AF159937 [GenBank] ; R. idahoensis U313, AF159943 [GenBank] ; R. melanophthalma U219, AF159929 [GenBank] ; R. melanophthalma U278, AF159934 [GenBank] ; R. melanophthalma U281, AF159935 [GenBank] ; R. peltata U198, AF159925 [GenBank] ; R. peltata U282, AF159936 [GenBank] ; R. subdiscrepans U387, AF163113 [GenBank] ; R. subidahoensis U314, AF159944 [GenBank] ) were aligned with DNAMAN4.0 (Lynnon Biosoft) without the flanking regions of the small subunit and large subunit rDNA. Phylogenetic analysis was executed with software MEGA2 (Kumar et al 2001Go) and Parmelia sulcata was selected as the outgroup taxon. We used the Kimura 2-parameter model and retained gaps initially but excluded them in the pairwise distance estimation. The neighbor joining (NJ) method was used in constructing the phylogenetic tree, and the reliability of the inferred tree was tested by 1000 bootstrap replications.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Morphology and chemistry.— – Two main morphotypes (TABLE IIGo) were identified in the specimens examined. Morphotype I was characterized by the gray thalli and narrow lobes or granules and morphotype II by the yellow-green thalli and wide lobes. However morphotype I was found only in the samples from Mount Wuling (FIG. 2Go) in the specimens we examined. A total of 12 chemotypes were detected in the present study. Neither morphotype correlated with the two major groups of chemotypes (TABLE IIIGo). The first major group of chemotypes (gray thallii and partial yellow-green thalli, lacks pseudoplacodiolic acid and unidentified substance 2, Group 1 in TABLE IIIGo) includes chemotypes I to X. The second major group of chemotypes (yellow-green thalli from Mount Wuling, contains pseudoplacodiolic acid and unidentified substance 2, Group 2 in TABLE IIIGo) includes chemotypes XI and XII. Chemotypes VII, V and XI correspond respectively to the chemotypes 1, 2 and 4 in Leuckert et al (1977)Go. Whereas chemotypes 3 (usnic + pseudoplacodiolic acids), 5 (usnic + placodiolic + psoromic acids) and 6 (usnic + pseudoplacodiolic + psoromic acids) detected by Leuckert et al (1977)Go for R. chrysoleuca were not found in our specimens examined. Chemotypes I, II, III, VI and XII were reported for the first time. It is interesting that almost all samples collected from Mount Wuling lack lecanoric acid, which is present in most samples from other areas.


Figure 2
View larger version (53K):
[in this window]
[in a new window]
 
FIG. 2. The NJ tree inferred from ITS data. The individuals that are marked with same tag were taken from the same area and these tags are corresponding to those in FIG. 1Go (both "{circ}" and "•" indicate that samples were taken from Mount Wuling, and "{circ}" refers to samples with gray thallus). The numbers in each node represent bootstrap support values, the absence of a number indicating <50% bootstrap support. (HB = Hebei, NM = Nei Mongol, XJ = Xinjiang, JL = Jilin). The morphological types, the chemotypes and the insertion distribution patterns are listed at the right columns.

 
ITS phylogenetic analysis.— – A total of 43 ITS sequences were determined from the mycobiont of Rhizoplaca. The ranges of the length of the ITS region in the samples are 506–512 bp. The samples of R. chrysoleuca from Mount Wuling formed two distinct well-supported clades (clade I and clade II, marked respectively with {circ} and • in FIG. 2Go). Little variations in ITS rDNA sequences were found within either group. The R. chrysoleuca samples from other areas did not form monophyletic groups according to their collection sites; R. chrysoleuca G8 was an exception because only one collection was obtained from Jimusaer, Xinjiang.

R. chrysoleuca samples14A and 15A grouped with U192, and this clade was supported by a 99% bootstrap value. With 97% support other R. chrysoleuca samples appeared as a sister branch to this clade. Although the nucleotides differences in ITS rDNA were obvious, all R. chrysoleuca individuals grouped together and this topology has 100% bootstrap support. R. melanophthalma, R. cerebriformis, R. subidahoensis, R. cylindrical, R. idahoensis and R. haydenii composed the core group of Rhizoplaca (Arup and Grube 2000Go) with 99% bootstrap support. R. subdiscrepans and R. chrysoleuca appeared as a sister taxon to the core group with 99% support. R. peltata, on the other hand, was on a separate branch.

Insertion distribution patterns in SSU rDNA.— – Four types of length polymorphism were observed in the SSU rDNA of R. chrysoleuca as revealed by the PCR reaction (TABLE IGo), and the result of sequencing indicated that it was caused by different numbers of insertions. Twelve PCR products were sequenced in this study. Together with R. chrysoleuca 144 whose SSU rDNA sequence had been determined in a previous study, partial SSU rDNA sequences were obtained for 13 individuals. Because the sequences of some individuals were identical only one was submitted to GenBank (TABLE IGo, the samples sharing the same accession number had identical sequences). The result of sequencing showed that the same length fragments were caused by homological insertions in the same positions. Four kinds of insertion patterns were observed corresponding to four types of length polymorphism (FIG. 3Go). The insertions in position 1462 (according to Saccharomyces cerevisiae with accession number Z75578 [GenBank] ) were spliceosomal introns according to their terminal-sequences (Cubero et al 2000Go); while the rest were group I introns according to their secondary structures (Cech 1988Go). The group I introns in position 1779 were determined in the ITS region sequencing because primer NS24UCB was ahead of this site. The 5' and 3' flank sequences of the PCR product from R. chrysoleuca 15A indicated that there was one spliceosomal intron in position 1462 and no group I introns in positions 999 or 1779 but there was at least one other insertion according to the length of the PCR product in the region not sequenced (FIG. 3DGo).


Figure 3
View larger version (11K):
[in this window]
[in a new window]
 
FIG. 3. Insertions in the different positions of the SSU rDNA in R. chrysoleuca. An inverted triangle indicates the presence of an intron in a particular position. The number below the horizontal line gives an indication of the nucleotide’s position (according to Saccharomyces cerevisiae with accession number Z75578). The region which was not sequenced is shown with the broken line.

 
The common parts of the SSU rDNA for these individuals (not including 15A) were compared with DNAMAN 4.0 (after excluding the introns) and there were only 14 variable positions among about 900 nucleotides, and the number of parsimony-informative characters was just three (data not shown). As a result the data from SSU rDNA was not used in constructing phylogenetic trees.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
After examining many specimens of R. chrysoleuca from different regions it was obvious to us that this species is quite variable in thallus morphology and chemistry. The morphological differences among R. chrysoleuca individuals were so obvious that we could readily divide them into two main groups merely on the basis of the color of the upper surface and the size of the lobe. Although Arup & Grube (1998Go, 2000)Go considered that the growth form of the thallus was highly variable within Rhizoplaca, the distinction between gray and green individuals of R. chrysoleuca is so marked that one would suspect that they are different species.

The insertion distribution patterns in R. chrysoleuca populations.— – Multiple insertions were found in the SSU rDNA of R. chrysoleuca. This phenomenon first was reported in lichenized fungi in the Cladonia chlorophaea complex (DePriest and Been 1992Go, DePriest 1993Go). Group I introns are considered to be gained by horizontal transfer between different organisms (e.g. Gargas et al 1995Go, Nishida and Sugiyama 1995Go, Hibbett 1996Go, Adams et al 1998Go, Perotto et al 2000Go), and spliceosomal introns also are considered to be of recent origin (e.g. Bhattacharya et al 2000Go). Consequently group I introns and spliceosomal introns generally are excluded in systematic analysis. However some studies revealed that group I introns had been immobile for the most part and were transmitted vertically during diversification at lower taxonomic levels (e.g. Nikoh and Fukatsu 2001Go, Bhattacharya et al 2002Go). Hence they have been used in some systematic analysis (e.g. Thell 1999Go, Myllys et al 2001Go). Because no insertions occur in the same position in all samples of R. chrysoleuca, we did not apply them in the phylogenetic analysis but simply compared them. The result shows that the insertions in the same position are high homologous.

There are five positions of insertion in the partial SSU rDNA (802–1800 in FIG. 3Go) of R. chrysoleuca, so there should be 25–1=31 type of insertion distribution patterns. However only four patterns were observed. Even if the positions of group I introns, which are generally accepted transposable elements, are only considered (and assuming that samples 5A and 17A had alternative patterns), just 5(3 + 2) of 24–1=15 combinations were found. Of interest, the samples from clades I and II have distinct insertion distribution patterns. It seems that individuals maintain the same insertion distribution pattern within a recombination population and that it is difficult to form new multiform patterns by intron acquisition and loss. This means that the mechanisms of intron acquisition and intron loss created equilibrium in these populations. Because group I introns without ORF would be eliminated rapidly by intron loss (Dujon 1989Go), we think that sexual reproduction is the main mechanism for intron gain. Further intronless loci would acquire introns from their intron-harboring alleles in zygotes, which suggest that the insertions distribution patterns in a recombination population would be the same. Thus the two different insertion distribution patterns found in gray and green individuals of R. chrysoleuca on Mount Wuling respectively suggests that no gene exchange occurs between the two populations and the isolation mechanism has established.

The phylogenetic analysis of R. chrysoleuca.— – Individuals from Mount Wuling with a gray thallus (morphotype I) assemble in a monophyletic clade (clade I) that is supported by a 97% bootstrap value and those with a green thallus assemble in a second clade (clade II) with 97% support (FIG. 2Go). Neither clade is well resolved, which is the result of an absence of informative characters within each. So it means that the variation in the ITS sequence is insufficient to discriminate the inner structures of R. chrysoleuca populations on Mount Wuling or on the other hand that reproduction in these small populations is clonal. If it implied exclusive clonal reproduction, one would expect that the genotype of all individuals would correspond to their chemotypes. However our data did not support such a result (FIG. 2Go). One genotype could share more than one chemotype, which means that the individuals with the same genotype are not necessarily clonal descendants and recombination must occur at some time. In fact the species of Rhizoplaca, including R. chrysoleuca, produce abundant sexual structures (ascomata) and sexual spores seem to provide the only means of dispersal because no symbiotic vegetative propagules are observed. Therefore recombination should occur frequently in Rhizoplaca.

The two groups of chemotypes do not cross the boundary between clade I and clade II on Mount Wuling, which demonstrates the evidence of genetic isolation. If only samples collected from Mount Wuling are considered, combining the data from the ITS sequences, chemotype analyses and the insertion distribution patterns indicates that there is no gene flow between clade I and clade II and the recombination occur just within each clade. It is surprising because the individuals were collected in a small geographic area so that the fungi certainly had the opportunity to interbreed (FIG. 1B, CGo). This indicates that sexual isolating mechanisms have been established between gray and green thalline groups on Mount Wuling, which implies that the specimens of R. chrysoleuca from Mount Wuling should be divided into two phylogenetic species.

However, when all the collections of R. chrysoleuca are taken into consideration, the conclusion is quite different. R. chrysoleuca is supported as a monophyletic group by a 100% bootstrap value in the ITS phylogenetic analysis. If we accepted that clade I and clade II represented different phylogenetic species, we would confront much more difficult problems in dealing with individuals outside those two clades. In fact R. chrysoleuca is a taxon with a broad distribution; therefore there are many separate geographical and ecological populations in the world. No evidence indicates that the potential for gene exchange among these different R. chrysoleuca populations has been lost. Only when populations differ enough to inhibit any gene exchange at all do we commonly view them as separate species. However a series of intermediate forms between clades I clade II can be found in other areas (FIG. 2Go) and some samples have close genetic relation with these two clades (for example, 17A and clade I, or 5A and clade II) and may be able to mate them. Furthermore the close genetic relationship among individuals of R. chrysoleuca is also reflected by the conservative SSU rDNA sequences.

Therefore we think that R. chrysoleuca as a monophyletic taxon is still accepted and that the graythalline and green-thalline individuals from Mount Wuling that represented different evolutional lineages could be treated as to represent two varieties.

Sexual isolation and speciation.— – In general there are two main modes of speciation: allopatric and sympatric. And speciation should occur in this sequence: (i) geographical or ecological isolation between populations; (ii) increased genetic divergence; (iii) selection for increased reproductive isolation; (iv) speciation completed. Our result indicates that on Mount Wuling R. chrysoleuca consist of two distinct populations based on morphological, chemical and genetic characters. Although reproductive isolation between these two populations exists we think speciation is still in the third stage and does not complete because a series of intermediate forms (such as 5A, 7A–10A, 12A, 17A, G8 and U302) that have close genetic relation with these two populations can be found in other areas and may be able to mate either of them. Thus it may provide a good example of speciation in lichenized fungi. We speculate that the population with a gray thallus eventually will evolve into a separate species. It has been shown that sexual reproduction predominates in lichen communities growing directly on rocks and at the frontiers of terrestrial life in polar regions. As a pioneer species R. chrysoleuca thus may provide a model for the speciation of saxicolous organisms and for the evolution of breeding systems in extreme environments.

The variation in ITS rDNA sequence patterns and secondary chemistry among individuals in R. chrysoleuca within a small area would suggest that most species mats do not represent a single genetic clone but multiple genetic clones or that these lichen thalli may have been composed of more than one sexually active fungal genotypes. Our results indicate that the relationships between genotypes and chemotypes are complicated and it should be avoided and that single individuals or few individuals from the same site should be used to represent the population or whole species in systematics study.

The interior branches within R. chrysoleuca and R. subdiscrepans are longer than the interspecific branches in the core group and those within R. peltata (FIG. 2Go). In fact, Arup and Grube had noticed the unusual longer branch of R. chrysoleuca in their study (2000); however, they did not make any interpretation. Such a large variety of R. chrysoleuca is the most likely to be an adoption to the multiple geographical and ecological zones. Another interesting phenomena is that the ITS sequences within R. chrysoleuca from the same site are sometimes more variable than those from different sites. It seems that it was caused by the movement of different R. chrysoleuca populations to the same geographical area. However the further conclusion can not be drawn until an extensive investigation of R. chrysoleuca on a large scale is undertaken.


    ACKNOWLEDGMENTS
 
We thank Mr Hai-ying Wang, Mr Jun Yang and Miss Xin-Li Wei for help with the field collections. The authors thank Dr J. Elix (Chemistry Department, Australian National University) for reviewing the manuscript. We would like to thank two anonymous reviewers for their helpful comments. This work was supported by grants from National Natural Science Foundation of China (grants 30230020, 30470007) and Chinese Academy of Sciences (Grant kscxz-sw-101C).


    FOOTNOTES
 
Accepted for publication December 7, 2005.

1 Corresponding author. E-mail: shouyug{at}yahoo.com


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Adams KL, Clements MJ, Vaughn JC. 1998. The Peperomia mitochondrial coxI group I intron: timing of horizontal transfer and subsequent evolution of the intron. J Mol Evol 46:689–696.[CrossRef][Medline]

Arup U, Grube M. 1998. Molecular systematics of Lecanora subgeneus Placodium. Lichenologist 30:415–425.

———, ———. 2000. Is Rhizoplaca (Lecanorales, lichenized Ascomycota) a monophyletic genus? Can J Bot 78:318–327.[CrossRef]

Bhattacharya D, Lutzoni F, Reed V, Simon D, Nason J, Fernandez F. 2000. Widespread occurrence of spliceosomal introns in the rDNA genes of Ascomycetes. Mol Biol Evol 17:1971–1984.[Abstract/Free Full Text]

———, Friedl T, Helms G. 2002. Vertical evolution and intragenetic spread of lichen-fungal group I introns. J Mol Evol 55:74–84.[CrossRef][Medline]

Carter DA, Taylor JW, Dechairo B, Burt A, Koenig GL, White TJ. 2001. Amplified single-nucleotide polymorphisms and a (GA)n microsatellite marker reveal genetic differentiation between populations of Histoplasma capsulatum from the Americas. Fungal Genet Biol 34:37–48.[CrossRef][Medline]

Cech TR. 1988. Conserved sequences and structures of group I introns: building an active site for RNA catalysis—a review. Gene 73:258–271.

Choisy M. 1929. Genres nouveaux pour la lichénologie dans le groupe des Lécanoracées. Bull Soc Bot Fr 76:521–527.

Cruse M, Telerant R, Gallagher T, Lee T, Taylor JW. 2002. Cryptic species in Stachybotrys chartarum. Mycologia 94:814–822.[Abstract/Free Full Text]

Cubero OF, Brideg PD, Crespo A. 2000. Terminal-sequence conservation identifies spliceosomal introns in Ascomycete 18S RNA genes. Mol Biol Evol 17:751–756.[Abstract/Free Full Text]

Culberson CF. 1972. Improved conditions and new data for the identification of lichen products by a standardized thin-layer chromatographic method. J Chromatogr 72:113–125.[CrossRef][Medline]

———, Armaleo D. 1992. Induction of complete secondary-product pathway in a cultured lichen fungus. Experimental Mycology 16:52–63.[CrossRef]

———, Culberson WL, Johnson A. 1988. Gene flow in lichens. Amer J Bot 75:1135–1139.[CrossRef]

Culberson WL. 1969. The use of chemistry in the systematics of lichens. Taxon 18:152–166.

———. 1986. Chemistry and sibling speciation in the lichen-forming fungi: ecological and biological considerations. The Bryologist 89:123–131.[CrossRef]

DePriest PT. 1992. Molecular-genetic analysis of ribosomal DNA polymorphism in the Cladonia chlorophaea complex. [Ph. D. dissertation]. Durham: Duck University.

———. 1993. Small subunit nrDNA variation in a population of lichen fungi due to optional group I introns. Gene 134:67–74.[CrossRef][Medline]

———. 1994. Variation in the Cladonia chlorophaea complex II: Ribosomal DNA variation in a Southern Appalachian population. The Bryologist 89:123–131.[CrossRef]

———, Been MD. 1992. Numberous group I introns with variable distributions in the rDNA of a lichen fungus. J Mol Biol 228:315–321.[CrossRef][Medline]

Dujon B. 1989. Group I introns as mobile genetic elements: facts and mechanistic speculations—a review. Gene 82:91–114.[CrossRef][Medline]

Gargas A, Taylor JW. 1992. Polymerase chain reaction (PCR) primers for amplifying and sequencing nuclear 18S rDNA from lichenized fungi. Mycologia 84:589–592.[CrossRef]

———, DePriest PT, Taylor JW. 1995. Positions of multiple insertions in SSU rDNA of lichen-forming fungi. Mol Biol Evol 12:208–218.[Abstract]

Goward T, Goffinet B. 2000. Peltigera chionophila, a new lichen (Ascomycetes) from the western cordillera of North American. The Bryologist 103:493–498.[CrossRef]

Groner U, LaGreca S. 1997. The ‘Mediterranean’ Ramalina panizzei north of the alps: morphological, chemical and rDNA sequence data. Lichenologist 29:582–586.

Hawksworth DL. 1976. Lichen chemotaxonomy. In: Brown DH, Hawksworth DL, Bailey RH, eds. Lichenology: progress and problems. London: Academic Press. p 139–184.

Hibbett DS. 1996. Phylogenetic evidence for horizontal transmission of group I intron in the nuclear ribosomal DNA of mushroom-forming fungi. Mol Biol Evol 13:903–917.[Abstract]

Högberg N, Kroken S, Thor G, Taylor JW. 2002. Reproductive mode and genetic variation suggest a North American origin of European Letharia vulpina. Mol Ecol 11:1191–1196.[CrossRef][Medline]

Kieft TL. 1988. Ice nucleation activity in lichens. Appl Environ Microbiol 54:1678–1681.[Abstract/Free Full Text]

———, Ruscetti T. 1990. Characterization of biological ice nuclei from a lichen. J Bacteriol 172:3519–3523.[Abstract/Free Full Text]

Kirk PM, Cannon PF, David JC, Stalpers JA. 2001. Dictionary of the Fungi. Wallingford: CABI publishing.

Koch MA, Dobes C, Mitchell-Olds T. 2003. Multiple hybrid formation in natural populations: concerted evolution of the internal transcribed spacer of nuclear ribosomal DNA (ITS) in North American Arabis divaricarpa (Brassicaceae). Mol Biol Evol 20:338–350.[Abstract/Free Full Text]

Koufopanou V, Burt A, Szaro T, Taylor JW. 2001. Gene genealogies, cryptic species, and molecular evolution in the human pathogen Coccidioides immitis and relatives (Ascomycota, Onygenales). Mol Biol Evol 18:1246–1258.[Abstract/Free Full Text]

Kroken S, Taylor JW. 2001. Outcrossing and recombination in the lichenized fungus Letharia. Fung Genet Biol 34:83–92.

Kumar S, Tamura K, Jakobsen IB, Nei M. 2001. MEGA2: Molecular evolutionary genetics analysis software. Tempe, Arizona, USA: Arizona State University.

Leuckert C, Schmitz KE, Feige GB. 1977. Zur chemotaxonomie der eurasischen arten der Flechtengattung Rhizoplaca. Nova Hedwigia 28:71–129.

Mehmann B, Brunner I, Braus GH. 1994. Nucleotide sequence variation of chitin synthase genes among ectomycorrhizal fungi and its potential use in taxonomy. Appl Environ Microb 60:3105–3111.[Abstract/Free Full Text]

Murtagh GJ, Dyer PS, Crittenden PD. 2000. Sex and the single lichen. Nature 404:564–564.[CrossRef][Medline]

Myllys L, Lohtander K, Tehler A. 2001. ß-tubulin, ITS and group I intron sequences challenge the species pair concept in Physcia aipolia and P. caesia. Mycologia 93:335–343.[CrossRef]

———, Stenroos S, Thell A. 2002. New genes for phylogenetic studies of lichenized fungi: glyceraldehyde-3-phosphate dehydrogenase and beta-tubulin genes. Lichenologist 34:237–246.[CrossRef]

Narui T, Culberson CF, Shibata S. 1996. A contribution to the chemistry of the lichen family Umbilicaricaeae (Ascomycotina). The Bryologist 99:199–211.[CrossRef]

Nikoh N, Fukatsu T. 2001. Evolutionary dynamics of multiple group I introns in nuclear ribosomal RNA genes of endoparasitic fungi of the genus Cordyceps. Mol Biol Evol 18:1631–1642.[Abstract/Free Full Text]

Nishida H, Sugiyama J. 1995. A common group I intron between a plant parasitic fungus and its host. Mol Biol Evol 12:883–886.[Abstract]

Perotto S, Nepote-Fus P, Saletta L, Bandi C, Peter J, Young W. 2000. A diverse population of introns in the nuclear ribosomal genes of ericoid mycorrhizal fungi includes elements with sequence similarity to endonuclease-coding genes. Mol Biol Evol 17:44–59.[Abstract/Free Full Text]

Presa P, Pardo GG, Martinez P, Bernatchez L. 2002. Phylogeographic congruence between mtDNA and rDNA ITS markers in brown trout. Mol Biol Evol 19:2161–2175.[Abstract/Free Full Text]

Rogers SO, Bendich AJ. 1988. Extraction of DNA from plant tissues. Plant Molecular Biology Manual A6:1–10. Kluwer Academic Publishers, Dordrecht.

Ryan BD, Nash III TH. 1997. Systematics of Lecanora subgenus Placodium (lichenizied Ascomycotina) in North America: an overview with keys. Nova Hedwigia 64:111–127.

Thell A. 1999. Group I intron versus ITS sequences in phylogeny of cetrariod lichens. Lichenologist 31:441–449.

Walser J-C, Holderegger R, Gugerli F, Hoebee SE, Scheidegger C. 2005. Microsatellites reveal regional population differentiation and isolation in Lobaria pulmonaria, an epiphytic lichen. Mol Ecol 14:457–467.[CrossRef][Medline]

Wei J. 1984. A preliminary study of lichen genus Rhizoplaca from China. Acta Mycol Sin 3:207–213.

White YJ, Bruns T, Lee S, Taylor JW. 1990. Amplication and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, eds. PCR protocols: a guide to methods and applications. San Diego, California: Academic Press. p 315–322.





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Services
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zhou, Q.-M.
Right arrow Articles by Wei, J.-C.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Zhou, Q.-M.
Right arrow Articles by Wei, J.-C.
Agricola
Right arrow Articles by Zhou, Q.-M.
Right arrow Articles by Wei, J.-C.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS