Natrix: a Snakemake-based workflow for processing, clustering, and taxonomically assigning amplicon sequencing reads.

BACKGROUND: Sequencing of marker genes amplified from environmental samples, known as amplicon sequencing, allows us to resolve some of the hidden diversity and elucidate evolutionary relationships and ecological processes among complex microbial communities. The analysis of large numbers of samples at high sequencing depths generated by high throughput sequencing technologies requires efficient, flexible, and reproducible bioinformatics pipelines. Only a few existing workflows can be run in a user-friendly, scalable, and reproducible manner on different computing devices using an efficient workflow management system. RESULTS: We present Natrix, an open-source bioinformatics workflow for preprocessing raw amplicon sequencing data. The workflow contains all analysis steps from quality assessment, read assembly, dereplication, chimera detection, split-sample merging, sequence representative assignment (OTUs or ASVs) to the taxonomic assignment of sequence representatives. The workflow is written using Snakemake, a workflow management engine for developing data analysis workflows. In addition, Conda is used for version control. Thus, Snakemake ensures reproducibility and Conda offers version control of the utilized programs. The encapsulation of rules and their dependencies support hassle-free sharing of rules between workflows and easy adaptation and extension of existing workflows. Natrix is freely available on GitHub ( https://github.com/MW55/Natrix ) or as a Docker container on DockerHub ( https://hub.docker.com/r/mw55/natrix ). CONCLUSION: Natrix is a user-friendly and highly extensible workflow for processing Illumina amplicon data.

Centers of endemism of freshwater protists deviate from pattern of taxon richness on a continental scale.

Here, we analyzed patterns of taxon richness and endemism of freshwater protists in Europe. Even though the significance of physicochemical parameters but also of geographic constraints for protist distribution is documented, it remains unclear where regional areas of high protist diversity are located and whether areas of high taxon richness harbor a high proportion of endemics. Further, patterns may be universal for protists or deviate between taxonomic groups. Based on amplicon sequencing campaigns targeting the SSU and ITS region of the rDNA we address these patterns at two different levels of phylogenetic resolution. Our analyses demonstrate that protists have restricted geographical distribution areas. For many taxonomic groups the regions of high taxon richness deviate from those having a high proportion of putative endemics. In particular, the diversity of high mountain lakes as azonal habitats deviated from surrounding lowlands, i.e. many taxa were found exclusively in high mountain lakes and several putatively endemic taxa occurred in mountain regions like the Alps, the Pyrenees or the Massif Central. Beyond that, taxonomic groups showed a pronounced accumulation of putative endemics in distinct regions, e.g. Dinophyceae along the Baltic Sea coastline, and Chrysophyceae in Scandinavia. Many other groups did not have pronounced areas of increased endemism but geographically restricted taxa were found across Europe.

Factors shaping community patterns of protists and bacteria on a European scale.

Factors shaping community patterns of microorganisms are controversially discussed. Physical and chemical factors certainly limit the survival of individual taxa and maintenance of diversity. In recent years, a contribution of geographic distance and dispersal barriers to distribution patterns of protists and bacteria has been demonstrated. Organismic interactions such as competition, predation and mutualism further modify community structure and maintenance of distinct taxa. Here, we address the relative importance of these different factors in shaping protists and bacterial communities on a European scale using high-throughput sequencing data obtained from lentic freshwater ecosystems. We show that community patterns of protists are similar to those of bacteria. Our results indicate that cross-domain organismic factors are important variables with a higher influence on protists as compared with bacteria. Abiotic physical and chemical factors also contributed significantly to community patterns. The contribution of these latter factors was higher for bacteria, which may reflect a stronger biogeochemical coupling. The contribution of geographical distance was similar for both microbial groups.

Synchrony of Eukaryotic and Prokaryotic Planktonic Communities in Three Seasonally Sampled Austrian Lakes.

Freshwater systems are characterized by an enormous diversity of eukaryotic protists and prokaryotic taxa. The community structures in different lakes are thereby influenced by factors such as habitat size, lake chemistry, biotic interactions, and seasonality. In our study, we used high throughput 454 sequencing to study the diversity and temporal changes of prokaryotic and eukaryotic planktonic communities in three Austrian lakes during the ice-free season. In the following year, one lake was sampled again with a reduced set of sampling dates to observe reoccurring patterns. Cluster analyses (based on SSU V9 (eukaryotic) and V4 (prokaryotic) OTU composition) grouped samples according to their origin followed by separation into seasonal clusters, indicating that each lake has a unique signature based on OTU composition. These results suggest a strong habitat-specificity of microbial communities and in particular of community patterns at the OTU level. A comparison of the prokaryotic and eukaryotic datasets via co-inertia analysis (CIA) showed a consistent clustering of prokaryotic and eukaryotic samples, probably reacting to the same environmental forces (e.g., pH, conductivity). In addition, the shifts in eukaryotic and bacterioplanktonic communities generally occurred at the same time and on the same scale. Regression analyses revealed a linear relationship between an increase in Bray-Curtis dissimilarities and elapsed time. Our study shows a pronounced coupling between bacteria and eukaryotes in seasonal samplings of the three analyzed lakes. However, our temporal resolution (biweekly sampling) and data on abiotic factors were insufficient to determine if this was caused by direct biotic interactions or by reacting to the same seasonally changing environmental forces.

Responses of stream microbes to multiple anthropogenic stressors in a mesocosm study.

Stream ecosystems are affected by multiple anthropogenic stressors worldwide. Even though effects of many single stressors are comparatively well studied, the effects of multiple stressors are difficult to predict. In particular bacteria and protists, which are responsible for the majority of ecosystem respiration and element flows, are infrequently studied with respect to multiple stressors responses. We conducted a stream mesocosm experiment to characterize the responses of single and multiple stressors on microbiota. Two functionally important stream habitats, leaf litter and benthic phototrophic rock biofilms, were exposed to three stressors in a full factorial design: fine sediment deposition, increased chloride concentration (salinization) and reduced flow velocity. We analyzed the microbial composition in the two habitat types of the mesocosms using an amplicon sequencing approach. Community analysis on different taxonomic levels as well as principle component analyses (PCoAs) based on realtive abundances of operational taxonomic units (OTUs) showed treatment specific shifts in the eukaryotic biofilm community. Analysis of variance (ANOVA) revealed that Bacillariophyta responded positively salinity and sediment increase, while the relative read abundance of chlorophyte taxa decreased. The combined effects of multiple stressors were mainly antagonistic. Therefore, the community composition in multiply stressed environments resembled the composition of the unstressed control community in terms of OTU occurrence and relative abundances.

Evolution of heterotrophy in chrysophytes as reflected by comparative transcriptomics.

Shifts in the nutritional mode between phototrophy, mixotrophy and heterotrophy are a widespread phenomenon in the evolution of eukaryotic diversity. The transition between nutritional modes is particularly pronounced in chrysophytes and occurred independently several times through parallel evolution. Thus, chrysophytes provide a unique opportunity for studying the molecular basis of nutritional diversification and of the accompanying pathway reduction and degradation of plastid structures. In order to analyze the succession in switching the nutritional mode from mixotrophy to heterotrophy, we compared the transcriptome of the mixotrophic Poterioochromonas malhamensis with the transcriptomes of three obligate heterotrophic species of Ochromonadales. We used the transcriptome of P. malhamensis as a reference for plastid reduction in the heterotrophic taxa. The analyzed heterotrophic taxa were in different stages of plastid reduction. We investigated the reduction of several photosynthesis related pathways e.g. the xanthophyll cycle, the mevalonate pathway, the shikimate pathway and the tryptophan biosynthesis as well as the reduction of plastid structures and postulate a presumable succession of pathway reduction and degradation of accompanying structures.

Cryptophyta as major bacterivores in freshwater summer plankton.

Small bacterivorous eukaryotes play a cardinal role in aquatic food webs and their taxonomic classification is currently a hot topic in aquatic microbial ecology. Despite increasing interest in their diversity, core questions regarding predator-prey specificity remain largely unanswered, e.g., which heterotrophic nanoflagellates (HNFs) are the main bacterivores in freshwaters and which prokaryotes support the growth of small HNFs. To answer these questions, we fed natural communities of HNFs from Rímov reservoir (Czech Republic) with five different bacterial strains of the ubiquitous betaproteobacterial genera Polynucleobacter and Limnohabitans. We combined amplicon sequencing and catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH) targeting eukaryotic 18 S rRNA genes to track specific responses of the natural HNF community to prey amendments. While amplicon sequencing provided valuable qualitative data and a basis for designing specific probes, the number of reads was insufficient to accurately quantify certain eukaryotic groups. We also applied a double-hybridization technique that allows simultaneous phylogenetic identification of both predator and prey. Our results show that community composition of HNFs is strongly dependent upon prey type. Surprisingly, Cryptophyta were the most abundant bacterivores, although this phylum has been so far assumed to be mainly autotrophic. Moreover, the growth of a small lineage of Cryptophyta (CRY1 clade) was strongly stimulated by one Limnohabitans strain in our experiment. Thus, our study is the first report that colorless Cryptophyta are major bacterivores in summer plankton samples and can play a key role in the carbon transfer from prokaryotes to higher trophic levels.

