Next Article in Journal
Just around the Corner: Advances in the Optimization of Yeasts and Filamentous Fungi for Lactic Acid Production
Next Article in Special Issue
The First Miniature, Small Foliose, Brown Xanthoparmelia in the Northern Hemisphere
Previous Article in Journal
Phylogenetic and Taxonomic Analyses of Five New Wood-Inhabiting Fungi of Botryobasidium, Coltricia and Coltriciella (Basidiomycota) from China
Previous Article in Special Issue
Contribution to the Study of Lichenicolous Fungi from Northwest Iberian Peninsula (León and Lugo Provinces)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Myco–Phycobiont Interactions within the “Ramalina farinacea Group”: A Geographical Survey over Europe and Macaronesia

1
Instituto Cavanilles de Biodiversidad y Biología Evolutiva (ICBiBE)—Departament de Botànica, Universitat de València, C/Dr. Moliner, 50, E-46100 Burjassot, València, Spain
2
Department of Life Sciences, University of Trieste, Via L. Giorgieri 10, 34127 Trieste, Italy
3
Instituto de Biotecnología y Biomedicina (BIOTECMED), Universitat de València, E-46100 Burjassot, Spain
4
Departament de Botànica i Geologia, Universitat de València, C/Dr. Moliner, 50, E-46100 Burjassot, València, Spain
*
Author to whom correspondence should be addressed.
J. Fungi 2024, 10(3), 206; https://doi.org/10.3390/jof10030206
Submission received: 31 January 2024 / Revised: 4 March 2024 / Accepted: 6 March 2024 / Published: 8 March 2024
(This article belongs to the Special Issue Lichen Forming Fungi—in Honour of Prof. Ana Rosa Burgaz)

Abstract

:
Ramalina farinacea is a widely distributed epiphytic lichen from the Macaronesian archipelagos to Mediterranean and Boreal Europe. Previous studies have indicated a specific association between R. farinacea and Trebouxia microalgae species. Here, we examined the symbiotic interactions in this lichen and its closest allies (the so-called “R. farinacea group”) across ten biogeographic subregions, spanning diverse macroclimates, analyzing the climatic niche of the primary phycobionts, and discussing the specificity of these associations across the studied area. The most common phycobionts in the “R. farinacea group” were T. jamesii and T. lynnae, which showed a preference for continentality and insularity, respectively. The Canarian endemic R. alisiosae associated exclusively with T. lynnae, while the other Ramalina mycobionts interacted with both microalgae. The two phycobionts exhibited extensive niche overlap in an area encompassing Mediterranean, temperate Europe, and Macaronesian localities. However, T. jamesii occurred in more diverse climate types, whereas T. lynnae preferred warmer and more humid climates, often close to the sea, which could be related to its tolerance to salinity. With the geographical perspective gained in this study, it was possible to show how the association with different phycobionts may shape the ecological adaptation of lichen symbioses.

1. Introduction

The lichen-forming fungus Ramalina farinacea (L.) Ach. is an iconic member of the family Ramalinaceae (Lecanorales, Ascomycota) [1] that develops pendant, whitish–greenish fruticose thalli composed of slender thallus ramifications (i.e., laciniae), which usually display well-developed asexual reproduction structures, i.e., soralia. These structures consist of soredia, minute clusters of fungal hyphae enwrapping microalgal cells, dispersed together by wind or animals [2,3]. Sexual reproduction involving meiotic spores produced in fungal apothecia is, however, exceptionally rare in this lichen. Ramalina farinacea is widely distributed across North America, Macaronesia, and Europe, from the climatically mild Mediterranean Basin to the cooler Boreal region [4], showing remarkable ecological versatility as it thrives in diverse habitats from tree trunks and branches in shaded deciduous forests to isolated trees exposed to sun and wind, as well as hedgerows, scrublands, rocks, and walls [5]. Although this lichen is relatively easily identifiable in the field, its populations also show variability in morphology and chemical composition, producing various secondary metabolites [6]. That variability was the basis for the description of the related species R. subfarinacea (Nyl. ex Cromb.) Nyl. and R. alisiosae by Pérez-Vargas and Pérez-Ortega [7,8]. In a recent paper, Moya et al. [4] showed that all these species together with an undescribed Ramalina sp. are indeed phylogenetically closely related, as they likely evolved recently and therefore still represent a species complex, the so-called “Ramalina farinacea group”. Indeed, phylogenetic analyses using a pair of genetic markers, including the fungal barcode nuclear ribosomal internal transcribed spacer (nrITS), suggested incomplete lineage sorting and conflicting phylogenetic signal [4].
The work by Moya et al. [4] also studied the phylogeography of the fungal partners (mycobionts) in the Ramalina farinacea group. In particular, they evaluated the hypothesis proposed by Krog and Østhagen [9], posing an origin for R. farinacea in the Macaronesian–Mediterranean region and its gradual geographic expansion to most temperate and Boreal regions of the Northern Hemisphere, with the Canary Islands as its probable southernmost limit in the Atlantic region. Aptroot and Schumm [10] have extended this boundary to include Cape Verde, another Macaronesian archipelago south of the Canary Islands. Preliminary phylogenetic studies by Del Campo et al. [11] and Molins et al. [12] had identified several clades in the mycobiont phylogeny, some with specimens restricted to the Iberian Peninsula or the Canary Islands, whereas at least one clade was found to be shared among these two regions and California. Finally, Moya et al. [4] extended the sampling to include other areas of the Mediterranean Basin and higher latitudes in Europe and found a great number of haplotypes restricted to the Macaronesian region. Altogether, their results suggest that these lichenized fungi probably originated in southern latitudes during relatively recent geological times (Pleistocene) and then expanded north, thus supporting the original hypothesis by Krog and Østhagen [9].
In recent decades, R. farinacea has also received considerable attention regarding the diversity and physiology of its associated microalgae (phycobionts) [12]. del Campo et al. [11] and Moya et al. [13] have found that Trebouxia Puymaly diversity and the composition in thalli of R. farinacea was strongly correlated with the geographic origin of the samples. Trebouxia jamesii (Hildreth and Ahmadjian) Gärtner was the main primary phycobiont in thalli from the Iberian Peninsula, while T. lynnae Barreno (former Trebouxia sp. TR9) was mainly present in thalli from the Canary Islands and Madeira. Furthermore, Casano et al. [14] and del Campo et al. [11] have demonstrated the coexistence of these two Trebouxia phycobionts within a single lichen thallus through a combination of microscopic techniques, isolation in axenic cultures, and molecular characterization [5,11]. Casano et al. [14] have found that these two phycobionts exhibited distinct responses to abiotic stresses, and they speculated that their concurrent presence in thalli is advantageous for the holobiont, a phenomenon that has been observed in other lichens with wide distributions and ecologies [15,16,17]. So far, no further studies have attempted to unveil the diversity of the associated Trebouxia in R. farinacea from areas without the influence of the Mediterranean climate.
The present work provides a global picture of phycobiont diversity and mycobiont–phycobiont interaction patterns in the R. farinacea group by considering an enlarged dataset consisting of samples included in Moya et al. [4], newly collected thalli, as well as publicly available sequence data in GenBank. Because previous studies have shown that lichen-forming fungi tend to preferentially select photobiont lineages which are better adapted to local environmental conditions [18,19,20,21,22,23], we will here test the following hypotheses: (i) mycobionts in the R. farinacea group associate differentially with phycobionts in Macaronesia and Mediterranean, temperate, and Boreal regions in Europe; and (ii) the association with different phycobionts shapes the ecological distribution of the mycobionts. In doing this, we will investigate the complexity of interactions among mycobiont and phycobiont lineages, the degree of specificity across their geographic distributions [24,25,26], and the climatic niche width of the primary phycobionts.

2. Materials and Methods

2.1. Sampling, Pretreatment of the Samples, and DNA Extraction

The present study considered 469 thalli of the “Ramalina farinacea group”, following the species concepts in Moya et al. [4]: 434 Ramalina farinacea, 25 Ramalina sp., 5 R. subfarinacea, and 5 R. alisiosae. Sampling included the Canary Islands and Madeira (Macaronesian archipelagos), the Mediterranean Basin (Algeria, the Iberian and Italian peninsulas, and the Balearic Islands), and central (France, Czech Republic, Austria, and Germany) and northern (Estonia, Finland, Sweden, and Norway) Europe (Figure 1; Table S1). The sampling localities can be grouped into ten biogeographic subregions, according to Rivas-Martínez et al. [27] (Figure S1): Boreal, Atlantic Europe, Central Europe, Alpine, Mediterranean West and Central Iberian Peninsula, Balearic-Catalonian-Provençal, Italo-Thyrrenian, Adriatic, and Macaronesian. Fresh collected medium size thalli (up to 5 cm) were air-dried and then stored at −20 °C. Before DNA extraction, the thalli were inspected under a stereomicroscope and cleaned with sterile water to avoid contamination by other fungi and epiphytic microalgae. Fragments from different parts of each thallus were arbitrarily excised and pooled together, put into an Eppendorf tube, and ground using a pestle and 400 μL of lysis buffer of the DNeasy Plant Mini kit (Qiagen, Hilden, Germany). Total genomic DNA was extracted and purified using the above kit following the manufacturer’s instructions, and DNA was finally eluted in a volume of 50 μL.

2.2. PCR Amplification and Sequencing

