Dadasnake, a Snakemake implementation of DADA2 to process amplicon sequencing data for microbial ecology

Modeling microbial impact on macrophyte debris decomposition in macrophyte-dominated eutrophic lakes

The decomposition of macrophytes plays a crucial role in the nutrient cycles of macrophyte-dominated eutrophication lakes. While research on plant decomposition mechanisms and microbial influences has rapid developed, it is curious that plant decomposition models have remained stagnant at the single-stage model from 50 years ago, without endeavor to consider any important factors. Our research conducted in-situ experiments and identified the optimal metrics for decomposition-related microbes, thereby establishing models for microbial impacts on decomposition rates (k_RDR). Using backward elimination in stepwise regression, we found that the optimal subset of independent variables-specifically Gammaproteobacteria-Q-L, Actinobacteriota-Q-L, and Ascomycota-Q-L-increased the adjusted R-squared (Ra2) to 0.93, providing the best modeling for decomposition rate (p = 0.002). Additionally, k_RDR can be modeled by synergic parameters of ACHB-Q-L, LDB-Q-L, and AB-Q-L for bacteria, and SFQ for fungi, albeit with a slightly lower Ra2 of 0.7-0.9 (p < 0.01). The primary contribution of our research lies in two key aspects. Firstly, we introduced optimal metrics for modeling microbes, opting for debris surface microbes over sediment microbes, and prioritizing absolute abundance over relative abundance. Secondly, our model represents a noteworthy advancement in debris modeling. Alongside elucidating the focus and innovative aspects of our work, we also addressed existing limitations and proposed directions for future research. SYNOPSIS: This study explores optimum metrics for decomposition-related microbes, offering precise microbial models for enhanced lake nutrient cycle simulation.

Structural properties and microbial diversity of the biofilm colonizing plastic substrates in Terra Nova Bay (Antarctica)

Microbial colonization on plastic polymers has been extensively explored, however the temporal dynamics of biofilm community in Antarctic environments are almost unknown. As a contribute to fill this knowledge gap, the structural characteristics and microbial diversity of the biofilm associated with polyvinyl chloride (PVC) and polyethylene (PE) panels submerged at 5 m of depth and collected after 3, 9 and 12 months were investigated in four coastal sites of the Ross Sea. Additional panels placed at 5 and 20 m were retrieved after 12 months. Chemical characterization was performed by FTIR-ATR and Raman (through Surface-Enhanced Raman Scattering, SERS) spectroscopy. Bacterial community composition was quantified at a single cell level by Catalyzed Reporter Deposition Fluorescence In Situ Hybridization (CARD-FISH) and Confocal Laser Scanning Microscopy (CLSM); microbial diversity was assessed by 16S rRNA gene sequencing. This multidisciplinary approach has provided new insights into microbial community dynamics during biofouling process, shedding light on the biofilm diversity and temporal succession on plastic substrates in the Ross Sea. Significant differences between free-living and microbial biofilm communities were found, with a more consolidated and structured community composition on PVC compared to PE. Spectral features ascribable to tyrosine, polysaccharides, nucleic acids and lipids characterized the PVC-associated biofilms. Pseudomonadota (among Gamma-proteobacteria) and Alpha-proteobacteria dominated the microbial biofilm community. Interestingly, in Road Bay, close to the Italian "Mario Zucchelli" research station, the biofilm growth - already observed during summer season, after 3 months of submersion - continued afterwards leading to a massive microbial abundance at the end of winter (after 12 months). After 3 months, higher percentages of Gamma-proteobacteria in Road Bay than in the not-impacted site were found. These observations lead us to hypothesize that in this site microbial fouling developed during the first 3 months could serve as a starter pioneering community stimulating the successive growth during winter.

A pile of pipelines: An overview of the bioinformatics software for metabarcoding data analyses

Environmental DNA (eDNA) metabarcoding has gained growing attention as a strategy for monitoring biodiversity in ecology. However, taxa identifications produced through metabarcoding require sophisticated processing of high-throughput sequencing data from taxonomically informative DNA barcodes. Various sets of universal and taxon-specific primers have been developed, extending the usability of metabarcoding across archaea, bacteria and eukaryotes. Accordingly, a multitude of metabarcoding data analysis tools and pipelines have also been developed. Often, several developed workflows are designed to process the same amplicon sequencing data, making it somewhat puzzling to choose one among the plethora of existing pipelines. However, each pipeline has its own specific philosophy, strengths and limitations, which should be considered depending on the aims of any specific study, as well as the bioinformatics expertise of the user. In this review, we outline the input data requirements, supported operating systems and particular attributes of thirty-two amplicon processing pipelines with the goal of helping users to select a pipeline for their metabarcoding projects.

Gut microbiome is not associated with mild cognitive impairment in Parkinson's disease

Gut microbiome differences between people with Parkinson's disease (PD) and control subjects without Parkinsonism are widely reported, but potential alterations related to PD with mild cognitive impairment (MCI) have yet to be comprehensively explored. We compared gut microbial features of PD with MCI (n = 58) to cognitively unimpaired PD (n = 60) and control subjects (n = 90) with normal cognition. Our results did not support a specific microbiome signature related to MCI in PD.

Chemical and microbiological insights into two littoral Antarctic demosponge species: Haliclona (Rhizoniera) dancoi (Topsent 1901) and Haliclona (Rhizoniera) scotti (Kirkpatrick 1907)

Introduction

Antarctic Porifera have gained increasing interest as hosts of diversified associated microbial communities that could provide interesting insights on the holobiome system and its relation with environmental parameters.

