UnCoVar: a reproducible and scalable workflow for transparent and robust virus variant calling and lineage assignment using SARS-CoV-2 as an example.

BACKGROUND: At a global scale, the SARS-CoV-2 virus did not remain in its initial genotype for a long period of time, with the first global reports of variants of concern (VOCs) in late 2020. Subsequently, genome sequencing has become an indispensable tool for characterizing the ongoing pandemic, particularly for typing SARS-CoV-2 samples obtained from patients or environmental surveillance. For such SARS-CoV-2 typing, various in vitro and in silico workflows exist, yet to date, no systematic cross-platform validation has been reported. RESULTS: In this work, we present the first comprehensive cross-platform evaluation and validation of in silico SARS-CoV-2 typing workflows. The evaluation relies on a dataset of 54 patient-derived samples sequenced with several different in vitro approaches on all relevant state-of-the-art sequencing platforms. Moreover, we present UnCoVar, a robust, production-grade reproducible SARS-CoV-2 typing workflow that outperforms all other tested approaches in terms of precision and recall. CONCLUSIONS: In many ways, the SARS-CoV-2 pandemic has accelerated the development of techniques and analytical approaches. We believe that this can serve as a blueprint for dealing with future pandemics. Accordingly, UnCoVar is easily generalizable towards other viral pathogens and future pandemics. The fully automated workflow assembles virus genomes from patient samples, identifies existing lineages, and provides high-resolution insights into individual mutations. UnCoVar includes extensive quality control and automatically generates interactive visual reports. UnCoVar is implemented as a Snakemake workflow. The open-source code is available under a BSD 2-clause license at github.com/IKIM-Essen/uncovar.

Small-scale wastewater-based epidemiology (WBE) for infectious diseases and antibiotic resistance: A scoping review.

Wastewater analysis can serve as a source of public health information. In recent years, wastewater-based epidemiology (WBE) has emerged and proven useful for the detection of infectious diseases. However, insights from the wastewater treatment plant do not allow for the small-scale differentiation within the sewer system that is needed to analyze the target population under study in more detail. Small-scale WBE offers several advantages, but there has been no systematic overview of its application. The aim of this scoping review is to provide a comprehensive overview of the current state of knowledge on small-scale WBE for infectious diseases, including methodological considerations for its application. A systematic database search was conducted, considering only peer-reviewed articles. Data analyses included quantitative summary and qualitative narrative synthesis. Of 2130 articles, we included 278, most of which were published since 2020. The studies analyzed wastewater at the building level (n = 203), especially healthcare (n = 110) and educational facilities (n = 80), and at the neighborhood scale (n = 86). The main analytical parameters were viruses (n = 178), notably SARS-CoV-2 (n = 161), and antibiotic resistance (ABR) biomarkers (n = 99), often analyzed by polymerase chain reaction (PCR), with DNA sequencing techniques being less common. In terms of sampling techniques, active sampling dominated. The frequent lack of detailed information on the specification of selection criteria and the characterization of the small-scale sampling sites was identified as a concern. In conclusion, based on the large number of studies, we identified several methodological considerations and overarching strategic aspects for small-scale WBE. An enabling environment for small-scale WBE requires inter- and transdisciplinary knowledge sharing across countries. Promoting the adoption of small-scale WBE will benefit from a common international conceptualization of the approach, including standardized and internationally accepted terminology. In particular, the development of good WBE practices for different aspects of small-scale WBE is warranted. This includes the establishment of guidelines for a comprehensive characterization of the local sewer system and its sub-sewersheds, and transparent reporting to ensure comparability of small-scale WBE results.

Long-term mitigation of drought changes the functional potential and life-strategies of the forest soil microbiome involved in organic matter decomposition.

