Molecular operational taxonomic units as approximations of species in the light of evolutionary models and empirical data from Fungi

Amanita satotamagotake sp. nov., a cryptic species formerly included in Amanita caesareoides

We evaluated the inclusion of a cryptic species in a Japanese Amanita caesareoides population. We sampled A. caesareoides specimens under various vegetation and climate conditions, and then conducted phylogenetic analyses on sequences from seven loci. The A. caesareoides specimens showed two distinct groups, except when the ITS phylogeny was considered. These two phylogroups showed different distributions: subalpine-cool temperate and temperate-subtropical areas. Although these two phylogroups overlapped in terms of basidiospore size, the latter tended to exhibit smaller basidiospores. In addition, only the former showed mycelial growth on nutrient agar. Based on these phylo-morpho-ecophysiological characteristics, we separated the specimens labeled with the name A. caesareoides into two species. As the lectotype of A. caesareoides showed similarity to the former by DNA analysis, the latter was described as a new species, namely A. satotamagotake. Based on the geographic patterns of the two species, A. satotamagotake may have invaded the natural habit of A. caesareoides because of global warming.

Pseudocryptic diversity and species boundaries in the sea cucumber Stichopus cf. horrens (Echinodermata: Stichopodidae) revealed by mitochondrial and microsatellite markers

Morphologically cryptic and pseudo-cryptic species pose a challenge to taxonomic identification and assessments of species diversity and distributions. Such is the case for the sea cucumber Stichopus horrens, commonly confused with Stichopus monotuberculatus. Here, we used mitochondrial cytochrome oxidase subunit I (COI) and microsatellite markers to examine genetic diversity in Stichopus cf. horrens throughout the Philippine archipelago, to aid species identification and clarify species boundaries. Phylogenetic analysis reveals two recently diverged COI lineages (Clade A and Clade B; c. 1.35-2.54 Mya) corresponding to sequence records for specimens identified as S. monotuberculatus and S. horrens, respectively. Microsatellite markers reveal two significantly differentiated genotype clusters broadly concordant with COI lineages (Cluster 1, Cluster 2). A small proportion of individuals were identified as later-generation hybrids indicating limited contemporary gene flow between genotype clusters, thus confirming species boundaries. Morphological differences in papillae distribution and form are observed for the two species, however tack-like spicules from the dorsal papillae are not a reliable diagnostic character. An additional putative cryptic species was detected within Clade B-Cluster 2 specimens warranting further examination. We propose that these lineages revealed by COI and genotype data be referred to as Stichopus cf. horrens species complex.

Does One Size Fit All? Variations in the DNA Barcode Gaps of Macrofungal Genera

The nuclear ribosomal internal transcribed spacer (nrITS) region has been widely used in fungal diversity studies. Environmental metabarcoding has increased the importance of the fungal DNA barcode in documenting fungal diversity and distribution. The DNA barcode gap is seen as the difference between intra- and inter-specific pairwise distances in a DNA barcode. The current understanding of the barcode gap in macrofungi is limited, inhibiting the development of best practices in applying the nrITS region toward research on fungal diversity. This study examined the barcode gap using 5146 sequences representing 717 species of macrofungi from eleven genera, eight orders and two phyla in datasets assembled by taxonomic experts. Intra- and inter-specific pairwise distances were measured from sequence and phylogenetic data. The results demonstrate that barcode gaps are influenced by differences in intra- and inter-specific variance in pairwise distances. In terms of DNA barcode behavior, variance is greater in the ITS1 than ITS2, and variance is greater in both relative to the combined nrITS region. Due to the difference in variance, the barcode gaps in the ITS2 region are greater than in the ITS1. Additionally, the taxonomic approach of "splitting" taxa into numerous taxonomic units produces greater barcode gaps when compared to "lumping". The results show variability in the barcode gaps between fungal taxa, demonstrating a need to understand the accuracy of DNA barcoding in quantifying species richness. For taxonomic studies, variability in nrITS sequence data supports the application of multiple molecular markers to corroborate the taxonomic and systematic delineation of species.

Effects of operational taxonomic unit inference methods on soil microeukaryote community analysis using long-read metabarcoding

Long amplicon metabarcoding has opened the door for phylogenetic analysis of the largely unknown communities of microeukaryotes in soil. Here, we amplified and sequenced the ITS and LSU regions of the rDNA operon (around 1500 bp) from grassland soils using PacBio SMRT sequencing. We tested how three different methods for generation of operational taxonomic units (OTUs) effected estimated richness and identified taxa, and how well large-scale ecological patterns associated with shifting environmental conditions were recovered in data from the three methods. The field site at Kungsängen Nature Reserve has drawn frequent visitors since Linnaeus's time, and its species rich vegetation includes the largest population of Fritillaria meleagris in Sweden. To test the effect of different OTU generation methods, we sampled soils across an abrupt moisture transition that divides the meadow community into a Carex acuta dominated plant community with low species richness in the wetter part, which is visually distinct from the mesic-dry part that has a species rich grass-dominated plant community including a high frequency of F. meleagris. We used the moisture and plant community transition as a framework to investigate how detected belowground microeukaryotic community composition was influenced by OTU generation methods. Soil communities in both moisture regimes were dominated by protists, a large fraction of which were taxonomically assigned to Ciliophora (Alveolata) while 30%-40% of all reads were assigned to kingdom Fungi. Ecological patterns were consistently recovered irrespective of OTU generation method used. However, different methods strongly affect richness estimates and the taxonomic and phylogenetic resolution of the characterized community with implications for how well members of the microeukaryotic communities can be recognized in the data.

DNA metabarcoding and morphological methods show complementary patterns in the metacommunity organization of lentic epiphytic diatoms

Diatoms are important organisms in freshwater ecosystems due to their position as primary producers and therefore, analyzing their assemblages provides relevant information on ecosystem functioning. Diatoms have historically been identified based on morphological traits, which is time-consuming and requires well-trained specialists. Nevertheless, DNA barcoding offers an alternative approach to overcome some limitations of the morphological method. Here, we assess if both approaches are comparable methods to study patterns and mechanisms (including environmental filtering and dispersal limitation) of epiphytic diatom metacommunities using a comprehensive dataset from 22 Mediterranean ponds at different taxonomic resolutions. We used a fragment of rbcL barcode gene combined with High-Throughput Sequencing to infer diatom community composition. The overall degree of correspondence between both approaches was assessed by Procrustean rotation analysis and Procrustean randomization tests, whereas the role of local environmental variables and geographical distances was studied using a comprehensive combination of BIOENV, Mantel tests and distance-based redundancy analysis. Our results showed a relatively poor correspondence in the compositional variation of diatom metacommunity between both approaches. We speculate that the incompleteness of the reference database and the bioinformatics processing are the biases most likely affecting the molecular approach, whereas the limited counting effort and the presence of cryptic species are presumably the major biases related with the morphological method. On the other hand, variation in diatom community composition detected with both approaches was strongly related to the environmental template, which may be related with the narrow community-environment relationships in diatoms. Nevertheless, we found no significant relationship between compositional variation and geographical distances. Overall, our work shows the complementary nature of both approaches and highlights the importance of DNA metabarcoding to address empirical research questions of community ecology in freshwaters, especially once the reference databases include most genotypes of occurring taxa and bioinformatics biases are overcome.

Enlightening the black and white: species delimitation and UNITE species hypothesis testing in the Russula albonigra species complex

Russula albonigra is considered a well-known species, morphologically delimited by the context of the basidiomata blackening without intermediate reddening, and the menthol-cooling taste of the lamellae. It is supposed to have a broad ecological range and a large distribution area. A thorough molecular analysis based on four nuclear markers (ITS, LSU, RPB2 and TEF1-α) shows this traditional concept of R. albonigra s. lat. represents a species complex consisting of at least five European, three North American, and one Chinese species. Morphological study shows traditional characters used to delimit R. albonigra are not always reliable. Therefore, a new delimitation of the R. albonigra complex is proposed and a key to the described European species of R. subgen. Compactae is presented. A lectotype and an epitype are designated for R. albonigra and three new European species are described: R. ambusta, R. nigrifacta, and R. ustulata. Different thresholds of UNITE species hypotheses were tested against the taxonomic data. The distance threshold of 0.5% gives a perfect match to the phylogenetically defined species within the R. albonigra complex. Publicly available sequence data can contribute to species delimitation and increase our knowledge on ecology and distribution, but the pitfalls are short and low quality sequences.

'Species' without species

Biological science uses multiple species concepts. Order can be brought to this diversity if we recognize two key features. First, any given species concept is likely to have a patchwork structure, generated by repeated application of the concept to new domains. We illustrate this by showing how two species concepts (biological and ecological) have been modified from their initial eukaryotic applications to apply to prokaryotes. Second, both within and between patches, distinct species concepts may interact and hybridize. We thus defend a semantic picture of the species concept as a collection of interacting patchwork structures. Thus, although not all uses of the term pick out the same kind of unit in nature, the diversity of uses reflects something more than mere polysemy. We suggest that the emphasis on the use of species to pick out natural units is itself problematic, because that is not the term's sole function. In particular, species concepts are used to manage inquiry into processes of speciation, even when these processes do not produce clearly delimited species.

Long- and short-read metabarcoding technologies reveal similar spatiotemporal structures in fungal communities

Fungi form diverse communities and play essential roles in many terrestrial ecosystems, yet there are methodological challenges in taxonomic and phylogenetic placement of fungi from environmental sequences. To address such challenges, we investigated spatiotemporal structure of a fungal community using soil metabarcoding with four different sequencing strategies: short-amplicon sequencing of the ITS2 region (300-400 bp) with Illumina MiSeq, Ion Torrent Ion S5 and PacBio RS II, all from the same PCR library, as well as long-amplicon sequencing of the full ITS and partial LSU regions (1200-1600 bp) with PacBio RS II. Resulting community structure and diversity depended more on statistical method than sequencing technology. The use of long-amplicon sequencing enables construction of a phylogenetic tree from metabarcoding reads, which facilitates taxonomic identification of sequences. However, long reads present issues for denoising algorithms in diverse communities. We present a solution that splits the reads into shorter homologous regions prior to denoising, and then reconstructs the full denoised reads. In the choice between short and long amplicons, we suggest a hybrid approach using short amplicons for sampling breadth and depth, and long amplicons to characterize the local species pool for improved identification and phylogenetic analyses.

Naming the untouchable - environmental sequences and niche partitioning as taxonomical evidence in fungi

Due to their submerged and cryptic lifestyle, the vast majority of fungal species are difficult to observe and describe morphologically, and many remain known to science only from sequences detected in environmental samples. The lack of practices to delimit and name most fungal species is a staggering limitation to communication and interpretation of ecology and evolution in kingdom Fungi. Here, we use environmental sequence data as taxonomical evidence and combine phylogenetic and ecological data to generate and test species hypotheses in the class Archaeorhizomycetes (Taphrinomycotina, Ascomycota). Based on environmental amplicon sequencing from a well-studied Swedish pine forest podzol soil, we generate 68 distinct species hypotheses of Archaeorhizomycetes, of which two correspond to the only described species in the class. Nine of the species hypotheses represent 78% of the sequenced Archaeorhizomycetes community, and are supported by long read data that form the backbone for delimiting species hypothesis based on phylogenetic branch lengths.

Soil fungal communities are shaped by environmental filtering and competitive exclusion so that closely related species are less likely to co-occur in a niche if adaptive traits are evolutionarily conserved. In soil profiles, distinct vertical horizons represent a testable niche dimension, and we found significantly differential distribution across samples for a well-supported pair of sister species hypotheses. Based on the combination of phylogenetic and ecological evidence, we identify two novel species for which we provide molecular diagnostics and propose names. While environmental sequences cannot be automatically translated to species, they can be used to generate phylogenetically distinct species hypotheses that can be further tested using sequences as ecological evidence. We conclude that in the case of abundantly and frequently observed species, environmental sequences can support species recognition in the absences of physical specimens, while rare taxa remain uncaptured at our sampling and sequencing intensity.

Brazilian fungal diversity represented by DNA markers generated over 20 years

Molecular techniques using fungal DNA barcoding (ITS) and other markers have been key to identifying the biodiversity of different geographic areas, mainly in megadiverse countries. Here, we provide an overview of the fungal diversity in Brazil based on DNA markers of phylogenetic importance generated since 1996. We retrieved fungal sequences of ITS, LSU, SSU, tef1-α, β-tubulin, rpb1, rpb2, actin, chitin synthase, and ATP6 from GenBank using different field keywords that indicated their origin in Brazil. A total of 19,440 sequences were recovered. ITS is the most representative marker (11,209 sequences), with 70.1% belonging to Ascomycota, 18.6% Basidiomycota, 10.2% unidentified, 1.1% Mucoromycota, two sequences of Olpidium bornovanus (Fungi incertae sedis), one sequence of Blastocladiomycota (Allomyces arbusculus), and one sequence of Chytridiomycota (Batrachochytrium dendrobatidis). Considering the sequences of all selected markers, only the phyla Cryptomycota and Entorrhizomycota were not represented. Based on ITS, using a cutoff of 98%, all sequences comprise 3047 OTUs, with the majority being Ascomycota (2088 OTUs) and Basidiomycota (681 OTUs). Previous numbers based mainly on morphological and bibliographical data revealed 5264 fungal species from Brazil, with a predominance of Basidiomycota (2741 spp.) and Ascomycota (1881 spp.). The unidentified ITS sequences not assigned to a higher taxonomic level represent 1.61% of all ITS sequences sampled and correspond to 38 unknown class-level lineages (75% cutoff). A maximum likelihood phylogeny based on LSU illustrates the fungal classes occurring in Brazil.

Coalescent-based delimitation and species-tree estimations reveal Appalachian origin and Neogene diversification in Russula subsection Roseinae

Numerous lineages of mushroom-forming fungi have been subject to bursts of diversification throughout their evolutionary history, events that can impact our ability to infer well-resolved phylogenies. However, groups that have undergone quick genetic change may have the highest adaptive potential. As the second largest genus of mushroom-forming fungi, Russula provides an excellent model for studying hyper-diversification and processes in evolution that drives it. This study focuses on the morphologically defined group - Russula subsection Roseinae. Species hypotheses based on morphological differentiation and multi-locus phylogenetic analyses are tested in the Roseinae using different applications of the multi-species coalescent model. Based on this combined approach, we recognize fourteen species in Roseinae including the Albida and wholly novel Magnarosea clades. Reconstruction of biogeographic and host association history suggest that parapatric speciation in refugia during glacial cycles of the Pleistocene drove diversification within the Roseinae, which is found to have a Laurasian distribution with an evolutionary origin in the Appalachian Mountains of eastern North America. Finally, we detect jump dispersal at a continental scale that has driven diversification since the most recent glacial cycles.

Tissue-Specific and Geographical Variation in Endophytic Fungi of Ageratina adenophora and Fungal Associations With the Environment

To understand the distribution of the cultivable fungal community in plant tissues and the associations of these fungi with their surrounding environments during the geographical expansion of an invasive plant, Ageratina adenophora, we isolated the cultivable fungi from 72 plant tissues, 12 soils, and 12 air samples collected from six areas in Yunnan Province, China. A total of 4066 isolates were investigated, including 1641 endophytic fungi, 233 withered leaf fungi, 1255 fungi from air, and 937 fungi from soil. These fungi were divided into 458 and 201 operational taxonomic units (OTUs) with unique and 97% ITS gene sequence identity, respectively. Phylogenetic analysis showed that the fungi belonged to four phyla, including Ascomycota (94.20%), Basidiomycota (2.71%), Mortierellomycota (3.03%), and Mucoromycota (0.07%). The dominant genera of cultivable endophytic fungi were Colletotrichum (34.61%), Diaporthe (17.24%), Allophoma (8.03%), and Fusarium (4.44%). Colletotrichum and Diaporthe were primarily isolated from mature leaves, Allophoma from stems, and Fusarium from roots, indicating that the enrichment of endophytic fungi is tissue-specific and fungi rarely grew systemically within A. adenophora. In the surrounding environment, Alternaria (21.46%), Allophoma (19.31%), Xylaria (18.45%), and Didymella (18.03%) were dominant in the withered leaves, Cladosporium (22.86%), Trichoderma (14.27%), and Epicoccum (9.83%) were dominant in the canopy air, and Trichoderma (27.27%) and Mortierella (20.46%) were dominant in the rhizosphere soils. Further analysis revealed that the cultivable endophytic fungi changed across geographic areas and showed a certain degree of variation in different tissues of A. adenophora. The cultivable fungi in mature and withered leaves fluctuated more than those in roots and stems. We also found that some cultivable endophytic fungi might undergo tissue-to-tissue migration and that the stem could be a transport tissue by which airborne fungi infect roots. Finally, we provided evidence that the fungal community within A. adenophora was partially shared with the contiguous environment. The data suggested a frequent interaction between fungi associated with A. adenophora and those in surrounding environments, reflecting a compromise driven by both functional requirements for plant growth and local environmental conditions.

Phylogenetic study documents different speciation mechanisms within the Russula globispora lineage in boreal and arctic environments of the Northern Hemisphere

The Russula globispora lineage is a morphologically and phylogenetically well-defined group of ectomycorrhizal fungi occurring in various climatic areas. In this study we performed a multi-locus phylogenetic study based on collections from boreal, alpine and arctic habitats of Europe and Western North America, subalpine collections from the southeast Himalayas and collections from subtropical coniferous forests of Pakistan. European and North American collections are nearly identical and probably represent a single species named R. dryadicola distributed from the Alps to the Rocky Mountains. Collections from the southeast Himalayas belong to two distinct species: R. abbottabadensis sp. nov. from subtropical monodominant forests of Pinus roxburghii and R. tengii sp. nov. from subalpine mixed forests of Abies and Betula. The results suggest that speciation in this group is driven by a climate disjunction and adaptation rather than a host switch and geographical distance.

Chromosomal polymorphism and molecular variability in the pearly razorfish Xyrichtys novacula (Labriformes, Labridae): taxonomic and biogeographic implications

The pearly razorfish Xyrichtys novacula (Linnaeus, 1758) is a sedentary benthic species distributed in both sides of the Atlantic Ocean and in the Mediterranean Sea. Previous cytogenetic analysis reported different diploid numbers in samples from Italy, Venezuela and Brazil. This research aims to test the hypothesis that samples from American Atlantic coast and Mediterranean Sea belong to the same single evolutionary lineage, characterized by intra-specific chromosome polymorphism. To this purpose a cytogenetic and molecular (mitochondrial COI sequences) survey was undertaken. Results revealed the existence of three different pearly razorfish molecular lineages: one present in Mediterranean Sea and two in the central and south American area, which are characterized by different karyotypes. One of these lineages shows substantial intra-population chromosomal polymorphism (2n = 45-48) determined by Robertsonian fusions that produce large metacentric chromosomes. On the whole data suggest that specimens morphologically identified as X. novacula correspond to three cryptic species.

Balsamia (Sequestrate Helvellaceae, Ascomycota) in western North America

Balsamia, a hypogeous, sequestrate genus in the Helvellaceae, has been characterized variously as having three to eight species in North America, and these have been considered either different from or conspecific with European species. No available modern systematic treatment of Balsamia exists to allow for accurate identification at the species level. We sequenced DNA from recent western North American Balsamia collections, assessed relationships by sequence similarity, and identified molecular taxonomic units. From these data, we determined which matched descriptions and types of named species. ITS sequences supported 12 Balsamia species in western North America, five originally described by Harkness and Fischer and seven new species that we describe here. No sequences from Balsamia collections in western North America were nested among those of European species. We found no clear evidence for separation of Balsamia into multiple genera.

Ectomycorrhizal fungal communities in high mountain conifer forests in central Mexico and their potential use in the assisted migration of Abies religiosa

Abies religiosa forests in central Mexico are the only overwinter refuge of the monarch butterfly and provide important ecosystem services. These forests have lost 55% of their original area and as a consequence, diversity and biotic interactions in these ecosystems are in risk. The aim of this study was to compare the soil fungal diversity and community structure in the Abies religiosa forests and surrounding Pinus montezumae, Pinus hartwegii, and coniferous mixed forest plant communities to provide data on ecology of mycorrhizal interactions for the assisted migration of A. religiosa. We sampled soil from five coniferous forests, extracted total soil DNA, and sequenced the ITS2 region by Illumina MiSeq. The soil fungi community was integrated by 1746 taxa with a species turnover ranging from 0.280 to 0.461 between sampling sites. In the whole community, the more abundant and frequent species were Russula sp. (aff. olivobrunnea), Mortierella sp.1, and Piloderma sp. (aff. olivacearum). The ectomycorrhizal fungi were the more frequent and abundant functional group. A total of 298 species (84 ectomycorrhizal) was shared in the five conifer forests; these widely distributed species were dominated by Russulaceae and Clavulinaceae. The fungal community composition was significantly influenced by altitude and the lowest species turnover happened between the two A. religiosa forests even though they have different soil types. As Pinus montezumae forests have a higher altitudinal distribution adjacent to A. religiosa and share the largest number of ectomycorrhizal fungi with it, we suggest these forests as a potential habitat for new A. religiosa populations.

Ten reasons why a sequence-based nomenclature is not useful for fungi anytime soon

The large number of species still to be discovered in fungi, together with an exponentially growing number of environmental sequences that cannot be linked to known taxa, has fuelled the idea that it might be necessary to formally name fungi on the basis of sequence data only. Here we object to this idea due to several shortcomings of the approach, ranging from concerns regarding reproducibility and the violation of general scientific principles to ethical issues. We come to the conclusion that sequence-based nomenclature is potentially harmful for mycology as a discipline. Additionally, a classification based on sequences as types is not within reach anytime soon, because there is a lack of consensus regarding common standards due to the fast pace at which sequencing technologies develop.

Sequence clustering threshold has little effect on the recovery of microbial community structure

Analysis of microbial community structure by multivariate ordination methods, using data obtained by high-throughput sequencing of amplified markers (i.e., DNA metabarcoding), often requires clustering of DNA sequences into operational taxonomic units (OTUs). Parameters for the clustering procedure tend not to be justified but are set by tradition rather than being based on explicit knowledge. In this study, we explore the extent to which ordination results are affected by variation in parameter settings for the clustering procedure. Amplicon sequence data from nine microbial community studies, representing different sampling designs, spatial scales and ecosystems, were subjected to clustering into OTUs at seven different similarity thresholds (clustering thresholds) ranging from 87% to 99% sequence similarity. The 63 data sets thus obtained were subjected to parallel DCA and GNMDS ordinations. The resulting community structures were highly similar across all clustering thresholds. We explain this pattern by the existence of strong ecological structuring gradients and phylogenetically diverse sets of abundant OTUs that are highly stable across clustering thresholds. Removing low-abundance, rare OTUs had negligible effects on community patterns. Our results indicate that microbial data sets with a clear gradient structure are highly robust to choice of sequence clustering threshold.

New light on names and naming of dark taxa

A growing proportion of fungal species and lineages are known only from sequence data and cannot be linked to any physical specimen or resolved taxonomic name. Such fungi are often referred to as “dark taxa” or “dark matter fungi”. As they lack a taxonomic identity in the form of a name, they are regularly ignored in many important contexts, for example in legalisation and species counts. It is therefore very urgent to find a system to also deal with these fungi. Here, issues relating to the taxonomy and nomenclature of dark taxa are discussed and a number of questions that the mycological community needs to consider before deciding on what system/s to implement are highlighted.

Genetic barcoding of dark-spored myxomycetes (Amoebozoa)-Identification, evaluation and application of a sequence similarity threshold for species differentiation in NGS studies

Unicellular, eukaryotic organisms (protists) play a key role in soil food webs as major predators of microorganisms. However, due to the polyphyletic nature of protists, no single universal barcode can be established for this group, and the structure of many protistean communities remains unresolved. Plasmodial slime moulds (Myxogastria or Myxomycetes) stand out among protists by their formation of fruit bodies, which allow for a morphological species concept. By Sanger sequencing of a large collection of morphospecies, this study presents the largest database to date of dark-spored myxomycetes and evaluate a partial 18S SSU gene marker for species annotation. We identify and discuss the use of an intraspecific sequence similarity threshold of 99.1% for species differentiation (OTU picking) in environmental PCR studies (ePCR) and estimate a hidden diversity of putative species, exceeding those of described morphospecies by 99%. When applying the identified threshold to an ePCR data set (including sequences from both NGS and cloning), we find 64 OTUs of which 21.9% had a direct match (>99.1% similarity) to the database and the remaining had on average 90.2 ± 0.8% similarity to their best match, thus thought to represent undiscovered diversity of dark-spored myxomycetes.

Molecular operational taxonomic units of mosquitoes (Diptera: Culicidae) collected in high Andean mountain ecosystems of Antioquia, Colombia

Accurate taxonomic identification of highland mosquito species may be complicated because of the lack of comprehensive regional morphological keys and taxonomic specialists, particularly for mosquitoes of medical or ecological importance. We applied a multi-locus approach to explore the diversity of genera/species collected, to define the Molecular Operational Taxonomic Units (MOTUs) and to perform phylogenetic clustering. Twenty MOTUs and three single sequences were revealed from 78 concatenated cox1 + ITS2 sequences, and the species name was allocated for five of these. This study provides molecular taxonomic information of culicid fauna present in high Andean mountain ecosystems in Antioquia, Colombia. However, future morphological and integrative taxonomic studies should be conducted to achieve the specific identity of all detected MOTUs.

Environmental DNA (eDNA) metabarcoding assays to detect invasive invertebrate species in the Great Lakes

Describing and monitoring biodiversity comprise integral parts of ecosystem management. Recent research coupling metabarcoding and environmental DNA (eDNA) demonstrate that these methods can serve as important tools for surveying biodiversity, while significantly decreasing the time, expense and resources spent on traditional survey methods. The literature emphasizes the importance of genetic marker development, as the markers dictate the applicability, sensitivity and resolution ability of an eDNA assay. The present study developed two metabarcoding eDNA assays using the mtDNA 16S RNA gene with Illumina MiSeq platform to detect invertebrate fauna in the Laurentian Great Lakes and surrounding waterways, with a focus for use on invasive bivalve and gastropod species monitoring. We employed careful primer design and in vitro testing with mock communities to assess ability of the markers to amplify and sequence targeted species DNA, while retaining rank abundance information. In our mock communities, read abundances reflected the initial input abundance, with regressions having significant slopes (p<0.05) and high coefficients of determination (R2) for all comparisons. Tests on field environmental samples revealed similar ability of our markers to measure relative abundance. Due to the limited reference sequence data available for these invertebrate species, care must be taken when analyzing results and identifying sequence reads to species level. These markers extend eDNA metabarcoding research for molluscs and appear relevant to other invertebrate taxa, such as rotifers and bryozoans. Furthermore, the sphaeriid mussel assay is group-specific, exclusively amplifying bivalves in the Sphaeridae family and providing species-level identification. Our assays provide useful tools for managers and conservation scientists, facilitating early detection of invasive species as well as improving resolution of mollusc diversity.

Increased phylogenetic resolution within the ecologically important Rhizopogon subgenus Amylopogon using 10 anonymous nuclear loci

Rhizopogon species are ecologically significant ectomycorrhizal fungi in conifer ecosystems. The importance of this system merits the development and utilization of a more robust set of molecular markers specifically designed to evaluate their evolutionary ecology. Anonymous nuclear loci (ANL) were developed for R. subgenus Amylopogon. Members of this subgenus occur throughout the United States and are exclusive fungal symbionts associated with Pterospora andromedea, a threatened mycoheterotrophic plant endemic to disjunct eastern and western regions of North America. Candidate ANL were developed from 454 shotgun pyrosequencing and assessed for positive amplification across targeted species, sequencing success, and recovery of phylogenetically informative sites. Ten ANL were successfully developed and were subsequently used to sequence representative taxa, herbaria holotype and paratype specimens in R. subgenus Amylopogon. Phylogenetic reconstructions were performed on individual and concatenated data sets by Bayesian inference and maximum likelihood methods. Phylogenetic analyses of these 10 ANL were compared with a phylogeny traditionally constructed using the universal fungal barcode nuc rDNA ITS1-5.8S-ITS2 region (ITS). The resulting ANL phylogeny was consistent with most of the species designations delineated by ITS. However, the ANL phylogeny provided much greater phylogenetic resolution, yielding new evidence for cryptic species within previously defined species of R. subgenus Amylopogon. Additionally, the rooted ANL phylogeny provided an alternate topology to the ITS phylogeny, which inferred a novel set of evolutionary relationships not identified in prior phylogenetic studies.

Scaling up discovery of hidden diversity in fungi: impacts of barcoding approaches

The fungal kingdom is a hyperdiverse group of multicellular eukaryotes with profound impacts on human society and ecosystem function. The challenge of documenting and describing fungal diversity is exacerbated by their typically cryptic nature, their ability to produce seemingly unrelated morphologies from a single individual and their similarity in appearance to distantly related taxa. This multiplicity of hurdles resulted in the early adoption of DNA-based comparisons to study fungal diversity, including linking curated DNA sequence data to expertly identified voucher specimens. DNA-barcoding approaches in fungi were first applied in specimen-based studies for identification and discovery of taxonomic diversity, but are now widely deployed for community characterization based on sequencing of environmental samples. Collectively, fungal barcoding approaches have yielded important advances across biological scales and research applications, from taxonomic, ecological, industrial and health perspectives. A major outstanding issue is the growing problem of 'sequences without names' that are somewhat uncoupled from the traditional framework of fungal classification based on morphology and preserved specimens. This review summarizes some of the most significant impacts of fungal barcoding, its limitations, and progress towards the challenge of effective utilization of the exponentially growing volume of data gathered from high-throughput sequencing technologies.This article is part of the themed issue 'From DNA barcodes to biomes'.

Phylogenetically Structured Differences in rRNA Gene Sequence Variation among Species of Arbuscular Mycorrhizal Fungi and Their Implications for Sequence Clustering

Arbuscular mycorrhizal (AM) fungi form mutualisms with plant roots that increase plant growth and shape plant communities. Each AM fungal cell contains a large amount of genetic diversity, but it is unclear if this diversity varies across evolutionary lineages. We found that sequence variation in the nuclear large-subunit (LSU) rRNA gene from 29 isolates representing 21 AM fungal species generally assorted into genus- and species-level clades, with the exception of species of the genera Claroideoglomus and Entrophospora However, there were significant differences in the levels of sequence variation across the phylogeny and between genera, indicating that it is an evolutionarily constrained trait in AM fungi. These consistent patterns of sequence variation across both phylogenetic and taxonomic groups pose challenges to interpreting operational taxonomic units (OTUs) as approximations of species-level groups of AM fungi. We demonstrate that the OTUs produced by five sequence clustering methods using 97% or equivalent sequence similarity thresholds failed to match the expected species of AM fungi, although OTUs from AbundantOTU, CD-HIT-OTU, and CROP corresponded better to species than did OTUs from mothur or UPARSE. This lack of OTU-to-species correspondence resulted both from sequences of one species being split into multiple OTUs and from sequences of multiple species being lumped into the same OTU. The OTU richness therefore will not reliably correspond to the AM fungal species richness in environmental samples. Conservatively, this error can overestimate species richness by 4-fold or underestimate richness by one-half, and the direction of this error will depend on the genera represented in the sample.

IMPORTANCE

Arbuscular mycorrhizal (AM) fungi form important mutualisms with the roots of most plant species. Individual AM fungi are genetically diverse, but it is unclear whether the level of this diversity differs among evolutionary lineages. We found that the amount of sequence variation in an rRNA gene that is commonly used to identify AM fungal species varied significantly between evolutionary groups that correspond to different genera, with the exception of two genera that are genetically indistinguishable from each other. When we clustered groups of similar sequences into operational taxonomic units (OTUs) using five different clustering methods, these patterns of sequence variation caused the number of OTUs to either over- or underestimate the actual number of AM fungal species, depending on the genus. Our results indicate that OTU-based inferences about AM fungal species composition from environmental sequences can be improved if they take these taxonomically structured patterns of sequence variation into account.

Determining threshold values for barcoding fungi: lessons from Cortinarius (Basidiomycota), a highly diverse and widespread ectomycorrhizal genus

Different distance-based threshold selection approaches were used to assess and compare use of the internal transcribed spacer (ITS) region to distinguish among 901 Cortinarius species represented by >3000 collections. Sources of error associated with genetic markers and selection approaches were explored and evaluated using MOTUs from genus and lineage based-alignments. Our study indicates that 1%-2% more species can be distinguished by using the full-length ITS barcode as compared to either the ITS1 or ITS2 regions alone. Optimal threshold values for different picking approaches and genetic marker lengths inferred from a subset of species containing major lineages ranged from 97.0% to 99.5% sequence similarity using clustering optimization and UNITE SH, and from 1% to 2% sequence dissimilarity with CROP. Errors for the optimal cutoff ranged from 0% to 70%, and these can be reduced to a maximum of 22% when excluding species lacking a barcode gap. A threshold value of 99% is suitable for distinguishing species in the majority of lineages in the genus using the entire ITS region but only 90% of the species could be identified using just the ITS1 or ITS2 region. Prior identification of species, lacking barcode gaps and their subsequent separate analyses, maximized the accuracy of threshold approaches.