Conserved genomic collinearity as a source of broadly applicable, fast evolving, markers to resolve species complexes: A case study using the lichen-forming genus Peltigera section Polydactylon

The mitogenomes characterization of two Peltigera species (Peltigera elisabethae and Peltigera polydactylon) and comparative mitogenomic analyses of six Peltigera

In the present study, the mitochondrial genomes of Peltigera elisabethae and P. polydactylon were sequenced and assembled. The two mitogenomes were composed of circular DNA molecules, with sizes of 64,034 bp and 59,208 bp, respectively. Comparative analysis showed that the genome size, GC content, GC skew, and AT skew varied between the two mitochondrial genomes. In codon analysis, phenylalanine (Phe), isoleucine (Ile), and leucine (Leu) were most frequently used in six Peltigera genomes. Evolutionary analysis showed that all 14 protein-coding genes (PCGs) were subject to purifying selection in the six Peltigera species. Regarding gene rearrangement, the PCGs of Peltigera had the same gene sequence and gene content, and a few intron sequences and spacer sequences were rearranged in Peltigera. In the phylogenetic analysis, we used Bayesian and ML methods to construct a phylogenetic tree. Two phylogenetic trees with consistent topology with high support indicate that mitochondrial genes were reliable molecular markers for analyzing the phylogenetic relationships. The present study enriches the mitochondrial genome data of Peltigera and promotes further understanding of the genetics and evolution of the Peltigera genus.

Isolation of starch and protein degrading strain Bacillus subtilis FYZ1-3 from tobacco waste and genomic analysis of its tolerance to nicotine and inhibition of fungal growth

Aerobic fermentation is an effective technique for the large-scale processing of tobacco waste. However, the specificity of the structure and composition of tobacco-derived organic matter and the toxic alkaloids in the material make it currently difficult to directly use microbial agents. In this study, a functional strain FYZ1-3 was isolated and screened from thermophilic phase samples of tobacco waste composting. This strain could withstand temperatures as high as 80°C and grow normally at 0.6% nicotine content. Furthermore, it had a strong decomposition capacity of tobacco-derived starch and protein, with amylase activity of 122.3  U/mL and protease activity and 52.3  U/mL, respectively. To further understand the mechanism of the metabolic transformation of the target, whole genome sequencing was used and the secondary metabolite gene cluster was predicted. The inhibitory effect of the strain on common tobacco fungi was verified using the plate confrontation and agar column methods. The results showed that the strain FYZ1-3 was Bacillus subtilis, with a genome size of 4.17  Mb and GC content of 43.68%; 4,338 coding genes were predicted. The genome was annotated and analyzed using multiple databases to determine its ability to efficiently degrade starch proteins at the molecular level. Moreover, 14 functional genes related to nicotine metabolism were identified, primarily located on the distinct genomic island of FYZ1-3, giving a speculation for its nicotine tolerance capability on the molecular mechanism. By mining the secondary metabolite gene cluster prediction, we found potential synthetic bacteriocin, antimicrobial peptide, and other gene clusters on its chromosome, which may have certain antibacterial properties. Further experiments confirmed that the FYZ1-3 strain was a potent growth inhibitor of Penicillium chrysogenum, Aspergillus sydowii, A. fumigatus, and Talaromyces funiculosus. The creation and industrial use of the functional strains obtained in this study provide a theoretical basis for its industrial use, where it would be of great significance to improve the utilization rate of tobacco waste.

High species richness in the lichen genus Peltigera (Ascomycota, Lecanoromycetes): 34 species in the dolichorhizoid and scabrosoid clades of section Polydactylon, including 24 new to science