Climate change can alter the flow of nutrients and energy through terrestrial ecosystems. Using an inverse climate change field experiment in the central European Alps, we explored how long-term irrigation of a naturally drought-stressed pine forest altered the metabolic potential of the soil microbiome and its ability to decompose lignocellulolytic compounds as a critical ecosystem function. Drought mitigation by a decade of irrigation stimulated profound changes in the functional capacity encoded in the soil microbiome, revealing alterations in carbon and nitrogen metabolism as well as regulatory processes protecting microorganisms from starvation and desiccation. Despite the structural and functional shifts from oligotrophic to copiotrophic microbial lifestyles under irrigation and the observation that different microbial taxa were involved in the degradation of cellulose and lignin as determined by a time-series stable-isotope probing incubation experiment with 13C-labeled substrates, degradation rates of these compounds were not affected by different water availabilities. These findings provide new insights into the impact of precipitation changes on the soil microbiome and associated ecosystem functioning in a drought-prone pine forest and will help to improve our understanding of alterations in biogeochemical cycling under a changing climate.

Analyzing community wastewater in sub-sewersheds for the small-scale detection of SARS-CoV-2 variants in a German metropolitan area.

Wastewater surveillance of SARS-CoV-2 proved useful, including for identifying the local appearance of newly identified virus variants. Previous studies focused on wastewater treatment plants (WWTP) with sewersheds of several hundred thousand people or at single building level, representing only a small number of people. Both approaches may prove inadequate for small-scale intra-urban inferences for early detection of emerging or novel virus variants. Our study aims (i) to analyze SARS-CoV-2 single nucleotide variants (SNVs) in wastewater of sub-sewersheds and WWTP using whole genome sequencing in order to (ii) investigate the potential of small-scale detection of novel known SARS-CoV-2 variants of concern (VOC) within a metropolitan wastewater system. We selected three sub-sewershed sampling sites, based on estimated population- and built environment-related indicators, and the inlet of the receiving WWTP in the Ruhr region, Germany. Untreated wastewater was sampled weekly between October and December 2021, with a total of 22 samples collected. SARS-CoV-2 RNA was analyzed by RT-qPCR and whole genome sequencing. For all samples, genome sequences were obtained, while only 13 samples were positive for RT-qPCR. We identified multiple specific SARS-CoV-2 SNVs in the wastewater samples of the sub-sewersheds and the WWTP. Identified SNVs reflected the dominance of VOC Delta at the time of sampling. Interestingly, we could identify an Omicron-specific SNV in one sub-sewershed. A concurrent wastewater study sampling the same WWTP detected the VOC Omicron one week later. Our observations suggest that the small-scale approach may prove particularly useful for the detection and description of spatially confined emerging or existing virus variants circulating in populations. Future studies applying small-scale sampling strategies taking into account the specific features of the wastewater system will be useful to analyze temporal and spatial variance in more detail.

Immune responses in COVID-19 patients during breakthrough infection with SARS-CoV-2 variants Delta, Omicron-BA.1 and Omicron-BA.5.

Background: Breakthrough infections with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants are increasingly observed in vaccinated individuals. Immune responses towards SARS-CoV-2 variants, particularly Omicron-BA.5, are poorly understood. We investigated the humoral and cellular immune responses of hospitalized COVID-19 patients during Delta and Omicron infection waves. Methods: The corresponding SARS-CoV-2 variant of the respective patients were identified by whole genome sequencing. Humoral immune responses were analyzed by ELISA and a cell culture-based neutralization assay against SARS-CoV-2 D614G isolate (wildtype), Alpha, Delta (AY.43) and Omicron (BA.1 and BA.5). Cellular immunity was evaluated with an IFN-γ ELISpot assay. Results: On a cellular level, patients showed a minor IFN-γ response after stimulating PBMCs with mutated regions of SARS-CoV-2 variants. Neutralizing antibody titers against Omicron-BA.1 and especially BA.5 were strongly reduced. Double-vaccinated patients with Delta breakthrough infection showed a significantly increased neutralizing antibody response against Delta compared to double-vaccinated uninfected controls (median complete neutralization titer (NT100) 640 versus 80, p<0.05). Omicron-BA.1 infection increased neutralization titers against BA.1 in double-vaccinated patients (median NT100 of 160 in patients versus 20 in controls, p=0.07) and patients that received booster vaccination (median NT100 of 50 in patients versus 20 in controls, p=0.68). For boosted patients with BA.5 breakthrough infection, we found no enhancing effect on humoral immunity against SARS-CoV-2 variants. Conclusion: Neutralizing antibody titers against Omicron-BA.1 and especially BA.5 were strongly reduced in SARS-CoV-2 breakthrough infections. Delta and Omicron-BA.1 but not Omicron-BA.5 infections boosted the humoral immunity in double-vaccinated patients and patients with booster vaccination. Despite BA.5 breakthrough infection, those patients may still be vulnerable for reinfections with BA.5 or other newly emerging variants of concern.

