The IMG/M data management and analysis system v.6.0: new tools and advanced capabilities

Bacteroides ovatus colonization influences the abundance of intestinal short chain fatty acids and neurotransmitters

Gut microbes can synthesize multiple neuro-active metabolites. We profiled neuro-active compounds produced by the gut commensal Bacteroides ovatus in vitro and in vivo by LC-MS/MS. We found that B. ovatus generates acetic acid, propionic acid, isobutyric acid, and isovaleric acid. In vitro, B. ovatus consumed tryptophan and glutamate and synthesized the neuro-active compounds glutamine and GABA. Consistent with our LC-MS/MS-based in vitro data, we observed elevated levels of acetic acid, propionic acid, isobutyric acid, and isovaleric acid in the intestines of B. ovatus mono-associated mice compared with germ-free controls. B. ovatus mono-association also increased the concentrations of intestinal GABA and decreased the concentrations of tryptophan and glutamine compared with germ-free controls. Computational network analysis revealed unique links between SCFAs, neuro-active compounds, and colonization status. These results highlight connections between microbial colonization and intestinal neurotransmitter concentrations, suggesting that B. ovatus selectively influences the presence of intestinal neurotransmitters.

Characterizing the mucin-degrading capacity of the human gut microbiota

Mucin-degrading microbes are known to harbor glycosyl hydrolases (GHs) which cleave specific glycan linkages. Although several microbial species have been identified as mucin degraders, there are likely many other members of the healthy gut community with the capacity to degrade mucins. The aim of the present study was to systematically examine the CAZyme mucin-degrading profiles of the human gut microbiota. Within the Verrucomicrobia phylum, all Akkermansia glycaniphila and muciniphila genomes harbored multiple gene copies of mucin-degrading GHs. The only representative of the Lentisphaerae phylum, Victivallales, harbored a GH profile that closely mirrored Akkermansia. In the Actinobacteria phylum, we found several Actinomadura, Actinomyces, Bifidobacterium, Streptacidiphilus and Streptomyces species with mucin-degrading GHs. Within the Bacteroidetes phylum, Alistipes, Alloprevotella, Bacteroides, Fermenitomonas Parabacteroides, Prevotella and Phocaeicola species had mucin degrading GHs. Firmicutes contained Abiotrophia, Blautia, Enterococcus, Paenibacillus, Ruminococcus, Streptococcus, and Viridibacillus species with mucin-degrading GHs. Interestingly, far fewer mucin-degrading GHs were observed in the Proteobacteria phylum and were found in Klebsiella, Mixta, Serratia and Enterobacter species. We confirmed the mucin-degrading capability of 23 representative gut microbes using a chemically defined media lacking glucose supplemented with porcine intestinal mucus. These data greatly expand our knowledge of microbial-mediated mucin degradation within the human gut microbiota.

Comparative genome analysis among Variovorax species and genome guided aromatic compound degradation analysis emphasizing 4-hydroxybenzoate degradation in Variovorax sp. PAMC26660

Background

While the genus Variovorax is known for its aromatic compound metabolism, no detailed study of the peripheral and central pathways of aromatic compound degradation has yet been reported. Variovorax sp. PAMC26660 is a lichen-associated bacterium isolated from Antarctica. The work presents the genome-based elucidation of peripheral and central catabolic pathways of aromatic compound degradation genes in Variovorax sp. PAMC26660. Additionally, the accessory, core and unique genes were identified among Variovorax species using the pan genome analysis tool. A detailed analysis of the genes related to xenobiotic metabolism revealed the potential roles of Variovorax sp. PAMC26660 and other species in bioremediation.

Results

TYGS analysis, dDDH, phylogenetic placement and average nucleotide identity (ANI) analysis identified the strain as Variovorax sp. Cell morphology was assessed using scanning electron microscopy (SEM). On analysis of the core, accessory, and unique genes, xenobiotic metabolism accounted only for the accessory and unique genes. On detailed analysis of the aromatic compound catabolic genes, peripheral pathway related to 4-hydroxybenzoate (4-HB) degradation was found among all species while phenylacetate and tyrosine degradation pathways were present in most of the species including PAMC26660. Likewise, central catabolic pathways, like protocatechuate, gentisate, homogentisate, and phenylacetyl-CoA, were also present. The peripheral pathway for 4-HB degradation was functionally tested using PAMC26660, which resulted in the growth using it as a sole source of carbon.

Conclusions

Computational tools for genome and pan genome analysis are important to understand the behavior of an organism. Xenobiotic metabolism-related genes, that only account for the accessory and unique genes infer evolution through events like lateral gene transfer, mutation and gene rearrangement. 4-HB, an aromatic compound present among lichen species is utilized by lichen-associated Variovorax sp. PAMC26660 as the sole source of carbon. The strain holds genes and pathways for its utilization. Overall, this study outlines the importance of Variovorax in bioremediation and presents the genomic information of the species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08589-3.

The Notable Achievements and the Prospects of Bacterial Pathogen Genomics

Throughout the entirety of human history, bacterial pathogens have played an important role and even shaped the fate of civilizations. The application of genomics within the last 27 years has radically changed the way we understand the biology and evolution of these pathogens. In this review, we discuss how the short- (Illumina) and long-read (PacBio, Oxford Nanopore) sequencing technologies have shaped the discipline of bacterial pathogen genomics, in terms of fundamental research (i.e., evolution of pathogenicity), forensics, food safety, and routine clinical microbiology. We have mined and discuss some of the most prominent data/bioinformatics resources such as NCBI pathogens, PATRIC, and Pathogenwatch. Based on this mining, we present some of the most popular sequencing technologies, hybrid approaches, assemblers, and annotation pipelines. A small number of bacterial pathogens are of very high importance, and we also present the wealth of the genomic data for these species (i.e., which ones they are, the number of antimicrobial resistance genes per genome, the number of virulence factors). Finally, we discuss how this discipline will probably be transformed in the near future, especially by transitioning into metagenome-assembled genomes (MAGs), thanks to long-read sequencing.

Metagenomic and Metatranscriptomic Insights into Population Diversity of Microcystis Blooms: Spatial and Temporal Dynamics of mcy Genotypes, Including a Partial Operon That Can Be Abundant and Expressed

Cyanobacterial harmful algal blooms (cyanoHABs) degrade freshwater ecosystems globally. Microcystis aeruginosa often dominates cyanoHABs and produces microcystin (MC), a class of hepatotoxins that poses threats to human and animal health. Microcystin toxicity is influenced by distinct structural elements across a diversity of related molecules encoded by variant mcy operons. However, the composition and distribution of mcy operon variants in natural blooms remain poorly understood. Here, we characterized the variant composition of mcy genes in western Lake Erie Microcystis blooms from 2014 and 2018. Sampling was conducted across several spatial and temporal scales, including different bloom phases within 2014, extensive spatial coverage on the same day (2018), and frequent, autonomous sampling over a 2-week period (2018). Mapping of metagenomic and metatranscriptomic sequences to reference sequences revealed three Microcystis mcy genotypes: complete (all genes present [mcyA-J]), partial (truncated mcyA, complete mcyBC, and missing mcyD-J), and absent (no mcy genes). We also detected two different variants of mcyB that may influence the production of microcystin congeners. The relative abundance of these genotypes was correlated with pH and nitrate concentrations. Metatranscriptomic analysis revealed that partial operons were, at times, the most abundant genotype and expressed in situ, suggesting the potential biosynthesis of truncated products. Quantification of genetic divergence between genotypes suggests that the observed strains are the result of preexisting heterogeneity rather than de novo mutation during the sampling period. Overall, our results show that natural Microcystis populations contain several cooccurring mcy genotypes that dynamically shift in abundance spatiotemporally via strain succession and likely influence the observed diversity of the produced congeners. IMPORTANCE Cyanobacteria are responsible for producing microcystins (MCs), a class of potent and structurally diverse toxins, in freshwater systems around the world. While microcystins have been studied for over 50 years, the diversity of their chemical forms and how this variation is encoded at the genetic level remain poorly understood, especially within natural populations of cyanobacterial harmful algal blooms (cyanoHABs). Here, we leverage community DNA and RNA sequences to track shifts in mcy genes responsible for producing microcystin, uncovering the relative abundance, expression, and variation of these genes. We studied this phenomenon in western Lake Erie, which suffers annually from cyanoHAB events, with impacts on drinking water, recreation, tourism, and commercial fishing.

