The ecology of lichenicolous lichens: a case-study in Italy

This paper, with Italy as a case-study, provides a general overview on the ecology of lichenicolous lichens, i.e. those which start their life-cycle on the thallus of other lichens. It aims at testing whether some ecological factors do exert a positive selective pressure on the lichenicolous lifestyle. The incidence of some biological traits (photobionts, growth-forms and reproductive strategies) in lichenicolous and non-lichenicolous lichens was compared, on a set of 3005 infrageneric taxa potentially occurring in Italy, 189 of which are lichenicolous. Lichenicolous lichens have a much higher incidence of coccoid (non-trentepohlioid) green algae, crustose growth-forms and sexual reproduction. A matrix of the 2762 species with phycobionts and some main ecological descriptors was subjected to ordination. Lichenicolous lichens occupy a well-defined portion of the ecological space, tending to grow on rocks in dry, well-lit habitats where a germinating spore is likely to have a short life-span, at all altitudes. This corroborates the hypothesis that at least some of them are not true "parasites", as they are often called, but gather the photobionts - which have already adapted to local ecological conditions - from their hosts, eventually developing an independent thallus.

Myrmecia, Not Asterochloris, Is the Main Photobiont of Cladonia subturgida (Cladoniaceae, Lecanoromycetes)

This study explores the diversity of photobionts associated with the Mediterranean lichen-forming fungus Cladonia subturgida. For this purpose, we sequenced the whole ITS rDNA region by Sanger using a metabarcoding method for ITS2. A total of 41 specimens from Greece, Italy, France, Portugal, and Spain were studied. Additionally, two specimens from Spain were used to generate four cultures. Our molecular studies showed that the genus Myrmecia is the main photobiont of C. subturgida throughout its geographic distribution. This result contrasts with previous studies, which indicated that the main photobiont for most Cladonia species is Asterochloris. The identity of Myrmecia was also confirmed by ultrastructural studies of photobionts within the lichen thalli and cultures. Photobiont cells showed a parietal chloroplast lacking a pyrenoid, which characterizes the species in this genus. Phylogenetic analyses indicate hidden diversity within this genus. The results of amplicon sequencing showed the presence of multiple ASVs in 58.3% of the specimens studied.

Regionally Varying Habitat Relationships in Lichens: The Concept and Evidence with an Emphasis on North-Temperate Ecosystems

Habitat ecology of lichens (lichen-forming fungi) involves diverse adaptations to stressful environments where lichens use specific habitat conditions. Field observations confirm that such habitat ‘preferences’ can vary significantly across species’ distribution ranges, sometimes revealing abrupt changes over short distances. We critically review and generalize such empirical evidence as broad ecological patterns, link these with the likely physiological mechanisms and evolutionary processes involved, and outline the implications for lichen conservation. Non-replicated correlative studies remain only suggestive because the data are frequently compromised by sampling bias and pervasive random errors; further noise is related to unrecognized cryptic species. Replicated evidence exists for three macroecological patterns: (a) regional limiting factors excluding a species from a part of its microhabitat range in suboptimal areas; (b) microhabitat shifts to buffer regionally adverse macroclimates; (c) substrate suitability changed by the chemical environment, notably air pollution. All these appear to be primarily buffering physiological challenges of the adverse conditions at the macrohabitat scale or, in favorable environments, coping with competition or predation. The roles of plasticity, adaptation, dispersal, and population-level stochasticity remain to be studied. Although lichens can inhabit various novel microhabitats, there is no evidence for a related adaptive change. A precautionary approach to lichen conservation is to maintain long-term structural heterogeneity in lichen habitats, and consider lichen ecotypes as potential evolutionarily significant units and a bet-hedging strategy for addressing the climate change-related challenges to biodiversity.

Phylogeography of Ramalina farinacea (Lichenized Fungi, Ascomycota) in the Mediterranean Basin, Europe, and Macaronesia

Ramalina farinacea is an epiphytic lichen-forming fungus with a broad geographic distribution, especially in the Northern Hemisphere. In the eighties of the last century, it was hypothesized that R. farinacea had originated in the Macaronesian–Mediterranean region, with the Canary Islands as its probable southernmost limit, and thereafter it would have increased its distribution area. In order to explore the phylogeography of this emblematic lichen, we analyzed 120 thalli of R. farinacea collected in 38 localities distributed in temperate and boreal Europe, the Western Mediterranean Basin, and several Macaronesian archipelagos in the Atlantic Ocean. Data from two nuclear markers (nrITS and uid70) of the mycobiont were obtained to calculate genetic diversity indices to infer the phylogenies and haplotype networks and to investigate population structure. In addition, dating analysis was conducted to provide a valuable hypothesis of the timing of the origin and diversification of R. farinacea and its close allies. Our results highlight that phylogenetic species circumscription in the “Ramalina farinacea group” is complex and suggests that incomplete lineage sorting is at the base of conflicting phylogenetic signals. The existence of a high number of haplotypes restricted to the Macaronesian region, together with the diversification of R. farinacea in the Pleistocene, suggests that this species and its closest relatives originated during relatively recent geological times and then expanded its range to higher latitudes. However, our data cannot rule out whether the species originated from the Macaronesian archipelagos exclusively or also from the Mediterranean Basin. In conclusion, the present work provides a valuable biogeographical hypothesis for disentangling the evolution of this epiphytic lichen in space and time.

Trebouxia lynnae sp. nov. (Former Trebouxia sp. TR9): Biology and Biogeography of an Epitome Lichen Symbiotic Microalga