Temporal disturbance of a model stream ecosystem by high microbial diversity from treated wastewater.

Microbial communities in freshwater streams play an essential role in ecosystem functioning via biogeochemical cycling. Yet, the impacts of treated wastewater influx into stream ecosystems on microbial strain diversity remain mostly unexplored. Here, we coupled full-length 16S ribosomal RNA gene Nanopore sequencing and strain-resolved metagenomics to investigate the impact of treated wastewater on a mesocosm system (AquaFlow) run with restored river water. Over 10 days, community Bray-Curtis dissimilarities between treated and control mesocosm decreased (0.57 ± 0.058 to 0.26 ± 0.046) based on ribosomal protein S3 gene clustering, finally converging to nearly identical communities. Similarly, strain-resolved metagenomics revealed a high diversity of bacteria and viruses after the introduction of treated wastewater; these microbes also decreased over time resulting in the same strain clusters in control and treatment at the end of the experiment. Specifically, 39.2% of viral strains detected in all samples were present after the introduction of treated wastewater only. Although bacteria present at low abundance in the treated wastewater introduced additional antibiotic resistance genes, signals of naturally occurring ARG-encoding organisms resembled the resistome at the endpoint. Our results suggest that the previously stressed freshwater stream and its microbial community are resilient to a substantial introduction of treated wastewater.

Identification of herbal teas and their compounds eliciting antiviral activity against SARS-CoV-2 in vitro.

BACKGROUND: The SARS-CoV-2/COVID-19 pandemic has inflicted medical and socioeconomic havoc, and despite the current availability of vaccines and broad implementation of vaccination programs, more easily accessible and cost-effective acute treatment options preventing morbidity and mortality are urgently needed. Herbal teas have historically and recurrently been applied as self-medication for prophylaxis, therapy, and symptom alleviation in diverse diseases, including those caused by respiratory viruses, and have provided sources of natural products as basis for the development of therapeutic agents. To identify affordable, ubiquitously available, and effective treatments, we tested herbs consumed worldwide as herbal teas regarding their antiviral activity against SARS-CoV-2. RESULTS: Aqueous infusions prepared by boiling leaves of the Lamiaceae perilla and sage elicit potent and sustained antiviral activity against SARS-CoV-2 when applied after infection as well as prior to infection of cells. The herbal infusions exerted in vitro antiviral effects comparable to interferon-ß and remdesivir but outperformed convalescent sera and interferon-α2 upon short-term treatment early after infection. Based on protein fractionation analyses, we identified caffeic acid, perilla aldehyde, and perillyl alcohol as antiviral compounds. Global mass spectrometry (MS) analyses performed comparatively in two different cell culture infection models revealed changes of the proteome upon treatment with herbal infusions and provided insights into the mode of action. As inferred by the MS data, induction of heme oxygenase 1 (HMOX-1) was confirmed as effector mechanism by the antiviral activity of the HMOX-1-inducing compounds sulforaphane and fraxetin. CONCLUSIONS: In conclusion, herbal teas based on perilla and sage exhibit antiviral activity against SARS-CoV-2 including variants of concern such as Alpha, Beta, Delta, and Omicron, and we identified HMOX-1 as potential therapeutic target. Given that perilla and sage have been suggested as treatment options for various diseases, our dataset may constitute a valuable resource also for future research beyond virology.

