One hundred new species of lichenized fungi: a signature of undiscovered global diversity

High‐throughput metabarcoding studies on fungi and other eukaryotic microorganisms are rapidly becoming more frequent and more complex, requiring researchers to handle ever increasing amounts of raw sequence data. Here, we provide a flexible pipeline for pruning and analyzing fungal barcode (ITS rDNA) data generated as paired‐end reads on Illumina MiSeq sequencers. The pipeline presented includes specific steps fine‐tuned for ITS, that are mostly missing from pipelines developed for prokaryotes. It (1) employs state of the art programs and follows best practices in fungal high‐throughput metabarcoding; (2) consists of modules and scripts easily modifiable by the user to ensure maximum flexibility with regard to specific needs of a project or future methodological developments; and (3) is straightforward to use, also in classroom settings. We provide detailed descriptions and revision techniques for each step, thus giving the user maximum control over data treatment and avoiding a black‐box approach. Employing this pipeline will improve and speed up the tedious and error‐prone process of cleaning fungal Illumina metabarcoding data.

Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution. A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting. The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.