Airborne DNA reveals predictable spatial and seasonal dynamics of fungi.

Fungi are among the most diverse and ecologically important kingdoms in life. However, the distributional ranges of fungi remain largely unknown as do the ecological mechanisms that shape their distributions1,2. To provide an integrated view of the spatial and seasonal dynamics of fungi, we implemented a globally distributed standardized aerial sampling of fungal spores3. The vast majority of operational taxonomic units were detected within only one climatic zone, and the spatiotemporal patterns of species richness and community composition were mostly explained by annual mean air temperature. Tropical regions hosted the highest fungal diversity except for lichenized, ericoid mycorrhizal and ectomycorrhizal fungi, which reached their peak diversity in temperate regions. The sensitivity in climatic responses was associated with phylogenetic relatedness, suggesting that large-scale distributions of some fungal groups are partially constrained by their ancestral niche. There was a strong phylogenetic signal in seasonal sensitivity, suggesting that some groups of fungi have retained their ancestral trait of sporulating for only a short period. Overall, our results show that the hyperdiverse kingdom of fungi follows globally highly predictable spatial and temporal dynamics, with seasonality in both species richness and community composition increasing with latitude. Our study reports patterns resembling those described for other major groups of organisms, thus making a major contribution to the long-standing debate on whether organisms with a microbial lifestyle follow the global biodiversity paradigms known for macroorganisms4,5.

Experimental evidence that root-associated fungi improve plant growth at high altitude.

Unravelling how species communities change along environmental gradients requires a dual understanding: the direct responses of the species to their abiotic surroundings and the indirect variation of these responses through biotic interactions. Here, we focus on the interactive relationships between plants and their symbiotic root-associated fungi (RAF) along stressful abiotic gradients. We investigate whether variations in RAF community composition along altitudinal gradients influence plant growth at high altitudes, where both plants and fungi face harsher abiotic conditions. We established a translocation experiment between pairs of Bistorta vivipara populations across altitudinal gradients. To separate the impact of shifting fungal communities from the overall influence of changing abiotic conditions, we used a root barrier to prevent new colonization by RAF following translocation. To characterize the RAF communities, we applied DNA barcoding to the root samples. Through the utilization of joint species distribution modelling, we assessed the relationship between changes in plant functional traits resulting from experimental treatments and the corresponding changes in the RAF communities. Our findings indicate that RAF communities influence plant responses to stressful abiotic conditions. Plants translocated from low to high altitudes grew more when they were able to associate with the resident high-altitude RAF compared to those plants that were not allowed to associate with the resident RAF. We conclude that interactions with RAF impact how plants respond to stressful abiotic conditions. Our results provide experimental support that interactions with RAF improve plant stress tolerance to altitudinal stressors such as colder temperatures and less nutrient availability.

PROTAX-GPU: a scalable probabilistic taxonomic classification system for DNA barcodes.

DNA-based identification is vital for classifying biological specimens, yet methods to quantify the uncertainty of sequence-based taxonomic assignments are scarce. Challenges arise from noisy reference databases, including mislabelled entries and missing taxa. PROTAX addresses these issues with a probabilistic approach to taxonomic classification, advancing on methods that rely solely on sequence similarity. It provides calibrated probabilistic assignments to a partially populated taxonomic hierarchy, accounting for taxa that lack references and incorrect taxonomic annotation. While effective on smaller scales, global application of PROTAX necessitates substantially larger reference libraries, a goal previously hindered by computational barriers. We introduce PROTAX-GPU, a scalable algorithm capable of leveraging the global Barcode of Life Data System (>14 million specimens) as a reference database. Using graphics processing units (GPU) to accelerate similarity and nearest-neighbour operations and the JAX library for Python integration, we achieve over a 1000 × speedup compared with the central processing unit (CPU)-based implementation without compromising PROTAX's key benefits. PROTAX-GPU marks a significant stride towards real-time DNA barcoding, enabling quicker and more efficient species identification in environmental assessments. This capability opens up new avenues for real-time monitoring and analysis of biodiversity, advancing our ability to understand and respond to ecological dynamics. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Global Spore Sampling Project: A global, standardized dataset of airborne fungal DNA.

