Multiple ITS Haplotypes in the Genome of the Lichenized Basidiomycete Cora inversa (Hygrophoraceae): Fact or Artifact?

Elucidation of intragenomic variation of ribosomal DNA sequences in the enigmatic fungal genus Ceraceosorus, including a newly described species Ceraceosorus americanus

Multicopy nuclear ribosomal DNA (rDNA) genes have been used as markers for fungal identification for three decades. The rDNA sequences in a genome are thought to be homogeneous due to concerted evolution. However, intragenomic variation of rDNA sequences has recently been observed in many fungi, which may make fungal identification and species abundance estimation based on these loci problematic. Ceraceosorus is an enigmatic genus in the smut lineage Ustilaginomycotina for which very limited distribution data exist. Our previous research demonstrated intragenomic variation in the internal transcribed spacer (ITS1-5.8S-ITS2) region of two Ceraceosorus species. In this study, we described the fourth known species of Ceraceosorus, C. americanus, isolated from an asymptomatic rosemary leaf collected in Louisiana, USA. This is the first report of this genus in the Americas. We then selected all four known Ceraceosorus species, plus exemplar smut fungi representing all major lineages of subphylum Ustilaginomycotina, to examine sequence heterogeneity in three regions of the rDNA repeat (partial 18S, ITS, and partial 28S regions). Three methods were used: PCR-cloning-Sanger sequencing, targeted amplicon high-throughput sequencing, and whole-genome shotgun high-throughput sequencing. Our results show that Ceraceosorus is the only sampled fungal genus in Ustilaginomycotina with significant intragenomic variation at the ITS, with up to 25 nucleotide variant sites in the ITS1-5.8S-ITS2 region and 2.6% divergence among analyzed ITS haplotypes. We found many conflicting patterns across the three detection methods, with up to 27 conflicting variant sites recorded from a single individual. At least 40% of the conflicting patterns are possibly due to PCR-cloning-sequencing errors, as the corresponding variant sites were not observed in the other detection methods. Based on our data and the literature, we evaluated the characteristics and advantages/disadvantages of each detection method. Finally, a model for how intragenomic variation in the rDNA copies within a genome may arise is presented.

Novel insights in cryptic diversity of snow and glacier ice algae communities combining 18S rRNA gene and ITS2 amplicon sequencing

Melting snow and glacier surfaces host microalgal blooms in polar and mountainous regions. The aim of this study was to determine the dominant taxa at the species level in the European Arctic and the Alps. A standardized protocol for amplicon metabarcoding using the 18S rRNA gene and ITS2 markers was developed. This is important because previous biodiversity studies have been hampered by the dominance of closely related algal taxa in snow and ice. Due to the limited resolution of partial 18S rRNA Illumina sequences, the hypervariable ITS2 region was used to further discriminate between the genotypes. Our results show that red snow was caused by the cosmopolitan Sanguina nivaloides (Chlamydomonadales, Chlorophyta) and two as of yet undescribed Sanguina species. Arctic orange snow was dominated by S. aurantia, which was not found in the Alps. On glaciers, at least three Ancylonema species (Zygnematales, Streptophyta) dominated. Golden-brown blooms consisted of Hydrurus spp. (Hydrurales, Stramenophiles) and these were mainly an Arctic phenomenon. For chrysophytes, only the 18S rRNA gene but not ITS2 sequences were amplified, showcasing how delicate the selection of eukaryotic 'universal' primers for community studies is and that primer specificity will affect diversity results dramatically. We propose our approach as a 'best practice'.

Extensive intragenomic variation in the internal transcribed spacer region of fungi

Fungi are among the most biodiverse organisms in the world. Accurate species identification is imperative for studies on fungal ecology and evolution. The internal transcribed spacer (ITS) rDNA region has been widely accepted as the universal barcode for fungi. However, several recent studies have uncovered intragenomic sequence variation within the ITS in multiple fungal species. Here, we mined the genome of 2414 fungal species to determine the prevalence of intragenomic variation and found that the genomes of 641 species, about one-quarter of the 2414 species examined, contained multiple ITS copies. Of those 641 species, 419 (∼65%) contained variation among copies revealing that intragenomic variation is common in fungi. We proceeded to show how these copies could result in the erroneous description of hundreds of fungal species and skew studies evaluating environmental DNA (eDNA) especially when making diversity estimates. Additionally, many genomes were found to be contaminated, especially those of unculturable fungi.

