Tor Carlsen

Fungal community analysis by high-throughput sequencing of amplified markers – a user's guide

Novel high-throughput sequencing methods outperform earlier approaches in terms of resolution and magnitude. They enable identification and relative quantification of community members and offer new insights into fungal community ecology. These methods are currently taking over as the primary tool to assess fungal communities of plant-associated endophytes, pathogens, and mycorrhizal symbionts, as well as free-living saprotrophs. Taking advantage of the collective experience of six research groups, we here review the different stages involved in fungal community analysis, from field sampling via laboratory procedures to bioinformatics and data interpretation. We discuss potential pitfalls, alternatives, and solutions. Highlighted topics are challenges involved in: obtaining representative DNA/RNA samples and replicates that encompass the targeted variation in community composition, selection of marker regions and primers, options for amplification and multiplexing, handling of sequencing errors, and taxonomic identification. Without awareness of methodological biases, limitations of markers, and bioinformatics challenges, large-scale sequencing projects risk yielding artificial results and misleading conclusions.

Changes in the root-associated fungal communities along a primary succession gradient analysed by 454 pyrosequencing

We investigated changes in the root‐associated fungal communities associated with the ectomycorrhizal herb Bistorta vivipara along a primary succession gradient using 454 amplicon sequencing. Our main objective was to assess the degree of variation in fungal richness and community composition as vegetation cover increases along the chronosequence. Sixty root systems of B. vivipara were sampled in vegetation zones delimited by dated moraines in front of a retreating glacier in Norway. We extracted DNA from rinsed root systems, amplified the ITS1 region using fungal‐specific primers and analysed the amplicons using 454 sequencing. Between 437 and 5063 sequences were obtained from each root system. Clustering analyses using a 98.5% sequence similarity cut‐off yielded a total of 470 operational taxonomic units (OTUs), excluding singletons. Between eight and 41 fungal OTUs were detected within each root system. Already in the first stage of succession, a high fungal diversity was present in the B. vivipara root systems. Total number of OTUs increased significantly along the gradient towards climax vegetation, but the average number of OTUs per root system stayed unchanged. There was a high patchiness in distribution of fungal OTUs across root systems, indicating that stochastic processes to a large extent structure the fungal communities. However, time since deglaciation had impact on the fungal community structure, as a systematic shift in the community composition was observed along the chronosequence. Ectomycorrhizal basidiomycetes were the dominant fungi in the roots of B. vivipara, when it comes to both number of OTUs and number of sequences.

High diversity of root associated fungi in both alpine and arctic Dryas octopetala.

BackgroundDryas octopetala is a widespread dwarf shrub in alpine and arctic regions that forms ectomycorrhizal (ECM) symbiotic relationships with fungi. In this study we investigated the fungal communities associated with roots of D. octopetala in alpine sites in Norway and in the High Arctic on Svalbard, where we aimed to reveal whether the fungal diversity and species composition varied across the Alpine and Arctic regions. The internal transcribed spacer (ITS) region of nuclear ribosomal DNA was used to identify the fungal communities from bulk root samples obtained from 24 plants.ResultsA total of 137 operational taxonomic units (OTUs) were detected (using 97% similarity cut off during sequence clustering) and well-known ECM genera such as Cenococcum, Cortinarius, Hebeloma, Inocybe and Tomentella occurred frequently. There was no decrease in fungal diversity with increasing latitude. The overall spatial heterogeneity was high, but a weak geographical structuring of the composition of OTUs in the root systems was observed. Calculated species accumulation curves did not level off.ConclusionsThis study indicates that the diversity of fungi associated with D. octopetala does not decrease in high latitude arctic regions, which contrasts observations made in a wide spectrum of other organism groups. A high degree of patchiness was observed across root systems, but the fungal communities were nevertheless weakly spatially structured. Non-asymptotical species accumulation curves and the occurrence of a high number of singletons indicated that only a small fraction of the fungal diversity was detected.

Towards standardization of the description and publication of next‐generation sequencing datasets of fungal communities

