Increasing ecological inference from high throughput sequencing of fungi in the environment through a tagging approach

Year-Long Microbial Succession on Microplastics in Wastewater: Chaotic Dynamics Outweigh Preferential Growth

Microplastics are a globally-ubiquitous aquatic pollutant and have been heavily studied over the last decade. Of particular interest are the interactions between microplastics and microorganisms, especially the pursuit to discover a plastic-specific biome, the so-called plastisphere. To follow this up, a year-long microcosm experimental setup was deployed to expose five different microplastic types (and silica beads control) to activated aerobic wastewater in controlled conditions, with microbial communities being measured four times over the course of the year using 16S rDNA (bacterial) and ITS (fungal) amplicon sequencing. The biofilm community shows no evidence of a specific plastisphere, even after a year of incubation. Indeed, the microbial communities (particularly bacterial) show a clear trend of increasing dissimilarity between plastic types as time increases. Despite little evidence for a plastic-specific community, there was a slight grouping observed for polyolefins (PE and PP) in 6–12-month biofilms. Additionally, an OTU assigned to the genus Devosia was identified on many plastics, increasing over time while showing no growth on silicate (natural particle) controls, suggesting this could be either a slow-growing plastic-specific taxon or a symbiont to such. Both substrate-associated findings were only possible to observe in samples incubated for 6–12 months, which highlights the importance of studying long-term microbial community dynamics on plastic surfaces.

Plant Identity Influences Foliar Fungal Symbionts More Than Elevation in the Colorado Rocky Mountains

Despite colonizing nearly every plant on Earth, foliar fungal symbionts have received little attention in studies on the biogeography of host-associated microbes. Evidence from regional scale studies suggests that foliar fungal symbiont distributions are influenced both by plant hosts and environmental variation in climate and soil resources. However, previous surveys have focused on either one plant host across an environmental gradient or one gradient and multiple plant hosts, making it difficult to disentangle the influence of host identity from the influence of the environment on foliar endophyte communities. We used a culture-based approach to survey fungal symbiont composition in the leaves of nine C₃ grass species along replicated elevation gradients in grasslands of the Colorado Rocky Mountains. In these ecosystems, the taxonomic richness and composition of foliar fungal symbionts were mostly structured by the taxonomic identity of the plant host rather than by variation in climate. Plant traits related to size (height and leaf length) were the best predictors of foliar fungal symbiont composition and diversity, and composition did not vary predictably with plant evolutionary history. The largest plants had the most diverse and distinctive fungal communities. These results suggest that across the ~ 300 m elevation range that we sampled, foliar fungal symbionts may indirectly experience climate change by tracking the shifting distributions of plant hosts rather than tracking climate directly.

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.

Three New Eremophilanes from a Ligularia Hybrid Collected in China

Three new eremophilanes were isolated from a sample of Ligularia collected in Yunnan province of China and their structures were determined. The chemical constituents including eleven eremophilanes (major constituents: 6-acyloxyfuranoeremophilan-10β-ols) were similar to those of a hybrid of L. cyathiceps and L. subspicata, however, the base sequence of the ITS1-5.8S-ITS2 region was typical of L. cyathiceps, indicating backcrossing.

Indexed PCR Primers Induce Template-Specific Bias in Large-Scale DNA Sequencing Studies

Massively parallel sequencing is rapidly emerging as an efficient way to quantify biodiversity at all levels, from genetic variation and expression to ecological community assemblage. However, the number of reads produced per sequencing run far exceeds the number required per sample for many applications, compelling researchers to sequence multiple samples per run in order to maximize efficiency. For studies that include a PCR step, this can be accomplished using primers that include an index sequence allowing sample origin to be determined after sequencing. The use of indexed primers assumes they behave no differently than standard primers; however, we found that indexed primers cause substantial template sequence-specific bias, resulting in radically different profiles of the same environmental sample. Likely the outcome of differential amplification efficiency due to primer-template mismatch, two indexed primer sets spuriously change the inferred sequence abundance from the same DNA extraction by up to 77.1%. We demonstrate that a double PCR approach alleviates these effects in applications where indexed primers are necessary.

