MicrobeAnnotator: a user-friendly, comprehensive functional annotation pipeline for microbial genomes

Unraveling radiation resistance strategies in two bacterial strains from the high background radiation area of Chavara-Neendakara: A comprehensive whole genome analysis

This paper reports the results of gamma irradiation experiments and whole genome sequencing (WGS) performed on vegetative cells of two radiation resistant bacterial strains, Metabacillus halosaccharovorans (VITHBRA001) and Bacillus paralicheniformis (VITHBRA024) (D10 values 2.32 kGy and 1.42 kGy, respectively), inhabiting the top-ranking high background radiation area (HBRA) of Chavara-Neendakara placer deposit (Kerala, India). The present investigation has been carried out in the context that information on strategies of bacteria having mid-range resistance for gamma radiation is inadequate. WGS, annotation, COG and KEGG analyses and manual curation of genes helped us address the possible pathways involved in the major domains of radiation resistance, involving recombination repair, base excision repair, nucleotide excision repair and mismatch repair, and the antioxidant genes, which the candidate could activate to survive under ionizing radiation. Additionally, with the help of these data, we could compare the candidate strains with that of the extremely radiation resistant model bacterium Deinococccus radiodurans, so as to find the commonalities existing in their strategies of resistance on the one hand, and also the rationale behind the difference in D10, on the other. Genomic analysis of VITHBRA001 and VITHBRA024 has further helped us ascertain the difference in capability of radiation resistance between the two strains. Significantly, the genes such as uvsE (NER), frnE (protein protection), ppk1 and ppx (non-enzymatic metabolite production) and those for carotenoid biosynthesis, are endogenous to VITHBRA001, but absent in VITHBRA024, which could explain the former's better radiation resistance. Further, this is the first-time study performed on any bacterial population inhabiting an HBRA. This study also brings forward the two species whose radiation resistance has not been reported thus far, and add to the knowledge on radiation resistant capabilities of the phylum Firmicutes which are abundantly observed in extreme environment.

Eupransor demetentiae gen. nov., sp. nov., a novel fructophilic lactic acid bacterium from bumble bees

Strain LMG 33000T was isolated from a Bombus lapidarius gut sample. It shared the highest percentage 16S rRNA sequence identity, average amino acid identity, and amino acid identity of conserved genes with Convivina intestini LMG 28291T (95.86 %, 69.9 and 76.2 %, respectively), and the highest percentage OrthoANIu value with Fructobacillus fructosus DSM 20349T (71.4 %). Phylogenomic analyses by means of 107 or 120 conserved genes consistently revealed Convivina as nearest neighbour genus. The draft genome of strain LMG 33000T was 1.44 Mbp in size and had a DNA G+C content of 46.1 mol%. Genomic and physiological analyses revealed that strain LMG 33000T was a typical obligately fructophilic lactic acid bacterium that lacked the adhE and aldh genes and that did not produce ethanol during glucose or fructose metabolism. In contrast, Convivina species have the adhE and aldh genes in their genomes and produced ethanol from glucose and fructose metabolism, which is typical for heterofermentative lactic acid bacteria. Moreover, strain LMG 33000T exhibited catalase activity, an unusual characteristic among lactic acid bacteria, that is not shared with Convivina species. Given its position in the phylogenomic trees, and the difference in genomic percentage G+C content and in physiological and metabolic characteristics between strain LMG 33000T and Convivina species, we considered it most appropriate to classify strain LMG 33000T into a novel genus and species within the Lactobacillaceae family for which we propose the name Eupransor demetentiae gen. nov., sp. nov., with LMG 33000T (=CECT 30958T) as the type strain.

Microbial adaptation to spaceflight is correlated with bacteriophage-encoded functions

Evidence from the International Space Station suggests microbial populations are rapidly adapting to the spacecraft environment; however, the mechanism of this adaptation is not understood. Bacteriophages are prolific mediators of bacterial adaptation on Earth. Here we survey 245 genomes sequenced from bacterial strains isolated on the International Space Station for dormant (lysogenic) bacteriophages. Our analysis indicates phage-associated genes are significantly different between spaceflight strains and their terrestrial counterparts. In addition, we identify 283 complete prophages, those that could initiate bacterial lysis and infect additional hosts, of which 21% are novel. These prophage regions encode functions that correlate with increased persistence in extreme environments, such as spaceflight, to include antimicrobial resistance and virulence, DNA damage repair, and dormancy. Our results correlate microbial adaptation in spaceflight to bacteriophage-encoded functions that may impact human health in spaceflight.

A comparative genomic analysis of Fructobacillus evanidus sp. nov. from bumble bees

The increase in studies on bee microbiomes is prompted by concerns over global pollinator declines. Bumble bees host core and non-core microbiota which may contribute to increased lifetime fitness. The presence of Fructobacillus in the gut microbiomes of bumble bee workers, or the replacement of core symbionts with Fructobacillus bacteria, has been considered a marker of dysbiosis. A phylogenomic analysis and functional genomic characterization of the genomes of 21 Fructobacillus isolates from bumble bees demonstrated that they represented four species, i.e. Fructobacillus cardui, Fructobacillus fructosus, Fructobacillus tropaeoli, and the novel species Fructobacillus evanidus sp. nov. Our results confirmed and substantiated the presence of two phylogenetically and functionally distinct Fructobacillus species clades that differ in genome size, percentage G + C content, the number of coding DNA sequences and metabolic characteristics. Clade 1 and clade 2 species differed in amino acid and, to a lesser extent, in carbohydrate metabolism, with F. evanidus and F. tropaeoli genomes featuring a higher number of complete metabolic pathways. While Fructobacillus genomes encoded genes that allow adhesion, biofilm formation, antibacterial activity and detoxification, other bacteria isolated from the bumble bee gut appeared better equipped to co-exist with the bumble bee host. The isolation and identification of multiple Fructobacillus species from several bumble bee gut samples in the present study also argued against a specific partnership between Fructobacillus species and their bumble bee hosts.

MerCat2: a versatile k-mer counter and diversity estimator for database-independent property analysis obtained from omics data

Motivation

MerCat2 ("Mer-Catenate2") is a versatile, parallel, scalable and modular property software package for robustly analyzing features in omics data. Using massively parallel sequencing raw reads, assembled contigs, and protein sequences from any platform as input, MerCat2 performs k-mer counting of any length k, resulting in feature abundance counts tables, quality control reports, protein feature metrics, and graphical representation (i.e. principal component analysis (PCA)).

Results

MerCat2 allows for direct analysis of data properties in a database-independent manner that initializes all data, which other profilers and assembly-based methods cannot perform. MerCat2 represents an integrated tool to illuminate omics data within a sample for rapid cross-examination and comparisons.

Availability and implementation

MerCat2 is written in Python and distributed under a BSD-3 license. The source code of MerCat2 is freely available at https://github.com/raw-lab/mercat2. MerCat2 is compatible with Python 3 on Mac OS X and Linux. MerCat2 can also be easily installed using bioconda: mamba create -n mercat2 -c conda-forge -c bioconda mercat2.

Unveiling the Microbial Realm with VEBA 2.0: A modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic, and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible platform that bridges the gap between genomics and biotechnological solutions.

MetaCerberus: distributed highly parallelized HMM-based processing for robust functional annotation across the tree of life

MOTIVATION

MetaCerberus is a massively parallel, fast, low memory, scalable annotation tool for inference gene function across genomes to metacommunities. MetaCerberus provides an elusive HMM/HMMER-based tool at a rapid scale with low memory. It offers scalable gene elucidation to major public databases, including KEGG (KO), COGs, CAZy, FOAM, and specific databases for viruses, including VOGs and PHROGs, from single genomes to metacommunities.

RESULTS

MetaCerberus is 1.3× as fast on a single node than eggNOG-mapper v2 on 5× less memory using an exclusively HMM/HMMER mode. In a direct comparison, MetaCerberus provides better annotation of viruses, phages, and archaeal viruses than DRAM, Prokka, or InterProScan. MetaCerberus annotates more KOs across domains when compared to DRAM, with a 186× smaller database, and with 63× less memory. MetaCerberus is fully integrated for automatic analysis of statistics and pathways using differential statistic tools (i.e. DESeq2 and edgeR), pathway enrichment (GAGE R), and pathview R. MetaCerberus provides a novel tool for unlocking the biosphere across the tree of life at scale.