Survey of Lichenized Fungi DNA Barcodes on King George Island (Antarctica): An Aid to Species Discovery

DNA barcoding is a powerful method for the identification of lichenized fungi groups for which the diversity is already well-represented in nucleotide databases, and an accurate, robust taxonomy has been established. However, the effectiveness of DNA barcoding for identification is expected to be limited for understudied taxa or regions. One such region is Antarctica, where, despite the importance of lichens and lichenized fungi identification, their genetic diversity is far from characterized. The aim of this exploratory study was to survey the lichenized fungi diversity of King George Island using a fungal barcode marker as an initial identification tool. Samples were collected unrestricted to specific taxa in coastal areas near Admiralty Bay. Most samples were identified using the barcode marker and verified up to the species or genus level with a high degree of similarity. A posterior morphological evaluation focused on samples with novel barcodes allowed for the identification of unknown Austrolecia, Buellia, and Lecidea s.l. species. These results contribute to better represent the lichenized fungi diversity in understudied regions such as Antarctica by increasing the richness of the nucleotide databases. Furthermore, the approach used in this study is valuable for exploratory surveys in understudied regions to guide taxonomic efforts towards species recognition and discovery.

A long-read amplicon approach to scaling up the metabarcoding of lichen herbarium specimens

Reference sequence databases are critical to the accurate detection and identification of fungi in the environment. As repositories of large numbers of well-curated specimens, herbaria and fungal culture collections have the material resources to generate sequence data for large number of taxa, and could therefore allow filling taxonomic gaps often present in reference sequence databases. Financial resources to do that are however often lacking, so that recent efforts have focused on decreasing sequencing cost by increasing the number of multiplexed samples per sequencing run while maintaining high sequence quality. Following a previous study that aimed at decreasing sequencing cost for lichen specimens by generating fungal ITS barcodes for 96 specimens using PacBio amplicon sequencing, we present a method that further decreases lichen specimen metabarcoding costs. A total of 384 mixed DNA extracts obtained from lichen herbarium specimens, mostly from the four genera Buellia, Catillaria, Endocarpon and Parmotrema, were used to generate new fungal ITS sequences using a Sequel I sequencing platform and the PacBio M13 barcoded primers. The average success rate across all taxa was high (86.5%), with particularly high rates for the crustose saxicolous taxa (Buellia, Catillaria and others; 93.3%) and the terricolous squamulose taxa (Endocarpon and others; 96.5%). On the other hand, the success rate for the foliose genus Parmotrema was lower (60.4%). With this taxon sampling, greater specimen age did not appear to impact sequencing success. In fact, the 1966–1980 collection date category showed the highest success rate (97.3%). Compared to the previous study, the abundance-based sequence denoising method showed some limitations, but the cost of generating ITS barcodes was further decreased thanks to the higher multiplexing level. In addition to contributing new ITS barcodes for specimens of four interesting lichen genera, this study further highlights the potential and challenges of using new sequencing technologies on collection specimens to generate DNA sequences for reference databases.

Extensive photobiont sharing in a rapidly radiating cyanolichen clade