The barcode nrITS region was selected to assess the taxonomic identity of the primary phycobiont as well as the mycobiont of the collected Ramalina thalli. The primers pairs used were nr-SSU-1780 [29] and ITS4T [30] for the phycobiont, and ITS1F [31] and ITS4a [32] for the mycobiont. PCR reactions were performed in a total volume of 25 μL using the EmeraldAmp GT PCR Master Mix (Takara, Shiga, Japan), which required the addition of the template DNA (1 μL), specific primers (1 μL each, 10 μM), and distilled water. The PCR amplification program comprised an initial denaturation step at 94 °C for 2 min, followed by 30 cycles at 94 °C for 30 s, 56 °C for 45 s, and 72 °C for 1 min, and a final elongation step at 72 °C for 5 min. Amplifications were carried out on 96-well lab cyclers, SensoQuest (Progen Scientific Ltd., Mexborough, UK). The PCR products were then electrophoresed in a 1% agarose gel and visualized using GelRed. The products were purified using the Gel Band Purification Kit (GE Healthcare Life Science, Piscataway, NJ, USA). The amplified PCR products were sequenced with an ABI 3730XL sequencer using the BigDye Terminator 3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA) at StabVida (Lisbon, Portugal). Raw electropherograms were manually checked, trimmed, and assembled using SeqmanII v.5.07© (DnaStar Inc., Madison, WI, USA). GenBank accession numbers are listed in bold in Table S1 (OR978676 to OR979068 for fungal nrITS and OR990615 to OR991088 for Trebouxia nrITS).

2.3. Inference of Genealogical Relationships among Phycobiont Haplotypes