Applying molecular methods to fungi establishing lichenized associations with green algae or cyanobacteria has repeatedly revealed the existence of numerous phylogenetic taxa overlooked by classical taxonomic approaches. Here, we report taxonomical conclusions based on multiple species delimitation and validation analyses performed on an eight-locus dataset that includes world-wide representatives of the dolichorhizoid and scabrosoid clades in section Polydactylon of the genus Peltigera. Following the recommendations resulting from a consensus species delimitation approach and additional species validation analysis (BPP) performed in this study, we present a total of 25 species in the dolichorhizoid clade and nine in the scabrosoid clade, including respectively 18 and six species that are new to science and formally described. Additionally, one combination and three varieties (including two new to science) are proposed in the dolichorhizoid clade. The following 24 new species are described: P. appalachiensis, P. asiatica, P. borealis, P. borinquensis, P. chabanenkoae, P. clathrata, P. elixii, P. esslingeri, P. flabellae, P. gallowayi, P. hawaiiensis, P. holtanhartwigii, P. itatiaiae, P. hokkaidoensis, P. kukwae, P. massonii, P. mikado, P. nigriventris, P. orientalis, P. rangiferina, P. sipmanii, P. stanleyensis, P. vitikainenii and P. willdenowii; the following new varieties are introduced: P. kukwae var. phyllidiata and P. truculenta var. austroscabrosa; and the following new combination is introduced: P. hymenina var. dissecta. Each species from the dolichorhizoid and scabrosoid clades is morphologically and chemically described, illustrated, and characterised with ITS sequences. Identification keys are provided for the main biogeographic regions where species from the two clades occur. Morphological and chemical characters that are commonly used for species identification in the genus Peltigera cannot be applied to unambiguously recognise most molecularly circumscribed species, due to high variation of thalli formed by individuals within a fungal species, including the presence of distinct morphs in some cases, or low interspecific variation in others. The four commonly recognised morphospecies: P. dolichorhiza, P. neopolydactyla, P. pulverulenta and P. scabrosa in the dolichorhizoid and scabrosoid clades represent species complexes spread across multiple and often phylogenetically distantly related lineages. Geographic origin of specimens is often helpful for species recognition; however, ITS sequences are frequently required for a reliable identification. Citation: Magain N, Miadlikowska J, Goffinet B, et al. 2023. High species richness in the lichen genus Peltigera (Ascomycota, Lecanoromycetes): 34 species in the dolichorhizoid and scabrosoid clades of section Polydactylon, including 24 new to science. Persoonia 51: 1-88. doi: 10.3767/persoonia.2023.51.01.

Phylogenetic structure of specialization: A new approach that integrates partner availability and phylogenetic diversity to quantify biotic specialization in ecological networks

Biotic specialization holds information about the assembly, evolution, and stability of biological communities. Partner availabilities can play an important role in enabling species interactions, where uneven partner availabilities can bias estimates of biotic specialization when using phylogenetic diversity indices. It is therefore important to account for partner availability when characterizing biotic specialization using phylogenies. We developed an index, phylogenetic structure of specialization (PSS), that avoids bias from uneven partner availabilities by uncoupling the null models for interaction frequency and phylogenetic distance. We incorporate the deviation between observed and random interaction frequencies as weights into the calculation of partner phylogenetic α‐diversity. To calculate the PSS index, we then compare observed partner phylogenetic α‐diversity to a null distribution generated by randomizing phylogenetic distances among the same number of partners. PSS quantifies the phylogenetic structure (i.e., clustered, overdispersed, or random) of the partners of a focal species. We show with simulations that the PSS index is not correlated with network properties, which allows comparisons across multiple systems. We also implemented PSS on empirical networks of host–parasite, avian seed‐dispersal, lichenized fungi–cyanobacteria, and hummingbird pollination interactions. Across these systems, a large proportion of taxa interact with phylogenetically random partners according to PSS, sometimes to a larger extent than detected with an existing method that does not account for partner availability. We also found that many taxa interact with phylogenetically clustered partners, while taxa with overdispersed partners were rare. We argue that species with phylogenetically overdispersed partners have often been misinterpreted as generalists when they should be considered specialists. Our results highlight the important role of randomness in shaping interaction networks, even in highly intimate symbioses, and provide a much‐needed quantitative framework to assess the role that evolutionary history and symbiotic specialization play in shaping patterns of biodiversity. PSS is available as an R package at https://github.com/cjpardodelahoz/pss.

Species in lichen-forming fungi: balancing between conceptual and practical considerations, and between phenotype and phylogenomics