Recent studies have uncovered remarkable diversity in Dictyonema s.lat. basidiolichens, here recognized as subtribe Dictyonemateae. This group includes five genera and 148 species, but hundreds more await description. The photobionts of these lichens belong to Rhizonema, a recently resurrected cyanobacterial genus known by a single species. To further investigate photobiont diversity within Dictyonemateae, we generated 765 new cyanobacterial sequences from 635 specimens collected from 18 countries. The ITS barcoding locus supported the recognition of 200 mycobiont (fungal) species among these samples, but the photobiont diversity was comparatively low. Our analyses revealed three main divisions of Rhizonema, with two repeatedly recovered as monophyletic (proposed as new species), and the third mostly paraphyletic. The paraphyletic lineage corresponds to R. interruptum and partnered with mycobionts from all five genera in Dictyonemateae. There was no evidence of photobiont-mycobiont co-speciation, but one of the monophyletic lineages of Rhizonema appears to partner predominantly with one of the two major clades of Cora (mycobiont) with samples collected largely from the northern Andes. Molecular clock estimations indicate the Rhizonema species are much older than the fungal species in the Dictyonemateae, suggesting that these basidiolichens obtained their photobionts from older ascolichen lineages and the photobiont variation in extant lineages of Dictyonemateae is the result of multiple photobiont switches. These results support the hypothesis of lichens representing "fungal farmers," in which diverse mycobiont lineages associate with a substantially lower diversity of photobionts by sharing those photobionts best suited for the lichen symbiosis among multiple and often unrelated mycobiont lineages.

Comparing Morphological and Molecular Estimates of Species Diversity in the Freshwater Genus Synura (Stramenopiles): A Model for Understanding Diversity of Eukaryotic Microorganisms

We performed a comparison of molecular and morphological diversity in a freshwater colonial genus Synura (Chrysophyceae, Stramenopiles), using the island of Newfoundland (Canada) as a case study. We examined the morphological species diversity in collections from 79 localities, and compared these findings to diversity based on molecular characters for 150 strains isolated from the same sites. Of 27 species or species-level lineages identified, only one third was recorded by both molecular and morphological techniques, showing both approaches are complementary in estimating species diversity within this genus. Eight taxa, each representing young evolutionary lineages, were recovered only by sequencing of isolated colonies, whereas ten species were recovered only microscopically. Our complex investigation, involving both morphological and molecular examinations, indicates that our knowledge of Synura diversity is still poor, limited only to a few well-studied areas. We revealed considerable cryptic diversity within the core S. petersenii and S. leptorrhabda lineages. We further resolved the phylogenetic position of two previously described taxa, S. kristiansenii and S. petersenii f. praefracta, propose species-level status for S. petersenii f. praefracta, and describe three new species, S. vinlandica, S. fluviatilis, and S. cornuta. Our findings add to the growing body of literature detailing distribution patterns observed in the genus, ranging from cosmopolitan species, to highly restricted taxa, to species such as S. hibernica found along coastal regions on multiple continents. Finally, our study illustrates the usefulness of combining detailed morphological information with gene sequence data to examine species diversity within chrysophyte algae.

Evidence of Intra-individual SSU Polymorphisms in Dark-spored Myxomycetes (Amoebozoa)

The nuclear small subunit rRNA gene (SSU or 18S) is a marker frequently used in phylogenetic and barcoding studies in Amoebozoa, including Myxomycetes. Despite its common usage and the confirmed existence of divergent copies of ribosomal genes in other protists, the potential presence of intra-individual SSU variability in Myxomycetes has never been studied before. Here we investigated the pattern of nucleotide polymorphism in the 5' end fragment of SSU by cloning and sequencing a total of 238 variants from eight specimens, each representing a species of the dark-spored orders Stemonitidales and Physarales. After excluding singletons, a relatively low SSU intra-individual variability was found but our data indicate that this might be a widely distributed phenomenon in Myxomycetes as all samples analyzed possessed various ribotypes. To determine if the occurrence of multiple SSU variants within a single specimen has a negative effect on the circumscription of species boundaries, we conducted phylogenetic analyses that revealed that clone variation may be detrimental for inferring phylogenetic relationships among some of the specimens analyzed. Despite that intra-individual variability should be assessed in additional taxa, our results indicate that special care should be taken for species identification when working with closely related species.

PacBio amplicon sequencing for metabarcoding of mixed DNA samples from lichen herbarium specimens