Diagnosis of Inflammatory Bowel Disease and Colorectal Cancer through Multi-View Stacked Generalization Applied on Gut Microbiome Data.

Most of the microbiome studies suggest that using ensemble models such as Random Forest results in best predictive power. In this study, we empirically evaluate a more powerful ensemble learning algorithm, multi-view stacked generalization, on pediatric inflammatory bowel disease and adult colorectal cancer patients' cohorts. We aim to check whether stacking would lead to better results compared to using a single best machine learning algorithm. Stacking achieves the best test set Average Precision (AP) on inflammatory bowel disease dataset reaching AP = 0.69, outperforming both the best base classifier (AP = 0.61) and the baseline meta learner built on top of base classifiers (AP = 0.63). On colorectal cancer dataset, the stacked classifier also outperforms (AP = 0.81) both the best base classifier (AP = 0.79) and the baseline meta learner (AP = 0.75). Stacking achieves best predictive performance on test set outperforming the best classifiers on both patient cohorts. Application of the stacking solves the issue of choosing the most appropriate machine learning algorithm by automating the model selection procedure. Clinical application of such a model is not limited to diagnosis task only, but it also can be extended to biomarker selection thanks to feature selection procedure.

Metadata harmonization-Standards are the key for a better usage of omics data for integrative microbiome analysis.

BACKGROUND: Tremendous amounts of data generated from microbiome research studies during the last decades require not only standards for sampling and preparation of omics data but also clear concepts of how the metadata is prepared to ensure re-use for integrative and interdisciplinary microbiome analysis. RESULTS: In this Commentary, we present our views on the key issues related to the current system for metadata submission in omics research, and propose the development of a global metadata system. Such a system should be easy to use, clearly structured in a hierarchical way, and should be compatible with all existing microbiome data repositories, following common standards for minimal required information and common ontology. Although minimum metadata requirements are essential for microbiome datasets, the immense technological progress requires a flexible system, which will have to be constantly improved and re-thought. While FAIR principles (Findable, Accessible, Interoperable, and Reusable) are already considered, international legal issues on genetic resource and sequence sharing provided by the Convention on Biological Diversity need more awareness and engagement of the scientific community. CONCLUSIONS: The suggested approach for metadata entries would strongly improve retrieving and re-using data as demonstrated in several representative use cases. These integrative analyses, in turn, would further advance the potential of microbiome research for novel scientific discoveries and the development of microbiome-derived products.

The MG-RAST API explorer: an on-ramp for RESTful query composition.

BACKGROUND: The MG-RAST API provides search capabilities and delivers organism and function data as well as raw or annotated sequence data via the web interface and its RESTful API. For casual users, however, RESTful APIs are hard to learn and work with. RESULTS: We created the graphical MG-RAST API explorer to help researchers more easily build and export API queries; understand the data abstractions and indices available in MG-RAST; and use the results presented in-browser for exploration, development, and debugging. CONCLUSIONS: The API explorer lowers the barrier to entry for occasional or first-time MG-RAST API users.

Federation in genomics pipelines: techniques and challenges.

Federation is a popular concept in building distributed cyberinfrastructures, whereby computational resources are provided by multiple organizations through a unified portal, decreasing the complexity of moving data back and forth among multiple organizations. Federation has been used in bioinformatics only to a limited extent, namely, federation of datastores, e.g. SBGrid Consortium for structural biology and Gene Expression Omnibus (GEO) for functional genomics. Here, we posit that it is important to federate both computational resources (CPU, GPU, FPGA, etc.) and datastores to support popular bioinformatics portals, with fast-increasing data volumes and increasing processing requirements. A prime example, and one that we discuss here, is in genomics and metagenomics. It is critical that the processing of the data be done without having to transport the data across large network distances. We exemplify our design and development through our experience with metagenomics-RAST (MG-RAST), the most popular metagenomics analysis pipeline. Currently, it is hosted completely at Argonne National Laboratory. However, through a recently started collaborative National Institutes of Health project, we are taking steps toward federating this infrastructure. Being a widely used resource, we have to move toward federation without disrupting 50 K annual users. In this article, we describe the computational tools that will be useful for federating a bioinformatics infrastructure and the open research challenges that we see in federating such infrastructures. It is hoped that our manuscript can serve to spur greater federation of bioinformatics infrastructures by showing the steps involved, and thus, allow them to scale to support larger user bases.