Quantitative Proteomics Reveals Ecophysiological Effects of Light and Silver Stress on the Mixotrophic Protist Poterioochromonas malhamensis.

Aquatic environments are heavily impacted by human activities including climate warming and the introduction of xenobiotics. Due to the application of silver nanoparticles as bactericidal agent the introduction of silver into the environment strongly has increased during the past years. Silver ions affect the primary metabolism of algae, in particular photosynthesis. Mixotrophic algae are an interesting test case as they do not exclusively rely on photosynthesis which may attenuate the harmful effect of silver. In order to study the effect of silver ions on mixotrophs, cultures of the chrysophyte Poterioochromonas malhamensis were treated in a replicate design in light and darkness with silver nitrate at a sub-lethal concentration. At five time points samples were taken for the identification and quantitation of proteins by mass spectrometry. In our analysis, relative quantitative protein mass spectrometry has shown to be a useful tool for functional analyses in conjunction with transcriptome reference sequences. A total of 3,952 proteins in 63 samples were identified and quantified, mapping to 4,829 transcripts of the sequenced and assembled transcriptome. Among them, 720 and 104 proteins performing various cellular functions were differentially expressed after eight days in light versus darkness and after three days of silver treatment, respectively. Specifically pathways of the energy and primary carbon metabolism were differentially affected by light and the utilization of expensive reactions hints to an energy surplus of P. malhamensis under light conditions. The excess energy is not invested in growth, but in the synthesis of storage metabolites. The effects of silver were less explicit, observable especially in the dark treatments where the light effect could not mask coinciding but weaker effects of silver. Photosynthesis, particularly the light harvesting complexes, and several sulphur containing enzymes were affected presumably due to a direct interference with the silver ions, mainly affecting energy supply.