Metabarcoding of fungal communities associated with bark beetles

Many species of fungi are closely allied with bark beetles, including many tree pathogens, but their species richness and patterns of distribution remain largely unknown. We established a protocol for metabarcoding of fungal communities directly from total genomic DNA extracted from individual beetles, showing that the ITS3/4 primer pair selectively amplifies the fungal ITS. Using three specimens of bark beetle from different species, we assess the fungal diversity associated with these specimens and the repeatability of these estimates in PCRs conducted with different primer tags. The combined replicates produced 727 fungal Operational Taxonomic Units (OTUs) for the specimen of Hylastes ater, 435 OTUs for Tomicus piniperda, and 294 OTUs for Trypodendron lineatum, while individual PCR reactions produced on average only 229, 54, and 31 OTUs for the three specimens, respectively. Yet, communities from PCR replicates were very similar in pairwise comparisons, in particular when considering species abundance, but differed greatly among the three beetle specimens. Different primer tags or the inclusion of amplicons in separate libraries did not impact the species composition. The ITS2 sequences were identified with the Lowest Common Ancestor approach and correspond to diverse lineages of fungi, including Ophiostomaceae and Leotiomycetes widely found to be tree pathogens. We conclude that Illumina MiSeq metabarcoding reliably captures fungal diversity associated with bark beetles, although numerous PCR replicates are recommended for an exhaustive sample. Direct PCR from beetle DNA extractions provides a rapid method for future surveys of fungal species diversity and their associations with bark beetles and environmental variables.

Rich and cold: diversity, distribution and drivers of fungal communities in patterned-ground ecosystems of the North American Arctic

Fungi are abundant and functionally important in the Arctic, yet comprehensive studies of their diversity in relation to geography and environment are not available. We sampled soils in paired plots along the North American Arctic Transect (NAAT), which spans all five bioclimatic subzones of the Arctic. Each pair of plots contrasted relatively bare, cryoturbated patterned-ground features (PGFs) and adjacent vegetated between patterned-ground features (bPGFs). Fungal communities were analysed via sequencing of 7834 ITS-LSU clones. We recorded 1834 OTUs - nearly half the fungal richness previously reported for the entire Arctic. These OTUs spanned eight phyla, 24 classes, 75 orders and 120 families, but were dominated by Ascomycota, with one-fifth belonging to lichens. Species richness did not decline with increasing latitude, although there was a decline in mycorrhizal taxa that was offset by an increase in lichen taxa. The dominant OTUs were widespread even beyond the Arctic, demonstrating no dispersal limitation. Yet fungal communities were distinct in each subzone and were correlated with soil pH, climate and vegetation. Communities in subzone E were distinct from the other subzones, but similar to those of the boreal forest. Fungal communities on disturbed PGFs differed significantly from those of paired stable areas in bPGFs. Indicator species for PGFs included lichens and saprotrophic fungi, while bPGFs were characterized by ectomycorrhizal and pathogenic fungi. Our results suggest that the Arctic does not host a unique mycoflora, while Arctic fungi are highly sensitive to climate and vegetation, with potential to migrate rapidly as global change unfolds.

An open source chimera checker for the fungal ITS region

