Hidden Genetic Diversity in an Asexually Reproducing Lichen Forming Fungal Group

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.

Reference-Based RADseq Unravels the Evolutionary History of Polar Species in 'the Crux Lichenologorum' Genus Usnea (Parmeliaceae, Ascomycota)

Nearly 90% of fungal diversity, one of the most speciose branches in the tree of life, remains undescribed. Lichenized fungi as symbiotic associations are still a challenge for species delimitation, and current species diversity is vastly underestimated. The ongoing democratization of Next-Generation Sequencing is turning the tables. Particularly, reference-based RADseq allows for metagenomic filtering of the symbiont sequence and yields robust phylogenomic trees of closely related species. We implemented reference-based RADseq to disentangle the evolution of neuropogonoid lichens, which inhabit harsh environments and belong to Usnea (Parmeliaceae, Ascomycota), one of the most taxonomically intriguing genera within lichenized fungi. Full taxon coverage of neuropogonoid lichens was sampled for the first time, coupled with phenotype characterizations. More than 20,000 loci of 126 specimens were analyzed through concatenated and coalescent-based methods, including time calibrations. Our analysis addressed the major taxonomic discussions over recent decades. Subsequently, two species are newly described, namely U. aymondiana and U. fibriloides, and three species names are resurrected. The late Miocene and Pliocene-Pleistocene boundary is inferred as the timeframe for neuropogonoid lichen diversification. Ultimately, this study helped fill the gap of fungal diversity by setting a solid backbone phylogeny which raises new questions about which factors may trigger complex evolutionary scenarios.

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.

Population genomics of a reindeer lichen species from North American lichen woodlands

PREMISE

Lichens are one of the main structural components of plant communities in the North American boreal biome. They play a pivotal role in lichen woodlands, a large ecosystem situated north of the closed-crown forest zone, and south of the forest-tundra zone. In Eastern Canada (Quebec), there is a remnant LW found 500 km south of its usual distribution range, in the Parc National des Grands-Jardins, originated mainly because of wildfires. We inferred the origin of the lichen Cladonia stellaris from this LW and assessed its genetic diversity in a postfire succession.

METHODS

We genotyped 122 individuals collected across a latitudinal gradient in Quebec. Using the software Stacks, we compared four different approaches of locus selection and single-nucleotide polymorphism calling. We identified the best fitting approach to investigate population structure and estimate genetic diversity of C. stellaris.

RESULTS

Populations in southern Quebec are not genetically different from those of northern LWs. The species consists of at least four phylogenetic lineages with elevated levels of genetic diversity and low co-ancestry. In Parc National des Grands-Jardins, we reported high values of genetic diversity not related with time since fire disturbance and low genetic differentiation among populations with different fire histories.

CONCLUSIONS

This first population genomic study of C. stellaris is an important step forward to understand the origin and biogeographic patterns of lichen woodlands in North America. Our findings also contribute to the understanding of the effect of postfire succession on the genetic structure of the species.

Fungal species boundaries in the genomics era

Genomic data has opened new possibilities to understand how organisms change over time, and could enable the discovery of previously undescribed species. Although taxonomy used to be based on phenotypes, molecular data has frequently revealed that morphological traits are insufficient to describe biodiversity. Genomics holds the promise of revealing even more genetic discontinuities, but the parameters on how to describe species from genomic data remain unclear. Fungi have been a successful case in which the use of molecular markers has uncovered the existence of genetic boundaries where no crosses are possible. In this minireview, we highlight recent advances, propose a set of standards to use genomic sequences to uncover species boundaries, point out potential pitfalls, and present possible future research directions.

Evaluating methodologies for species delimitation: the mismatch between phenotypes and genotypes in lichenized fungi (Bryoria sect. Implexae, Parmeliaceae)