A phycobiont haplotype network for the R. farinacea group was built using an alignment of nrITS sequences. The DNA dataset was assembled with the sequences obtained in the present study and those obtained from the specimens studied in Moya et al. [4] and a selection of 58 sequences from GenBank (available at http://www.ncbi.nlm.nih.gov/; accessed on 19 December 2023) (Table S1). The search criterion in GenBank was “Trebouxia AND Ramalina farinacea”, and only those accessions with accompanying information about sampling locality were considered [11,33,34,35]. The final sequence dataset consisted of 527 sequences, which were subsequently aligned with MAFFT v.7.308 [36,37] using the algorithm FFT-NS-I x1000, the 200PAM/k = 2 scoring matrix, a gap open penalty of 1.5, and an offset value of 0.123. The resulting nrITS alignment was manually optimized in Geneious v.9.0.2 to trim the ends of longer sequences that included fragments of 18S–26S ribosomal subunits. Genealogical relationships among haplotypes were then calculated under a statistical parsimony framework in PopART v.1.7 [38] using the TCS method under a 95% parsimony probability criterion, with gaps treated as a 5th character state [39,40,41]. Because the inference of haplotype networks is sensitive to ambiguous base calls and missing data [42], the alignment was first pruned and reduced to 415 complete sequences. The network was artistically edited in Adobe Illustrator 2023 and haplotypes were labelled according to the above-mentioned ten biogeographic subregions [27].

2.4. Interaction Networks

We built bipartite interaction networks among phycobiont haplotypes of the R. farinacea group and (a) different mycobionts, considering species circumscriptions in Moya et al. [4], (b) mycobiont genetic clusters without any consideration of species boundaries, and (c) biogeographic areas. Hierarchical Bayesian genetic clustering was performed using fastbaps v.1.0.8 [43] in R v.4.3.1 [44]. The fast BAPS algorithm is based on applying the hierarchical Bayesian clustering algorithm [45] to the problem of clustering genetic sequences using the same likelihood as BAPS [46]. We employed optimized.baps for the Dirichlet prior hyperparameter, as it generally outperforms other prior options [43], and the value of k.init was set to 20. This analysis was seeded with a mycobiont nrITS alignment that was built using newly obtained sequences (n = 395) and a selection from GenBank (n = 227) with the following search criteria: “Ramalina farinacea AND internal”, “Ramalina alisiosae AND internal”, and “Ramalina subfarinacea AND internal”. The final dataset consisted of 622 sequences, which were subsequently aligned with MAFFT v.7.308 [36,37] as described above. The three bipartite networks were constructed using the function plotweb in the R package bipartite [47].

2.5. Niche Hypervolumes

Because the two associated microalgae species corresponded to Trebouxia jamesii and T. lynnae (see Section 3 below), we represented their climatic niche, along with those of the four predominant haplotypes of T. jamesii, using the Hutchinsonian niche concept, defining a species’ niche as an n-dimensional hypervolume with dimensions corresponding to environmental variables [48]. The climatic hypervolumes were constructed using multivariate kernel density estimation [49]. A PCA based on 19 bioclimatic variables was conducted [50], and the first two axes, which explained 78% of the total variance, were selected to calculate hypervolumes for each species-level lineage and haplotype. The boundaries of the kernel density estimates were delineated by the probability threshold using the 0.85 quantile value. Hypervolume contours were plotted to project niche spaces based on 5000 random background points, using the alphahull contour type with the alpha smoothing value set to 0.55. For ease of comprehension, differences in climatic variables (BIO1–BIO19) between species were visualized by using box-plots. The significance of a difference between species was tested by using Wilcoxon tests, and Bonferroni correction was applied to adjust the p-value. All analyses were performed in R studio v.4.3.1 [44] using base functions and the packages hypervolume [49] and alphahull [51].

2.6. Trebouxia jamesii Geographical Occurrence

The GenBank sequence data of the nrITS were considered to evaluate the global distribution of Trebouxia jamesii and the taxonomical breadth of associated mycobionts. BLASTn searches [52] against that database using complete or partial nrITS sequences and following the search criterion “Trebouxia jamesii AND internal” were performed; consequently, 1154 hits with a 99–100% of nucleotide identity were retrieved, and their corresponding sequences were downloaded. The final dataset also included four T. jamesii haplotypes obtained in the present study and a representative sequence of a strain of T. jamesii (FJ626733 UTEX 2233) located at the Culture Collection of Algae of the University of Göttingen, Germany (SAG; available at https://sagdb.uni-goettingen.de/, accessed on 13 January 2024). Sequences of T. lynnae were also included in the dataset, and Asterochloris mediterranea was used to root the phylogenetic tree. Then, a phylogeny was built to explore within-species phylodiversity and to corroborate the validity of the taxonomic labels of GenBank sequences using an alignment obtained with MAFFT employing the same parameter setting as above. The online version of RAxML-HPC2 hosted at the CIPRES Science Gateway [53,54,55] was used to estimate phylogeny under a maximum likelihood (ML). The analysis used the GTRGAMMA substitution model and one thousand rapid bootstrap pseudoreplicates were conducted to evaluate nodal support. The resulting ML phylogenetic tree was visualized with the iTOL web tool [56], and Adobe Illustrator 2023 was used for artwork. Tree nodes with bootstrap support (BS) values equal to or greater than 70% were regarded as significantly supported. Based on phylogenetic clade clustering, only 211 sequences that matched with representative T. jamesii lineages were considered for assessing the geographic distribution and range of associated fungal partners of this microalga (Table S2) [11,12,33,35,57,58,59,60,61,62,63,64,65,66,67].

3. Results

3.1. Phycobiont Haplotype Network

The number of nrITS phycobiont haplotypes inferred from the newly sequenced specimens and the selection of phycobiont sequences from the R. farinacea group from GenBank was 52. The haplotype network showed the presence of two different Trebouxia species, T. jamesii and T. lynnae, separated by at least 28 different nucleotides (Figure 2). While T. lynnae occurred in the Macaronesian archipelagos and the European continent, T. jamesii was restricted to continental localities and archipelagos in the Mediterranean Basin. At the level of the 76 sampling localities considered in the present work, T. jamesii was the only phycobiont in 60 of them, whereas T. lynnae was the sole microalga found in eight localities in the Macaronesian archipelagos and the Iberian and Italian peninsulas (Table S3). In the eight remaining localities, both phycobiont species co-existed; the predominance of each microalga in these localities varied. Furthermore, considering the distance to the sea from each sampling locality in the straight shortest line (kilometres), T. jamesii occurred in both coastal and inland localities (up to 400 km), while T. lynnae preferred coastal localities, with a maximum inland occurrence at 58 km in the locality of El Toro in the eastern Iberian Peninsula (Figure S2).
The two most frequent T. jamesii haplotypes, H3 and H9, had 21 other haplotypes differing from them by a single nucleotide polymorphism (SNP) (Figure 2). Haplotypes H3 and H9 included specimens from all over the sampling area, except for the Macaronesian archipelagos, while most of the other 21 connected haplotypes were individually represented by a single biogeographical subregion. Haplotypes H1 and H2 were also occurred very frequently in nearly all considered localities. Central Europe and the Mediterranean Central Iberian Peninsula were the two areas hosting a higher number of exclusive T. jamesii haplotypes. On the other hand, T. lynnae co-occurred with T. jamesii in certain localities of the Iberian Peninsula, Balearic Islands, and Italy, which corresponded to the following biogeographic subregions: Atlantic Europe, Mediterranean Central Iberian Peninsula, Balearic-Catalonian-Provençal, and Italo-Thyrrenian. The latter two areas showed the highest number of exclusive minor haplotypes.

3.2. fastbaps Cluster Assignment and Interaction Networks

The four species forming the R. farinacea group sensu Moya et al. [4] showed differential associations to phycobiont species and haplotypes (Figure 3A). Ramalina farinacea, Ramalina sp., and R. subfarinacea were associated with both phycobionts, T. jamesii and T. lynnae. Only R. alisiosae was found to be exclusively associated with T. lynnae in Macaronesia (Figure 3A,C). Ramalina farinacea specimens (the most abundant in our assembled dataset) formed the highest number of associations with different T. jamesii haplotypes. From the phycobiont perspective, the two most-frequent T. jamesii haplotypes, H3 and H9, interacted with all mycobiont species except for R. alisiosae. The less frequent T. jamesii haplotypes (e.g., H6, H19, H20, etc.) only interacted with R. farinacea. At a geographical level, we found that both microalgae were not exclusive to a single biogeographic subregion (Figure 3C). Although T. jamesii was widespread in continental Europe and present in the Balearic Islands, this species was not detected in Macaronesian samples of the R. farinacea group. On the other hand, T. lynnae occurred predominantly in biogeographic areas that are often characterized by displaying milder climatic conditions, e.g., Macaronesia, the eastern Iberian Peninsula, and Balearic Islands, as well as localities in Atlantic Europe rather close to the coast. The Trebouxia jamesii haplotypes H1, H3, and H9 were shown to be widespread in subregions with markedly different climatic conditions, such as Boreal, Central Europe, Atlantic Europe, and the climatically milder Mediterranean Basin. H3 was also found in North California, which indicates a trans-Atlantic distribution.
Furthermore, the fast BAPS algorithm inferred eight clusters based on the mycobiont nrITS dataset. Figure 3B shows interactions of these clusters (each represented by a Roman number) with phycobiont nrITS haplotypes. For a better comprehension only clusters I, IV, VI, VII, and VIII were represented (Figure 3B). This was because either the remaining clusters derived from mycobiont GenBank nrITS sequences for which there was no information about the associated phycobiont, or they were associated with phycobiont haplotypes consisting of a single sequence/individual. Cluster VIII encompassed the bulk of R. farinacea samples and was associated with the highest number of different T. jamesii haplotypes. This cluster and the less frequent clusters I, IV, and VI were associated with both Trebouxia microalgae. On the opposite side, cluster VII (corresponding to the Macaronesian samples of R. alisiosae) was only associated with the three haplotypes of T. lynnae.

3.3. Climatic Niches

Two-dimensional hypervolumes for T. jamesii and T. lynnae and for the four most abundant T. jamesii haplotypes were built, and PC1 and PC2 explained 78% of the variation in climatic variables (Figure 4). Due to considerable differences in climatic traits among the Macaronesian archipelagos, Mediterranean, and Northern Europe, the hypervolumes of each species were split into two isolated areas, one broader than the other (Figure 4B,C). A discernible gradient of climatic niches along PC2 was observed, explaining variations in temperature and precipitation. Both species exhibited large niches overlapping (Figure 4B) in an area representative of localities placed in the Mediterranean Basin and temperate Europe (especially the Atlantic Europe and Central Europe biogeographic subregions). Trebouxia jamesii occurred in more diverse climate types than T. lynnae, including an isolated zone in the PCA representing the Boreal European region that shows a wide temperature seasonality. On the other hand, T. lynnae showed preference for warmer and drier climates (BIO1, 6, 16, 17, Wilcox tests p < 0.001 for all comparisons; Figure 5), and an isolated area related to the Macaronesian region (Figure 4B). The hypervolumes of the most frequent T. jamesii haplotypes H1, H2, H3, and H9 revealed rather negligible differences regarding their climatic niches (Figure 4C). H1 and H3 emerged as the most widespread in the analysis, with two distinct areas in the PCA, the smaller representing Boreal Europe.

3.4. Phycobiont Geographical Occurrence

A total of 220 sequences were clustered with T. jamesii (Table S2 and Figure S3). According to the information reported in GenBank accession metadata, T. jamesii has been detected as a primary phycobiont in the following lichen species: Amandinea punctata, Anaptychia runcinata, Biatora sp., Candelariella vitellina, Evernia prunastri, Lecanographa amylacea, Lecanora argentata, L. bicincta, L. chlarotera, L. frustulosa, L. glabrata, L. rupicola, L. sulphurea, Lecidea roseotincta, Lecidella elaeochroma, Lepra amara, Loxospora elatina, Melanelixia fuliginosa, M. glabratula, Myriolecis hagenii, Pertusaria amara, Pe. coccodes, Pe. leioplaca, Phlyctis argena, Protoparmelia badia, Pr. montagnei, Pr. psarophana, Ramalina calicaris, R. capitata, R. farinacea, R. fastigiata, R. fraxinea, R. lusitanica, R. pollinaria, R. pontica, Rhizocarpon geographicum, Rhizoplaca sp., Tephromela atra, T. grumosa, and Umbilicaria grisea. According to the analyzed GenBank metadata and our own dataset, the geographic distribution of T. jamesii would encompass continental Europe, the British Islands, and islands in the Mediterranean Basin (Austria, Crete, Crimean Peninsula, Czech Republic, Denmark, Estonia, Finland, France, Germany, Greece, Italy, Norway, Poland, Portugal, Slovakia, Slovenia, Spain, Sweden, United Kingdom, and Turkey); Peru in South America and the US (e.g., California) in North America; Algeria in northern Africa; and India in Asia.
Barreno et al. [68] provided a comprehensive analysis of the geographic distribution of T. lynnae (see Supplementary Table S2 in [68]), and according to their work and the results here obtained, T. lynnae associates as a primary phycobiont with Lecanographa amylacea, Protoparmelia montagnei, the Ramalina decipiens group, R. farinacea, R. fastigata, R. menziesii, and an undescribed Ramalina sp. that belongs to the R. farinacea group. Geographically, this microalga is distributed in continental Europe (Italy, Poland, Sweden, and Spain), the Balearic Islands, the Macaronesian archipelagos (Canary Islands and Cape Verde), the island of Madeira, North America, and New Zealand.

4. Discussion

Lichens are usually defined as paradigms of mutualistic symbioses, in which fungal mycobionts, populations of photosynthetic microalgae, and other microorganisms interact extracellularly and form morphologically and physiologically complex thalli (holobiomes) [69,70,71,72,73]. A question that still raises considerable curiosity among lichenologists is whether (and why) a given lichenized microalga can be replaced by equally compatible photobionts throughout the lichen’s distribution range [29,64,74,75,76,77,78,79]. This question has been already addressed at various taxonomic, geographic and ecological scales, generating valuable data about the evolutionary history of lichen symbioses [20,21,22,33,80,81,82,83], although, fungal–algal association patterns are still unknown for the vast majority of lichens [73]. In the present study, we explored patterns of symbiont associations in species of the R. farinacea group, which comprises R. subfarinacea, R. alisiosae, and Ramalina sp. [4], through an exhaustive European and Macaronesian sampling of specimens spanning ten biogeographic subregions with highly diverse climatic conditions.
Our results revealed that in the sampled regions, the most frequent phycobionts in the R. farinacea group were Trebouxia jamesii and T. lynnae, two phylogenetically closely related species [68]. Associations of closely related mycobionts with single or multiple phycobiont(s) which are also among themselves phylogenetically closely related have been already detected in other lichen symbioses, such as in the Macaronesian Ramalina decipiens group [22], the amphitropically distributed Pseudephebe spp. [84], or the Mediterranean and temperate epiphytic Parmelia sulcata and P. saxatilis groups [78]. A pattern is also observed at the level of the single mycobiont species in Ramalina menziesii in western North America [33] and the Mediterranean Seirophora villosa [85]. These findings do not necessarily mean that strict parallel diversification between lichenized fungi and their green algal partners took place [29,86,87]. Instead, it might suggest that fungal–algal cell recognition mechanisms are rather conserved and that only some interactions succeed during the first lichenization steps. Furthermore, the composition of the “pool” of locally available microalgae may exert a strong influence on the range of possible interactions as well [88,89,90,91]. Singh et al. [35] have found R. farinacea in association only with T. jamesii in the emblematic Białowieża Forest (Poland). This microalga was the most predominant in the epiphytic lichen community, being present in at least 53% of the sampled species, most of which were phylogenetically unrelated to R. farinacea and reproduced either sexually or asexually. Certainly, future research is needed to understand the composition of the “pool” of phycobionts in lichen communities thriving in Mediterranean regions, providing context for our results showing the high diversity of haplotypes and interactions in the R. farinacea group.
Insular vs. continental distribution of the phycobionts—We found that phycobiont association is not balanced amongst the studied mycobionts. The Canarian R. alisiosae showed a strict association with T. lynnae, whereas the other three mycobionts interacted with both T. jamesii and T. lynnae. Although it might be questioned that the results are biased by the low number of samples of R. alisiosae against those of R. farinacea, our results are well in line with those of Blázquez et al. [22], who have used high-throughput sequencing techniques to evaluate the phycobiont diversity in fifteen species of the Ramalina decipiens group in the Canary Islands, Madeira, and Cape Verde. Although the authors did not find remarkable differences in the identity of associated phycobionts, T. lynnae turned out to be the most frequent microalgal partner. The less frequent phycobionts were T. jamesii, T. aggregata, T. decolorans, T. australis, and some still undescribed Trebouxia species. The R. decipiens group represents a putative radiation of endemic lichenized fungi inhabiting the Macaronesian region, whose diversification does not seem to be linked with differing associations to phycobiont lineages [22]. The R. farinacea group likely diversified in the Macaronesian–Mediterranean region as well [4,9]. The strict interaction of R. alisiosae with T. lynnae is in line with the observation that this microalga is better adapted to the environment conditions of oceanic islands than other phycobionts [12,14,18,22]. This is further supported by the lack of T. jamesii lineages in the Macaronesian thalli of R. farinacea. However, the phylogeny of the R. farinacea group is still not well resolved [4], and due to the low number of studied R. alisiosae specimens we cannot assess whether phycobiont switch could have contributed to diversification in our target species.
The relationship of the lichen R. farinacea with T. jamesii breaks in the Canary Islands, Madeira, Sicily, and the Balearic Islands, where the mycobiont also (or only) associates with T. lynnae. del Campo et al. [11] and Molins et al. [12] have previously suggested a preference for continentality for T. jamesii and insularity for T. lynnae, which is further evidenced by our ecological niche analysis. The preference for warmer and more humid climates of T. lynnae aligns with the results of Casano et al. [14], who have demonstrated that this microalga has better photosynthetic performance at higher temperature and irradiance, while T. jamesii thrives at moderate temperature and irradiance. The better physiological performance of T. lynnae under more oxidative conditions than T. jamesii may reflect its greater capacity to undertake key metabolic adjustments, including increased non-photochemical quenching, higher antioxidant protection, and the induction of repair mechanisms [14]. These considerations align well with the previous physiological analyses performed on axenic isolates in vitro of the two phycobionts.
Furthermore, Hinojosa-Vidal et al. [92] and Pérez-Rodrigo et al. [93] have studied the effects of prolonged exposure to high salt concentrations on both T. lynnae and T. jamesii species and have demonstrated the extraordinarily higher tolerance to osmotic and saline stress of these two Trebouxia species compared to other microalga genera, such as Asterochloris and Chlorella. Their results suggest that the two Trebouxia phycobionts could cope with highly saline environments. Thus, a clear relationship was observed in the present study between T. lynnae occurrence and the distance to the sea, which might be related to its tolerance to salinity by water sprays and aerosols. Although T. jamesii also inhabits coastal areas (which could also be explained by its tolerance to salinity), it is able to colonize inland areas where T. lynnae does not spread into. Coastal zones usually display a warmer climate due to the isothermal influence of the air’s relative humidity, and temperatures tend to be more stable throughout the year. Contrary, inland zones experience more pronounced seasonal thermal variations. It is likely that a combined effect of the influence of temperature and distance from the sea controls the distribution of these microalgae. Even within thalli, the tolerances and physiological roles of photobionts vary; for example, two different photobionts provide different polysaccharides under changing seawater immersion in the marine lichen Lichina pygmaea [79]. Thus, it is a good example for testing hypotheses of how complex photobiont communities might react to environmental changes in other lichen symbioses, such as the cases studied in this work.
The multiplicity of haplotypes—Previous studies have already described a specific association of T. jamesii and T. lynnae with R. farinacea based only on Iberian Peninsula and Canary Islands specimens [11,12,14]. Here, we found that the association between R. farinacea and T. jamesii spanned from 36 to 63 degrees of latitude, i.e., from northern Africa (Algeria) to Central, Alpine, and Boreal regions in continental Europe. Some phycobiont haplotypes were restricted to localities in Central and Boreal Europe (e.g., H1, H6, H19, and H35), which might point towards fine-tuned fungal–algal interactions in these regions, assuming the existence of “pools” of potentially better-adapted microalgae [88,89,90,94]. In fact, phycobiont switches have been suggested for various epiphytic, terricolous, and saxicolous lichens throughout their distribution ranges, often aligning with notable changes in environmental conditions [18,20,61,81,84,85,94,95,96]. Alternatively, rare haplotypes might just represent the genetic differentiation in the phycobiont by intrathalline somatic mutations, as suggested for the microalgae occurring in the lichens Xanthoria parietina and Anaptychia ciliaris [97]. However, there were other T. jamesii haplotypes, such as H3 and H9, which showed a widespread distribution across different biogeographic regions. This finding is not strange, as identical nrITS phycobiont haplotypes have been reported in symbiosis with the same mycobiont over large geographic distances, even encompassing temperate, north, and south Polar Regions [84,94]. Given that R. farinacea reproduces mainly asexually using soredia, in which both symbionts are dispersed together, we might hypothesize that this lichen achieved its current broad geographic distribution by a successful first co-dispersion. Secondly, once the soredium developed and hyphae grew further, the mycobiont would have settled in that area by enwrapping those locally adapted phycobionts, i.e., those that were found to be less frequent and more geographically restricted haplotypes. This switch might have been achieved in at least two ways in R. farinacea: by the association of fungal hyphae projecting from soredia with a new phycobiont, as suggested by Nelsen and Gargas [86] in the strictly asexually reproducing Lepraria; alternatively, by a de novo lichenization processes in which case the mycobiont would reproduce sexually through the dispersion of meiotic spores. The number of different geographically adapted photobiont haplotypes may also be influenced by the species richness of the lichen communities present at a certain site. Events of photobiont switches are known to occur in lichens [17,23,29,65], especially if mycobionts and photobionts present low specificity and selectivity toward themselves, and thus many phycobionts are available for lichenization on growth substrates (i.e., wood, rocks, soil).
The present study has extended our knowledge about the geographic distribution and range of possible mycobiont partners of the two studied phycobionts. Trebouxia lynnae associates with lichen fungi belonging to phylogenetically distant families and from different areas of the planet, including continental Europe (Iberian Peninsula, Poland, and Sweden), the Balearic Islands, Macaronesian archipelagos (Canary Islands, Madeira, and Cape Verde), North America, and New Zealand [22,68]. Trebouxia jamesii is also a cosmopolitan phycobiont, which associates with several lichenized fungi [35,98].

5. Conclusions

In this research, we faced certain shortcomings when using the GenBank database to retrieve phycobionts sequences using nrITS BLASTn searches. Indeed, we found that out of the 1154 nrITS sequences labelled as T. jamesii in GenBank, only 220 corresponded to the authentic T. jamesii. This confusion might have stemmed from having formerly considered T. jamesii and T. simplex conspecific taxa due to the lack of morphological differences [99]. Nonetheless, differences between these two phycobionts were later found and the distinction of both species was finally accepted [63,89,100,101]. It is currently acknowledged that the authentic (type) strain of T. jamesii (UTEX 2233) belongs to the Trebouxia clade A [102], whereas T. simplex falls within the Trebouxia clade S. It is thus crucial to emphasize that species identification of Trebouxia cannot rely solely on BLASTn matches. In any study involving symbiotic microalgae, it is a necessary condition that identifications should be based on phylogenies encompassing species from all clades, including the 27 species officially described to date and indicated as references for phylogenetic datasets [103].
Future studies must reveal whether the apparently strict association of mycobiont members in the R. farinacea group with T. jamesii and T. lynnae is maintained in other areas of the planet, especially in North America. It would also be essential to examine whether the different haplotypes of T. jamesii found across its European distribution perform differently at the physiological level, and if the differences are related to varying environmental conditions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof10030206/s1. Figure S1: Biogeographic map of Europe [27], Figure S2: Bubble chart indicating the distance to the sea (in Km) from each sampling locality. The series coloured in light green represents Trebouxia lynnae, whereas T. jamesii series is coloured in dark green. The size of each bubble is proportional to the number of lichen thalli (n) associated with a given Trebouxia. The increasing scale of n, represented by the smaller bubble to the bigger one, is displayed above the bubble chart, Figure S3: Single circular phylograms estimated with RAxML showing the relationships among nrITS sequences downloaded from GenBank (n = 1154) to evaluate the global distribution of Trebouxia jamesii. The dataset included four T. jamesii haplotypes obtained in the present study and a representative sequence of a strain of T. jamesii located at the Culture Collection of Algae of the University of Göttingen, Germany (SAG). Sequences of T. lynnae were also included in the dataset and Asterochloris mediterranea was used to root phylogenetic trees; for the sake of clarity, tree representation ignored branch lengths; minute green dots on nodes denote nodal support (BS ≥ 70%); Table S1: Specimens analyzed in this study with GenBank accession numbers (in bold) of fungal and algal nrITS and details of the collection location. Including newly sequences generated in this study, those obtained from specimens studied in Moya et al. [4] and a selection of 58 sequences from the GenBank selected under the search criterion “Trebouxia AND Ramalina farinacea”, and only those accessions with accompanying information about sampling locality were considered, Table S2: GenBank accession numbers, host, locality, and reference for the Trebouxia species that matched with Trebouxia jamesii in GenBank, Table S3: Trebouxia composition, number of sequence and distance to the sea in a straight line (kilometres) from each sampling locality included in this study.

Author Contributions

Conceptualization, P.M., I.G.-B., S.C. and E.B.; methodology, I.G.-B., P.M. and S.C.; software, I.G.-B., P.M. and T.P.; sampling, P.M., I.G-B., S.C. and L.M.; validation, I.G.-B., P.M. and S.C.; formal analysis, I.G.-B., P.M. and S.C.; investigation, P.M., I.G.-B., S.C. and E.B.; resources, E.B., P.C. and I.G.-B.; writing—original draft preparation, P.M. and I.G.-B.; writing—review and editing, P.M., I.G-B., S.C., L.M., E.B., P.C. and T.P.; project administration, P.C., I.G.-B. and E.B.; funding acquisition, E.B., P.C. and I.G.-B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by PID2021-127087NB-I00 (Spanish Ministry of Science and Innovation to I.G.-B. and P.C.); PROMETEO/2021/005 (Prometeo Excellence Research Program, Generalitat Valenciana, Spain, to E.B. and P.C.); and a postdoctoral contract (Next generation EU, MS21-058) by Ministerio de Universidades—Spain (S.C.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The dataset generated during the current study is available in the GenBank (see Table S1): OR978676 to OR979068 for fungal nrITS and OR990615 to OR991088 for Trebouxia nrITS.

Acknowledgments

We would like to thank Violeta Atienza (UV), Wolfgang von Brackel (Röttenbach, GE), Per-Anders Essen (UMU, SE), Francisco Gasulla (UAH), Tarek Hamel (UBMA, Annaba), Carlos Lobo (MADJ, Madeira), Kristiina Mark (EMU, Tartu), Arantxa Molins (INAGEA), Leena Myllis (LUOMUS, FI), Ondrej Peksa (ZCM, CZ), Sergio Prats (Málaga), Domenico Puntillo (UniCal, IT), Arnoldo Santos (Tenerife), Gabriel Segarra (UV), Jana Steinová (CUNI, CZ), Tor Tonsberg (UIB, NO), Martin Westberg (UU, SE), Juan Carlos Zamora (CJB, CH), for contributing to the collection, andLucie Vančurová for the ecological niche script.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Acharius, E. Lichenographia Universalis; Apud Iust. Frid. Danckwerts: Göttingen, Germany, 1810. [Google Scholar]
  2. Armstrong, R.A. Dispersal of soredia from individual soralia of the lichen Hypogymnia physodes (L.) Nyl. in a simple wind tunnel. Environ. Exp. Bot. 1994, 34, 39–45. [Google Scholar] [CrossRef]
  3. Walser, J.C.; Zoller, S.; Büchler, U.; Scheidegger, C. Species-specific detection of Lobaria pulmonaria (lichenized ascomycete) diaspores in litter samples trapped in snow cover. Mol. Ecol. 2001, 10, 2129–2138. [Google Scholar] [CrossRef]
  4. Moya, P.; Garrido-Benavent, I.; Chiva, S.; Pérez-Ortega, S.; Blázquez, M.; Pazos, T.; Hamel, T.; Myllys, L.; Tønsberg, T.; Esseen, P.-A.; et al. Phylogeography of Ramalina farinacea (Lichenized Fungi, Ascomycota) in the Mediterranean Basin, Europe, and Macaronesia. Diversity 2023, 15, 310. [Google Scholar] [CrossRef]
  5. del Campo, E.; Gimeno, J.; Casano, L.M.; Gasulla, F.; García-Breijo, F.; Reig-Armiñana, J.; Barreno, E. South European populations of Ramalina farinacea (L.) Ach. share different Trebouxia algae. Bibl. Lichenol. 2010, 105, 247–256. [Google Scholar]
  6. Stocker-Wörgötter, E.; Elix, J.A.; Grube, M. Secondary chemistry of lichen-forming fungi: Chemosyndromic variation and DNA-analyses of cultures and chemotypes in the Ramalina farinacea complex. Bryologist 2004, 107, 152–162. [Google Scholar] [CrossRef]
  7. Nylander, W. Observata lichenologia in Pyranaeis orientalibus. Bull. Soc. Linnéenne Normandie 1872, 6, 256–328. [Google Scholar]
  8. Pérez-Vargas, I.; Pérez-Ortega, S. A new endemic Ramalina species from the Canary Islands (Ascomycota, Lecanorales). Phytotaxa 2014, 159, 269–278. [Google Scholar] [CrossRef]
  9. Krog, H.; Østhagen, H. The genus Ramalina in the Canary Islands. Norw. J. Bot. 1980, 27, 255–296. [Google Scholar]
  10. Aptroot, A.; Schumm, F. Key to Ramalina species known from Atlantic islands, with two new species from the Azores. Sauteria 2008, 15, 21–57. [Google Scholar]
  11. del Campo, E.M.; Català, S.; Casano, L.M.; Gimeno, J.; del Hoyo, A.; Martínez-Alberola, F.; Casano, L.M.; Grube, M.; Barreno, E. The genetic structure of the cosmopolitan three-partner lichen Ramalina farinacea evidences the concerted diversification of symbionts. FEMS Microbiol. Ecol. 2013, 83, 310–323. [Google Scholar] [CrossRef] [PubMed]
  12. Molins, A.; Moya, P.; Muggia, L.; Barreno, E. Thallus growth stage and geographic origin shape microalgal diversity in Ramalina farinacea lichen holobionts. J. Phycol. 2021, 57, 975–987. [Google Scholar] [CrossRef]
  13. Moya, P.; Molins, A.; Martínez-Alberola, F.; Muggia, L.; Barreno, E. Unexpected associated microalgal diversity in the lichen Ramalina farinacea is uncovered by pyrosequencing analyses. PLoS ONE 2017, 12, e0175091. [Google Scholar] [CrossRef]
  14. Casano, L.M.; del Campo, E.M.; García-Breijo, F.J.; Reig-Armiñana, J.; Gasulla, F.; Del Hoyo, A.; Guéra, A.; Barreno, E. Two Trebouxia algae with different physiological performances are ever-present in lichen thalli of Ramalina farinacea. Coexistence versus competition? Environ. Microbiol. 2011, 13, 806–818. [Google Scholar] [CrossRef]
  15. del Hoyo, A.; Álvarez, R.; del Campo, E.; Gasulla, F.; Barreno, E.; Casano, L.M. Oxidative stress induces distinct physiological responses in the two Trebouxia phycobionts of the lichen Ramalina farinacea. Ann. Bot. 2011, 107, 109–118. [Google Scholar] [CrossRef]
  16. Muggia, L.; Vancurová, L.; Škaloud, P.; Peksa, O.; Wedin, M.; Grube, M. The symbiotic playground of lichen thalli–a highly flexible photobiont association in rock-inhabiting lichens. FEMS Microbiol. Ecol. 2013, 85, 313–323. [Google Scholar] [CrossRef]
  17. Ohmura, Y.; Takeshita, S.; Kawachi, M. Photobiont diversity within populations of a vegetatively reproducing lichen, Parmotrema tinctorum, can be generated by photobiont switching. Symbiosis 2019, 77, 59–72. [Google Scholar] [CrossRef]
  18. Peksa, O.; Škaloud, P. Do photobionts influence the ecology of lichens? A case study of environmental preferences in symbiotic green alga Asterochloris (Trebouxiophyceae). Mol. Ecol. 2011, 20, 3936–3948. [Google Scholar] [CrossRef]
  19. Muggia, L.; Pérez-Ortega, S.; Fryday, A.; Spribille, T.; Grube, M. Global assessment of genetic variation and phenotypic plasticity in the lichen-forming species Tephromela atra. Fungal Divers. 2014, 64, 233–251. [Google Scholar] [CrossRef]
  20. Dal Grande, F.; Rolshausen, G.; Divakar, P.K.; Crespo, A.; Otte, J.; Schleuning, M.; Schmitt, I. Environment and host identity structure communities of green algal symbionts in lichens. New Phytol. 2018, 217, 277–289. [Google Scholar] [CrossRef] [PubMed]
  21. Vančurová, L.; Muggia, L.; Peksa, O.; Řídká, T.; Škaloud, P. The complexity of symbiotic interactions influences the ecological amplitude of the host: A case study in Stereocaulon (lichenized ascomycota). Mol. Ecol. 2018, 27, 3016–3033. [Google Scholar] [CrossRef] [PubMed]
  22. Blázquez, M.; Hernández-Moreno, L.S.; Gasulla, F.; Pérez-Vargas, I.; Pérez-Ortega, S. The role of photobionts as drivers of diversification in an island radiation of lichen-forming fungi. Front. Microbiol. 2022, 12, 784182. [Google Scholar] [CrossRef] [PubMed]
  23. Škvorová, Z.; Černajová, I.; Steinová, J.; Peksa, O.; Moya, P.; Škaloud, P. Promiscuity in lichens follows clear rules: Partner switching in Cladonia is regulated by climatic factors and soil chemistry. Front. Microbiol. 2022, 12, 781585. [Google Scholar] [CrossRef] [PubMed]
  24. Steinová, J.; Holien, H.; Košuthová, A.; Škaloud, P. An exception to the rule? Could photobiont identity be a better predictor of lichen phenotype than mycobiont identity? J. Fungi 2022, 8, 275. [Google Scholar] [CrossRef] [PubMed]
  25. Garrido-Benavent, I.; Mora-Rodríguez, M.R.; Chiva, S.; Fos, S.; Barreno, E. Punctelia borreri and P. subrudecta (Parmeliaceae) associate with a partially overlapping pool of Trebouxia gelatinosa lineages. Lichenologist 2023, 55, 389–399. [Google Scholar] [CrossRef]
  26. Rodríguez-Arribas, C.; Martínez, I.; Aragón, G.; Zamorano-Elgueta, C.; Cavieres, L.; Prieto, M. Specialization patterns in symbiotic associations: A community perspective over spatial scales. Ecol. Evol. 2023, 13, e10296. [Google Scholar] [CrossRef] [PubMed]
  27. Rivas-Martínez, S.; Penas, Á.; del Río, S.; Díaz González, T.E.; Rivas-Sáenz, S. Bioclimatology of the Iberian Peninsula and the Balearic Islands. In The Vegetation of the Iberian Peninsula; Loidi, J., Ed.; Springer: Cham, Switzerland, 2017; Volume I, pp. 29–80. [Google Scholar]
  28. Wickham, H. Data analysis. In ggplot2; Springer: Cham, Switzerland, 2016; pp. 189–201. [Google Scholar]
  29. Piercey-Normore, M.D.; DePriest, P.T. Algal switching among lichen symbioses. Am. J. Bot. 2001, 88, 1490–1498. [Google Scholar] [CrossRef] [PubMed]
  30. Kroken, S.; Taylor, J.W. Phylogenetic species, reproductive mode, and specificity of the green alga Trebouxia forming lichens with the fungal genus Letharia. Bryologist 2000, 103, 645–660. [Google Scholar] [CrossRef]
  31. Gardes, M.; Bruns, T.D. ITS primers with enhanced specificity for basidiomycetes-application to the identification of mycorrhizae and rusts. Mol. Ecol. 1993, 2, 113–118. [Google Scholar] [CrossRef]
  32. Larena, I.; Salazar, O.; González, V.; Julián, M.C.; Rubio, V. Design of a primer for ribosomal DNA internal transcribed spacer with enhanced specificity for ascomycetes. J. Biotechnol. 1999, 75, 187–194. [Google Scholar] [CrossRef]
  33. Werth, S.; Sork, V.L. Ecological specialization in Trebouxia (Trebouxiophyceae) photobionts of Ramalina menziesii (Ramalinaceae) across six range-covering ecoregions of western North America. Am. J. Bot. 2014, 101, 1127–1140. [Google Scholar] [CrossRef]
  34. Álvarez, R.; Del Hoyo, A.; Díaz-Rodríguez, C.; Coello, A.J.; del Campo, E.M.; Barreno, E.; Casano, L.M. Lichen rehydration in heavy metal-polluted environments: Pb modulates the oxidative response of both Ramalina farinacea thalli and its isolated microalgae. Microbial. Ecol. 2015, 69, 698–709. [Google Scholar] [CrossRef] [PubMed]
  35. Singh, G.; Kukwa, M.; Dal Grande, F.; Łubek, A.; Otte, J.; Schmitt, I. A glimpse into genetic diversity and symbiont interaction patterns in lichen communities from areas with different disturbance histories in Białowieża forest, Poland. Microorganisms 2019, 7, 335. [Google Scholar] [CrossRef] [PubMed]
  36. Katoh, K.; Misawa, K.; Kuma, K.I.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [PubMed]
  37. Katoh, K.; Standley, D.M. Mafft multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed]
  38. Leigh, J.W.; Bryant, D. PopART: Full feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  39. Templeton, A.R.; Crandall, K.A.; Sing, C.F. A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics 1992, 132, 619–633. [Google Scholar] [CrossRef]
  40. Clement, M.; Snell, Q.; Walker, P.; Posada, D.; Crandall, K. TCS: A computer program to estimate gene genealogies. Mol. Ecol. 2000, 9, 1657–1660. [Google Scholar] [CrossRef] [PubMed]
  41. Shi, W.; Zhu, C. Coding and applying gaps in phylogenetic analyses. Acta Zootax Sin. 2007, 32, 1–5. [Google Scholar]
  42. Joly, S.; Stevens, M.I.; van Vuuren, B.J. Haplotype Networks Can Be Misleading in the Presence of Missing Data. Syst. Biol. 2007, 56, 857–862. [Google Scholar] [CrossRef]
  43. Tonkin-Hill, G.; Lees, J.A.; Bentley, S.D.; Frost, S.D.; Corander, J. Fast hierarchical Bayesian analysis of population structure. Nucleic Acids Res. 2019, 47, 5539–5549. [Google Scholar] [CrossRef]
  44. R Studio Team. R Studio: Integrated Development for R; R Studio; PBC: Boston, MA, USA, 2023; Available online: http://www.rstudio.com/ (accessed on 17 January 2024).
  45. Heller, K.A.; Ghahramani, Z. Bayesian Hierarchical Clustering. In Proceedings of the 22nd International Conference on Machine Learning, ICML ’05, Bonn, Germany, 7–11 August 2005; ACM: New York, NY, USA, 2005; pp. 297–304. [Google Scholar]
  46. Cheng, L.; Connor, T.R.; Sirén, J.; Aanensen, D.M.; Corander, J. Hierarchical and spatially explicit clustering of DNA sequences with BAPS software. Mol. Biol. Ecol. 2013, 30, 1224–1228. [Google Scholar] [CrossRef] [PubMed]
  47. Dormann, C.F.; Gruber, B.; Fründ, J. Introducing the bipartite package: Analysing ecological networks. R News 2008, 8, 8–11. [Google Scholar]
  48. Hutchinson, G.E. Concluding remarks. Cold Spring Harb. Symp. Quant. Biol. 1957, 22, 415–427. [Google Scholar] [CrossRef]
  49. Blonder, B.; Lamanna, C.; Violle, C.; Enquist, B.J. The n-dimensional hypervolume. Glob. Ecol. Biogeogr. 2014, 23, 595–609. [Google Scholar] [CrossRef]
  50. Karger, D.N.; Conrad, O.; Böhner, J.; Kawohl, T.; Kreft, H.; Soria-Auza, R.W.; Zimmermann, N.E.; Linder, H.P.; Kessler, M. Climatologies at high resolution for the earth’s land surface areas. Sci. Data 2017, 4, 170122. [Google Scholar] [CrossRef] [PubMed]
  51. Pateiro-Lopez, B.; Rodriguez-Casal, A. Alphahull: Generalization of the Convex Hull of a Sample of Points in the Plane Manual. 2016. Available online: https://cran.r-project.org/package=alphahull (accessed on 13 January 2024).
  52. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef] [PubMed]
  53. Stamatakis, A. RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006, 22, 2688–2690. [Google Scholar] [CrossRef]
  54. Stamatakis, A.; Hoover, P.; Rougemont, J. A rapid boostrap algorithm for the RAxML web server. Syst. Biol. 2008, 57, 758–771. [Google Scholar] [CrossRef]
  55. Miller, M.A.; Pfeiffer, W.; Schwartz, T. The CIPRES science gateway: A community resource for phylogenetic analyses. In Proceedings of the Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010. [Google Scholar]
  56. Letunic, I.; Bork, P. Interactive Tree Of Life (iTOL) v5: An online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef]
  57. Helms, G.; Friedl, T.; Rambold, G.; Mayrhofer, H. Identification of photobionts from the lichen family Physciaceae using algal-specific ITS rDNA sequencing. Lichenologist 2001, 33, 73–86. [Google Scholar] [CrossRef]
  58. Cordeiro, L.M.; Reis, R.A.; Cruz, L.M.; Stocker-Wörgötter, E.; Grube, M.; Iacomini, M. Molecular studies of photobionts of selected lichens from the coastal vegetation of Brazil. FEMS Microbiol. Ecol. 2005, 54, 381–390. [Google Scholar] [CrossRef]
  59. Blaha, J.; Baloch, E.; Grube, M. High photobiont diversity associated with the euryoecious lichen-forming ascomycete Lecanora rupicola (Lecanoraceae, Ascomycota). Biol. J. Linnean Soc. 2006, 88, 283–293. [Google Scholar] [CrossRef]
  60. Schmull, M.; Miadlikowska, J.; Pelzer, M.; Stocker-Wörgötter, E.; Hofstetter, V.; Fraker, E.; Hodkinson, B.P.; Reeb, V.; Kukwa, M.; Lumbsch, H.T.; et al. Phylogenetic affiliations of members of the heterogeneous lichen-forming fungi of the genus Lecidea sensu Zahlbruckner (Lecanoromycetes, Ascomycota). Mycologia 2011, 103, 983–1003. [Google Scholar] [CrossRef] [PubMed]
  61. Muggia, L.; Pérez-Ortega, S.; Kopun, T.; Zellnig, G.; Grube, M. Photobiont selectivity leads to ecological tolerance and evolutionary divergence in a polymorphic complex of lichenized fungi. Ann. Bot. 2014, 114, 463–475. [Google Scholar] [CrossRef] [PubMed]
  62. Catalá, S.; Del Campo, E.M.; Barreno, E.; García-Breijo, F.J.; Reig-Armiñana, J.; Casano, L.M. Coordinated ultrastructural and phylogenomic analyses shed light on the hidden phycobiont diversity of Trebouxia microalgae in Ramalina fraxinea. Mol. Phyl. Evol. 2016, 94, 765–777. [Google Scholar] [CrossRef] [PubMed]
  63. Voytsekhovich, A.; Beck, A. Lichen photobionts of the rocky outcrops of Karadag massif (Crimean Peninsula). Symbiosis 2016, 68, 9–24. [Google Scholar] [CrossRef]
  64. Singh, G.; Dal Grande, F.; Divakar, P.K.; Otte, J.; Crespo, A.; Schmitt, I. Fungal–algal association patterns in lichen symbiosis linked to macroclimate. New Phytol. 2017, 214, 317–329. [Google Scholar] [CrossRef]
  65. Ertz, D.; Guzow-Krzemińska, B.; Thor, G.; Łubek, A.; Kukwa, M. Photobiont switching causes changes in the reproduction strategy and phenotypic dimorphism in the Arthoniomycetes. Sci. Rep. 2018, 8, 4952. [Google Scholar] [CrossRef]
  66. Mark, K.; Laanisto, L.; Bueno, C.G.; Niinemets, Ü.; Keller, C.; Scheidegger, C. Contrasting co-occurrence patterns of photobiont and cystobasidiomycete yeast associated with common epiphytic lichen species. New Phytol. 2020, 227, 1362–1375. [Google Scholar] [CrossRef]
  67. Bordenave, C.D.; Muggia, L.; Chiva, S.; Leavitt, S.D.; Carrasco, P.; Barreno, E. Chloroplast morphology and pyrenoid ultrastructural analyses reappraise the diversity of the lichen phycobiont genus Trebouxia (Chlorophyta). Algal. Res. 2022, 61, 102561. [Google Scholar] [CrossRef]
  68. Barreno, E.; Muggia, L.; Chiva, S.; Molins, A.; Bordenave, C.; García-Breijo, F.; Moya, P. Trebouxia lynnae sp. nov. (Former Trebouxia sp. TR9): Biology and Biogeography of an Epitome Lichen Symbiotic Microalga. Biology 2022, 11, 1196. [Google Scholar] [CrossRef]
  69. Guerrero, R.; Margulis, L.; Berlanga, M. Symbiogenesis: The Holobiont as a Unit of Evolution. Int. Microbiol. 2013, 16, 133–143. [Google Scholar]
  70. Hawksworth, D.L.; Grube, M. Lichens redefined as complex ecosystems. New Phytol. 2020, 227, 1281–1283. [Google Scholar] [CrossRef] [PubMed]
  71. Allen, J.L.; Lendemer, J.C. A call to reconceptualize lichen symbioses. Trends Ecol. Evol. 2022, 37, 582–589. [Google Scholar] [CrossRef] [PubMed]
  72. Kranner, I.; Pichler, G.; Grube, M. The lichen market place. New Phytol. 2022, 234, 541–543. [Google Scholar] [CrossRef]
  73. Grube, M. Lichens. In Fungal Associations. The Mycota; Hsueh, Y.P., Blackwell, M., Eds.; Springer: Cham, Switzerlands, 2024; Volume 9. [Google Scholar]
  74. Werth, S. Fungal-algal interactions in Ramalina menziesii and its associated epiphytic lichen community. Lichenologist 2012, 44, 543–560. [Google Scholar] [CrossRef]
  75. Moya, P.; Škaloud, P.; Chiva, S.; García-Breijo, F.J.; Reig-Arminana, J.; Vančurová, L.; Barreno, E. Molecular phylogeny and ultrastructure of the lichen microalga Asterochloris mediterranea sp. nov. from Mediterranean and Canary Islands ecosystems. Int. J. Syst. Evol. Microbiol. 2015, 65, 1838–1854. [Google Scholar] [CrossRef] [PubMed]
  76. Moya, P.; Chiva, S.; Molins, A.; Jadrná, I.; Škaloud, P.; Peksa, O.; Barreno, E. Myrmecia israeliensis as the primary symbiotic microalga in squamulose lichens growing in European and Canary Island terricolous communities. Fottea 2018, 18, 72–85. [Google Scholar] [CrossRef]
  77. Moya, P.; Molins, A.; Chiva, S.; Bastida, J.; Barreno, E. Symbiotic microalgal diversity within lichenicolous lichens and crustose hosts on Iberian Peninsula gypsum biocrusts. Sci. Rep. 2020, 10, 14060. [Google Scholar] [CrossRef]
  78. Moya, P.; Molins, A.; Škaloud, P.; Divakar, P.K.; Chiva, S.; Dumitru, C.; Molina, M.C.; Crespo, A.; Barreno, E. Biodiversity patterns and ecological preferences of the photobionts associated with the lichen-forming genus Parmelia. Front. Microbiol. 2021, 12, 765310. [Google Scholar] [CrossRef]
  79. Chrismas, N.; Tindall-Jones, B.; Jenkins, H.; Harley, J.; Bird, K.; Cunliffe, M. Metatranscriptomics reveals diversity of symbiotic interaction and mechanisms of carbon exchange in the marine cyanolichen Lichina pygmaea. New Phytol. 2023, 241, 2243–2257. [Google Scholar] [CrossRef]
  80. Helms, G. Taxonomy and Symbiosis in Associations of Physciaceae and Trebouxia. Ph.D. Thesis, University of Göttingen, Göttingen, Germany, 2003. [Google Scholar]
  81. Yahr, R.; Vilgalys, R.; Depriest, P.T. Strong fungal specificity and selectivity for algal symbionts in Florida scrub Cladonia lichens. Mol. Ecol. 2004, 13, 3367–3378. [Google Scholar] [CrossRef]
  82. Castillo, R.V.; Beck, A. Photobiont selectivity and specificity in Caloplaca species in a fog-induced community in the Atacama Desert, northern Chile. Fung Biol. 2012, 116, 665–676. [Google Scholar] [CrossRef]
  83. Nyati, S.; Scherrer, S.; Werth, S.; Honegger, R. Green-algal photobiont diversity (Trebouxia spp.) in representatives of Teloschistaceae (Lecanoromycetes, lichen-forming ascomycetes). Lichenologist 2014, 46, 189–212. [Google Scholar] [CrossRef]
  84. Garrido-Benavent, I.; Pérez-Ortega, S.; de los Ríos, A.; Fernández-Mendoza, F. Amphitropical variation of the algal partners of Pseudephebe (Parmeliaceae, lichenized fungi). Symbiosis 2020, 82, 35–48. [Google Scholar] [CrossRef]
  85. Garrido-Benavent, I.; Molins, A.; Barreno, E. Genetic variation in the symbiont partners in the endangered macrolichen Seirophora villosa (Teloschistaceae: Ascomycota). Bot. J. Linnean Soc. 2022, 199, 816–829. [Google Scholar] [CrossRef]
  86. Nelsen, M.P.; Gargas, A. Dissociation and horizontal transmission of codispersing lichen symbionts in the genus Lepraria (Lecanorales: Stereocaulaceae). New Phytol. 2008, 177, 264–275. [Google Scholar] [CrossRef]
  87. Pino-Bodas, R.; Stenroos, S. Global biodiversity patterns of the photobionts associated with the genus Cladonia (Lecanorales, Ascomycota). Microb. Ecol. 2021, 82, 173–187. [Google Scholar] [CrossRef]
  88. Beck, A. Photobiont inventory of a lichen community growing on heavy-metal-rich rock. Lichenologist 1999, 31, 501–510. [Google Scholar] [CrossRef]
  89. Beck, A. Selektivität der Symbionten Schwermetalltoleranter Flechten. Ph.D. Thesis, Ludwig-Maximilians-Universität München, Munich, Germany, 2002. [Google Scholar]
  90. Doering, M.; Piercey-Normore, M.D. Genetically divergent algae shape an epiphytic lichen community on Jack Pine in Manitoba. Lichenologist 2009, 41, 69–80. [Google Scholar] [CrossRef]
  91. Berlinches de Gea, A.; Verdú, M.; Villar-de Pablo, M.; Pérez-Ortega, S. The combined effects of habitat fragmentation and life history traits on specialisation in lichen symbioses. J. Ecol. 2023, 112, 200–216. [Google Scholar] [CrossRef]
  92. Hinojosa-Vidal, E.; Marco, F.; Martínez-Alberola, F.; Escaray, F.J.; García-Breijo, F.J.; Reig-Armiñana, J.; Carrasco, P.; Barreno, E. Characterization of the responses to saline stress in the symbiotic green microalga Trebouxia sp. TR9. Planta 2018, 248, 1473–1486. [Google Scholar] [CrossRef]
  93. Pérez-Rodrigo, M.; Moya, P.; Marco, F.; Carrasco, P.; Barreno, E.A. Symbiotic Trebouxia sp. TR9, Asterochloris erici and free-living Chlorella vulgaris green microalgae respond differentially to osmotic and saline stresses. In Proceedings of the IAL9 Program & Abstracts Book, Universidade Federal de Sergipe, Brasil, International Association for Lichenology 9th Symposium “Unlocking the Inner Lichen” Virtual-IAL9, Bonito, Brazil, 1–6 August 2021. [Google Scholar]
  94. Fernández-Mendoza, F.; Domaschke, S.; García, M.A.; Jordan, P.; Martín, M.P.; Printzen, C. Population structure of mycobionts and photobionts of the widespread lichen Cetraria aculeata. Mol. Ecol. 2011, 20, 1208–1232. [Google Scholar] [CrossRef]
  95. Yahr, R.; Vilgalys, R.; DePriest, P.T. Geographic variation in algal partners of Cladonia subtenuis (Cladoniaceae) highlights the dynamic nature of a lichen symbiosis. New Phytol. 2006, 171, 847–860. [Google Scholar] [CrossRef]
  96. Piercey-Normore, M.D. The lichen-forming ascomycete Evernia mesomorpha associates with multiple genotypes of Trebouxia jamesii. Phytologist 2006, 169, 331–344. [Google Scholar] [CrossRef]
  97. Dal Grande, F.; Alors, D.; Divakar, P.K.; Bálint, M.; Crespo, A.; Schmitt, I. Insights into intrathalline genetic diversity of the cosmopolitan lichen symbiotic green alga Trebouxia decolorans Ahmadjian using microsatellite markers. Mol. Phyl. Evol. 2014, 72, 54–60. [Google Scholar] [CrossRef]
  98. Sanders, W.B.; Masumoto, H. Lichen algae: The photosynthetic partners in lichen symbioses. Lichenologist 2021, 53, 347–393. [Google Scholar] [CrossRef]
  99. Friedl, T. Comparative ultrastructure of pyrenoids in Trebouxia (Microthamniales, Chlorophyta). Plant Syst. Evol. 1989, 164, 145–159. [Google Scholar] [CrossRef]
  100. Peršoh, D.; Beck, A.; Rambold, G. The distribution of ascus types and photobiontal selection in Lecanoromycetes (Ascomycota) against the background of a revised SSU nrDNA phylogeny. Mycol. Prog. 2004, 3, 103–121. [Google Scholar] [CrossRef]
  101. Opanowicz, M.; Grube, M. Photobiont genetic variation in Flavocetraria nivalis from Poland (Parmeliaceae, lichenized Ascomycota). Lichenologist 2004, 36, 125–131. [Google Scholar] [CrossRef]
  102. Muggia, L.; Nelsen, M.P.; Kirika, P.M.; Barreno, E.; Beck, A.; Lindgren, H.; Lumbsch, H.T.; Leavitt, S.D.; Trebouxia Working Group. Formally described species woefully underrepresent phylogenetic diversity in the common lichen photobiont genus Trebouxia (Trebouxiophyceae, Chlorophyta): An impetus for developing an integrated taxonomy. Mol. Phyl. Evol. 2020, 149, 106821. [Google Scholar] [CrossRef] [PubMed]
  103. Guiry, M.D.; Guiry, G.M. AlgaeBase. World-Wide Electronic Publication, National University of Ireland, Galway. 2024. Available online: https://www.algaebase.org (accessed on 17 January 2024).
