There is increasing interest in employing shotgun sequencing, rather than amplicon sequencing, to analyze microbiome samples. Typical projects may involve hundreds of samples and billions of sequencing reads. The comparison of such samples against a protein reference database generates billions of alignments and the analysis of such data is computationally challenging. To address this, we have substantially rewritten and extended our widely-used microbiome analysis tool MEGAN so as to facilitate the interactive analysis of the taxonomic and functional content of very large microbiome datasets. Other new features include a functional classifier called InterPro2GO, gene-centric read assembly, principal coordinate analysis of taxonomy and function, and support for metadata. The new program is called MEGAN Community Edition (CE) and is open source. By integrating MEGAN CE with our high-throughput DNA-to-protein alignment tool DIAMOND and by providing a new program MeganServer that allows access to metagenome analysis files hosted on a server, we provide a straightforward, yet powerful and complete pipeline for the analysis of metagenome shotgun sequences. We illustrate how to perform a full-scale computational analysis of a metagenomic sequencing project, involving 12 samples and 800 million reads, in less than three days on a single server. All source code is available here: https://github.com/danielhuson/megan-ce Microbiome sequencing projects continue to grow rapidly, both in the number of samples considered and sequencing reads collected. With MEGAN Community Edition (CE), we provide a highly efficient program for interactive analysis and comparison of such data, allowing one to explore hundreds of samples and billions of reads. While taxonomic profiling is performed based on the NCBI taxonomy, we provide a number of different functional profiling approaches such as SEED, eggNOG, KEGG, and a new InterPro2GO classification scheme. MEGAN CE also supports the use of metadata in the context of principal coordinate analysis and clustering analysis.

There are numerous computational tools for taxonomic or functional analysis of microbiome samples, optimized to run on hundreds of millions of short, high quality sequencing reads. Programs such as MEGAN allow the user to interactively navigate these large datasets. Long read sequencing technologies continue to improve and produce increasing numbers of longer reads (of varying lengths in the range of 10k-1M bps, say), but of low quality. There is an increasing interest in using long reads in microbiome sequencing, and there is a need to adapt short read tools to long read datasets. We describe a new LCA-based algorithm for taxonomic binning, and an interval-tree based algorithm for functional binning, that are explicitly designed for long reads and assembled contigs. We provide a new interactive tool for investigating the alignment of long reads against reference sequences. For taxonomic and functional binning, we propose to use LAST to compare long reads against the NCBI-nr protein reference database so as to obtain frame-shift aware alignments, and then to process the results using our new methods. All presented methods are implemented in the open source edition of MEGAN, and we refer to this new extension as MEGAN-LR (MEGAN long read). We evaluate the LAST+MEGAN-LR approach in a simulation study, and on a number of mock community datasets consisting of Nanopore reads, PacBio reads and assembled PacBio reads. We also illustrate the practical application on a Nanopore dataset that we sequenced from an anammox bio-rector community. This article was reviewed by Nicola Segata together with Moreno Zolfo, Pete James Lockhart and Serghei Mangul. This work extends the applicability of the widely-used metagenomic analysis software MEGAN to long reads. Our study suggests that the presented LAST+MEGAN-LR pipeline is sufficiently fast and accurate.