AVAILABILITY AND IMPLEMENTATION

MetaCerberus is written in Python and distributed under a BSD-3 license. The source code of MetaCerberus is freely available at https://github.com/raw-lab/metacerberus compatible with Python 3 and works on both Mac OS X and Linux. MetaCerberus can also be easily installed using bioconda: mamba create -n metacerberus -c bioconda -c conda-forge metacerberus.

Investigating the viral ecology and contribution to the microbial ecology in full-scale mesophilic anaerobic digesters

In an attempt to assess the diversity of viruses and their potential to modulate the metabolism of functional microorganisms in anaerobic digesters, we collected digestate from three mesophilic anaerobic digesters in full-scale wastewater treatment plants treating real municipal wastewater. The reads were analyzed using bioinformatics algorithms to elucidate viral diversity, identify their potential role in modulating the metabolism of functional microorganisms, and provide essential genomic information for the potential use of virus-mediated treatment in controlling the anaerobic digester microbiome. We found that Siphoviridae was the dominant family in mesophilic anaerobic digesters, followed by Myoviridae and Podoviridae. Lysogeny was prevalent in mesophilic anaerobic digesters as the majority of metagenome-assembled genomes contained at least one viral genome within them. One virus within the genome of an acetoclastic methanogen (Methanothrix soehngenii) was observed with a gene (fwdE) acquired via lateral transfer from hydrogenotrophic methanogens. The virus-mediated acquisition of fwdE gene enables possibility of mixotrophic methanogenesis in Methanothrix soehngenii. This evidence highlighted that lysogeny provides fitness advantage to methanogens in anaerobic digesters by adding flexibility to changing substrates. Similarly, we found auxiliary metabolic genes, such as cellulase and alpha glucosidase, of bacterial origin responsible for sludge hydrolysis in viruses. Additionally, we discovered novel viral genomes and provided genomic information on viruses infecting acidogenic, acetogenic, and pathogenic bacteria that can potentially be used for virus-mediated treatment to deal with the souring problem in anaerobic digesters and remove pathogens from biosolids before land application. Collectively, our study provides a genome-level understanding of virome in conjunction with the microbiome in anaerobic digesters that can be used to optimize the anaerobic digestion process for efficient biogas generation.

Pseudogenomic insights into the evolution of Mycobacterium ulcerans

BACKGROUND

Buruli ulcer (BU) disease, caused by Mycobacterium ulcerans (MU), and characterized by necrotic ulcers is still a health problem in Africa and Australia. The genome of the bacterium has several pseudogenes due to recent evolutionary events and environmental pressures. Pseudogenes are genetic elements regarded as nonessential in bacteria, however, they are less studied due to limited available tools to provide understanding of their evolution and roles in MU pathogenicity.

RESULTS

This study developed a bioinformatic pipeline to profile the pseudogenomes of sequenced MU clinical isolates from different countries. One hundred and seventy-two MU genomes analyzed revealed that pseudogenomes of African strains corresponded to the two African lineages 1 and 2. Pseudogenomes were lineage and location specific and African lineage 1 was further divided into A and B. Lineage 2 had less relaxation in positive selection than lineage 1 which may signify different evolutionary points. Based on the Gil-Latorre model, African MU strains may be in the latter stages of evolutionary adaption and are adapting to an environment rich in metabolic resources with a lower temperature and decreased UV radiation. The environment fosters oxidative metabolism and MU may be less reliant on some secondary metabolites. In-house pseudogenomes from Ghana and Cote d'Ivoire were different from other African strains, however, they were identified as African strains.

CONCLUSION

Our bioinformatic pipeline provides pseudogenomic insights to complement other whole genome analyses, providing a better view of the evolution of the genome of MU and suggest an adaptation model which is important in understanding transmission. MU pseudogene profiles vary based on lineage and country, and an apparent reduction in insertion sequences used for the detection of MU which may adversely affect the sensitivity of diagnosis.

Adaptation of Anaerobic Digestion Microbial Communities to High Ammonium Levels: Insights from Strain-Resolved Metagenomics

Ammonia release from proteinaceous feedstocks represents the main inhibitor of the anaerobic digestion (AD) process, which can result in a decreased biomethane yield or even complete failure of the process. The present study focused on the adaptation of mesophilic AD communities to a stepwise increase in the concentration of ammonium chloride in synthetic medium with casein used as the carbon source. An adaptation process occurring over more than 20 months allowed batch reactors to reach up to 20 g of NH4+ N/L without collapsing in acidification nor ceasing methane production. To decipher the microbial dynamics occurring during the adaptation and determine the genes mostly exposed to selective pressure, a combination of biochemical and metagenomics analyses was performed, reconstructing the strains of key species and tracking them over time. Subsequently, the adaptive metabolic mechanisms were delineated by following the single nucleotide variants (SNVs) characterizing the strains and prioritizing the associated genes according to their function. An in-depth exploration of the archaeon Methanoculleus bourgensis vb3066 and the putative syntrophic acetate-oxidizing bacteria Acetomicrobium sp. ma133 identified positively selected SNVs on genes involved in stress adaptation. The intraspecies diversity with multiple coexisting strains in a temporal succession pattern allows us to detect the presence of an additional level of diversity within the microbial community beyond the species level.

Decoding Microbial Responses to Ammonia Shock Loads in Biogas Reactors through Metagenomics and Metatranscriptomics

The presence of elevated ammonia levels is widely recognized as a significant contributor to process inhibition in biogas production, posing a common challenge for biogas plant operators. The present study employed a combination of biochemical, genome-centric metagenomic and metatranscriptomic data to investigate the response of the biogas microbiome to two shock loads induced by single pulses of elevated ammonia concentrations (i.e., 1.5 g NH4+/LR and 5 g NH4+/LR). The analysis revealed a microbial community of high complexity consisting of 364 Metagenome Assembled Genomes (MAGs). The hydrogenotrophic pathway was the primary route for methane production during the entire experiment, confirming its efficiency even at high ammonia concentrations. Additionally, metatranscriptomic analysis uncovered a metabolic shift in the methanogens Methanothrix sp. MA6 and Methanosarcina flavescens MX5, which switched their metabolism from the acetoclastic to the CO2 reduction route during the second shock. Furthermore, multiple genes associated with mechanisms for maintaining osmotic balance in the cell were upregulated, emphasizing the critical role of osmoprotection in the rapid response to the presence of ammonia. Finally, this study offers insights into the transcriptional response of an anaerobic digestion community, specifically focusing on the mechanisms involved in recovering from ammonia-induced stress.

Dissemination of antibiotic resistance genes through soil-plant-earthworm continuum in the food production environment

Treated municipal wastewater (TMW) can provide a reliable source of irrigation water for crops, which is especially important in arid areas where water resources are limited or prone to drought. Nonetheless, TMW may contain residual antibiotics, potentially exposing the crops to these substances. The goal of this study was to investigate the dissemination of antibiotics resistance genes (ARGs) in the soil-plant-earthworm continuum after irrigation of spinach and radish plants with TMW containing trimethoprim, sulfamethoxazole, and sulfapyridine in a greenhouse experiment, followed by feeding of earthworms with harvested plant materials. Our results showed that antibiotic resistance genes (ARGs) were enriched in the soil-plant-earthworm microbiomes irrigated with TMW and TMW spiked with higher concentrations of antibiotics. The number of ARGs and antibiotic-resistant bacteria (ARB) enrichment varied with plant type, with spinach harboring a significantly higher amount of ARGs and ARB compared to radish. Our data showed that bulk and rhizosphere soils of spinach and radish plants irrigated with MilliQ water, TMW, TMW10, or TMW100 had significant differences in bacterial community (p < 0.001), ARG (p < 0.001), and virulence factor gene (VFG) (p < 0.001) diversities. The abundance of ARGs significantly decreased from bulk soil to rhizosphere to phyllosphere and endosphere. Using metagenome assembled genomes (MAGs), we recovered many bacterial MAGs and a near complete genome (>90 %) of bacterial MAG of genus Leclercia adecarboxylata B from the fecal microbiome of earthworm that was fed harvested radish tubers and spinach leaves grown on TMW10 irrigated waters, and this bacterium has been shown to be an emerging pathogen causing infection in immunocompromised patients that may lead to health complications and death. Therefore, crops irrigated with TMW containing residual antibiotics and ARGs may lead to increased incidences of enrichment of ARB in the soil-plant-earthworm continuum.