In many lichen-forming fungi, molecular phylogenetic analyses lead to the discovery of cryptic species within traditional morphospecies. However, in some cases, molecular sequence data also questions the separation of phenotypically characterised species. Here we apply an integrative taxonomy approach - including morphological, chemical, molecular, and distributional characters - to re-assess species boundaries in a traditionally speciose group of hair lichens, Bryoria sect. Implexae. We sampled multilocus sequence and microsatellite data from 142 specimens from a broad intercontinental distribution. Molecular data included DNA sequences of the standard fungal markers ITS, IGS, GAPDH, two newly tested loci (FRBi15 and FRBi16), and SSR frequencies from 18 microsatellite markers. Datasets were analysed with Bayesian and maximum likelihood phylogenetic reconstruction, phenogram reconstruction, STRUCTURE Bayesian clustering, principal coordinate analysis, haplotype network, and several different species delimitation analyses (ABGD, PTP, GMYC, and DISSECT). Additionally, past population demography and divergence times are estimated. The different approaches to species recognition do not support the monophyly of the 11 currently accepted morphospecies, and rather suggest the reduction of these to four phylogenetic species. Moreover, three of these are relatively recent in origin and cryptic, including phenotypically and chemically variable specimens. Issues regarding the integration of an evolutionary perspective into taxonomic conclusions in species complexes, which have undergone recent diversification, are discussed. The four accepted species, all epitypified by sequenced material, are Bryoria fuscescens, B. glabra, B. kockiana, and B. pseudofuscescens. Ten species rank names are reduced to synonymy. In the absence of molecular data, they can be recorded as the B. fuscescens complex. Intraspecific phenotype plasticity and factors affecting the speciation of different morphospecies in this group of Bryoria are outlined.

Discovering cryptic species in the Aspiciliella intermutans complex (Megasporaceae, Ascomycota) - First results using gene concatenation and coalescent-based species tree approaches

Taxonomic identifications in some groups of lichen-forming fungi have been challenge largely due to the scarcity of taxonomically relevant features and limitations of morphological and chemical characters traditionally used to distinguish closely related taxa. Delineating species boundaries in closely related species or species complexes often requires a range of multisource data sets and comprehensive analytical methods. Here we aim to examine species boundaries in a group of saxicolous lichen forming fungi, the Aspiciliella intermutans complex (Megasporaceae), widespread mainly in the Mediterranean. We gathered DNA sequences of the nuclear ribosomal internal transcribed spacer (nuITS), the nuclear large subunit (nuLSU), the mitochondrial small subunit (mtSSU) ribosomal DNA, and the DNA replication licensing factor MCM7 from 80 samples mostly from Iran, Caucasia, Greece and eastern Europe. We used a combination of phylogenetic strategies and a variety of empirical, sequence-based species delimitation approaches to infer species boundaries in this group. The latter included: the automatic barcode gap discovery (ABGD), the multispecies coalescent approach *BEAST and Bayesian Phylogenetics and Phylogeography (BPP) program. Different species delimitation scenarios were compared using Bayes factors species delimitation analysis. Furthermore, morphological, chemical, ecological and geographical features of the sampled specimens were examined. Our study uncovered cryptic species diversity in A. intermutans and showed that morphology-based taxonomy may be unreliable, underestimating species diversity in this group of lichens. We identified a total of six species-level lineages in the A. intermutans complex using inferences from multiple empirical operational criteria. We found little corroboration between morphological and ecological features with our proposed candidate species, while secondary metabolite data do not corroborate tree topology. The present study on the A. intermutans species-complex indicates that the genus Aspiciliella, as currently circumscribed, is more diverse in Eurasia than previously expected.

Formal description of sequence-based voucherless Fungi: promises and pitfalls, and how to resolve them