Novel methods for sampling and characterizing biodiversity hold great promise for re-evaluating patterns of life across the planet. The sampling of airborne spores with a cyclone sampler, and the sequencing of their DNA, have been suggested as an efficient and well-calibrated tool for surveying fungal diversity across various environments. Here we present data originating from the Global Spore Sampling Project, comprising 2,768 samples collected during two years at 47 outdoor locations across the world. Each sample represents fungal DNA extracted from 24 m3 of air. We applied a conservative bioinformatics pipeline that filtered out sequences that did not show strong evidence of representing a fungal species. The pipeline yielded 27,954 species-level operational taxonomic units (OTUs). Each OTU is accompanied by a probabilistic taxonomic classification, validated through comparison with expert evaluations. To examine the potential of the data for ecological analyses, we partitioned the variation in species distributions into spatial and seasonal components, showing a strong effect of the annual mean temperature on community composition.

Effects of landscape and distance in automatic audio based bird species identification.

The present work focuses on how the landscape and distance between a bird and an audio recording unit affect automatic species identification. Moreover, it is shown that automatic species identification can be improved by taking into account the effects of landscape and distance. The proposed method uses measurements of impulse responses between the sound source and the recorder. These impulse responses, characterizing the effect of a landscape, can be measured in the real environment, after which they can be convolved with any number of recorded bird sounds to modify an existing set of bird sound recordings. The method is demonstrated using autonomous recording units on an open field and in two different types of forests, varying the distance between the sound source and the recorder. Species identification accuracy improves significantly when the landscape and distance effect is taken into account when building the classification model. The method is demonstrated using bird sounds, but the approach is applicable to other animal and non-animal vocalizations as well.

MycoPins: a metabarcoding-based method to monitor fungal colonization of fine woody debris.

The MycoPins method described here is a rapid and affordable protocol to monitor early colonization events in communities of wood-inhabiting fungi in fine woody debris. It includes easy to implement field sampling techniques and sample processing, followed by data processing, and analysis of the development of early dead wood fungal communities. The method is based on fieldwork from a time series experiment on standard sterilized colonization targets followed by the metabarcoding analysis and automated molecular identification of species. This new monitoring method through its simplicity, moderate costs, and scalability paves a way for a broader and scalable project pipeline. MycoPins establishes a standard routine for research stations or regularly visited field sites for monitoring of fungal colonization of woody substrates. The routine uses widely available consumables and therefore presents a unifying method for monitoring of fungi of this type.

Spatial cumulant models enable spatially informed treatment strategies and analysis of local interactions in cancer systems.

Theoretical and applied cancer studies that use individual-based models (IBMs) have been limited by the lack of a mathematical formulation that enables rigorous analysis of these models. However, spatial cumulant models (SCMs), which have arisen from theoretical ecology, describe population dynamics generated by a specific family of IBMs, namely spatio-temporal point processes (STPPs). SCMs are spatially resolved population models formulated by a system of differential equations that approximate the dynamics of two STPP-generated summary statistics: first-order spatial cumulants (densities), and second-order spatial cumulants (spatial covariances). We exemplify how SCMs can be used in mathematical oncology by modelling theoretical cancer cell populations comprising interacting growth factor-producing and non-producing cells. To formulate model equations, we use computational tools that enable the generation of STPPs, SCMs and mean-field population models (MFPMs) from user-defined model descriptions (Cornell et al. Nat Commun 10:4716, 2019). To calculate and compare STPP, SCM and MFPM-generated summary statistics, we develop an application-agnostic computational pipeline. Our results demonstrate that SCMs can capture STPP-generated population density dynamics, even when MFPMs fail to do so. From both MFPM and SCM equations, we derive treatment-induced death rates required to achieve non-growing cell populations. When testing these treatment strategies in STPP-generated cell populations, our results demonstrate that SCM-informed strategies outperform MFPM-informed strategies in terms of inhibiting population growths. We thus demonstrate that SCMs provide a new framework in which to study cell-cell interactions, and can be used to describe and perturb STPP-generated cell population dynamics. We, therefore, argue that SCMs can be used to increase IBMs' applicability in cancer research.

A molecular-based identification resource for the arthropods of Finland.