Lichens are symbiotic associations resulting from interactions among fungi (primary and secondary mycobionts), algae and/or cyanobacteria (primary and secondary photobionts), and specific elements of the bacterial microbiome associated with the lichen thallus. The question of what is a species, both concerning the lichen as a whole and its main fungal component, the primary mycobiont, has faced many challenges throughout history and has reached new dimensions with the advent of molecular phylogenetics and phylogenomics. In this paper, we briefly revise the definition of lichens and the scientific and vernacular naming conventions, concluding that the scientific, Latinized name usually associated with lichens invariably refers to the primary mycobiont, whereas the vernacular name encompasses the entire lichen. Although the same lichen mycobiont may produce different phenotypes when associating with different photobionts or growing in axenic culture, this discrete variation does not warrant the application of different scientific names, but must follow the principle "one fungus = one name". Instead, broadly agreed informal designations should be used for such discrete morphologies, such as chloromorph and cyanomorph for lichens formed by the same mycobiont but with either green algae or cyanobacteria. The taxonomic recognition of species in lichen-forming fungi is not different from other fungi and conceptual and nomenclatural approaches follow the same principles. We identify a number of current challenges and provide recommendations to address these. Species delimitation in lichen-forming fungi should not be tailored to particular species concepts but instead be derived from empirical evidence, applying one or several of the following principles in what we call the LPR approach: lineage (L) coherence vs. divergence (phylogenetic component), phenotype (P) coherence vs. divergence (morphological component), and/or reproductive (R) compatibility vs. isolation (biological component). Species hypotheses can be established based on either L or P, then using either P or L (plus R) to corroborate them. The reliability of species hypotheses depends not only on the nature and number of characters but also on the context: the closer the relationship and/or similarity between species, the higher the number of characters and/or specimens that should be analyzed to provide reliable delimitations. Alpha taxonomy should follow scientific evidence and an evolutionary framework but should also offer alternative practical solutions, as long as these are scientifically defendable. Taxa that are delimited phylogenetically but not readily identifiable in the field, or are genuinely cryptic, should not be rejected due to the inaccessibility of proper tools. Instead, they can be provisionally treated as undifferentiated complexes for purposes that do not require precise determinations. The application of infraspecific (gamma) taxonomy should be restricted to cases where there is a biological rationale, i.e., lineages of a species complex that show limited phylogenetic divergence but no evidence of reproductive isolation. Gamma taxonomy should not be used to denote discrete phenotypical variation or ecotypes not warranting the distinction at species level. We revise the species pair concept in lichen-forming fungi, which recognizes sexually and asexually reproducing morphs with the same underlying phenotype as different species. We conclude that in most cases this concept does not hold, but the actual situation is complex and not necessarily correlated with reproductive strategy. In cases where no molecular data are available or where single or multi-marker approaches do not provide resolution, we recommend maintaining species pairs until molecular or phylogenomic data are available. This recommendation is based on the example of the species pair Usnea aurantiacoatra vs. U. antarctica, which can only be resolved with phylogenomic approaches, such as microsatellites or RADseq. Overall, we consider that species delimitation in lichen-forming fungi has advanced dramatically over the past three decades, resulting in a solid framework, but that empirical evidence is still missing for many taxa. Therefore, while phylogenomic approaches focusing on particular examples will be increasingly employed to resolve difficult species complexes, broad screening using single barcoding markers will aid in placing as many taxa as possible into a molecular matrix. We provide a practical protocol how to assess and formally treat taxonomic novelties. While this paper focuses on lichen fungi, many of the aspects discussed herein apply generally to fungal taxonomy. The new combination Arthonia minor (Lücking) Lücking comb. et stat. nov. (Bas.: Arthonia cyanea f. minor Lücking) is proposed.

Exploring syntenic conservation across genomes for phylogenetic studies of organisms subjected to horizontal gene transfers: A case study with Cyanobacteria and cyanolichens

