Fungi Sailing the Arctic Ocean: Speciose Communities in North Atlantic Driftwood as Revealed by High-Throughput Amplicon Sequencing.

High amounts of driftwood sail across the oceans and provide habitat for organisms tolerating the rough and saline environment. Fungi have adapted to the extremely cold and saline conditions which driftwood faces in the high north. For the first time, we applied high-throughput sequencing to fungi residing in driftwood to reveal their taxonomic richness, community composition, and ecology in the North Atlantic. Using pyrosequencing of ITS2 amplicons obtained from 49 marine logs, we found 807 fungal operational taxonomic units (OTUs) based on clustering at 97 % sequence similarity cut-off level. The phylum Ascomycota comprised 74 % of the OTUs and 20 % belonged to Basidiomycota. The richness of basidiomycetes decreased with prolonged submersion in the sea, supporting the general view of ascomycetes being more extremotolerant. However, more than one fourth of the fungal OTUs remained unassigned to any fungal class, emphasising the need for better DNA reference data from the marine habitat. Different fungal communities were detected in coniferous and deciduous logs. Our results highlight that driftwood hosts a considerably higher fungal diversity than currently known. The driftwood fungal community is not a terrestrial relic but a speciose assemblage of fungi adapted to the stressful marine environment and different kinds of wooden substrates found in it.

Primary succession of Bistorta vivipara (L.) Delabre (Polygonaceae) root-associated fungi mirrors plant succession in two glacial chronosequences.

Glacier chronosequences are important sites for primary succession studies and have yielded well-defined primary succession models for plants that identify environmental resistance as an important determinant of the successional trajectory. Whether plant-associated fungal communities follow those same successional trajectories and also respond to environmental resistance is an open question. In this study, 454 amplicon pyrosequencing was used to compare the root-associated fungal communities of the ectomycorrhizal (ECM) herb Bistorta vivipara along two primary succession gradients with different environmental resistance (alpine versus arctic) and different successional trajectories in the vascular plant communities (directional replacement versus directional non-replacement). At both sites, the root-associated fungal communities were dominated by ECM basidiomycetes and community composition shifted with increasing time since deglaciation. However, the fungal community's successional trajectory mirrored the pattern observed in the surrounding plant community at both sites: the alpine site displayed a directional-replacement successional trajectory, and the arctic site displayed a directional-non-replacement successional trajectory. This suggests that, like in plant communities, environmental resistance is key in determining succession patterns in root-associated fungi. The need for further replicated study, including in other host species, is emphasized.

Mycorrhizal fungal communities in coastal sand dunes and heaths investigated by pyrosequencing analyses.

Maritime sand dunes and coastal ericaceous heaths are unstable and dynamic habitats for mycorrhizal fungi. Creeping willow (Salix repens) is an important host plant in these habitats in parts of Europe. In this study, we wanted to assess which mycorrhizal fungi are associated with S. repens in four different coastal vegetation types in Southern Norway, three types from sand dunes and one from heaths. Moreover, we investigated which ecological factors are important for the fungal community structure in these vegetation types. Mycorrhizal fungi on S. repens root samples were identified by 454 pyrosequencing of tag-encoded internal transcribed spacer 1 (ITS1) amplicons. Significantly higher fungal richness was observed in hummock dunes and dune slacks compared to eroded dune vegetation. The compositional variation was mainly accounted for by location (plot) and vegetation type and was significantly correlated to content of carbon, nitrogen and phosphorus in soil. The investigated maritime sand dunes and coastal ericaceous heaths hosted mycorrhizal taxa mainly associated with Helotiales, Sebacinales, Thelephorales and Agaricales.

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.

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.

The devil in the details: interactions between the branch-length prior and likelihood model affect node support and branch lengths in the phylogeny of the Psoraceae.

In popular use of Bayesian phylogenetics, a default branch-length prior is almost universally applied without knowing how a different prior would have affected the outcome. We performed Bayesian and maximum likelihood (ML) inference of phylogeny based on empirical nucleotide sequence data from a family of lichenized ascomycetes, the Psoraceae, the morphological delimitation of which has been controversial. We specifically assessed the influence of the combination of Bayesian branch-length prior and likelihood model on the properties of the Markov chain Monte Carlo tree sample, including node support, branch lengths, and taxon stability. Data included two regions of the mitochondrial ribosomal RNA gene, the internal transcribed spacer region of the nuclear ribosomal RNA gene, and the protein-coding largest subunit of RNA polymerase II. Data partitioning was performed using Bayes' factors, whereas the best-fitting model of each partition was selected using the Bayesian information criterion (BIC). Given the data and model, short Bayesian branch-length priors generate higher numbers of strongly supported nodes as well as short and topologically similar trees sampled from parts of tree space that are largely unexplored by the ML bootstrap. Long branch-length priors generate fewer strongly supported nodes and longer and more dissimilar trees that are sampled mostly from inside the range of tree space sampled by the ML bootstrap. Priors near the ML distribution of branch lengths generate the best marginal likelihood and the highest frequency of "rogue" (unstable) taxa. The branch-length prior was shown to interact with the likelihood model. Trees inferred under complex partitioned models are more affected by the stretching effect of the branch-length prior. Fewer nodes are strongly supported under a complex model given the same branch-length prior. Irrespective of model, internal branches make up a larger proportion of total tree length under the shortest branch-length priors compared with longer priors. Relative effects on branch lengths caused by the branch-length prior can be problematic to downstream phylogenetic comparative methods making use of the branch lengths. Furthermore, given the same branch-length prior, trees are on average more dissimilar under a simple unpartitioned model compared with a more complex partitioned models. The distribution of ML branch lengths was shown to better fit a gamma or Pareto distribution than an exponential one. Model adequacy tests indicate that the best-fitting model selected by the BIC is insufficient for describing data patterns in 5 of 8 partitions. More general substitution models are required to explain the data in three of these partitions, one of which also requires nonstationarity. The two mitochondrial ribosomal RNA gene partitions need heterotachous models. We found no significant correlations between, on the one hand, the amount of ambiguous data or the smallest branch-length distance to another taxon and, on the other hand, the topological stability of individual taxa. Integrating over several exponentially distributed means under the best-fitting model, node support for the family Psoraceae, including Psora, Protoblastenia, and the Micarea sylvicola group, is approximately 0.96. Support for the genus Psora is distinctly lower, but we found no evidence to contradict the current classification.