To associate specimens identified by molecular characters to other biological knowledge, we need reference sequences annotated by Linnaean taxonomy. In this study, we (1) report the creation of a comprehensive reference library of DNA barcodes for the arthropods of an entire country (Finland), (2) publish this library, and (3) deliver a new identification tool for insects and spiders, as based on this resource. The reference library contains mtDNA COI barcodes for 11,275 (43%) of 26,437 arthropod species known from Finland, including 10,811 (45%) of 23,956 insect species. To quantify the improvement in identification accuracy enabled by the current reference library, we ran 1000 Finnish insect and spider species through the Barcode of Life Data system (BOLD) identification engine. Of these, 91% were correctly assigned to a unique species when compared to the new reference library alone, 85% were correctly identified when compared to BOLD with the new material included, and 75% with the new material excluded. To capitalize on this resource, we used the new reference material to train a probabilistic taxonomic assignment tool, FinPROTAX, scoring high success. For the full-length barcode region, the accuracy of taxonomic assignments at the level of classes, orders, families, subfamilies, tribes, genera, and species reached 99.9%, 99.9%, 99.8%, 99.7%, 99.4%, 96.8%, and 88.5%, respectively. The FinBOL arthropod reference library and FinPROTAX are available through the Finnish Biodiversity Information Facility (www.laji.fi) at https://laji.fi/en/theme/protax. Overall, the FinBOL investment represents a massive capacity-transfer from the taxonomic community of Finland to all sectors of society.

Alternative developmental and transcriptomic responses to host plant water limitation in a butterfly metapopulation.

Predicting how climate change affects biotic interactions poses a challenge. Plant-insect herbivore interactions are particularly sensitive to climate change, as climate-induced changes in plant quality cascade into the performance of insect herbivores. Whereas the immediate survival of herbivore individuals depends on plastic responses to climate change-induced nutritional stress, long-term population persistence via evolutionary adaptation requires genetic variation for these responses. To assess the prospects for population persistence under climate change, it is therefore crucial to characterize response mechanisms to climate change-induced stressors, and quantify their variability in natural populations. Here, we test developmental and transcriptomic responses to water limitation-induced host plant quality change in a Glanville fritillary butterfly (Melitaea cinxia) metapopulation. We combine nuclear magnetic resonance spectroscopy on the plant metabolome, larval developmental assays and an RNA sequencing analysis of the larval transcriptome. We observed that responses to feeding on water-limited plants, in which amino acids and aromatic compounds are enriched, showed marked variation within the metapopulation, with individuals of some families performing better on control and others on water-limited plants. The transcriptomic responses were concordant with the developmental responses: families exhibiting opposite developmental responses also produced opposite transcriptomic responses (e.g. in growth-associated transcripts). The divergent responses in both larval development and transcriptome are associated with differences between families in amino acid catabolism and storage protein production. The results reveal intrapopulation variability in plasticity, suggesting that the Finnish M. cinxia metapopulation harbours potential for buffering against drought-induced changes in host plant quality.

High resolution DNA barcode library for European butterflies reveals continental patterns of mitochondrial genetic diversity.

The study of global biodiversity will greatly benefit from access to comprehensive DNA barcode libraries at continental scale, but such datasets are still very rare. Here, we assemble the first high-resolution reference library for European butterflies that provides 97% taxon coverage (459 species) and 22,306 COI sequences. We estimate that we captured 62% of the total haplotype diversity and show that most species possess a few very common haplotypes and many rare ones. Specimens in the dataset have an average 95.3% probability of being correctly identified. Mitochondrial diversity displayed elevated haplotype richness in southern European refugia, establishing the generality of this key biogeographic pattern for an entire taxonomic group. Fifteen percent of the species are involved in barcode sharing, but two thirds of these cases may reflect the need for further taxonomic research. This dataset provides a unique resource for conservation and for studying evolutionary processes, cryptic species, phylogeography, and ecology.

Transcriptomic and Phenotypic Analyses of the Sigma B-Dependent Characteristics and the Synergism between Sigma B and Sigma L in Listeria monocytogenes EGD-e.

Numerous gene expression and stress adaptation responses in L. monocytogenes are regulated through alternative sigma factors σB and σL. Stress response phenotypes and transcriptomes were compared between L. monocytogenes EGD-e and its ΔsigB and ΔsigBL mutants. Targeted growth phenotypic analysis revealed that the ΔsigB and ΔsigBL mutants are impaired during growth under cold and organic-acid stress conditions. Phenotypic microarrays revealed increased sensitivity in both mutants to various antimicrobial compounds. Genes de-regulated in these two mutants were identified by genome-wide transcriptome analysis during exponential growth in BHI. The ΔsigB and ΔsigBL strains repressed 198 and 254 genes, respectively, compared to the parent EGD-e strain at 3 °C, whereas 86 and 139 genes, respectively, were repressed in these mutants during growth at 37 °C. Genes repressed in these mutants are involved in various cellular functions including transcription regulation, energy metabolism and nutrient transport functions, and viral-associated processes. Exposure to cold stress induced a significant increase in σB and σL co-dependent genes of L. monocytogenes EGD-e since most (62%) of the down-regulated genes uncovered at 3 °C were detected in the ΔsigBL double-deletion mutant but not in ΔsigB or ΔsigL single-deletion mutants. Overall, the current study provides an expanded insight into σB and σL phenotypic roles and functional interactions in L. monocytogenes. Besides previously known σB- and σL-dependent genes, the transcriptomes defined in ΔsigB and ΔsigBL mutants reveal several new genes that are positively regulated by σB alone, as well as those co-regulated through σB- and σL-dependent mechanisms during L. monocytogenes growth under optimal and cold-stress temperature conditions.

A general mathematical method for predicting spatio-temporal correlations emerging from agent-based models.

Agent-based models are used to study complex phenomena in many fields of science. While simulating agent-based models is often straightforward, predicting their behaviour mathematically has remained a key challenge. Recently developed mathematical methods allow the prediction of the emerging spatial patterns for a general class of agent-based models, whereas the prediction of spatio-temporal pattern has been thus far achieved only for special cases. We present a general and mathematically rigorous methodology that allows deriving the spatio-temporal correlation structure for a general class of individual-based models. To do so, we define an auxiliary model, in which each agent type of the primary model expands to three types, called the original, the past and the new agents. In this way, the auxiliary model keeps track of both the initial and current state of the primary model, and hence the spatio-temporal correlations of the primary model can be derived from the spatial correlations of the auxiliary model. We illustrate the agreement between analytical predictions and agent-based simulations using two example models from theoretical ecology. In particular, we show that the methodology is able to correctly predict the dynamical behaviour of a host-parasite model that shows spatially localized oscillations.

Higher host plant specialization of root-associated endophytes than mycorrhizal fungi along an arctic elevational gradient.

How community-level specialization differs among groups of organisms, and changes along environmental gradients, is fundamental to understanding the mechanisms influencing ecological communities. In this paper, we investigate the specialization of root-associated fungi for plant species, asking whether the level of specialization varies with elevation. For this, we applied DNA barcoding based on the ITS region to root samples of five plant species equivalently sampled along an elevational gradient at a high arctic site. To assess whether the level of specialization changed with elevation and whether the observed patterns varied between mycorrhizal and endophytic fungi, we applied a joint species distribution modeling approach. Our results show that host plant specialization is not environmentally constrained in arctic root-associated fungal communities, since there was no evidence for changing specialization with elevation, even if the composition of root-associated fungal communities changed substantially. However, the level of specialization for particular plant species differed among fungal groups, root-associated endophytic fungal communities being highly specialized on particular host species, and mycorrhizal fungi showing almost no signs of specialization. Our results suggest that plant identity affects associated mycorrhizal and endophytic fungi differently, highlighting the need of considering both endophytic and mycorrhizal fungi when studying specialization in root-associated fungal communities.

Fungal communities decline with urbanization-more in air than in soil.

Increasing evidence suggests that degradation of biodiversity in human populated areas is a threat for the ecosystem processes that are relevant for human well-being. Fungi are a megadiverse kingdom that plays a key role in ecosystem processes and affects human well-being. How urbanization influences fungi has remained poorly understood, partially due to the methodological difficulties in comprehensively surveying fungi. Here we show that both aerial and soil fungal communities are greatly poorer in urban than in natural areas. Strikingly, a fivefold reduction in fungal DNA abundance took place in both air and soil samples already at 1 km scale when crossing the edge from natural to urban habitats. Furthermore, in the air, fungal diversity decreased with urbanization even more than in the soil. This result is counterintuitive as fungal spores are known to disperse over large distances. A large proportion of the fungi detectable in the air are specialized to natural habitats, whereas soil fungal communities comprise a large proportion of habitat generalists. The sensitivity of the aerial fungal community to anthropogenic disturbance makes this method a reliable and efficient bioindicator of ecosystem health in urban areas.