Understanding the evolutionary history of symbiotic Cyanobacteria at a fine scale is essential to unveil patterns of associations with their hosts and factors driving their spatiotemporal interactions. As for bacteria in general, Horizontal Gene Transfers (HGT) are expected to be rampant throughout their evolution, which justified the use of single-locus phylogenies in macroevolutionary studies of these photoautotrophic bacteria. Genomic approaches have greatly increased the amount of molecular data available, but the selection of orthologous, congruent genes that are more likely to reflect bacterial macroevolutionary histories remains problematic. In this study, we developed a synteny-based approach and searched for Collinear Orthologous Regions (COR), under the assumption that genes that are present in the same order and orientation across a wide monophyletic clade are less likely to have undergone HGT. We searched sixteen reference Nostocales genomes and identified 99 genes, part of 28 COR comprising three to eight genes each. We then developed a bioinformatic pipeline, designed to minimize inter-genome contamination and processed twelve Nostoc-associated lichen metagenomes. This reduced our original dataset to 90 genes representing 25 COR, which were used to infer phylogenetic relationships within Nostocales and among lichenized Cyanobacteria. This dataset was narrowed down further to 71 genes representing 22 COR by selecting only genes part of one (largest) operon per COR. We found a relatively high level of congruence among trees derived from the 90-gene dataset, but congruence was only slightly higher among genes within a COR compared to genes across COR. However, topological congruence was significantly higher among the 71 genes part of one operon per COR. Nostocales phylogenies resulting from concatenation and species tree approaches based on the 90- and 71-gene datasets were highly congruent, but the most highly supported result was obtained when using synteny, collinearity, and operon information (i.e., 71-gene dataset) as gene selection criteria, which outperformed larger datasets with more genes.

Towards Genomic Criteria for Delineating Fungal Species

The discussion of fungal species delineation has yet to reach a consensus, despite the advancements in technology, which helped modernise traditional approaches. In particular, the phylogenetic species concept was one of the tools that has been used with considerable success across the fungal kingdom. The fast rise of fungal genomics provides an unprecedented opportunity to expand measuring the relatedness of fungal strains to the level of whole genomes. However, the use of genomic information in taxonomy has only just begun, and few methodological guidelines have been suggested so far. Here, a simple approach of computationally measuring genomic distances and their use as a standard for species delineation is investigated. A fixed threshold genomic distance calculated by the quick and easy-to-use tools Mash and Dashing proved to be an unexpectedly widely applicable and robust criterion for determining whether two genomes belong to the same or to different species. The accuracy of species delineation in an uncurated dataset of GenBank fungal genomes was close to 90%—and exceeded 90% with minimal curation. As expected, the discriminative power of this approach was lower at higher taxonomic ranks, but still significantly larger than zero. Simple instructions for calculation of a genomic distance between two genomes and species similarity thresholds at different k-mer sizes are suggested. The calculation of genomic distance is identified as a powerful approach for delineating fungal species and is proposed—not as the only criterion—but as an additional tool in the versatile toolbox of fungal taxonomy.

The fruticose genera in the Ramalinaceae (Ascomycota, Lecanoromycetes): their diversity and evolutionary history

We present phylogenetic analyses of the fruticose Ramalinaceae based on extensive collections from many parts of the world, with a special focus on the Vizcaíno deserts in north-western Mexico and the coastal desert in Namibia. We generate a four-locus DNA sequence dataset for accessions of Ramalina and two additional loci for Niebla and Vermilacinia. Four genera are strongly supported: the subcosmopolitan Ramalina, the new genus Namibialina endemic to SW Africa, and a duo formed by Niebla and Vermilacinia, endemic to the New World except the sorediate V. zebrina that disjunctly occurs in Namibia. The latter three genera are restricted to coastal desert and chaparral where vegetation depends on moisture from ocean fog. Ramalina is subcosmopolitan and much more diverse in its ecology. We show that Ramalina and its sister genus Namibialina diverged from each other at c. 48 Myrs, whereas Vermilacinia and Niebla split at c. 30 Myrs. The phylogeny of the fruticose genera remains unresolved to their ancestral crustose genera. Species delimitation within Namibialina and Ramalina is rather straightforward. The phylogeny and taxonomy of Vermilacinia are fully resolved, except for the two youngest clades of corticolous taxa, and support current taxonomy, including four new taxa described here. Secondary metabolite variation in Niebla generally coincides with major clades which are comprised of species complexes with still unresolved phylogenetic relationships. A micro-endemism pattern of allopatric species is strongly suspected for both genera, except for the corticolous taxa within Vermilacinia. Both Niebla and saxicolous Vermilacinia have chemotypes unique to species clades that are largely endemic to the Vizcaíno deserts. The following new taxa are described: Namibialina gen. nov. with N. melanothrix (comb. nov.) as type species, a single new species of Ramalina (R. krogiae) and four new species of Vermilacinia (V. breviloba, V. lacunosa, V. pustulata and V. reticulata). The new combination V. granulans is introduced. Two epithets are re-introduced for European Ramalina species: R. crispans (= R. peruviana auct. eur.) and R. rosacea (= R. bourgeana auct. p.p). A lectotype is designated for Vermilacinia procera. A key to saxicolous species of Vermilacinia is presented.