Bonfante P, Genre A. 2008. Plants and arbuscular mycorrhizal fungi: an evolutionary-developmental perspective. Trends in Plant Science 13: 492–498. Bonfante P, Selosse MA. 2010. A glimpse into the past of land plants and of their mycorrhizal affairs: from fossils to evo-devo. New Phytologist 186: 267–270. Boullard B. 1979. Considerations sur la symbiose fongique chez les pteridophytes. Syllogeus 19: 1–58. Brachmann A, Parniske M. 2006. The most widespread symbiosis on Earth. PLoS Biology 4: e239. Brundrett MC. 2002. Coevolution of roots and mycorrhizas of land plants. New Phytologist 154: 275–304. Brundrett MC. 2004. Diversity and classification of mycorrhizal associations. Biological Reviews 78: 473–495. Koltai H, Kapulnik Y, Eds. 2010. Arbuscular mycorrhizas: physiology and function, 2nd edn. Dordrecht, the Netherlands: Springer. Parniske M. 2008. Arbuscular mycorrhiza: the mother of plant root endosymbioses. Nature Reviews: Microbiology 6: 763–775. Phillips TL, DiMichele WA. 1992. Comparative ecology and life-history biology of arborescent lycopsids in Late Carboniferous swamps of Euramerica. Annals of the Missouri Botanical Garden 79: 560–588. Pressel S, Bidartondo MI, Ligrone R, Duckett JG. 2010. Fungal symbioses in bryophytes: new insights in the twenty first century. Phytotaxa 9: 238–253. Rothwell GW, Erwin DM. 1985. The rhizomorphic apex of Paurodendron; implications for homologies among the rooting organs of Lycopsida. American Journal of Botany 72: 86–98. Stewart WN. 1947. A comparative study of stigmarian appendages and Isoetes roots. American Journal of Botany 34: 315–324. Strullu-Derrien C, Rioult JP, Strullu DG. 2009. Mycorrhizas in Upper Carboniferous Radiculites-type cordaitalean rootlets. New Phytologist 182: 561–564. Strullu-Derrien C, Strullu DG. 2007. Mycorrhization of fossil and living plants. Comptes Rendus Palevol 6: 483–494. Stubblefield SP, Rothwell GW. 1981. Embryology and reproductive biology of Bothrodendrostrobus mundus (Lycopsida). American Journal of Botany 68: 625–634. Sudová R, Rydlová J, Ctvrtlı́ková M, Havránek P, Adamec L 2011. The incidence of arbuscular mycorrhiza in two submerged Isoëtes species. Aquatic Botany 94: 183–187. Taber RA, Trappe JM. 1982. Vesicular-arbuscular mycorrhiza in rhizomes, scale-like leaves, roots, and xylem of ginger. Mycologia 74: 156–161. Taylor TN, Taylor EL, Krings M. 2009. Paleobotany. The biology and evolution of fossil plants, 2nd edn. New York, NY, USA: Elsevier ⁄ Academic Press. Wagner CA, Taylor TN. 1981. Evidence for endomycorrhizae in Pennsylvanian age plant fossils. Science 212: 562–563. Winther JL, Friedman WE. 2008. Arbuscular mycorrhizal associations in Lycopodiaceae. New Phytologist 177: 790–801.

CLOTU: An online pipeline for processing and clustering of 454 amplicon reads into OTUs followed by taxonomic annotation

BackgroundThe implementation of high throughput sequencing for exploring biodiversity poses high demands on bioinformatics applications for automated data processing. Here we introduce CLOTU, an online and open access pipeline for processing 454 amplicon reads. CLOTU has been constructed to be highly user-friendly and flexible, since different types of analyses are needed for different datasets.ResultsIn CLOTU, the user can filter out low quality sequences, trim tags, primers, adaptors, perform clustering of sequence reads, and run BLAST against NCBInr or a customized database in a high performance computing environment. The resulting data may be browsed in a user-friendly manner and easily forwarded to downstream analyses. Although CLOTU is specifically designed for analyzing 454 amplicon reads, other types of DNA sequence data can also be processed. A fungal ITS sequence dataset generated by 454 sequencing of environmental samples is used to demonstrate the utility of CLOTU.ConclusionsCLOTU is a flexible and easy to use bioinformatics pipeline that includes different options for filtering, trimming, clustering and taxonomic annotation of high throughput sequence reads. Some of these options are not included in comparable pipelines. CLOTU is implemented in a Linux computer cluster and is freely accessible to academic users through the Bioportal web-based bioinformatics service (http://www.bioportal.uio.no).

High consistency between replicate 454 pyrosequencing analyses of ectomycorrhizal plant root samples.

In this methodological study, we compare 454 sequencing and a conventional cloning and Sanger sequencing approach in their ability to characterize fungal communities PCR amplified from four root systems of the ectomycorrhizal plant Bistorta vivipara. To examine variation introduced by stochastic processes during the laboratory work, we replicated all analyses using two independently obtained DNA extractions from the same root systems. The ITS1 region was used as DNA barcode and the sequences were clustered into OTUs as proxies for species using single linkage clustering (BLASTClust) and 97% sequence similarity cut-off. A relatively low overlap in fungal OTUs was observed between the 454 and the clone library datasets — even among the most abundant OTUs. In a non-metric multidimensional scaling analysis, the samples grouped more according to methodology compared to plant. Some OTUs frequently detected by 454, most notably those OTUs with taxonomic affinity to Glomales, were not detected in the Sanger dataset. Likewise, a few OTUs, including Cenococcum sp., only appeared in the clone libraries. Surprisingly, we observed a significant relationship between GC/AT content of the OTUs and their proportional abundances in the 454 versus the clone library datasets. Reassuringly, a very good consistency in OTU recovery was observed between replicate runs of both sequencing methods. This indicates that stochastic processes had little impact when applying the same sequencing technique on replicate samples.

Employing 454 amplicon pyrosequencing to reveal intragenomic divergence in the internal transcribed spacer rDNA region in fungi

The rDNA internal transcribed spacer (ITS) region has been accepted as a DNA barcoding marker for fungi and is widely used in phylogenetic studies; however, intragenomic ITS variability has been observed in a broad range of taxa, including prokaryotes, plants, animals, and fungi, and this variability has the potential to inflate species richness estimates in molecular investigations of environmental samples. In this study 454 amplicon pyrosequencing of the ITS1 region was applied to 99 phylogenetically diverse axenic single-spore cultures of fungi (Dikarya: Ascomycota and Basidiomycota) to investigate levels of intragenomic variation. Three species (one Basidiomycota and two Ascomycota), in addition to a positive control species known to contain ITS paralogs, displayed levels of molecular variation indicative of intragenomic variation; taxon inflation due to presumed intragenomic variation was ≈9%. Intragenomic variability in the ITS region appears to be widespread but relatively rare in fungi (≈3–5% of species investigated in this study), suggesting this problem may have minor impacts on species richness estimates relative to PCR and/or pyrosequencing errors. Our results indicate that 454 amplicon pyrosequencing represents a powerful tool for investigating levels of ITS intragenomic variability across taxa, which may be valuable for better understanding the fundamental mechanisms underlying concerted evolution of repetitive DNA regions.

Arctic root-associated fungal community composition reflects environmental filtering.

There is growing evidence that root‐associated fungi have important roles in Arctic ecosystems. Here, we assess the diversity of fungal communities associated with roots of the ectomycorrhizal perennial herb Bistorta vivipara on the Arctic archipelago of Svalbard and investigate whether spatial separation and bioclimatic variation are important structuring factors of fungal community composition. We sampled 160 plants of B. vivipara from 32 localities across Svalbard. DNA was extracted from entire root systems, and 454 pyrosequencing of ITS1 amplicons was used to profile the fungal communities. The fungal communities were predominantly composed of Basidiomycota (55% of reads) and Ascomycota (35%), with the orders Thelephorales (24%), Agaricales (13.8%), Pezizales (12.6%) and Sebacinales (11.3%) accounting for most of the reads. Plants from the same site or region had more similar fungal communities to one another than plants from other sites or regions, and sites clustered together along a weak latitudinal gradient. Furthermore, a decrease in per‐plant OTU richness with increasing latitude was observed. However, no statistically significant spatial autocorrelation between sites was detected, suggesting that environmental filtering, not dispersal limitation, causes the observed patterns. Our analyses suggest that while latitudinal patterns in community composition and richness might reflect bioclimatic influences at global spatial scales, at the smaller spatial scale of the Svalbard archipelago, these changes more likely reflect varied bedrock composition and associated edaphic factors. The need for further studies focusing on identifying those specific bioclimatic and edaphic factors structuring root‐associated fungal community composition at both global and local scales is emphasized.

Biogeography and Phylogeny of Cardamine (Brassicaceae)1

Abstract The biogeography and phylogeny of Cardamine L. were inferred based on sequences of the nuclear ribosomal ITS regions and the plastid trnL intron and trnL-F spacer regions. This genus is one of the largest and polyploid-rich genera of the Brassicaceae and has its center of diversity in Eurasia. Species were included from all populated continents, representing all sections except two monotypic ones. The results support a hypothesis of recent and rapid speciation in the genus. The traditional sectional classification was not supported. We found evidence for several extremely long-distance dispersal events. Colonization of the Southern Hemisphere and the Arctic has occurred repeatedly; we identified at least three phylogenetically distinct Arctic lineages, two distinct Oceanian lineages, and four distinct South American lineages. Polyploidization has occurred independently many times during the evolution of Cardamine. Recent divergence combined with widespread polyploidization offer an explanation for the complex taxonomy of the genus.

Different bacterial communities in ectomycorrhizae and surrounding soil

Several eukaryotic symbioses have shown to host a rich diversity of prokaryotes that interact with their hosts. Here, we study bacterial communities associated with ectomycorrhizal root systems of Bistorta vivipara compared to bacterial communities in bulk soil using pyrosequencing of 16S rRNA amplicons. A high richness of Operational Taxonomic Units (OTUs) was found in plant roots (3,571 OTUs) and surrounding soil (3,476 OTUs). The community composition differed markedly between these two environments. Actinobacteria, Armatimonadetes, Chloroflexi and OTUs unclassified at phylum level were significantly more abundant in plant roots than in soil. A large proportion of the OTUs, especially those in plant roots, presented low similarity to Sanger 16S rRNA reference sequences, suggesting novel bacterial diversity in ectomycorrhizae. Furthermore, the bacterial communities of the plant roots were spatially structured up to a distance of 60 cm, which may be explained by bacteria using fungal hyphae as a transport vector. The analyzed ectomycorrhizae presents a distinct microbiome, which likely influence the functioning of the plant-fungus symbiosis.