The internal transcribed spacer (ITS) region of the nuclear ribosomal repeat unit holds a central position in the pursuit of the taxonomic affiliation of fungi recovered through environmental sampling. Newly generated fungal ITS sequences are typically compared against the International Nucleotide Sequence Databases for a species or genus name using the sequence similarity software suite blast. Such searches are not without complications however, and one of them is the presence of chimeric entries among the query or reference sequences. Chimeras are artificial sequences, generated unintentionally during the polymerase chain reaction step, that feature sequence data from two (or possibly more) distinct species. Available software solutions for chimera control do not readily target the fungal ITS region, but the present study introduces a blast-based open source software package (available at http://www.emerencia.org/chimerachecker.html) to examine newly generated fungal ITS sequences for the presence of potentially chimeric elements in batch mode. We used the software package on a random set of 12 300 environmental fungal ITS sequences in the public sequence databases and found 1.5% of the entries to be chimeric at the ordinal level after manual verification of the results. The proportion of chimeras in the sequence databases can be hypothesized to increase as emerging sequencing technologies drawing from pooled DNA samples are becoming important tools in molecular ecology research.

Dandruff is associated with disequilibrium in the proportion of the major bacterial and fungal populations colonizing the scalp

The bacterial and fungal communities associated with dandruff were investigated using culture-independent methodologies in the French subjects. The major bacterial and fungal species inhabiting the scalp subject’s were identified by cloning and sequencing of the conserved ribosomal unit regions (16S for bacterial and 28S-ITS for fungal) and were further quantified by quantitative PCR. The two main bacterial species found on the scalp surface were Propionibacterium acnes and Staphylococcus epidermidis, while Malassezia restricta was the main fungal inhabitant. Dandruff was correlated with a higher incidence of M. restricta and S. epidermidis and a lower incidence of P. acnes compared to the control population (p<0.05). These results suggested for the first time using molecular methods, that dandruff is linked to the balance between bacteria and fungi of the host scalp surface.

Evaluation of the authenticity of a highly novel environmental sequence from boreal forest soil using ribosomal RNA secondary structure modeling

The number of sequences from both formally described taxa and uncultured environmental DNA deposited in the International Nucleotide Sequence Databases has increased substantially over the last two decades. Although the majority of these sequences represent authentic gene copies, there is evidence of DNA artifacts in these databases as well. These include lab artifacts, such as PCR chimeras, and biological artifacts such as pseudogenes or other paralogous sequences. Sequences that fall in basal positions in phylogenetic trees and appear distant from known sequences are particularly suspect. Phylogenetic analyses suggest that a novel sequence type (NS1) found in two boreal forest soil clone libraries belongs to the fungal kingdom but does not fall unambiguously within any known phylum. We have evaluated this sequence type using an array of secondary-structure analyses. To our knowledge, such analyses have never been used on environmental ribosomal sequences. Ribosomal secondary structure was modeled for four rRNA loci (ITS1, 5.8S, ITS2, 5' LSU). These models were analyzed for the presence of conserved domains, conserved nucleotide motifs, and compensatory base changes. Minimal free energy (MFE) foldings and GC contents of sequences representing the major fungal clades, as well as NS1, were also compared. NS1 displays secondary rRNA structures consistent with other fungi and many, but not all, conserved nucleotide motifs found across eukaryotes. However, our analyses show that many other authentic sequences from basal fungi lack more of these conserved motifs than does NS1. Together our findings suggest that NS1 represents an authentic gene copy. The methods described here can be used on any rRNA-coding sequence, not just environmental fungal sequences. As new-generation sequencing methods that yield shorter sequences become more widely implemented, methods that evaluate sequence authenticity should also be more widely implemented. For fungi, the adjacent 5.8S and ITS2 loci should be prioritized. This region is not only suited to distinguishing between closely related species, but it is also more informative in terms of expected secondary structure.

Mycorrhizal response to experimental pH and P manipulation in acidic hardwood forests

Many temperate forests of the Northeastern United States and Europe have received significant anthropogenic acid and nitrogen (N) deposition over the last century. Although temperate hardwood forests are generally thought to be N-limited, anthropogenic deposition increases the possibility of phosphorus (P) limiting productivity in these forest ecosystems. Moreover, inorganic P availability is largely controlled by soil pH and biogeochemical theory suggests that forests with acidic soils (i.e., Collapse

Key Words Collapse

MESH Headings

• Biomass
• Ecosystem
• Hydrogen-Ion Concentration
• Mycorrhizae/metabolism
• Ohio
• Phosphorus/metabolism
• Plant Roots/metabolism
• Plant Roots/microbiology
• Soil
• Soil Microbiology
• Trees/metabolism
• Trees/microbiology

Collapse

Grants Collapse

Affiliation(s)

Laurel A Kluber Department of Biology, Case Western Reserve University, Cleveland, OH, USA.
Sarah R Carrino-Kyker
Kaitlin P Coyle
Jared L DeForest
Charlotte R Hewins
Alanna N Shaw
Kurt A Smemo
David J Burke

Collapse

Back to basics: an evaluation of NaOH and alternative rapid DNA extraction protocols for DNA barcoding, genotyping, and disease diagnostics from fungal and oomycete samples

The ubiquity, high diversity and often-cryptic manifestations of fungi and oomycetes frequently necessitate molecular tools for detecting and identifying them in the environment. In applications including DNA barcoding, pathogen detection from plant samples, and genotyping for population genetics and epidemiology, rapid and dependable DNA extraction methods scalable from one to hundreds of samples are desirable. We evaluated several rapid extraction methods (NaOH, Rapid one-step extraction (ROSE), Chelex 100, proteinase K) for their ability to obtain DNA of quantity and quality suitable for the following applications: PCR amplification of the multicopy barcoding locus ITS1/5.8S/ITS2 from various fungal cultures and sporocarps; single-copy microsatellite amplification from cultures of the phytopathogenic oomycete Phytophthora ramorum; probe-based P. ramorum detection from leaves. Several methods were effective for most of the applications, with NaOH extraction favored in terms of success rate, cost, speed and simplicity. Frozen dilutions of ROSE and NaOH extracts maintained PCR viability for over 32 months. DNA from rapid extractions performed poorly compared to CTAB/phenol-chloroform extracts for TaqMan diagnostics from tanoak leaves, suggesting that incomplete removal of PCR inhibitors is an issue for sensitive diagnostic procedures, especially from plants with recalcitrant leaf chemistry. NaOH extracts exhibited lower yield and size than CTAB/phenol-chloroform extracts; however, NaOH extraction facilitated obtaining clean sequence data from sporocarps contaminated by other fungi, perhaps due to dilution resulting from low DNA yield. We conclude that conventional extractions are often unnecessary for routine DNA sequencing or genotyping of fungi and oomycetes, and recommend simpler strategies where source materials and intended applications warrant such use.

How well do ITS rDNA sequences differentiate species of true morels (Morchella)?

Arguably more mycophiles hunt true morels (Morchella) during their brief fruiting season each spring in the northern hemisphere than any other wild edible fungus. Concerns about overharvesting by individual collectors and commercial enterprises make it essential that science-based management practices and conservation policies are developed to ensure the sustainability of commercial harvests and to protect and preserve morel species diversity. Therefore, the primary objectives of the present study were to: (i) investigate the utility of the ITS rDNA locus for identifying Morchella species, using phylogenetic species previously inferred from multilocus DNA sequence data as a reference; and (ii) clarify insufficiently identified sequences and determine whether the named sequences in GenBank were identified correctly. To this end, we generated 553 Morchella ITS rDNA sequences and downloaded 312 additional ones generated by other researchers from GenBank using emerencia and analyzed them phylogenetically. Three major findings emerged: (i) ITS rDNA sequences were useful in identifying 48/62 (77.4%) of the known phylospecies; however, they failed to identify 12 of the 22 species within the species-rich Elata Subclade and two closely related species in the Esculenta Clade; (ii) at least 66% of the named Morchella sequences in GenBank are misidentified; and (iii) ITS rDNA sequences of up to six putatively novel Morchella species were represented in GenBank. Recognizing the need for a dedicated Web-accessible reference database to facilitate the rapid identification of known and novel species, we constructed Morchella MLST (http://www.cbs.knaw.nl/morchella/), which can be queried with ITS rDNA sequences and those of the four other genes used in our prior multilocus molecular systematic studies of this charismatic genus.

Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi

Six DNA regions were evaluated as potential DNA barcodes for Fungi, the second largest kingdom of eukaryotic life, by a multinational, multilaboratory consortium. The region of the mitochondrial cytochrome c oxidase subunit 1 used as the animal barcode was excluded as a potential marker, because it is difficult to amplify in fungi, often includes large introns, and can be insufficiently variable. Three subunits from the nuclear ribosomal RNA cistron were compared together with regions of three representative protein-coding genes (largest subunit of RNA polymerase II, second largest subunit of RNA polymerase II, and minichromosome maintenance protein). Although the protein-coding gene regions often had a higher percent of correct identification compared with ribosomal markers, low PCR amplification and sequencing success eliminated them as candidates for a universal fungal barcode. Among the regions of the ribosomal cistron, the internal transcribed spacer (ITS) region has the highest probability of successful identification for the broadest range of fungi, with the most clearly defined barcode gap between inter- and intraspecific variation. The nuclear ribosomal large subunit, a popular phylogenetic marker in certain groups, had superior species resolution in some taxonomic groups, such as the early diverging lineages and the ascomycete yeasts, but was otherwise slightly inferior to the ITS. The nuclear ribosomal small subunit has poor species-level resolution in fungi. ITS will be formally proposed for adoption as the primary fungal barcode marker to the Consortium for the Barcode of Life, with the possibility that supplementary barcodes may be developed for particular narrowly circumscribed taxonomic groups.

Archaeorhizomycetes: unearthing an ancient class of ubiquitous soil fungi

Estimates suggest that only one-tenth of the true fungal diversity has been described. Among numerous fungal lineages known only from environmental DNA sequences, Soil Clone Group 1 is the most ubiquitous. These globally distributed fungi may dominate below-ground fungal communities, but their placement in the fungal tree of life has been uncertain. Here, we report cultures of this group and describe the class, Archaeorhizomycetes, phylogenetically placed within subphylum Taphrinomycotina in the Ascomycota. Archaeorhizomycetes comprises hundreds of cryptically reproducing filamentous species that do not form recognizable mycorrhizal structures and have saprotrophic potential, yet are omnipresent in roots and rhizosphere soil and show ecosystem and host root habitat specificity.

A bioinformatics pipeline for sequence-based analyses of fungal biodiversity

The internal transcribed spacer (ITS) is the locus of choice with which to characterize fungal diversity in environmental samples. However, methods to analyze ITS datasets have lagged behind the capacity to generate large amounts of sequence information. Here, we describe our bioinformatics pipeline to process large fungal ITS sequence datasets, from raw chromatograms to a spreadsheet of operational taxonomic unit (OTU) abundances across samples. Steps include assembling of reads originating from one clone, identifying primer "barcodes" or "tags," trimming vectors and primers, marking low-quality base calls and removing low-quality sequences, orienting sequences, extracting the ITS region from longer amplicons, and grouping sequences into OTUs. We expect that the principles and tools presented here are relevant to datasets arising from ever-evolving new technologies.

Phylogenetic and ecological analyses of soil and sporocarp DNA sequences reveal high diversity and strong habitat partitioning in the boreal ectomycorrhizal genus Russula (Russulales; Basidiomycota)

*Although critical for the functioning of ecosystems, fungi are poorly known in high-latitude regions. Here, we provide the first genetic diversity assessment of one of the most diverse and abundant ectomycorrhizal genera in Alaska: Russula. *We analyzed internal transcribed spacer rDNA sequences from sporocarps and soil samples using phylogenetic methods, operational taxonomic unit (OTU) delimitations and ordinations to compare species composition in various types of boreal forest. *The genus Russula is highly diverse in Alaska, with at least 42 nonsingleton OTUs (soil) and 50 phylogroups (soil + sporocarp). Russula taxa showed strong habitat preference to one of the two major forest types in the sampled regions (black spruce and birch-aspen-white spruce), and some preference for soil horizon. *Our results show that the vast majority of Russula species are present in the soil samples, although some additional taxa are expected to be found with extended sampling. OTU diversity in black spruce forests was only one-third of the diversity observed in mixed upland forests. Our findings suggest that some of the diversity is niche based, especially along host and successional axes, because most OTUs predictably occurred in specific habitats, regardless of geographical location.

Stable isotope fingerprinting: a novel method for identifying plant, fungal, or bacterial origins of amino acids

Amino acids play an important role in ecology as essential nutrients for animals and as currencies in symbiotic associations. Here we present a new approach to tracing the origins of amino acids by identifying unique patterns of carbon isotope signatures generated by amino acid synthesis in plants, fungi, and bacteria ("13C fingerprints"). We measured amino acid delta 13C from 10 C3 plants, 13 fungi, and 10 bacteria collected and isolated from a boreal forest in interior Alaska, USA, using gas chromatography-combustion-isotope ratio mass spectrometry (GC-C-IRMS). Microorganisms were cultured under amino-acid-free conditions and identified based on DNA sequences. Bacteria, fungi, and plants generated consistent, unique 13C fingerprints based on the more complex amino acids (five or more biosynthetic steps) that are classified as essential for animals. Linear discriminant analysis classified all samples correctly with >99% certainty and correctly classified nearly all insect samples from a previous study by diet. Our results suggest that 13C fingerprints of amino acids could provide a powerful in situ assay of the biosynthetic sources of amino acids and a potential new tool for understanding nutritional linkages in food webs.

Top-down systems biology modeling of host metabotype-microbiome associations in obese rodents

Covariation in the structural composition of the gut microbiome and the spectroscopically derived metabolic phenotype (metabotype) of a rodent model for obesity were investigated using a range of multivariate statistical tools. Urine and plasma samples from three strains of 10-week-old male Zucker rats (obese (fa/fa, n=8), lean (fa/-, n=8) and lean (-/-, n=8)) were characterized via high-resolution 1H NMR spectroscopy, and in parallel, the fecal microbial composition was investigated using fluorescence in situ hydridization (FISH) and denaturing gradient gel electrophoresis (DGGE) methods. All three Zucker strains had different relative abundances of the dominant members of their intestinal microbiota (FISH), with the novel observation of a Halomonas and a Sphingomonas species being present in the (fa/fa) obese strain on the basis of DGGE data. The two functionally and phenotypically normal Zucker strains (fa/- and -/-) were readily distinguished from the (fa/fa) obese rats on the basis of their metabotypes with relatively lower urinary hippurate and creatinine, relatively higher levels of urinary isoleucine, leucine and acetate and higher plasma LDL and VLDL levels typifying the (fa/fa) obese strain. Collectively, these data suggest a conditional host genetic involvement in selection of the microbial species in each host strain, and that both lean and obese animals could have specific metabolic phenotypes that are linked to their individual microbiomes.

Molecular phylogenetic biodiversity assessment of arctic and boreal ectomycorrhizal Lactarius Pers. (Russulales; Basidiomycota) in Alaska, based on soil and sporocarp DNA

Despite the critical roles fungi play in the functioning of ecosystems, especially as symbionts of plants and recyclers of organic matter, their biodiversity is poorly known in high-latitude regions. In this paper, we discuss the molecular diversity of one of the most diverse and abundant groups of ectomycorrhizal fungi: the genus Lactarius Pers. We analysed internal transcribed spacer rDNA sequences from both curated sporocarp collections and soil polymerase chain reaction clone libraries sampled in the arctic tundra and boreal forests of Alaska. Our genetic diversity assessment, based on various phylogenetic methods and operational taxonomic unit (OTU) delimitations, suggests that the genus Lactarius is diverse in Alaska, with at least 43 putative phylogroups, and 24 and 38 distinct OTUs based on 95% and 97% internal transcribed spacer sequence similarity, respectively. Some OTUs were identified to known species, while others were novel, previously unsequenced groups. Non-asymptotic species accumulation curves, the disparity between observed and estimated richness, and the high number of singleton OTUs indicated that many Lactarius species remain to be found and identified in Alaska. Many Lactarius taxa show strong habitat preference to one of the three major vegetation types in the sampled regions (arctic tundra, black spruce forests, and mixed birch-aspen-white spruce forests), as supported by statistical tests of UniFrac distances and principal coordinates analyses (PCoA). Together, our data robustly demonstrate great diversity and nonrandom ecological partitioning in an important boreal ectomycorrhizal genus within a relatively small geographical region. The observed diversity of Lactarius was much higher in either type of boreal forest than in the arctic tundra, supporting the widely recognized pattern of decreasing species richness with increasing latitude.