Two microalgal species, Trebouxia jamesii and Trebouxia sp. TR9, were detected as the main photobionts coexisting in the thalli of the lichen Ramalina farinacea. Trebouxia sp. TR9 emerged as a new taxon in lichen symbioses and was successfully isolated and propagated in in vitro culture and thoroughly investigated. Several years of research have confirmed the taxon Trebouxia sp. TR9 to be a model/reference organism for studying mycobiont−photobiont association patterns in lichen symbioses. Trebouxia sp. TR9 is the first symbiotic, lichen-forming microalga for which an exhaustive characterization of cellular ultrastructure, physiological traits, genetic and genomic diversity is available. The cellular ultrastructure was studied by light, electron and confocal microscopy; physiological traits were studied as responses to different abiotic stresses. The genetic diversity was previously analyzed at both the nuclear and organelle levels by using chloroplast, mitochondrial, and nuclear genome data, and a multiplicity of phylogenetic analyses were carried out to study its intraspecific diversity at a biogeographical level and its specificity association patterns with the mycobiont. Here, Trebouxia sp. TR9 is formally described by applying an integrative taxonomic approach and is presented to science as Trebouxia lynnae, in honor of Lynn Margulis, who was the primary modern proponent for the significance of symbiosis in evolution. The complete set of analyses that were carried out for its characterization is provided.

The Role of Photobionts as Drivers of Diversification in an Island Radiation of Lichen-Forming Fungi

Speciation in oceanic islands has attracted the interest of scientists since the 19th century. One of the most striking evolutionary phenomena that can be studied in islands is adaptive radiation, that is, when a lineage gives rise to different species by means of ecological speciation. Some of the best-known examples of adaptive radiation are charismatic organisms like the Darwin finches of the Galapagos and the cichlid fishes of the great African lakes. In these and many other examples, a segregation of the trophic niche has been shown to be an important diversification driver. Radiations are known in other groups of organisms, such as lichen-forming fungi. However, very few studies have investigated their adaptive nature, and none have focused on the trophic niche. In this study, we explore the role of the trophic niche in a putative radiation of endemic species from the Macaronesian Region, the Ramalina decipiens group. The photobiont diversity was studied by Illumina MiSeq sequencing of the ITS2 region of 197 specimens spanning the phylogenetic breadth and geographic range of the group. A total of 66 amplicon sequence variants belonging to the four main clades of the algal genus Trebouxia were found. Approximately half of the examined thalli showed algal coexistence, but in most of them, a single main photobiont amounted to more than 90% of the reads. However, there were no significant differences in photobiont identity and in the abundance of ITS2 reads across the species of the group. We conclude that a segregation of the trophic niche has not occurred in the R. decipiens radiation.

Chemical Characterization of the Lichen-Symbiont Microalga Asterochloris erici and Study of Its Cytostatic Effect on the L929 Murine Fibrosarcoma Cell Line

New resources of food, pharmaceuticals or biotechnological products are needed. The huge biodiversity of aero-terrestrial lichen-symbiont microalgae belonging to the Chlorophyta group remains unexplored despite they present interesting features such as extreme stress tolerance and growth in water shortage. Appropriateness for human consumption demands the demonstration of the absence of toxic effects. In vitro biocompatibility of crude homogenates of axenic microalga Asterochloris erici, isolated from the lichen Cladonia cristatella, was analyzed after treatment of cultured L929 fibroblasts with different concentrations of microalgal homogenates. The microalgal protein content (37%) was similar to spirulina or soybean. Antioxidant capacity (10.6 ± 0.6 µmol TE/g WW) or phenolic content (7.5 ± 0.5 mg GAE/g DW) were high compared to Chlorella. The results show that crude homogenates of A. erici do not induce cytotoxicity but seem to have some cytostatic effect inducing slight cell cycle alterations and intracellular reactive oxygen species (ROS) increase at the highest concentration. Carotenoid analysis demonstrates high contents of lutein (1211 µg/g microalga DW), a xanthophyll with antioxidant and cytostatic properties in vivo and high commercial added value. These findings confirm that Asterochloris erici can be suitable for the development of alimentary or pharmaceutical applications and further in vivo animal testing. The cytostatic effects should be further investigated for antitumor agents.

Macroecological diversification and convergence in a clade of keystone symbionts

Lichens are classic models of symbiosis, and one of the most frequent nutritional modes among fungi. The ecologically and geographically widespread lichen-forming algal (LFA) genus Trebouxia is one of the best-studied groups of LFA and associates with over 7000 fungal species. Despite its importance, little is known about its diversification. We synthesized twenty years of publicly available data by characterizing the ecological preferences of this group and testing for time-variant shifts in climatic regimes over a distribution of trees. We found evidence for limited shifts among regimes, but that disparate lineages convergently evolved similar ecological tolerances. Early Trebouxia lineages were largely forest specialists or habitat generalists that occupied a regime whose extant members occur in moderate climates. Trebouxia then convergently diversified in non-forested habitats and expanded into regimes whose modern representatives occupy wet-warm and cool-dry climates. We rejected models in which climatic diversification slowed through time, suggesting climatic diversification is inconsistent with that expected under an adaptive radiation. In addition, we found that climatic and vegetative regime shifts broadly coincided with the evolution of biomes and associated or similar taxa. Together, our work illustrates how this keystone symbiont from an iconic symbiosis evolved to occupy diverse habitats across the globe.

Thallus Growth Stage and Geographic Origin Shape Microalgal Diversity in Ramalina farinacea Lichen Holobionts

Lichen symbioses are microecosystems hosting many other living organisms besides the two major lichen symbionts (i.e., lichenized fungi [the mycobiont] and green microalgae or cyanobacteria [the photobiont]). Recent investigations evidenced that other fungi, non-photosynthetic bacteria, and microalgae co-inhabit within the lichen thalli, but their diversity and their roles are still underinvestigated. Here we present an ad hoc stratified sampling design and in-depth Illumina paired-end metabarcoding approach to explore microalgal diversity in lichen thalli of the model species Ramalina farinacea from different ecologies. Lichen thalli were surveyed according to three different sizes, and different thallus parts were considered for molecular, bioinformatics, and community diversity analyses. The results revealed that microalgal diversity strongly depends on the growth stage of the thalli, the geographic area, and the habitat type. The results also show that microalgal diversity does not vary along the thallus branches (lacinias)-that is, it does not correlate with the apical growth and founder effects-and that there is no balanced co-presence of two main photobionts as previously established in R. farinacea. The sampling design performed here minimizes bias in the assessment of photobiont diversity in lichens and is proposed to be reliable and applicable to further study microalgal diversity in lichen symbioses.

