Cascabel: A Scalable and Versatile Amplicon Sequence Data Analysis Pipeline Delivering Reproducible and Documented Results.

Marker gene sequencing of the rRNA operon (16S, 18S, ITS) or cytochrome c oxidase I (CO1) is a popular means to assess microbial communities of the environment, microbiomes associated with plants and animals, as well as communities of multicellular organisms via environmental DNA sequencing. Since this technique is based on sequencing a single gene, or even only parts of a single gene rather than the entire genome, the number of reads needed per sample to assess the microbial community structure is lower than that required for metagenome sequencing. This makes marker gene sequencing affordable to nearly any laboratory. Despite the relative ease and cost-efficiency of data generation, analyzing the resulting sequence data requires computational skills that may go beyond the standard repertoire of a current molecular biologist/ecologist. We have developed Cascabel, a scalable, flexible, and easy-to-use amplicon sequence data analysis pipeline, which uses Snakemake and a combination of existing and newly developed solutions for its computational steps. Cascabel takes the raw data as input and delivers a table of operational taxonomic units (OTUs) or Amplicon Sequence Variants (ASVs) in BIOM and text format and representative sequences. Cascabel is a highly versatile software that allows users to customize several steps of the pipeline, such as selecting from a set of OTU clustering methods or performing ASV analysis. In addition, we designed Cascabel to run in any linux/unix computing environment from desktop computers to computing servers making use of parallel processing if possible. The analyses and results are fully reproducible and documented in an HTML and optional pdf report. Cascabel is freely available at Github: https://github.com/AlejandroAb/CASCABEL.

The absence of intact polar lipid-derived GDGTs in marine waters dominated by Marine Group II: Implications for lipid biosynthesis in Archaea.

The marine pelagic archaeal community is dominated by three major groups, the marine group I (MGI) Thaumarchaeota, and the marine groups II and III (MGII and MGIII) Euryarchaeota. Studies of both MGI cultures and the environment have shown that the MGI core membrane lipids are predominantly composed of glycerol dibiphytanyl glycerol tetraether (GDGT) lipids and the diether lipid archaeol. However, there are no cultured representatives of MGII and III archaea and, therefore, both their membrane lipid composition and potential contribution to the marine archaeal lipid pool remain unknown. Here, we show that GDGTs present in suspended particulate matter of the (sub)surface waters of the North Atlantic Ocean and the coastal North Sea are derived from MGI archaea, and that MGII archaea do not significantly contribute to the pool of GDGTs and archaeol. This implies, in contrast to previous suggestions, that their lipids do not affect the widely used sea surface temperature proxy TEX86. These findings also indicate that MGII archaea are not able to produce any known archaeal lipids, implying that our understanding of the evolution of membrane lipid biosynthesis in Archaea is far from complete.