Figure 1. (A) An epiphytic thallus of Ramalina farinacea collected in Bocairent (eastern Spain; photo: IGB) showing soredia (yellow asterisks) and apothecia (red asterisk). (B) Georeferenced map showing sampling localities of specimens of the “R. farinacea group” in the ten studied biogeographic subregions: Boreal (n = 25), Atlantic Europe (n = 56), Central (n = 51) Europe, Alpine (n = 6), Mediterranean West (n = 11), Central (n = 167) Iberian Peninsula, Balearic-Catalonian-Provençal (n = 35), Italo-Thyrrenian (n = 52), Adriatic (n = 14), and Macaronesian (n = 6). The map was created using the function map_data in the R package ggplot2 [28].
Figure 1. (A) An epiphytic thallus of Ramalina farinacea collected in Bocairent (eastern Spain; photo: IGB) showing soredia (yellow asterisks) and apothecia (red asterisk). (B) Georeferenced map showing sampling localities of specimens of the “R. farinacea group” in the ten studied biogeographic subregions: Boreal (n = 25), Atlantic Europe (n = 56), Central (n = 51) Europe, Alpine (n = 6), Mediterranean West (n = 11), Central (n = 167) Iberian Peninsula, Balearic-Catalonian-Provençal (n = 35), Italo-Thyrrenian (n = 52), Adriatic (n = 14), and Macaronesian (n = 6). The map was created using the function map_data in the R package ggplot2 [28].
Jof 10 00206 g001
Figure 2. Statistical parsimony network of phycobiont haplotypes calculated using the extended nrITS sequence sampling that included available GenBank data of Trebouxia from Ramalina farinacea; circle colors indicate the biogeographic subregions where specimens were collected (see legend below); the sizes of circles in each network are proportional to the numbers of individuals bearing the haplotype; circles may represent two or more haplotypes when these are separated only by indels; black-filled circles indicate missing haplotypes, and hatch marks indicate mutations. The geographic distributions of the two phycobionts of the R. farinacea group are represented over the map.
Figure 2. Statistical parsimony network of phycobiont haplotypes calculated using the extended nrITS sequence sampling that included available GenBank data of Trebouxia from Ramalina farinacea; circle colors indicate the biogeographic subregions where specimens were collected (see legend below); the sizes of circles in each network are proportional to the numbers of individuals bearing the haplotype; circles may represent two or more haplotypes when these are separated only by indels; black-filled circles indicate missing haplotypes, and hatch marks indicate mutations. The geographic distributions of the two phycobionts of the R. farinacea group are represented over the map.
Jof 10 00206 g002
Figure 3. Interaction networks among the phycobionts Trebouxia jamesii and T. lynnae and (A) the different mycobiont species in the “R. farinacea group”; (B) mycobiont genetic clusters delimited with fastbaps (I, IV, VI, VII and VIII); and (C) biogeographic subregions spanning from Macaronesian archipelagos to continental Europe. The width of the links is proportional to the number of specimens forming the association. The size of rectangles depicting each phycobiont haplotype is proportional to the number of obtained sequences; however, the sizes are not proportional among the three networks. Roman numerals represent different genetic clusters inferred with fastbaps. The biogeographic subregions [27] where specimens were collected are indicated as BE: Boreal Europe, CE: Central Europe, NCal: North California, Unkn: Unknown; AE: Atlantic Europe, Adr: Adriatic, Alp: Alpine, MWIP: Mediterranean West Iberian Peninsula, MCIP: Mediterranean Central Iberian Peninsula, IT: Italo-Thyrrenian, BCP: Balearic-Catalonian-Provençal, Mac: Macaronesian.
Figure 3. Interaction networks among the phycobionts Trebouxia jamesii and T. lynnae and (A) the different mycobiont species in the “R. farinacea group”; (B) mycobiont genetic clusters delimited with fastbaps (I, IV, VI, VII and VIII); and (C) biogeographic subregions spanning from Macaronesian archipelagos to continental Europe. The width of the links is proportional to the number of specimens forming the association. The size of rectangles depicting each phycobiont haplotype is proportional to the number of obtained sequences; however, the sizes are not proportional among the three networks. Roman numerals represent different genetic clusters inferred with fastbaps. The biogeographic subregions [27] where specimens were collected are indicated as BE: Boreal Europe, CE: Central Europe, NCal: North California, Unkn: Unknown; AE: Atlantic Europe, Adr: Adriatic, Alp: Alpine, MWIP: Mediterranean West Iberian Peninsula, MCIP: Mediterranean Central Iberian Peninsula, IT: Italo-Thyrrenian, BCP: Balearic-Catalonian-Provençal, Mac: Macaronesian.
Jof 10 00206 g003
Figure 4. (A) Principal coordinate analysis (PCA) of the 19 BioClim variables: 1 = annual mean temperature, 2 = mean diurnal range, 3 = isothermality, 4 = temperature seasonality, 5 = max temperature of warmest month, 6 = min temperature of coldest month, 7 = temperature annual range, 8 = mean temperature of wettest quarter, 9 = mean temperature of driest quarter, 10 = mean temperature of warmest quarter, 11 = mean temperature of coldest quarter, 12 = annual precipitation, 13 = precipitation of wettest month, 14 = precipitation of driest month, 15 = precipitation seasonality, 16 = precipitation of wettest quarter, 17 = precipitation of driest quarter, 18 = precipitation of warmest quarter, 19 = precipitation of coldest quarter. Climatic niche hypervolumes based on climatic PC1–PC2 axes (78% of variation explained) for (B) phycobiont species T. jamesii and T. lynnae and for (C) the four most abundant T. jamesii haplotypes.
Figure 4. (A) Principal coordinate analysis (PCA) of the 19 BioClim variables: 1 = annual mean temperature, 2 = mean diurnal range, 3 = isothermality, 4 = temperature seasonality, 5 = max temperature of warmest month, 6 = min temperature of coldest month, 7 = temperature annual range, 8 = mean temperature of wettest quarter, 9 = mean temperature of driest quarter, 10 = mean temperature of warmest quarter, 11 = mean temperature of coldest quarter, 12 = annual precipitation, 13 = precipitation of wettest month, 14 = precipitation of driest month, 15 = precipitation seasonality, 16 = precipitation of wettest quarter, 17 = precipitation of driest quarter, 18 = precipitation of warmest quarter, 19 = precipitation of coldest quarter. Climatic niche hypervolumes based on climatic PC1–PC2 axes (78% of variation explained) for (B) phycobiont species T. jamesii and T. lynnae and for (C) the four most abundant T. jamesii haplotypes.
Jof 10 00206 g004
Figure 5. Box-plot diagram representing differences in climate preferences in Trebouxia jamesii and T. lynnae. Climatic data were obtained from the Global Climate Data—WorldClim. (A) BIO1 = annual mean temperature (°C), (B) BIO6 = min temperature of coldest month (°C), (C) BIO16 = precipitation of wettest quarter and (D) BIO17 = precipitation of driest quarter.
Figure 5. Box-plot diagram representing differences in climate preferences in Trebouxia jamesii and T. lynnae. Climatic data were obtained from the Global Climate Data—WorldClim. (A) BIO1 = annual mean temperature (°C), (B) BIO6 = min temperature of coldest month (°C), (C) BIO16 = precipitation of wettest quarter and (D) BIO17 = precipitation of driest quarter.
Jof 10 00206 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Moya, P.; Chiva, S.; Pazos, T.; Barreno, E.; Carrasco, P.; Muggia, L.; Garrido-Benavent, I. Myco–Phycobiont Interactions within the “Ramalina farinacea Group”: A Geographical Survey over Europe and Macaronesia. J. Fungi 2024, 10, 206. https://doi.org/10.3390/jof10030206

AMA Style

Moya P, Chiva S, Pazos T, Barreno E, Carrasco P, Muggia L, Garrido-Benavent I. Myco–Phycobiont Interactions within the “Ramalina farinacea Group”: A Geographical Survey over Europe and Macaronesia. Journal of Fungi. 2024; 10(3):206. https://doi.org/10.3390/jof10030206

Chicago/Turabian Style

Moya, Patricia, Salvador Chiva, Tamara Pazos, Eva Barreno, Pedro Carrasco, Lucia Muggia, and Isaac Garrido-Benavent. 2024. "Myco–Phycobiont Interactions within the “Ramalina farinacea Group”: A Geographical Survey over Europe and Macaronesia" Journal of Fungi 10, no. 3: 206. https://doi.org/10.3390/jof10030206

APA Style

Moya, P., Chiva, S., Pazos, T., Barreno, E., Carrasco, P., Muggia, L., & Garrido-Benavent, I. (2024). Myco–Phycobiont Interactions within the “Ramalina farinacea Group”: A Geographical Survey over Europe and Macaronesia. Journal of Fungi, 10(3), 206. https://doi.org/10.3390/jof10030206

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop