| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Institute of Plant Sciences, Holteigasse 6, Karl-Franzens-University Graz, A-8010 Graz, Austria
| ABSTRACT |
|---|
|
|
|---|
A phylogenetic analysis of the Lecanora rupicola group based on combined nITS rDNA and ß-tubulin sequences and a combined dataset of ITS, ß-tubulin and partial sequences of polyketide synthase genes reveals a previously unrecognized species, which here is introduced under the name Lecanora bicinctoidea. The new species is a sister group of the L. swartzii complex (including L. swartzii and L. lojkaeana), which is characterized by eucorticate ascomata, and a morphological diversity that includes also a dwarf-fruticose lineage. The preferential occurrence on vertical to overhanging siliceous rocks corresponds more closely to L. swartzii. A detailed investigation of phenotypic characters reveals that the new species differs from the superficially similar morphospecies L. bicincta in several ways, such as a thallus of comparatively small areoles and broadly sessile ascomata and the development of an amphithecial cortex devoid of algal remnants (i.e. an eucortex). L. bicinctoidea contains methyl 3
-hydroxy-4-O-demethylbarbatate, a chemical compound not known from other members of the L. rupicola group. We also discuss the importance of eucortex formation as one of several factors that are involved in the evolution of substrate-detached thallus structures.
Key words: ß-tubulin, ITS, KS-domain, Lecanora rupicola group, lichenized ascomycota, phylogeny, secondary metabolites
| INTRODUCTION |
|---|
|
|
|---|
Grube et al (2004)
also confirmed that the Lecanora rupicola group includes not only the saxicolous taxa monographed by Leuckert and Poelt (1989)
but also corticolous taxa of the Lecanora carpinea group. Moreover an additional sterile species, Lecanora rouxii, previously treated under a different name in the genus Lepraria, was added. The study however also questioned previous concepts of species recognition in the genus and showed that the polymorphic species such as L. rupicola and L. bicincta are not supported by ITS data alone.
L. rupicola and L. bicincta are distinguished by the presence of an apically pigmented parathecium that is formed in the latter. Both morphospecies are otherwise morphologically diverse, and therefore Leuckert and Poelt (1989)
attempted an infraspecific classification primarily using characters of secondary chemistry. Only in cases of clearly deviating morphotypes, subspecific entities were based on morphology (e.g. L. swartzii ssp. nuorensis). Similar chemotypes were found in both morphospecies and numerous specimens were labeled by Josef Poelt as L. bicincta s.l. and L. rupicola s.l. in the herbarium at the Institute of Plant Sciences in Graz, where the major part of his substantial collections of this group is curated. One of these samples deviated from typical L. bicincta by somewhat smaller thallus areoles with narrowly sinuous margins, and we were interested whether this morphotype could be shown to be distinct by molecular data. Using the information on the label we located this morphotype from the original collecting locality and subjected this material to a molecular analysis. In this study we show that this morphotype is indeed distinct from the unresolved bulk of L. rupicola/bicincta species complex and we describe it as an additional and distinct species in the L. rupicola group.
| MATERIAL AND METHODS |
|---|
|
|
|---|
|
DNA extraction, amplification and sequencing.
For molecular analyses the thalli first were checked for well developed areoles that did not show any externally visible contaminations by other organisms. Total DNA was extracted from individual thalli, according to a modified CTAB method (Cubero et al 1999
) or using the DNeasy Plant Mini Kit (QIAGEN, Vienna) following the manufacturers instructions, with the exception that the DNA was eluted in sterile water. Primers used for PCR of the nuclear ribosomal ITS region were ITS1F (Gardes and Bruns 1993
) and ITS4 (White et al 1990
), those for ß-tubulin were bt2a and bt2b (Glass and Donaldson 1995
) and those for the KS-domain of polyketide synthase genes were LC1 and LC2C (Bingle et al 1999
). Fifty µL PCR mix (10 mM Tris pH8.3, 50 mM KCl/1.5 mM MgCl2/50 µg gelatine) contained 1.25 units of Taq DNA Polymerase (Amersham Pharmacia Biotech Inc.), 0.2 mM of each of the four dNTPs, 0.5 µM of each primer and ca. 1050 ng genomic DNA. Products were cleaned with QIAquick PCR Purification Kit (QIAGEN, Vienna). Both complementary strands were sequenced with the BigDyeTerminator Cycle Sequencing Ready Reaction Kit (Applera, Vienna) according to the manufacturers instructions. Sequences were run on an ABI 310 automated sequencer (Applera, Vienna) and raw data were assembled with AutoAssembler (Applera, Vienna). To assess the relationship within the Lecanora rupicola group, we constructed an alignment that included selected species of the group by using the Clustal algorithm as implemented in BioEdit (http://jwbrown.mbio.ncsu.edu/BioEdit/bioedit.html). The alignment was improved manually and group I intron sequences at the end of the SSU rDNA gene (at position 1506, Gargas et al 1995
) as well as ambiguously aligned positions were removed. The sequences are submitted to EMBL/GenBank, and a list of sequenced taxa and sequences retrieved from GenBank for comparisons together with their GenBank accession numbers is presented (TABLE I
).
Phylogenetic analysis.
The phylogenetic hypothesis was constructed with a Bayesian approach as implemented in the program MrBayes v3.1. (Huelsenbeck and Ronquist 2001
). Because the efficiency of a Bayesian analysis is sensitive to the model parameters we obtained the appropriate model of nucleotide substitution for each gene with the MrModeltest 2.2. (by J.A.A. Nylander, available at http://www.ebc.uu.se/systzoo/staff/nylander.html). According to the Akaike information criterion (AIC), the symmetrical model SYM +
(Zharkihk 1994) with estimation of site specific rates and assuming a discrete gamma distribution within a class was suggested as the optimal model for ß-tubulin and the KS-domain of putative polyketide synthase (PKS) genes and the general time reversible substitution model GTR +
(Rodriguez et al 1990
) assuming a discrete gamma distribution for ITS. To improve the resolution we combined ITS and ß-tubulin for a larger dataset, and ITS, ß-tubulin and PKS for selected taxa. Before combining, separately analysed datasets revealed no conflicts in the tree topologies. We also performed a partition homogeneity test as implemented in PAUP*4.08b (Swofford 2002
), which showed that the partitions are not incongruent (P value = 0.01). The appropriate models (see above) were set for each partition in the combined dataset. In both phylogenies the Markov chain Monte Carlo (MCMC) analysis was run 3 000 000 generations, with four chains starting from a random tree, and using the temperature of 0.1. Every 100th tree was sampled, while the first 10 000 generations were discarded as burn-in. The parsimony analyses were performed with PAUP*4.08b (Swafford 2002), using 1000 replicates. Bootstrap values higher than 75% are indicated as numbers above branches (FIGS. 1
, 2
). As outgroup we chose Lecanora macrocyclos and L. allophana (FIG. 1
) and L. carpinea (FIG. 2
), respectively.
|
|
| RESULTS |
|---|
|
|
|---|
The majority rule consensus tree of the combined ITS/ß-tubulin dataset (FIG. 1
) was calculated from 580 002 trees, the average log-likelihood of the tree sample is 4264.91 (arithmetic mean). The likelihood parameters in the sample had these average values (variance) for the ITS partition: rate matrix r(GT) = 0.083 (0), r(CT) = 0.393 (± 0.001), r(CG) = 0.113 (0), r(AT) = 0.128 (0), r(AG) = 0.188 (0), r(AC) = 0.094 (± 0.003), gamma shape parameter
= 0.413 (± 0.005); and these for the ß-tubulin partition: rate matrix r(GT) = 0.092 (0), r(CT) = 0.397 (± 0.001), r(CG) = 0.084 (0), r(AT) = 0.164 (0), r(AG) = 0.180 (0), r(AC) = 0.092 (0), and the gamma shape parameter
= 2.431 (± 0.030). Estimated base frequencies of the GTR model for the ITS dataset are:
(A) = 0.221 (0),
(C) = 0.272 (0),
(G) = 0.260 (0),
(T) = 0.246 (0).
The majority rule consensus tree of the combined ITS, ß-tubulin and KS dataset (FIG. 2
) was calculated from 580 002 trees, the average log-likelihood of the tree sample is 4023.08 (arithmetic mean). The likelihood parameters in the sample had these average values (variance) for the ITS partition: rate matrix r(GT) = 0.067 (0), r(CT) = 0.426 (± 0.001), r(CG) = 0.127 (0), r(AT) = 0.050 (0), r(AG) = 0.214 (± 0,001), r(AC) = 0.117 (0), gamma shape parameter
= 1.280 (± 0.004); and these for the ß-tubulin partition: rate matrix r(GT) = 0.075 (± 0.001), r(CT) = 0.412 (± 0.004), r(CG) = 0.163 (± 0.002), r(AT) = 0.095 (± 0.001), r(AG) = 0.201 (± 0.002), r(AC) = 0.053 (± 0.001), gamma shape parameter
= 1.368 (± 0.009); and these for the KS partition: rate matrix r(GT) = 0.048 (0), r(CT) = 0.513 (± 0.003), r(CG) = 0.068 (0), r(AT) = 0.019 (0), r(AG) = 0.277 (± 0.002), r(AC) = 0.075 (0) and the gamma shape parameter
= 0.240 (± 0.086). Estimated base frequencies of the GTR model for the ITS dataset are:
(A) = 0.223 (0),
(C) = 0.262 (0),
(G) = 0.262 (0),
(T) = 0.251 (0).
The combined Bayesian analysis of ITS and ß-tubulin (FIG. 1
) resulted in four distinct, well supported clades (posterior probability 100%) within the Lecanora rupicola group. Clade A comprises the two closely related species L. rupicola and L. bicincta, yet the two species are not resolved. Representing the Lecanora swartzii subgroup, Clade B includes the new species L. bicinctoidea and its sister Clade C comprises L. swartzii as well as the sorediate L. lojkaeana. The L. carpinea subgroup with the corticolous L. carpinea and L. subcarpinea (Clade D) is the sister group of the saxicolous species. The phylogeny of three combined genes (ITS, ß-tubulin, KS-domain of putative polyketide synthase genes) of selected taxa (FIG. 2
), with L. carpinea as outgroup, shows a similar topology.
In addition to the known secondary compounds of the members of the group (i.e. atranorin, arthothelin, eugenitol, psoromic acid, sordidone and thiophanic acid) the chemical analysis of L. bicinctoidea revealed a new compound in the group, the depside methyl 3
-hydroxy-4-O-demethylbarbatate.
| TAXONOMY |
|---|
|
|
|---|
|
-hydroxy-4-O-demethylbarbatatum continentibus. Thallus crustose, cream-white to white, areoles small 0.40.60.8(1.0) mm diam, flat to slightly convex, areole edges sinuous, surface ± smooth; soredia absent; prothallus distinct, black;
Ascomata 0.7 1.0 1.3(1.7) mm diam, round, broadly sessile, arising singly to juxtaposed, thalline exciple entire and slightly sinuous; disks blue-gray to blue, flat to convex, densely gray to white pruinose; epithecium dark olivaceous brown, interspersed with crystals dissolving in K; hymenium hyaline, 60 75 90(100) µm tall; paraphyses 2 µm thick, sparsely branched, apices not or slightly thickened, 23 µm. Asci 8-spored, ca. 15 µm broad, Lecanora-type; ascospores hyaline, simple, ellipsoid, wall thick, 10.211.613(15) ± 5.76.67.5(8); conidiomata with dark brown to black ostiolum, conidia filiform, (18)25(29) ± 1 µm.
Chemistry.
Thallus K ± pale yellow, C , apothecial disk Pd, K, C+ yellow to orange. Arthothelin, atranorin, chloroatranorin, eugenitol, isoarthothelin, methyl 3
-hydroxy-4-O-demethylbarbatate, psoromic acid (facultative), sordidone, thiophanic acid (facultative).
Etymology. The epithet bicinctoidea is based on the similarity of the new species and Lecanora bicincta; both share a more or less distinct black margin in the upper part of the ascomata between the hymenium and the thalline excipulum.
Ecology and distribution. Lecanora bicinctoidea grows on vertical to overhanging siliceous rocks. So far the species is known from various mountain ranges in Styria, Austria, (Fischbacher Alpen, Gurktaler Alpen, Seckauer Tauern, Seetaler Alpen, Steirisches Randgebirge) and Italy (Nimis personal communication).
HOLOTYPE: AUSTRIA. STYRIA:. Steirisches Randgebirge, Koralpe, Handalpe, hiking trail from the Weinebene to the Moserkogel, NW Handalpe, ca. 1800 m, siliceous rocks, Öfen, 28-IX-2002; J. Blaha (GZU).
ADDITIONAL SPECIMENS EXAMINED: AUSTRIA. STYRIA:. Seetaler Alpen, Zirbitzkogel, 200 m SW of Großer Winterleitensee; J. Blaha B277 (GZU). Seetaler Alpen, Zirbitzkogel, steep rocky faces on way to the Ochsenboden; J. Blaha B282, B284 (GZU). Seckauer Tauern, Gaalgraben, Krugspitze, 50 m N of the Krugsee; J. Blaha B321 (GZU). Fischbacher Alpen, Stuhleck; M. Grube EB26 (GZU). Steirisches Randgebirge, Koralpe, Handalpe; M. Grube EB 31 (GZU). Steirisches Randgebirge, Koralpe, Handalpe; E. Baloch EB111 (GZU). Gurktaler Alpen, Frauenalpe, J. Blaha B451 (GZU). Seckauer Alpen, Hochreichart, B. Emmerer, H. Komposch, J. Blaha B437, B438 (GZU).
| DISCUSSION |
|---|
|
|
|---|
Lecanora bicinctoidea (FIGS. 3, 4
) is clearly distinguished from the poorly resolved Lecanora bicincta/rupicola complex by molecular, morphological and chemical data. L. bicinctoidea contains, in addition to sordidone, eugenitol and atranorin, the depside methyl 3
-hydroxy-4-O-demethylbarbatate, which is not known from other members of the L. rupicola group. This ß-orcinol derivate is biosynthetically closely related to atranorin, methyl barbatate and methyl 4-O-demethylbarbatate. The reduction of atranorin at
-C3 results in methyl 3
-hydroxy-4-O-demethylbarbatate. It is not common in lichenized ascomycetes and first was documented in Oropogon loxensis (Culberson and Culberson 1978
). The compound was detected both by HP-TLC (according to Arup et al 1993
, Rf-classes 5/5/5) and corroborated with HPLC by J. Elix. The ecological requirements of Lecanora bicinctoidea are vertical to overhanging siliceous rocks, where it consistently associates with Trebouxia simplex (as assessed by algal ITS sequencing, data not shown). This humid and shaded environment corresponds to that of L. swartzii (Leuckert and Poelt 1989
), which usually co-occurs at sampling sites in several mountain ranges in Styria investigated so far.
In this study we also included the ketoacyl synthase (KS) domain of polyketide synthase genes (fungal type I PKS) as an additional locus to confirm the distinctiveness of the new species. Polyketide synthase genes represent a gene family of considerable functional importance in fungi because the translated proteins are essential in the production of secondary metabolites. Grube and Blaha (2003)
suggested that paralogy of PKS genes can be detected by phylogenetic incongruency with rDNA data, and Schmitt et al (2005)
found that numerous paralogs for nonreducing polyketide synthases exist in lichens. Within Lecanora three such paralogs have been discovered (Blaha et al unpublished data). We therefore have taken care to include only orthologous KS-sequences from our samples in this analysis. The PKS paralog included here belongs to a subgroup in Clade V of Schmitt et al (2005)
, in which all paralogs detected directly from the Lecanora rupicola group with the primers LC1 and LC2C (Bingle et al 1999
) form a monophyletic group. In addition the presence of an mRNA intron at the same position further corroborates the orthology of this PKS gene. Whereas the functional assignment of the analysed PKS genes to a particular biosynthetic pathway is difficult and would require extensive experimental studies, we are convinced that orthologous PKS genes represent valuable markers for phylogenetic studies within genera.
Lecanora swartzii ssp. caulescens was discussed before as an example for the evolution of fruticose growth form among crustose relatives. Poelt (1989)
found that an eucortex is developed at the upper surface in this dwarf-fruticose lichen. He hypothesized that the early formation of an eucortex is usually connected to the ascomatal ontogeny in this subspecies. Subsequently Poelt somehow generalized this observation and suggested that the proleptic formation of ascomatal eucortices and the delay of ascoma formation could be a principal phenomenon in the evolution of foliose thalli; in this context the eucortex was functionally interpreted as a requirement to stabilize the thallus structure ("exoskeleton", Poelt 1991
). However this hypothesis was not developed further (e.g. by comparative ontogenetic studies).
We assume that the exclusion of algal remnants in otherwise gelatinized fungal plectenchyma could promote the evolution of substrate-detached thallus types. Layers in heteromerous thalli may differ in the capacity to increase in size with water uptake. Layers, which consist of a matrix of extracellular polysaccharides (e.g. hymenia, eucortices), usually swell more pronouncedly in water than layers without such a matrix (Blaha unpublished data). As a consequence the repeated swelling and shrinking of thalli under natural conditions of fluctuating hydration likely acts as a mechanic force in a physiologically active and growing thallus, especially when the thallus is fixed with the lower surface to a rigid substrate (such as rocks). Under these conditions strictly crustose thalli with an eucortex might be regarded as exceptional. It in fact is interesting that growth-form transitions observed in other groups of rock-inhabiting Lecanora also correlate with increasingly well developed eucortices, e.g. in the Lecanora polytropa group (crustose L. polytropa vs. lobate L. concolor, Arup and Grube 1998
) and in the L. novomexicana group (lobate L. novomexicana vs. foliose Rhizoplaca melanophthalma, Arup and Grube 2000
). Similar hygroscopic forces also could influence apothecial functions. It might be hypothesized that the eucorticate margin could compensate for the hygroscopic movements of hymenial parts. Preliminary observations in L. swartzii suggest that the eucorticate amphithecial cortex acts like a clamp that keeps hydrated ascomatal sections in shape, whereas removal of the eucortex causes the hymenium to disintegrate from subtending structures (Blaha unpublished data).
Our data show that the Lecanora swartzii group including L. bicinctoidea contains only species that develop a eucorticate amphithecium, whereas the related strictly crustose L. rupicola/bicincta complex is characterized by ascomatal margins with a phenocortex (see FIGS. 58
). L. bicinctoidea is thus a basal crustose lineage in a species group that exhibits an interesting morphological variation, from broadly sessile to distinctly elevated and basally constricted ascomata, and from crustose to dwarf-fruticose thalli.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 Corresponding author. Current address: Department of Biology, University of Crete, P.O. Box 2208, 71409 Heraklion, Greece. E-mail: julia.blaha{at}uni-graz.at
| LITERATURE CITED |
|---|
|
|
|---|
, . 2000. Is Rhizoplaca (Lecanorales, lichenized Ascomycota) a monophyletic genus? Can J Bot 78(3):318327.[CrossRef]
, Ekman S, Lindblom L, Mattsson JE. 1993. High performance thin layer chromatography (HPTLC), an improved technique for screening lichen substances. Lichenologist 25:6171.
Bingle LEH, Simpson TJ, Lazarus CM. 1999. Ketosynthase domain probes identify two subclasses of fungal polyketide synthase genes. Fung Gen Biol 26:209223.[CrossRef]
Cubero OF, Crespo A, Fatehi J, Bridge PD. 1999. DNA Extraction and PCR amplification method suitable for fresh, herbarium stored and lichenized fungi. Plant Syst Evol 216:243249.[CrossRef]
Culberson CF, Culberson WL. 1978. ß-Orcinol derivatives in lichens: biogenetic evidence from Oropogon loxensis. Exp Mycol 2:245257.[CrossRef]
Gardes M, Bruns TD. 1993. ITS primers with enhanced specificity for basidiomycetesapplication for the identification of mycorrhizae and rusts. Mol Ecol 2: 113118.[Medline]
Gargas A, DePriest PT, Taylor JW. 1995. Positions of multiple insertions in SSU rDNA of lichen-forming fungi. Mol Biol Evol 12:208218.[Abstract]
Glass NL, Donaldson GC. 1995. Development of primer sets designed for use with the PCR to amplify conserved genes from filamentous ascomycetes. Appl Environ Microbiol 61:13231330.
Grube M, Blaha J. 2003. On the phylogeny of some polyketide synthase genes in the lichenized genus Lecanora. Mycol Res 107:14191426.[CrossRef][Medline]
, Baloch E, Arup U. 2004. A phylogenetic study of the Lecanora rupicola group (Lecanoraceae, Ascomycota). Mycol Res 108:506514.[CrossRef][Medline]
Huelsenbeck JP, Ronquist F. 2001. MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754755.
Kirk PM, Cannon PF, David JC, Stalpers JA, eds. 2001. Ainsworth & Bisbys dictionary of the Fungi. 9th ed. Wallingford: CABI Publishing. 655 p.
Leuckert C, Poelt J. 1989. Studien über die Lecanora rupicola-Gruppe in Europa (Lecanoraceae). Nov Hedwig 49:121167.
Poelt J. 1989. Die Entstehung einer Stauchflechte aus einem Formenkreis krustiger Verwandter. The origin of a fruticose lichen from a related group of crustose taxa. Flora 183:6572.
. 1991. Homologies and analogies in the evolution of lichens. In: Hawksworth DL, ed. Frontiers in Mycology. Wallingford: CAB International. p 8597.
Rodriguez F, Oliver JL, Marin A, Medina JR. 1990. The general stochastic model of nucleotide substitution. J Theor Biol 142:485501.[Medline]
Schmitt I, Martin MP, Kautz S, Lumbsch HT. 2005. Diversity of non-reducing polyketide synthase genes in the Pertusariales (lichenized Ascomycota): A phylogenetic perspective. Phytochemistry 66:1241253.[CrossRef][Medline]
Swofford DL. 2002. PAUP*: phylogenetic analysis using parsimony (*and other methods). Version 4. Sunderland, Massachusetts: Sinauer Associates.
White TJ, Bruns TD, Lee S, Taylor J. 1990. Amplification and direct sequencing of fungal ribosomal DNA genes for phylogenies. In: Innis MA, Gelfand DH, Sninsky JJ, White TJ, eds. PCR protocols: a guide to methods and applications. San Diego: Academic Press. p 315322.
Zharkikh A. 1994. Estimation of evolutionary distances between nucleotide sequences. J Mol Evol 39:31529.[CrossRef][Medline]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |