Next Article in Journal
Diversity, Abundance, and Ecological Roles of Planktonic Fungi in Marine Environments
Next Article in Special Issue
Antibiotic-Induced Treatments Reveal Stress-Responsive Gene Expression in the Endangered Lichen Lobaria pulmonaria
Previous Article in Journal
Phylogeny of Crepidotus applanatus Look-Alikes Reveals a Convergent Morphology Evolution and a New Species C. pini
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JoF
Volume 8
Issue 5
10.3390/jof8050490
jof-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Providing Scale to a Known Taxonomic Unknown—At Least a 70-Fold Increase in Species Diversity in a Cosmopolitan Nominal Taxon of Lichen-Forming Fungi

by
Yanyun Zhang
1,2,†,
Jeffrey Clancy
3,†,
Jacob Jensen
3,
Richard Troy McMullin
4,
Lisong Wang
1,* and
Steven D. Leavitt
5,*
1
Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Heilongtan, Kunming 650201, China
2
College of Life Science, Anhui Normal University, Wuhu 241000, China
3
Department of Biology, Brigham Young University, 4102 Life Science Building, Provo, UT 84602, USA
4
Canadian Museum of Nature, Research and Collections, Ottawa, ON K1P 6P4, Canada
5
Department of Biology, M. L. Bean Life Science Museum, Brigham Young University, 4102 Life Science Building, Provo, UT 84602, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
J. Fungi 2022, 8(5), 490; https://doi.org/10.3390/jof8050490
Submission received: 3 March 2022 / Revised: 3 May 2022 / Accepted: 4 May 2022 / Published: 8 May 2022
(This article belongs to the Special Issue Species Diversity, Ecology and Evolutionary History of Lichenized Fungi)
Browse Figures
Review Reports Versions Notes

Abstract

:
Robust species delimitations provide a foundation for investigating speciation, phylogeography, and conservation. Here we attempted to elucidate species boundaries in the cosmopolitan lichen-forming fungal taxon Lecanora polytropa. This nominal taxon is morphologically variable, with distinct populations occurring on all seven continents. To delimit candidate species, we compiled ITS sequence data from populations worldwide. For a subset of the samples, we also generated alignments for 1209 single-copy nuclear genes and an alignment spanning most of the mitochondrial genome to assess concordance among the ITS, nuclear, and mitochondrial inferences. Species partitions were empirically delimited from the ITS alignment using ASAP and bPTP. We also inferred a phylogeny for the L. polytropa clade using a four-marker dataset. ASAP species delimitations revealed up to 103 species in the L. polytropa clade, with 75 corresponding to the nominal taxon L. polytropa. Inferences from phylogenomic alignments generally supported that these represent evolutionarily independent lineages or species. Less than 10% of the candidate species were comprised of specimens from multiple continents. High levels of candidate species were recovered at local scales but generally with limited overlap across regions. Lecanora polytropa likely ranks as one of the largest species complexes of lichen-forming fungi known to date.

1. Introduction

Quantifying the full scope of species diversity is perhaps one of the largest undertakings of modern taxonomy [1]. Form and function have historically played fundamental roles in inferring species boundaries and relationships among species [2]. More recently, the incorporation of genetic information has highlighted limitations in past attempts to characterize the earth’s biological diversity based exclusively on phenotypes [3]. While the science of naming, describing, and classifying organisms has fallen largely to taxonomists, the implications of taxonomy extend to questions of the origin of novelty, speciation, symbiosis, conservation, and ecology [4,5,6]. Complementary, modern taxonomists use an ever-expanding toolbox from other disciplines for systematic and taxonomic research, integrating traditional approaches with genetic and ecological data, computational modeling, and empirical species delimitation [7,8,9,10].
As with other organismal groups, this integrative taxonomic approach has transformed our understanding of fungal diversity—challenging current taxonomy at multiple levels and highlighting rampant unrecognized diversity [11,12]. Though underestimates and challenges in establishing accurate inventories of fungal diversity may be expected due to the mostly cryptic lifestyle of fungi, reconsidering diversity in well-known, conspicuous fungi has also led to the realization that species diversity in these groups may also be mischaracterized [13,14,15,16]. Traditionally, and often due to lack of available genetic data, fungal species boundaries have been delimited based largely on morphological characteristics as they can be incredibly practical, especially in the field. However, robust phenotype-based species delimitations in fungi are often confounded by the existence of morphologically similar species in which intra- and inter-specific variations overlap regarding some characters that have been traditionally used to separate taxa [17].
Factors relating to fungi in obligate symbioses add another layer of complexity in delimiting species-level lineages. In lichen-forming fungi, stable interactions among multiple symbionts—the fungal nutritional specialists that acquire fixed carbon from photoautotrophs and are supported by other microbes [18,19,20,21]—result in a persistent thallus which can be considered a symbiotic phenotype that results only from the interaction of unrelated organisms [22]. While the symbiotic phenotype of lichens can provide important insight into species boundaries of lichen-forming fungi—the mycobiont, interactions of the mycobiont with differing suites of microorganisms may result in diverse symbiotic phenotypes of the same fungal species [23,24]. In other cases, carefully considering the lichen-associated symbionts can provide crucial insight for resolving otherwise apparently cryptic fungal species [25]. While the complexities of lichen symbioses and their implications for interpreting lichen-forming fungal diversity have recently been more fully recognized, commensurate advances in practical approaches to empirically delimiting lichen-forming fungi have also occurred [26].
Ongoing advances in sequencing technologies have been key in revealing the presence of previously unknown species diversity in fungi [12], opening the gates for research into lichen symbioses that would have previously been impossible [27,28,29,30]. Research focusing on long-standing, well-known problematic groups has been shown to lead to crucial resolution. For example, detailed investigations of a widespread montane, tropical macrolichen representing a single nominal fungal species highlighted that the species diversity of the name-giving mycobiont had been underestimated by orders of magnitude [16,31]. The question remains, is this vast underestimate of species-level diversity in lichen-forming fungi an anomaly or a more widespread phenomenon across lichen fungi?
Lecanora polytropa (Hoffm.) Rabenh. (Figure 1)—the name-giving mycobiont for the associated symbiotic phenotype—is distributed across largely disjunct, intercontinental populations and occurs on siliceous rocks (especially granite) in montane, alpine, and arctic/Antarctic habitats (Figure 2). As the name suggests, L. polytropa, and related species, exhibit a wide range of morphological diversity (Figure 1). Morphological diversity across the name-giving nominal taxon has been recognized over the course of two centuries of research, with dozens of described forms or varieties. At the same time, limited studies have explored the relationship of L. polytropa with other Lecanora species [32,33,34,35,36], and the full extent of diversity within the ‘L. polytropa group’ itself remains largely uncharted. Delimiting species boundaries is made more difficult by the potential role that the environment plays in shaping morphology. Specimens of L. polytropa that grow near copper mines, for instance, often have their color change from the normal yellow-green to turquoise blue [37]. While careful phenotypic assessments have led to the description of some species within the ‘L. polytropa group’ [32,33], the notorious variability of this species group, coupled with the occurrence of intermediate morphotypes within the ‘L. polytropa group’ [38] highlights the pressing need for a thorough survey of this group.
To better understand the scale of unrecognized species-level diversity in this common, cosmopolitan lichen-forming fungal species, here we circumscribe candidate species within this complex using DNA sequence data. We sampled over 300 specimens collected across multiple, intercontinental populations to meet this aim. From these specimens, we generated DNA sequence data from the standard DNA barcode marker—the internal transcribed spacer region [39]. A subset of specimens was then selected to represent the genetic diversity observed from the ITS data to generate multi-locus and genome-scale datasets. Based on these data, we provide compelling evidence that species diversity within the ‘L. polytropa group’ is vastly underestimated, with the nominal taxon L. polytropa representing at least 70 candidate species based on current, limited sampling and additional unrecognized species diversity in other species within the ‘L. polytropa group’. In so doing, this study lays a transformative framework for future studies to characterize species diversity more fully within the ‘L. polytropa’ species complex.

2. Materials and Methods

2.1. Taxon Sampling

Our sampling targeted specimens in the ‘L. polytropa group’ [35,40]. Recently, the ‘L. polytropa group’ was found to be a major lineage within the provisionally named ‘MPRPS clade’ sensu [35] within Lecanoraceae (comprising Myriolecis, Protoparmeliopsis, Rhizoplaca, the ‘L. polytropa group’, Bryonora, and the “Lecanora” saligna group). For this study, sampling efforts focused on L. polytropa sensu lato populations in Asia, Europe, and North America, supplemented with ITS sequences from GenBank, including Lecanora chlorophaeodes Nyl. (n = 3), L. dispersoareolata (Schaer.) Lamy (1), L. fuscobrunnea Dodge & Baker (31), L. intricata (Ach.) Ach. (4), L. polytropa (22), L. solaris Yakovchenko & Davydov (7), L. cf. subcinctula (Nyl.) Th. Fr. (2), L. subintricata (Nyl.) Th. Fr. (19), Rhizoplaca aspidophora (Vain.) Redón (1), and unidentified sequences also recovered within the ‘polytropa group’ clade (5) (Figure 2). The genus Carbonea was shown to be closely related to the ‘polytropa group’ [40], and all 21 ITS sequences currently available on GenBank were included here. Ultimately, a total of 340 specimens were included (Supplementary File S1). We note that in several cases, multiple lichen thalli were selected from the same vouchered collection when the voucher included multiple, distinct L. polytropa sensu lato thalli. To assess the range of diversity at a limited geographic scale in this nominal taxon, ‘L. polytropa group’ specimens were relatively densely sampled from the La Sal Mountains—a sky island on the Colorado Plateau, Utah, USA. Other geographic regions were not sampled as densely, aiming to characterize broader-scale patterns of distributions of putative species-level lineages. We could not obtain fresh material from alpine/subalpine habitats in Africa, Antarctica, Australia, Central America, or New Zealand, and only limited sequences from South America were available [35].
To explore the potential for diagnostic phenotypic traits separating specimens for distinct candidate species-level lineages circumscribed in this study, we characterized (i) general growth forms, (ii) spore sizes, and (iii) secondary metabolites. Morphological characters—e.g., thallus characters, surface color/texture, apothecial disk color, apothecia margin, etc.—were assessed using an Olympus SZH dissecting microscope. Observations and measurements of ascospores were made in water with an Olympus BH-2 microscope, with multiple ascospores measured from at least two apothecia on each specimen. Chemical constituents were identified using thin-layer chromatography (TLC), following standard methods with solvent systems ‘C’ and ‘G’ [41,42].