Lack of specific definitions of clinical characteristics, disease severity, and risk and preventive factors of post-COVID-19 syndrome (PCS) severely impacts research and discovery of new preventive and therapeutics drugs. This prospective multicenter cohort study was conducted from February 2020 to June 2022 in 5 countries, enrolling SARS-CoV-2 out- and in-patients followed at 3-, 6-, and 12-month from diagnosis, with assessment of clinical and biochemical features, antibody (Ab) response, Variant of Concern (VoC), and physical and mental quality of life (QoL). Outcome of interest was identification of risk and protective factors of PCS by clinical phenotype, setting, severity of disease, treatment, and vaccination status. We used SF-36 questionnaire to assess evolution in QoL index during follow-up and unsupervised machine learning algorithms (principal component analysis, PCA) to explore symptom clusters. Severity of PCS was defined by clinical phenotype and QoL. We also used generalized linear models to analyse the impact of PCS on QoL and associated risk and preventive factors. CT registration number: NCT05097677. Among 1796 patients enrolled, 1030 (57%) suffered from at least one symptom at 12-month. PCA identified 4 clinical phenotypes: chronic fatigue-like syndrome (CFs: fatigue, headache and memory loss, 757 patients, 42%), respiratory syndrome (REs: cough and dyspnoea, 502, 23%); chronic pain syndrome (CPs: arthralgia and myalgia, 399, 22%); and neurosensorial syndrome (NSs: alteration in taste and smell, 197, 11%). Determinants of clinical phenotypes were different (all comparisons p < 0.05): being female increased risk of CPs, NSs, and CFs; chronic pulmonary diseases of REs; neurological symptoms at SARS-CoV-2 diagnosis of REs, NSs, and CFs; oxygen therapy of CFs and REs; and gastrointestinal symptoms at SARS-CoV-2 diagnosis of CFs. Early treatment of SARS-CoV-2 infection with monoclonal Ab (all clinical phenotypes), corticosteroids therapy for mild/severe cases (NSs), and SARS-CoV-2 vaccination (CPs) were less likely to be associated to PCS (all comparisons p < 0.05). Highest reduction in QoL was detected in REs and CPs (43.57 and 43.86 vs 57.32 in PCS-negative controls, p < 0.001). Female sex (p < 0.001), gastrointestinal symptoms (p = 0.034) and renal complications (p = 0.002) during the acute infection were likely to increase risk of severe PCS (QoL <50). Vaccination and early treatment with monoclonal Ab reduced the risk of severe PCS (p = 0.01 and p = 0.03, respectively). Our study provides new evidence suggesting that PCS can be classified by clinical phenotypes with different impact on QoL, underlying possible different pathogenic mechanisms. We identified factors associated to each clinical phenotype and to severe PCS. These results might help in designing pathogenesis studies and in selecting high-risk patients for inclusion in therapeutic and management clinical trials. The study received funding from the Horizon 2020 ORCHESTRA project, grant 101016167; from the 10.13039/501100001826 Netherlands Organisation for Health Research and Development (ZonMw), grant 10430012010023; from 10.13039/501100001677 Inserm , REACTing (REsearch & ACtion emergING infectious diseases) consortium and the French Ministry of Health, grant PHRC 20-0424.

Short-read sequencing technologies have long been the work-horse of microbiome analysis. Continuing technological advances are making the application of long-read sequencing to metagenomic samples increasingly feasible. We demonstrate that whole bacterial chromosomes can be obtained from an enriched community, by application of MinION sequencing to a sample from an EBPR bioreactor, producing 6 Gb of sequence that assembles into multiple closed bacterial chromosomes. We provide a simple pipeline for processing such data, which includes a new approach to correcting erroneous frame-shifts. Advances in long-read sequencing technology and corresponding algorithms will allow the routine extraction of whole chromosomes from environmental samples, providing a more detailed picture of individual members of a microbiome.

Structural biology experiments and structure prediction tools have provided many high-resolution three-dimensional structures of nucleic acids. Also, molecular dynamics force field parameters have been adapted to simulating charged and flexible nucleic acid structures on microsecond time scales. Therefore, we can generate the dynamics of DNA or RNA molecules, but we still lack adequate tools for the analysis of the resulting huge amounts of data. We present MINT (Motif Identifier for Nucleic acids Trajectory) — an automatic tool for analyzing three-dimensional structures of RNA and DNA, and their full-atom molecular dynamics trajectories or other conformation sets (e.g. X-ray or nuclear magnetic resonance-derived structures). For each RNA or DNA conformation MINT determines the hydrogen bonding network resolving the base pairing patterns, identifies secondary structure motifs (helices, junctions, loops, etc.) and pseudoknots. MINT also estimates the energy of stacking and phosphate anion-base interactions. For many conformations, as in a molecular dynamics trajectory, MINT provides averages of the above structural and energetic features and their evolution. We show MINT functionality based on all-atom explicit solvent molecular dynamics trajectory of the 30S ribosomal subunit.