Methods

The Antarctic demosponge species Haliclona dancoi and Haliclona scotti were targeted for the determination of persistent organic pollutant (i. e., polychlorobiphenyls, PCBs, and polycyclic aromatic hydrocarbons, PAHs) and trace metal concentrations, along with the characterization of the associated prokaryotic communities by the 16S rRNA next generation sequencing, to evaluate possible relationships between pollutant accumulation (e.g., as a stress factor) and prokaryotic community composition in Antarctic sponges. To the best of our knowledge, this approach has been never applied before.

Results

Notably, both chemical and microbiological data on H. scotti (a quite rare species in the Ross Sea) are here reported for the first time, as well as the determination of PAHs in Antarctic Porifera. Both sponge species generally contained higher amounts of pollutants than the surrounding sediment and seawater, thus demonstrating their accumulation capability. The structure of the associated prokaryotic communities, even if differing at order and genus levels between the two sponge species, was dominated by Proteobacteria and Bacteroidota (with Archaea abundances that were negligible) and appeared in sharp contrast to communities inhabiting the bulk environment.

Discussions

Results suggested that some bacterial groups associated with H. dancoi and H. scotti were significantly (positively or negatively) correlated to the occurrence of certain contaminants.

Enzymatic machinery of wood-inhabiting fungi that degrade temperate tree species

Deadwood provides habitat for fungi and serves diverse ecological functions in forests. We already have profound knowledge of fungal assembly processes, physiological and enzymatic activities, and resulting physico-chemical changes during deadwood decay. However, in situ detection and identification methods, fungal origins, and a mechanistic understanding of the main lignocellulolytic enzymes are lacking. This study used metaproteomics to detect the main extracellular lignocellulolytic enzymes in 12 tree species in a temperate forest that have decomposed for 8 ½ years. Mainly white-rot (and few brown-rot) Basidiomycota were identified as the main wood decomposers, with Armillaria as the dominant genus; additionally, several soft-rot xylariaceous Ascomycota were identified. The key enzymes involved in lignocellulolysis included manganese peroxidase, peroxide-producing alcohol oxidases, laccase, diverse glycoside hydrolases (cellulase, glucosidase, xylanase), esterases, and lytic polysaccharide monooxygenases. The fungal community and enzyme composition differed among the 12 tree species. Ascomycota species were more prevalent in angiosperm logs than in gymnosperm logs. Regarding lignocellulolysis as a function, the extracellular enzyme toolbox acted simultaneously and was interrelated (e.g. peroxidases and peroxide-producing enzymes were strongly correlated), highly functionally redundant, and present in all logs. In summary, our in situ study provides comprehensive and detailed insight into the enzymatic machinery of wood-inhabiting fungi in temperate tree species. These findings will allow us to relate changes in environmental factors to lignocellulolysis as an ecosystem function in the future.

Education as Risk Factor of Mild Cognitive Impairment: The Link to the Gut Microbiome

BACKGROUND

With differences apparent in the gut microbiome in mild cognitive impairment (MCI) and dementia, and risk factors of dementia linked to alterations of the gut microbiome, the question remains if gut microbiome characteristics may mediate associations of education with MCI.

OBJECTIVES

We sought to examine potential mediation of the association of education and MCI by gut microbiome diversity or composition.

DESIGN

Cross-sectional study.

SETTING

Luxembourg, the Greater Region (surrounding areas in Belgium, France, Germany).

PARTICIPANTS

Control participants of the Luxembourg Parkinson's Study.

MEASUREMENTS

Gut microbiome composition, ascertained with 16S rRNA gene amplicon sequencing. Differential abundance, assessed across education groups (0-10, 11-16, 16+ years of education). Alpha diversity (Chao1, Shannon and inverse Simpson indices). Mediation analysis with effect decomposition was conducted with education as exposure, MCI as outcome and gut microbiome metrics as mediators.

RESULTS