Host genetic regulation of human gut microbial structural variation

Although the impact of host genetics on gut microbial diversity and the abundance of specific taxa is well established1-6, little is known about how host genetics regulates the genetic diversity of gut microorganisms. Here we conducted a meta-analysis of associations between human genetic variation and gut microbial structural variation in 9,015 individuals from four Dutch cohorts. Strikingly, the presence rate of a structural variation segment in Faecalibacterium prausnitzii that harbours an N-acetylgalactosamine (GalNAc) utilization gene cluster is higher in individuals who secrete the type A oligosaccharide antigen terminating in GalNAc, a feature that is jointly determined by human ABO and FUT2 genotypes, and we could replicate this association in a Tanzanian cohort. In vitro experiments demonstrated that GalNAc can be used as the sole carbohydrate source for F. prausnitzii strains that carry the GalNAc-metabolizing pathway. Further in silico and in vitro studies demonstrated that other ABO-associated species can also utilize GalNAc, particularly Collinsella aerofaciens. The GalNAc utilization genes are also associated with the host's cardiometabolic health, particularly in individuals with mucosal A-antigen. Together, the findings of our study demonstrate that genetic associations across the human genome and bacterial metagenome can provide functional insights into the reciprocal host-microbiome relationship.

Functional and genomic characterization of a novel probiotic Lactobacillus johnsonii KD1 against shrimp WSSV infection

White Spot syndrome virus (WSSV) causes rapid shrimp mortality and production loss worldwide. This study demonstrates potential use of Lactobacillus johnsonii KD1 as an anti-WSSV agent for post larva shrimp cultivation and explores some potential mechanisms behind the anti-WSSV properties. Treatment of Penaeus vannamei shrimps with L. johnsonii KD1 prior to oral challenge with WSSV-infected tissues showed a significantly reduced mortality. In addition, WSSV copy numbers were not detected and shrimp immune genes were upregulated. Genomic analysis of L. johnsonii KD1 based on Illumina and Nanopore platforms revealed a 1.87 Mb chromosome and one 15.4 Kb plasmid. Only one antimicrobial resistance gene (ermB) in the chromosome was identified. Phylogenetic analysis comparing L. johnsonii KD1 to other L. johnsonii isolates revealed that L. johnsonii KD1 is closely related to L. johnsonii GHZ10a isolated from wild pigs. Interestingly, L. johnsonii KD1 contains isolate-specific genes such as genes involved in a type I restriction-modification system and CAZymes belonging to the GT8 family. Furthermore, genes coding for probiotic survival and potential antimicrobial/anti-viral metabolites such as a homolog of the bacteriocin helveticin-J were found. Protein-protein docking modelling suggests the helveticin-J homolog may be able to block VP28-PmRab7 interactions and interrupt WSSV infection.

Why are so many enteric pathogen infections asymptomatic? Pathogen and gut microbiome characteristics associated with diarrhea symptoms and carriage of diarrheagenic E. coli in northern Ecuador

A high proportion of enteric infections, including those caused by diarrheagenic Escherichia coli (DEC), are asymptomatic for diarrhea. The factors responsible for the development of diarrhea symptoms, or lack thereof, remain unclear. Here, we used DEC isolate genome and whole stool microbiome data from a case-control study of diarrhea in Ecuador to examine factors associated with diarrhea symptoms accompanying DEC carriage. We investigated i) pathogen abundance, ii) gut microbiome characteristics, and iii) strain-level pathogen characteristics from DEC infections with diarrhea symptoms (symptomatic infections) and without diarrhea symptoms (asymptomatic infections). We also included data from individuals with and without diarrhea who were not infected with DEC (uninfected cases and controls). i) E. coli relative abundance in the gut microbiome was highly variable, but higher on-average in individuals with symptomatic compared to asymptomatic DEC infections. Similarly, the number and relative abundances of virulence genes in the gut were higher in symptomatic than asymptomatic DEC infections. ii) Measures of microbiome diversity were similar regardless of diarrhea symptoms or DEC carriage. Proteobacterial families that have been described as pathobionts were enriched in symptomatic infections and uninfected cases, whereas potentially beneficial taxa, including the Bacteroidaceae and Bifidobacteriaceae, were more abundant in individuals without diarrhea. An analysis of high-level gene functions recovered in metagenomes revealed that genes that were differentially abundant by diarrhea and DEC infection status were more abundant in symptomatic than asymptomatic DEC infections. iii) DEC isolates from symptomatic versus asymptomatic individuals showed no significant differences in virulence or accessory gene content, and there was no phylogenetic signal associated with diarrhea symptoms. Together, these data suggest signals that distinguish symptomatic from asymptomatic DEC infections. In particular, the abundance of E. coli, the virulence gene content of the gut microbiome, and the taxa present in the gut microbiome have an apparent role.

Applications of genome-scale metabolic models to investigate microbial metabolic adaptations in response to genetic or environmental perturbations

Environmental perturbations are encountered by microorganisms regularly and will require metabolic adaptations to ensure an organism can survive in the newly presenting conditions. In order to study the mechanisms of metabolic adaptation in such conditions, various experimental and computational approaches have been used. Genome-scale metabolic models (GEMs) are one of the most powerful approaches to study metabolism, providing a platform to study the systems level adaptations of an organism to different environments which could otherwise be infeasible experimentally. In this review, we are describing the application of GEMs in understanding how microbes reprogram their metabolic system as a result of environmental variation. In particular, we provide the details of metabolic model reconstruction approaches, various algorithms and tools for model simulation, consequences of genetic perturbations, integration of '-omics' datasets for creating context-specific models and their application in studying metabolic adaptation due to the change in environmental conditions.

Comparative genomic analysis of Periweissella and the characterization of novel motile species

The genus Periweissella was proposed as a novel genus in the Lactobacillaceae in 2022. However, the phylogenetic relationship between Periweissella and other heterofermentative lactobacilli, and the genetic and physiological properties of this genus remain unclear. This study aimed to determine the phylogenetic relationship between Periweissella and the two closest genera, Weissella and Furfurilactobacillus, by the phylogenetic analysis and calculation of (core gene) pairwise average amino acid identity. Targeted genomic analysis showed that fructose bisphosphate aldolase was only present in the genome of Pw. cryptocerci. Mannitol dehydrogenase was found in genomes of Pw. beninensis, Pw. fabaria, and Pw. fabalis. Untargeted genomic analysis identified the presence of flagellar genes in Periweissella but not in other closely related genera. Phenotypes related to carbohydrate fermentation and motility matched the genotypes. Motility genes were organized in a single operon and the proteins shared a high amino acid similarity in the genus Periweissella. The relatively low similarity of motility operons between Periweissella and other motile lactobacilli indicated the acquisition of motility by the ancestral species. Our findings facilitate the phylogenetic, genetic, and phenotypic understanding of the genus Periweissella.ImportanceThe genus Periweissella is a heterofermentative genus in the Lactobacillaceae which includes predominantly isolates from cocoa fermentations in tropical climates. Despite the relevance of the genus in food fermentations, genetic and physiological properties of the genus are poorly characterized and genome sequences became available only after 2020. This study characterized strains of the genus by functional genomic analysis, and by determination of metabolic and physiological traits. Phylogenetic analysis revealed that Periweissella is the evolutionary link between rod-shaped heterofermentative lactobacilli and the coccoid Leuconostoc clade with the genera Weissella and Furfurilactobacillus as closest relatives. Periweissella is the only heterofermentative genus in the Lactobacillaceae which comprises predominantly motile strains. The genomic, physiological, and metabolic characterization of Periweissella may facilitate the potential use of strains of the genus as starter culture in traditional or novel food fermentations.

Stringent Response of Cyanobacteria and Other Bacterioplankton during Different Stages of a Harmful Cyanobacterial Bloom

We conducted a field study to investigate the role of stringent response in cyanobacteria and coexisting bacterioplankton during nutrient-deprived periods at various stages of bloom in a freshwater lake (Utah Lake) for the first time. Using metagenomics and metatranscriptomics analyses, we examined the cyanobacterial ecology and expression of important functional genes related to stringent response, N and P metabolism, and regulation. Our findings mark a significant advancement in understanding the mechanisms by which toxic cyanobacteria survive and proliferate during nitrogen (N) and phosphorus (P) limitations. We successfully identified and analyzed the metagenome-assembled genomes (MAGs) of the dominant bloom-forming cyanobacteria, namely, Dolichospermum circinale, Aphanizomenon flos-aquae UKL13-PB, Planktothrix agardhii, and Microcystis aeruginosa. By mapping RNA-seq data to the coding sequences of the MAGs, we observed that these four prevalent cyanobacteria species activated multiple functions to adapt to the depletion of inorganic nutrients. During and after the blooms, the four dominant cyanobacteria species expressed high levels of transcripts related to toxin production, such as microcystins (mcy), anatoxins (ana), and cylindrospermopsins (cyr). Additionally, genes associated with polyphosphate (poly-P) storage and the stringent response alarmone (p)ppGpp synthesis/hydrolysis, including ppk, relA, and spoT, were highly activated in both cyanobacteria and bacterioplankton. Under N deficiency, the main N pathways shifted from denitrification and dissimilatory nitrate reduction in bacterioplankton toward N2-fixing and assimilatory nitrate reduction in certain cyanobacteria with a corresponding shift in the community composition. P deprivation triggered a stringent response mediated by spoT-dependent (p)ppGpp accumulation and activation of the Pho regulon in both cyanobacteria and bacterioplankton, facilitating inorganic and organic P uptake. The dominant cyanobacterial MAGs exhibited the presence of multiple alkaline phosphatase (APase) transcripts (e.g., phoA in Dolichospermum, phoX in Planktothrix, and Microcystis), suggesting their ability to synthesize and release APase enzymes to convert ambient organic P into bioavailable forms. Conversely, transcripts associated with bacterioplankton-dominated pathways like denitrification were low and did not align with the occurrence of intense cyanoHABs. The strong correlations observed among N, P, stringent response metabolisms and the succession of blooms caused by dominant cyanobacterial species provide evidence that the stringent response, induced by nutrient limitation, may activate unique N and P functions in toxin-producing cyanobacteria, thereby sustaining cyanoHABs.

A Comprehensive Bioinformatics Resource Guide for Genome-Based Antimicrobial Resistance Studies

The use of high-throughput sequencing technologies and bioinformatic tools has greatly transformed microbial genome research. With the help of sophisticated computational tools, it has become easier to perform whole genome assembly, identify and compare different species based on their genomes, and predict the presence of genes responsible for proteins, antimicrobial resistance, and toxins. These bioinformatics resources are likely to continuously improve in quality, become more user-friendly to analyze the multiple genomic data, efficient in generating information and translating it into meaningful knowledge, and enhance our understanding of the genetic mechanism of AMR. In this manuscript, we provide an essential guide for selecting the popular resources for microbial research, such as genome assembly and annotation, antibiotic resistance gene profiling, identification of virulence factors, and drug interaction studies. In addition, we discuss the best practices in computer-oriented microbial genome research, emerging trends in microbial genomic data analysis, integration of multi-omics data, the appropriate use of machine-learning algorithms, and open-source bioinformatics resources for genome data analytics.

Scalable and versatile container-based pipelines for de novo genome assembly and bacterial annotation

Background: Advancements in DNA sequencing technology have transformed the field of bacterial genomics, allowing for faster and more cost effective chromosome level assemblies compared to a decade ago. However, transforming raw reads into a complete genome model is a significant computational challenge due to the varying quality and quantity of data obtained from different sequencing instruments, as well as intrinsic characteristics of the genome and desired analyses. To address this issue, we have developed a set of container-based pipelines using Nextflow, offering both common workflows for inexperienced users and high levels of customization for experienced ones. Their processing strategies are adaptable based on the sequencing data type, and their modularity enables the incorporation of new components to address the community's evolving needs. Methods: These pipelines consist of three parts: quality control, de novo genome assembly, and bacterial genome annotation. In particular, the genome annotation pipeline provides a comprehensive overview of the genome, including standard gene prediction and functional inference, as well as predictions relevant to clinical applications such as virulence and resistance gene annotation, secondary metabolite detection, prophage and plasmid prediction, and more. Results: The annotation results are presented in reports, genome browsers, and a web-based application that enables users to explore and interact with the genome annotation results. Conclusions: Overall, our user-friendly pipelines offer a seamless integration of computational tools to facilitate routine bacterial genomics research. The effectiveness of these is illustrated by examining the sequencing data of a clinical sample of Klebsiella pneumoniae.

Brytella acorum gen. nov., sp. nov., a novel acetic acid bacterium from sour beverages

Polyphasic taxonomic and comparative genomic analyses revealed that a series of lambic beer isolates including strain LMG 32668T and the kombucha isolate LMG 32879 represent a novel species among the acetic acid bacteria, with Acidomonas methanolica as the nearest phylogenomic neighbor with a valid name. Overall genomic relatedness indices and phylogenomic and physiological analyses revealed that this novel species was best classified in a novel genus for which we propose the name Brytella acorum gen. nov., sp. nov., with LMG 32668T (=CECT 30723T) as the type strain. The B. acorum genomes encode a complete but modified tricarboxylic acid cycle, and complete pentose phosphate, pyruvate oxidation and gluconeogenesis pathways. The absence of 6-phosphofructokinase which rendered the glycolysis pathway non-functional, and an energy metabolism that included both aerobic respiration and oxidative fermentation are typical metabolic characteristics of acetic acid bacteria. Neither genome encodes nitrogen fixation or nitrate reduction genes, but both genomes encode genes for the biosynthesis of a broad range of amino acids. Antibiotic resistance genes or virulence factors are absent.

Predictive modeling of antibiotic eradication therapy success for new-onset Pseudomonas aeruginosa pulmonary infections in children with cystic fibrosis

Chronic Pseudomonas aeruginosa (Pa) lung infections are the leading cause of mortality among cystic fibrosis (CF) patients; therefore, the eradication of new-onset Pa lung infections is an important therapeutic goal that can have long-term health benefits. The use of early antibiotic eradication therapy (AET) has been shown to clear the majority of new-onset Pa infections, and it is hoped that identifying the underlying basis for AET failure will further improve treatment outcomes. Here we generated machine learning models to predict AET outcomes based on pathogen genomic data. We used a nested cross validation design, population structure control, and recursive feature selection to improve model performance and showed that incorporating population structure control was crucial for improving model interpretation and generalizability. Our best model, controlling for population structure and using only 30 recursively selected features, had an area under the curve of 0.87 for a holdout test dataset. The top-ranked features were generally associated with motility, adhesion, and biofilm formation.

Arcobacteraceae comparative genome analysis demonstrates genome heterogeneity and reduction in species isolated from animals and associated with human illness

The Arcobacteraceae family groups Gram-negative bacterial species previously included in the family Campylobacteraceae. These species of which some are considered foodborne pathogens, have been isolated from different environmental niches and hosts. They have been isolated from various types of foods, though predominantly from food of animal origin, as well as from stool of humans with enteritis. Their different abilities to survive in different hosts and environments suggest an evolutionary pressure with consequent variation in their genome content. Moreover, their different physiological and genomic characteristics led to the recent proposal to create new genera within this family, which is however criticized due to the lack of discriminatory features and biological and clinical relevance. Aims of the present study were to assess the Arcobacteraceae pangenome, and to characterize existing similarities and differences in 20 validly described species. For this, analysis has been conducted on the genomes of the corresponding type strains obtained by Illumina sequencing, applying several bioinformatic tools. Results of the present study do not support the proposed division into different genera and revealed the presence of pangenome partitions with numbers comparable to other Gram-negative bacteria genera, such as Campylobacter. Different gene class compositions in animal and human-associated species are present, including a higher percentage of virulence-related gene classes such as cell motility genes. The adaptation to environmental and/or host conditions of some species was identified by the presence of specific genes. Furthermore, a division into pathogenic and non-pathogenic species is suggested, which can support future research on food safety and public health.

StandEnA: a customizable workflow for standardized annotation and generating a presence-absence matrix of proteins

Motivation

Several genome annotation tools standardize annotation outputs for comparability. During standardization, these tools do not allow user-friendly customization of annotation databases; limiting their flexibility and applicability in downstream analysis.

Results

StandEnA is a user-friendly command-line tool for Linux that facilitates the generation of custom databases by retrieving protein sequences from multiple databases. Directed by a user-defined list of standard names, StandEnA retrieves synonyms to search for corresponding sequences in a set of public databases. Custom databases are used in prokaryotic genome annotation to generate standardized presence-absence matrices and reference files containing standard database identifiers. To showcase StandEnA, we applied it to six metagenome-assembled genomes to analyze three different pathways.

Availability and implementation

StandEnA is an open-source software available at https://github.com/mdsufz/StandEnA.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Potential reservoirs of antimicrobial resistance in livestock waste and treated wastewater that can be disseminated to agricultural land

Livestock manure, dairy lagoon effluent, and treated wastewater are known reservoirs of antibiotic resistance genes (ARGs), antibiotic-resistant bacteria (ARB), and virulence factor genes (VFGs), and their application to agricultural farmland could be a serious public health threat. However, their dissemination to agricultural lands and impact on important geochemical pathways such as the nitrogen (N) cycle have not been jointly explored. In this study, shotgun metagenomic sequencing and analyses were performed to examine the diversity and composition of microbial communities, ARGs, VFGs, and N cycling genes in different livestock manure/lagoon and treated wastewater collected from concentrated animal feeding operations (CAFOs) and a municipal wastewater treatment plant along the west coast of the United States. Multivariate analysis showed that diversity indices of bacterial taxa from the different microbiomes were not significantly different based on InvSimpson (P = 0.05), but differences in ARG mechanisms were observed between swine manure and other microbiome sources. Comparative resistome profiling showed that ARGs in microbiome samples belonged to four core resistance classes: aminoglycosides (40-55 %), tetracyclines (30-45 %), beta-lactam-resistance (20-35 %), macrolides (18-30 %), and >50 % of the VFGs that the 24 microbiomes harbored were phyletically affiliated with two bacteria, Bacteroidetes fragilis and Enterobacter aerogenes. Network analysis based on Spearman correlation showed co-occurrence patterns between several genes such as transporter-gene and regulator, efflux pump and involved-in-polymyxin- resistance, aminoglycoside, beta-lactam, and macrolide with VFGs and bacterial taxa such as Firmicutes, Candidatus Themoplasmatota, Actinobacteria, and Bacteroidetes. Metabolic reconstruction of metagenome-assembled genome (MAGs) analysis showed that the most prevalent drug resistance mechanisms were associated with carbapenem resistance, multidrug resistance (MDR), and efflux pump. Bacteroidales was the main taxa involved in dissimilatory nitrate reduction (DNRA) in dairy lagoon effluent. This study demonstrates that the dissemination of waste from these sources can increase the spread of ARGs, ARB, and VFGs into agricultural lands, negatively impacting both soil and human health.

MEGAnnotator2: a pipeline for the assembly and annotation of microbial genomes

The reconstruction of microbial genome sequences by bioinformatic pipelines and the consequent functional annotation of their genes' repertoire are fundamental activities aiming at unveiling their biological mechanisms, such as metabolism, virulence factors, and antimicrobial resistances. Here, we describe the development of the MEGAnnotator2 pipeline able to manage all next-generation sequencing methodologies producing short- and long-read DNA sequences. Starting from raw sequencing data, the updated pipeline can manage multiple analyses leading to the assembly of high-quality genome sequences and the functional classification of their genetic repertoire, providing the user with a useful report constituting features and statistics related to the microbial genome. The updated pipeline is fully automated from the installation to the delivery of the output, thus requiring minimal bioinformatics knowledge to be executed.

A phylogenomic and comparative genomic analysis of Commensalibacter, a versatile insect symbiont

BACKGROUND

To understand mechanisms of adaptation and plasticity of pollinators and other insects a better understanding of diversity and function of their key symbionts is required. Commensalibacter is a genus of acetic acid bacterial symbionts in the gut of honey bees and other insect species, yet little information is available on the diversity and function of Commensalibacter bacteria. In the present study, whole-genome sequences of 12 Commensalibacter isolates from bumble bees, butterflies, Asian hornets and rowan berries were determined, and publicly available genome assemblies of 14 Commensalibacter strains were used in a phylogenomic and comparative genomic analysis.

RESULTS

The phylogenomic analysis revealed that the 26 Commensalibacter isolates represented four species, i.e. Commensalibacter intestini and three novel species for which we propose the names Commensalibacter melissae sp. nov., Commensalibacter communis sp. nov. and Commensalibacter papalotli sp. nov. Comparative genomic analysis revealed that the four Commensalibacter species had similar genetic pathways for central metabolism characterized by a complete tricarboxylic acid cycle and pentose phosphate pathway, but their genomes differed in size, G + C content, amino acid metabolism and carbohydrate-utilizing enzymes. The reduced genome size, the large number of species-specific gene clusters, and the small number of gene clusters shared between C. melissae and other Commensalibacter species suggested a unique evolutionary process in C. melissae, the Western honey bee symbiont.

CONCLUSION

The genus Commensalibacter is a widely distributed insect symbiont that consists of multiple species, each contributing in a species specific manner to the physiology of the holobiont host.

The Skin Microbiome: Current Landscape and Future Opportunities

Our skin is the largest organ of the body, serving as an important barrier against the harsh extrinsic environment. Alongside preventing desiccation, chemical damage and hypothermia, this barrier protects the body from invading pathogens through a sophisticated innate immune response and co-adapted consortium of commensal microorganisms, collectively termed the microbiota. These microorganisms inhabit distinct biogeographical regions dictated by skin physiology. Thus, it follows that perturbations to normal skin homeostasis, as occurs with ageing, diabetes and skin disease, can cause microbial dysbiosis and increase infection risk. In this review, we discuss emerging concepts in skin microbiome research, highlighting pertinent links between skin ageing, the microbiome and cutaneous repair. Moreover, we address gaps in current knowledge and highlight key areas requiring further exploration. Future advances in this field could revolutionise the way we treat microbial dysbiosis associated with skin ageing and other pathologies.

A nontuberculous mycobacterium could solve the mystery of the lady from the Franciscan church in Basel, Switzerland

BACKGROUND

In 1975, the mummified body of a female has been found in the Franciscan church in Basel, Switzerland. Molecular and genealogic analyses unveiled her identity as Anna Catharina Bischoff (ACB), a member of the upper class of post-reformed Basel, who died at the age of 68 years, in 1787. The reason behind her death is still a mystery, especially that toxicological analyses revealed high levels of mercury, a common treatment against infections at that time, in different body organs. The computed tomography (CT) and histological analysis showed bone lesions in the femurs, the rib cage, and the skull, which refers to a potential syphilis case.

RESULTS

Although we could not detect any molecular signs of the syphilis-causing pathogen Treponema pallidum subsp. pallidum, we realized high prevalence of a nontuberculous mycobacterium (NTM) species in brain tissue sample. The genome analysis of this NTM displayed richness of virulence genes and toxins, and similarity to other infectious NTM, known to infect immunocompromised patients. In addition, it displayed potential resistance to mercury compounds, which might indicate a selective advantage against the applied treatment. This suggests that ACB might have suffered from an atypical mycobacteriosis during her life, which could explain the mummy's bone lesion and high mercury concentrations.

CONCLUSIONS

The study of this mummy exemplifies the importance of employing differential diagnostic approaches in paleopathological analysis, by combining classical anthropological, radiological, histological, and toxicological observations with molecular analysis. It represents a proof-of-concept for the discovery of not-yet-described ancient pathogens in well-preserved specimens, using de novo metagenomic assembly.

