Variation in genomic traits of microbial communities among ecosystems

Free-living bacteria in nutrient limited environments often exhibit traits which may reduce the cost of reproduction, such as smaller genome size, low GC content, and fewer sigma (σ) factor and 16S rRNA gene copies. Despite the potential utility of these traits to detect relationships between microbial communities and ecosystem-scale properties, few studies have assessed these traits on a community-scale. Here, we analyzed these traits from publicly available metagenomes derived from marine, soil, host-associated, and thermophilic communities. In marine and thermophilic communities, genome size and GC content declined in parallel, consistent with genomic streamlining, with GC content in thermophilic communities generally higher than in marine systems. In contrast, soil communities averaging smaller genomes featured higher GC content and were often from low-carbon environments, suggesting unique selection pressures in soil bacteria. The abundance of specific σ-factors varied with average genome size and ecosystem type. In oceans, abundance of fliA, a σ-factor controlling flagella biosynthesis, was positively correlated with community average genome size – reflecting known trade-offs between nutrient conservation and chemotaxis. In soils, a high abundance of the stress response σ-factor gene rpoS was associated with smaller average genome size and often located in harsh and/or carbon-limited environments – a result which tracks features observed in culture and indicates an increased capacity for stress response in nutrient-poor soils. This work shows how ecosystem-specific constraints are associated with trade-offs which are embedded in the genomic features of bacteria in microbial communities, and which can be detected at the community level, highlighting the importance of genomic features in microbial community analysis.

A general model of the dynamics of genes with dual σ factor preference

Escherichia coli uses the ability of σ factors to recognize specific DNA sequences in order to quickly control large gene cohorts. While most genes respond to only one s factor, approximately 5% have dual s factor preference. The ones in significant numbers are ‘σ 70+38 genes’, responsive to σ 70 , which controls housekeeping genes, as well as to σ 38 , which activates genes during stationary growth and stresses. We show that σ 70+38 genes are almost as upregulated in stationary growth as genes responsive to σ 38 alone. Also, their response strengths to σ 38 are predictable from their promoter sequences. Next, we propose a sequence- and σ 38 level-dependent, analytical model of σ 70+38 genes applicable in the exponential, stationary, and in the transition period between the two growth phases. Finally, we propose a general model, applicable to other σ factors as well. This model can guide the design of synthetic circuits with sequence-dependent sensitivity and plasticity to transitions between the exponential and stationary growth phases.

Selection for constrained peptides that bind to a single target protein

Peptide secondary metabolites are common in nature and have diverse pharmacologically-relevant functions, from antibiotics to cross-kingdom signaling. Here, we present a method to design large libraries of modified peptides in Escherichia coli and screen them in vivo to identify those that bind to a single target-of-interest. Constrained peptide scaffolds were produced using modified enzymes gleaned from microbial RiPP (ribosomally synthesized and post-translationally modified peptide) pathways and diversified to build large libraries. The binding of a RiPP to a protein target leads to the intein-catalyzed release of an RNA polymerase σ factor, which drives the expression of selectable markers. As a proof-of-concept, a selection was performed for binding to the SARS-CoV-2 Spike receptor binding domain. A 1625 Da constrained peptide (AMK-1057) was found that binds with similar affinity (990 ± 5 nM) as an ACE2-derived peptide. This demonstrates a generalizable method to identify constrained peptides that adhere to a single protein target, as a step towards “molecular glues” for therapeutics and diagnostics.

Listeria monocytogenes 10403S Alternative Sigma-54 Factor σL Has a Negative Role on Survival Ability Under Bile Exposure