After exclusion of participants below 50, or with missing data, n=258 participants (n=58 MCI) were included (M [SD] Age=64.6 [8.3] years). Higher education (16+ years) was associated with MCI (Odds ratio natural direct effect=0.35 [95% CI 0.15-0.81]. Streptococcus and Lachnospiraceae-UCG-001 genera were more abundant in higher education.

CONCLUSIONS

Education is associated with gut microbiome composition and MCI risk without clear evidence for mediation. However, our results suggest signatures of the gut microbiome that have been identified previously in AD and MCI to be reflected in lower education and suggest education as important covariate in microbiome studies.

Responses of rhizosphere fungi to the root economics space in grassland monocultures of different age

Recent studies on root traits have shown that there are two axes explaining trait variation belowground: the collaboration axis with mycorrhizal partners and the conservation ('fast - slow') axis. However, it is yet unknown whether these trait axes affect the assembly of soilborne fungi. We expect saprotrophic fungi to link to the conservation axis of root traits, whereas pathogenic and arbuscular mycorrhizal fungi link to the collaboration axis, but in opposite directions, as arbuscular mycorrhizal fungi might provide pathogen protection. To test these hypotheses, we sequenced rhizosphere fungal communities and measured root traits in monocultures of 25 grassland plant species, differing in age. Within the fungal guilds, we evaluated fungal species richness, relative abundance and community composition. Contrary to our hypotheses, fungal diversity and relative abundance were not strongly related to the root trait axes. However, saprotrophic fungal community composition was affected by the conservation gradient and pathogenic community composition by the collaboration gradient. The rhizosphere AMF community composition did not change along the collaboration gradient, even though the root trait axis was in line with the root mycorrhizal colonization rate. Overall, our results indicate that in the long term, the root trait axes are linked with fungal community composition.

Nitrogen addition increases mass loss of gymnosperm but not of angiosperm deadwood without changing microbial communities

Enhanced nitrogen (N) deposition due to combustion of fossil fuels and agricultural fertilization is a global phenomenon which has severely altered carbon (C) and N cycling in temperate forest ecosystems in the northern hemisphere. Although deadwood holds a substantial amount of C in forest ecosystems and thus plays a crucial role in nutrient cycling, the effect of increased N deposition on microbial processes and communities, wood chemical traits and deadwood mass loss remains unclear. Here, we simulated high N deposition rates by adding reactive N in form of ammonium-nitrate (40 kg N ha-1 yr-1) to deadwood of 13 temperate tree species over nine years in a field experiment in Germany. Non-treated deadwood from the same logs served as control with background N deposition. Our results show that chronically elevated N levels alters deadwood mass loss alongside respiration, enzymatic activities and wood chemistry depending on tree clade and species. In gymnosperm deadwood, elevated N increased mass loss by +38 %, respiration by +37 % and increased laccase activity 12-fold alongside increases of white-rot fungal abundance +89 % (p = 0.03). Furthermore, we observed marginally significant (p = 0.06) shifts of bacterial communities in gymnosperm deadwood. In angiosperm deadwood, we did not detect consistent effects on mass loss, physico-chemical properties, extracellular enzymatic activity or changes in microbial communities except for changes in abundance of 10 fungal OTUs in seven tree species and 28 bacterial OTUs in 10 tree species. We conclude that N deposition alters decomposition processes exclusively in N limited gymnosperm deadwood in the long term by enhancing fungal activity as expressed by increases in respiration rate and extracellular enzyme activity with minor shifts in decomposing microbial communities. By contrast, deadwood of angiosperm tree species had higher N concentrations and mass loss as well as community composition did not respond to N addition.

A Deep Insight into the Diversity of Microfungal Communities in Arctic and Antarctic Lakes

We assessed fungal diversity in water and sediment samples obtained from five Arctic lakes in Ny-Ålesund (Svalbard Islands, High Arctic) and five Antarctic lakes on Livingston and Deception Islands (South Shetland Islands), using DNA metabarcoding. A total of 1,639,074 fungal DNA reads were detected and assigned to 5980 ASVs amplicon sequence variants (ASVs), with only 102 (1.7%) that were shared between the two Polar regions. For Arctic lakes, unknown fungal taxa dominated the sequence assemblages, suggesting the dominance of possibly undescribed fungi. The phylum Chytridiomycota was the most represented in the majority of Arctic and Antarctic samples, followed by Rozellomycota, Ascomycota, Basidiomycota, and the less frequent Monoblepharomycota, Aphelidiomycota, Mortierellomycota, Mucoromycota, and Neocallimastigomycota. At the genus level, the most abundant genera included psychrotolerant and cosmopolitan cold-adapted fungi including Alternaria, Cladosporium, Cadophora, Ulvella (Ascomycota), Leucosporidium, Vishniacozyma (Basidiomycota), and Betamyces (Chytridiomycota). The assemblages displayed high diversity and richness. The assigned diversity was composed mainly of taxa recognized as saprophytic fungi, followed by pathogenic and symbiotic fungi.

FungalTraits vs. FUNGuild: Comparison of Ecological Functional Assignments of Leaf- and Needle-Associated Fungi Across 12 Temperate Tree Species

Recently, a new annotation tool "FungalTraits" was created based on the previous FUNGuild and Fun^Fun databases, which has attracted high attention in the scientific community. These databases were widely used to gain more information from fungal sequencing datasets by assigning fungal functional traits. More than 1500 publications so far employed FUNGuild and the aim of this study is to compare this successful database with the recent FungalTraits database. Quality and quantity of the assignment by FUNGuild and FungalTraits to a fungal internal transcribed spacer (ITS)-based amplicon sequencing dataset on amplicon sequence variants (ASVs) were addressed. Sequencing dataset was derived from leaves and needles of 12 temperate broadleaved and coniferous tree species. We found that FungalTraits assigned more functional traits than FUNGuild, and especially the coverage of saprotrophs, plant pathogens, and endophytes was higher while lichenized fungi revealed similar findings. Moreover, ASVs derived from leaves and needles of each tree species were better assigned to all available fungal traits as well as to saprotrophs by FungalTraits compared to FUNGuild in particular for broadleaved tree species. Assigned ASV richness as well as fungal functional community composition was higher and more diverse after analyses with FungalTraits compared to FUNGuild. Moreover, datasets of both databases showed similar effect of environmental factors for saprotrophs but for endophytes, unidentical patterns of significant corresponding factors were obtained. As a conclusion, FungalTraits is superior to FUNGuild in assigning a higher quantity and quality of ASVs as well as a higher frequency of significant correlations with environmental factors.

Future climate conditions accelerate wheat straw decomposition alongside altered microbial community composition, assembly patterns, and interaction networks

Although microbial decomposition of plant litter plays a crucial role in nutrient cycling and soil fertility, we know less about likely links of specific microbial traits and decomposition, especially in relation to climate change. We study here wheat straw decomposition under ambient and manipulated conditions simulating a future climate scenario (next 80 years) in agroecosystems, including decay rates, macronutrient dynamics, enzyme activity, and microbial communities. We show that future climate will accelerate straw decay rates only during the early phase of the decomposition process. Additionally, the projected climate change will increase the relative abundance of saprotrophic fungi in decomposing wheat straw. Moreover, the impact of future climate on microbial community assembly and molecular ecological networks of both bacteria and fungi will strongly depend on the decomposition phase. During the early phase of straw decomposition, stochastic processes dominated microbial assembly under ambient climate conditions, whereas deterministic processes highly dominated bacterial and fungal communities under simulated future climate conditions. In the later decomposition phase, similar assembly processes shaped the microbial communities under both climate scenarios. Furthermore, over the early phases of decomposition, simulated future climate enhanced the complexity of microbial interaction networks. We concluded that the impact of future climate on straw decay rate and associated microbial traits like assembly processes and inter-community interactions is restricted to the early phase of decomposition.

Bark-inhabiting fungal communities of European chestnut undergo substantial alteration by canker formation following chestnut blight infection

Background

Chestnut forests are severely threatened by chestnut blight caused by the fungal pathogen Cryphonectria parasitica and the infected trees exhibit bark canker in the later stage of the disease. European chestnut (Castanea sativa) is further infected by Gnomoniopsis smithogilvyi, another canker-causing fungal pathogen. We explored whether and how chestnut blight is reflected in bark-inhabiting fungal communities of European chestnut and also assessed the co-occurrence of C. parasitica and G. smithogilvyi.

Materials and methods

We initially investigated the fungal communities of European chestnut bark tissues and further monitored changes in these fungal communities with regard to disease progression from infection to canker formation by analyzing bark samples from asymptomatic trees, asymptomatic trees with latent C. parasitica infection, and infected trees with canker tissues, using amplicon sequencing of the ITS2 region of rDNA.

Results

The results showed that fungal community composition and diversity differed between the sample types. The fungal community composition was substantially reshaped by canker formation, whereas latent C. parasitica infection and more specifically pre-canker infection period per se had a weak effect. Fungal communities of canker samples was less diverse and more dissimilar to those of other sample types. C. parasitica dominated the mycobiome of canker samples, whereas G. smithogilvyi was found in only 9% of canker samples at very low abundances. However, G. smithogilvyi was a dominant fungus in the bark of healthy plants.

Conclusion

This study highlights that canker formation is the principal driver of decreasing diversity and altered composition of the mycobiome in bark tissues of European chestnut infected by C. parasitica infection. It additionally emphasizes the scarce co-occurrence of C. parasitica and G. smithogilvyi on European chestnut.

Amplicon Sequencing Pipelines in Metagenomics

Taxonomic profiling among a large number of samples is a fundamental task during amplicon sequencing analysis. The heterogeneity and technical noises in the sample handling, library preparation, and sequencing present a major challenge to how the biological conclusions are drawn from the data analysis, and accordingly, many tools have been developed to address specific issues related to each step of the data analysis. Nowadays, several sophisticated computational pipelines with flexible parameters are made available to provide one-stop comprehensive solutions by integrating various tools, which significantly mitigate the burden imposed by the complexity of the metagenomics data analysis. This chapter discusses the best practices related to the data generation and describes bioinformatics approaches to achieving greater accuracy from data processing. It offers two independent stepwise pipelines using mothur and DADA2 in a parallel way, presents the basic principles in the key steps of the analysis, and enables the comparisons between the two pipelines straightforwardly.

Evaluation of primers for the detection of deadwood-inhabiting archaea via amplicon sequencing

Archaea have been reported from deadwood of a few different tree species in temperate and boreal forest ecosystems in the past. However, while one of their functions is well linked to methane production any additional contribution to wood decomposition is not understood and underexplored which may be also attributed to lacking investigations on their diversity in this substrate. With this current work, we aim at encouraging further investigations by providing aid in primer choice for DNA metabarcoding using Illumina amplicon sequencing. We tested 16S primer pairs on genomic DNA extracted from woody tissue of four temperate deciduous tree species. Three primer pairs were specific to archaea and one prokaryotic primer pair theoretically amplifies both, bacterial and archaeal DNA. Methanobacteriales and Methanomassiliicoccales have been consistently identified as dominant orders across all datasets but significant variability in ASV richness was observed using different primer combinations. Nitrososphaerales have only been identified when using archaea-specific primer sets. In addition, the most commonly applied primer combination targeting prokaryotes in general yielded the lowest relative proportion of archaeal sequences per sample, which underlines the fact, that using target specific primers unraveled a yet unknown diversity of archaea in deadwood. Hence, archaea seem to be an important group of the deadwood-inhabiting community and further research is needed to explore their role during the decomposition process.

A poisoned apple: First insights into community assembly and networks of the fungal pathobiome of healthy-looking senescing leaves of temperate trees in mixed forest ecosystem

Despite the abundance of observations of foliar pathogens, our knowledge is severely lacking regarding how the potential fungal pathobiome is structured and which processes determine community assembly. In this study, we addressed these questions by analysing the potential fungal pathobiome associated with the senescing leaves and needles of 12 temperate tree species. We compared fungal plant pathogen load in the senescing leaves/needles and demonstrated that healthy-looking leaves/needles are inhabited by diverse and distinct fungal plant pathogens. We detected 400 fungal plant pathogenic ASVs belonging to 130 genera. The fungal plant pathogenic generalist, Mycosphaerella, was found to be the potential most significant contributor to foliar disease in seedlings. The analyses of assembly process and co-occurrence network showed that the fungal plant pathogenic communities in different tree types are mainly determined by stochastic processes. However, the homogenising dispersal highly contributes in broadleaf trees, whereas ecological drift plays an important role in coniferious trees. The deterministic assembly processes (dominated by variable selection) contributed more in broadleaf trees as compared to coniferous trees. We found that pH and P level significantly corresponded with fungal plant pathogenic community compositions in both tree types. Our study provides the first insight and mechanistic understanding into the community assembly, networks, and complete taxonomy of the foliar fungal pathobiome in senescing leaves and needles.

Dark side of a bio-based and biodegradable plastic? Assessment of pathogenic microbes associated with poly(butylene succinate-co-adipate) under ambient and future climates using next-generation sequencing

Bio-based and biodegradable plastic mulching films have been proposed to replace the non-biodegradable plastic mulch films to solve plastic pollution problems in agricultural soils. However, the impact of bio-based and biodegradable plastics on plant and human health remains largely unexplored. Here, we aimed to assess the risk under field conditions of a bio-based and biodegradable poly(butylene succinate-co-adipate; PBSA), a widely used mulching film as carrier of potential pathogenic microorganisms (bacteria and fungi) at ambient and future climate conditions. Overall, we affiliated 64 fungal and 11 bacterial operational taxonomic units (OTUs) as pathogens by using Next-Generation Sequencing approach. Our results revealed that PBSA hosted at least 53 plant pathogens, of which 51 were classified as fungi, while the other two were bacteria. Most fungal plant pathogens were able to withstand the anticipated future climate changes. We detected 13 fungal and eight bacterial OTUs, which were classified as opportunistic human pathogens. Only one bacterial OTU (Enterococcus faecium) was assigned to a human pathogen. While future climate conditions only significantly impacted on the presence and frequency of detection of few pathogens, incubation time was found to significantly impacted on nine pathogens. This result demonstrates the temporal dynamics of pathogens associated with PBSA. The threats to plant and human health were discussed. We emphasize that the risks to human health are relatively low because we mainly found opportunistic pathogens associated with PBSA and the amount are comparable to the plant debris. However, the risks to plant health may be considered as moderate because many plant pathogens were discovered and/or enriched in PBSA. Furthermore, in soil environments, the pathogenic risk of plastic is highly depending on the surrounding soil pathobiome where plastic is being decomposed.

More than you can see: Unraveling the ecology and biodiversity of lichenized fungi associated with leaves and needles of 12 temperate tree species using high-throughput sequencing

Currently, lichen surveys are generally based on the examination of fruiting bodies. Lichens in the mycelial stage, in spores, or awaiting conditions for fruiting body formation are usually overlooked, even though they are important for maintaining biodiversity and ecosystem functions. This study aimed to explore the lichenized fungal community composition and richness associated with leaves and needles of 12 temperate tree species using Illumina MiSeq-based amplicon sequencing of the internal transcribed spacer (ITS) 2 region. Picea abies harbored the highest richness and number of lichenized fungal species. We found that the lichenized fungus Physcia adscendens dominated the leaves and needles of the most temperate tree species. Eleven lichenized fungal species detected in this study were recorded for the first time on leaves and needles. In addition, we identified Athallia cerinella, Fellhanera bouteillei, and Melanohalea exasperata that are on the German national red lists. Lichenized fungal richness was higher in conifer compared to broadleaf trees. Overall, tree species (within coniferous trees) and tree types (broadleaved vs. coniferous trees) harbored significantly different lichenized fungal community compositions pointing out the importance of host species. Diversity and community composition patterns of lichenized fungi were correlated mainly with tree species. Our study demonstrates that the diversity of foliicolous lichens associated with leaves and needles of 12 temperate tree species can be appropriately analyzed and functionally assigned using the ITS-based high-throughput sequencing. We highlighted the importance of conifers for maintaining the biodiversity of foliicolous lichens. Based on the discovery of many red list lichens, our methodological approach and results are important contributions to subsequent actions in the bio-conversation approaches.

Antarctic Salt-Cones: An Oasis of Microbial Life? The Example of Boulder Clay Glacier (Northern Victoria Land)

The evaporation of a localized, highly saline water body of the Boulder Clay debris-covered glacier, in the Northern Victoria Land, probably generated the accumulation of mirabilite (Na₂SO₄ × 10H₂O) and thenardite (Na₂SO₄) in a glacier salt-cone. Such an extremely cold and salty environment resembles the conditions on Mars, so it can be considered a terrestrial analog. The study was aimed at gaining a first glimpse at the prokaryotic community associated with Antarctic mirabilite and thenardite minerals and also to find clues about the origin of the salts. For this purpose, samples were analyzed by a next generation approach to investigate the prokaryotic (Bacteria and Archaea) diversity. Phylogenetic analysis allowed the identification of Bacteroidota, Actinobacteriota, Firmicutes, and Gammaproteobacteria as the main bacterial lineages, in addition to Archaea in the phylum Halobacterota. The genera Arthrobacter, Rhodoglobus, Gillisia, Marinobacter and Psychrobacter were particularly abundant. Interestingly, several bacterial and archaeal sequences were related to halotolerant and halophilic genera, previously reported in a variety of marine environments and saline habitats, also in Antarctica. The analyzed salt community also included members that are believed to play a major role in the sulfur cycle.

Tourmaline: A containerized workflow for rapid and iterable amplicon sequence analysis using QIIME 2 and Snakemake

Background

Amplicon sequencing (metabarcoding) is a common method to survey diversity of environmental communities whereby a single genetic locus is amplified and sequenced from the DNA of whole or partial organisms, organismal traces (e.g., skin, mucus, feces), or microbes in an environmental sample. Several software packages exist for analyzing amplicon data, among which QIIME 2 has emerged as a popular option because of its broad functionality, plugin architecture, provenance tracking, and interactive visualizations. However, each new analysis requires the user to keep track of input and output file names, parameters, and commands; this lack of automation and standardization is inefficient and creates barriers to meta-analysis and sharing of results.

Findings

We developed Tourmaline, a Python-based workflow that implements QIIME 2 and is built using the Snakemake workflow management system. Starting from a configuration file that defines parameters and input files—a reference database, a sample metadata file, and a manifest or archive of FASTQ sequences—it uses QIIME 2 to run either the DADA2 or Deblur denoising algorithm; assigns taxonomy to the resulting representative sequences; performs analyses of taxonomic, alpha, and beta diversity; and generates an HTML report summarizing and linking to the output files. Features include support for multiple cores, automatic determination of trimming parameters using quality scores, representative sequence filtering (taxonomy, length, abundance, prevalence, or ID), support for multiple taxonomic classification and sequence alignment methods, outlier detection, and automated initialization of a new analysis using previous settings. The workflow runs natively on Linux and macOS or via a Docker container. We ran Tourmaline on a 16S ribosomal RNA amplicon data set from Lake Erie surface water, showing its utility for parameter optimization and the ability to easily view interactive visualizations through the HTML report, QIIME 2 viewer, and R- and Python-based Jupyter notebooks.

Conclusion

Automated workflows like Tourmaline enable rapid analysis of environmental amplicon data, decreasing the time from data generation to actionable results. Tourmaline is available for download at github.com/aomlomics/tourmaline.

Effects of Tree Composition and Soil Depth on Structure and Functionality of Belowground Microbial Communities in Temperate European Forests

Depending on their tree species composition, forests recruit different soil microbial communities. Likewise, the vertical nutrient gradient along soil profiles impacts these communities and their activities. In forest soils, bacteria and fungi commonly compete, coexist, and interact, which is challenging for understanding the complex mechanisms behind microbial structuring. Using amplicon sequencing, we analyzed bacterial and fungal diversity in relation to forest composition and soil depth. Moreover, employing random forest models, we identified microbial indicator taxa of forest plots composed of either deciduous or evergreen trees, or their mixtures, as well as of three soil depths. We expected that forest composition and soil depth affect bacterial and fungal diversity and community structure differently. Indeed, relative abundances of microbial communities changed more across soil depths than in relation to forest composition. The microbial Shannon diversity was particularly affected by soil depth and by the proportion of evergreen trees. Our results also reflected that bacterial communities are primarily shaped by soil depth, while fungi were influenced by forest tree species composition. An increasing proportion of evergreen trees did not provoke differences in main bacterial metabolic functions, e.g., carbon fixation, degradation, or photosynthesis. However, significant responses related to specialized bacterial metabolisms were detected. Saprotrophic, arbuscular mycorrhizal, and plant pathogenic fungi were related to the proportion of evergreen trees, particularly in topsoil. Prominent microbial indicator taxa in the deciduous forests were characterized to be r-strategists, whereas K-strategists dominated evergreen plots. Considering simultaneously forest composition and soil depth to unravel differences in microbial communities, metabolic pathways and functional guilds have the potential to enlighten mechanisms that maintain forest soil functionality and provide resistance against disturbances.

Snaq: A Dynamic Snakemake Pipeline for Microbiome Data Analysis With QIIME2

Optimizing and automating a protocol for 16S microbiome data analysis with QIIME2 is a challenging task. It involves a multi-step process, and multiple parameters and options that need to be tested and determined. In this article, we describe Snaq, a snakemake pipeline that helps automate and optimize 16S data analysis using QIIME2. Snaq offers an informative file naming system and automatically performs the analysis of a data set by downloading and installing the required databases and classifiers, all through a single command-line instruction. It works natively on Linux and Mac and on Windows through the use of containers, and is potentially extendable by adding new rules. This pipeline will substantially reduce the efforts in sending commands and prevent the confusion caused by the accumulation of analysis results due to testing multiple parameters.

Water Deficit History Selects Plant Beneficial Soil Bacteria Differently Under Conventional and Organic Farming

Water deficit tolerance is critical for plant fitness and survival, especially when successive drought events happen. Specific soil microorganisms are however able to improve plant tolerance to stresses, such as those displaying a 1-aminocyclopropane-1-carboxylate (ACC) deaminase activity. Microorganisms adapted to dry conditions can be selected by plants over time because of properties such as sporulation, substrate preference, or cell-wall thickness. However, the complexity and interconnection between abiotic factors, like drought or soil management, and biotic factors, like plant species identity, make it difficult to elucidate the general selection processes of such microorganisms. Using a pot experiment in which wheat and barley were grown on conventional and organic farming soils, we determined the effect of water deficit history on soil microorganisms by comparing single and successive events of water limitation. The analysis showed that water deficit strongly impacts the composition of both the total microbial community (16S rRNA genes) and one of ACC deaminase-positive (acdS+) microorganisms in the rhizosphere. In contrast, successive dry conditions moderately influence the abundance and diversity of both communities compared to a single dry event. We revealed interactive effects of the farming soil type and the water deficit conditioning treatment. Indeed, possibly due to better nutrient status, plants grown on soils from conventional farming showed higher growth and were able to select more adapted microbial taxa. Some of them are already known for their plant-beneficial properties like the Actinobacteria Streptomyces, but interestingly, some Proteobacteria were also enriched after a water deficit history under conventional farming. Our approach allowed us to identify key microbial taxa promoting drought adaptation of cereals, thus improving our understanding of drought effects on plant-microbe interactions.

Impact of Soil Inoculation with Bacillus amyloliquefaciens FZB42 on the Phytoaccumulation of Germanium, Rare Earth Elements, and Potentially Toxic Elements

Bioaugmentation promises benefits for agricultural production as well as for remediation and phytomining approaches. Thus, this study investigated the effect of soil inoculation with the commercially available product RhizoVital^®42, which contains Bacillus amyloliquefaciens FZB42, on nutrient uptake and plant biomass production as well as on the phytoaccumulation of potentially toxic elements, germanium, and rare earth elements (REEs). Zea mays and Fagopyrum esculentum were selected as model plants, and after harvest, the element uptake was compared between plants grown on inoculated versus reference soil. The results indicate an enrichment of B. amyloliquefaciens in inoculated soils as well as no significant impact on the inherent bacterial community composition. For F. esculentum, inoculation increased the accumulation of most nutrients and As, Cu, Pb, Co, and REEs (significant for Ca, Cu, and Co with 40%, 2042%, and 383%, respectively), while it slightly decreased the uptake of Ge, Cr, and Fe. For Z. mays, soil inoculation decreased the accumulation of Cr, Pb, Co, Ge, and REEs (significant for Co with 57%) but showed an insignificant increased uptake of Cu, As, and nutrient elements. Summarily, the results suggest that bioaugmentation with B. amyloliquefaciens is safe and has the potential to enhance/reduce the phytoaccumulation of some elements and the effects of inoculation are plant specific.

Molecular Screening of Microorganisms Associated with Discolored Wood in Dead European Beech Trees Suffered from Extreme Drought Event Using Next Generation Sequencing

Drought events weaken trees and make them vulnerable to attacks by diverse plant pathogens. Here, we propose a molecular method for fast screening of microorganisms associated with European beech decline after an extreme drought period (2018) in a forest of Thuringia, Germany. We used Illumina sequencing with a recent bioinformatics approach based on DADA2 to identify archaeal, bacterial, and fungal ASVs (amplicon sequence variants) based on bacterial and archaeal 16S and fungal ITS genes. We show that symptomatic beech trees are associated with both bacterial and fungal plant pathogens. Although the plant pathogen sequences were detected in both discolored and non-discolored wood areas, they were highly enriched in the discolored wood areas. We show that almost each individual tree was associated with a different combination of pathogens. Cytospora spp. and Neonectria coccinea were among the most frequently detected fungal pathogens, whereas Erwinia spp. and Pseudomonas spp. were the dominant bacterial plant pathogens. We demonstrate that bacterial plant pathogens may be of major importance in beech decline.

Distinct effects of host and neighbour tree identity on arbuscular and ectomycorrhizal fungi along a tree diversity gradient

Plant diversity and plant-related ecosystem functions have been important in biodiversity-ecosystem functioning studies. However, biotic interactions with mycorrhizal fungi have been understudied although they are crucial for plant-resource acquisition. Here, we investigated the effects of tree species richness and tree mycorrhizal type on arbuscular (AMF) and ectomycorrhizal fungal (EMF) communities. We aimed to understand how dissimilarities in taxa composition and beta-diversity are related to target trees and neighbours of the same or different mycorrhizal type. We sampled a tree diversity experiment with saplings (~7 years old), where tree species richness (monocultures, 2-species, and 4-species mixtures) and mycorrhizal type were manipulated. AMF and EMF richness significantly increased with increasing tree species richness. AMF richness of mixture plots resembled that of the sum of the respective monocultures, whereas EMF richness of mixture plots was lower compared to the sum of the respective monocultures. Specialisation scores revealed significantly more specialised AMF than EMF suggesting that, in contrast to previous studies, AMF were more specialised, whereas EMF were not. We further found that AMF communities were little driven by the surrounding trees, whereas EMF communities were. Our study revealed drivers of mycorrhizal fungal communities and further highlights the distinct strategies of AMF and EMF.

First Evidence That Nematode Communities in Deadwood Are Related to Tree Species Identity and to Co-Occurring Fungi and Prokaryotes

Nematodes represent a diverse and ubiquitous group of metazoans in terrestrial environments. They feed on bacteria, fungi, plants, other nematodes or parasitize a variety of animals and hence may be considered as active members of many food webs. Deadwood is a structural component of forest ecosystems which harbors many niches for diverse biota. As fungi and bacteria are among the most prominent decomposing colonizers of deadwood, we anticipated frequent and diverse nematode populations to co-occur in such ecosystems. However, knowledge about their ability to colonize this habitat is still limited. We applied DNA-based amplicon sequencing (metabarcoding) of the 18S rRNA gene to analyze nematode communities in sapwood and heartwood of decaying logs from 13 different tree species. We identified 247 nematode ASVs (amplicon sequence variants) from 27 families. Most of these identified families represent bacterial and fungal feeders. Their composition strongly depended on tree species identity in both wood compartments. While pH and water content were the only wood properties that contributed to nematodes' distribution, co-occurring fungal and prokaryotic (bacteria and archaea) α- and β-diversities were significantly related to nematode communities. By exploring thirteen different tree species, which exhibit a broad range of wood characteristics, this study provides first and comprehensive insights into nematode diversity in deadwood of temperate forests and indicates connectivity to other wood-inhabiting organisms.

Interactions between soil properties, agricultural management and cultivar type drive structural and functional adaptations of the wheat rhizosphere microbiome to drought

Rhizosphere microbial communities adapt their structural and functional compositions to water scarcity and have the potential to substantially mitigate drought stress of crops. To unlock this potential, it is crucial to understand community responses to drought in the complex interplay between soil properties, agricultural management and crop species. Two winter wheat cultivars, demanding and non-demanding, were exposed to drought stress in loamy Chernozem and sandy Luvisol soils under conventional or organic farming management. Structural and functional adaptations of the rhizosphere bacteria were assessed by 16S amplicon sequencing, the predicted abundance of drought-related functional genes in the bacterial community based on 16S amplicon sequences (Tax4Fun) and the activity potentials of extracellular enzymes involved in the carbon cycle. Bacterial community composition was strongly driven by drought and soil type. Under drought conditions, Gram-positive phyla became relatively more abundant, but either less or more diverse in Luvisol and Chernozem soil respectively. Enzyme activities and functional gene abundances related to carbon degradation were increased under drought in the rhizosphere of the demanding wheat cultivar in organic farming. We demonstrate that soil type, farming system and wheat cultivar each constitute important factors during the structural and/or functional adaptation of rhizobacterial communities in response to drought.

Dadaist2: A Toolkit to Automate and Simplify Statistical Analysis and Plotting of Metabarcoding Experiments

The taxonomic composition of microbial communities can be assessed using universal marker amplicon sequencing. The most common taxonomic markers are the 16S rDNA for bacterial communities and the internal transcribed spacer (ITS) region for fungal communities, but various other markers are used for barcoding eukaryotes. A crucial step in the bioinformatic analysis of amplicon sequences is the identification of representative sequences. This can be achieved using a clustering approach or by denoising raw sequencing reads. DADA2 is a widely adopted algorithm, released as an R library, that denoises marker-specific amplicons from next-generation sequencing and produces a set of representative sequences referred to as 'Amplicon Sequence Variants' (ASV). Here, we present Dadaist2, a modular pipeline, providing a complete suite for the analysis that ranges from raw sequencing reads to the statistics of numerical ecology. Dadaist2 implements a new approach that is specifically optimised for amplicons with variable lengths, such as the fungal ITS. The pipeline focuses on streamlining the data flow from the command line to R, with multiple options for statistical analysis and plotting, both interactive and automatic.

Organic agricultural practice enhances arbuscular mycorrhizal symbiosis in correspondence to soil warming and altered precipitation patterns

Climate and agricultural practice interact to influence both crop production and soil microbes in agroecosystems. Here, we carried out a unique experiment in Central Germany to simultaneously investigate the effects of climates (ambient climate vs. future climate expected in 50-70 years), agricultural practices (conventional vs. organic farming), and their interaction on arbuscular mycorrhizal fungi (AMF) inside wheat (Triticum aestivum L.) roots. AMF communities were characterized using Illumina sequencing of 18S rRNA gene amplicons. We showed that climatic conditions and agricultural practices significantly altered total AMF community composition. Conventional farming significantly affected the AMF community and caused a decline in AMF richness. Factors shaping AMF community composition and richness at family level differed greatly among Glomeraceae, Gigasporaceae and Diversisporaceae. An interactive impact of climate and agricultural practices was detected in the community composition of Diversisporaceae. Organic farming mitigated the negative effect of future climate and promoted total AMF and Gigasporaceae richness. AMF richness was significantly linked with nutrient content of wheat grains under both agricultural practices.

Soil Texture, Sampling Depth and Root Hairs Shape the Structure of ACC Deaminase Bacterial Community Composition in Maize Rhizosphere

Preservation of the phytostimulatory functions of plant growth-promoting bacteria relies on the adaptation of their community to the rhizosphere environment. Here, an amplicon sequencing approach was implemented to specifically target microorganisms with 1-aminocyclopropane-1-carboxylate deaminase activity, carrying the acdS gene. We stated the hypothesis that the relative phylogenetic distribution of acdS carrying microorganisms is affected by the presence or absence of root hairs, soil type, and depth. To this end, a standardized soil column experiment was conducted with maize wild type and root hair defective rth3 mutant in the substrates loam and sand, and harvest was implemented from three depths. Most acdS sequences (99%) were affiliated to Actinobacteria and Proteobacteria, and the strongest influence on the relative abundances of sequences were exerted by the substrate. Variovorax, Acidovorax, and Ralstonia sequences dominated in loam, whereas Streptomyces and Agromyces were more abundant in sand. Soil depth caused strong variations in acdS sequence distribution, with differential levels in the relative abundances of acdS sequences affiliated to Tetrasphaera, Amycolatopsis, and Streptomyces in loam, but Burkholderia, Paraburkholderia, and Variovorax in sand. Maize genotype influenced the distribution of acdS sequences mainly in loam and only in the uppermost depth. Variovorax acdS sequences were more abundant in WT, but Streptomyces, Microbacterium, and Modestobacter in rth3 rhizosphere. Substrate and soil depth were strong and plant genotype a further significant single and interacting drivers of acdS carrying microbial community composition in the rhizosphere of maize. This suggests that maize rhizosphere acdS carrying bacterial community establishes according to the environmental constraints, and that root hairs possess a minor but significant impact on acdS carrying bacterial populations.