Extensive photobiont sharing in a rapidly radiating cyanolichen clade

Recent studies have uncovered remarkable diversity in Dictyonema s.lat. basidiolichens, here recognized as subtribe Dictyonemateae. This group includes five genera and 148 species, but hundreds more await description. The photobionts of these lichens belong to Rhizonema, a recently resurrected cyanobacterial genus known by a single species. To further investigate photobiont diversity within Dictyonemateae, we generated 765 new cyanobacterial sequences from 635 specimens collected from 18 countries. The ITS barcoding locus supported the recognition of 200 mycobiont (fungal) species among these samples, but the photobiont diversity was comparatively low. Our analyses revealed three main divisions of Rhizonema, with two repeatedly recovered as monophyletic (proposed as new species), and the third mostly paraphyletic. The paraphyletic lineage corresponds to R. interruptum and partnered with mycobionts from all five genera in Dictyonemateae. There was no evidence of photobiont-mycobiont co-speciation, but one of the monophyletic lineages of Rhizonema appears to partner predominantly with one of the two major clades of Cora (mycobiont) with samples collected largely from the northern Andes. Molecular clock estimations indicate the Rhizonema species are much older than the fungal species in the Dictyonemateae, suggesting that these basidiolichens obtained their photobionts from older ascolichen lineages and the photobiont variation in extant lineages of Dictyonemateae is the result of multiple photobiont switches. These results support the hypothesis of lichens representing "fungal farmers," in which diverse mycobiont lineages associate with a substantially lower diversity of photobionts by sharing those photobionts best suited for the lichen symbiosis among multiple and often unrelated mycobiont lineages.

Symbiotic microalgal diversity within lichenicolous lichens and crustose hosts on Iberian Peninsula gypsum biocrusts

This study analyses the interactions among crustose and lichenicolous lichens growing on gypsum biocrusts. The selected community was composed of Acarospora nodulosa, Acarospora placodiiformis, Diploschistes diacapsis, Rhizocarpon malenconianum and Diplotomma rivas-martinezii. These species represent an optimal system for investigating the strategies used to share phycobionts because Acarospora spp. are parasites of D. diacapsis during their first growth stages, while in mature stages, they can develop independently. R. malenconianum is an obligate lichenicolous lichen on D. diacapsis, and D. rivas-martinezii occurs physically close to D. diacapsis. Microalgal diversity was studied by Sanger sequencing and 454-pyrosequencing of the nrITS region, and the microalgae were characterized ultrastructurally. Mycobionts were studied by performing phylogenetic analyses. Mineralogical and macro- and micro-element patterns were analysed to evaluate their influence on the microalgal pool available in the substrate. The intrathalline coexistence of various microalgal lineages was confirmed in all mycobionts. D. diacapsis was confirmed as an algal donor, and the associated lichenicolous lichens acquired their phycobionts in two ways: maintenance of the hosts' microalgae and algal switching. Fe and Sr were the most abundant microelements in the substrates but no significant relationship was found with the microalgal diversity. The range of associated phycobionts are influenced by thallus morphology.

Formally described species woefully underrepresent phylogenetic diversity in the common lichen photobiont genus Trebouxia (Trebouxiophyceae, Chlorophyta): An impetus for developing an integrated taxonomy

Lichens provide valuable systems for studying symbiotic interactions. In lichens, these interactions are frequently described in terms of availability, selectivity and specificity of the mycobionts and photobionts towards one another. The lichen-forming, green algal genus Trebouxia Puymaly is among the most widespread photobiont, associating with a broad range of lichen-forming fungi. To date, 29 species have been described, but studies consistently indicate that the vast majority of species-level lineages still lack formal description, and new, previously unrecognized lineages are frequently reported. To reappraise the diversity and the evolutionary relationships of species-level lineages in Trebouxia, we assembled DNA sequence data from over 1600 specimens, compiled from a range of sequences from previously published studies, axenic algal cultures, and lichens collected from poorly sampled regions. From these samples, we selected representatives of the currently known genetic diversity in the lichenized Trebouxia and inferred a phylogeny from multi-locus sequence data (ITS, rbcL, cox2). We demonstrate that the current formally described species woefully underrepresent overall species-level diversity in this important lichen-forming algal genus. We anticipate that an integrative taxonomic approach, incorporating morphological and physiological data from axenic cultures with genetic data, will be required to establish a robust, comprehensive taxonomy for Trebouxia. The data presented here provide an important impetus and reference dataset for more reliably characterizing diversity in lichenized algae and in using lichens to investigate the evolution of symbioses and holobionts.

The chloroplast genome of the lichen-symbiont microalga Trebouxia sp. Tr9 (Trebouxiophyceae, Chlorophyta) shows short inverted repeats with a single gene and loss of the rps4 gene, which is encoded by the nucleus