There is urgent need for a formal nomenclature of sequence-based, voucherless Fungi, given that environmental sequencing has accumulated more than one billion fungal ITS reads in the Sequence Read Archive, about 1,000 times as many as fungal ITS sequences in GenBank. These unnamed Fungi could help to bridge the gap between 115,000 to 140,000 currently accepted and 2.2 to 3.8 million predicted species, a gap that cannot realistically be filled using specimen or culture-based inventories. The Code never aimed at placing restrictions on the nature of characters chosen for taxonomy, and the requirement for physical types is now becoming a constraint on the advancement of science. We elaborate on the promises and pitfalls of sequence-based nomenclature and provide potential solutions to major concerns of the mycological community. Types of sequence-based taxa, which by default lack a physical specimen or culture, could be designated in four alternative ways: (1) the underlying sample ('bag' type), (2) the DNA extract, (3) fluorescent in situ hybridization (FISH), or (4) the type sequence itself. Only (4) would require changes to the Code and the latter would be the most straightforward approach, complying with three of the five principal functions of types better than physical specimens. A fifth way, representation of the sequence in an illustration, has been ruled as unacceptable in the Code. Potential flaws in sequence data are analogous to flaws in physical types, and artifacts are manageable if a stringent analytical approach is applied. Conceptual errors such as homoplasy, intragenomic variation, gene duplication, hybridization, and horizontal gene transfer, apply to all molecular approaches and cannot be used as a specific argument against sequence-based nomenclature. The potential impact of these phenomena is manageable, as phylogenetic species delimitation has worked satisfactorily in Fungi. The most serious shortcoming of sequence-based nomenclature is the likelihood of parallel classifications, either by describing taxa that already have names based on physical types, or by using different markers to delimit species within the same lineage. The probability of inadvertently establishing sequence-based species that have names available is between 20.4 % and 1.5 % depending on the number of globally predicted fungal species. This compares favourably to a historical error rate of about 30 % based on physical types, and this rate could be reduced to practically zero by adding specific provisions to this approach in the Code. To avoid parallel classifications based on different markers, sequence-based nomenclature should be limited to a single marker, preferably the fungal ITS barcoding marker; this is possible since sequence-based nomenclature does not aim at accurate species delimitation but at naming lineages to generate a reference database, independent of whether these lineages represent species, closely related species complexes, or infraspecies. We argue that clustering methods are inappropriate for sequence-based nomenclature; this approach must instead use phylogenetic methods based on multiple alignments, combined with quantitative species recognition methods. We outline strategies to obtain higher-level phylogenies for ITS-based, voucherless species, including phylogenetic binning, 'hijacking' species delimitation methods, and temporal banding. We conclude that voucherless, sequence-based nomenclature is not a threat to specimen and culture-based fungal taxonomy, but a complementary approach capable of substantially closing the gap between known and predicted fungal diversity, an approach that requires careful work and high skill levels.

Sequence-based nomenclature: a reply to Thines et al. and Zamora et al. and provisions for an amended proposal "from the floor" to allow DNA sequences as types of names

We reply to two recently published, multi-authored opinion papers by opponents of sequence-based nomenclature, namely Zamora et al. (IMA Fungus9: 167-175,2018) and Thines et al. (IMA Fungus9: 177-183, 2018). While we agree with some of the principal arguments brought forward by these authors, we address misconceptions and demonstrate that some of the presumed evidence presented in these papers has been wrongly interpreted. We disagree that allowing sequences as types would fundamentally alter the nature of types, since a similar nature of abstracted features as type is already allowed in the Code (Art. 40.5), namely an illustration. We also disagree that there is a high risk of introducing artifactual taxa, as this risk can be quantified at well below 5 %, considering the various types of high-throughput sequencing errors. Contrary to apparently widespread misconceptions, sequence-based nomenclature cannot be based on similarity-derived OTUs and their consensus sequences, but must be derived from rigorous, multiple alignment-based phylogenetic methods and quantitative, single-marker species recognition algorithms, using original sequence reads; it is therefore identical in its approach to single-marker studies based on physical types, an approach allowed by the Code. We recognize the limitations of the ITS as a single fungal barcoding marker, but point out that these result in a conservative approach, with "false negatives" surpassing "false positives"; a desirable feature of sequence-based nomenclature. Sequence-based nomenclature does not aim at accurately resolving species, but at naming sequences that represent unknown fungal lineages so that these can serve as a means of communication, so ending the untenable situation of an exponentially growing number of unlabeled fungal sequences that fill online repositories. The risks are outweighed by the gains obtained by a reference library of named sequences spanning the full array of fungal diversity. Finally, we elaborate provisions in addition to our original proposal to amend the Code that would take care of the issues brought forward by opponents to this approach. In particular, taking up the idea of the Candidatus status of invalid, provisional names in prokaryote nomenclature, we propose a compromise that would allow valid publication of voucherless, sequence-based names in a consistent manner, but with the obligate designation as "nom. seq." (nomen sequentiae). Such names would not have priority over specimen- or culture-based names unless either epitypified with a physical type or adopted for protection on the recommendation of a committee of the International Commission on the Taxonomy of Fungi following evaluation based on strict quality control of the underlying studies based on established rules or recommendations.