The detection and identification of species of fungi in the environment using molecular methods heavily depends on reliable reference sequence databases. However, these databases are largely incomplete in terms of taxon coverage, and a significant effort is required from herbaria and living fungal collections for the mass-barcoding of well-identified and well-curated fungal specimens or strains. Here, a PacBio amplicon sequencing approach is applied to recent lichen herbarium specimens for the sequencing of the fungal ITS barcode, allowing a higher throughput sample processing than Sanger sequencing, which often required the use of cloning. Out of 96 multiplexed samples, a full-length ITS sequence of the target lichenised fungal species was recovered for 85 specimens. In addition, sequences obtained for co-amplified fungi gave an interesting insight into the diversity of endolichenic fungi. Challenges encountered at both the laboratory and bioinformatic stages are discussed, and cost and quality are compared with Sanger sequencing. With increasing data output and reducing sequencing cost, PacBio amplicon sequencing is seen as a promising approach for the generation of reference sequences for lichenised fungi as well as the characterisation of lichen-associated fungal communities.

First insights into the living groundwater mycobiome of the terrestrial biogeosphere

Although fungi play important roles in biogeochemical cycling in aquatic ecosystems and have received a great deal of attention, much remains unknown about the living fractions of fungal communities in aquifers of the terrestrial subsurface in terms of diversity, community dynamics, functional roles, the impact of environmental factors and presence of fungal pathogens. Here we address this gap in knowledge by using RNA-based high throughput pair-end illumina sequencing analysis of fungal internal transcribed spacer (ITS) gene markers, to target the living fractions of groundwater fungal communities from fractured alternating carbonate-/siliciclastic-rock aquifers of the Hainich Critical Zone Exploratory. The probed levels of the hillslope multi-storey aquifer system differ primarily in their oxygen and nitrogen content due to their different connections to the surface. We discovered highly diverse living fungal communities (384 Operational Taxonomic Units, OTUs) with different taxonomic affiliations and ecological functions. The observed fungal communities primarily belonged to three phyla: Ascomycota, Basidiomycota and Chytridiomycota. Perceived dynamics in the composition of living fungal communities were significantly shaped by the concentration of ammonium in the moderately agriculturally impacted aquifer system. Apart from fungal saprotrophs, we also detected living plant and animal pathogens for the first time in this aquifer system. This work also demonstrates that the RNA-based high throughput pair-end illumina sequencing method can be used in future for water quality monitoring in terms of living fungal load and subsequent risk assessments. In general, this study contributes towards the growing knowledge of aquatic fungi in terrestrial subsurface biogeosphere.

TBC1D21 Potentially Interacts with and Regulates Rap1 during Murine Spermatogenesis

Few papers have focused on small guanosine triphosphate (GTP)-binding proteins and their regulation during spermatogenesis. TBC1D21 genes (also known as male germ cell RAB GTPase-activating protein MGCRABGAP) are related to sterility, as determined through cDNA microarray testing of human testicular tissues exhibiting spermatogenic defects. TBC1D21 is a protein specifically expressed in the testes that exhibits specific localizations of elongating and elongated spermatids during mammalian spermiogenesis. Furthermore, through co-immunoprecipitation (co-IP) and nano liquid chromatography–tandem mass spectrometry (nano LC–MS/MS), Rap1 has been recognized as a potential TBC1D21 interactor. This study determined the possible roles of Rap1 and TBC1D21 during mammalian spermiogenesis. First, the binding ability between Rap1 and TBC1D21 was verified using co-IP. Second, the stronger signals of Rap1 expressed in elongating and elongated murine spermatids extracted from testicular sections, namely spermatogonia, spermatocytes, and round spermatids, were compared. Third, Rap1 and TBC1D21 exhibited similar localizations at postacrosomal regions of spermatids and at the midpieces of mature sperms, through isolated male germ cells. Fourth, the results of an activating Rap1 pull-down assay indicated that TBC1D21 overexpression inactivates Rap1 activity in cell models. In conclusion, TBC1D21 may interact with and potentially regulate Rap1 during murine spermatogenesis.

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.

rRNA Pseudogenes in Filamentous Ascomycetes as Revealed by Genome Data