The Trebouxiophyceae is the class of Chlorophyta algae from which the highest number of chloroplast genome (cpDNA) sequences has been obtained. Several species in this class participate in symbioses with fungi to form lichens. However, no cpDNA has been obtained from any Trebouxia lichen-symbiont microalgae, which are present in approximately half of all lichens. Here, we report the sequence of the completely assembled cpDNA from Trebouxia sp. TR9 and a comparative study with other Trebouxio-phyceae. The organization of the chloroplast genome of Trebouxia sp. TR9 has certain features that are unusual in the Trebouxiophyceae and other green algae. The most remarkable characteristics are the presence of long intergenic spacers, a quadripartite structure with short inverted repeated sequences (IRs), and the loss of the rps4 gene. The presence of long intergenic spacers accounts for a larger cpDNA size in comparison to other closely related Trebouxiophyceae. The IRs, which were thought to be lost in the Trebouxiales, are distinct from most of cpDNAs since they lack the rRNA operon and uniquely includes the rbcL gene. The functional transfer of the rps4 gene to the nuclear genome has been confirmed by sequencing and examination of the gene architecture, which includes three spliceosomal introns as well as the verification of the presence of the corresponding transcript. This is the first documented transfer of the rps4 gene from the chloroplast to the nucleus among Viridiplantae. Additionally, a fairly well-resolved phylogenetic reconstruction, including Trebouxia sp. TR9 along with other Trebouxiophyceae, was obtained based on a set of conserved chloroplast genes.

Holobiont chronobiology: mycorrhiza may be a key to linking aboveground and underground rhythms

Circadian clocks are nearly ubiquitous timing mechanisms that can orchestrate rhythmic behavior and gene expression in a wide range of organisms. Clock mechanisms are becoming well understood in fungal, animal, and plant model systems, yet many of these organisms are surrounded by a complex and diverse microbiota which should be taken into account when examining their biology. Of particular interest are the symbiotic relationships between organisms that have coevolved over time, forming a unit called a holobiont. Several studies have now shown linkages between the circadian rhythms of symbiotic partners. Interrelated regulation of holobiont circadian rhythms seems thus important to coordinate shifts in activity over the day for all the partners. Therefore, we suggest that the classical view of "chronobiological individuals" should include "a holobiont" rather than an organism. Unfortunately, mechanisms that may regulate interspecies temporal acclimation and the evolution of the circadian clock in holobionts are far from being understood. For the plant holobiont, our understanding is particularly limited. In this case, the holobiont encompasses two different ecosystems, one above and the other below the ground, with the two potentially receiving timing information from different synchronizing signals (Zeitgebers). The arbuscular mycorrhizal (AM) symbiosis, formed by plant roots and fungi, is one of the oldest and most widespread associations between organisms. By mediating the nutritional flux between the plant and the many microbes in the soil, AM symbiosis constitutes the backbone of the plant holobiont. Even though the importance of the AM symbiosis has been well recognized in agricultural and environmental sciences, its circadian chronobiology remains almost completely unknown. We have begun to study the circadian clock of arbuscular mycorrhizal fungi, and we compile and here discuss the available information on the subject. We propose that analyzing the interrelated temporal organization of the AM symbiosis and determining its underlying mechanisms will advance our understanding of the role and coordination of circadian clocks in holobionts in general.

Dynamic evolution of mitochondrial genomes in Trebouxiophyceae, including the first completely assembled mtDNA from a lichen-symbiont microalga (Trebouxia sp. TR9)

Trebouxiophyceae (Chlorophyta) is a species-rich class of green algae with a remarkable morphological and ecological diversity. Currently, there are a few completely sequenced mitochondrial genomes (mtDNA) from diverse Trebouxiophyceae but none from lichen symbionts. Here, we report the mitochondrial genome sequence of Trebouxia sp. TR9 as the first complete mtDNA sequence available for a lichen-symbiont microalga. A comparative study of the mitochondrial genome of Trebouxia sp. TR9 with other chlorophytes showed important organizational changes, even between closely related taxa. The most remarkable change is the enlargement of the genome in certain Trebouxiophyceae, which is principally due to larger intergenic spacers and seems to be related to a high number of large tandem repeats. Another noticeable change is the presence of a relatively large number of group II introns interrupting a variety of tRNA genes in a single group of Trebouxiophyceae, which includes Trebouxiales and Prasiolales. In addition, a fairly well-resolved phylogeny of Trebouxiophyceae, along with other Chlorophyta lineages, was obtained based on a set of seven well-conserved mitochondrial genes.

Reproductive and dispersal strategies shape the diversity of mycobiont-photobiont association in Cladonia lichens

Ecological preferences, partner compatibility, or partner availability are known to be important factors shaping obligate and intimate lichen symbioses. We considered a complex of Cladonia species, traditionally differentiated by the extent of sexual reproduction and the type of vegetative propagules, to assess if the reproductive and dispersal strategies affect mycobiont-photobiont association patterns. In total 85 lichen thalli from 72 European localities were studied, two genetic markers for both Cladonia mycobionts and Asterochloris photobionts were analyzed. Variance partitioning analysis by multiple regression on distance matrices was performed to describe and partition variance in photobiont genetic diversity. Asexually reproducing Cladonia in our study were found to be strongly specific to their photobionts, associating with only two closely related Asterochloris species. In contrast, sexually reproducing lichens associated with seven unrelated Asterochloris lineages, thus being photobiont generalists. The reproductive mode had the largest explanatory power, explaining 44% of the total photobiont variability. Reproductive and dispersal strategies are the key factors shaping photobiont diversity in this group of Cladonia lichens. A strict photobiont specialisation observed in two studied species may steer both evolutionary flexibility and responses to ecological changes of these organisms, and considerably limit their distribution ranges.

Comparing Sanger sequencing and high-throughput metabarcoding for inferring photobiont diversity in lichens