Oligocene origin and drivers of diversification in the genus Sticta (Lobariaceae, Ascomycota)

A major challenge to evolutionary biologists is to understand how biodiversity is distributed through space and time and across the tree of life. Diversification of organisms is influenced by many factors that act at different times and geographic locations but it is still not clear which have a significant impact and how drivers interact. To study diversification, we chose the lichen genus Sticta, by sampling through most of the global range and producing a time tree. We estimate that Sticta originated about 30 million years ago, but biogoegraphic analysis was unclear in estimating the origin of the genus. Furthermore, we investigated the effect of dispersal ability finding that Sticta has a high dispersal rate, as collections from Hawaii showed that divergent lineages colonized the islands at least four times. Symbiont interactions were investigated using BiSSE to understand if green-algal or cyanobacterial symbiont interactions influenced diversification, only to find that the positive results were driven almost completely by Type I error. On the other hand, another BiSSE analysis found that an association with Andean tectonic activity increases the speciation rate of species.

Fungal Diversity Revisited: 2.2 to 3.8 Million Species

The question of how many species of Fungi there are has occasioned much speculation, with figures mostly posited from around half a million to 10 million, and in one extreme case even a sizable portion of the spectacular number of 1 trillion. Here we examine new evidence from various sources to derive an updated estimate of global fungal diversity. The rates and patterns in the description of new species from the 1750s show no sign of approaching an asymptote and even accelerated in the 2010s after the advent of molecular approaches to species delimitation. Species recognition studies of (semi-)cryptic species hidden in morpho-species complexes suggest a weighted average ratio of about an order of magnitude for the number of species recognized after and before such studies. New evidence also comes from extrapolations of plant:fungus ratios, with information now being generated from environmental sequence studies, including comparisons of molecular and fieldwork data from the same sites. We further draw attention to undescribed species awaiting discovery in biodiversity hot spots in the tropics, little-explored habitats (such as lichen-inhabiting fungi), and material in collections awaiting study. We conclude that the commonly cited estimate of 1.5 million species is conservative and that the actual range is properly estimated at 2.2 to 3.8 million. With 120,000 currently accepted species, it appears that at best just 8%, and in the worst case scenario just 3%, are named so far. Improved estimates hinge particularly on reliable statistical and phylogenetic approaches to analyze the rapidly increasing amount of environmental sequence data.

Importance of Resolving Fungal Nomenclature: the Case of Multiple Pathogenic Species in the Cryptococcus Genus

Cryptococcosis is a major fungal disease caused by members of the Cryptococcus gattii and Cryptococcus neoformans species complexes. After more than 15 years of molecular genetic and phenotypic studies and much debate, a proposal for a taxonomic revision was made. The two varieties within C. neoformans were raised to species level, and the same was done for five genotypes within C. gattii. In a recent perspective (K. J. Kwon-Chung et al., mSphere 2:e00357-16, 2017, https://doi.org/10.1128/mSphere.00357-16), it was argued that this taxonomic proposal was premature and without consensus in the community. Although the authors of the perspective recognized the existence of genetic diversity, they preferred the use of the informal nomenclature "C. neoformans species complex" and "C. gattii species complex." Here we highlight the advantage of recognizing these seven species, as ignoring these species will impede deciphering further biologically and clinically relevant differences between them, which may in turn delay future clinical advances.