Evolutionary investigations of the biosynthetic diversity in the skin microbiome using lsaBGC

BGCFlow: systematic pangenome workflow for the analysis of biosynthetic gene clusters across large genomic datasets

Genome mining is revolutionizing natural products discovery efforts. The rapid increase in available genomes demands comprehensive computational platforms to effectively extract biosynthetic knowledge encoded across bacterial pangenomes. Here, we present BGCFlow, a novel systematic workflow integrating analytics for large-scale genome mining of bacterial pangenomes. BGCFlow incorporates several genome analytics and mining tools grouped into five common stages of analysis such as: (i) data selection, (ii) functional annotation, (iii) phylogenetic analysis, (iv) genome mining, and (v) comparative analysis. Furthermore, BGCFlow provides easy configuration of different projects, parallel distribution, scheduled job monitoring, an interactive database to visualize tables, exploratory Jupyter Notebooks, and customized reports. Here, we demonstrate the application of BGCFlow by investigating the phylogenetic distribution of various biosynthetic gene clusters detected across 42 genomes of the Saccharopolyspora genus, known to produce industrially important secondary/specialized metabolites. The BGCFlow-guided analysis predicted more accurate dereplication of BGCs and guided the targeted comparative analysis of selected RiPPs. The scalable, interoperable, adaptable, re-entrant, and reproducible nature of the BGCFlow will provide an effective novel way to extract the biosynthetic knowledge from the ever-growing genomic datasets of biotechnologically relevant bacterial species.

The porcine skin microbiome exhibits broad fungal antagonism

The skin and its microbiome function to protect the host from pathogen colonization and environmental stressors. In this study, using the Wisconsin Miniature Swine™ model, we characterize the porcine skin fungal and bacterial microbiomes, identify bacterial isolates displaying antifungal activity, and use whole-genome sequencing to identify biosynthetic gene clusters encoding for secondary metabolites that may be responsible for the antagonistic effects on fungi. Through this comprehensive approach of paired microbiome sequencing with culturomics, we report the discovery of novel species of Corynebacterium and Rothia. Further, this study represents the first comprehensive evaluation of the porcine skin mycobiome and the evaluation of bacterial-fungal interactions on this surface. Several diverse bacterial isolates exhibit potent antifungal properties against opportunistic fungal pathogens in vitro. Genomic analysis of inhibitory species revealed a diverse repertoire of uncharacterized biosynthetic gene clusters suggesting a reservoir of novel chemical and biological diversity. Collectively, the porcine skin microbiome represents a potential unique source of novel antifungals.

Comprehensive analysis of genomic variation, pan-genome and biosynthetic potential of Corynebacterium glutamicum strains

Corynebacterium glutamicum is a non-pathogenic species of the Corynebacteriaceae family. It has been broadly used in industrial biotechnology for the production of valuable products. Though it is widely accepted at the industrial level, knowledge about the genomic diversity of the strains is limited. Here, we investigated the comparative genomic features of the strains and pan-genomic characteristics. We also observed phylogenetic relationships among the strains based on average nucleotide identity (ANI). We found diversity between strains at the genomic and pan-genomic levels. Less than one-third of the C. glutamicum pan-genome consists of core genes and soft-core genes. Whereas, a large number of strain-specific genes covered about half of the total pan-genome. Besides, C. glutamicum pan-genome is open and expanding, which indicates the possible addition of new gene families to the pan-genome. We also investigated the distribution of biosynthetic gene clusters (BGCs) among the strains. We discovered slight variations of BGCs at the strain level. Several BGCs with the potential to express novel bioactive secondary metabolites have been identified. Therefore, by utilizing the characteristic advantages of C. glutamicum, different strains can be potential applicants for natural drug discovery.

Phylogenomics and genetic analysis of solvent-producing Clostridium species

The genus Clostridium is a large and diverse group within the Bacillota (formerly Firmicutes), whose members can encode useful complex traits such as solvent production, gas-fermentation, and lignocellulose breakdown. We describe 270 genome sequences of solventogenic clostridia from a comprehensive industrial strain collection assembled by Professor David Jones that includes 194 C. beijerinckii, 57 C. saccharobutylicum, 4 C. saccharoperbutylacetonicum, 5 C. butyricum, 7 C. acetobutylicum, and 3 C. tetanomorphum genomes. We report methods, analyses and characterization for phylogeny, key attributes, core biosynthetic genes, secondary metabolites, plasmids, prophage/CRISPR diversity, cellulosomes and quorum sensing for the 6 species. The expanded genomic data described here will facilitate engineering of solvent-producing clostridia as well as non-model microorganisms with innately desirable traits. Sequences could be applied in conventional platform biocatalysts such as yeast or Escherichia coli for enhanced chemical production. Recently, gene sequences from this collection were used to engineer Clostridium autoethanogenum, a gas-fermenting autotrophic acetogen, for continuous acetone or isopropanol production, as well as butanol, butanoic acid, hexanol and hexanoic acid production.

Skin-associated Corynebacterium amycolatum shares cobamides