The implementation of HTS (high-throughput sequencing) approaches is rapidly changing our understanding of the lichen symbiosis, by uncovering high bacterial and fungal diversity, which is often host-specific. Recently, HTS methods revealed the presence of multiple photobionts inside a single thallus in several lichen species. This differs from Sanger technology, which typically yields a single, unambiguous algal sequence per individual. Here we compared HTS and Sanger methods for estimating the diversity of green algal symbionts within lichen thalli using 240 lichen individuals belonging to two species of lichen-forming fungi. According to HTS data, Sanger technology consistently yielded the most abundant photobiont sequence in the sample. However, if the second most abundant photobiont exceeded 30% of the total HTS reads in a sample, Sanger sequencing generally failed. Our results suggest that most lichen individuals in the two analyzed species, Lasallia hispanica and L. pustulata, indeed contain a single, predominant green algal photobiont. We conclude that Sanger sequencing is a valid approach to detect the dominant photobionts in lichen individuals and populations. We discuss which research areas in lichen ecology and evolution will continue to benefit from Sanger sequencing, and which areas will profit from HTS approaches to assessing symbiont diversity.

Secondary metabolism in the lichen symbiosis

Lichens, which are defined by a core symbiosis between a mycobiont (fungal partner) and a photobiont (photoautotrophic partner), are in fact complex assemblages of microorganisms that constitute a largely untapped source of bioactive secondary metabolites. Historically, compounds isolated from lichens have predominantly been those produced by the dominant fungal partner, and these continue to be of great interest for their unique chemistry and biotechnological potential. In recent years it has become apparent that many photobionts and lichen-associated bacteria also produce a range of potentially valuable molecules. There is evidence to suggest that the unique nature of the symbiosis has played a substantial role in shaping many aspects of lichen chemistry, for example driving bacteria to produce metabolites that do not bring them direct benefit but are useful to the lichen as a whole. This is most evident in studies of cyanobacterial photobionts, which produce compounds that differ from free living cyanobacteria and are unique to symbiotic organisms. The roles that these and other lichen-derived molecules may play in communication and maintaining the symbiosis are poorly understood at present. Nonetheless, advances in genomics, mass spectrometry and other analytical technologies are continuing to illuminate the wealth of biological and chemical diversity present within the lichen holobiome. Implementation of novel biodiscovery strategies such as metagenomic screening, coupled with synthetic biology approaches to reconstitute, re-engineer and heterologously express lichen-derived biosynthetic gene clusters in a cultivable host, offer a promising means for tapping into this hitherto inaccessible wealth of natural products.

Environment and host identity structure communities of green algal symbionts in lichens

An understanding of how biotic interactions shape species' distributions is central to predicting host-symbiont responses under climate change. Switches to locally adapted algae have been proposed to be an adaptive strategy of lichen-forming fungi to cope with environmental change. However, it is unclear how lichen photobionts respond to environmental gradients, and whether they play a role in determining the fungal host's upper and lower elevational limits. Deep-coverage Illumina DNA metabarcoding was used to track changes in the community composition of Trebouxia algae associated with two phylogenetically closely related, but ecologically divergent fungal hosts along a steep altitudinal gradient in the Mediterranean region. We detected the presence of multiple Trebouxia species in the majority of thalli. Both altitude and host genetic identity were strong predictors of photobiont community assembly in these two species. The predominantly clonally dispersing fungus showed stronger altitudinal structuring of photobiont communities than the sexually reproducing host. Elevation ranges of the host were not limited by the lack of compatible photobionts. Our study sheds light on the processes guiding the formation and distribution of specific fungal-algal combinations in the lichen symbiosis. The effect of environmental filtering acting on both symbiotic partners appears to shape the distribution of lichens.

ITS1 metabarcoding highlights low specificity of lichen mycobiomes at a local scale

As self-supporting and long-living symbiotic structures, lichens provide a habitat for many other organisms beside the traditionally considered lichen symbionts-the myco- and the photobionts. The lichen-inhabiting fungi either develop diagnostic phenotypes or occur asymptomatically. Because the degree of specificity towards the lichen host is poorly known, we studied the diversity of these fungi among neighbouring lichens on rocks in an alpine habitat. Using a sequencing metabarcoding approach, we show that lichen mycobiomes clearly reflect the overlap of multiple ecological sets of taxa, which differ in their trophic association with lichen thalli. The lack of specificity to the lichen mycobiome is further supported by the lack of community structure observed using clustering and ordination methods. The communities encountered across samples largely result from the subsampling of a shared species pool, in which we identify three major ecological components: (i) a generalist environmental pool, (ii) a lichenicolous/endolichenic pool and (iii) a pool of transient species. These taxa majorly belong to the fungal classes Dothideomycetes, Eurotiomycetes and Tremellomycetes with close relatives in adjacent ecological niches. We found no significant evidence that the phenotypically recognized lichenicolous fungi influence the occurrence of the other asymptomatic fungi in the host thalli. We claim that lichens work as suboptimal habitats or as a complex spore and mycelium bank, which modulate and allow the regeneration of local fungal communities. By performing an approach that minimizes ambiguities in the taxonomic assignments of fungi, we present how lichen mycobiomes are also suitable targets for improving bioinformatic analyses of fungal metabarcoding.

Unexpected associated microalgal diversity in the lichen Ramalina farinacea is uncovered by pyrosequencing analyses

The current literature reveals that the intrathalline coexistence of multiple microalgal taxa in lichens is more common than previously thought, and additional complexity is supported by the coexistence of bacteria and basidiomycete yeasts in lichen thalli. This replaces the old paradigm that lichen symbiosis occurs between a fungus and a single photobiont. The lichen Ramalina farinacea has proven to be a suitable model to study the multiplicity of microalgae in lichen thalli due to the constant coexistence of Trebouxia sp. TR9 and T. jamesii in long-distance populations. To date, studies involving phycobiont diversity within entire thalli are based on Sanger sequencing, but this method seems to underestimate the diversity. Here, we aim to analyze both the microalgal diversity and its community structure in a single thallus of the lichen R. farinacea by applying a 454 pyrosequencing approach coupled with a careful ad hoc-performed protocol for lichen sample processing prior to DNA extraction. To ascertain the reliability of the pyrosequencing results and the applied bioinformatics pipeline results, the thalli were divided into three sections (apical, middle and basal zones), and a mock community sample was used. The developed methodology allowed 40448 filtered algal reads to be obtained from a single lichen thallus, which encompassed 31 OTUs representative of different microalgae genera. In addition to corroborating the coexistence of the two Trebouxia sp. TR9 and T. jamesii taxa in the same thallus, this study showed a much higher microalgal diversity associated with the lichen. Along the thallus ramifications, we also detected variations in phycobiont distribution that might correlate with different microenvironmental conditions. These results highlight R. farinacea as a suitable material for studying microalgal diversity and further strengthen the concept of lichens as multispecies microecosystems. Future analyses will be relevant to ecophysiological and evolutionary studies to understand the roles of the multiple photobionts in lichen symbioses.