Accounting for environmental variation in co-occurrence modelling reveals the importance of positive interactions in root-associated fungal communities.

Understanding the role of interspecific interactions in shaping ecological communities is one of the central goals in community ecology. In fungal communities, measuring interspecific interactions directly is challenging because these communities are composed of large numbers of species, many of which are unculturable. An indirect way of assessing the role of interspecific interactions in determining community structure is to identify the species co-occurrences that are not constrained by environmental conditions. In this study, we investigated co-occurrences among root-associated fungi, asking whether fungi co-occur more or less strongly than expected based on the environmental conditions and the host plant species examined. We generated molecular data on root-associated fungi of five plant species evenly sampled along an elevational gradient at a high arctic site. We analysed the data using a joint species distribution modelling approach that allowed us to identify those co-occurrences that could be explained by the environmental conditions and the host plant species, as well as those co-occurrences that remained unexplained and thus more probably reflect interactive associations. Our results indicate that not only negative but also positive interactions play an important role in shaping microbial communities in arctic plant roots. In particular, we found that mycorrhizal fungi are especially prone to positively co-occur with other fungal species. Our results bring new understanding to the structure of arctic interaction networks by suggesting that interactions among root-associated fungi are predominantly positive.

A unified framework for analysis of individual-based models in ecology and beyond.

Individual-based models, 'IBMs', describe naturally the dynamics of interacting organisms or social or financial agents. They are considered too complex for mathematical analysis, but computer simulations of them cannot give the general insights required. Here, we resolve this problem with a general mathematical framework for IBMs containing interactions of an unlimited level of complexity, and derive equations that reliably approximate the effects of space and stochasticity. We provide software, specified in an accessible and intuitive graphical way, so any researcher can obtain analytical and simulation results for any particular IBM without algebraic manipulation. We illustrate the framework with examples from movement ecology, conservation biology, and evolutionary ecology. This framework will provide unprecedented insights into a hitherto intractable panoply of complex models across many scientific fields.

Effect of wet storage conditions on potato tuber transcriptome, phytohormones and growth.

BACKGROUND: Stored potato (Solanum tuberosum L.) tubers are sensitive to wet conditions that can cause rotting in long-term storage. To study the effect of water on the tuber surface during storage, microarray analysis, RNA-Seq profiling, qRT-PCR and phytohormone measurements were performed to study gene expression and hormone content in wet tubers incubated at two temperatures: 4 °C and 15 °C. The growth of the plants was also observed in a greenhouse after the incubation of tubers in wet conditions. RESULTS: Wet conditions induced a low-oxygen response, suggesting reduced oxygen availability in wet tubers at both temperatures when compared to that in the corresponding dry samples. Wet conditions induced genes coding for heat shock proteins, as well as proteins involved in fermentative energy production and defense against reactive oxygen species (ROS), which are transcripts that have been previously associated with low-oxygen stress in hypoxic or anoxic conditions. Wet treatment also induced senescence-related gene expression and genes involved in cell wall loosening, but downregulated genes encoding protease inhibitors and proteins involved in chloroplast functions and in the biosynthesis of secondary metabolites. Many genes involved in the production of phytohormones and signaling were also affected by wet conditions, suggesting altered regulation of growth by wet conditions. Hormone measurements after incubation showed increased salicylic acid (SA), abscisic acid (ABA) and auxin (IAA) concentrations as well as reduced production of jasmonate 12-oxo-phytodienoic acid (OPDA) in wet tubers. After incubation in wet conditions, the tubers produced fewer stems and more roots compared to controls incubated in dry conditions. CONCLUSIONS: In wet conditions, tubers invest in ROS protection and defense against the abiotic stress caused by reduced oxygen due to excessive water. Changes in ABA, SA and IAA that are antagonistic to jasmonates affect growth and defenses, causing induction of root growth and rendering tubers susceptible to necrotrophic pathogens. Water on the tuber surface may function as a signal for growth, similar to germination of seeds.

Scaling up the effects of inbreeding depression from individuals to metapopulations.

Inbreeding is common in nature, and many laboratory studies have documented that inbreeding depression can reduce the fitness of individuals. Demonstrating the consequences of inbreeding depression on the growth and persistence of populations is more challenging because populations are often regulated by density- or frequency-dependent selection and influenced by demographic and environmental stochasticity. A few empirical studies have shown that inbreeding depression can increase extinction risk of local populations. The importance of inbreeding depression at the metapopulation level has been conjectured based on population-level studies but has not been evaluated. We quantified the impact of inbreeding depression affecting the fitness of individuals on metapopulation persistence in heterogeneous habitat networks of different sizes and habitat configuration in a context of natural butterfly metapopulations. We developed a spatial individual-based simulation model of metapopulations with explicit genetics. We used Approximate Bayesian Computation to fit the model to extensive demographic, genetic and life-history data available for the well-studied Glanville fritillary butterfly (Melitaea cinxia) metapopulations in the Åland islands in SW Finland. We compared 18 semi-independent habitat networks differing in size and fragmentation. The results show that inbreeding is more frequent in small habitat networks, and consequently, inbreeding depression elevates extinction risks in small metapopulations. Metapopulation persistence and neutral genetic diversity maintained in the metapopulations increase with the total habitat amount in and mean patch size of habitat networks. Dispersal and mating behaviour interact with landscape structure to determine how likely it is to encounter kin while looking for mates. Inbreeding depression can decrease the viability of small metapopulations even when they are strongly influenced by stochastic extinction-colonization dynamics and density-dependent selection. The findings from this study support that genetic factors, in addition to demographic factors, can contribute to extinctions of small local populations and also of metapopulations.

BARCOSEL: a tool for selecting an optimal barcode set for high-throughput sequencing.

BACKGROUND: Current high-throughput sequencing platforms provide capacity to sequence multiple samples in parallel. Different samples are labeled by attaching a short sample specific nucleotide sequence, barcode, to each DNA molecule prior pooling them into a mix containing a number of libraries to be sequenced simultaneously. After sequencing, the samples are binned by identifying the barcode sequence within each sequence read. In order to tolerate sequencing errors, barcodes should be sufficiently apart from each other in sequence space. An additional constraint due to both nucleotide usage and basecalling accuracy is that the proportion of different nucleotides should be in balance in each barcode position. The number of samples to be mixed in each sequencing run may vary and this introduces a problem how to select the best subset of available barcodes at sequencing core facility for each sequencing run. There are plenty of tools available for de novo barcode design, but they are not suitable for subset selection. RESULTS: We have developed a tool which can be used for three different tasks: 1) selecting an optimal barcode set from a larger set of candidates, 2) checking the compatibility of user-defined set of barcodes, e.g. whether two or more libraries with existing barcodes can be combined in a single sequencing pool, and 3) augmenting an existing set of barcodes. In our approach the selection process is formulated as a minimization problem. We define the cost function and a set of constraints and use integer programming to solve the resulting combinatorial problem. Based on the desired number of barcodes to be selected and the set of candidate sequences given by user, the necessary constraints are automatically generated and the optimal solution can be found. The method is implemented in C programming language and web interface is available at http://ekhidna2.biocenter.helsinki.fi/barcosel . CONCLUSIONS: Increasing capacity of sequencing platforms raises the challenge of mixing barcodes. Our method allows the user to select a given number of barcodes among the larger existing barcode set so that both sequencing errors are tolerated and the nucleotide balance is optimized. The tool is easy to access via web browser.

Protax-fungi: a web-based tool for probabilistic taxonomic placement of fungal internal transcribed spacer sequences.

Incompleteness of reference sequence databases and unresolved taxonomic relationships complicates taxonomic placement of fungal sequences. We developed Protax-fungi, a general tool for taxonomic placement of fungal internal transcribed spacer (ITS) sequences, and implemented it into the PlutoF platform of the UNITE database for molecular identification of fungi. With empirical data on root- and wood-associated fungi, Protax-fungi reliably identified (with at least 90% identification probability) the majority of sequences to the order level but only around one-fifth of them to the species level, reflecting the current limited coverage of the databases. Protax-fungi outperformed the Sintax and Rdb classifiers in terms of increased accuracy and decreased calibration error when applied to data on mock communities representing species groups with poor sequence database coverage. We applied Protax-fungi to examine the internal consistencies of the Index Fungorum and UNITE databases. This revealed inconsistencies in the taxonomy database as well as mislabelling and sequence quality problems in the reference database. The according improvements were implemented in both databases. Protax-fungi provides a robust tool for performing statistically reliable identifications of fungi in spite of the incompleteness of extant reference sequence databases and unresolved taxonomic relationships.