Phenazines and toxoflavin act as interspecies modulators of resilience to diverse antibiotics

Bacterial opportunistic pathogens make diverse secondary metabolites both in the natural environment and when causing infections, yet how these molecules mediate microbial interactions and their consequences for antibiotic treatment are still poorly understood. Here, we explore the role of three redox-active secondary metabolites, pyocyanin, phenazine-1-carboxylic acid and toxoflavin, as interspecies modulators of antibiotic resilience. We find that these molecules dramatically change susceptibility levels of diverse bacteria to clinical antibiotics. Pyocyanin and phenazine-1-carboxylic acid are made by Pseudomonas aeruginosa, while toxoflavin is made by Burkholderia gladioli, organisms that infect cystic fibrosis and other immunocompromised patients. All molecules alter the susceptibility profile of pathogenic species within the "Burkholderia cepacia complex" to different antibiotics, either antagonizing or potentiating their effects, depending on the drug's class. Defense responses regulated by the redox-sensitive transcription factor SoxR potentiate the antagonistic effects these metabolites have against fluoroquinolones, and the presence of genes encoding SoxR and the efflux systems it regulates can be used to predict how these metabolites will affect antibiotic susceptibility of different bacteria. Finally, we demonstrate that inclusion of secondary metabolites in standard protocols used to assess antibiotic resistance can dramatically alter the results, motivating the development of new tests for more accurate clinical assessment.

The Genome of the Acid Soil-Adapted Strain Rhizobium favelukesii OR191 Encodes Determinants for Effective Symbiotic Interaction With Both an Inverted Repeat Lacking Clade and a Phaseoloid Legume Host

Although Medicago sativa forms highly effective symbioses with the comparatively acid-sensitive genus Ensifer, its introduction into acid soils appears to have selected for symbiotic interactions with acid-tolerant R. favelukesii strains. Rhizobium favelukesii has the unusual ability of being able to nodulate and fix nitrogen, albeit sub-optimally, not only with M. sativa but also with the promiscuous host Phaseolus vulgaris. Here we describe the genome of R. favelukesii OR191 and genomic features important for the symbiotic interaction with both of these hosts. The OR191 draft genome contained acid adaptation loci, including the highly acid-inducible lpiA/acvB operon and olsC, required for production of lysine- and ornithine-containing membrane lipids, respectively. The olsC gene was also present in other acid-tolerant Rhizobium strains but absent from the more acid-sensitive Ensifer microsymbionts. The OR191 symbiotic genes were in general more closely related to those found in Medicago microsymbionts. OR191 contained the nodA, nodEF, nodHPQ, and nodL genes for synthesis of polyunsaturated, sulfated and acetylated Nod factors that are important for symbiosis with Medicago, but contained a truncated nodG, which may decrease nodulation efficiency with M. sativa. OR191 contained an E. meliloti type BacA, which has been shown to specifically protect Ensifer microsymbionts from Medicago nodule-specific cysteine-rich peptides. The nitrogen fixation genes nifQWZS were present in OR191 and P. vulgaris microsymbionts but absent from E. meliloti-Medicago microsymbionts. The ability of OR191 to nodulate and fix nitrogen symbiotically with P. vulgaris indicates that this host has less stringent requirements for nodulation than M. sativa but may need rhizobial strains that possess nifQWZS for N₂-fixation to occur. OR191 possessed the exo genes required for the biosynthesis of succinoglycan, which is required for the Ensifer-Medicago symbiosis. However, ¹H-NMR spectra revealed that, in the conditions tested, OR191 exopolysaccharide did not contain a succinyl substituent but instead contained a 3-hydroxybutyrate moiety, which may affect its symbiotic performance with Medicago hosts. These findings provide a foundation for the genetic basis of nodulation requirements and symbiotic effectiveness with different hosts.

Identifying candidate structured RNAs in CRISPR operons

Noncoding RNAs with secondary structures play important roles in CRISPR-Cas systems. Many of these structures likely remain undiscovered. We used a large-scale comparative genomics approach to predict 156 novel candidate structured RNAs from 36,111 CRISPR-Cas systems. A number of these were found to overlap with coding genes, including palindromic candidates that overlapped with a variety of Cas genes in type I and III systems. Among these 156 candidates, we identified 46 new models of CRISPR direct repeats and 1 tracrRNA. This tracrRNA model occasionally overlapped with predicted cas9 coding regions, emphasizing the importance of expanding our search windows for novel structure RNAs in coding regions. We also demonstrated that the antirepeat sequence in this tracrRNA model can be used to accurately assign thousands of predicted CRISPR arrays to type II-C systems. This study highlights the importance of unbiased identification of candidate structured RNAs across CRISPR-Cas systems.

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.

Comparative Genomics of Xylella fastidiosa Explores Candidate Host-Specificity Determinants and Expands the Known Repertoire of Mobile Genetic Elements and Immunity Systems

Xylella fastidiosa causes diseases in many plant species. Originally confined to the Americas, infecting mainly grapevine, citrus, and coffee, X. fastidiosa has spread to several plant species in Europe causing devastating diseases. Many pathogenicity and virulence factors have been identified, which enable the various X. fastidiosa strains to successfully colonize the xylem tissue and cause disease in specific plant hosts, but the mechanisms by which this happens have not been fully elucidated. Here we present thorough comparative analyses of 94 whole-genome sequences of X. fastidiosa strains from diverse plant hosts and geographic regions. Core-genome phylogeny revealed clades with members sharing mostly a geographic region rather than a host plant of origin. Phylogenetic trees for 1605 orthologous CDSs were explored for potential candidates related to host specificity using a score of mapping metrics. However, no candidate host-specificity determinants were strongly supported using this approach. We also show that X. fastidiosa accessory genome is represented by an abundant and heterogeneous mobilome, including a diversity of prophage regions. Our findings provide a better understanding of the diversity of phylogenetically close genomes and expand the knowledge of X. fastidiosa mobile genetic elements and immunity systems.

A functional Wood-Ljungdahl pathway devoid of a formate dehydrogenase in the gut acetogens Blautia wexlerae, Blautia luti and beyond

Species of the genus Blautia are typical inhabitants of the human gut and considered as beneficial gut microbes. However, their role in the gut microbiome and their metabolic features are poorly understood. Blautia schinkii was described as an acetogenic bacterium, characterized by a functional Wood-Ljungdahl pathway (WLP) of acetogenesis from H₂ + CO₂ . Here we report that two relatives, Blautia luti and Blautia wexlerae do not grow on H₂ + CO₂ . Inspection of the genome sequence revealed all genes of the WLP except genes encoding a formate dehydrogenase and an electron-bifurcating hydrogenase. Enzyme assays confirmed this prediction. Accordingly, resting cells neither converted H₂ + CO₂ nor H₂ + HCOOH + CO₂ to acetate. Carbon monoxide is an intermediate of the WLP and substrate for many acetogens. B. luti and B. wexlerae had an active CO dehydrogenase and resting cells performed acetogenesis from HCOOH + CO₂ + CO, demonstrating a functional WLP. Bioinformatic analyses revealed that many Blautia strains as well as other gut acetogens lack formate dehydrogenases and hydrogenases. Thus, the use of formate instead of H₂ + CO₂ as an interspecies hydrogen and electron carrier seems to be more common in the gut microbiome. This article is protected by copyright. All rights reserved.

Leveraging big data bioinformatics approaches to extract knowledge from Staphylococcus aureus public omics data

Staphylococcus aureus is a notorious pathogen posing challenges in the medical industry due to drug resistance and biofilm formation. The horizon of knowledge on S. aureus pathogenesis has expanded with the advancement of data-driven bioinformatics techniques. Mining information from sequenced genomes and their expression data is an economic approach that alleviates wastage of resources and redundancy in experiments. The current review covers how big data bioinformatics has been used in the analysis of S. aureus from publicly available -omics data to uncover mechanisms of infection and inhibition. Particularly, advances in the past two decades in biomarker discovery, host responses, phenotype identification, consolidation of information, and drug development are discussed highlighting the challenges and shortcomings. Overall, the review summarizes the diverse aspects of scrupulous re-analysis of S. aureus proteomic and transcriptomic expression datasets retrieved from public repositories in terms of the efforts taken, benefits offered, and follow-up actions. The detailed review thus serves as a reference and aid for (i) Computational biologists by briefing the approaches utilized for bacterial omics re-analysis concerning S. aureus and (ii) Experimental biologists by elucidating the potential of bioinformatics in biological research to generate reliable postulates in a prompt and economical manner.