Fungal and cyanobacterial gene expression in a lichen symbiosis: Effect of temperature and location

Organisms have evolved different cellular mechanisms to deal with environmental stress, primarily through complex molecular mechanisms including protein refolding and DNA repair. As mutualistic symbioses, lichens offer the possibility of analyzing molecular stress responses in a particularly tight interspecific relationship. We study the widespread cyanolichen Peltigera membranacea, a key player in carbon and nitrogen cycling in terrestrial ecosystems at northern latitudes. We ask whether increased temperature is reflected in mRNA levels of selected damage control genes, and do the response patterns show geographical associations? Using real-time PCR quantification of 38 transcripts, differential expression was demonstrated for nine cyanobacterial and nine fungal stress response genes (plus the fungal symbiosis-related lec2 gene) when the temperature was increased from 5 °C to 15 °C and 25 °C. Principle component analysis (PCA) revealed two gene groups with different response patterns. Whereas a set of cyanobacterial DNA repair genes and the fungal lec2 (PC1 group) showed an expression drop at 15 °C vs. 5 °C, most fungal candidates (PC2 group) showed increased expression at 25 °C vs. 5 °C. PC1 responses also correlated with elevation. The correlated downregulation of lec2 and cyanobacterial DNA repair genes suggests a possible interplay between the symbionts warranting further studies.

Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi

We studied the evolutionary history of the Parmeliaceae (Lecanoromycetes, Ascomycota), one of the largest families of lichen-forming fungi with complex and variable morphologies, also including several lichenicolous fungi. We assembled a six-locus data set including nuclear, mitochondrial and low-copy protein-coding genes from 293 operational taxonomic units (OTUs). The lichenicolous lifestyle originated independently three times in lichenized ancestors within Parmeliaceae, and a new generic name is introduced for one of these fungi. In all cases, the independent origins occurred c. 24 million yr ago. Further, we show that the Paleocene, Eocene and Oligocene were key periods when diversification of major lineages within Parmeliaceae occurred, with subsequent radiations occurring primarily during the Oligocene and Miocene. Our phylogenetic hypothesis supports the independent origin of lichenicolous fungi associated with climatic shifts at the Oligocene-Miocene boundary. Moreover, diversification bursts at different times may be crucial factors driving the diversification of Parmeliaceae. Additionally, our study provides novel insight into evolutionary relationships in this large and diverse family of lichen-forming ascomycetes.

Differences in the cell walls and extracellular polymers of the two Trebouxia microalgae coexisting in the lichen Ramalina farinacea are consistent with their distinct capacity to immobilize extracellular Pb

Trebouxia TR1 and T. TR9 are the two microalgae that coexist within Ramalina farinacea thalli. In the presence of Pb, TR9 formed extracellular aggregates, while TR1 showed a lower wall lead retention capability. Herein, we studied the cell walls and extracellular polymers (EPS) of TR1 and TR9, and their possible implication in the different Pb retention capacity of these microalgae. The proportion of cell walls on the overall cell biomass was 2.6 times higher in TR9 than in TR1. Glycosyl linkage analysis indicated the presence of hot-water soluble β-galactofuranan(s) in both cell walls, distinct from that previously described in Asterochloris, with increased rhamnose content in TR9 and a higher substitution by side chains in TR1. In the EPS of both microalgae, monosaccharide composition resembled their cell walls; however, TR9 produced two times more EPS than TR1. EPS also contained proteins, displaying distinct polypeptide patterns, modulated by Pb in TR9. Uronic acid content in EPS was higher in TR9, and increased in this phycobiont after the addition of Pb. Our findings suggest that the differential extracellular Pb retention capabilities of TR1 and TR9 microalgae can be related to differences in the main features of their cell walls and EPS.

Molecular phylogeny and ultrastructure of the lichen microalga Asterochloris mediterranea sp. nov. from Mediterranean and Canary Islands ecosystems

The microalgae of the genus Asterochloris are the preferential phycobionts in Cladonia, Lepraria and Stereocaulon lichens. Recent studies have highlighted the hidden diversity of the genus, even though phycobionts hosting species of the genus Cladonia in Mediterranean and Canarian ecosystems have been poorly explored. Phylogenetic analyses were made by concatenation of the sequences obtained with a plastid – LSU rDNA – and two nuclear – internal transcribed spacer (ITS) rDNA and actin – molecular markers of the phycobionts living in several populations of the Cladonia convoluta-Cladonia foliacea complex, Cladonia rangiformis and Cladonia cervicornis s. str. widely distributed in these areas in a great variety of substrata and habitats. A new strongly supported clade was obtained in relation to the previously published Asterochloris phylogenies. Minimum genetic variation was detected between our haplotypes and other sequences available in the GenBank database. The correct identification of the fungal partners was corroborated by the ITS rDNA barcode. In this study we provide a detailed characterization comprising chloroplast morphology, and ultrastructural and phylogenetic analyses of a novel phycobiont species, here described as Asterochloris mediterranea sp. nov. Barreno, Chiva, Moya et Škaloud. A cryopreserved holotype specimen has been deposited in the Culture Collection of Algae of Charles University in Prague, Czech Republic (CAUP) as CAUP H 1015. We suggest the use of a combination of several nuclear and plastid molecular markers, as well as ultrastructural (transmission electron and confocal microscopy) techniques, both in culture and in the symbiotic state, to improve novel species delimitation of phycobionts in lichens.