Trade-off between taxon diversity and functional diversity in European lake ecosystems.

Inferring ecosystem functioning and ecosystem services through inspections of the species inventory is a major aspect of ecological field studies. Ecosystem functions are often stable despite considerable species turnover. Using metatranscriptome analyses, we analyse a thus-far unparalleled freshwater data set which comprises 21 mainland European freshwater lakes from the Sierra Nevada (Spain) to the Carpathian Mountains (Romania) and from northern Germany to the Apennines (Italy) and covers an altitudinal range from 38 m above sea level (a.s.l) to 3110 m a.s.l. The dominant taxa were Chlorophyta and streptophytic algae, Ciliophora, Bacillariophyta and Chrysophyta. Metatranscriptomics provided insights into differences in community composition and into functional diversity via the relative share of taxa to the overall read abundance of distinct functional genes on the ecosystem level. The dominant metabolic pathways in terms of the fraction of expressed sequences in the cDNA libraries were affiliated with primary metabolism, specifically oxidative phosphorylation, photosynthesis and the TCA cycle. Our analyses indicate that community composition is a good first proxy for the analysis of ecosystem functions. However, differential gene regulation modifies the relative importance of taxa in distinct pathways. Whereas taxon composition varies considerably between lakes, the relative importance of distinct metabolic pathways is much more stable, indicating that ecosystem functioning is buffered against shifts in community composition through a functional redundancy of taxa.

Protistan community analysis: key findings of a large-scale molecular sampling.

Protists are perhaps the most lineage-rich of microbial lifeforms, but remain largely unknown. High-throughput sequencing technologies provide opportunities to screen whole habitats in depth and enable detailed comparisons of different habitats to measure, compare and map protistan diversity. Such comparisons are often limited by low sample numbers within single studies and a lack of standardisation between studies. Here, we analysed 232 samples from 10 sampling campaigns using a standardised PCR protocol and bioinformatics pipeline. We show that protistan community patterns are highly consistent within habitat types and geographic regions, provided that sample processing is standardised. Community profiles are only weakly affected by fluctuations of the abundances of the most abundant taxa and, therefore, provide a sound basis for habitat comparison beyond random short-term fluctuations in the community composition. Further, we provide evidence that distribution patterns are not solely resulting from random processes. Distinct habitat types and distinct taxonomic groups are dominated by taxa with distinct distribution patterns that reflect their ecology with respect to dispersal and habitat colonisation. However, there is no systematic shift of the distribution pattern with taxon abundance.

Biodiversity and new records of microfungi in the Ruhrarea (north Rhine Westfalia), Germany.

During our investigations of the microflora in NRW (Duisburg, Düsseldorf and Essen incl. the greenhouse of the Botanical Garden) in 2007 and 2008, we were able to collect and identify about 55 species on trees, bushes and ornamental plants as parasites and saprophytes. Some of these species are new for Germany or have been only rarely found until now. Most of the species belong the Ascomycotina, Basidiomycotina and Deuteromycotina for example Arthrocladiella mougeotii (Lév.) Vassilkov. on Lycium barbarum L., Caudospora taleola (Fr.) Starb on Quercus robur L., Colletotrichum coffeanum F. Noak on Coffea arabica L. (new for Germany) Colletotrichum trichellum (Fr.) Duke on Hedera helix L., Erysiphe buhrii U. Braun on Lychnis cf. coronaria (L.) Desr. (Anamorph. Oidium dianthi Jacz.), Erysiphe spec. on Acer opalus Mill (new host), Erysiphe flexuosa (Peck) U. Braun & S. Takam. on Aesculus spec. (new for Europe)), Erysiphe heraclei DC. on Tinguarra montana (Webb ex Christ ) A.Hansen & G.Kunkel, Erysiphe necator Schwein. = Uncinula necator (Schwein.) Burrill on Cissus cf. rhombifolia Vahl. (new for NRW), Erysphe trifolii Grev. on Trigonella caerulea (L.) Ser., Golovinomyces cichoracearum (DC.) V.P.Gelyuta (Oidium spec.) on Argyranthemum pinnatifidum (L.f.) R.T. Lowe (new host), Lobatopedis foliicola P.M. Kirk on Quercus robur L. (new for NRW), Lophodermium juniperinum (Fr.) de Not. on Juniperus communis L., Mamiania coryli De Not. on Corylus avellana L., Marssonina juglandis (Lib.) Magnus on Juglans regia L., Oidium hortensia Jørst on Philadelphus coronarius L., Oidium spec. on Dahlia variabilis (Willd.) Desf. (new for Germany), Oidium longipes Noordeloos & Loerak on Petunia hybrida Vilm., Oidium pedilanthi M. Yen on Pedilanthus titymaloides (L.) Poit, Oidium pedaliacearum H.D. Shin sp. nov. (= Oidium sesami H.D. Shin) on Ibicella lutea (Lindl.) van Eselt. (= Martynia lutea Lindl.), Passalora pastinacae (Sacc.) U. Braun = Pseudocercosporella pastinacae (P. Karst.) U. Braun (new for Germany) on Pastinaca sativa L., Podosphaera tridactyla (WalIr.) de Bary on Prunus laurocerasus L., Septoria cornicola Desm. on Cornus sanguinea L., Stigmina tinea (Sacc.) M.B.Ellis on Viburnum opulus L., Torula herbarum (Pers.) Link on Potentilla argentea L., etc. All species are located in the herbarium Mycotheca parva collection G.B. Feige and N. Ale-Agha.

Quantitative variability of direct chemical defense in primary and secondary leaves of lima bean (Phaseolus lunatus) and consequences for a natural herbivore.

Ontogenetic variability in chemical plant defenses against herbivores is a common phenomenon, but the effects of this variability on herbivore-plant interactions are little understood. In a previous study on lima bean (Phaseolus lunatus), we found a trade-off between cyanogenesis, a direct defense, and the release of herbivore-induced volatile organic compounds (VOCs; mainly functioning as an indirect defense). Moreover, the expression of these two defenses could change during plant ontogeny. The present study aimed at elucidating whether such ontogenetic changes in plant defense can affect herbivore-plant interactions. We quantified feeding rates of a natural insect herbivore, the Mexican bean beetle (Epilachna varivestis), on primary and secondary leaves of individual lima bean plants. These insects strongly preferred low cyanogenic primary leaves over high cyanogenic secondary leaves. Although weakly defended by cyanogenesis, lima beans' primary leaves showed protein concentrations and photosynthetic activities that did not differ significantly from secondary leaves at the time of analysis. Based on our findings, we suggest that lima beans' long-lived primary leaves function as efficient source organs, even beyond the stage of seedlings. This hypothesis may explain why primary leaves express a strong indirect defense by the release of herbivore induced VOCs.

Lateral diffusion of CO2 from shaded to illuminated leaf parts affects photosynthesis inside homobaric leaves.

Gas exchange is generally regarded to occur between the leaf interior and ambient air, i.e. in vertical (anticlinal) directions of leaf blades. However, inside homobaric leaves, gas movement occurs also in lateral directions. The aim of the present study was to ascertain whether lateral CO2 diffusion affects leaf photosynthesis when illuminated leaves are partially shaded. Measurements using gas exchange and chlorophyll fluorescence imaging techniques were performed on homobaric leaves of Vicia faba and Nicotiana tabacum or on heterobaric leaves of Glycine max and Phaseolus vulgaris. For homobaric leaves, gas exchange inside a clamp-on leaf chamber was affected by shading the leaf outside the chamber. The quantum yield of photosystem II (Phi(PSII)) was highest directly adjacent to a light/shade border (LSB). Phi(PSII) decreased in the illuminated leaf parts with distance from the LSB, while the opposite was observed for nonphotochemical quenching. These effects became most pronounced at low stomatal conductance. They were not observed in heterobaric leaves. The results suggest that plants with homobaric leaves can benefit from lateral CO2 flux, in particular when stomata are closed (e.g. under drought stress). This may enhance photosynthetic, instead of nonphotochemical, processes near LSBs in such leaves and reduce the photoinhibitory effects of excess light.