2.2. DNA Extraction and Sequencing

Total genomic DNA was extracted from specimens collected for this study using the ZR fungal/bacterial DNA miniprep kit (Zymo Research, Irvine, CA, USA), the Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA), or the DNAsecure Plant Kit (Tiangen Biotech, Beijing, China). For all specimens, we attempted to generate sequence data from the internal transcribed spacer region (ITS)—the standard DNA barcoding marker for fungi [39]. For selected specimens representing the genetic diversity observed from the ITS dataset, for example, attempting to represent as many candidate species-level lineages as possible, including multiple representatives for candidate species, if available (see Results), we targeted four additional loci traditionally used in phylogenetic analyses of Lecanoraceae [40]—a portion of the nuclear large-subunit (nuLSU), a fragment of the gene encoding the mitochondrial small subunit (mtSSU), and fragments from two nuclear protein-coding loci, the RNA polymerase II subunit 1 (RPB1) and RNA polymerase II subunit 2 (RPB2). Temperature profiles for polymerase chain reaction (PCR) amplification for all loci follow previous studies [40]. PCR amplifications were performed using Ready-To-Go PCR Beads (GE Healthcare, Pittsburgh, PA, USA); or alternatively, in 25 μL reactions containing 12.5 μL 2 × Taq PCR Mix (Tiangen Biotech, Beijing, China), 0.5 μL of each primer, 10.5 μL ddH2 O and 1 μL of DNA. PCR products were visualized on 1% agarose gel and cleaned using ExoSAP-IT (USB, Cleveland, OH, USA), following the manufacturer’s recommendations. We sequenced complementary strands with the same primers used for PCR amplification, and sequencing reactions were performed using BigDye 3.1 (Applied Biosystems, Foster City, CA, USA). Products were run on an ABI 3730 automated sequencer (Applied Biosystems) at the DNA Sequencing Center at Brigham Young University, Provo, UT, USA.
Single marker and multi-locus approaches may be insufficient to robustly delimit species boundaries, particularly among closely related species [43,44]. Genome-scale data provide unprecedented insight into species boundaries and testing concordance among independent loci [45,46,47]. Therefore, in addition to Sanger sequencing of traditional phylogenetic markers, metagenomic reads were newly generated from 32 specimens from the ‘L. polytropa group’ using short-read shotgun sequencing [48]. Using metagenomic reads from the selected specimens, independent DNA datasets were assembled to investigate concordance among the species partitions inferred from the ITS markers with clades inferred from genome-scale nuclear and mitochondrial datasets [45]. Specimens for Illumina sequencing were selected based on genetic diversity initially observed from sampling in western North America; as the project expanded, we could not include specimens representing additional diversity observed in subsequent broad geographic sampling for Illumina sequencing. For the specimens selected for metagenomic high-throughput sequencing, total genomic DNA was extracted from a small portion of lichen thalli (comprised of the mycobiont, photobiont, and other associated microbes) using the E.Z.N.A. Plant DNA DS Mini Kit (Omega Bio-Tek, Inc., Norcross, GA, USA) and following the manufacturers’ recommendations. Total genomic DNA was prepared following the standard Illumina whole genome sequencing (WGS) library preparation process using Adaptive Focused Acoustics for shearing (Covaris, Sydney, Australia), followed by an AMPure cleanup step. The DNA was then processed with the NEBNext Ultra™ II End Repair/dA-Tailing Module end-repair and the NEBNext Ultra™ II Ligation Module (New England Biolabs, Ipswich, MA, USA) while using standard Illumina index primers. Libraries were pooled and sequenced with the HiSeq 2500 sequencer in high output mode at the DNA Sequencing Center, Brigham Young University, Provo, UT, USA, using 250 cycle paired-end (PE) reads.

2.3. Short-Read Processing and Data Assembly

Raw reads were trimmed using Trimmomatic v0.39 [49] to remove adapter and primer sequences and low-quality reads. Bases at the start and end of reads were trimmed when they had a quality below 3 and 10, respectively, and when the quality of 5-bp sliding windows was <20. All trimmed reads <36 bp were filtered out. We performed a de novo genome assembly using PE reads from L. polytropa specimen “Leavitt 16-650” using SPAdes [50]. To identify single-copy nuclear genes for phylogenomic reconstructions from the assembled mycobiont contigs, we used Benchmarking Universal Single-Copy Orthologs to extract up to 1438 gene regions (BUSCO; [51]. Assembled contigs were analyzed using the BUSCO pipeline implemented in the Cyverse.org Discovery Environment [52,53]. The Fungi Odb10 dataset was used to identify BUSCO genes from the assembled L. polytropa contigs. Exploratory BLAST searches and assessments of relative sequencing coverage were used to infer that the extracted BUSCO genes likely originated from the L. polytropa genome and not other co-occurring fungi. Partial and multi-copy BUSCOs were excluded, and the remaining filtered, single-copy BUSCO genes were used as targets for bait sequence capture using HybPiper [54] to extract these genes regions from each ‘L. polytropa group’ metagenomic sample (e.g., [48]). MAFFT [55] was used to generate alignments for individual BUSCO genes using the default parameters, and the alignment algorithm for each locus was chosen automatically by MAFFT. For each alignment, any sample which had an average completion (assembly length/target length) of <0.20 was removed. Genes with average coverage <75% across all BUSCO genes were removed. Phylogenetic gene trees were reconstructed using maximum likelihood (ML) as implemented by IQ-TREE [56]. The substitution model used for each tree was selected using ModelFinder [57]. To assemble genome-scale data from the mycobiont mitochondrial genome, we identified mitochondrial contigs from the SPAdes assembly using BLAST comparisons [58]. The three longest mitochondrial contigs—representing a total of 87.9 Kb—were used as targets in RealPhy v1.12 [59]. For the RealPhy assembly, we used the following parameters to generate the mitochondrial genome alignment: -readLength 100; –perBaseCov 5; –gapThreshold 0.2.

2.4. Candidate Species Delimitation Using the Standard Fungal DNA Barcode

Initial candidate species partitions for the L. polotropa group were inferred using Assemble Species by Automatic Partitioning (ASAP) [60] based on the multiple sequence alignment of the standard fungal DNA barcode—ITS [39]. ASAP circumscribes species partitions using an implementation of a hierarchal clustering algorithm based on pairwise genetic distances from single-locus sequence alignments [60]. The pairwise genetic distances are used to build a list of partitions ranked by a score, computed using the probability of groups to define panmictic species. ASAP provides an objective approach to circumscribe relevant species hypotheses as a first step in the process of species delimitation. ITS sequences generated for this study were combined with those GenBank and aligned using the program MAFFT v7 [55,61]. We implemented the G-INS-i alignment algorithm and ‘1PAM/K = 2’ scoring matrix with an offset value of 0.1, the ‘unalignlevel’ = 0.2, and the remaining parameters were set to default values. The multiple sequence alignment was analyzed using the ASAP Web Server (https://bioinfo.mnhn.fr/abi/public/asap/, accessed on 27 January 2022), with the ‘asap-score’ considered to select the optimal number of species partitions [60]. In addition to ASAP species delimitations, we also implemented a single-locus tree-based species delimitation method—the Bayesian implementation of the Poisson tree process model (bPTP) [62]. A maximum-likelihood (ML) tree was inferred from the ITS alignment using IQ-TREE [56], which was subsequently analyzed using bPTP with 1,000,000 generations and a burn-in of 10%.

2.5. Phylogenetic Analyses and Tests of Genomic Concordance

To assess the monophyly of the ‘L. polytropa group’ within Lecanoraceae, ITS, and mtSSU sequences from the present study were combined with the relatively comprehensive “2-locus dataset”, including 251 OTUs representing 150 species, originally reported in [40]. Sequences were aligned using MAFFT v7 [55,61], implementing the same parameters described above for ITS. To minimize ambiguities in multiple sequence alignments (MSA), subsequent MSA and phylogenetic reconstructions were restricted to members of the ‘L. polytropa group’ and specimens representing the two putative sister lineages inferred from the family-wide, “2-locus dataset”—Carbonea species and those representing L. subintricata (see Results). For the ‘L. polytropa group’ and putative sister groups, the ITS, nuLSU, mtSSU, RPB1, and RPB2 sequences, were aligned using MAFFT v7 as described above. For both the nuLSU and mtSSU alignments, ambiguously aligned regions were excluded using the Gblocks webserver [63], implementing the options for a less stringent selection. A ML tree for the ‘L. polytropa group’ was inferred from the concatenated five-marker alignment using IQ-TREE [56]. The concatenated alignment was partitioned by loci, with substitution models selected using ModelFinder [57] and nodal support assessed using 2000 ultra-fast bootstrap replicates [64].
To assess concordance between the candidate species inferred from the standard fungal barcode (ITS) and genome-scale data, we compared ASAP partitions with (i) clades inferred from concatenated single-copy nuclear markers spanning 2.28 Mb, (ii) clades inferred from a mitochondrial alignment spanning 65.5 Kb, and (iii) a phylogenomic approach to species delimitation. Our genome-scale sampling only represented a subset of the candidate species inferred from the ITS data, and comparisons were limited to the 12 ASAP partitions that also had representative samples with short-read data (32 specimens). We reconstructed a phylogeny from a supermatrix comprised of the 1209 single-copy BUSCO markers. Concatenation approaches provide accurate inferences under a range of conditions [65]. We used IQ-TREE v1.6.9 to generate a ML tree from the concatenated BUSCO supermatrix, with nodal support assessed using 1000 ultra-fast bootstrap replicates [64]. We also used SODA [66], an ultra-fast and relatively accurate method for species delimitation, to assess candidate species boundaries inferred from the BUSCO data. SODA uses frequencies of quartet topologies to determine if each branch in a guide tree inferred from gene trees (1209 BUSCO topologies) is likely to have a positive length. It uses the results to infer a new species tree that defines species boundaries. We ran SODA, implemented in ASTRAL [66], with a p-value cut-off of 0.001. The mitochondrial topology was reconstructed from the RealPhy alignment using IQ-TREE [56], with the substitution model selected using ModelFinder [57], and nodal support assessed using 1000 ultra-fast bootstrap replicates [64].

3. Results

3.1. Sequence Data

Newly generated ITS, nuLSU, RPB1, RPB2, and mtSSU sequences are deposited in GenBank under accession numbers ON179980–ON180462 and ON217582–ON217793. Illumina short reads from the 32 ‘L. polytropa group’ specimens are available in the NCBI Short Read Archive (PRJNA823672). Of the 1438 BUSCO genes searched, 1209 complete, single-copy BUSCO genes were recovered (93.5% of all BUSCO groups). The concatenated alignments of the 1209 nuclear BUSCO markers spanned a total of 2.28 Mb. The REALPHY genome skimming approach for generating mitochondrial data resulted in an alignment of 65456 bp. The concatenated five-marker dataset comprised 3891 aligned nucleotide position characters—ITS (n = 380; 603 bp MSA), nuLSU (n = 122; 843 bp MSA [ambiguous sites removed]), RPB1 (n = 117; 816 bp MSA), RPB2 (n = 105; 828 bp MSA), and mtSSU (n = 112; 801 bp MSA [ambiguous sites removed]).

3.2. Candidate Species Inferred Using the Standard DNA Barcode (ITS)

The best-scoring ASAP species partitions delimited between 62–103 candidate species within the ‘L. polytropa group’—the 103- and 102-species models had the best asap-scores (Supplementary File S2 and S3). Specimens identified as L. polytropa represented up to 75 distinct species partitions in the ASAP analyses of the ITS alignment. Multiple candidate species were also recovered within the nominal taxa L. concolor (up to 2 species partitions), L. dispersoareolata (2), L. intricata (6), L. solaris (2), L. sommervelli (2), and L. subintricata (5). Six species partitions comprised Carbonea sequences—the currently unsettled sister-clade to the ‘L. polytropa group’. Based on current sampling, most candidate species partitions are found in geographically limited areas, with only seven species partitions including samples from multiple continents (Supplementary File S1). The tree-based bPTP species delimitation model resulted in 73 candidate species within the ‘L. polytropa group’, which were concordant with most of the ASAP species partitions, with ten of the 73 bPTP candidate species combining multiple ASAP partitions from the best-supported models (Supplementary File S4).

3.3. Multi-Locus Phylogenetic Inference and Phylogenetic Concordance among Data Sets

In the combined ITS/mtSSU topology, the ‘L. polytropa group’ was recovered as monophyletic with strong bootstrap (BS) support within the ‘MPRPS clade’ sensu [35]. In the combined ITS/mtSSU topology, Carbonea specimens were recovered as monophyletic within the ‘L. polytropa group’; and L. subintricata specimens were recovered as sisters to the entire ‘L. polytropa group’ with weak support (58% BS; Supplementary Files S5 and S6). In the 5-marker ‘L. polytropa group’ topology, including both Carbonea and L. subintricata specimens, a midpoint rooted topology showed Carbonea specimens as sisters to the ‘L. polytropa group’, with L. subintricata specimens nested within. Given the uncertainty of the sister clade to the ‘L. polytropa group’, we opted to show the midpoint root topology (Figure 3), rather than rooted with L. subintricata specimens as inferred from the “2-locus dataset”.
In the 5-marker ‘L. polytropa group’ topology, the vast majority of candidate species partitions inferred from the ITS MSA using ASAP were recovered as reciprocally monophyletic (Figure 3; Supplementary Files S7 and S8). However, the candidate species partitions inferred from the ITS MSA were not recovered as reciprocally monophyletic in any of the topologies inferred from four traditional markers individually (nuLSU, RPB1, RPB2, and mtSSU), although some species partitions were monophyletic in some single-gene topologies. Specimens representing previously described species that occur within the L. polytropa clade were generally recovered as monophyletic—L. concolor, L. dispersoareolata, L. fuscobrunnea, L. solaris, L. sommervelli, and L. subintricata; although multiple candidate species were inferred from ITS sequence data in most of these nominal taxa (Figure 3; Supplementary File S7). In contrast, sequences representing L. chlorophaeodes and L. intricata were not recovered as monophyletic.
The limited nuclear phylogenomic sampling comprising specimens from the western USA unambiguously supported the ASAP candidate species partitions represented by short-read data (Figure 4A). The mitochondrial phylogenomic data also generally supported the ASAP species partitions with one instance of putative mitochondrial introgression (Figure 4B). The SODA species delimitation analyses based on 1209 BUSCO gene topologies consistently further subdivided the candidate species partitions delimited from ITS sequence data using ASAP, delimiting 20 species from the specimens represented by genome-scale data using SODA, in contrast to the 12 ASAP species partitions (Figure 4A).

3.4. Assessing Morphological Concordance with Species Partitions

ASAP partitions represented by multiple specimens often comprised polymorphic morphologies. No phenotypic characters were observed that consistently diagnosed any of the ASAP partitions comprised of three or more specimens, although some characters or combinations of characters (e.g., spore size, growth form, secondary metabolite variation, and morphology) loosely corresponded with some distinct species partitions. In a limited number of cases, secondary metabolite concentrations qualitatively varied among specimens in different ASAP partitions. However, given the uneven representation of specimens per candidate species, with many candidate species represented by a very limited number of species, no quantitative comparisons were made. Lecanora stenotropa Nyl. is morphologically similar to L. polytropa but with a more brownish-green thallus and smaller, narrowly ellipsoid spores. However, based on anatomical observations, no L. stenotropa specimens from western North America were sampled, and the relationship of this taxon to other species-level clades in the L. polytropa group remains unknown.

4. Discussion

Here, molecular sequence data reveal that the well-known, problematic nominal taxon Lecanora polytropa likely ranks as one of the largest species complexes of lichen-forming fungi known to date. Based on species delimitations using alignments of ITS sequences, as many as 103 species partitions were circumscribed with the L. polytropa clade using ASAP, including ca. 75 candidate species identified as L. polytropa and multiple candidate species in other formally described species in the clade. For comparison, perhaps the most species-rich group is the nominal taxon Dictyonema glabratum, likely comprising at least 400 distinct species-level lineages [16]. Most other nominal taxa known to mask multiple species-level lineages comprise perhaps up to ten distinct species-level lineages [29,67,68,69,70,71,72]. Overall, these findings corroborated the perspective that some currently circumscribed conspicuous, well-known lichens may harbor spectacular levels of unrecognized species-level diversity [16]. While phenotypic variation within L. poytropa s. lat. is well documented and has historically been interpreted to circumscribe dozens of forms or varieties, a widely accepted workable taxonomy has been elusive for this group over the past two centuries. Below we discuss important implications relative to the likely extent of the species-level diversity in the L. polytropa clade.
We provide the most comprehensive sampling to date for members within the L. polytropa clade, with over 300 sampled specimens collected from populations across the globe (Figure 2). However, despite the more than 70-fold increase in putative species in the nominal taxon L. polytropa, the present sampling effort is likely insufficient to fully capture species-level diversity in this group. For this study, our densest sampling effort targeted a limited number of mountain ranges in western China and the southwestern USA, complemented by opportunistic sampling in other regions. At local scales, for example, individual mountain ranges, high levels of candidate species were recovered (Supplementary File S1). For example, in the La Sal Mountains, in southeastern Utah, USA, up to 14 candidate species were observed. However, candidate species observed on other nearby mountains, e.g., Uinta Mountains (UT, USA), Beartooth Plateau (MT/WY, USA), Saguache Range (CO, USA), had only partly overlapping suites of candidate species (Figure 3). Hence, we speculate that a substantial portion of species diversity in the L. polytropa clade remains undiscovered even in the regions most densely sampled for this study. The speculated underestimate of species diversity in the group is further confounded by the fact that vast regions remain poorly sampled, for example, Europe, Northern Asia, and isolated populations such as those in Australasia (Figure 2). Predictive models for estimating species richness, such as that implemented by Lucking et al. [16], can help in providing quantitative estimates of diversity and directing future sampling.
Given the apparent species richness in the group and the unsettled species boundaries, effective strategies must be developed to robustly delimit species boundaries [26,46]. Based on our results, it appears that the ITS coupled with sequence-based species delimitation approaches, such as ASAP [60], is appropriate as a first pass investigation into species boundaries in the L. polytropa clade. Following an initial screening of ITS diversity, targeted genome-scale sequencing will be essential to robustly test evolutionary independence among candidate species-level lineages [46,73]. While our multi-locus phylogenetic reconstruction generally resulted in distinct, well-supported clades coinciding with candidate species inferred from ITS sequence data (Figure 3), the multi-locus dataset was insufficient to resolve many backbone-level relationships nor fully corroborate the independence of candidate species. The genome skimming method implemented here provided a wide range of phylogenomic data that also supported rampant unrecognized diversity in the L. polytropa clade. By investigating concordance among different DNA datasets [45], e.g., the ITS (standard fungal barcode) and inferences from the mitochondrial and nuclear genomes, we observed a general pattern that suggests long-term, evolutionary independence among the candidate species within the L. polytropa clade. Specifically, data from whole-genome sequencing of a subset of the L. polytropa group samples support the inference of rampant, unrecognized species-level diversity (Figure 3 and Figure 4), largely congruent with ASAP delimitations. These results suggest that the candidate species in the L. polytropa group are not an artifact of limitations due to single-locus species delimitation methods or intragenomic variation within the multi-copy ribosomal tandem repeat [74]. Rather, there are most likely legitimate reproductive barriers resulting in the observed phylogenetic structure and inferred species partitions.
Although genome skimming approaches, such as those implemented here, can provide crucial genome-scale data [48,75,76], more cost- and time-efficient methods, such as target enrichment sequencing [77,78] or restriction-site associated sequencing (RADseq) [26,47,79] might be more appropriate for this clade given the high number of samples that will likely need to be included. Ultimately, genome-scale data can also be used to reconstruct a robust phylogeny for this group to explore additional evolutionary, phylogeographic, and taxonomic questions. Arguably, genetic data alone is not sufficient for robust delimitations of species boundaries and is best when evaluated in conjunction with other information [80]. While these candidate species provide a powerful framework for taxonomic hypothesis testing, we ultimately support a hypothesis-based, integrative approach to species delimitation [26]. In this study, biased specimen sampling, small sample sizes for most candidate species, and pending detailed anatomical assessments all limit the robustness of the species boundaries inferred here. We hope that these species hypotheses serve as a starting point for integrated approaches required for robust species boundaries and good taxonomy [80,81].
The spectacular species-level diversity and unexpected biogeographic patterns highlight the complex speciation and phylogeographic history of the L. polytropa clade (Figure 3). Based on current, limited sampling, we found a mix of a few cosmopolitan lineages and many local/regional endemics. Other lichen-forming fungi have shown similar patterns of complex intercontinental species distributions [82,83,84,85]. Presently, little is known of evolutionary processes that may give rise to reproductive isolation and ultimately speciation among populations in the L. polytropa clade. Our results indicate frequent, long-distance dispersal throughout the diversification history of this clade (Figure 3). While most candidate species inferred here were not found across broad, intercontinental distributions, ca. 10% comprised specimens from multiple continents. Future work should attempt to elucidate species distributions and different ecological niches among candidate species and factors influencing ecological specialization, dispersal limitations, etc. For example, some morphologically similar species, such as L. microloba, appear to have specialized ecological niches and distributions [36], while our data suggest that others have wide ecological amplitude, resulting in more broadly distributed species. Our hope is that this study provides the impetus for using members of the L. polytropa clade as a model to explore phylogeography and speciation in symbiotic fungi.
The results of this study have complex taxonomic implications for L. polytropa and closely related species. Given the long, unsettled taxonomic history of L. polytropa and associated forms and varieties, a careful revision of the complex taxonomy and available type material will be required to identify which candidate species coincide with L. polytropa and other formally described species before naming new species can proceed [26]. This effort may be confounded by the lack of consistent diagnostic features separating evolutionarily distinct lineages, potentially due, in part, to symbiotic interactions [23,86]. Initial phenotypic investigations here failed to reveal consistent taxonomically diagnostic traits, e.g., spore size (Supplementary File S5), corroborating a subset of the distinct candidate species. In some cases, intraspecific phenotype variation/plasticity relating to the thallus morphology, secondary metabolite variation, and spore size was observed. Lecanora fuscobrunnea, an Antarctic endemic with distinct apothecia morphology, was found to be closely related to several L. polytropa specimens from western North America (Figure 3) and were inferred to be conspecific in several ASAP species partition models (Supplementary File S2). However, given the sparse sample size for most candidate species circumscribed in this study and the limited number of specimens examined, this conclusion should be held only tentatively. In contrast, our results show/corroborate that several morphologically distinct taxa also belong to the L. polytropa clade, including L. chlorophaeodes, L. dispersoareolata, L. fuscobrunnea, L. intricata, L. solaris, L. cf. subcinctula, and L. subintricata. While multiple candidate species also occur in a number of these other nominal taxa, the morphologically distinct groups are generally separated from specimens representing L. polytropa (Figure 3). Importantly, we propose that future species descriptions include genetic information in the formal description, as the highly variable morphology is often non-diagnostic, making classification based simply on morphology largely unreliable, as evidenced by two centuries of unresolved species boundaries in the L. polytropa clade. Despite the indisputable merits of molecular sequence data in taxonomy, these data are only scarcely used for the formal description of taxa. Novel approaches for DNA-based diagnoses, including diagnostic nucleotide combinations in DNA sequence alignments [87], can be used to provide formal diagnoses for these challenging groups.
Our study lays the groundwork for elucidating important evolutionary insight and formally recognizing undescribed species diversity and taxonomic diversity for this important cosmopolitan clade of symbiotic fungi.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof8050490/s1. Supplementary File S1: List of specimens from the L. polytropa group included in this study. Data include voucher number and herbarium for specimens collected for this study and GenBank accession numbers for sequences from GenBank; DNA number used for labels in multiple sequences alignments; taxonomic assignments for all specimens, including both candidate ASAP species partition codes (shown in Figure 3) and morphology-based specimen identifications; geographic information for specimens (note: geographic origin for sequences from GenBank is incomplete; and in a number of cases, latitude and longitude were inferred [shown in red text]); and GenBank accession numbers for the internal transcribed spacer region (ITS), the nuclear ribosomal large subunit (nuLSU), RPB1, RPB2, mitochondrial small subunit (mtSSU), and short reads for selected specimens. Supplementary File S2: ASAP scores and rankings for the ‘L. polytropa group’ inferred from a multiple sequence alignment of the internal transcribed spacer region (ITS, standard fungal barcode marker) and comprising 380 sequences. Ten best partitions and coinciding asap-scores and ranking (panel on left); and graphical representation of all asap-scores (center and right panels), highlighting the ten best partitions shown in the panel. Supplementary File S3: Comparison of sample assignments to candidate species across the ten best asap partitions. Supplementary File S4: Results of the Bayesian implementation of the Poisson tree process (bPTP) species delimitation model. A maximum-likelihood (ML) tree was inferred from the ITS alignment using IQ-TREE, which was subsequently analyzed using bPTP with 1,000,000 generations and a burn-in of 10%—see trace file below. Supplementary File S5: Two-marker (ITS + mtSSU) maximum likelihood topology of Lecanoraceae (summarized in Figure 2). Specimen sampling includes L. polytropa samples compiled for this study combined with a broad, Lecanoraceae-wide sampling reported in [40]. Bootstrap support is indicated at each node and the ‘L. polytropa group’ is shown in red branches, with L. subintricata specimens and Carbonea specimens with blue branches. Supplementary File S6: Concatenated two-marker (ITS + mtSSU) alignment representing Lecanoraceae to assess monophyly of the L. polytropa group. The two-marked dataset comprised 630 specimens and spanned 1337 aligned nucleotide position characters: ITS, 781 bp; mtSSU, 556 bp. Supplementary File S7: Midpoint-rooted, five-marker (ITS, nuLSU, RPB1, RPB2, mtSSU) maximum likelihood topology of the ‘L. polytropa group’, including Carbonea species and L. subintricata. Bootstrap values greater than 50% are shown for each node. This topology is summarized in Figure 3, and locality info for each specimen can be found in Supplementary File S1. Supplementary File S8: Concatenated five-marker (ITS, nuLSU, RPB1, RPB2, mtSSU) alignment representing the L. polytropa group and putative sister clades. The five-marker dataset comprised 382 specimens and spanned 3891 aligned nucleotide position characters—ITS (n = 380; 603 bp MSA), nuLSU (n = 122; 843 bp MSA [ambiguous sites removed]), RPB1 (n = 117; 816 bp MSA), RPB2 (n = 105; 828 bp MSA), and mtSSU (n = 112; 801 bp MSA [ambiguous sites removed]). The ITS portion of the concatenated alignment (comprising alignment positions 1–603) was used for the ASAP and BPTP species delimitations. Supplementary File S9: Concatenated alignment of 1209 nuclear BUSCO markers, spanning a total of 2.28 Mb. Of the 1438 BUSCO genes searched, 1209 complete, single-copy BUSCO genes were recovered (93.5% of all BUSCO groups). Supplementary File S10: The REALPHY genome skimming approach for generating mitochondrial data resulted in an alignment of 65,456 bp.

Author Contributions

Conceptualization, Y.Z., J.C. and S.D.L.; methodology, Y.Z., J.C., J.J., R.T.M., L.W. and S.D.L.; software, J.J.; validation, Y.Z., J.C., J.J., R.T.M., L.W. and S.D.L.; formal analysis, J.C., J.J. and S.D.L.; investigation, Y.Z., J.C., J.J., R.T.M., L.W. and S.D.L.; resources, Y.Z., R.T.M., L.W. and S.D.L.; data curation, Y.Z. and S.D.L.; writing—original draft preparation, J.C. and S.D.L.; writing—review and editing, Y.Z., J.C., J.J., R.T.M., L.W. and S.D.L.; visualization, J.C. and S.D.L.; supervision, L.W. and S.D.L.; project administration, Y.Z. and S.D.L.; funding acquisition, Y.Z., R.T.M., L.W. and S.D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Nos 31970022, 31750001), Canyonlands Natural History Association (Moab, UT, USA), Second Tibetan Plateau Scientific Expedition and Research (STEP) program [No. 2019QZKK0503], and the M.L. Bean Museum of Life Sciences at Brigham Young University (Provo, UT, USA).

Data Availability Statement

Newly generated ITS, nuLSU, RPB1, RPB2, and mtSSU sequences are deposited in GenBank under accession numbers ON179980–ON180462 and ON217582–ON217793. Illumina short reads from the 32 ‘L. polytropa group’ specimens are available in the NCBI Short Read Archive (PRJNA823672). All alignments are available provided as Supplementary Files.

Acknowledgments

We also thank Barb Smith (Moab District, Manti–La Sal National Forest, Moab, UT, USA), Larry St. Clair (Brigham Young University, Provo, UT, USA), Taylor Hawes, and the Leavitt family for assistance with fieldwork. We thank two anonymous referees and Silke Werth (Academic Editor) for the invaluable feedback that greatly improved this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Blaxter, M.; Archibald, J.M.; Childers, A.K.; Coddington, J.A.; Crandall, K.A.; Di Palma, F.; Durbin, R.; Edwards, S.V.; Graves, J.A.M.; Hackett, K.J.; et al. Why sequence all eukaryotes? Proc. Natl. Acad. Sci. USA 2022, 119, e2115636118. [Google Scholar] [CrossRef] [PubMed]
  2. Ogle, W. (Ed.) Aristotle: On the Parts of Animals; Kegan Paul, Trench & Co.: London, UK, 1882. [Google Scholar]
  3. Burki, F.; Roger, A.J.; Brown, M.W.; Simpson, A.G. The New Tree of Eukaryotes. Trends Ecol. Evol. 2020, 35, 43–55. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Bennett, B.C.; Balick, M.J. Does the name really matter? The importance of botanical nomenclature and plant taxonomy in biomedical research. J. Ethnopharmacol. 2014, 152, 387–392. [Google Scholar] [CrossRef] [PubMed]
  5. Borowitzka, M.A. Systematics, taxonomy and species names: Do They Matter? In The Physiology of Microalgae; Borowitzka, M.A., Beardall, J., Raven, J.A., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 655–681. [Google Scholar] [CrossRef]
  6. Barron, E.S.; Sthultz, C.M.; Hurley, D.; Pringle, A. Names matter: Interdisciplinary research on taxonomy and nomenclature for ecosystem management. Prog. Phys. Geogr. Earth Environ. 2015, 39, 640–660. [Google Scholar] [CrossRef] [Green Version]
  7. Schlick-Steiner, B.C.; Steiner, F.M.; Seifert, B.; Stauffer, C.; Christian, E.; Crozier, R.H. Integrative Taxonomy: A Multisource Approach to Exploring Biodiversity. Annu. Rev. Èntomol. 2010, 55, 421–438. [Google Scholar] [CrossRef]
  8. Dayrat, B. Towards integrative taxonomy. Biol. J. Linn. Soc. 2005, 85, 407–415. [Google Scholar] [CrossRef]
  9. Smith, M.L.; Carstens, B.C. Process-based species delimitation leads to identification of more biologically relevant species. Evolution 2020, 74, 216–229. [Google Scholar] [CrossRef]
  10. Camargo, A.; Sites, J. Species Delimitation: A Decade After the Renaissance. In The Species Problem-Ongoing Issues; Pavlinov, I., Ed.; InTech: London, UK, 2013. [Google Scholar] [CrossRef] [Green Version]
  11. Nilsson, R.H.; Larsson, K.-H.; Taylor, A.F.S.; Bengtsson-Palme, J.; Jeppesen, T.S.; Schigel, D.; Kennedy, P.; Picard, K.; Glöckner, F.O.; Tedersoo, L.; et al. The UNITE database for molecular identification of fungi: Handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 2019, 47, D259–D264. [Google Scholar] [CrossRef]
  12. Baldrian, P.; Větrovský, T.; Lepinay, C.; Kohout, P. High-throughput sequencing view on the magnitude of global fungal diversity. Fungal Divers. 2021. [Google Scholar] [CrossRef]
  13. Salgado-Salazar, C.; Rossman, A.Y.; Chaverri, P. Not as Ubiquitous as We Thought: Taxonomic Crypsis, Hidden Diversity and Cryptic Speciation in the Cosmopolitan Fungus Thelonectria discophora (Nectriaceae, Hypocreales, Ascomycota). PLoS ONE 2013, 8, e76737. [Google Scholar] [CrossRef] [Green Version]
  14. Chethana, K.W.T.; Manawasinghe, I.S.; Hurdeal, V.G.; Bhunjun, C.S.; Appadoo, M.A.; Gentekaki, E.; Raspé, O.; Promputtha, I.; Hyde, K.D. What are fungal species and how to delineate them? Fungal Divers. 2021, 109, 1–25. [Google Scholar] [CrossRef]
  15. Kobmoo, N.; Mongkolsamrit, S.; Arnamnart, N.; Luangsa-Ard, J.J.; Giraud, T. Population genomics revealed cryptic species within host-specific zombie-ant fungi (Ophiocordyceps unilateralis). Mol. Phylogenet. Evol. 2019, 140, 106580. [Google Scholar] [CrossRef] [PubMed]
  16. Lücking, R.; Dal-Forno, M.; Sikaroodi, M.; Gillevet, P.M.; Bungartz, F.; Moncada, B.; Yánez-Ayabaca, A.; Chaves, J.L.; Coca, L.F.; Lawrey, J.D. A single macrolichen constitutes hundreds of unrecognized species. Proc. Natl. Acad. Sci. USA 2014, 111, 11091–11096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Cao, B.; Haelewaters, D.; Schoutteten, N.; Begerow, D.; Boekhout, T.; Giachini, A.J.; Gorjón, S.P.; Gunde-Cimerman, N.; Hyde, K.D.; Kemler, M.; et al. Delimiting species in Basidiomycota: A review. Fungal Divers. 2021, 109, 181–237. [Google Scholar] [CrossRef]
  18. Peksa, O.; Gebouská, T.; Škvorová, Z.; Vančurová, L.; Škaloud, P. The guilds in green algal lichens—An insight into the life of terrestrial symbiotic communities. FEMS Microbiol. Ecol. 2022, 98, fiac008. [Google Scholar] [CrossRef]
  19. Spribille, T.; Tuovinen, V.; Resl, P.; Vanderpool, D.; Wolinski, H.; Aime, M.C.; Schneider, K.; Stabentheiner, E.; Toome-Heller, M.; Thor, G.; et al. Basidiomycete yeasts in the cortex of ascomycete macrolichens. Science 2016, 353, 488–492. [Google Scholar] [CrossRef] [Green Version]
  20. 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]
  21. Rolshausen, G.; Hallman, U.; Grande, F.D.; Otte, J.; Knudsen, K.; Schmitt, I. Expanding the mutualistic niche: Parallel symbiont turnover along climatic gradients. Proc. R. Soc. B 2020, 287, 20192311. [Google Scholar] [CrossRef]
  22. Honegger, R. The Symbiotic Phenotype of Lichen-Forming Ascomycetes. In Fungal Associations; Hock, B., Ed.; Springer: Berlin/Heidelberg, Germany, 2001; Volume 9, pp. 165–188. [Google Scholar]
  23. Spribille, T. Relative symbiont input and the lichen symbiotic outcome. Curr. Opin. Plant Biol. 2018, 44, 57–63. [Google Scholar] [CrossRef]
  24. Armaleo, D.; Clerc, P. Lichen chimeras: DNA analysis suggests that one fungus forms two morphotypes. Exp. Mycol. 1991, 15, 1–10. [Google Scholar] [CrossRef]
  25. Schneider, K.; Resl, P.; Spribille, T. Escape from the cryptic species trap: Lichen evolution on both sides of a cyanobacterial acquisition event. Mol. Ecol. 2016, 25, 3453–3468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lücking, R.; Leavitt, S.D.; Hawksworth, D.L. Species in lichen-forming fungi: Balancing between conceptual and practical considerations, and between phenotype and phylogenomics. Fungal Divers. 2021, 109, 99–154. [Google Scholar] [CrossRef]
  27. Gilbert, S.F.; Sapp, J.; Tauber, A.I. A Symbiotic View of Life: We Have Never Been Individuals. Q. Rev. Biol. 2012, 87, 325–341. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Tuovinen, V.; Ekman, S.; Thor, G.; Vanderpool, D.; Spribille, T.; Johannesson, H. Two Basidiomycete Fungi in the Cortex of Wolf Lichens. Curr. Biol. 2019, 29, 476–483.e5. [Google Scholar] [CrossRef] [Green Version]
  29. Magain, N.; Tniong, C.; Goward, T.; Niu, D.; Goffinet, B.; Sérusiaux, E.; Vitikainen, O.; Lutzoni, F.; Miadlikowska, J. Species delimitation at a global scale reveals high species richness with complex biogeography and patterns of symbiont association in Peltigera section Peltigera (lichenized Ascomycota: Lecanoromycetes). TAXON 2018, 67, 836–870. [Google Scholar] [CrossRef] [Green Version]
  30. McKenzie, S.K.; Walston, R.F.; Allen, J.L. Complete, high-quality genomes from long-read metagenomic sequencing of two wolf lichen thalli reveals enigmatic genome architecture. Genomics 2020, 112, 3150–3156. [Google Scholar] [CrossRef]
  31. Lücking, R.; Forno, M.D.; Moncada, B.; Coca, L.F.; Vargas-Mendoza, L.Y.; Aptroot, A.; Arias, L.J.; Besal, B.; Bungartz, F.; Cabrera-Amaya, D.M.; et al. Turbo-taxonomy to assemble a megadiverse lichen genus: Seventy new species of Cora (Basidiomycota: Agaricales: Hygrophoraceae), honouring David Leslie Hawksworth’s seventieth birthday. Fungal Divers. 2017, 84, 139–207. [Google Scholar] [CrossRef]
  32. Lagreca, S.; Lumbsch, H.T. Three Species of Lecanora New to North America, With Notes on Other Poorly Known Lecanoroid Lichens. Bryologist 2001, 104, 204–211. [Google Scholar] [CrossRef]
  33. Yakovchenko, L.S.; Davydov, E.A.; Ohmura, Y.; Printzen, C. The phylogenetic position of species of Lecanora s. l. containing calycin and usnic acid, with the description of Lecanora solaris Yakovchenko & Davydov sp. nov. Lichenologist 2019, 51, 147–156. [Google Scholar] [CrossRef]
  34. Obermayer, W.; Poelt, J. Contributions to the Knowledge of the Lichen Flora of the Himalayas III. On Lecanora Somervellii Paulson (Lichenized Ascomycotina, Lecanoraceae). Lichenologist 1992, 24, 111–117. [Google Scholar] [CrossRef]
  35. Medeiros, I.D.; Mazur, E.; Miadlikowska, J.; Flakus, A.; Rodriguez-Flakus, P.; Pardo-De la Hoz, C.J.; Cieślak, E.; Śliwa, L.; Lutzoni, F. Turnover of Lecanoroid Mycobionts and Their Trebouxia Photobionts Along an Elevation Gradient in Bolivia Highlights the Role of Environment in Structuring the Lichen Symbiosis. Front. Microbiol. 2021, 12, 774839. [Google Scholar] [CrossRef] [PubMed]
  36. Śliwa, L.; Flakus, A. Lecanora microloba, a new saxicolous species from Poland. Lichenologist 2011, 43, 1–6. [Google Scholar] [CrossRef]
  37. Purvis, O.W.; Pawlik-Skowrońska, B.; Cressey, G.; Jones, G.C.; Kearsley, A.; Spratt, J. Mineral phases and element composition of the copper hyperaccumulator lichen Lecanora polytropa. Miner. Mag. 2008, 72, 607–616. [Google Scholar] [CrossRef]
  38. De La Rosa, I.N.; Passo, A.; Rodríguez, J.M.; Chiapella, J.O.; Messuti, M.I. A new species and new records of Lecanora (Lecanoraceae, lichenized Ascomycota) with usnic acid from the Antarctic region. Phytotaxa 2016, 261, 185–193. [Google Scholar] [CrossRef]
  39. Schoch, C.L.; Seifert, K.A.; Huhndorf, S.; Robert, V.; Spouge, J.L.; Levesque, C.A.; Chen, W.; Fungal Barcoding Consortium; Fungal Barcoding Consortium Author List; Bolchacova, E.; et al. Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc. Natl. Acad. Sci. USA 2012, 109, 6241–6246. [Google Scholar] [CrossRef] [Green Version]
  40. Zhao, X.; Leavitt, S.D.; Zhao, Z.T.; Zhang, L.L.; Arup, U.; Grube, M.; Pérez-Ortega, S.; Printzen, C.; Śliwa, L.; Kraichak, E.; et al. Towards a revised generic classification of lecanoroid lichens (Lecanoraceae, Ascomycota) based on molecular, morphological and chemical evidence. Fungal Divers. 2016, 78, 293–304. [Google Scholar] [CrossRef]
  41. Culberson, C.F. Improved conditions and new data for identification of lichen products by standardized thin-layer chromatographic method. J. Chromatogr. A 1972, 72, 113–125. [Google Scholar] [CrossRef]
  42. Orange, A.; James, P.W.; White, F.J. Microchemical Methods for the Identification of Lichens; British Lichen Society: London, UK, 2001; p. 101. [Google Scholar]
  43. Dupuis, J.R.; Roe, A.D.; Sperling, F.A.H. Multi-locus species delimitation in closely related animals and fungi: One marker is not enough. Mol. Ecol. 2012, 21, 4422–4436. [Google Scholar] [CrossRef]
  44. Carstens, B.C.; Pelletier, T.A.; Reid, N.M.; Satler, J.D. How to fail at species delimitation. Mol. Ecol. 2013, 22, 4369–4383. [Google Scholar] [CrossRef]
  45. Avise, J.; Ball, R. Principles of genealogical concordance in species concepts and biological taxonomy. Oxf. Surv. Evol. Biol. 1990, 7, 45–67. [Google Scholar]
  46. Maharachchikumbura, S.S.N.; Chen, Y.; Ariyawansa, H.A.; Hyde, K.D.; Haelewaters, D.; Perera, R.H.; Samarakoon, M.C.; Wanasinghe, D.N.; Bustamante, D.E.; Liu, J.-K.; et al. Integrative approaches for species delimitation in Ascomycota. Fungal Divers. 2021, 109, 155–179. [Google Scholar] [CrossRef]
  47. Grewe, F.; Lagostina, E.; Wu, H.; Printzen, C.; Lumbsch, T. Population genomic analyses of RAD sequences resolves the phylogenetic relationship of the lichen-forming fungal species Usneaantarctica and Usneaaurantiacoatra. MycoKeys 2018, 43, 91–113. [Google Scholar] [CrossRef] [PubMed]
  48. Keuler, R.; Garretson, A.; Saunders, T.; Erickson, R.J.; Andre, N.S.; Grewe, F.; Smith, R.; Lumbsch, H.T.; Huang, J.-P.; Clair, L.L.S.; et al. Genome-scale data reveal the role of hybridization in lichen-forming fungi. Sci. Rep. 2020, 10, 1497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  50. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  51. Seppey, M.; Manni, M.; Zdobnov, E.M. BUSCO: Assessing Genome Assembly and Annotation Completeness. Methods Mol. Biol. 2019, 1962, 227–245. [Google Scholar] [CrossRef]
  52. Goff, S.A.; Vaughn, M.; McKay, S.; Lyons, E.; Stapleton, A.E.; Gessler, D.; Matasci, N.; Wang, L.; Hanlon, M.; Lenards, A.; et al. The iPlant Collaborative: Cyberinfrastructure for Plant Biology. Front. Plant Sci. 2011, 2, 34. [Google Scholar] [CrossRef] [Green Version]
  53. Merchant, N.; Lyons, E.; Goff, S.; Vaughn, M.; Ware, D.; Micklos, D.; Antin, P. The iPlant Collaborative: Cyberinfrastructure for Enabling Data to Discovery for the Life Sciences. PLoS Biol. 2016, 14, e1002342. [Google Scholar] [CrossRef] [Green Version]
  54. Johnson, M.G.; Gardner, E.M.; Liu, Y.; Medina, R.; Goffinet, B.; Shaw, A.J.; Zerega, N.; Wickett, N. HybPiper: Extracting coding sequence and introns for phylogenetics from high-throughput sequencing reads using target enrichment. Appl. Plant Sci. 2016, 4, 1600016. [Google Scholar] [CrossRef] [Green Version]
  55. 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] [Green Version]
  56. Nguyen, L.-T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  57. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [Green Version]
  58. 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]
  59. Bertels, F.; Silander, O.; Pachkov, M.; Rainey, P.B.; van Nimwegen, E. Automated Reconstruction of Whole-Genome Phylogenies from Short-Sequence Reads. Mol. Biol. Evol. 2014, 31, 1077–1088. [Google Scholar] [CrossRef] [PubMed]
  60. Puillandre, N.; Brouillet, S.; Achaz, G. ASAP: Assemble species by automatic partitioning. Mol. Ecol. Resour. 2021, 21, 609–620. [Google Scholar] [CrossRef] [PubMed]
  61. Katoh, K.; Toh, H. Recent developments in the MAFFT multiple sequence alignment program. Brief. Bioinform. 2008, 9, 286–298. [Google Scholar] [CrossRef] [Green Version]
  62. Zhang, J.; Kapli, P.; Pavlidis, P.; Stamatakis, A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics 2013, 29, 2869–2876. [Google Scholar] [CrossRef] [Green Version]
  63. Talavera, G.; Castresana, J. Improvement of Phylogenies after Removing Divergent and Ambiguously Aligned Blocks from Protein Sequence Alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef] [Green Version]
  64. Hoang, D.T.; Chernomor, O.; Von Haeseler, A.; Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the Ultrafast Bootstrap Approximation. Mol. Biol. Evol. 2018, 35, 518–522. [Google Scholar] [CrossRef]
  65. Tonini, J.; Moore, A.; Stern, D.; Shcheglovitova, M.; Ortí, G. Concatenation and Species Tree Methods Exhibit Statistically Indistinguishable Accuracy under a Range of Simulated conditions. PLoS Curr. 2015, 7. [Google Scholar] [CrossRef]
  66. Rabiee, M.; Mirarab, S. SODA: Multi-locus species delimitation using quartet frequencies. Bioinformatics 2021, 36, 5623–5631. [Google Scholar] [CrossRef] [PubMed]
  67. Kroken, S.; Taylor, J.W. A gene genealogical approach to recognize phylogenetic species boundaries in the lichenized fungus Letharia. Mycologia 2001, 93, 38–53. [Google Scholar] [CrossRef]
  68. Del-Prado, R.; Buaruang, K.; Lumbsch, H.T.; Crespo, A.; Divakar, P.K. DNA sequence-based identification and barcoding of a morphologically highly plastic lichen forming fungal genus (Parmotrema, Parmeliaceae) from the tropics. Bryologist 2019, 122, 281–291. [Google Scholar] [CrossRef]
  69. Leavitt, S.; Fernández-Mendoza, F.; Pérez-Ortega, S.; Sohrabi, M.; Divakar, P.K.; Lumbsch, T.; Clair, L.S. DNA barcode identification of lichen-forming fungal species in the Rhizoplaca melanophthalma species-complex (Lecanorales, Lecanoraceae), including five new species. MycoKeys 2013, 7, 1–22. [Google Scholar] [CrossRef] [Green Version]
  70. Divakar, P.K.; Leavitt, S.D.; Molina, M.C.; Del-Prado, R.; Lumbsch, H.T.; Crespo, A. A DNA barcoding approach for identification of hidden diversity in Parmeliaceae (Ascomycota): Parmelia sensu strictoas a case study. Bot. J. Linn. Soc. 2015, 180, 21–29. [Google Scholar] [CrossRef] [Green Version]
  71. Leavitt, S.D.; Kraichak, E.; Vondrak, J.; Nelsen, M.; Sohrabi, M.; Perez-Ortega, S.; Clair, L.L.S.; Lumbsch, T. Cryptic diversity and symbiont interactions in rock-posy lichens. Mol. Phylogenet. Evol. 2016, 99, 261–274. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Leavitt, S.D.; Johnson, L.A.; Goward, T.; Clair, L.L.S. Species delimitation in taxonomically difficult lichen-forming fungi: An example from morphologically and chemically diverse Xanthoparmelia (Parmeliaceae) in North America. Mol. Phylogenet. Evol. 2011, 60, 317–332. [Google Scholar] [CrossRef]
  73. Stanton, D.W.G.; Frandsen, P.; Waples, R.K.; Heller, R.; Russo, I.-R.M.; Orozco-Terwengel, P.A.; Pedersen, C.-E.T.; Siegismund, H.R.; Bruford, M.W. More grist for the mill? Species delimitation in the genomic era and its implications for conservation. Conserv. Genet. 2019, 20, 101–113. [Google Scholar] [CrossRef] [Green Version]
  74. Bradshaw, M.; Grewe, F.; Thomas, A.; Harrison, C.H.; Lindgren, H.; Muggia, L.; Clair, L.L.S.; Lumbsch, H.T.; Leavitt, S.D. Characterizing the ribosomal tandem repeat and its utility as a DNA barcode in lichen-forming fungi. BMC Evol. Biol. 2020, 20, 2. [Google Scholar] [CrossRef] [Green Version]
  75. Fu, C.; Mo, Z.; Yang, J.; Cai, J.; Ye, L.; Zou, J.; Qin, H.; Zheng, W.; Hollingsworth, P.M.; Li, D.; et al. Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron. Mol. Ecol. Resour. 2022, 22, 404–414. [Google Scholar] [CrossRef]
  76. Leavitt, S.D.; Grewe, F.; Widhelm, T.; Muggia, L.; Wray, B.; Lumbsch, H.T. Resolving evolutionary relationships in lichen-forming fungi using diverse phylogenomic datasets and analytical approaches. Sci. Rep. 2016, 6, 22262. [Google Scholar] [CrossRef] [PubMed]
  77. Widhelm, T.J.; Grewe, F.; Goffinet, B.; Wedin, M.; Goward, T.; Coca, L.F.; Distefano, I.; Košuthová, A.; Lumbsch, H.T. Phylogenomic reconstruction addressing the Peltigeralean backbone (Lecanoromycetes, Ascomycota). Fungal Divers. 2021, 110, 59–73. [Google Scholar] [CrossRef]
  78. Grewe, F.; Ametrano, C.; Widhelm, T.J.; Leavitt, S.; Distefano, I.; Polyiam, W.; Pizarro, D.; Wedin, M.; Crespo, A.; Divakar, P.K.; et al. Using target enrichment sequencing to study the higher-level phylogeny of the largest lichen-forming fungi family: Parmeliaceae (Ascomycota). IMA Fungus 2020, 11, 27. [Google Scholar] [CrossRef] [PubMed]
  79. Leaché, A.D.; Fujita, M.K.; Minin, V.N.; Bouckaert, R.R. Species Delimitation using Genome-Wide SNP Data. Syst. Biol. 2014, 63, 534–542. [Google Scholar] [CrossRef] [Green Version]
  80. Sukumaran, J.; Holder, M.T.; Knowles, L.L. Incorporating the speciation process into species delimitation. PLoS Comput. Biol. 2021, 17, e1008924. [Google Scholar] [CrossRef]
  81. Fujita, M.K.; Leaché, A.D.; Burbrink, F.T.; McGuire, J.A.; Moritz, C. Coalescent-based species delimitation in an integrative taxonomy. Trends Ecol. Evol. 2012, 27, 480–488. [Google Scholar] [CrossRef]
  82. Leavitt, S.D.; Westberg, M.; Nelsen, M.; Elix, J.A.; Timdal, E.; Sohrabi, M.; Clair, L.L.S.; Williams, L.; Wedin, M.; Lumbsch, H.T. Multiple, Distinct Intercontinental Lineages but Isolation of Australian Populations in a Cosmopolitan Lichen-Forming Fungal Taxon, Psora decipiens (Psoraceae, Ascomycota). Front. Microbiol. 2018, 9, 283. [Google Scholar] [CrossRef] [Green Version]
  83. Lagostina, E.; Andreev, M.; Grande, F.D.; Grewe, F.; Lorenz, A.; Lumbsch, H.T.; Rozzi, R.; Ruprecht, U.; Sancho, L.G.; Søchting, U.; et al. Effects of dispersal strategy and migration history on genetic diversity and population structure of Antarctic lichens. J. Biogeogr. 2021, 48, 1635–1653. [Google Scholar] [CrossRef]
  84. Divakar, P.K.; Wei, X.; McCune, B.; Cubas, P.; Boluda, C.G.; Leavitt, S.D.; Crespo, A.; Tchabanenko, S.; Lumbsch, H.T.; Crespo, H.A. Parallel Miocene dispersal events explain the cosmopolitan distribution of the Hypogymnioid lichens. J. Biogeogr. 2019, 46, 945–955. [Google Scholar] [CrossRef]
  85. Garrido-Benavent, I.; Pérez-Ortega, S. Past, present, and future research in bipolar lichen-forming fungi and their photobionts. Am. J. Bot. 2017, 104, 1660–1674. [Google Scholar] [CrossRef] [Green Version]
  86. Boluda, C.; Rico, V.; Divakar, P.; Nadyeina, O.; Myllys, L.; McMullin, R.; Zamora, J.; Scheidegger, C.; Hawksworth, D. Evaluating methodologies for species delimitation: The mismatch between phenotypes and genotypes in lichenized fungi (Bryoria sect. Implexae, Parmeliaceae). Pers. Mol. Phylogeny Evol. Fungi 2019, 42, 75–100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Fedosov, A.; Achaz, G.; Gontchar, A.; Puillandre, N. MOLD, a novel software to compile accurate and reliable DNA diagnoses for taxonomic descriptions. Mol. Ecol. Resour. 2022, in press. [Google Scholar] [CrossRef]