Lichen rehydration in heavy metal-polluted environments: Pb modulates the oxidative response of both Ramalina farinacea thalli and its isolated microalgae

Lichens are adapted to desiccation/rehydration and accumulate heavy metals, which induce ROS especially from the photobiont photosynthetic pigments. Although their mechanisms of abiotic stress tolerance are still to be unravelled, they seem related to symbionts' reciprocal upregulation of antioxidant systems. With the aim to study the effect of Pb on oxidative status during rehydration, the kinetics of intracellular ROS, lipid peroxidation and chlorophyll autofluorescence of whole Ramalina farinacea thalli and its isolated microalgae (Trebouxia TR1 and T. TR9) was recorded. A genetic characterization of the microalgae present in the thalli used was also carried out in order to assess possible correlations among the relative abundance of each phycobiont, their individual physiological responses and that of the entire thallus. Unexpectedly, Pb decreased ROS and lipid peroxidation in thalli and its phycobionts, associated with a lower chlorophyll autofluorescence. Each phycobiont showed a particular pattern, but the oxidative response of the thallus paralleled the TR1's, agreeing with the genetic identification of this strain as the predominant phycobiont. We conclude that: (1) the lichen oxidative behaviour seems to be modulated by the predominant phycobiont and (2) Pb evokes in R. farinacea and its phycobionts strong mechanisms to neutralize its own oxidant effects along with those of rehydration.

Lichen symbiosis: nature's high yielding machines for induced hydrogen production

Hydrogen is a promising future energy source. Although the ability of green algae to produce hydrogen has long been recognized (since 1939) and several biotechnological applications have been attempted, the greatest obstacle, being the O2-sensitivity of the hydrogenase enzyme, has not yet been overcome. In the present contribution, 75 years after the first report on algal hydrogen production, taking advantage of a natural mechanism of oxygen balance, we demonstrate high hydrogen yields by lichens. Lichens have been selected as the ideal organisms in nature for hydrogen production, since they consist of a mycobiont and a photobiont in symbiosis. It has been hypothesized that the mycobiont's and photobiont's consumption of oxygen (increase of COX and AOX proteins of mitochondrial respiratory pathways and PTOX protein of chrolorespiration) establishes the required anoxic conditions for the activation of the phycobiont's hydrogenase in a closed system. Our results clearly supported the above hypothesis, showing that lichens have the ability to activate appropriate bioenergetic pathways depending on the specific incubation conditions. Under light conditions, they successfully use the PSII-dependent and the PSII-independent pathways (decrease of D1 protein and parallel increase of PSaA protein) to transfer electrons to hydrogenase, while under dark conditions, lichens use the PFOR enzyme and the dark fermentative pathway to supply electrons to hydrogenase. These advantages of lichen symbiosis in combination with their ability to survive in extreme environments (while in a dry state) constitute them as unique and valuable hydrogen producing natural factories and pave the way for future biotechnological applications.

Phylogenetic congruence of lichenised fungi and algae is affected by spatial scale and taxonomic diversity

The role of species' interactions in structuring biological communities remains unclear. Mutualistic symbioses, involving close positive interactions between two distinct organismal lineages, provide an excellent means to explore the roles of both evolutionary and ecological processes in determining how positive interactions affect community structure. In this study, we investigate patterns of co-diversification between fungi and algae for a range of New Zealand lichens at the community, genus, and species levels and explore explanations for possible patterns related to spatial scale and pattern, taxonomic diversity of the lichens considered, and the level sampling replication. We assembled six independent datasets to compare patterns in phylogenetic congruence with varied spatial extent of sampling, taxonomic diversity and level of specimen replication. For each dataset, we used the DNA sequences from the ITS regions of both the fungal and algal genomes from lichen specimens to produce genetic distance matrices. Phylogenetic congruence between fungi and algae was quantified using distance-based redundancy analysis and we used geographic distance matrices in Moran's eigenvector mapping and variance partitioning to evaluate the effects of spatial variation on the quantification of phylogenetic congruence. Phylogenetic congruence was highly significant for all datasets and a large proportion of variance in both algal and fungal genetic distances was explained by partner genetic variation. Spatial variables, primarily at large and intermediate scales, were also important for explaining genetic diversity patterns in all datasets. Interestingly, spatial structuring was stronger for fungal than algal genetic variation. As the spatial extent of the samples increased, so too did the proportion of explained variation that was shared between the spatial variables and the partners' genetic variation. Different lichen taxa showed some variation in their phylogenetic congruence and spatial genetic patterns and where greater sample replication was used, the amount of variation explained by partner genetic variation increased. Our results suggest that the phylogenetic congruence pattern, at least at small spatial scales, is likely due to reciprocal co-adaptation or co-dispersal. However, the detection of these patterns varies among different lichen taxa, across spatial scales and with different levels of sample replication. This work provides insight into the complexities faced in determining how evolutionary and ecological processes may interact to generate diversity in symbiotic association patterns at the population and community levels. Further, it highlights the critical importance of considering sample replication, taxonomic diversity and spatial scale in designing studies of co-diversification.

Photobiont selectivity leads to ecological tolerance and evolutionary divergence in a polymorphic complex of lichenized fungi

BACKGROUND AND AIMS

The integrity and evolution of lichen symbioses depend on a fine-tuned combination of algal and fungal genotypes. Geographically widespread species complexes of lichenized fungi can occur in habitats with slightly varying ecological conditions, and it remains unclear how this variation correlates with symbiont selectivity patterns in lichens. In an attempt to address this question, >300 samples were taken of the globally distributed and ecologically variable lichen-forming species complex Tephromela atra, together with closely allied species, in order to study genetic diversity and the selectivity patterns of their photobionts.