MG-RAST version 4-lessons learned from a decade of low-budget ultra-high-throughput metagenome analysis.

As technologies change, MG-RAST is adapting. Newly available software is being included to improve accuracy and performance. As a computational service constantly running large volume scientific workflows, MG-RAST is the right location to perform benchmarking and implement algorithmic or platform improvements, in many cases involving trade-offs between specificity, sensitivity and run-time cost. The work in [Glass EM, Dribinsky Y, Yilmaz P, et al. ISME J 2014;8:1-3] is an example; we use existing well-studied data sets as gold standards representing different environments and different technologies to evaluate any changes to the pipeline. Currently, we use well-understood data sets in MG-RAST as platform for benchmarking. The use of artificial data sets for pipeline performance optimization has not added value, as these data sets are not presenting the same challenges as real-world data sets. In addition, the MG-RAST team welcomes suggestions for improvements of the workflow. We are currently working on versions 4.02 and 4.1, both of which contain significant input from the community and our partners that will enable double barcoding, stronger inferences supported by longer-read technologies, and will increase throughput while maintaining sensitivity by using Diamond and SortMeRNA. On the technical platform side, the MG-RAST team intends to support the Common Workflow Language as a standard to specify bioinformatics workflows, both to facilitate development and efficient high-performance implementation of the community's data analysis tasks.

Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.

The metagenomic data life-cycle: standards and best practices.

Metagenomics data analyses from independent studies can only be compared if the analysis workflows are described in a harmonized way. In this overview, we have mapped the landscape of data standards available for the description of essential steps in metagenomics: (i) material sampling, (ii) material sequencing, (iii) data analysis, and (iv) data archiving and publishing. Taking examples from marine research, we summarize essential variables used to describe material sampling processes and sequencing procedures in a metagenomics experiment. These aspects of metagenomics dataset generation have been to some extent addressed by the scientific community, but greater awareness and adoption is still needed. We emphasize the lack of standards relating to reporting how metagenomics datasets are analysed and how the metagenomics data analysis outputs should be archived and published. We propose best practice as a foundation for a community standard to enable reproducibility and better sharing of metagenomics datasets, leading ultimately to greater metagenomics data reuse and repurposing.

Identification of the Core Set of Carbon-Associated Genes in a Bioenergy Grassland Soil.

Despite the central role of soil microbial communities in global carbon (C) cycling, little is known about soil microbial community structure and even less about their metabolic pathways. Efforts to characterize soil communities often focus on identifying differences in gene content across environmental gradients, but an alternative question is what genes are similar in soils. These genes may indicate critical species or potential functions that are required in all soils. Here we identified the "core" set of C cycling sequences widely present in multiple soil metagenomes from a fertilized prairie (FP). Of 226,887 sequences associated with known enzymes involved in the synthesis, metabolism, and transport of carbohydrates, 843 were identified to be consistently prevalent across four replicate soil metagenomes. This core metagenome was functionally and taxonomically diverse, representing five enzyme classes and 99 enzyme families within the CAZy database. Though it only comprised 0.4% of all CAZy-associated genes identified in FP metagenomes, the core was found to be comprised of functions similar to those within cumulative soils. The FP CAZy-associated core sequences were present in multiple publicly available soil metagenomes and most similar to soils sharing geographic proximity. In soil ecosystems, where high diversity remains a key challenge for metagenomic investigations, these core genes represent a subset of critical functions necessary for carbohydrate metabolism, which can be targeted to evaluate important C fluxes in these and other similar soils.