Jof 08 00490 g001 550
Figure 1. Morphological diversity within the Lecanora polytropa group. (A) L. intricata and L. polytropa (McMullin 13834 CANL); (B) L. polytropa (McMullin 17629 CANL); (C) L. polytropa (McMullin 17695 CANL); (D) L. polytropa (McMullin 17811 CANL); (E) L. polytropa (McMullin 13266 CANL); (F) L. polytropa (McMullin 22536 CANL); (G) L. polytropa (McMullin 8820CANL); (H) L. polytropa (McMullin 22539 CANL); (I) L. polytropa (McMullin 21121 CANL); (J) L. polytropa (Leavitt SL18256 BRY-C); (K) L. polytropa (Leavitt SL18280 BRY-C); (L) L. polytropa (Leavitt SL18356 BRY-C); (M) L. polytropa (Leavitt SL18454 BRY-C); (N) L. polytropa (Leavitt SL18455 BRY-C); and (O) L. polytropa (Leavitt SL18653 BRY-C). Scale bar = 1 mm.
Figure 1. Morphological diversity within the Lecanora polytropa group. (A) L. intricata and L. polytropa (McMullin 13834 CANL); (B) L. polytropa (McMullin 17629 CANL); (C) L. polytropa (McMullin 17695 CANL); (D) L. polytropa (McMullin 17811 CANL); (E) L. polytropa (McMullin 13266 CANL); (F) L. polytropa (McMullin 22536 CANL); (G) L. polytropa (McMullin 8820CANL); (H) L. polytropa (McMullin 22539 CANL); (I) L. polytropa (McMullin 21121 CANL); (J) L. polytropa (Leavitt SL18256 BRY-C); (K) L. polytropa (Leavitt SL18280 BRY-C); (L) L. polytropa (Leavitt SL18356 BRY-C); (M) L. polytropa (Leavitt SL18454 BRY-C); (N) L. polytropa (Leavitt SL18455 BRY-C); and (O) L. polytropa (Leavitt SL18653 BRY-C). Scale bar = 1 mm.
Jof 08 00490 g001
Jof 08 00490 g002 550
Figure 2. Geographic distribution of species within the Lecanora polytropa group based on records available from GBIF (https://www.gbif.org; accessed on 3 February 2022). Sampling sites are indicated with a ‘star’. The locations of specimens representing L. fuscobrunnea (Antarctic endemic, downloaded from GenBank) are not shown. The map was procuded using SimpleMappr (https://www.simplemappr.net; accessed on 3 February 2022).
Figure 2. Geographic distribution of species within the Lecanora polytropa group based on records available from GBIF (https://www.gbif.org; accessed on 3 February 2022). Sampling sites are indicated with a ‘star’. The locations of specimens representing L. fuscobrunnea (Antarctic endemic, downloaded from GenBank) are not shown. The map was procuded using SimpleMappr (https://www.simplemappr.net; accessed on 3 February 2022).
Jof 08 00490 g002
Jof 08 00490 g003 550
Figure 3. (previous page). Midpoint-rooted, five-marker (ITS, nuLSU, RPB1, RPB2, mtSSU) maximum likelihood topology of the L. polytropa group. Clades are labeled as the candidate species inferred from ITS data using ASAP, with the number of specimens representing each candidate species in parentheses. Thickened black branches indicate bootstrap (BS) values between 95–100%, dashed lines indicate BS values below 50%, and otherwise, bootstrap values are indicated at nodes. Clades highlighted in pink represent other taxa recovered in the L. polytropa clade—‘ch’, L. chlorophaeodes; ‘co’, L. concolor; ‘di’, L. dispersoareolata; ‘fu’, L. fuscobrunnea; ‘in’, L. intricata; ‘sb’, L. cf. subcintula’; ‘sm’, L. somervelii; ‘so’, L. solaris; ‘su’, L. subintricata—tips labels coincide with names in Supplementary File S1 in the ‘Candidate ASAP species partition code’ column. The geographic distribution of each candidate species is shown following the tip label: Antarctica (no distinction for different regions); Arctic (‘N’ = Nunavik; ‘S’ = Svalbard; ‘M’ = multiple locations); China (‘Q’= Qinghai; ‘S’= Sichuan; ‘T’ = Tibet; ‘M’ = multiple locations); Russia (no distinction for different regions); Japan (no distinction for different regions); Taiwan (no distinction for different regions); ‘Western NA’ is western North America (‘B’ = Beartooth Plateau, MT, USA; ‘E’ = Escalante Region, UT, USA; ‘L’ = La Sal Mountains, UT, USA; ‘P’ = Transverse Range, CA, USA; ‘S’= Saguache Range, CO, USA; ‘U’ = Uinta Mountains, UT, USA; ‘M’ = multiple locations); ‘Eastern NA’ is eastern North America (‘M’ = multiple locations); Bolivia, South America (no distinction for different regions). Carbonea specimens were recovered as sisters to the remaining samples and not shown—the complete tree is available as Supplementary File S3.
Figure 3. (previous page). Midpoint-rooted, five-marker (ITS, nuLSU, RPB1, RPB2, mtSSU) maximum likelihood topology of the L. polytropa group. Clades are labeled as the candidate species inferred from ITS data using ASAP, with the number of specimens representing each candidate species in parentheses. Thickened black branches indicate bootstrap (BS) values between 95–100%, dashed lines indicate BS values below 50%, and otherwise, bootstrap values are indicated at nodes. Clades highlighted in pink represent other taxa recovered in the L. polytropa clade—‘ch’, L. chlorophaeodes; ‘co’, L. concolor; ‘di’, L. dispersoareolata; ‘fu’, L. fuscobrunnea; ‘in’, L. intricata; ‘sb’, L. cf. subcintula’; ‘sm’, L. somervelii; ‘so’, L. solaris; ‘su’, L. subintricata—tips labels coincide with names in Supplementary File S1 in the ‘Candidate ASAP species partition code’ column. The geographic distribution of each candidate species is shown following the tip label: Antarctica (no distinction for different regions); Arctic (‘N’ = Nunavik; ‘S’ = Svalbard; ‘M’ = multiple locations); China (‘Q’= Qinghai; ‘S’= Sichuan; ‘T’ = Tibet; ‘M’ = multiple locations); Russia (no distinction for different regions); Japan (no distinction for different regions); Taiwan (no distinction for different regions); ‘Western NA’ is western North America (‘B’ = Beartooth Plateau, MT, USA; ‘E’ = Escalante Region, UT, USA; ‘L’ = La Sal Mountains, UT, USA; ‘P’ = Transverse Range, CA, USA; ‘S’= Saguache Range, CO, USA; ‘U’ = Uinta Mountains, UT, USA; ‘M’ = multiple locations); ‘Eastern NA’ is eastern North America (‘M’ = multiple locations); Bolivia, South America (no distinction for different regions). Carbonea specimens were recovered as sisters to the remaining samples and not shown—the complete tree is available as Supplementary File S3.
Jof 08 00490 g003
Jof 08 00490 g004 550
Figure 4. Maximum likelihood (ML) topologies inferred from nuclear (A) and mitochondrial (B) phylogenomic datasets. (A). A ML topology from the concatenated alignments of 1209 BUSCO single copy nuclear markers spanning 2.28 MBps (Supplementary File S9). Colored clades and clade names are linked to the candidate species inferred from the standard DNA barcoding marker for fungi (ITS), and candidate species inferred from the 1209 BUSCO gene trees using SODA are shown with dashed boxes. Only bootstrap values below 100% are shown. (B). A ML topology from the 66.5 KBp mitochondrial alignment generated using RealPhy (Supplementary File S10). Only bootstrap values below 100% are shown.
Figure 4. Maximum likelihood (ML) topologies inferred from nuclear (A) and mitochondrial (B) phylogenomic datasets. (A). A ML topology from the concatenated alignments of 1209 BUSCO single copy nuclear markers spanning 2.28 MBps (Supplementary File S9). Colored clades and clade names are linked to the candidate species inferred from the standard DNA barcoding marker for fungi (ITS), and candidate species inferred from the 1209 BUSCO gene trees using SODA are shown with dashed boxes. Only bootstrap values below 100% are shown. (B). A ML topology from the 66.5 KBp mitochondrial alignment generated using RealPhy (Supplementary File S10). Only bootstrap values below 100% are shown.
Jof 08 00490 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, Y.; Clancy, J.; Jensen, J.; McMullin, R.T.; Wang, L.; Leavitt, S.D. Providing Scale to a Known Taxonomic Unknown—At Least a 70-Fold Increase in Species Diversity in a Cosmopolitan Nominal Taxon of Lichen-Forming Fungi. J. Fungi 2022, 8, 490. https://doi.org/10.3390/jof8050490

AMA Style

Zhang Y, Clancy J, Jensen J, McMullin RT, Wang L, Leavitt SD. Providing Scale to a Known Taxonomic Unknown—At Least a 70-Fold Increase in Species Diversity in a Cosmopolitan Nominal Taxon of Lichen-Forming Fungi. Journal of Fungi. 2022; 8(5):490. https://doi.org/10.3390/jof8050490

Chicago/Turabian Style

Zhang, Yanyun, Jeffrey Clancy, Jacob Jensen, Richard Troy McMullin, Lisong Wang, and Steven D. Leavitt. 2022. "Providing Scale to a Known Taxonomic Unknown—At Least a 70-Fold Increase in Species Diversity in a Cosmopolitan Nominal Taxon of Lichen-Forming Fungi" Journal of Fungi 8, no. 5: 490. https://doi.org/10.3390/jof8050490

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