The nuclear ribosomal DNA (rDNA) is considered as a paradigm of concerted evolution. Components of the rDNA tandem repeats (45S) are widely used in phylogenetic studies of different organisms and the internal transcribed spacer (ITS) region was recently selected as a fungal DNA bar code. However, rRNA pseudogenes, as one kind of escape from concerted evolution, were reported in a wide range of organisms, especially in plants and animals. Moreover, large numbers of 5S rRNA pseudogenes were identified in several filamentous ascomycetes. To study whether rDNA evolves in a strict concerted manner and test whether rRNA pseudogenes exist in more species of ascomycetes, intragenomic rDNA polymorphisms were analyzed using whole genome sequences. Divergent rDNA paralogs were found to coexist within a single genome in seven filamentous ascomycetes examined. A great number of paralogs were identified as pseudogenes according to the mutation and secondary structure analyses. Phylogenetic analyses of the three rRNA coding regions of the 45S rDNA repeats, i.e., 18S, 5.8S, and 28S, revealed an interspecies clustering pattern of those different rDNA paralogs. The identified rRNA pseudogenic sequences were validated using specific primers designed. Mutation analyses revealed that the repeat-induced point (RIP) mutation was probably responsible for the formation of those rRNA pseudogenes.

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.

Characterization of Unexplored Deadwood Mycobiome in Highly Diverse Subtropical Forests Using Culture-independent Molecular Technique

The deadwood mycobiome, also known as wood-inhabiting fungi (WIF), are among the key players in wood decomposition, having a large impact on nutrient cycling in forest soils. However, our knowledge of WIF richness and distribution patterns in different forest biomes is limited. Here, we used pyrotag sequencing of the fungal internal transcribed spacer (ITS2) region to characterize the deadwood mycobiome of two tree species with greatly different wood characteristics (Schima superba and Pinus massoniana) in a Chinese subtropical forest ecosystem. Specifically, we tested (i) the effects of tree species and wood quality properties on WIF OTU richness and community composition; (ii) the role of biotic and abiotic factors in shaping the WIF communities; and (iii) the relationship between WIF OTU richness, community composition and decomposition rates. Due to different wood chemical properties, we hypothesized that the WIF communities derived from the two tree species would be correlated differently with biotic and abiotic factors. Our results show that deadwood in subtropical forests harbors diverse fungal communities comprising six ecological functional groups. We found interesting colonization patterns for this subtropical biome, where Resinicium spp. were highly detected in both broadleaved and coniferous deadwood. In addition, the members of Xylariales were frequently found in Schima. The two deadwood species differed significantly in WIF OTU richness (Pinus > Schima) and community composition (P < 0.001). Variations in WIF community composition of both tree species were significantly explained by wood pH and ecological factors (biotic: deadwood species, basal area and abiotic: soil pH), but the WIF communities derived from each tree species correlated differently with abiotic factors. Interestingly, we found that deadwood decomposition rate significantly correlated with WIF communities and negatively correlated with WIF OTU richness. We conclude that the pattern of WIF OTU richness and community composition are controlled by multiple interacting biotic and abiotic factors. Overall, our study provides an in-depth picture of the deadwood mycobiome in this subtropical forest. Furthermore, by comparing our results to results from temperate and boreal forests we contribute to a better understanding of patterns of WIF communities across different biomes and geographic locations.

Barcoding lichen-forming fungi using 454 pyrosequencing is challenged by artifactual and biological sequence variation

Although lichens (lichen-forming fungi) play an important role in the ecological integrity of many vulnerable landscapes, only a minority of lichen-forming fungi have been barcoded out of the currently accepted ∼18 000 species. Regular Sanger sequencing can be problematic when analyzing lichens since saprophytic, endophytic, and parasitic fungi live intimately admixed, resulting in low-quality sequencing reads. Here, high-throughput, long-read 454 pyrosequencing in a GS FLX+ System was tested to barcode the fungal partner of 100 epiphytic lichen species from Switzerland using fungal-specific primers when amplifying the full internal transcribed spacer region (ITS). The present study shows the potential of DNA barcoding using pyrosequencing, in that the expected lichen fungus was successfully sequenced for all samples except one. Alignment solutions such as BLAST were found to be largely adequate for the generated long reads. In addition, the NCBI nucleotide database—currently the most complete database for lichen-forming fungi—can be used as a reference database when identifying common species, since the majority of analyzed lichens were identified correctly to the species or at least to the genus level. However, several issues were encountered, including a high sequencing error rate, multiple ITS versions in a genome (incomplete concerted evolution), and in some samples the presence of mixed lichen-forming fungi (possible lichen chimeras).