Listeria monocytogenes is a Gram-positive bacterium causing listeriosis in animals and humans. To initiate a foodborne infection, L. monocytogenes has to pass through the host gastrointestinal tract (GIT). In this study, we evaluated survival abilities of L. monocytogenes 10403S wild type (WT) and its isogenic mutants in alternative sigma (σ) factor genes (i.e., sigB, sigC, sigH, and sigL) under simulated gastric, duodenal, and bile fluids. Within 10min of exposures, only bile fluid was able to significantly reduce survival ability of L. monocytogenes WT by 2 logs CFU/ml. Loss of sigL showed the greatest bile resistance among 16 strains tested, p<0.0001, (i.e., WT, four single alternative σ factor mutants, six double mutants, four triple mutants, and one quadruple mutant). To further investigate the role of σL in bile response, RNA-seq was conducted to compare the transcriptional profiles among L. monocytogenes 10403S ΔBCH triple mutant (lacking sigB, sigC, and sigH genes; expressing housekeeping σA and σL) and ΔBCHL quadruple mutant (lacking all alternative sigma factor genes; expressing only σA) strains under BHI and 1% bile conditions. A total of 216 and 176 differentially expressed genes (DEGs) were identified in BHI and bile, respectively. We confirmed that mpt operon was shown to be strongly activated by σL. Interestingly, more than 80% of DEGs were found to be negatively regulated in the presence of σL. This includes PrfA regulon and its mediated genes (i.e., hly, hpt, inlB, clpP, clpE, groL, and inlC) which were downregulated in response to bile in the presence of σL. This result suggests the potential negative role of σL on bile survival, and the roles of σL and σB might be in a seesaw model prior to host cell invasion.

Transcriptional organization and regulation of the Pseudomonas putida flagellar system

A single region of the Pseudomonas putida genome, designated the flagellar cluster, includes 59 genes potentially involved in the biogenesis and function of the flagellar system. Here we combine bioinformatics and in vivo gene expression analyses to clarify the transcriptional organization and regulation of the flagellar genes in the cluster. We have identified eleven flagellar operons and characterized twenty-one primary and internal promoter regions. Our results indicate that synthesis of the flagellar apparatus and core chemotaxis machinery is regulated by a three-tier cascade in which fleQ is the sole Class I gene, standing at the top of the transcriptional hierarchy. FleQ- and σ54-dependent Class II genes encode most components of the flagellar structure, part of the chemotaxis machinery and multiple regulatory elements, including the flagellar σ factor FliA. FliA activation of Class III genes enables synthesis of the filament, one stator complex and completion of the chemotaxis apparatus. Accessory regulatory proteins and an intricate operon architecture add complexity to the regulation by providing feedback and feed-forward loops to the main circuit. Because of the high conservation of the gene arrangement and promoter motifs, we believe that the regulatory circuit presented here may also apply to other environmental Pseudomonas.

A Codon Constrained Method for Both Eliminating and Creating Intragenic Bacterial Promoters

Future applications of synthetic biology will require refactored genetic sequences devoid of internal regulatory elements within coding sequences. These regulatory elements include cryptic and intragenic promoters which may constitute up to a third of predicted Escherichia coli promoters. Promoter activity is dependent on the structural interaction of core bases with a σ factor. Rational engineering can be used to alter key promoter element nucleotides interacting with σ factors and eliminate downstream transcriptional activity. In this paper, we present COdon Restrained Promoter SilEncing (CORPSE), a system for removing intragenic promoters. CORPSE exploits the DNA-σ factor structural relationship to disrupt σ70 promoters embedded within gene coding sequences, with a minimum of synonymous codon changes. Additionally, we present an inverted CORPSE system, iCORPSE, which can create highly active promoters within a gene sequence while not perturbing the function of the modified gene.

A novel mycothiol-dependent thiol–disulfide reductase in Corynebacterium glutamicum involving oxidative stress resistance

Abstractncgl2478 gene from Corynebacterium glutamicum encodes a thiol–disulfide oxidoreductase enzyme annotated as dithiol–disulfide isomerase DsbA. It preserves a Cys–Pro–Phe–Cys active-site motif, which is presumed to be an exclusive characteristic of the novel DsbA–mycoredoxin 1 (Mrx1) cluster. However, the real mode of action, the nature of the electron donor pathway and biological functions of NCgl2478 in C. glutamicum have remained enigmatic so far. Herein, we report that NCgl2478 plays an important role in stress resistance. Deletion of the ncgl2478 gene increases the size of growth inhibition zones. The ncgl2478 expression is induced in the stress-responsive extra-cytoplasmic function-sigma (ECF-σ) factor SigH-dependent manner by stress. It receives electrons preferentially from the mycothiol (MSH)/mycothione reductase (Mtr)/NADPH pathway. Further, NCgl2478 reduces S-mycothiolated mixed disulfides and intramolecular disulfides via a monothiol–disulfide and a dithiol–disulfide exchange mechanism, respectively. NCgl2478 lacks oxidase activity; kinetic properties of its demycothiolation are different from those of Mrx1. Site-directed mutagenesis confirms Cys24 is the resolving Cys residue, while Cys21 is the nucleophilic cysteine that is oxidized to a sulfenic acid and then forms an intramolecular disulfide bond with Cys24 or a mixed disulfide with MSH under oxidative stress. In conclusion, our study presents the first evidence that NCgl2478 protects against various stresses by acting as an MSH-dependent thiol–disulfide reductase, belonging to a novel DsbA–Mrx1 cluster.