Habitat Protection Approaches Facilitate Conservation of Overlooked Fungal Diversity - A Case Study From the Norwegian Coastal Heathland System.

European coastal heathlands are distinct ecosystems shaped by land use tradition and they have experienced an 80% area reduction from their historical maximum. These mosaics of mires and wind exposed patches have ericaceous shrub dominated vegetation, and soils within coastal heathlands are characterized by low pH and high levels of recalcitrant debris. Using a culture-based approach with molecular identification of isolates, we characterized root-associated fungal communities of six ericaceous species in eight heathland localities along Norway's western coast. Site-level alpha diversity ranged from 21-38 OTUs, while the total estimated gamma diversity for culturable heathland root fungi was 190-231 OTUs. Most species recovered are previously reported at low abundance in Norway, suggesting the biodiversity in this community is underreported, rather than novel for science. The fungi recovered were primarily Ascomycota, specifically endophytic Phialocephala, and Pezicula, and no host specificity was observed in the communities. The fungal communities exhibited high turnover and low nestedness, both between ericaceous hosts and across heathland sites. We observed no spatial patterns in fungal betadiversity, and this heterogeneity may be a product of the unique historic land use practices at each locality creating a distinct mycofloral "fingerprint". Robust diversity estimates will be key for managing fungal biodiversity in coastal heathlands. Our results indicate that sampling schemes that maximize the number of host plants sampled per site, rather than the number of cultures per plant yield improved alpha diversity estimates. Similarly, gamma diversity estimates are improved by maximizing the total number of localities sampled, rather than increasing the number of plants sampled per locality. We argue that while the current protected status of coastal heathland habitats and restoration efforts have knock-on effects for the conservation of fungal biodiversity, fungi have a vital functional role in the ecosystem and holistic conservation plans that consider fungal biodiversity would be beneficial.

Directed endozoochorous dispersal by scavengers facilitate sexual reproduction in otherwise clonal plants at cadaver sites.

The regeneration niche of many plant species involves spatially and temporally unpredictable disturbances, called recruitment windows of opportunity. However, even species with clear dispersal adaptations such as fleshy berries may not successfully reach such elusive regeneration microsites. Ericaceous, berry-producing species in the northern hemisphere demonstrate this dispersal limitation. They are said to display a reproductive paradox owing to their lack of regeneration in apparently suitable microsites despite considerable investment in producing large quantities of berries.Cadavers generate vegetation-denuded and nutrient-rich disturbances termed cadaver decomposition islands (CDIs). Cadavers attract facultative scavengers with considerable capacity for endozoochorous seed dispersal. We hypothesize that CDIs facilitate recruitment in berry-producing ericaceous species due to endozoochorous dispersal directed toward favorable microsites with low competition.We examined seedling establishment within a permanent, semi-regular 10 × 10 m grid across an ungulate mass die-off on the Hardangervidda plateau in southeastern Norway. Competing models regarding the relative importance of factors governing recruitment were evaluated, specifically cadaver location (elevated seed rain) and microsite conditions (competition).We found that CDIs did facilitate seedling establishment, as cadaver density was the best predictor of seedling distribution. Other important factors governing seedling establishment such as percentage cover of soil and vascular plants alone were inadequate to explain seedling establishment. Synthesis: This study provides a novel understanding of sexual reproduction in species with cryptic generative reproduction. The directed nature of endozoochorous dispersal combined with long-distance dispersal abilities of medium to large vertebrate scavengers toward cadavers allows plants to exploit the advantageous but ephemeral resource provided by CDIs.

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

BACKGROUND: The 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. RESULTS: In 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. CONCLUSIONS: CLOTU 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).

Fine-scale diversity patterns in belowground microbial communities are consistent across kingdoms.

The belowground environment is heterogeneous and complex at fine spatial scales. Physical structures, biotic components and abiotic conditions create a patchwork mosaic of potential niches for microbes. Questions remain about mechanisms and patterns of community assembly belowground, including: Do fungal and bacterial communities assemble differently? How do microbes reach the roots of host plants? Within a 4 m2 plot in alpine vegetation, high throughput sequencing of the 16S (bacteria) and ITS1 (fungal) ribosomal RNA genes was used to characterise microbial community composition in roots and adjacent soil of a viviparous host plant (Bistorta vivipara). At fine spatial scales, beta-diversity patterns in belowground bacterial and fungal communities were consistent, although compositional change was greater in bacteria than fungi. Spatial structure and distance-decay relationships were also similar for bacteria and fungi, with significant spatial structure detected at <50 cm among root- but not soil-associated microbes. Recruitment of root microbes from the soil community appeared limited at this sampling and sequencing depth. Possible explanations for this include recruitment from low-abundance populations of soil microbes, active recruitment from neighbouring plants and/or vertical transmission of symbionts to new clones, suggesting varied methods of microbial community assembly for viviparous plants. Our results suggest that even at relatively small spatial scales, deterministic processes play a significant role in belowground microbial community structure and assembly.

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.