The underlying interactions that occur to maintain skin microbiome composition, function, and overall skin health are largely unknown. Often, these types of interactions are mediated by microbial metabolites. Cobamides, the vitamin B12 family of cofactors, are essential for metabolism in many bacteria, but are only synthesized by a small fraction of prokaryotes, including certain skin-associated species. Therefore, we hypothesize that cobamide sharing mediates skin community dynamics. Preliminary work predicts that several skin-associated Corynebacterium species encode de novo cobamide biosynthesis and that their abundance is associated with skin microbiome diversity. Here, we show that commensal Corynebacterium amycolatum produces cobamides and that this synthesis can be tuned by cobalt limitation. To demonstrate cobamide sharing by C. amycolatum, we employed a co-culture assay using an E. coli cobamide auxotroph and show that C. amycolatum produces sufficient cobamides to support E. coli growth, both in liquid co-culture and when separated spatially on solid medium. We also generated a C. amycolatum non-cobamide-producing strain (cob-) using UV mutagenesis that contains mutated cobamide biosynthesis genes cobK and cobO and confirm that disruption of cobamide biosynthesis abolishes support of E. coli growth through cobamide sharing. Our study provides a unique model to study metabolite sharing by microorganisms, which will be critical for understanding the fundamental interactions that occur within complex microbiomes and for developing approaches to target the human microbiota for health advances.

Genomic characterization of the C. tuberculostearicum species complex, a prominent member of the human skin microbiome

IMPORTANCE

Amplicon sequencing data combined with isolate whole genome sequencing have expanded our understanding of Corynebacterium on the skin. Healthy human skin is colonized by a diverse collection of Corynebacterium species, but Corynebacterium tuberculostearicum predominates on many skin sites. Our work supports the emerging idea that C. tuberculostearicum is a species complex encompassing several distinct species. We produced a collection of genomes that help define this complex, including a potentially new species we term Corynebacterium hallux based on a preference for sites on the feet, whole-genome average nucleotide identity, pangenomic analysis, and growth in skin-like media. This isolate collection and high-quality genome resource set the stage for developing engineered strains for both basic and translational clinical studies.

Identification of staphyloxanthin and derivates in yellow-pigmented Staphylococcus capitis subsp. capitis

Introduction

Staphylococcus capitis naturally colonizes the human skin but as an opportunistic pathogen, it can also cause biofilm-associated infections and bloodstream infections in newborns. Previously, we found that two strains from the subspecies S. capitis subsp. capitis produce yellow carotenoids despite the initial species description, reporting this subspecies as non-pigmented. In Staphylococcus aureus, the golden pigment staphyloxanthin is an important virulence factor, protecting cells against reactive oxygen species and modulating membrane fluidity.

Methods

In this study, we used two pigmented (DSM 111179 and DSM 113836) and two non-pigmented S. capitis subsp. capitis strains (DSM 20326T and DSM 31028) to identify the pigment, determine conditions under which pigment-production occurs and investigate whether pigmented strains show increased resistance to ROS and temperature stress.

Results

We found that the non-pigmented strains remained colorless regardless of the type of medium, whereas intensity of pigmentation in the two pigmented strains increased under low nutrient conditions and with longer incubation times. We were able to detect and identify staphyloxanthin and its derivates in the two pigmented strains but found that methanol cell extracts from all four strains showed ROS scavenging activity regardless of staphyloxanthin production. Increased survival to cold temperatures (-20°C) was detected in the two pigmented strains only after long-term storage compared to the non-pigmented strains.

Conclusion

The identification of staphyloxanthin in S. capitis is of clinical relevance and could be used, in the same way as in S. aureus, as a possible target for anti-virulence drug design.

zol & fai: large-scale targeted detection and evolutionary investigation of gene clusters

Many universally and conditionally important genes are genomically aggregated within clusters. Here, we introduce fai and zol, which together enable large-scale comparative analysis of different types of gene clusters and mobile-genetic elements (MGEs), such as biosynthetic gene clusters (BGCs) or viruses. Fundamentally, they overcome a current bottleneck to reliably perform comprehensive orthology inference at large scale across broad taxonomic contexts and thousands of genomes. First, fai allows the identification of orthologous or homologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of protein-encoding ortholog groups for individual genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of statistics for each inferred ortholog group. These programs are showcased through application to: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-genetic insights of two common BGCs in a fungal species, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Metabolic capabilities are highly conserved among human nasal-associated Corynebacterium species in pangenomic analyses

Corynebacterium species are globally ubiquitous in human nasal microbiota across the lifespan. Moreover, nasal microbiota profiles typified by higher relative abundances of Corynebacterium are often positively associated with health. Among the most common human nasal Corynebacterium species are C. propinquum, C. pseudodiphtheriticum, C. accolens, and C. tuberculostearicum. Based on the prevalence of these species, at least two likely coexist in the nasal microbiota of 82% of adults. To gain insight into the functions of these four species, we identified genomic, phylogenomic, and pangenomic properties and estimated the functional protein repertoire and metabolic capabilities of 87 distinct human nasal Corynebacterium strain genomes: 31 from Botswana and 56 from the U.S. C. pseudodiphtheriticum had geographically distinct clades consistent with localized strain circulation, whereas some strains from the other species had wide geographic distribution across Africa and North America. All four species had similar genomic and pangenomic structures. Gene clusters assigned to all COG metabolic categories were overrepresented in the persistent (core) compared to the accessory genome of each species indicating limited strain-level variability in metabolic capacity. Moreover, core metabolic capabilities were highly conserved among the four species indicating limited species-level metabolic variation. Strikingly, strains in the U.S. clade of C. pseudodiphtheriticum lacked genes for assimilatory sulfate reduction present in the Botswanan clade and in the other studied species, indicating a recent, geographically related loss of assimilatory sulfate reduction. Overall, the minimal species and strain variability in metabolic capacity implies coexisting strains might have limited ability to occupy distinct metabolic niches.