Unambiguous identification of fungi: where do we stand and how accurate and precise is fungal DNA barcoding?

True fungi (Fungi) and fungus-like organisms (e.g. Mycetozoa, Oomycota) constitute the second largest group of organisms based on global richness estimates, with around 3 million predicted species. Compared to plants and animals, fungi have simple body plans with often morphologically and ecologically obscure structures. This poses challenges for accurate and precise identifications. Here we provide a conceptual framework for the identification of fungi, encouraging the approach of integrative (polyphasic) taxonomy for species delimitation, i.e. the combination of genealogy (phylogeny), phenotype (including autecology), and reproductive biology (when feasible). This allows objective evaluation of diagnostic characters, either phenotypic or molecular or both. Verification of identifications is crucial but often neglected. Because of clade-specific evolutionary histories, there is currently no single tool for the identification of fungi, although DNA barcoding using the internal transcribed spacer (ITS) remains a first diagnosis, particularly in metabarcoding studies. Secondary DNA barcodes are increasingly implemented for groups where ITS does not provide sufficient precision. Issues of pairwise sequence similarity-based identifications and OTU clustering are discussed, and multiple sequence alignment-based phylogenetic approaches with subsequent verification are recommended as more accurate alternatives. In metabarcoding approaches, the trade-off between speed and accuracy and precision of molecular identifications must be carefully considered. Intragenomic variation of the ITS and other barcoding markers should be properly documented, as phylotype diversity is not necessarily a proxy of species richness. Important strategies to improve molecular identification of fungi are: (1) broadly document intraspecific and intragenomic variation of barcoding markers; (2) substantially expand sequence repositories, focusing on undersampled clades and missing taxa; (3) improve curation of sequence labels in primary repositories and substantially increase the number of sequences based on verified material; (4) link sequence data to digital information of voucher specimens including imagery. In parallel, technological improvements to genome sequencing offer promising alternatives to DNA barcoding in the future. Despite the prevalence of DNA-based fungal taxonomy, phenotype-based approaches remain an important strategy to catalog the global diversity of fungi and establish initial species hypotheses.

Optimal markers for the identification of Colletotrichum species

Colletotrichum is among the most important genera of fungal plant pathogens. Molecular phylogenetic studies over the last decade have resulted in a much better understanding of the evolutionary relationships and species boundaries within the genus. There are now approximately 200 species accepted, most of which are distributed among 13 species complexes. Given their prominence on agricultural crops around the world, rapid identification of a large collection of Colletotrichum isolates is routinely needed by plant pathologists, regulatory officials, and fungal biologists. However, there is no agreement on the best molecular markers to discriminate species in each species complex. Here we calculate the barcode gap distance and intra/inter-specific distance overlap to evaluate each of the most commonly applied molecular markers for their utility as a barcode for species identification. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), histone-3 (HIS3), DNA lyase (APN2), intergenic spacer between DNA lyase and the mating-type locus MAT1-2-1 (APN2/MAT-IGS), and intergenic spacer between GAPDH and a hypothetical protein (GAP2-IGS) have the properties of good barcodes, whereas sequences of actin (ACT), chitin synthase (CHS-1) and nuclear rDNA internal transcribed spacers (nrITS) are not able to distinguish most species. Finally, we assessed the utility of these markers for phylogenetic studies using phylogenetic informativeness profiling, the genealogical sorting index (GSI), and Bayesian concordance analyses (BCA). Although GAPDH, HIS3 and β-tubulin (TUB2) were frequently among the best markers, there was not a single set of markers that were best for all species complexes. Eliminating markers with low phylogenetic signal tends to decrease uncertainty in the topology, regardless of species complex, and leads to a larger proportion of markers that support each lineage in the Bayesian concordance analyses. Finally, we reconstruct the phylogeny of each species complex using a minimal set of phylogenetic markers with the strongest phylogenetic signal and find the majority of species are strongly supported as monophyletic.