Soil, ocean, hot spring, and host-associated environments reveal unique selection pressures on genomic features of bacteria in microbial communities

Free-living bacteria in nutrient limited environments often exhibit small genomes which curb the cost of reproduction - a phenomenon known as genomic streamlining. Streamlining has been associated with a suite of traits such as reduced GC content, fewer 16S rRNA copies, and a lower abundance of regulatory genes, such as sigma (σ)-factors. Here, we analyzed these traits from 116 publicly available metagenomes derived from marine, soil, host associated, and thermophilic communities. In marine and thermophilic communities, genome size and GC content declined in parallel, but GC content was higher in thermophilic communities. In soils, the relationship between genome size and GC content was negative, suggesting a different selection pressure on genome size and GC content in soil bacteria. The abundance of σ-factors varied with average genome size, ecosystem type, and the specific functions regulated by the sigma factor. In marine environments, housekeeping and heat-shock σ-factor genes (rpoD and rpoH respectively) increased as genome size declined, and σ-factor responsible for flagella biosynthesis (fliA) decreased, suggesting a trade-off between nutrient conservation and chemotaxis. In soils, a high abundance of fliA and the stress response σ-factor gene (rpoS) was associated with smaller average genome size and often located in harsh and/or carbon-limited environments such as deserts or agricultural fields - suggesting an increased capacity for stress response and mobility in nutrient-poor soils. This work showcases how ecosystem-specific environmental constraints force trade-offs which are then embedded in the genomic features of bacteria in microbial communities, specifically genome size, GC content, and regulatory genes, and further highlights the importance of considering these features in microbial community analysis.

Alternative σ Factors Regulate Overlapping as Well as Distinct Stress Response and Metabolic Functions in Listeria monocytogenes under Stationary Phase Stress Condition

Listeria monocytogenes can regulate and fine-tune gene expression, to adapt to diverse stress conditions encountered during foodborne transmission. To further understand the contributions of alternative sigma (σ) factors to the regulation of L. monocytogenes gene expression, RNA-Seq was performed on L. monocytogenes strain 10403S and five isogenic mutants (four strains bearing in-frame null mutations in three out of four alternative σ factor genes, ΔCHL, ΔBHL, ΔBCL, and ΔBCH, and one strain bearing null mutations in all four genes, ΔBCHL), grown to stationary phase. Our data showed that 184, 35, 34, and 20 genes were positively regulated by σB, σL, σH, and σC (posterior probability > 0.9 and Fold Change (FC) > 5.0), respectively. Moreover, σB-dependent genes showed the highest FC (based on comparisons between the ΔCHL and the ΔBCHL strain), with 44 genes showing an FC > 100; only four σL-dependent, and no σH- or σC-dependent genes showed FC >100. While σB-regulated genes identified in this study are involved in stress-associated functions and metabolic pathways, σL appears to largely regulate genes involved in a few specific metabolic pathways, including positive regulation of operons encoding phosphoenolpyruvate (PEP)-dependent phosphotransferase systems (PTSs). Overall, our data show that (i) σB and σL directly and indirectly regulate genes involved in several energy metabolism-related functions; (ii) alternative σ factors are involved in complex regulatory networks and appear to have epistatic effects in stationary phase cells; and (iii) σB regulates multiple stress response pathways, while σL and σH positively regulate a smaller number of specific pathways.