Contrasting evolutionary patterns of 28S and ITS rRNA genes reveal high intragenomic variation in Cephalenchus (Nematoda): Implications for species delimitation

Concerted evolution is often assumed to be the evolutionary force driving multi-family genes, including those from ribosomal DNA (rDNA) repeat, to complete homogenization within a species, although cases of non-concerted evolution have been also documented. In this study, sequence variation of 28S and ITS ribosomal RNA (rRNA) genes in the genus Cephalenchus is assessed at three different levels, intragenomic, intraspecific, and interspecific. The findings suggest that not all Cephalenchus species undergo concerted evolution. High levels of intraspecific polymorphism, mostly due to intragenomic variation, are found in Cephalenchus sp1 (BRA-01). Secondary structure analyses of both rRNA genes and across different species show a similar substitution pattern, including mostly compensatory (CBC) and semi-compensatory (SBC) base changes, thus suggesting the functionality of these rRNA copies despite the variation found in some species. This view is also supported by low sequence variation in the 5.8S gene in relation to the flanking ITS-1 and ITS-2 as well as by the existence of conserved motifs in the former gene. It is suggested that potential cross-fertilization in some Cephalenchus species, based on inspection of female reproductive system, might contribute to both intragenomic and intraspecific polymorphism of their rRNA genes. These results reinforce the potential implications of intragenomic and intraspecific genetic diversity on species delimitation, especially in biodiversity studies based solely on metagenetic approaches. Knowledge of sequence variation will be crucial for accurate species diversity estimation using molecular methods.

Bioinformatic Amplicon Read Processing Strategies Strongly Affect Eukaryotic Diversity and the Taxonomic Composition of Communities

Amplicon read sequencing has revolutionized the field of microbial diversity studies. The technique has been developed for bacterial assemblages and has undergone rigorous testing with mock communities. However, due to the great complexity of eukaryotes and the numbers of different rDNA copies, analyzing eukaryotic diversity is more demanding than analyzing bacterial or mock communities, so studies are needed that test the methods of analyses on taxonomically diverse natural communities. In this study, we used 20 samples collected from the Baltic Sea ice, slush and under-ice water to investigate three program packages (UPARSE, mothur and QIIME) and 18 different bioinformatic strategies implemented in them. Our aim was to assess the impact of the initial steps of bioinformatic strategies on the results when analyzing natural eukaryotic communities. We found significant differences among the strategies in resulting read length, number of OTUs and estimates of diversity as well as clear differences in the taxonomic composition of communities. The differences arose mainly because of the variable number of chimeric reads that passed the pre-processing steps. Singleton removal and denoising substantially lowered the number of errors. Our study showed that the initial steps of the bioinformatic amplicon read processing strategies require careful consideration before applying them to eukaryotic communities.

A single macrolichen constitutes hundreds of unrecognized species

The number of Fungi is estimated at between 1.5 and 3 million. Lichenized species are thought to make up a comparatively small portion of this figure, with unrecognized species richness hidden among little-studied, tropical microlichens. Recent findings, however, suggest that some macrolichens contain a large number of unrecognized taxa, increasing known species richness by an order of magnitude or more. Here we report the existence of at least 126 species in what until recently was believed to be a single taxon: the basidiolichen fungus Dictyonema glabratum, also known as Cora pavonia. Notably, these species are not cryptic but morphologically distinct. A predictive model suggests an even larger number, with more than 400 species. These results call into question species concepts in presumably well-known macrolichens and demonstrate the need for accurately documenting such species richness, given the importance of these lichens in endangered ecosystems such as paramos and the alarming potential for species losses throughout the tropics.