Long-Read Sequencing Improves Recovery of Picoeukaryotic Genomes and Zooplankton Marker Genes from Marine Metagenomes

Long-read sequencing offers the potential to improve metagenome assemblies and provide more robust assessments of microbial community composition and function than short-read sequencing. We applied Pacific Biosciences (PacBio) CCS (circular consensus sequencing) HiFi shotgun sequencing to 14 marine water column samples and compared the results with those for short-read metagenomes from the corresponding environmental DNA samples. We found that long-read metagenomes varied widely in quality and biological information. The community compositions of the corresponding long- and short-read metagenomes were frequently dissimilar, suggesting higher stochasticity and/or bias associated with PacBio sequencing. Long reads provided few improvements to the assembly qualities, gene annotations, and prokaryotic metagenome-assembled genome (MAG) binning results. However, only long reads produced high-quality eukaryotic MAGs and contigs containing complete zooplankton marker gene sequences. These results suggest that high-quality long-read metagenomes can improve marine community composition analyses and provide important insight into eukaryotic phyto- and zooplankton genetics, but the benefits may be outweighed by the inconsistent data quality. IMPORTANCE Ocean microbes provide critical ecosystem services, but most remain uncultivated. Their communities can be studied through shotgun metagenomic sequencing and bioinformatic analyses, including binning draft microbial genomes. However, most sequencing to date has been done using short-read technology, which rarely yields genome sequences of key microbes like SAR11. Long-read sequencing can improve metagenome assemblies but is hampered by technological shortcomings and high costs. In this study, we compared long- and short-read sequencing of marine metagenomes. We found a wide range of long-read metagenome qualities and minimal improvements to microbiome analyses. However, long reads generated draft genomes of eukaryotic algal species and provided full-length marker gene sequences of zooplankton species, including krill and copepods. These results suggest that long-read sequencing can provide greater genetic insight into the wide diversity of eukaryotic phyto- and zooplankton that interact as part of and with the marine microbiome.

Metagenomic views on taxonomic and functional profiles of the Himalayan Tsomgo cold lake and unveiling its deterzome potential

Cold habitat is considered a potential source for detergent industry enzymes. This study aims at the metagenomic investigation of Tsomgo lake for taxonomic and functional annotation, unveiling the deterzome potential of the residing microbiota at this site. The present investigation revealed molecular profiling of microbial community structure and functional potential of the high-altitude Tsomgo lake samples of two different temperatures, harvested during March and August. Bacteria were found to be the most dominant phyla, with traces of genomic pieces of evidence belonging to archaea, viruses, and eukaryotes. Proteobacteria and Actinobacteria were noted to be the most abundant bacterial phyla in the cold lake. In-depth metagenomic investigation of the cold aquatic habitat revealed novel genes encoding detergent enzymes, amylase, protease, and lipase. Further, metagenome-assembled genomes (MAGs) belonging to the psychrophilic bacterium, Arthrobacter alpinus, were constructed from the metagenomic data. The annotation depicted the presence of detergent enzymes and genes for low-temperature adaptation in Arthrobacter alpinus. Psychrophilic microbial isolates were screened for lipase, protease, and amylase activities to further strengthen the metagenomic findings. A novel strain of Acinetobacter sp. was identified with the dual enzymatic activity of protease and amylase. The bacterial isolates exhibited hydrolyzing activity at low temperatures. This metagenomic study divulged novel genomic resources for detergent industry enzymes, and the bacterial isolates secreting cold-active amylase, lipase, and protease enzymes. The findings manifest that Tsomgo lake is a potential bioresource of cold-active enzymes, vital for various industrial applications.

Effect of pressure on biomethanation process and spatial stratification of microbial communities in trickle bed reactors under decreasing gas retention time

The current study investigated the effect of elevating gas pressure on biomethanation in trickle-bed reactors (TBRs). The increased pressure led to successful biomethanation (CH₄ > 90 %) at a gas retention time (GRT) of 21 min, due to the improved transfer rates of H₂ and CO₂. On the contrary, the non-pressurized TBR performance was reduced at GRTs shorter than 40 min. Metagenomic analysis revealed that the microbial populations collected from the lower and middle parts of the reactor under the same GRT were more homogeneous compared with those developed in the upper layer. Comparison with previous experiments suggest that microbial stratification is mainly driven by the nutrient provision strategy. Methanobacterium species was the most dominant methanogen and it was mainly associated with the bottom and middle parts of TBRs. Overall, the increased pressure did not affect markedly the microbial composition, while the GRT was the most important parameter shaping the microbiomes.

Dissecting the role of the human microbiome in COVID-19 via metagenome-assembled genomes

Coronavirus disease 2019 (COVID-19), primarily a respiratory disease caused by infection with Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), is often accompanied by gastrointestinal symptoms. However, little is known about the relation between the human microbiome and COVID-19, largely due to the fact that most previous studies fail to provide high taxonomic resolution to identify microbes that likely interact with SARS-CoV-2 infection. Here we used whole-metagenome shotgun sequencing data together with assembly and binning strategies to reconstruct metagenome-assembled genomes (MAGs) from 514 COVID-19 related nasopharyngeal and fecal samples in six independent cohorts. We reconstructed a total of 11,584 medium-and high-quality microbial MAGs and obtained 5403 non-redundant MAGs (nrMAGs) with strain-level resolution. We found that there is a significant reduction of strain richness for many species in the gut microbiome of COVID-19 patients. The gut microbiome signatures can accurately distinguish COVID-19 cases from healthy controls and predict the progression of COVID-19. Moreover, we identified a set of nrMAGs with a putative causal role in the clinical manifestations of COVID-19 and revealed their functional pathways that potentially interact with SARS-CoV-2 infection. Finally, we demonstrated that the main findings of our study can be largely validated in three independent cohorts. The presented results highlight the importance of incorporating the human gut microbiome in our understanding of SARS-CoV-2 infection and disease progression.

Here, using whole-metagenome shotgun sequencing data from patients with COVID-19 and controls, the authors reconstruct 11,584 microbial metagenome-assembled genomes (MAGs) including 5,403 non-redundant MAGs, revealing microbiota and metabolic pathways associations with SARS-CoV-2 infection at strain-level resolution.

Selection of Anabaena sp. PCC 7938 as a Cyanobacterium Model for Biological ISRU on Mars

Crewed missions to Mars are expected to take place in the coming decades. After short-term stays, a permanent presence will be desirable to enable a wealth of scientific discoveries. This will require providing crews with life-support consumables in amounts that are too large to be imported from Earth. Part of these consumables could be produced on site with bioprocesses, but the feedstock should not have to be imported. A solution under consideration lies in using diazotrophic, rock-weathering cyanobacteria as primary producers: fed with materials naturally available on site, they would provide the nutrients required by other organisms. This concept has recently gained momentum but progress is slowed by a lack of consistency across contributing teams, and notably of a shared model organism. With the hope to address this issue, we present the work performed to select our current model. We started with preselected strains from the Nostocaceae family. After sequencing the genome of Anabaena sp. PCC 7938-the only one not yet available-we compared the strains' genomic data to determine their relatedness and provide insights into their physiology. We then assessed and compared relevant features: chiefly, their abilities to utilize nutrients from Martian regolith, their resistance to perchlorates (toxic compounds present in the regolith), and their suitability as feedstock for secondary producers (here a heterotrophic bacterium and a higher plant). This led to the selection of Anabaena sp. PCC 7938, which we propose as a model cyanobacterium for the development of bioprocesses based on Mars's natural resources. IMPORTANCE The sustainability of crewed missions to Mars could be increased by biotechnologies which are connected to resources available on site via primary producers: diazotrophic, rock-leaching cyanobacteria. Indeed, this could greatly reduce the mass of payloads to be imported from Earth. The concept is gaining momentum but progress is hindered by a lack of consistency across research teams. We consequently describe the selection process that led to the choice of our model strain, demonstrate its relevance to the field, and propose it as a shared model organism. We expect this contribution to support the development of cyanobacterium-based biotechnologies on Mars.

Comparative Genomic Analysis of Antarctic Pseudomonas Isolates with 2,4,6-Trinitrotoluene Transformation Capabilities Reveals Their Unique Features for Xenobiotics Degradation

The nitroaromatic explosive 2,4,6-trinitrotoluene (TNT) is a highly toxic and persistent environmental pollutant. Since physicochemical methods for remediation are poorly effective, the use of microorganisms has gained interest as an alternative to restore TNT-contaminated sites. We previously demonstrated the high TNT-transforming capability of three novel Pseudomonas spp. isolated from Deception Island, Antarctica, which exceeded that of the well-characterized TNT-degrading bacterium Pseudomonas putida KT2440. In this study, a comparative genomic analysis was performed to search for the metabolic functions encoded in the genomes of these isolates that might explain their TNT-transforming phenotype, and also to look for differences with 21 other selected pseudomonads, including xenobiotics-degrading species. Comparative analysis of xenobiotic degradation pathways revealed that our isolates have the highest abundance of key enzymes related to the degradation of fluorobenzoate, TNT, and bisphenol A. Further comparisons considering only TNT-transforming pseudomonads revealed the presence of unique genes in these isolates that would likely participate directly in TNT-transformation, and others involved in the β-ketoadipate pathway for aromatic compound degradation. Lastly, the phylogenomic analysis suggested that these Antarctic isolates likely represent novel species of the genus Pseudomonas, which emphasizes their relevance as potential agents for the bioremediation of TNT and other xenobiotics.

Quinone Moieties Link the Microbial Respiration of Natural Organic Matter to the Chemical Reduction of Diverse Nitroaromatic Compounds

Insensitive munitions compounds (IMCs) are emerging nitroaromatic contaminants developed by the military as safer-to-handle alternatives to conventional explosives. Biotransformation of nitroaromatics via microbial respiration has only been reported for a limited number of substrates. Important soil microorganisms can respire natural organic matter (NOM) by reducing its quinone moieties to hydroquinones. Thus, we investigated the NOM respiration combined with the abiotic reduction of nitroaromatics by the hydroquinones formed. First, we established nitroaromatic concentration ranges that were nontoxic to the quinone respiration. Then, an enrichment culture dominated by Geobacter anodireducens could indirectly reduce a broad array of nitroaromatics by first respiring NOM components or the NOM surrogate anthraquinone-2,6-disulfonate (AQDS). Without quinones, no nitroaromatic tested was reduced except for the IMC 3-nitro-1,2,4-triazol-5-one (NTO). Thus, the quinone respiration expanded the spectrum of nitroaromatics susceptible to transformation. The system functioned with very low quinone concentrations because NOM was recycled by the nitroaromatic reduction. A metatranscriptomic analysis demonstrated that the microorganisms obtained energy from quinone or NTO reduction since respiratory genes were upregulated when AQDS or NTO was the electron acceptor. The results indicated microbial NOM respiration sustained by the nitroaromatic-dependent cycling of quinones. This process can be applied as a nitroaromatic remediation strategy, provided that a quinone pool is available for microorganisms.

Toward shotgun metagenomic approaches for microbial source tracking sewage spills based on laboratory mesocosms

Little is known about the genomic diversity of the microbial communities associated with raw municipal wastewater (sewage), including whether microbial populations specific to sewage exist and how such populations could be used to improve source attribution and apportioning in contaminated waters. Herein, we used the influent of three wastewater treatment plants in Atlanta, Georgia (USA) to perturb laboratory freshwater mesocosms, simulating sewage contamination events, and followed these mesocosms with shotgun metagenomics over a 7-day observational period. We describe 15 abundant non-redundant bacterial metagenome-assembled genomes (MAGs) ubiquitous within all sewage inocula yet absent from the unperturbed freshwater control at our analytical limit of detection. Tracking the dynamics of the populations represented by these MAGs revealed varied decay kinetics, depending on (inferred) phenotypes, e.g., anaerobes decayed faster than aerobes under the well-aerated incubation conditions. Notably, a portion of these populations showed decay patterns similar to those of common markers, Enterococcus and HF183. Despite the apparent decay of these populations, the abundance of β-lactamase encoding genes remained high throughout incubation relative to the control. Lastly, we constructed genomic libraries representing several different fecal sources and outline a bioinformatic approach which leverages these libraries for identifying and apportioning contamination signal among multiple probable sources using shotgun metagenomic data.

Reconstruction of Archaeal Genomes from Short-Read Metagenomes

As the majority of biological diversity remains unexplored and uncultured, investigating it requires culture-independent approaches. Archaea in particular suffer from a multitude of issues that make their culturing problematic, from them being frequently members of the rare biosphere, to low growth rates, to them thriving under very specific and often extreme environmental and community conditions that are difficult to replicate. OMICs techniques are state of the art approaches that allow direct high-throughput investigations of environmental samples at all levels from nucleic acids to proteins, lipids, and secondary metabolites. Metagenomics, as the foundation for other OMICs techniques, facilitates the identification and functional characterization of the microbial community members and can be combined with other methods to provide insights into the microbial activities, both on the RNA and protein levels. In this chapter, we provide a step-by-step workflow for the recovery of archaeal genomes from metagenomes, starting from raw short-read sequences. This workflow can be applied to recover bacterial genomes as well.

Considerations for Constructing a Protein Sequence Database for Metaproteomics

Mass spectrometry-based metaproteomics has emerged as a prominent technique for interrogating the functions of specific organisms in microbial communities, in addition to total community function. Identifying proteins by mass spectrometry requires matching mass spectra of fragmented peptide ions to a database of protein sequences corresponding to the proteins in the sample. This sequence database determines which protein sequences can be identified from the measurement, and as such the taxonomic and functional information that can be inferred from a metaproteomics measurement. Thus, the construction of the protein sequence database directly impacts the outcome of any metaproteomics study. Several factors, such as source of sequence information and database curation, need to be considered during database construction to maximize accurate protein identifications traceable to the species of origin. In this review, we provide an overview of existing strategies for database construction and the relevant studies that have sought to test and validate these strategies. Based on this review of the literature and our experience we provide a decision tree and best practices for choosing and implementing database construction strategies.

Small and Equipped: the Rich Repertoire of Antibiotic Resistance Genes in Candidate Phyla Radiation Genomes

Microbes belonging to Candidate Phyla Radiation (CPR) have joined the tree of life as a new branch, thanks to the intensive application of metagenomics and sequencing technologies. CPR have been eventually identified by 16S rRNA analysis, and they represent more than 26% of microbial diversity. Despite their ultrasmall size, reduced genome, and metabolic pathways which mainly depend on exosymbiotic or exoparasitic relationships with the bacterial host, CPR microbes were found to be abundant in almost all environments. They can be considered survivors in highly competitive circumstances within microbial communities. However, their defense mechanisms and phenotypic characteristic remain poorly explored. Here, we conducted a thorough in silico analysis on 4,062 CPR genomes to search for antibiotic resistance (AR)-like enzymes using BLASTp and functional domain predictions against an exhaustive consensus AR database and conserved domain database (CDD), respectively. Our findings showed that a rich reservoir of divergent AR-like genes (n = 30,545 hits, mean = 7.5 hits/genome [0 to 41]) were distributed across the 13 CPR superphyla. These AR-like genes encode 89 different enzymes that are associated with 14 different chemical classes of antimicrobials. Most hits found (93.6%) were linked to glycopeptide, beta-lactam, macrolide-lincosamide-streptogramin (MLS), tetracycline, and aminoglycoside resistance. Moreover, two AR profiles were discerned for the Microgenomates group and "Candidatus Parcubacteria," which were distinct between them and differed from all other CPR superphyla. CPR cells seem to be active players during microbial competitive interactions; they are well equipped for microbial combat in different habitats, which ensures their natural survival and continued existence. IMPORTANCE To our knowledge, this study is one of the few studies that characterize the defense systems in the CPR group and describes the first repertoire of antibiotic resistance (AR) genes. The use of a BLAST approach with lenient criteria followed by a careful examination of the functional domains has yielded a variety of enzymes that mainly give three different mechanisms of action of resistance. Our genome analysis showed the existence of a rich reservoir of CPR resistome, which is associated with different antibiotic families. Moreover, this analysis revealed the hidden face of the reduced-genome CPR, particularly their weaponry with AR genes. These data suggest that CPR are competitive players in the microbial war, and they can be distinguished by specific AR profiles.

eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale

Even though automated functional annotation of genes represents a fundamental step in most genomic and metagenomic workflows, it remains challenging at large scales. Here, we describe a major upgrade to eggNOG-mapper, a tool for functional annotation based on precomputed orthology assignments, now optimized for vast (meta)genomic data sets. Improvements in version 2 include a full update of both the genomes and functional databases to those from eggNOG v5, as well as several efficiency enhancements and new features. Most notably, eggNOG-mapper v2 now allows for: 1) de novo gene prediction from raw contigs, 2) built-in pairwise orthology prediction, 3) fast protein domain discovery, and 4) automated GFF decoration. eggNOG-mapper v2 is available as a standalone tool or as an online service at http://eggnog-mapper.embl.de.

eggNOG-mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale

Even though automated functional annotation of genes represents a fundamental step in most genomic and metagenomic workflows, it remains challenging at large scales. Here, we describe a major upgrade to eggNOG-mapper, a tool for functional annotation based on precomputed orthology assignments, now optimized for vast (meta)genomic data sets. Improvements in version 2 include a full update of both the genomes and functional databases to those from eggNOG v5, as well as several efficiency enhancements and new features. Most notably, eggNOG-mapper v2 now allows for: 1) de novo gene prediction from raw contigs, 2) built-in pairwise orthology prediction, 3) fast protein domain discovery, and 4) automated GFF decoration. eggNOG-mapper v2 is available as a standalone tool or as an online service at http://eggnog-mapper.embl.de.

Comparative pangenomic analyses and biotechnological potential of cocoa-related Acetobacter senegalensis strains

Acetobacter senegalensis belongs to the group of acetic acid bacteria (AAB) that present potential biotechnological applications, for production of D-gluconate, cellulose and acetic acid. AAB can overcome heat and acid stresses by using strategies involving the overexpression of heat-shock proteins and enzymes from the complex pyrroquinoline-ADH, besides alcohol dehydrogenases (ADH). Nonetheless, the isolation of A. senegalensis and other AAB from food may be challenging due to presence of viable but non-culturable (VBNC) cells and due to uncertainties about nutritional requirements. To contribute for a better understanding of the ecology of AAB, this paper reports on the pangenome analysis of five strains of A. senegalensis recently isolated from a Brazilian spontaneous cocoa fermentation. The results showed biosynthetic clusters exclusively found in some cocoa-related AAB, such as those related to terpene pathways, which are important for flavour development. Genes related to oxidative stress were conserved in all the genomes, with multiple clusters. Moreover, there were genes coding for ADH and putative ABC transporters distributed in core, shell and cloud genomes, while chaperonin-encoding genes were present only in the core and soft-core genomes. Regarding quorum sensing, a response regulator gene was in the shell genome, and the gene encoding for acyl-homoserine lactone efflux protein was in the soft-core genome. There were quorum quenching-related genes, mainly encoding for lactonases, but also for acylases. Moreover, A. senegalensis did not have determinants of virulence or antibiotic resistance, which are good traits for strains intended to be applied in food fermentation.

FA-nf: A Functional Annotation Pipeline for Proteins from Non-Model Organisms Implemented in Nextflow

Functional annotation allows adding biologically relevant information to predicted features in genomic sequences, and it is, therefore, an important procedure of any de novo genome sequencing project. It is also useful for proofreading and improving gene structural annotation. Here, we introduce FA-nf, a pipeline implemented in Nextflow, a versatile computational workflow management engine. The pipeline integrates different annotation approaches, such as NCBI BLAST+, DIAMOND, InterProScan, and KEGG. It starts from a protein sequence FASTA file and, optionally, a structural annotation file in GFF format, and produces several files, such as GO assignments, output summaries of the abovementioned programs and final annotation reports. The pipeline can be broken easily into smaller processes for the purpose of parallelization and easily deployed in a Linux computational environment, thanks to software containerization, thus helping to ensure full reproducibility.

Viruses and Their Interactions With Bacteria and Archaea of Hypersaline Great Salt Lake

Viruses play vital biogeochemical and ecological roles by (a) expressing auxiliary metabolic genes during infection, (b) enhancing the lateral transfer of host genes, and (c) inducing host mortality. Even in harsh and extreme environments, viruses are major players in carbon and nutrient recycling from organic matter. However, there is much that we do not yet understand about viruses and the processes mediated by them in the extreme environments such as hypersaline habitats. The Great Salt Lake (GSL) in Utah, United States is a hypersaline ecosystem where the biogeochemical role of viruses is poorly understood. This study elucidates the diversity of viruses and describes virus–host interactions in GSL sediments along a salinity gradient. The GSL sediment virosphere consisted of Haloviruses (32.07 ± 19.33%) and members of families Siphoviridae (39.12 ± 19.8%), Myoviridae (13.7 ± 6.6%), and Podoviridae (5.43 ± 0.64%). Our results demonstrate that salinity alongside the concentration of organic carbon and inorganic nutrients (nitrogen and phosphorus) governs the viral, bacteria, and archaeal diversity in this habitat. Computational host predictions for the GSL viruses revealed a wide host range with a dominance of viruses that infect Proteobacteria, Actinobacteria, and Firmicutes. Identification of auxiliary metabolic genes for photosynthesis (psbA), carbon fixation (rbcL, cbbL), formaldehyde assimilation (SHMT), and nitric oxide reduction (NorQ) shed light on the roles played by GSL viruses in biogeochemical cycles of global relevance.

Accurate annotation of protein coding sequences with IDTAXA

The observed diversity of protein coding sequences continues to increase far more rapidly than knowledge of their functions, making classification algorithms essential for assigning a function to proteins using only their sequence. Most pipelines for annotating proteins rely on searches for homologous sequences in databases of previously annotated proteins using BLAST or HMMER. Here, we develop a new approach for classifying proteins into a taxonomy of functions and demonstrate its utility for genome annotation. Our algorithm, IDTAXA, was more accurate than BLAST or HMMER at assigning sequences to KEGG ortholog groups. Moreover, IDTAXA correctly avoided classifying sequences with novel functions to existing groups, which is a common error mode for classification approaches that rely on E-values as a proxy for confidence. We demonstrate IDTAXA's utility for annotating eukaryotic and prokaryotic genomes by assigning functions to proteins within a multi-level ontology and applied IDTAXA to detect genome contamination in eukaryotic genomes. Finally, we re-annotated 8604 microbial genomes with known antibiotic resistance phenotypes to discover two novel associations between proteins and antibiotic resistance. IDTAXA is available as a web tool (http://DECIPHER.codes/Classification.html) or as part of the open source DECIPHER R package from Bioconductor.

Insights on the Evolutionary Genomics of the Blautia Genus: Potential New Species and Genetic Content Among Lineages

Blautia, a genus established in 2008, is a relevantly abundant taxonomic group present in the microbiome of human and other mammalian gastrointestinal (GI) tracts. Several described (or proposed) Blautia species are available at this date. However, despite the increasing level of knowledge about Blautia, its diversity is still poorly understood. The increasing availability of Blautia genomic sequences in the public databases opens the possibility to study this genus from a genomic perspective. Here we report the pangenome analysis and the phylogenomic study of 225 Blautia genomes available in RefSeq. We found 33 different potential species at the genomic level, 17 of them previously undescribed; we also confirmed by genomic standards the status of 4 previously proposed new Blautia species. Comparative genomic analyses suggest that the Blautia pangenome is open, with a relatively small core genome (∼ 700-800 gene families). Utilizing a set of representative genomes, we performed a gene family gain/loss model for the genus, showing that despite terminal nodes suffered more massive gene gain events than internal nodes (i.e., predicted ancestors), some ancestors were predicted to have gained an important number of gene families, some of them associated with the possible acquisition of metabolic abilities. Gene loss events remained lower than gain events in most cases. General aspects regarding pangenome composition and gene gain/loss events are discussed, as well as the proposition of changes in the taxonomic assignment of B. coccoides TY and the proposition of a new species, "B. pseudococcoides.".