Species in section Peltidea (aphthosa group) of the genus Peltigera remain cryptic after molecular phylogenetic revision

Closely related lichen-forming fungal species circumscribed using phenotypic traits (morphospecies) do not always align well with phylogenetic inferences based on molecular data. Using multilocus data obtained from a worldwide sampling, we inferred phylogenetic relationships among five currently accepted morphospecies of Peltigera section Peltidea (P. aphthosa group). Monophyletic circumscription of all currently recognized morphospecies (P. britannica, P. chionophila, P. frippii and P. malacea) except P. aphthosa, which contained P. britannica, was confirmed with high bootstrap support. Following their re-delimitation using bGMYC and Structurama, BPP validated 14 putative species including nine previously unrecognized potential species (five within P. malacea, five within P. aphthosa, and two within P. britannica). Because none of the undescribed potential species are corroborated morphologically, chemically, geographically or ecologically, we concluded that these monophyletic entities represent intraspecific phylogenetic structure, and, therefore, should not be recognized as new species. Cyanobionts associated with Peltidea mycobionts (51 individuals) represented 22 unique rbcLX haplotypes from five phylogroups in Clade II subclades 2 and 3. With rare exceptions, Nostoc taxa involved in trimembered and bimembered associations are phylogenetically closely related (subclade 2) or identical, suggesting a mostly shared cyanobiont pool with infrequent switches. Based on a broad geographical sampling, we confirm a high specificity of Nostoc subclade 2 with their mycobionts, including a mutualistically exclusive association between phylogroup III and specific lineages of P. malacea.

Contrasting Symbiotic Patterns in Two Closely Related Lineages of Trimembered Lichens of the Genus Peltigera

Species circumscription is key to the characterization of patterns of specificity in symbiotic systems at a macroevolutionary scale. Here, a worldwide phylogenetic framework was used to assess the biodiversity and symbiotic patterns of association among partners in trimembered lichens from the genus Peltigera, section Chloropeltigera. We sequenced six loci of the main fungal partner and performed species discovery and validation analyses to establish putative species boundaries. Single locus phylogenies were used to establish the identity of both photobionts, Nostoc (cyanobacterium) and Coccomyxa (green alga). Distribution and specificity patterns were compared to the closely related clade, section Peltidea, which includes mainly Peltigera species with trimembered thalli. For section Chloropeltigera, eight fungal species (including five newly delimited putative species) were found in association with nine Nostoc phylogroups and two Coccomyxa species. In contrast, eight fungal species (including three newly delimited putative species) in section Peltidea were found in association with only four Nostoc phylogroups and the same two Coccomyxa species as for section Chloropeltigera. This difference in cyanobiont biodiversity between these two sections can potentially be explained by a significantly higher frequency of sexual reproductive structures in species from section Chloropeltigera compared to section Peltidea. Therefore, horizontal transmission of the cyanobiont might be more prevalent in Chloropeltigera species, while vertical transmission might be more common in Peltidea species. All Peltigera species in section Chloropeltigera are generalists in their association with Nostoc compared to more specialized Peltigera species in section Peltidea. Constrained distributions of Peltigera species that associate strictly with one species of green algae (Coccomyxa subellipsoidea) indicate that the availability of the green alga and the specificity of the interaction might be important factors limiting geographic ranges of trimembered Peltigera, in addition to constraints imposed by their interaction with Nostoc partners and by climatic factors.