Novel Viral DNA Polymerases From Metagenomes Suggest Genomic Sources of Strand-Displacing Biochemical Phenotypes

Viruses are the most abundant and diverse biological entities on the planet and constitute a significant proportion of Earth's genetic diversity. Most of this diversity is not represented by isolated viral-host systems and has only been observed through sequencing of viral metagenomes (viromes) from environmental samples. Viromes provide snapshots of viral genetic potential, and a wealth of information on viral community ecology. These data also provide opportunities for exploring the biochemistry of novel viral enzymes. The in vitro biochemical characteristics of novel viral DNA polymerases were explored, testing hypothesized differences in polymerase biochemistry according to protein sequence phylogeny. Forty-eight viral DNA Polymerase I (PolA) proteins from estuarine viromes, hot spring metagenomes, and reference viruses, encompassing a broad representation of currently known diversity, were synthesized, expressed, and purified. Novel functionality was shown in multiple PolAs. Intriguingly, some of the estuarine viral polymerases demonstrated moderate to strong innate DNA strand displacement activity at high enzyme concentration. Strand-displacing polymerases have important technological applications where isothermal reactions are desirable. Bioinformatic investigation of genes neighboring these strand displacing polymerases found associations with SNF2 helicase-associated proteins. The specific function of SNF2 family enzymes is unknown for prokaryotes and viruses. In eukaryotes, SNF2 enzymes have chromatin remodeling functions but do not separate nucleic acid strands. This suggests the strand separation function may be fulfilled by the DNA polymerase for viruses carrying SNF2 helicase-associated proteins. Biochemical data elucidated from this study expands understanding of the biology and ecological behavior of unknown viruses. Moreover, given the numerous biotechnological applications of viral DNA polymerases, novel viral polymerases discovered within viromes may be a rich source of biological material for further in vitro DNA amplification advancements.

An Unprecedented Number of Cytochrome P450s Are Involved in Secondary Metabolism in Salinispora Species

Cytochrome P450 monooxygenases (CYPs/P450s) are heme thiolate proteins present in species across the biological kingdoms. By virtue of their broad substrate promiscuity and regio- and stereo-selectivity, these enzymes enhance or attribute diversity to secondary metabolites. Actinomycetes species are well-known producers of secondary metabolites, especially Salinispora species. Despite the importance of P450s, a comprehensive comparative analysis of P450s and their role in secondary metabolism in Salinispora species is not reported. We therefore analyzed P450s in 126 strains from three different species Salinispora arenicola, S. pacifica, and S. tropica. The study revealed the presence of 2643 P450s that can be grouped into 45 families and 103 subfamilies. CYP107 and CYP125 families are conserved, and CYP105 and CYP107 families are bloomed (a P450 family with many members) across Salinispora species. Analysis of P450s that are part of secondary metabolite biosynthetic gene clusters (smBGCs) revealed Salinispora species have an unprecedented number of P450s (1236 P450s-47%) part of smBGCs compared to other bacterial species belonging to the genera Streptomyces (23%) and Mycobacterium (11%), phyla Cyanobacteria (8%) and Firmicutes (18%) and the classes Alphaproteobacteria (2%) and Gammaproteobacteria (18%). A peculiar characteristic of up to six P450s in smBGCs was observed in Salinispora species. Future characterization Salinispora species P450s and their smBGCs have the potential for discovering novel secondary metabolites.

Phylogenetic Distribution and Evolution of Type VI Secretion System in the Genus Xanthomonas

The type VI secretion system (T6SS) present in many Gram-negative bacteria is a contact-dependent apparatus that can directly deliver secreted effectors or toxins into diverse neighboring cellular targets including both prokaryotic and eukaryotic organisms. Recent reverse genetics studies with T6 core gene loci have indicated the importance of functional T6SS toward overall competitive fitness in various pathogenic Xanthomonas spp. To understand the contribution of T6SS toward ecology and evolution of Xanthomonas spp., we explored the distribution of the three distinguishable T6SS clusters, i3*, i3***, and i4, in approximately 1,740 Xanthomonas genomes, along with their conservation, genetic organization, and their evolutionary patterns in this genus. Screening genomes for core genes of each T6 cluster indicated that 40% of the sequenced strains possess two T6 clusters, with combinations of i3*** and i3* or i3*** and i4. A few strains of Xanthomonas citri, Xanthomonas phaseoli, and Xanthomonas cissicola were the exception, possessing a unique combination of i3* and i4. The findings also indicated clade-specific distribution of T6SS clusters. Phylogenetic analysis demonstrated that T6SS clusters i3* and i3*** were probably acquired by the ancestor of the genus Xanthomonas, followed by gain or loss of individual clusters upon diversification into subsequent clades. T6 i4 cluster has been acquired in recent independent events by group 2 xanthomonads followed by its spread via horizontal dissemination across distinct clades across groups 1 and 2 xanthomonads. We also noted reshuffling of the entire core T6 loci, as well as T6SS spike complex components, hcp and vgrG, among different species. Our findings indicate that gain or loss events of specific T6SS clusters across Xanthomonas phylogeny have not been random.

New insights into the energy metabolism and taxonomy of Deferribacteres revealed by the characterization of a new isolate from a hypersaline microbial mat

Strain L21-Ace-BES^T , isolated from a lithifying cyanobacterial mat, could be assigned to a novel species and genus within the Deferribacteres. It is an important model organism for the study of anaerobic acetate degradation under hypersaline conditions. The metabolism of strain L21-Ace-BES^T was characterized by biochemical studies, comparative genome analyses, and the evaluation of gene expression patterns. The central metabolic pathway is the citric acid cycle, which is mainly controlled by the enzyme succinyl-CoA:acetate-CoA transferase. The potential use of a reversed oxidative citric acid cycle to fix CO₂ has been revealed through genome analysis. However, no autotrophic growth was detected in this strain, whereas sulfide and H₂ can be used mixotrophically. Preferred electron acceptors for the anaerobic oxidation of acetate are nitrate, fumarate and DMSO, while oxygen can be utilized only under microoxic conditions. Aerotolerant growth by fermentation was observed at higher oxygen concentrations. The redox cycling of sulfur/sulfide enables the generation of reducing power for the assimilation of acetate during growth and could prevent the over-reduction of cells in stationary phase. Extracellular electron transfer appears to be an essential component of the respiratory metabolism in this clade of Deferribacteres and may be involved in the reduction of nitrite to ammonium. This article is protected by copyright. All rights reserved.

Scoring protein sequence alignments using deep Learning

BACKGROUND

A high-quality sequence alignment (SA) is the most important input feature for accurate protein structure prediction. For a protein sequence, there are many methods to generate a SA. However, when given a choice of more than one SA for a protein sequence, there are no methods to predict which SA may lead to more accurate models without actually building the models. In this work, we describe a method to predict the quality of a protein's SA.

METHODS

We created our own dataset by generating a variety of SAs for a set of 1,351 representative proteins and investigated various deep learning architectures to predict the local distance difference test (lDDT) scores of distance maps predicted with SAs as the input. These lDDT scores serve as indicators of the quality of the SAs.

RESULTS

Using two independent test datasets consisting of CASP13 and CASP14 targets, we show that our method is effective for scoring and ranking SAs when a pool of SAs is available for a protein sequence. With an example, we further discuss that SA selection using our method can lead to improved structure prediction.

AVAILABILITY

Code and datasets are available at https://github.com/ba-lab/Alignment-Score/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mechanistic insights into consumption of the food additive xanthan gum by the human gut microbiota

Processed foods often include food additives such as xanthan gum, a complex polysaccharide with unique rheological properties, that has established widespread use as a stabilizer and thickening agent. Xanthan gum's chemical structure is distinct from those of host and dietary polysaccharides that are more commonly expected to transit the gastrointestinal tract, and little is known about its direct interaction with the gut microbiota, which plays a central role in digestion of other dietary fibre polysaccharides. Here we show that the ability to digest xanthan gum is common in human gut microbiomes from industrialized countries and appears contingent on a single uncultured bacterium in the family Ruminococcaceae. Our data reveal that this primary degrader cleaves the xanthan gum backbone before processing the released oligosaccharides using additional enzymes. Some individuals harbour Bacteroides intestinalis that is incapable of consuming polymeric xanthan gum but grows on oligosaccharide products generated by the Ruminococcaceae. Feeding xanthan gum to germfree mice colonized with a human microbiota containing the uncultured Ruminococcaceae supports the idea that the additive xanthan gum can drive expansion of the primary degrader Ruminococcaceae, along with exogenously introduced B. intestinalis. Our work demonstrates the existence of a potential xanthan gum food chain involving at least two members of different phyla of gut bacteria and provides an initial framework for understanding how widespread consumption of a recently introduced food additive influences human microbiomes.

Metagenomics-enabled microbial surveillance

Lessons learnt from the COVID-19 pandemic include increased awareness of the potential for zoonoses and emerging infectious diseases that can adversely affect human health. Although emergent viruses are currently in the spotlight, we must not forget the ongoing toll of morbidity and mortality owing to antimicrobial resistance in bacterial pathogens and to vector-borne, foodborne and waterborne diseases. Population growth, planetary change, international travel and medical tourism all contribute to the increasing frequency of infectious disease outbreaks. Surveillance is therefore of crucial importance, but the diversity of microbial pathogens, coupled with resource-intensive methods, compromises our ability to scale-up such efforts. Innovative technologies that are both easy to use and able to simultaneously identify diverse microorganisms (viral, bacterial or fungal) with precision are necessary to enable informed public health decisions. Metagenomics-enabled surveillance methods offer the opportunity to improve detection of both known and yet-to-emerge pathogens.

The role of nitrogen-responsive regulators in controlling inorganic polyphosphate synthesis in Escherichia coli

Inorganic polyphosphate (polyP) is synthesized by bacteria under stressful environmental conditions and acts by a variety of mechanisms to promote cell survival. While the kinase that synthesizes polyP (PPK, encoded by the ppk gene) is well known, ppk transcription is not activated by environmental stress and little is understood about how environmental stress signals lead to polyP accumulation. Previous work has shown that the transcriptional regulators DksA, RpoN (σ54) and RpoE (σ24) positively regulate polyP production, but not ppk transcription, in Escherichia coli. In this work, we examine the role of the alternative sigma factor RpoN and nitrogen starvation stress response pathways in controlling polyP synthesis. We show that the RpoN enhancer binding proteins GlnG and GlrR impact polyP production, and uncover a new role for the nitrogen phosphotransferase regulator PtsN (EIIANtr) as a positive regulator of polyP production, acting upstream of DksA, downstream of RpoN and apparently independently of RpoE. However, neither these regulatory proteins nor common nitrogen metabolites appear to act directly on PPK, and the precise mechanism(s) by which polyP production is modulated after stress remain(s) unclear. Unexpectedly, we also found that the genes that impact polyP production vary depending on the composition of the rich media in which the cells were grown before exposure to polyP-inducing stress. These results constitute progress towards deciphering the regulatory networks driving polyP production under stress, and highlight the remarkable complexity of this regulation and its connections to a broad range of stress-sensing pathways.

Diverse ecophysiological adaptations of subsurface Thaumarchaeota in floodplain sediments revealed through genome-resolved metagenomics

The terrestrial subsurface microbiome contains vastly underexplored phylogenetic diversity and metabolic novelty, with critical implications for global biogeochemical cycling. Among the key microbial inhabitants of subsurface soils and sediments are Thaumarchaeota, an archaeal phylum that encompasses ammonia-oxidizing archaea (AOA) as well as non-ammonia-oxidizing basal lineages. Thaumarchaeal ecology in terrestrial systems has been extensively characterized, particularly in the case of AOA. However, there is little knowledge on the diversity and ecophysiology of Thaumarchaeota in deeper soils, as most lineages, particularly basal groups, remain uncultivated and underexplored. Here we use genome-resolved metagenomics to examine the phylogenetic and metabolic diversity of Thaumarchaeota along a 234 cm depth profile of hydrologically variable riparian floodplain sediments in the Wind River Basin near Riverton, Wyoming. Phylogenomic analysis of the metagenome-assembled genomes (MAGs) indicates a shift in AOA population structure from the dominance of the terrestrial Nitrososphaerales lineage in the well-drained top ~100 cm of the profile to the typically marine Nitrosopumilales in deeper, moister, more energy-limited sediment layers. We also describe two deeply rooting non-AOA MAGs with numerous unexpected metabolic features, including the reductive acetyl-CoA (Wood-Ljungdahl) pathway, tetrathionate respiration, a form III RuBisCO, and the potential for extracellular electron transfer. These MAGs also harbor tungsten-containing aldehyde:ferredoxin oxidoreductase, group 4f [NiFe]-hydrogenases and a canonical heme catalase, typically not found in Thaumarchaeota. Our results suggest that hydrological variables, particularly proximity to the water table, impart a strong control on the ecophysiology of Thaumarchaeota in alluvial sediments.

Analysis of laccase‐like enzymes secreted by fungi isolated from a cave in northern Spain

Laccases belong to a family of multicopper enzymes able to oxidize a broad spectrum of organic compounds. Despite the well‐known property of laccases to carry out bleaching and degradation of industrial dyes and polyphenolic compounds, their industrial use is often limited by the high cost, low efficiency, or instability of these enzymes. To look for new microorganisms which produce laccases that are potentially suitable for industrial applications, we have isolated several fungal strains from a cave in northern Spain. Their phenotypic analysis on agar plates supplemented with ABTS (2,2′‐azino‐bis(3‐ethylbenzothiazoline‐6‐sulfonic acid)) disclosed two laccase‐positive strains. Further genotyping revealed that they belonged to the Gliomastix murorum and Conidiobolus thromboides species. The secretion of G. murorum and C. thromboides laccase‐like enzymes was then confirmed by zymography. Further identification of these polypeptides by mass‐spectroscopy revealed the nature of the laccases and made it possible to predict their functional domains and other features. In addition, plate assays revealed that the laccases secreted by both G. murorum and C. thromboides were capable of degrading industrial dyes (Congo Red, Indigo, and Eriochrome Black T). Homology modeling and substrate docking predicted the putative structure of the currently uncrystallized G. murorum enzyme as well as its amino acid residues potentially involved in interactions with these dyes. In summary, new biochemical and structural insights into decolorization mediated by G. murorum laccase as well as identification of laccase‐like oxidase in C. thromboides point to a promising future for these enzymes in biotechnology.

Metagenomes and Metagenome-Assembled Genomes from Substrate-Amended Hot Spring Sediment Incubations from Yellowstone National Park

Here, we report on eight sediment metagenomes obtained from an alkaline hot spring, with their corresponding metagenome-assembled genomes. Samples had been incubated for 48 h with various substrate amendments in conjunction with the amino acid analog l-homopropargylglycine in a study targeted at identifying anabolicly active uncultured thermophilic archaea and bacteria.

deepNEC: a novel alignment-free tool for the identification and classification of nitrogen biochemical network-related enzymes using deep learning

Nitrogen is essential for life and its transformations are an important part of the global biogeochemical cycle. Being an essential nutrient, nitrogen exists in a range of oxidation states from +5 (nitrate) to -3 (ammonium and amino-nitrogen), and its oxidation and reduction reactions catalyzed by microbial enzymes determine its environmental fate. The functional annotation of the genes encoding the core nitrogen network enzymes has a broad range of applications in metagenomics, agriculture, wastewater treatment and industrial biotechnology. This study developed an alignment-free computational approach to determine the predicted nitrogen biochemical network-related enzymes from the sequence itself. We propose deepNEC, a novel end-to-end feature selection and classification model training approach for nitrogen biochemical network-related enzyme prediction. The algorithm was developed using Deep Learning, a class of machine learning algorithms that uses multiple layers to extract higher-level features from the raw input data. The derived protein sequence is used as an input, extracting sequential and convolutional features from raw encoded protein sequences based on classification rather than traditional alignment-based methods for enzyme prediction. Two large datasets of protein sequences, enzymes and non-enzymes were used to train the models with protein sequence features like amino acid composition, dipeptide composition (DPC), conformation transition and distribution, normalized Moreau-Broto (NMBroto), conjoint and quasi order, etc. The k-fold cross-validation and independent testing were performed to validate our model training. deepNEC uses a four-tier approach for prediction; in the first phase, it will predict a query sequence as enzyme or non-enzyme; in the second phase, it will further predict and classify enzymes into nitrogen biochemical network-related enzymes or non-nitrogen metabolism enzymes; in the third phase, it classifies predicted enzymes into nine nitrogen metabolism classes; and in the fourth phase, it predicts the enzyme commission number out of 20 classes for nitrogen metabolism. Among all, the DPC + NMBroto hybrid feature gave the best prediction performance (accuracy of 96.15% in k-fold training and 93.43% in independent testing) with an Matthews correlation coefficient (0.92 training and 0.87 independent testing) in phase I; phase II (accuracy of 99.71% in k-fold training and 98.30% in independent testing); phase III (overall accuracy of 99.03% in k-fold training and 98.98% in independent testing); phase IV (overall accuracy of 99.05% in k-fold training and 98.18% in independent testing), the DPC feature gave the best prediction performance. We have also implemented a homology-based method to remove false negatives. All the models have been implemented on a web server (prediction tool), which is freely available at http://bioinfo.usu.edu/deepNEC/.

Assessment of global health risk of antibiotic resistance genes

Antibiotic resistance genes (ARGs) have accelerated microbial threats to human health in the last decade. Many genes can confer resistance, but evaluating the relative health risks of ARGs is complex. Factors such as the abundance, propensity for lateral transmission and ability of ARGs to be expressed in pathogens are all important. Here, an analysis at the metagenomic level from various habitats (6 types of habitats, 4572 samples) detects 2561 ARGs that collectively conferred resistance to 24 classes of antibiotics. We quantitatively evaluate the health risk to humans, defined as the risk that ARGs will confound the clinical treatment for pathogens, of these 2561 ARGs by integrating human accessibility, mobility, pathogenicity and clinical availability. Our results demonstrate that 23.78% of the ARGs pose a health risk, especially those which confer multidrug resistance. We also calculate the antibiotic resistance risks of all samples in four main habitats, and with machine learning, successfully map the antibiotic resistance threats in global marine habitats with over 75% accuracy. Our novel method for quantitatively surveilling the health risk of ARGs will help to manage one of the most important threats to human and animal health.

Antibiotic resistance genes (ARGs) have accelerated microbial threats to human health. Here, Zhang et al. analyze 4572 metagenomic samples to illustrate the global patterns of ARG distribution in diverse habitats. They quantitatively evaluate the health risk to humans of 2561 ARGs by integrating human accessibility, mobility, pathogenicity and clinical availability. With the machine learning, they map the antibiotic resistance threats in global marine habitats.

Increased Butyrate Production in Clostridium saccharoperbutylacetonicum from Lignocellulose-Derived Sugars

Butyrate is produced by chemical synthesis based on crude oil, produced by microbial fermentation, or extracted from animal fats (M. Dwidar, J.-Y. Park, R. J. Mitchell, and B.-I. Sang, The Scientific World Journal, 2012:471417, 2012, https://doi.org/10.1100/2012/471417). Butyrate production by anaerobic bacteria is highly favorable since waste or sustainable resources can be used as the substrates. For this purpose, the native hyper-butanol producer Clostridium saccharoperbutylacetonicum N1-4(HMT) was used as a chassis strain due to its broad substrate spectrum. BLASTp analysis of the predicted proteome of C. saccharoperbutylacetonicum N1-4(HMT) resulted in the identification of gene products potentially involved in acetone-butanol-ethanol (ABE) fermentation. Their participation in ABE fermentation was either confirmed or disproven by the parallel production of acids or solvents and the respective transcript levels obtained by transcriptome analysis of this strain. The genes encoding phosphotransacetylase (pta) and butyraldehyde dehydrogenase (bld) were deleted to reduce acetate and alcohol formation. The genes located in the butyryl-CoA synthesis (bcs) operon encoding crotonase, butyryl-CoA dehydrogenase with electron-transferring protein subunits α and β, and 3-hydroxybutyryl-CoA dehydrogenase were overexpressed to channel the flux further towards butyrate formation. Thereby, the native hyper-butanol producer C. saccharoperbutylacetonicum N1-4(HMT) was converted into the hyper-butyrate producer C. saccharoperbutylacetonicum ΔbldΔpta [pMTL83151_BCS_P_bgaL]. The transcription pattern following deletion and overexpression was characterized by a second transcriptomic study, revealing partial compensation for the deletion. Furthermore, this strain was characterized in pH-controlled fermentations with either glucose or Excello, a substrate yielded from spruce biomass. Butyrate was the main product, with maximum butyrate concentrations of 11.7 g·L^-1 and 14.3 g·L^-1, respectively. Minimal amounts of by-products were detected. IMPORTANCE Platform chemicals such as butyrate are usually produced chemically from crude oil, resulting in the carry-over of harmful compounds. The selective production of butyrate using sustainable resources or waste without harmful by-products can be achieved by bacteria such as clostridia. The hyper-butanol producer Clostridium saccharoperbutylacetonicum N1-4(HMT) was converted into a hyper-butyrate producer. Butyrate production with very small amounts of by-products was established with glucose and the sustainable lignocellulosic sugar substrate Excello extracted from spruce biomass by the biorefinery Borregaard (Sarpsborg, Norway).

Biochemical Basis of Xylooligosaccharide Utilisation by Gut Bacteria

Xylan is one of the major structural components of the plant cell wall. Xylan present in the human diet reaches the large intestine undigested and becomes a substrate to species of the gut microbiota. Here, we characterised the capacity of Limosilactobacillus reuteri and Blautia producta strains to utilise xylan derivatives. We showed that L. reuteri ATCC 53608 and B. producta ATCC 27340 produced β-D-xylosidases, enabling growth on xylooligosaccharide (XOS). The recombinant enzymes were highly active on artificial (p-nitrophenyl β-D-xylopyranoside) and natural (xylobiose, xylotriose, and xylotetraose) substrates, and showed transxylosylation activity and tolerance to xylose inhibition. The enzymes belong to glycoside hydrolase family 120 with Asp as nucleophile and Glu as proton donor, as shown by homology modelling and confirmed by site-directed mutagenesis. In silico analysis revealed that these enzymes were part of a gene cluster in L. reuteri but not in Blautia strains, and quantitative proteomics identified other enzymes and transporters involved in B. producta XOS utilisation. Based on these findings, we proposed a model for an XOS metabolism pathway in L. reuteri and B. producta strains. Together with phylogenetic analyses, the data also revealed the extended xylanolytic potential of the gut microbiota.

MetaDecoder: a novel method for clustering metagenomic contigs

BACKGROUND

Clustering the metagenomic contigs into potential genomes is a key step to investigate the functional roles of microbial populations. Existing algorithms have achieved considerable success with simulated or real sequencing datasets. However, accurately classifying contigs from complex metagenomes is still a challenge.

RESULTS

We introduced a novel clustering algorithm, MetaDecoder, which can classify metagenomic contigs based on the frequencies of k-mers and coverages. MetaDecoder was built as a two-layer model with the first layer being a GPU-based modified Dirichlet process Gaussian mixture model (DPGMM), which controls the weight of each DPGMM cluster to avoid over-segmentation by dynamically dissolving contigs in small clusters and reassigning them to the remaining clusters. The second layer comprises a semi-supervised k-mer frequency probabilistic model and a modified Gaussian mixture model for modeling the coverage based on single copy marker genes. Benchmarks on simulated and real-world datasets demonstrated that MetaDecoder can be served as a promising approach for effectively clustering metagenomic contigs.

CONCLUSIONS

In conclusion, we developed the GPU-based MetaDecoder for effectively clustering metagenomic contigs and reconstructing microbial communities from microbial data. Applying MetaDecoder on both simulated and real-world datasets demonstrated that it could generate more complete clusters with lower contamination. Using MetaDecoder, we identified novel high-quality genomes and expanded the existing catalog of bacterial genomes. Video Abstract.

Diverse events have transferred genes for edible seaweed digestion from marine to human gut bacteria

Humans harbor numerous species of colonic bacteria that digest fiber polysaccharides in commonly consumed terrestrial plants. More recently in history, regional populations have consumed edible macroalgae seaweeds containing unique polysaccharides. It remains unclear how extensively gut bacteria have adapted to digest these nutrients. Here, we show that the ability of gut bacteria to digest seaweed polysaccharides is more pervasive than previously appreciated. Enrichment-cultured Bacteroides harbor previously discovered genes for seaweed degradation, which have mobilized into several members of this genus. Additionally, other examples of marine bacteria-derived genes, and their mobile DNA elements, are involved in gut microbial degradation of seaweed polysaccharides, including genes in gut-resident Firmicutes. Collectively, these results uncover multiple separate events that have mobilized the genes encoding seaweed-degrading-enzymes into gut bacteria. This work further underscores the metabolic plasticity of the human gut microbiome and global exchange of genes in the context of dietary selective pressures.

Genome-Based Classification of Strain 16-SW-7, a Marine Bacterium Capable of Converting B Red Blood Cells, as Pseudoalteromonas distincta and Proposal to Reclassify Pseudoalteromonas paragorgicola as a Later Heterotypic Synonym of Pseudoalteromonas distincta

A strictly aerobic, Gram-stain-negative, rod-shaped, and motile bacterium, designated strain 16-SW-7, isolated from a seawater sample, was investigated in detail due to its ability to produce a unique α-galactosidase converting B red blood cells into the universal type blood cells. The phylogenetic analysis based on 16S rRNA gene sequences revealed that the strain 16-SW-7 is a member of the Gammaproteobacteria genus Pseudoalteromonas. The closest relatives of the environmental isolate were Pseudoalteromonas distincta KMM 638^T and Pseudoalteromonas paragorgicola KMM 3548^T, with the plural paralogous 16S rRNA genes of 99.87-100% similarity. The strain 16-SW-7 grew with 1-10% NaCl and at 4-34°C, and hydrolyzed casein, gelatin, tyrosine, and DNA. The genomic DNA G+C content was 39.3 mol%. The prevalent fatty acids were C_16:1 ω7c, C_16:0, C_17:1 ω8c, C_18:1 ω7c, C_17:0, and C_12:0 3-OH. The polar lipid profile was characterized by the presence of phosphatidylethanolamine, phosphatidylglycerol, two unidentified amino lipids, and three unidentified lipids. The major respiratory quinone was Q-8. The finished genome of the strain 16-SW-7 (GenBank assembly accession number: GCA_005877035.1) has a size of 4,531,445 bp and comprises two circular chromosomes L1 and S1, deposited in the GenBank under the accession numbers CP040558 and CP040559, respectively. The strain 16-SW-7 has the ANI values of 98.2% with KMM 638^T and KMM 3548^T and the DDH values of 84.4 and 83.5%, respectively, indicating clearly that the three strains belonged to a single species. According to phylogenetic evidence and similarity for the chemotaxonomic and genotypic properties, the strain 16-SW-7 (= KCTC 52772 = KMM 701) represents a novel member of the species Pseudoalteromonas distincta. Also, we have proposed to reclassify Pseudoalteromonas paragorgicola as a later heterotypic synonym of P. distincta based on the rules of priority with the emendation of the species.

Scaling laws in enzyme function reveal a new kind of biochemical universality

Known examples of life all share the same core biochemistry going back to the last universal common ancestor (LUCA), but whether this feature is universal to other examples, including at the origin of life or alien life, is unknown. We show how a physics-inspired statistical approach identifies universal scaling laws across biochemical reactions that are not defined by common chemical components but instead, as macroscale patterns in the reaction functions used by life. The identified scaling relations can be used to predict statistical features of LUCA, and network analyses reveal some of the functional principles that underlie them. They are, therefore, prime candidates for developing new theory on the “laws of life” that might apply to all possible biochemistries.

All life on Earth is unified by its use of a shared set of component chemical compounds and reactions, providing a detailed model for universal biochemistry. However, this notion of universality is specific to known biochemistry and does not allow quantitative predictions about examples not yet observed. Here, we introduce a more generalizable concept of biochemical universality that is more akin to the kind of universality found in physics. Using annotated genomic datasets including an ensemble of 11,955 metagenomes, 1,282 archaea, 11,759 bacteria, and 200 eukaryotic taxa, we show how enzyme functions form universality classes with common scaling behavior in their relative abundances across the datasets. We verify that these scaling laws are not explained by the presence of compounds, reactions, and enzyme functions shared across known examples of life. We demonstrate how these scaling laws can be used as a tool for inferring properties of ancient life by comparing their predictions with a consensus model for the last universal common ancestor (LUCA). We also illustrate how network analyses shed light on the functional principles underlying the observed scaling behaviors. Together, our results establish the existence of a new kind of biochemical universality, independent of the details of life on Earth’s component chemistry, with implications for guiding our search for missing biochemical diversity on Earth or for biochemistries that might deviate from the exact chemical makeup of life as we know it, such as at the origins of life, in alien environments, or in the design of synthetic life.

Structural basis for an unprecedented enzymatic alkylation in cylindrocyclophane biosynthesis

The cyanobacterial enzyme CylK assembles the cylindrocyclophane natural products by performing two unusual alkylation reactions, forming new carbon–carbon bonds between aromatic rings and secondary alkyl halide substrates. This transformation is unprecedented in biology, and the structure and mechanism of CylK are unknown. Here, we report X-ray crystal structures of CylK, revealing a distinctive fusion of a Ca²⁺-binding domain and a β-propeller fold. We use a mutagenic screening approach to locate CylK’s active site at its domain interface, identifying two residues, Arg105 and Tyr473, that are required for catalysis. Anomalous diffraction datasets collected with bound bromide ions, a product analog, suggest that these residues interact with the alkyl halide electrophile. Additional mutagenesis and molecular dynamics simulations implicate Asp440 in activating the nucleophilic aromatic ring. Bioinformatic analysis of CylK homologs from other cyanobacteria establishes that they conserve these key catalytic amino acids, but they are likely associated with divergent reactivity and altered secondary metabolism. By gaining a molecular understanding of this unusual biosynthetic transformation, this work fills a gap in our understanding of how alkyl halides are activated and used by enzymes as biosynthetic intermediates, informing enzyme engineering, catalyst design, and natural product discovery.

Investigation and Alteration of Organic Acid Synthesis Pathways in the Mammalian Gut Symbiont Bacteroides thetaiotaomicron

Members of the gut-dwelling Bacteroides genus have remarkable abilities in degrading a diverse set of fiber polysaccharide structures, most of which are found in the mammalian diet. As part of their metabolism, they convert these fibers to organic acids that can in turn provide energy to their host. While many studies have identified and characterized the genes and corresponding proteins involved in polysaccharide degradation, relatively little is known about Bacteroides genes involved in downstream metabolic pathways. Bacteroides thetaiotaomicron is one of the most studied species from the genus and is representative of this group in producing multiple organic acids as part of its metabolism. We focused here on several organic acid synthesis pathways in B. thetaiotaomicron, including those involved in formate, lactate, propionate, and acetate production. We identified potential genes involved in each pathway and characterized these through gene deletions coupled to growth assays and organic acid quantification. In addition, we developed and employed a Golden Gate-compatible plasmid system to simplify alteration of native gene expression levels. Our work both validates and contradicts previous bioinformatic gene annotations, and we develop a model on which to base future efforts. A clearer understanding of Bacteroides metabolic pathways can inform and facilitate efforts to employ these bacteria for improved human health or other utilization strategies.

IMPORTANCE Both humans and animals host a large community of bacteria and other microorganisms in their gastrointestinal tracts. This community breaks down dietary fiber and produces organic acids that are used as an energy source by the body and can also help the host resist infection by various pathogens. While the Bacteroides genus is one of the most common in the gut microbiota, it is only distantly related to bacteria with well-characterized metabolic pathways and it is therefore unclear whether research insights on organic acid production in those species can also be directly applied to the Bacteroides. By investigating multiple genetic pathways for organic acid production in Bacteroides thetaiotaomicron, we provide a basis for deeper understanding of these pathways. The work further enables greater understanding of Bacteroides–host relationships, as well as inter-species relationships in the microbiota, which are of importance for both human and animal gut health.

Protease activities of vaginal Porphyromonas species disrupt coagulation and extracellular matrix in the cervicovaginal niche

Porphyromonas asaccharolytica and Porphyromonas uenonis are common inhabitants of the vaginal microbiome, but their presence has been linked to adverse health outcomes for women, including bacterial vaginosis and preterm birth. However, little is known about the pathogenesis mechanisms of these bacteria. The related oral opportunistic pathogen, Porphyromonas gingivalis, is comparatively well-studied and known to secrete numerous extracellular matrix-targeting proteases. Among these are the gingipain family of cysteine proteases that drive periodontal disease progression and hematogenic transmission to the placenta. In this study, we demonstrate that vaginal Porphyromonas species secrete broad-acting proteases capable of freely diffusing within the cervicovaginal niche. These proteases degrade collagens that are enriched within the cervix (type I) and chorioamniotic membranes (type IV), as well as fibrinogen, which inhibits clot formation. Bioinformatic queries confirmed the absence of gingipain orthologs and identified five serine, cysteine, and metalloprotease candidates in each species. Inhibition assays revealed that each species’ proteolytic activity can be partially attributed to a secreted metalloprotease with broad substrate specificity that is distantly related to the P. gingivalis endopeptidase PepO. This characterization of virulence activities in vaginal Porphyromonas species highlights their potential to alter the homeostasis of reproductive tissues and harm human pregnancy through clotting disruption, fetal membrane weakening, and premature cervical remodeling.

Similarities in Virulence and Extended Spectrum Beta-Lactamase Gene Profiles among Cefotaxime-Resistant Escherichia coli Wastewater and Clinical Isolates

The World Health Organization has identified antibiotic resistance as one of the largest threats to human health and food security. In this study, we compared antibiotic resistance patterns between ESBL-producing Escherichia coli from human clinical diseases and cefotaxime-resistant environmental strains, as well as their potential to be pathogenic. Antibiotic susceptibility was tested amongst clinical isolates (n = 11), hospital wastewater (n = 22), and urban wastewater (n = 36, both influent and treated effluents). Multi-drug resistance predominated (>70%) among hospitalwastewater and urban wastewater influent isolates. Interestingly, isolates from clinical and urban treated effluents showed similar multi-drug resistance rates (~50%). Most hospital wastewater isolates were Phylogroup A, while clinical isolates were predominately B2, with a more diverse phylogroup population in urban wastewater. ESBL characterization of cefotaxime-resistant populations identified blaCTX-M-1 subgroup as the most common, whereby blaKPC was more associated with ceftazidime and ertapenem resistance. Whole-genome sequencing of a carbapenemase-producing hospital wastewater E. coli strain revealed plasmid-mediated blaKPC-2. Among cefotaxime-resistant populations, over 60% of clinical and 30% of treated effluent E. coli encoded three or more virulence genes exhibiting a pathogenic potential. Together, the similarity among treated effluent E. coli populations and clinical strains suggest effluents could serve as a reservoir for future multi-drug resistant E. coli clinical infections.

Response of Coastal Shewanella and Duganella Bacteria to Planktonic and Terrestrial Food Substrates

Global warming scenarios indicate that in subarctic regions, the precipitation will increase in the future. Coastal bacteria will thus receive increasing organic carbon sources from land runoff. How such changes will affect the function and taxonomic composition of coastal bacteria is poorly known. We performed a 10-day experiment with two isolated bacteria: Shewanella baltica from a seaside location and Duganella sp. from a river mouth, and provided them with a plankton and a river extract as food substrate. The bacterial growth and carbon consumption were monitored over the experimental period. Shewanella and Duganella consumed 40% and 30% of the plankton extract, respectively, while the consumption of the river extract was low for both bacteria, ∼1%. Shewanella showed the highest bacterial growth efficiency (BGE) (12%) when grown on plankton extract, while when grown on river extract, the BGE was only 1%. Duganella showed low BGE when grown on plankton extract (< 1%) and slightly higher BGE when grown on river extract (2%). The cell growth yield of Duganella was higher than that of Shewanella when grown on river extract. These results indicate that Duganella is more adapted to terrestrial organic substrates with low nutritional availability, while Shewanella is adapted to eutrophied conditions. The different growth performance of the bacteria could be traced to genomic variations. A closely related genome of Shewanella was shown to harbor genes for the sequestration of autochthonously produced carbon substrates, while Duganella contained genes for the degradation of relatively refractive terrestrial organic matter. The results may reflect the influence of environmental drivers on bacterial community composition in natural aquatic environments. Elevated inflows of terrestrial organic matter to coastal areas in subarctic regions would lead to increased occurrence of bacteria adapted to the degradation of complex terrestrial compounds with a low bioavailability.

Exploring the roles of microbes in facilitating plant adaptation to climate change

Plants benefit from their close association with soil microbes which assist in their response to abiotic and biotic stressors. Yet much of what we know about plant stress responses is based on studies where the microbial partners were uncontrolled and unknown. Under climate change, the soil microbial community will also be sensitive to and respond to abiotic and biotic stressors. Thus, facilitating plant adaptation to climate change will require a systems-based approach that accounts for the multi-dimensional nature of plant-microbe-environment interactions. In this perspective, we highlight some of the key factors influencing plant-microbe interactions under stress as well as new tools to facilitate the controlled study of their molecular complexity, such as fabricated ecosystems and synthetic communities. When paired with genomic and biochemical methods, these tools provide researchers with more precision, reproducibility, and manipulability for exploring plant-microbe-environment interactions under a changing climate.

A positive correlation between GC content and growth temperature in prokaryotes

Background

GC pairs are generally more stable than AT pairs; GC-rich genomes were proposed to be more adapted to high temperatures than AT-rich genomes. Previous studies consistently showed positive correlations between growth temperature and the GC contents of structural RNA genes. However, for the whole genome sequences and the silent sites of the codons in protein-coding genes, the relationship between GC content and growth temperature is in a long-lasting debate.

Results

With a dataset much larger than previous studies (681 bacteria and 155 archaea with completely assembled genomes), our phylogenetic comparative analyses showed positive correlations between optimal growth temperature (Topt) and GC content both in bacterial and archaeal structural RNA genes and in bacterial whole genome sequences, chromosomal sequences, plasmid sequences, core genes, and accessory genes. However, in the 155 archaea, we did not observe a significant positive correlation of Topt with whole-genome GC content (GC_w) or GC content at four-fold degenerate sites. We randomly drew 155 samples from the 681 bacteria for 1000 rounds. In most cases (> 95%), the positive correlations between Topt and genomic GC contents became statistically nonsignificant (P > 0.05). This result suggested that the small sample sizes might account for the lack of positive correlations between growth temperature and genomic GC content in the 155 archaea and the bacterial samples of previous studies. Comparing the GC content among four categories (psychrophiles/psychrotrophiles, mesophiles, thermophiles, and hyperthermophiles) also revealed a positive correlation between GC_w and growth temperature in bacteria. By including the GC_w of incompletely assembled genomes, we expanded the sample size of archaea to 303. Positive correlations between GC_w and Topt appear especially after excluding the halophilic archaea whose GC contents might be strongly shaped by intense UV radiation.

Conclusions

This study explains the previous contradictory observations and ends a long debate. Prokaryotes growing in high temperatures have higher GC contents. Thermal adaptation is one possible explanation for the positive association. Meanwhile, we propose that the elevated efficiency of DNA repair in response to heat mutagenesis might have the by-product of increasing GC content like that happens in intracellular symbionts and marine bacterioplankton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08353-7.

A Deep-Sea Bacterium Senses Blue Light via a BLUF-Dependent Pathway

Light is a ubiquitous energy source and environmental signal that broadly impacts the lifestyle of a large number of photosynthetic/nonphotosynthetic microorganisms living in the euphotic layer. However, the responses of deep-sea microbes to light are largely unknown, even though blue light is proposed to be distributed in the deep ocean. Here, we successfully cultured a novel bacterial species, named Spongiibacter nanhainus CSC3.9, from deep-sea cold seep samples by a blue light induction approach. The growth of strain CSC3.9 was obviously promoted by the illumination of blue light. We next determined BLUF (a typical blue light photoreceptor) was the most essential factor directing light sensing of strain CSC3.9 through a combined proteomic and genetic method. The function of light sensing mediated by BLUF was further confirmed by the in vitro-synthesized protein. Notably, homologs of BLUF widely existed across the marine microorganisms (containing Spongiibacter species) derived from different environments, including cold seeps. This strongly indicates that the distribution of light utilization by the nonphototrophic bacteria living in the ocean is broad and has been substantially underestimated.

IMPORTANCE Extensive studies have been conducted to explore the mechanisms of light sensing and utilization by microorganisms that live in the photic zone. Strikingly, accumulated evidence shows that light is distributed in the deep biosphere. However, the existence and process of light sensing and utilization by microbes inhabiting the deep ocean have been seldom reported. In the present study, a novel bacterial strain, Spongiibacter nanhainus CSC3.9, was enriched and purified from a deep-sea cold seep sample through a blue light induction method. Combined with genomic, proteomic, genetic, and biochemical approaches, the mechanism of this novel strain sensing blue light through a BLUF-dependent pathway was detailed. Our study provides a good model to study the mechanisms of light sensing mediated by deep-sea nonphototrophic bacteria.

Petabase-scale sequence alignment catalyses viral discovery

Public databases contain a planetary collection of nucleic acid sequences, but their systematic exploration has been inhibited by a lack of efficient methods for searching this corpus, which (at the time of writing) exceeds 20 petabases and is growing exponentially1. Here we developed a cloud computing infrastructure, Serratus, to enable ultra-high-throughput sequence alignment at the petabase scale. We searched 5.7 million biologically diverse samples (10.2 petabases) for the hallmark gene RNA-dependent RNA polymerase and identified well over 105 novel RNA viruses, thereby expanding the number of known species by roughly an order of magnitude. We characterized novel viruses related to coronaviruses, hepatitis delta virus and huge phages, respectively, and analysed their environmental reservoirs. To catalyse the ongoing revolution of viral discovery, we established a free and comprehensive database of these data and tools. Expanding the known sequence diversity of viruses can reveal the evolutionary origins of emerging pathogens and improve pathogen surveillance for the anticipation and mitigation of future pandemics.

Small investments with big returns: environmental genomic bioprospecting of microbial life

Microorganisms and their natural products are major drivers of ecological processes and industrial applications. Microbial bioprospecting has been critical for the advancement in various fields such as pharmaceuticals, sustainable industries, food security and bioremediation. Next generation sequencing has been paramount in the exploration of diverse environmental microbiomes. It presents a culture-independent approach to investigating hitherto uncultured taxa, resulting in the creation of massive sequence databases, which are available in the public domain. Genome mining searches available (meta)genomic data for target biosynthetic genes, and combined with the large-scale public data, this in-silico bioprospecting method presents an efficient and extensive way to uncover microbial bioproducts. Bioinformatic tools have progressed to a stage where we can recover genomes from the environment; these metagenome-assembled genomes present a way to understand the metabolic capacity of microorganisms in a physiological and ecological context. Environmental sampling been extensive across various ecological settings, including microbiomes with unique physicochemical properties that could influence the discovery of novel functions and metabolic pathways. Although in-silico methods cannot completely substitute in-vitro studies, the contextual information it provides is invaluable for understanding the ecological and taxonomic distribution of microbial genotypes and to form effective strategies for future microbial bioprospecting efforts.

GAL08, an Uncultivated Group of Acidobacteria, Is a Dominant Bacterial Clade in a Neutral Hot Spring

GAL08 are bacteria belonging to an uncultivated phylogenetic cluster within the phylum Acidobacteria. We detected a natural population of the GAL08 clade in sediment from a pH-neutral hot spring located in British Columbia, Canada. To shed light on the abundance and genomic potential of this clade, we collected and analyzed hot spring sediment samples over a temperature range of 24.2–79.8°C. Illumina sequencing of 16S rRNA gene amplicons and qPCR using a primer set developed specifically to detect the GAL08 16S rRNA gene revealed that absolute and relative abundances of GAL08 peaked at 65°C along three temperature gradients. Analysis of sediment collected over multiple years and locations revealed that the GAL08 group was consistently a dominant clade, comprising up to 29.2% of the microbial community based on relative read abundance and up to 4.7 × 10⁵ 16S rRNA gene copy numbers per gram of sediment based on qPCR. Using a medium quality threshold, 25 single amplified genomes (SAGs) representing these bacteria were generated from samples taken at 65 and 77°C, and seven metagenome-assembled genomes (MAGs) were reconstructed from samples collected at 45–77°C. Based on average nucleotide identity (ANI), these SAGs and MAGs represented three separate species, with an estimated average genome size of 3.17 Mb and GC content of 62.8%. Phylogenetic trees constructed from 16S rRNA gene sequences and a set of 56 concatenated phylogenetic marker genes both placed the three GAL08 bacteria as a distinct subgroup of the phylum Acidobacteria, representing a candidate order (Ca. Frugalibacteriales) within the class Blastocatellia. Metabolic reconstructions from genome data predicted a heterotrophic metabolism, with potential capability for aerobic respiration, as well as incomplete denitrification and fermentation. In laboratory cultivation efforts, GAL08 counts based on qPCR declined rapidly under atmospheric levels of oxygen but increased slightly at 1% (v/v) O₂, suggesting a microaerophilic lifestyle.

PREGO: A Literature and Data-Mining Resource to Associate Microorganisms, Biological Processes, and Environment Types

To elucidate ecosystem functioning, it is fundamental to recognize what processes occur in which environments (where) and which microorganisms carry them out (who). Here, we present PREGO, a one-stop-shop knowledge base providing such associations. PREGO combines text mining and data integration techniques to mine such what-where-who associations from data and metadata scattered in the scientific literature and in public omics repositories. Microorganisms, biological processes, and environment types are identified and mapped to ontology terms from established community resources. Analyses of comentions in text and co-occurrences in metagenomics data/metadata are performed to extract associations and a level of confidence is assigned to each of them thanks to a scoring scheme. The PREGO knowledge base contains associations for 364,508 microbial taxa, 1090 environmental types, 15,091 biological processes, and 7971 molecular functions with a total of almost 58 million associations. These associations are available through a web portal, an Application Programming Interface (API), and bulk download. By exploring environments and/or processes associated with each other or with microbes, PREGO aims to assist researchers in design and interpretation of experiments and their results. To demonstrate PREGO’s capabilities, a thorough presentation of its web interface is given along with a meta-analysis of experimental results from a lagoon-sediment study of sulfur-cycle related microbes.

Comparative Genomics of Clostridium baratii Reveals Strain-Level Diversity in Toxin Abundance

Clostridium baratii strains are rare opportunistic pathogens associated with botulism intoxication. They have been isolated from foods, soil and be carried asymptomatically or cause botulism outbreaks. Is not taxonomically related to Clostridium botulinum, but some strains are equipped with BoNT/F7 cluster. Despite their relationship with diseases, our knowledge regarding the genomic features and phylogenetic characteristics is limited. We analyzed the pangenome of C. baratii to understand the diversity and genomic features of this species. We compared existing genomes in public databases, metagenomes, and one newly sequenced strain isolated from an asymptomatic subject. The pangenome was open, indicating it comprises genetically diverse organisms. The core genome contained 28.49% of the total genes of the pangenome. Profiling virulence factors confirmed the presence of phospholipase C in some strains, a toxin capable of disrupting eukaryotic cell membranes. Furthermore, the genomic analysis indicated significant horizontal gene transfer (HGT) events as defined by the presence of prophage genomes. Seven strains were equipped with BoNT/F7 cluster. The active site was conserved in all strains, identifying a missing 7-aa region upstream of the active site in C. baratii genomes. This analysis could be important to advance our knowledge regarding opportunistic clostridia and better understand their contribution to disease.

Challenges in Bioinformatics Workflows for Processing Microbiome Omics Data at Scale

The nascent field of microbiome science is transitioning from a descriptive approach of cataloging taxa and functions present in an environment to applying multi-omics methods to investigate microbiome dynamics and function. A large number of new tools and algorithms have been designed and used for very specific purposes on samples collected by individual investigators or groups. While these developments have been quite instructive, the ability to compare microbiome data generated by many groups of researchers is impeded by the lack of standardized application of bioinformatics methods. Additionally, there are few examples of broad bioinformatics workflows that can process metagenome, metatranscriptome, metaproteome and metabolomic data at scale, and no central hub that allows processing, or provides varied omics data that are findable, accessible, interoperable and reusable (FAIR). Here, we review some of the challenges that exist in analyzing omics data within the microbiome research sphere, and provide context on how the National Microbiome Data Collaborative has adopted a standardized and open access approach to address such challenges.

Removal of antibiotic thiamphenicol by bacterium Aeromonas hydrophila HS01

Thiamphenicol (TAP) is an amphenicol antibiotic, which has a broad-spectrum inhibitory effect on both gram-positive and gram-negative bacteria. Since it is widely used in animals and aquaculture, its residues in environment may bring potential risk for human health and ecosystems. While TAP can be removed through conventional physical or chemical methods, its bioremediation using microorganisms is less studied. Here, we report the removal of TAP by a bacterial strain, Aeromonas hydrophila HS01, which can remove more than 90.0% of TAP in a living cell-dependent manner. Our results indicated that its removal efficiency can be greatly affected by the growth condition. Proteomics studies revealed a number of differentially expressed proteins of HS01 in the presence of TAP, which may play critical roles in the transportation and degradation of TAP. All these results indicate bacterial strain A. hydrophila HS01 is a new microbial resource for efficiently removing TAP, and may shed new insights in developing bioremediation approaches for TAP pollution.