Search

Abstract Background The possibility of recovering metagenome-assembled genomes (MAGs) from sequence reads allows for further insights into microbial communities and their members, possibly even analyzing such sequences with tools designed for single-isolate genomes. As result quality depends on sequence quality, performance of tools for single-isolate genomes on MAGs should be tested beforehand. Bioinformatics can be leveraged to quickly create varied synthetic test sets with known composition for this purpose. Results We present MAGICIAN, a flexible, user-friendly pipeline for the simulation of MAGs. MAGICIAN combines a synthetic metagenome simulator with a metagenomic assembly and binning pipeline to simulate MAGs based on user-supplied input genomes, allowing users to test performance of tools on MAGs while having a ground truth to compare results to. Using MAGICIAN, we found that even very slight (1%) changes in depth of coverage can drastically affect whether a genome can be recovered. We also demonstrate the use of simulated MAGs by evaluating the suitability of such genomes obtained with MAGICIAN’s current default pipeline for analysis with the antimicrobial resistance gene identification tool ResFinder. Conclusions Using MAGICIAN, it is possible to simulate MAGs which, while generally high in quality, reflect issues encountered with real-world data, thus providing realistic best-case data. Evaluating the results of ResFinder analysis of these genomes revealed a risk for plausible-looking false positives, which underlines the need for pipeline validation so that researchers are aware of the potential issues when interpreting real-world data. Furthermore, the effects of fluctuations in depth of coverage on genome recovery in our simulated “random sequencing” warrant further investigation and indicate random subsampling of reads may affect discovery of more genomes.

Immunocompromised patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can have a longer duration of viral shedding and persistence of symptoms. The optimal treatment strategy for these patients remains to be established. This case describes a male in his late sixties with follicular lymphoma and persistent symptoms of infection with SARS-CoV-2 variant BA.2 who was treated with remdesivir five times over a period of six months. The clinical effect of remdesivir treatment decreased over time, and further viral sequencing revealed the emergence of mutations across the SARS-CoV-2 genome. Due to the lack of other treatment options, the patient was treated with a combination of remdesivir and molnupiravir for 10 days, and epcoritamab was discontinued, which led to the cessation of symptoms. This case illustrates the risk of a diminished effect of remdesivir with prolonged use and the need for treatment guidelines for immunocompromised patients with persistent COVID-19.

ABSTRACT Microbes produce a plethora of secondary (or specialized) metabolites that, although not essential for primary metabolism, benefit them to survive in the environment, communicate, and influence cell differentiation. Biosynthetic gene clusters (BGCs), responsible for the production of these secondary metabolites, are readily identifiable on bacterial genome sequences. Understanding the phylogeny and distribution of BGCs helps us to predict the natural product synthesis ability of new isolates. Here, we examined 310 genomes from the Bacillus subtilis group, determined the inter- and intraspecies patterns of absence/presence for all BGCs, and assigned them to defined gene cluster families (GCFs). This allowed us to establish patterns in the distribution of both known and unknown products. Further, we analyzed variations in the BGC structures of particular families encoding natural products, such as plipastatin, fengycin, iturin, mycosubtilin, and bacillomycin. Our detailed analysis revealed multiple GCFs that are species or clade specific and a few others that are scattered within or between species, which will guide exploration of the chemodiversity within the B. subtilis group. Surprisingly, we discovered that partial deletion of BGCs and frameshift mutations in selected biosynthetic genes are conserved within phylogenetically related isolates, although isolated from around the globe. Our results highlight the importance of detailed genomic analysis of BGCs and the remarkable phylogenetically conserved erosion of secondary metabolite biosynthetic potential in the B. subtilis group. IMPORTANCE Members of the B. subtilis species complex are commonly recognized producers of secondary metabolites, among those, the production of antifungals, which makes them promising biocontrol strains. While there are studies examining the distribution of well-known secondary metabolites in Bacilli, intraspecies clade-specific distribution has not been systematically reported for the B. subtilis group. Here, we report the complete biosynthetic potential within the B. subtilis group to explore the distribution of the biosynthetic gene clusters and to reveal an exhaustive phylogenetic conservation of secondary metabolite production within Bacillus that supports the chemodiversity within this species complex. We identify that certain gene clusters acquired deletions of genes and particular frameshift mutations, rendering them inactive for secondary metabolite biosynthesis, a conserved genetic trait within phylogenetically conserved clades of certain species. The overview guides the assignment of the secondary metabolite production potential of newly isolated Bacillus strains based on genome sequence and phylogenetic relatedness.