Bioclimatic factors at an intrabiome scale are more limiting than cyanobiont availability for the lichen-forming genus Peltigera

PREMISE OF THE STUDY

Factors shaping spatiotemporal patterns of associations in mutualistic systems are poorly understood. We used the lichen-forming fungi Peltigera and their cyanobacterial partners Nostoc to investigate the spatial structure of this symbiosis at an intrabiome scale and to identify potential factors shaping these associations.

METHODS

Ninety-three thalli were sampled in Québec, Canada, along a south-north and an east-west transect of ~1300 km each. We identified the two main partners (Peltigera species and Nostoc phylogroups) using molecular markers and modeled the effects of environmental variables and partner occurrence on Peltigera-Nostoc distributions.

KEY RESULTS

Peltigera species showed a high degree of specialization toward cyanobionts, whereas two Nostoc phylogroups dominated both transects by associating with several Peltigera species. Peltigera species had narrower ranges than these two main cyanobionts. Distributions of three Peltigera species were highly associated with precipitation and temperature variables, which was not detected for Nostoc phylogroups at this spatial scale.

CONCLUSIONS

For these cyanolichens, factors driving patterns of symbiotic associations are scale dependent. Contrary to global-scale findings, generalist Peltigera species were not more widespread within the boreal biome than specialists. Nostoc availability was not the only driver of Peltigera species' geographic ranges; environmental factors also contributed to their intrabiome distributions. Climatic conditions (especially precipitation) limited the range of some Peltigera species more than the range of their cyanobacterial partners at an intrabiome (boreal) scale.

Strong specificity and network modularity at a very fine phylogenetic scale in the lichen genus Peltigera

Identifying the drivers and evolutionary consequences of species interactions is a major goal of community ecology. Network-based analyses can provide mathematical tools to detect non-random patterns of interactions, and potentially help predicting the consequences of such patterns on evolutionary dynamics of symbiotic systems. Here, we characterize the structure of a lichen network at a very fine phylogenetic scale, by identifying the photosynthetic partners (i.e., cyanobacteria of the genus Nostoc) of lichenized fungi belonging to a monophyletic section of a single genus (i.e., section Polydactylon of the genus Peltigera), worldwide. Even at such a fine phylogenetic scale, we found that interactions were highly modular and anti-nested, indicating strong preferences in interactions. When considering local Peltigera communities, i.e., datasets at small spatial scales with only a slightly broader phylogenetic range, interactions remained modular but were asymmetric, with generalist Nostoc partners interacting with specialized Peltigera species. This asymmetry was not detected with our global spatial scale dataset. We discuss these results in the light of lichen community assembly, and explore how such interaction patterns may influence coevolution in lichens and the evolutionary stability of the mutualism in general.

High diversity, high insular endemism and recent origin in the lichen genus Sticta (lichenized Ascomycota, Peltigerales) in Madagascar and the Mascarenes

Lichen biodiversity and its generative evolutionary processes are practically unknown in the MIOI (Madagascar and Indian Ocean Islands) biodiversity hotspot. We sought to test the hypothesis that lichenized fungi in this region have undergone a rapid radiation, following a single colonization event, giving rise to narrow endemics, as is characteristic of other lineages of plants. We extensively sampled specimens of the lichen genus Sticta in the Mascarene archipelago (mainly Réunion) and in Madagascar, mainly in the northern range (Amber Mt and Marojejy Mt) and produced the fungal ITS barcode sequence for 148 thalli. We further produced a four-loci data matrix for 68 of them, representing the diversity and geographical distribution of ITS haplotypes. We reconstructed the phylogenetic relationships within this group, established species boundaries with morphological context, and estimated the date of the most recent common ancestor. Our inferences resolve a robust clade comprising 31 endemic species of Sticta that arose from the diversification following a single recent (c. 11 Mya) colonization event. All but three species have a very restricted range, endemic to either the Mascarene archipelago or a single massif in Madagascar. The first genus of lichens to be studied with molecular data in this region underwent a recent radiation, exhibits micro-endemism, and thus exemplifies the biodiversity characteristics found in other taxa in Madagascar and the Mascarenes.