METHODS

Lichen thalli of T. atra and of closely related species T. grumosa, T. nashii and T. atrocaesia were collected from six continents, across 24 countries and 62 localities representing a wide range of habitats. Analyses of genetic diversity and phylogenetic relationships were carried out both for photobionts amplified directly from the lichen thalli and from those isolated in axenic cultures. Morphological and anatomical traits were studied with light and transmission electron microscopy in the isolated algal strains.

KEY RESULTS

Tephromela fungal species were found to associate with 12 lineages of Trebouxia. Five new clades demonstrate the still-unrecognized genetic diversity of lichen algae. Culturable, undescribed lineages were also characterized by phenotypic traits. Strong selectivity of the mycobionts for the photobionts was observed in six monophyletic Tephromela clades. Seven Trebouxia lineages were detected in the poorly resolved lineage T. atra sensu lato, where co-occurrence of multiple photobiont lineages in single thalli was repeatedly observed.

CONCLUSIONS

Low selectivity apparently allows widespread lichen-forming fungi to establish successful symbioses with locally adapted photobionts in a broader range of habitats. This flexibility might correlate with both lower phylogenetic resolution and evolutionary divergence in species complexes of crustose lichen-forming fungi.

Integrating coalescent and phylogenetic approaches to delimit species in the lichen photobiont Trebouxia

The accurate assessment of species boundaries in symbiotic systems is a prerequisite for the study of speciation, co-evolution and selectivity. Many studies have shown the high genetic diversity of green algae from the genus Trebouxia, the most common photobiont of lichen-forming fungi. However, the phylogenetic relationships, and the amount of cryptic diversity of these algae are still poorly understood, and an adequate species concept for trebouxiophycean algae is still missing. In this study we used a multifaceted approach based on coalescence (GMYC, STEM) and phylogenetic relationships to assess species boundaries in the trebouxioid photobionts of the lichen-forming fungus Lasallia pustulata. We further investigated whether putative species of Trebouxia found in L. pustulata are shared with other lichen-forming fungi. We found that L. pustulata is associated with at least five species of Trebouxia and most of them are shared with other lichen-forming fungi, showing different patterns of species-to-species and species-to-community interactions. We also show that one of the putative Trebouxia species is found exclusively in association with L. pustulata and is restricted to thalli from localities with Mediterranean microclimate. We suggest that the species delimitation method presented in this study is a promising tool to address species boundaries within the heterogeneous genus Trebouxia.

Do photobiont switch and cephalodia emancipation act as evolutionary drivers in the lichen symbiosis? A case study in the Pannariaceae (Peltigerales)

Lichen symbioses in the Pannariaceae associate an ascomycete and either cyanobacteria alone (usually Nostoc; bipartite thalli) or green algae and cyanobacteria (cyanobacteria being located in dedicated structures called cephalodia; tripartite thalli) as photosynthetic partners (photobionts). In bipartite thalli, cyanobacteria can either be restricted to a well-delimited layer within the thallus ('pannarioid' thalli) or spread over the thallus that becomes gelatinous when wet ('collematoid' thalli). We studied the collematoid genera Kroswia and Physma and an undescribed tripartite species along with representatives of the pannarioid genera Fuscopannaria, Pannaria and Parmeliella. Molecular inferences from 4 loci for the fungus and 1 locus for the photobiont and statistical analyses within a phylogenetic framework support the following: (a) several switches from pannarioid to collematoid thalli occured and are correlated with photobiont switches; the collematoid genus Kroswia is nested within the pannarioid genus Fuscopannaria and the collematoid genus Physma is sister to the pannarioid Parmeliella mariana group; (b) Nostoc associated with collematoid thalli in the Pannariaceae are related to that of the Collemataceae (which contains only collematoid thalli), and never associated with pannarioid thalli; Nostoc associated with pannarioid thalli also associate in other families with similar morphology; (c) ancestors of several lineages in the Pannariaceae developed tripartite thalli, bipartite thalli probably resulting from cephalodia emancipation from tripartite thalli which eventually evolved and diverged, as suggested by the same Nostoc present in the collematoid genus Physma and in the cephalodia of a closely related tripartite species; Photobiont switches and cephalodia emancipation followed by divergence are thus suspected to act as evolutionary drivers in the family Pannariaceae.

Evolutionary implications of intron-exon distribution and the properties and sequences of the RPL10A gene in eukaryotes

The RPL10A gene encodes the RPL10 protein, required for joining 40S and 60S subunits into a functional 80S ribosome. This highly conserved gene, ubiquitous across all eukaryotic super-groups, is characterized by a variable number of spliceosomal introns, present in most organisms. These properties facilitate the recognition of orthologs among distant taxa and thus comparative studies of sequences as well as the distribution and properties of introns in taxonomically distant groups of eukaryotes. The present study examined the multiple ways in which RPL10A conservation vs. sequence changes in the gene over the course of evolution, including in exons, introns, and the encoded proteins, can be exploited for evolutionary analysis at different taxonomic levels. At least 25 different positions harboring introns within the RPL10A gene were determined in different taxa, including animals, plants, fungi, and alveolates. Generally, intron positions were found to be well conserved even across different kingdoms. However, certain introns seemed to be restricted to specific groups of organisms. Analyses of several properties of introns, including insertion site, phase, and length, along with exon and intron GC content and exon-intron boundaries, suggested biases within different groups of organisms. The use of a standard primer pair to analyze a portion of the intron-containing RPL10A gene in 12 genera of green algae within Chlorophyta is presented as a case study for evolutionary analyses of introns at intermediate and low taxonomic levels. Our study shows that phylogenetic reconstructions at different depths can be achieved using RPL10A nucleotide sequences from both exons and introns as well as the amino acid sequences of the encoded protein.