The primary transcriptome of the major human pathogen Helicobacter pylori

Molecular insights into the fine-tuning of pH-dependent ArsR-mediated regulation of the SabA adhesin in Helicobacter pylori

Adaptation to variations in pH is crucial for the ability of Helicobacter pylori to persist in the human stomach. The acid responsive two-component system ArsRS, constitutes the global regulon that responds to acidic conditions, but molecular details of how transcription is affected by the ArsR response regulator remains poorly understood. Using a combination of DNA-binding studies, in vitro transcription assays, and H. pylori mutants, we demonstrate that phosphorylated ArsR (ArsR-P) forms an active protein complex that binds DNA with high specificity in order to affect transcription. Our data showed that DNA topology is key for DNA binding. We found that AT-rich DNA sequences direct ArsR-P to specific sites and that DNA-bending proteins are important for the effect of ArsR-P on transcription regulation. The repression of sabA transcription is mediated by ArsR-P with the support of Hup and is affected by simple sequence repeats located upstream of the sabA promoter. Here stochastic events clearly contribute to the fine-tuning of pH-dependent gene regulation. Our results reveal important molecular aspects for how ArsR-P acts to repress transcription in response to acidic conditions. Such transcriptional control likely mediates shifts in bacterial positioning in the gastric mucus layer.

A simple BLASTn-based approach generates novel insights into the regulation and biological function of type I toxin-antitoxins

Bacterial chromosomal type I toxin-antitoxin systems consist of a small protein, typically under 60 amino acids, and a small RNA (sRNA) that represses toxin translation. These gene pairs have gained attention over the last decade for their contribution to antibiotic persistence and phage tolerance in bacteria. However, biological functions for many remain elusive as gene deletions often fail to produce an observable phenotype. For many pairs, it is still unknown when the toxin and/or antitoxin gene are natively expressed within the bacterium. We examined sequence conservation of three type I toxin-antitoxin systems, tisB/istR-1, shoB/ohsC, and zor/orz, in over 2,000 Escherichia coli strains, including pathogenic and commensal isolates. Using our custom database, we found that these gene pairs are widespread across E. coli and have expression potential via BLASTn. We identified an alternative, dominant sequence variant of TisB and confirmed that it is toxic upon overproduction. Additionally, analyses revealed a highly conserved sequence in the zorO mRNA untranslated region that is required for full toxicity. We further noted that over 30% of E. coli genomes contain an orz antitoxin gene only and confirmed its expression in a representative strain: the first confirmed report of a type I antitoxin without its cognate toxin. Our results add to our understanding of these systems, and our methodology is applicable for other type I loci to identify critical regulatory and functional features.IMPORTANCEChromosomal type I toxin-antitoxins are a class of genes that have gained increasing attention over the last decade for their roles in antibiotic persistence which may contribute to therapeutic failures. However, the control of many of these genes and when they function have remained elusive. We demonstrate that a simple genetic conservation-based approach utilizing free, publicly available data yields known and novel insights into the regulation and function of three chromosomal type I toxin-antitoxins in Escherichia coli. This study also provides a framework for how this approach could be applied to other genes of interest.

Type I toxin-antitoxin systems in bacteria: from regulation to biological functions

Toxin-antitoxin systems are ubiquitous in the prokaryotic world and widely distributed among chromosomes and mobile genetic elements. Several different toxin-antitoxin system types exist, but what they all have in common is that toxin activity is prevented by the cognate antitoxin. In type I toxin-antitoxin systems, toxin production is controlled by an RNA antitoxin and by structural features inherent to the toxin messenger RNA. Most type I toxins are small membrane proteins that display a variety of cellular effects. While originally discovered as modules that stabilize plasmids, chromosomal type I toxin-antitoxin systems may also stabilize prophages, or serve important functions upon certain stress conditions and contribute to population-wide survival strategies. Here, we will describe the intricate RNA-based regulation of type I toxin-antitoxin systems and discuss their potential biological functions.

(Re)-definition of the holo- and apo-Fur direct regulons of Helicobacter pylori

Iron homeostasis is a critical process for living organisms because this metal is an essential co-factor for fundamental biochemical activities, like energy production and detoxification, albeit its excess quickly leads to cell intoxication. The protein Fur (ferric uptake regulator) controls iron homeostasis in bacteria by switching from its apo- to holo-form as a function of the cytoplasmic level of ferrous ions, thereby modulating gene expression. The Helicobacter pylori HpFur protein has the rare ability to operate as a transcriptional commutator; apo- and holo-HpFur function as two different repressors with distinct DNA binding recognition properties for specific sets of target genes. Although the regulation of apo- and holo-HpFur in this bacterium has been extensively investigated, we propose a genome-wide redefinition of holo-HpFur direct regulon in H. pylori by integration of RNA-seq and ChIP-seq data, and a large extension of the apo-HpFur direct regulon. We show that in response to iron availability, new coding sequences, non-coding RNAs, toxin-antitoxin systems, and transcripts within open reading frames are directly regulated by apo- or holo-HpFur. These new targets and the more thorough validation and deeper characterization of those already known provide a complete and updated picture of the direct regulons of this two-faced transcriptional regulator.

Helicobacter pylori cheV1 mutants recover semisolid agar migration due to loss of a previously uncharacterized Type IV filament membrane alignment complex homolog

The bacterial chemotaxis system is a well-understood signaling pathway that promotes bacterial success. Chemotaxis systems comprise chemoreceptors and the CheA kinase, linked by CheW or CheV scaffold proteins. Scaffold proteins provide connections between chemoreceptors and CheA and also between chemoreceptors to create macromolecular arrays. Chemotaxis is required for host colonization by many microbes, including the stomach pathogen Helicobacter pylori. This bacterium builds chemoreceptor-CheA contacts with two distinct scaffold proteins, CheW and CheV1. H. pylori cheW or cheV1 deletion mutants both lose chemoreceptor array formation, but show differing semisolid agar chemotaxis assay behaviors: ∆cheW mutants exhibit total migration failure, whereas ∆cheV1::cat mutants display a 50% reduction. On investigating these varied responses, we found that both mutants initially struggle with migration. However, over time, ∆cheV1::cat mutants develop a stable, enhanced migration capability, termed "migration-able" (Mig+). Whole-genome sequencing analysis of four distinct ∆cheV1::cat Mig+ strains identified single-nucleotide polymorphisms (SNPs) in hpg27_252 (hp0273) that were predicted to truncate the encoded protein. Computational analysis of the hpg27_252-encoded protein revealed it encoded a hypothetical protein that was a remote homolog of the PilO Type IV filament membrane alignment complex protein. Although H. pylori lacks Type IV filaments, our analysis showed it retains an operon of genes for homologs of PilO, PilN, and PilM. Deleting hpg27_252 in the ∆cheV1::cat or wild type strain resulted in enhanced migration in semisolid agar. Our study thus reveals that while cheV1 mutants initially have significant migration defects, they can recover the migration ability through genetic suppressors, highlighting a complex regulatory mechanism in bacterial migration.

IMPORTANCE

Chemotactic motility, present in over half of bacteria, depends on chemotaxis signaling systems comprising receptors, kinases, and scaffold proteins. In Helicobacter pylori, a stomach pathogen, chemotaxis is crucial for colonization, with CheV1 and CheW as key scaffold proteins. While both scaffolds are essential for building chemoreceptor complexes, their roles vary in other assays. Our research reexamines cheV1 mutants' behavior in semisolid agar, a standard chemotaxis test. Initially, cheV1 mutants exhibited defects similar to those of cheW mutants, but they evolved genetic suppressors that enhanced migration. These suppressors involve mutations in a previously uncharacterized gene, unknown in motility behavior. Our findings highlight the significant chemotaxis defects in cheV1 mutants and identify new elements influencing bacterial motility.

Transcriptome fine-mapping in Fusobacterium nucleatum reveals FoxJ, a new σE-dependent small RNA with unusual mRNA activation activity

The oral commensal Fusobacterium nucleatum can spread to extra-oral sites, where it is associated with diverse pathologies, including pre-term birth and cancer. Due to the evolutionary distance of F. nucleatum to other model bacteria, we lack a deeper understanding of the RNA regulatory networks that allow this bacterium to adapt to its various niches. As a first step in that direction, we recently showed that F. nucleatum harbors a global stress response governed by the extracytoplasmic function sigma factor, σE, which displays a striking functional conservation with Proteobacteria and includes a noncoding arm in the form of a regulatory small RNA (sRNA), FoxI. To search for putative additional σE-dependent sRNAs, we comprehensively mapped the 5' and 3' ends of transcripts in the model strain ATCC 23726. This enabled the discovery of FoxJ, a ~156-nucleotide sRNA previously misannotated as the 5' untranslated region (UTR) of ylmH. FoxJ is tightly controlled by σE and activated by the same stress conditions as is FoxI. Both sRNAs act as mRNA repressors of the abundant porin FomA, but FoxJ also regulates genes that are distinct from the target suite of FoxI. Moreover, FoxJ differs from other σE-dependent sRNAs in that it also positively regulates genes at the post-transcriptional level. We provide preliminary evidence for a new mode of sRNA-mediated mRNA activation, which involves the targeting of intra-operonic terminators. Overall, our study provides an important resource through the comprehensive annotation of 5' and 3' UTRs in F. nucleatum and expands our understanding of the σE response in this evolutionarily distant bacterium.IMPORTANCEThe oral microbe Fusobacterium nucleatum can colonize secondary sites, including cancer tissue, and likely deploys complex regulatory systems to adapt to these new environments. These systems are largely unknown, partly due to the phylogenetic distance of F. nucleatum to other model organisms. Previously, we identified a global stress response mediated by σE that displays functional conservation with the envelope stress response in Proteobacteria, comprising a coding and noncoding regulatory arm. Through global identification of transcriptional start and stop sites, we uncovered the small RNA (sRNA) FoxJ as a novel component of the noncoding arm of the σE response in F. nucleatum. Together with its companion sRNA FoxI, FoxJ post-transcriptionally modulates the synthesis of envelope proteins, revealing a conserved function for σE-dependent sRNAs between Fusobacteriota and Proteobacteria. Moreover, FoxJ activates the gene expression for several targets, which is a mode of regulation previously unseen in the noncoding arm of the σE response.

MLDSPP: Bacterial Promoter Prediction Tool Using DNA Structural Properties with Machine Learning and Explainable AI

Bacterial promoters play a crucial role in gene expression by serving as docking sites for the transcription initiation machinery. However, accurately identifying promoter regions in bacterial genomes remains a challenge due to their diverse architecture and variations. In this study, we propose MLDSPP (Machine Learning and Duplex Stability based Promoter prediction in Prokaryotes), a machine learning-based promoter prediction tool, to comprehensively screen bacterial promoter regions in 12 diverse genomes. We leveraged biologically relevant and informative DNA structural properties, such as DNA duplex stability and base stacking, and state-of-the-art machine learning (ML) strategies to gain insights into promoter characteristics. We evaluated several machine learning models, including Support Vector Machines, Random Forests, and XGBoost, and assessed their performance using accuracy, precision, recall, specificity, F1 score, and MCC metrics. Our findings reveal that XGBoost outperformed other models and current state-of-the-art promoter prediction tools, namely Sigma70pred and iPromoter2L, achieving F1-scores >95% in most systems. Significantly, the use of one-hot encoding for representing nucleotide sequences complements these structural features, enhancing our XGBoost model's predictive capabilities. To address the challenge of model interpretability, we incorporated explainable AI techniques using Shapley values. This enhancement allows for a better understanding and interpretation of the predictions of our model. In conclusion, our study presents MLDSPP as a novel, generic tool for predicting promoter regions in bacteria, utilizing original downstream sequences as nonpromoter controls. This tool has the potential to significantly advance the field of bacterial genomics and contribute to our understanding of gene regulation in diverse bacterial systems.

Genome-wide annotation of transcript boundaries using bacterial Rend-seq datasets

Accurate annotation to single-nucleotide resolution of the transcribed regions in genomes is key to optimally analyse RNA-seq data, understand regulatory events and for the design of experiments. However, currently most genome annotations provided by GenBank generally lack information about untranslated regions. Additionally, information regarding genomic locations of non-coding RNAs, such as sRNAs, or anti-sense RNAs is frequently missing. To provide such information, diverse RNA-seq technologies, such as Rend-seq, have been developed and applied to many bacterial species. However, incorporating this vast amount of information into annotation files has been limited and is bioinformatically challenging, resulting in UTRs and other non-coding elements being overlooked or misrepresented. To overcome this problem, we present pyRAP (python Rend-seq Annotation Pipeline), a software package that analyses Rend-seq datasets to accurately resolve transcript boundaries genome-wide. We report the use of pyRAP to find novel transcripts, transcript isoforms, and RNase-dependent sRNA processing events. In Bacillus subtilis we uncovered 63 novel transcripts and provide genomic coordinates with single-nucleotide resolution for 2218 5'UTRs, 1864 3'UTRs and 161 non-coding RNAs. In Escherichia coli, we report 117 novel transcripts, 2429 5'UTRs, 1619 3'UTRs and 91 non-coding RNAs, and in Staphylococcus aureus, 16 novel transcripts, 664 5'UTRs, 696 3'UTRs, and 81 non-coding RNAs. Finally, we use pyRAP to produce updated annotation files for B. subtilis 168, E. coli K-12 MG1655, and S. aureus 8325 for use in the wider microbial genomics research community.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Environmental pollutants and phosphoinositide signaling in autoimmunity

Environmental pollution stands as one of the most critical challenges affecting human health, with an estimated mortality rate linked to pollution-induced non-communicable diseases projected to range from 20% to 25%. These pollutants not only disrupt immune responses but can also trigger immunotoxicity. Phosphoinositide signaling, a pivotal regulator of immune responses, plays a central role in the development of autoimmune diseases and exhibits high sensitivity to environmental stressors. Among these stressors, environmental pollutants have become increasingly prevalent in our society, contributing to the initiation and exacerbation of autoimmune conditions. In this review, we summarize the intricate interplay between phosphoinositide signaling and autoimmune diseases within the context of environmental pollutants and contaminants. We provide an up-to-date overview of stress-induced phosphoinositide signaling, discuss 14 selected examples categorized into three groups of environmental pollutants and their connections to immune diseases, and shed light on the associated phosphoinositide signaling pathways. Through these discussions, this review advances our understanding of how phosphoinositide signaling influences the coordinated immune response to environmental stressors at a biological level. Furthermore, it offers valuable insights into potential research directions and therapeutic targets aimed at mitigating the impact of environmental pollutants on the pathogenesis of autoimmune diseases. SYNOPSIS: Phosphoinositide signaling at the intersection of environmental pollutants and autoimmunity provides novel insights for managing autoimmune diseases aggravated by pollutants.

Refining the transcriptional landscapes for distinct clades of virulent phages infecting Pseudomonas aeruginosa

The introduction of high-throughput sequencing has resulted in a surge of available bacteriophage genomes, unveiling their tremendous genomic diversity. However, our current understanding of the complex transcriptional mechanisms that dictate their gene expression during infection is limited to a handful of model phages. Here, we applied ONT-cappable-seq to reveal the transcriptional architecture of six different clades of virulent phages infecting Pseudomonas aeruginosa. This long-read microbial transcriptomics approach is tailored to globally map transcription start and termination sites, transcription units, and putative RNA-based regulators on dense phage genomes. Specifically, the full-length transcriptomes of LUZ19, LUZ24, 14-1, YuA, PAK_P3, and giant phage phiKZ during early, middle, and late infection were collectively charted. Beyond pinpointing traditional promoter and terminator elements and transcription units, these transcriptional profiles provide insights in transcriptional attenuation and splicing events and allow straightforward validation of Group I intron activity. In addition, ONT-cappable-seq data can guide genome-wide discovery of novel regulatory element candidates, including noncoding RNAs and riboswitches. This work substantially expands the number of annotated phage-encoded transcriptional elements identified to date, shedding light on the intricate and diverse gene expression regulation mechanisms in Pseudomonas phages, which can ultimately be sourced as tools for biotechnological applications in phage and bacterial engineering.

The global RNA-RNA interactome of Klebsiella pneumoniae unveils a small RNA regulator of cell division

The ubiquitous RNA chaperone Hfq is involved in the regulation of key biological processes in many species across the bacterial kingdom. In the opportunistic human pathogen Klebsiella pneumoniae, deletion of the hfq gene affects the global transcriptome, virulence, and stress resistance; however, the ligands of the major RNA-binding protein in this species have remained elusive. In this study, we have combined transcriptomic, co-immunoprecipitation, and global RNA interactome analyses to compile an inventory of conserved and species-specific RNAs bound by Hfq and to monitor Hfq-mediated RNA-RNA interactions. In addition to dozens of RNA-RNA pairs, our study revealed an Hfq-dependent small regulatory RNA (sRNA), DinR, which is processed from the 3' terminal portion of dinI mRNA. Transcription of dinI is controlled by the master regulator of the SOS response, LexA. As DinR accumulates in K. pneumoniae in response to DNA damage, the sRNA represses translation of the ftsZ transcript by occupation of the ribosome binding site. Ectopic overexpression of DinR causes depletion of ftsZ mRNA and inhibition of cell division, while deletion of dinR antagonizes cell elongation in the presence of DNA damage. Collectively, our work highlights the important role of RNA-based gene regulation in K. pneumoniae and uncovers the central role of DinR in LexA-controlled division inhibition during the SOS response.

Exploring the transcriptional landscape of phage-host interactions using novel high-throughput approaches

In the last decade, powerful high-throughput sequencing approaches have emerged to analyse microbial transcriptomes at a global scale. However, to date, applications of these approaches to microbial viruses such as phages remain scarce. Tailoring these techniques to virus-infected bacteria promises to obtain a detailed picture of the underexplored RNA biology and molecular processes during infection. In addition, transcriptome study of stress and perturbations induced by phages in their infected bacterial hosts is likely to reveal new fundamental mechanisms of bacterial metabolism and gene regulation. Here, we provide references and blueprints to implement emerging transcriptomic approaches towards addressing transcriptome architecture, RNA-RNA and RNA-protein interactions, RNA modifications, structures and heterogeneity of transcription profiles in infected cells that will provide guides for future directions in phage-centric therapeutic applications and microbial synthetic biology.

The shapeshifting Helicobacter pylori: From a corkscrew to a ball

There is growing evidence that bacterial morphology is closely related to their lifestyle. The helical Helicobacter pylori relies on its unique shape for survival and efficient colonization of the human stomach. Yet, they have been observed to transform into another distinctive morphology, the spherical coccoid. Despite being hypothesized to be involved in the persistence and transmission of this species, years of effort in deciphering the roles of the coccoid form remain fruitless since contrasting observations regarding its lifestyle were reported. Here, we discuss the two forms of H. pylori with a focus on the coccoid form, the molecular mechanism behind its morphological transformation, and experimental approaches to further develop our understanding of this phenomenon. We also propose a putative mechanism of the coccoid formation in H. pylori through induction of a type-I toxin-antitoxin (TA) system recently shown to influence the morphology of this species.

Unravelling the transcriptome response of Enterobacter sp. S-33 under varying temperature

Enterobacter genus includes the bacteria occupying every aspect of environment, however, their adaptability at varying temperature is not clear. In the present study, we analyzed the transcriptome response of Enterobacter sp. S-33 and their functional genes under various temperatures (30-45 ℃) that were expressed and accumulated in cells under temperature-stress. During a temperature shift from 37 to 45 ℃, 165 genes showed differential expression including 112 up-regulated and 53 down-regulated. In particular, heat-shock genes such as CspA, 16 kDa heat shock protein A/B and transcriptional regulators such as LysR, TetR, and LuxR were differentially expressed, indicating the role of complex molecular mechanism of Enterobacter adapting to temperature stress. Similarly, genes associated to signal transduction, ABC transporters, iron homeostasis, and quorum sensing were also induced. The Gene ontology enrichment analysis of differentially expressed genes (DEGs) were categorized into "transmembrane transport", "tRNA binding", "hydrogen sulfide biosynthetic process" and "sulfate assimilation" terms. In addition, Kyoto Encyclopedia of Genes and Genomes pathways showed that ABC transporter as well as quorum sensing pathways were significantly enriched. Overall, current study has contributed to explore the adaptive molecular mechanisms of Enterobacter spp. upon temperature change, which further opens new avenues for future in-depth functional studies.

The Helicobacter pylori methylome is acid-responsive due to regulation by the two-component system ArsRS and the type I DNA methyltransferase HsdM1 (HP0463)

In addition to its role in genome protection, DNA methylation can regulate gene expression. In this study, we characterized the impact of acidity, phase variation, and the ArsRS TCS on the expression of the Type I m6A DNA methyltransferase HsdM1 (HP0463) of Helicobacter pylori 26695 and their subsequent effects on the methylome. Transcription of hsdM1 increases at least fourfold in the absence of the sensory histidine kinase ArsS, the major acid-sensing protein of H. pylori. hsdM1 exists in the phase-variable operon hsdR1-hsdM1. Phase-locking hsdR1 (HP0464), the restriction endonuclease gene, has significant impacts on the transcription of hsdM1. To determine the impacts of methyltransferase transcription patterns on the methylome, we conducted methylome sequencing on samples cultured at pH 7 or pH 5. We found differentially methylated motifs between these growth conditions and that deletions of arsS and/or hsdM1 interfere with the epigenetic acid response. Deletion of arsS leads to altered activity of HsdM1 and multiple other methyltransferases under both pH conditions indicating that the ArsRS TCS, in addition to direct effects on regulon transcription during acid acclimation, may also indirectly impact gene expression via regulation of the methylome. We determined the target motif of HsdM1 (HP0463) to be the complementary bipartite sequence pair 5'-TCAm6AVN6TGY-3' and 3'-AGTN6GAm6ACA-5'. This complex regulation of DNA methyltransferases, and thus differential methylation patterns, may have implications for the decades-long persistent infection by H. pylori. IMPORTANCE This study expands the possibilities for complex, epigenomic regulation in Helicobacter pylori. We demonstrate that the H. pylori methylome is plastic and acid sensitive via the two-component system ArsRS and the DNA methyltransferase HsdM1. The control of a methyltransferase by ArsRS may allow for a layered response to changing acidity. Likely, an early response whereby ArsR~P affects regulon expression, including the methyltransferase hsdM1. Then, a somewhat later effect as the altered methylome, due to altered HsdM1 expression, subsequently alters the expression of other genes involved in acclimation. The intermediate methylation of certain motifs supports the hypothesis that methyltransferases play a regulatory role. Untangling this additional web of regulation could play a key role in understanding H. pylori colonization and persistence.

Bacterial flagella hijack type IV pili proteins to control motility

Bacterial flagella and type IV pili (TFP) are surface appendages that enable motility and mechanosensing through distinct mechanisms. These structures were previously thought to have no components in common. Here, we report that TFP and some flagella share proteins PilO, PilN, and PilM, which we identified as part of the Helicobacter pylori flagellar motor. H. pylori mutants lacking PilO or PilN migrated better than wild type in semisolid agar because they continued swimming rather than aggregated into microcolonies, mimicking the TFP-regulated surface response. Like their TFP homologs, flagellar PilO/PilN heterodimers formed a peripheral cage that encircled the flagellar motor. These results indicate that PilO and PilN act similarly in flagella and TFP by differentially regulating motility and microcolony formation when bacteria encounter surfaces.

Experimental promoter identification of a foodborne pathogen Salmonella enterica subsp. enterica serovar Typhimurium with near single base-pair resolution

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a common foodborne pathogen which is frequently used as the reference strain for Salmonella. Investigating the sigma factor network and protomers is crucial to understand the genomic and transcriptomic properties of the bacterium. Its promoters were identified using various methods such as dRNA-seq, ChIP-chip, or ChIP-Seq. However, validation using ChIP-exo, which exhibits higher-resolution performance compared to conventional ChIP, has not been conducted to date. In this study, using the representative strain S. Typhimurium LT2 (LT2), the ChIP-exo experiment was conducted to accurately determine the binding sites of catalytic RNA polymerase subunit RpoB and major sigma factors (RpoD, RpoN, RpoS, and RpoE) during exponential phase. Integrated with the results of RNA-Seq, promoters and sigmulons for the sigma factors and their association with RpoB have been discovered. Notably, the overlapping regions among binding sites of each alternative sigma factor were found. Furthermore, comparative analysis with Escherichia coli str. K-12 substr. MG1655 (MG1655) revealed conserved binding sites of RpoD and RpoN across different species. In the case of small RNAs (sRNAs), 50 sRNAs observed their expression during the exponential growth of LT2. Collectively, the integration of ChIP-exo and RNA-Seq enables genome-scale promoter mapping with high resolution and facilitates the characterization of binding events of alternative sigma factors, enabling a comprehensive understanding of the bacterial sigma factor network and condition-specific active promoters.

Precision run-on sequencing (PRO-seq) for microbiome transcriptomics

Bacteria respond to environmental stimuli through precise regulation of transcription initiation and elongation. Bulk RNA sequencing primarily characterizes mature transcripts, so to identify actively transcribed loci we need to capture RNA polymerase (RNAP) complexed with nascent RNA. However, such capture methods have only previously been applied to culturable, genetically tractable organisms such as E. coli and B. subtilis. Here we apply precision run-on sequencing (PRO-seq) to profile nascent transcription in cultured E. coli and diverse uncultured bacteria. We demonstrate that PRO-seq can characterize the transcription of small, structured, or post-transcriptionally modified RNAs, which are often absent from bulk RNA-seq libraries. Applying PRO-seq to the human microbiome highlights taxon-specific RNAP pause motifs and pause-site distributions across non-coding RNA loci that reflect structure-coincident pausing. We also uncover concurrent transcription and cleavage of CRISPR guide RNAs and transfer RNAs. We demonstrate the utility of PRO-seq for exploring transcriptional dynamics in diverse microbial communities.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

Genome-Scale Analysis of the Structure and Function of RNA Pathways and Networks in Pseudomonas aeruginosa

In recent years, several genome-wide approaches based on RNA sequencing (RNA-seq) have been developed. These methods allow a comprehensive and dynamic view of the structure and function of the multi-layered RNA pathways and networks. Many of these approaches, including the promising one of single-cell transcriptome analysis, have been successfully applied to Pseudomonas aeruginosa. However, we are only at the beginning because only a few surrounding conditions have been considered. Here, we aim to illustrate the different types of approaches based on RNA-seq that will lead us in the future to a better understanding of the dynamics of RNA biology in P. aeruginosa.

Acetylation regulates the oligomerization state and activity of RNase J, the Helicobacter pylori major ribonuclease

In the gastric pathogen Helicobacter pylori, post-transcriptional regulation relies strongly on the activity of the essential ribonuclease RNase J. Here, we elucidated the crystal and cryo-EM structures of RNase J and determined that it assembles into dimers and tetramers in vitro. We found that RNase J extracted from H. pylori is acetylated on multiple lysine residues. Alanine substitution of several of these residues impacts on H. pylori morphology, and thus on RNase J function in vivo. Mutations of Lysine 649 modulates RNase J oligomerization in vitro, which in turn influences ribonuclease activity in vitro. Our structural analyses of RNase J reveal loops that gate access to the active site and rationalizes how acetylation state of K649 can influence activity. We propose acetylation as a regulatory level controlling the activity of RNase J and its potential cooperation with other enzymes of RNA metabolism in H. pylori.

Small proteins in bacteria - Big challenges in prediction and identification

Proteins with up to 100 amino acids have been largely overlooked due to the challenges associated with predicting and identifying them using traditional methods. Recent advances in bioinformatics and machine learning, DNA sequencing, RNA and Ribo-seq technologies, and mass spectrometry (MS) have greatly facilitated the detection and characterisation of these elusive proteins in recent years. This has revealed their crucial role in various cellular processes including regulation, signalling and transport, as toxins and as folding helpers for protein complexes. Consequently, the systematic identification and characterisation of these proteins in bacteria have emerged as a prominent field of interest within the microbial research community. This review provides an overview of different strategies for predicting and identifying these proteins on a large scale, leveraging the power of these advanced technologies. Furthermore, the review offers insights into the future developments that may be expected in this field.

TRS: a method for determining transcript termini from RNAtag-seq sequencing data

In bacteria, determination of the 3' termini of transcripts plays an essential role in regulation of gene expression, affecting the functionality and stability of the transcript. Several experimental approaches were developed to identify the 3' termini of transcripts, however, these were applied only to a limited number of bacteria and growth conditions. Here we present a straightforward approach to identify 3' termini from widely available RNA-seq data without the need for additional experiments. Our approach relies on the observation that the RNAtag-seq sequencing protocol results in overabundance of reads mapped to transcript 3' termini. We present TRS (Termini by Read Starts), a computational pipeline exploiting this property to identify 3' termini in RNAtag-seq data, and show that the identified 3' termini are highly reliable. Since RNAtag-seq data are widely available for many bacteria and growth conditions, our approach paves the way for studying bacterial transcription termination in an unprecedented scope.

FlgV forms a flagellar motor ring that is required for optimal motility of Helicobacter pylori

Flagella-driven motility is essential for Helicobacter pylori to colonize the human stomach, where it causes a variety of diseases, including chronic gastritis, peptic ulcer disease, and gastric cancer. H. pylori has evolved a high-torque-generating flagellar motor that possesses several accessories not found in the archetypical Escherichia coli motor. FlgV was one of the first flagellar accessory proteins identified in Campylobacter jejuni, but its structure and function remain poorly understood. Here, we confirm that deletion of flgV in H. pylori B128 and a highly motile variant of H. pylori G27 (G27M) results in reduced motility in soft agar medium. Comparative analyses of in-situ flagellar motor structures of wild-type, ΔflgV, and a strain expressing FlgV-YFP showed that FlgV forms a ring-like structure closely associated with the junction of two highly conserved flagellar components: the MS and C rings. The results of our studies suggest that the FlgV ring has adapted specifically in Campylobacterota to support the assembly and efficient function of the high-torque-generating motors.

Selective RNA Processing and Stabilization are Multi-Layer and Stoichiometric Regulators of Gene Expression in Escherichia coli

Selective RNA processing and stabilization (SRPS) facilitates the differential expression of multiple genes in polycistronic operons. However, how the coordinated actions of SRPS-related enzymes affect stoichiometric regulation remains unclear. In the present study, the first genome-wide targetome analysis is reported of these enzymes in Escherichia coli, at a single-nucleotide resolution. A strictly linear relationship is observed between the RNA pyrophosphohydrolase processing ratio and scores assigned to the first three nucleotides of the primary transcript. Stem-loops associated with PNPase targetomes exhibit a folding free energy that is negatively correlated with the termination ratio of PNPase at the 3' end. More than one-tenth of the RNase E processing sites in the 5'-untranslated regions（UTR） form different stem-loops that affect ribosome-binding and translation efficiency. The effectiveness of the SRPS elements is validated using a dual-fluorescence reporter system. The findings highlight a multi-layer and quantitative regulatory method for optimizing the stoichiometric expression of genes in bacteria and promoting the application of SRPS in synthetic biology.

Profiling of the Helicobacter pylori redox switch HP1021 regulon using a multi-omics approach

The gastric human pathogen Helicobacter pylori has developed mechanisms to combat stress factors, including reactive oxygen species (ROS). Here, we present a comprehensive study on the redox switch protein HP1021 regulon combining transcriptomic, proteomic and DNA-protein interactions analyses. Our results indicate that HP1021 modulates H. pylori's response to oxidative stress. HP1021 controls the transcription of 497 genes, including 407 genes related to response to oxidative stress. 79 proteins are differently expressed in the HP1021 deletion mutant. HP1021 controls typical ROS response pathways (katA, rocF) and less canonical ones, particularly DNA uptake and central carbohydrate metabolism. HP1021 is a molecular regulator of competence in H. pylori, as HP1021-dependent repression of the comB DNA uptake genes is relieved under oxidative conditions, increasing natural competence. Furthermore, HP1021 controls glucose consumption by directly regulating the gluP transporter and has an important impact on maintaining the energetic balance in the cell.

Systems and synthetic biology-driven engineering of live bacterial therapeutics

The past decade has seen growing interest in bacterial engineering for therapeutically relevant applications. While early efforts focused on repurposing genetically tractable model strains, such as Escherichia coli, engineering gut commensals is gaining traction owing to their innate capacity to survive and stably propagate in the intestine for an extended duration. Although limited genetic tractability has been a major roadblock, recent advances in systems and synthetic biology have unlocked our ability to effectively harness native gut commensals for therapeutic and diagnostic purposes, ranging from the rational design of synthetic microbial consortia to the construction of synthetic cells that execute "sense-and-respond" logic operations that allow real-time detection and therapeutic payload delivery in response to specific signals in the intestine. In this review, we outline the current progress and latest updates on microbial therapeutics, with particular emphasis on gut commensal engineering driven by synthetic biology and systems understanding of their molecular phenotypes. Finally, the challenges and prospects of engineering gut commensals for therapeutic applications are discussed.

Current and Future Technologies for the Detection of Antibiotic-Resistant Bacteria

Antibiotic resistance is a global public health concern, posing a significant threat to the effectiveness of antibiotics in treating bacterial infections. The accurate and timely detection of antibiotic-resistant bacteria is crucial for implementing appropriate treatment strategies and preventing the spread of resistant strains. This manuscript provides an overview of the current and emerging technologies used for the detection of antibiotic-resistant bacteria. We discuss traditional culture-based methods, molecular techniques, and innovative approaches, highlighting their advantages, limitations, and potential future applications. By understanding the strengths and limitations of these technologies, researchers and healthcare professionals can make informed decisions in combating antibiotic resistance and improving patient outcomes.

Methodologies for bacterial ribonuclease characterization using RNA-seq

Bacteria adjust gene expression at the post-transcriptional level through an intricate network of small regulatory RNAs and RNA-binding proteins, including ribonucleases (RNases). RNases play an essential role in RNA metabolism, regulating RNA stability, decay, and activation. These enzymes exhibit species-specific effects on gene expression, bacterial physiology, and different strategies of target recognition. Recent advances in high-throughput RNA sequencing (RNA-seq) approaches have provided a better understanding of the roles and modes of action of bacterial RNases. Global studies aiming to identify direct targets of RNases have highlighted the diversity of RNase activity and RNA-based mechanisms of gene expression regulation. Here, we review recent RNA-seq approaches used to study bacterial RNases, with a focus on the methods for identifying direct RNase targets.

Effects of C-Terminal Lys-Arg Residue of AapA1 Protein on Toxicity and Structural Mechanism

Previous experimental investigations have established the indispensability of the C-terminal Lys-Arg residues in the toxic activity of the AapA1 toxin protein. AapA1 is classified as a type I toxin-antitoxin (TA) bacterial toxin, and the precise impact of the C-terminal Lys-Arg residues on its structure and mechanism of action remains elusive. To address this knowledge gap, the present study employed molecular dynamics (MD) and enhanced sampling Well-tempered Two-dimensional Metadynamics (2D-MetaD) simulations to examine the behavior of the C-terminal Lys-Arg residues of truncated AapA1 toxin (AapA1-28) within the inner membrane of Escherichia coli. Specifically, the study focused on the elucidation of possible conformation states of AapA1-28 protein in POPE/POPG (3:1) bilayers and their interactions between the protein and POPE/POPG (3:1) bilayers. The findings of our investigation indicate that the AapA1-28 protein does not adopt a vertical orientation upon membrane insertion; rather, it assumes an angled conformation, with the side chain of Lys-23 directed toward the upper layer of the membrane. This non-transmembrane conformation of AapA1-28 protein impedes its ability to form pores within the membrane, resulting in reduced toxicity towards Escherichia coli. These results suggest that C-Terminal positively charged residues are essential for electrostatic binding to the negatively charged head group of bottom bilayer membrane, which stabilize the transmembrane conformation. These outcomes contribute to our comprehension of the impact of C-terminal charged residues on the structure and functionality of membrane-associated proteins, and provide an improved understanding of how protein sequence influences the antimicrobial effect.

Non-canonical transcriptional start sites in E. coli O157:H7 EDL933 are regulated and appear in surprisingly high numbers

Analysis of genome wide transcription start sites (TSSs) revealed an unexpected complexity since not only canonical TSS of annotated genes are recognized by RNA polymerase. Non-canonical TSS were detected antisense to, or within, annotated genes as well new intergenic (orphan) TSS, not associated with known genes. Previously, it was hypothesized that many such signals represent noise or pervasive transcription, not associated with a biological function. Here, a modified Cappable-seq protocol allows determining the primary transcriptome of the enterohemorrhagic E. coli O157:H7 EDL933 (EHEC). We used four different growth media, both in exponential and stationary growth phase, replicated each thrice. This yielded 19,975 EHEC canonical and non-canonical TSS, which reproducibly occurring in three biological replicates. This questions the hypothesis of experimental noise or pervasive transcription. Accordingly, conserved promoter motifs were found upstream indicating proper TSSs. More than 50% of 5,567 canonical and between 32% and 47% of 10,355 non-canonical TSS were differentially expressed in different media and growth phases, providing evidence for a potential biological function also of non-canonical TSS. Thus, reproducible and environmentally regulated expression suggests that a substantial number of the non-canonical TSSs may be of unknown function rather than being the result of noise or pervasive transcription.

Transcriptome architecture of the three main lineages of agrobacteria

Agrobacteria are a diverse, polyphyletic group of prokaryotes with multipartite genomes capable of transferring DNA into the genomes of host plants, making them an essential tool in plant biotechnology. Despite their utility in plant transformation, genome-wide transcriptional regulation is not well understood across the three main lineages of agrobacteria. Transcription start sites (TSSs) are a necessary component of gene expression and regulation. In this study, we used differential RNA-seq and a TSS identification algorithm optimized on manually annotated TSS, then validated with existing TSS to identify thousands of TSS with nucleotide resolution for representatives of each lineage. We extend upon the 356 TSSs previously reported in Agrobacterium fabrum C58 by identifying 1,916 TSSs. In addition, we completed genomes and phenotyping of Rhizobium rhizogenes C16/80 and Allorhizobium vitis T60/94, identifying 2,650 and 2,432 TSSs, respectively. Parameter optimization was crucial for an accurate, high-resolution view of genome and transcriptional dynamics, highlighting the importance of algorithm optimization in genome-wide TSS identification and genomics at large. The optimized algorithm reduced the number of TSSs identified internal and antisense to the coding sequence on average by 90.5% and 91.9%, respectively. Comparison of TSS conservation between orthologs of the three lineages revealed differences in cell cycle regulation of ctrA as well as divergence of transcriptional regulation of chemotaxis-related genes when grown in conditions that simulate the plant environment. These results provide a framework to elucidate the mechanistic basis and evolution of pathology across the three main lineages of agrobacteria. IMPORTANCE Transcription start sites (TSSs) are fundamental for understanding gene expression and regulation. Agrobacteria, a group of prokaryotes with the ability to transfer DNA into the genomes of host plants, are widely used in plant biotechnology. However, the genome-wide transcriptional regulation of agrobacteria is not well understood, especially in less-studied lineages. Differential RNA-seq and an optimized algorithm enabled identification of thousands of TSSs with nucleotide resolution for representatives of each lineage. The results of this study provide a framework for elucidating the mechanistic basis and evolution of pathology across the three main lineages of agrobacteria. The optimized algorithm also highlights the importance of parameter optimization in genome-wide TSS identification and genomics at large.

Global Trend on Machine Learning in Helicobacter within One Decade: A Scientometric Study

Purpose

This study aims to create a science map, provide structural analysis, investigate evolution, and identify new trends in Helicobacter pylori (H. pylori) research articles.

Methods

All Helicobacter publications were gathered from the Web of Science (WoS) database from August 2010 to 2021. The data were required for bibliometric analysis. The bibliometric analysis was performed with Bibliometrix R Tool. Bibliometric data were analyzed using the Bibliometrix Biblioshiny R-package software.

Results

A total of 17,413 articles were reviewed and analyzed, with descriptive characteristics of the H. pylori literature included. In journals, 21,102 keywords plus and 20,490 author keywords were reported. These articles were also written by 56,106 different authors, with 262 being single-author articles. Most authors' abstracts, titles, and keywords included "Helicobacter-pylori." Since 2010, the total number of H. pylori-related publications has been decreasing. Gut, PLOS ONE, and Gastroenterology are the most influential H. pylori journals, according to source impact. China, the United States, and Japan are the countries with most affiliations and subjects. In addition, Seoul National University has published the most articles about H. pylori. According to the cloud word plot, the authors' most frequently used keywords are gastric cancer (GC), H. pylori, gastritis, eradication, and inflammation. "Helicobacter pylori" and "infection" have the steepest slopes in terms of the upward trend of words used in articles from 2010 to 2021. Subjects such as GC, intestinal metaplasia, epidemiology, peptic ulcer, eradication, and clarithromycin are included in the diagram's motor theme section, according to strategic diagrams. According to the thematic evolution map, topics such as Helicobacter pylori infection, B-cell lymphoma, CagA, Helicobacter pylori, and infection were largely discussed between 2010 and 2015. From 2016 to 2021, the top topics covered included Helicobacter pylori, H. pylori infection, and infection.

Improved RNA stability estimation indicates that transcriptional interference is frequent in diverse bacteria

We used stochastic simulations and experimental data from E. coli, K. aerogenes, Synechococcus PCC 7002 and Synechocystis PCC 6803 to provide evidence that transcriptional interference via the collision mechanism is likely a prevalent mechanism for bacterial gene regulation. Rifampicin time-series data can be used to globally monitor and quantify collision between sense and antisense transcription-complexes. Our findings also highlight that transcriptional events, such as differential RNA decay, partial termination, and internal transcriptional start sites often deviate from gene annotations. Consequently, within a single gene annotation, there exist transcript segments with varying half-lives and transcriptional properties. To address these complexities, we introduce 'rifi', an R-package that analyzes transcriptomic data from rifampicin time series. 'rifi' employs a dynamic programming-based segmentation approach to identify individual transcripts, enabling accurate assessment of RNA stability and detection of diverse transcriptional events.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen Borrelia burgdorferi

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. Here, we use several RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. We identify complex gene arrangements and operons, untranslated regions and small RNAs. We predict intrinsic terminators and experimentally test examples of Rho-dependent transcription termination. Remarkably, 63% of RNA 3' ends map upstream of or internal to open reading frames (ORFs), including genes involved in the unique infectious cycle of B. burgdorferi. We suggest these RNAs result from premature termination, processing and regulatory events such as cis-acting regulation. Furthermore, the polyamine spermidine globally influences the generation of truncated mRNAs. Collectively, our findings provide insights into transcription termination and uncover an abundance of potential RNA regulators in B. burgdorferi.

Divide and conquer: genetics, mechanism, and evolution of the ferrous iron transporter Feo in Helicobacter pylori

Introduction

Feo is the most widespread and conserved system for ferrous iron uptake in bacteria, and it is important for virulence in several gastrointestinal pathogens. However, its mechanism remains poorly understood. Hitherto, most studies regarding the Feo system were focused on Gammaproteobacterial models, which possess three feo genes (feoA, B, and C) clustered in an operon. We found that the human pathogen Helicobacter pylori possesses a unique arrangement of the feo genes, in which only feoA and feoB are present and encoded in distant loci. In this study, we examined the functional significance of this arrangement.

Methods

Requirement and regulation of the individual H. pylori feo genes were assessed through in vivo assays and gene expression profiling. The evolutionary history of feo was inferred via phylogenetic reconstruction, and AlphaFold was used for predicting the FeoA-FeoB interaction.

Results and Discussion

Both feoA and feoB are required for Feo function, and feoB is likely subjected to tight regulation in response to iron and nickel by Fur and NikR, respectively. Also, we established that feoA is encoded in an operon that emerged in the common ancestor of most, but not all, helicobacters, and this resulted in feoA transcription being controlled by two independent promoters. The H. pylori Feo system offers a new model to understand ferrous iron transport in bacterial pathogens.

Genetic diversity of Helicobacter pylori type IV secretion system cagI and cagN genes and their association with clinical diseases

A number of cagPAI genes in the Helicobacter pylori genome are considered the most evolved genes under a diversifying selection and evolutionary pressure. Among them, cagI and cagN are described as a part of the two different-operon of cagPAI that are involved in the T4SS machinery, but the definite association of these factors with clinical manifestations is still unclear. A total of 70 H. pylori isolates were obtained from different gastroduodenal patients. All isolates were examined for the presence of primary H. pylori virulence genes by PCR analysis. Direct DNA sequence analysis was performed for the cagI and cagN genes. The results were compared with the reference strain. The cagI, cagN, cagA, cagL, vacA s1m1, vacA s1m2, vacA s2m2, babA2, sabA, and dupA genotypes were detected in 80, 91.4, 84, 91.4, 32.8, 42.8, 24.4, 97.1, 84.3, and 84.3% of the total isolates, respectively. The most variable codon usage in cagI was observed at residues 20-25, 55-60, 94, 181-199, 213-221, 241-268, and 319-320, while the most variable codon usage in CagN hypervariable motif (CagNHM) was observed at residues 53 to 63. Sequencing data analysis of cagN revealed a hypothetical hexapeptide motif (EAKDEN/K) in residues of 278-283 among six H. pylori isolates, which needs further studies to evaluate its putative function. The present study demonstrated a high prevalence of cagI and cagN genes among Iranian H. pylori isolates with gastroduodenal diseases. Furthermore, no significant correlation between cagI and cagN variants and clinical diseases was observed in the present study. However, all patients had a high prevalence of cagPAI genes including cagI, cagN, cagA, and cagL, which indicates more potential role of these genes in disease outcome.

Helicobacter pylori Modulates Heptose Metabolite Biosynthesis and Heptose-Dependent Innate Immune Host Cell Activation by Multiple Mechanisms

Heptose metabolites including ADP-d-glycero-β-d-manno-heptose (ADP-heptose) are involved in bacterial lipopolysaccharide and cell envelope biosynthesis. Recently, heptoses were also identified to have potent proinflammatory activity on human cells as novel microbe-associated molecular patterns. The gastric pathogenic bacterium Helicobacter pylori produces heptose metabolites, which it transports into human cells through its Cag type 4 secretion system. Using H. pylori as a model, we have addressed the question of how proinflammatory ADP-heptose biosynthesis can be regulated by bacteria. We have characterized the interstrain variability and regulation of heptose biosynthesis genes and the modulation of heptose metabolite production by H. pylori, which impact cell-autonomous proinflammatory human cell activation. HldE, a central enzyme of heptose metabolite biosynthesis, showed strong sequence variability between strains and was also variably expressed between strains. Amounts of gene transcripts in the hldE gene cluster displayed intrastrain and interstrain differences, were modulated by host cell contact and the presence of the cag pathogenicity island, and were affected by carbon starvation regulator A (CsrA). We reconstituted four steps of the H. pylori lipopolysaccharide (LPS) heptose biosynthetic pathway in vitro using recombinant purified GmhA, HldE, and GmhB proteins. On the basis of one- and two-dimensional nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry, the structures of major reaction products were identified as β-d-ADP-heptose and β-heptose-1-monophosphate. A proinflammatory heptose-monophosphate variant was also identified for the first time as a novel cell-active product in H. pylori bacteria. Separate purified HldE subdomains and variant HldE allowed us to uncover additional strain variation in generating heptose metabolites. IMPORTANCE Bacterial heptose metabolites, intermediates of lipopolysaccharide (LPS) biosynthesis, are novel microbe-associated molecular patterns (MAMPs) that activate proinflammatory signaling. In the gastric pathogen Helicobacter pylori, heptoses are transferred into host cells by the Cag type IV secretion system, which is also involved in carcinogenesis. Little is known about how H. pylori, which is highly strain variable, regulates heptose biosynthesis and downstream host cell activation. We report here that the regulation of proinflammatory heptose production by H. pylori is strain specific. Heptose gene cluster activity is modulated by the presence of an active cag pathogenicity island (cagPAI), contact with human cells, and the carbon starvation regulator A. Reconstitution with purified biosynthesis enzymes and purified bacterial lysates allowed us to biochemically characterize heptose pathway products, identifying a heptose-monophosphate variant as a novel proinflammatory metabolite. These findings emphasize that the bacteria use heptose biosynthesis to fine-tune inflammation and also highlight opportunities to mine the heptose biosynthesis pathway as a potential therapeutic target against infection, inflammation, and cancer.

Nitrogen-regulated antisense transcription in the adaptation to nitrogen deficiency in Nostoc sp. PCC 7120

Transcriptomic analyses using high-throughput methods have revealed abundant antisense transcription in bacteria. Antisense transcription is often due to the overlap of mRNAs with long 5' or 3' regions that extend beyond the coding sequence. In addition, antisense RNAs that do not contain any coding sequence are also observed. Nostoc sp. PCC 7120 is a filamentous cyanobacterium that, under nitrogen limitation, behaves as a multicellular organism with division of labor among two different cell types that depend on each other, the vegetative CO₂-fixing cells and the nitrogen-fixing heterocysts. The differentiation of heterocysts depends on the global nitrogen regulator NtcA and requires the specific regulator HetR. To identify antisense RNAs potentially involved in heterocyst differentiation, we assembled the Nostoc transcriptome using RNA-seq analysis of cells subjected to nitrogen limitation (9 or 24 h after nitrogen removal) in combination with a genome-wide set of transcriptional start sites and a prediction of transcriptional terminators. Our analysis resulted in the definition of a transcriptional map that includes >4,000 transcripts, 65% of which contain regions in antisense orientation to other transcripts. In addition to overlapping mRNAs, we identified nitrogen-regulated noncoding antisense RNAs transcribed from NtcA- or HetR-dependent promoters. As an example of this last category, we further analyzed an antisense (as_gltA) of the gene-encoding citrate synthase and showed that transcription of as_gltA takes place specifically in heterocysts. Since the overexpression of as_gltA reduces citrate synthase activity, this antisense RNA could eventually contribute to the metabolic remodeling that occurs during the differentiation of vegetative cells into heterocysts.

Characterization of Bacterial Transcriptional Regulatory Networks in Escherichia coli through Genome-Wide In Vitro Run-Off Transcription/RNA-seq (ROSE)

The global characterization of transcriptional regulatory networks almost exclusively uses in vivo conditions, thereby providing a snapshot on multiple regulatory interactions at the same time. To complement these approaches, we developed and applied a method for characterizing bacterial promoters genome-wide by in vitro transcription coupled to transcriptome sequencing specific for native 5'-ends of transcripts. This method, called ROSE (run-off transcription/RNA-sequencing), only requires chromosomal DNA, ribonucleotides, RNA polymerase (RNAP) core enzyme, and a specific sigma factor, recognizing the corresponding promoters, which have to be analyzed. ROSE was performed on E. coli K-12 MG1655 genomic DNA using Escherichia coli RNAP holoenzyme (including σ70) and yielded 3226 transcription start sites, 2167 of which were also identified in in vivo studies, and 598 were new. Many new promoters not yet identified by in vivo experiments might be repressed under the tested conditions. Complementary in vivo experiments with E. coli K-12 strain BW25113 and isogenic transcription factor gene knockout mutants of fis, fur, and hns were used to test this hypothesis. Comparative transcriptome analysis demonstrated that ROSE could identify bona fide promoters that were apparently repressed in vivo. In this sense, ROSE is well-suited as a bottom-up approach for characterizing transcriptional networks in bacteria and ideally complementary to top-down in vivo transcriptome studies.

Discovery of Type IV filament membrane alignment complex homologs in H. pylori that promote soft-agar migration

The stomach pathogen Helicobacter pylori utilizes two scaffold proteins, CheW and CheV1, to build critical chemotaxis arrays. Chemotaxis helps bacteria establish and maintain infection. Mutants lacking either of these chemotaxis proteins have different soft agar phenotypes: deletion of cheW creates non-chemotactic strains, while deletion of cheV1 results in 50% loss of chemotaxis. In this work, we characterized the cheV1 deletion mutant phenotype in detail. cheV1 deletion mutants had poor soft-agar migration initially, but regained migration ability over time. This improved bacterial migration was stable, suggesting a genetic suppressor phenotype, termed Che+. Whole-genome sequencing analysis of four distinct cheV1 Che+ strains revealed single nucleotide polymorphisms (SNPs) in a common gene, HPG27_252 (HP0273). These SNPs were predicted to truncate the encoded protein. To confirm the role of HPG27_252 in the cheV1 phenotype, we created a targeted deletion of HPG27_252 and found that loss of HPG27_252 enhanced soft-agar migration. HPG27_252 and CheV1 appear to interact directly, based on bacterial two-hybrid analysis. HPG27_252 is predicted to encode a 179 amino acid, 21 kDa protein annotated as a hypothetical protein. Computational analysis revealed this protein to be a remote homolog of the PilO Type IV filament membrane alignment complex protein. Although H. pylori is not known to possess Type IV filaments, our analysis showed it retains an operon of genes for homologs of PilO, PilN, and PilM, but does not possess other Type IV pili genes. Our data suggest the PilO homolog plays a role in regulating H. pylori chemotaxis and motility, suggesting new ideas about evolutionary steps for controlling migration through semi-solid media.

Assessing the Orthogonality of Phage-Encoded RNA Polymerases for Tailored Synthetic Biology Applications in Pseudomonas Species

The phage T7 RNA polymerase (RNAP) and lysozyme form the basis of the widely used pET expression system for recombinant expression in the biotechnology field and as a tool in microbial synthetic biology. Attempts to transfer this genetic circuitry from Escherichia coli to non-model bacterial organisms with high potential have been restricted by the cytotoxicity of the T7 RNAP in the receiving hosts. We here explore the diversity of T7-like RNAPs mined directly from Pseudomonas phages for implementation in Pseudomonas species, thus relying on the co-evolution and natural adaptation of the system towards its host. By screening and characterizing different viral transcription machinery using a vector-based system in P. putida., we identified a set of four non-toxic phage RNAPs from phages phi15, PPPL-1, Pf-10, and 67PfluR64PP, showing a broad activity range and orthogonality to each other and the T7 RNAP. In addition, we confirmed the transcription start sites of their predicted promoters and improved the stringency of the phage RNAP expression systems by introducing and optimizing phage lysozymes for RNAP inhibition. This set of viral RNAPs expands the adaption of T7-inspired circuitry towards Pseudomonas species and highlights the potential of mining tailored genetic parts and tools from phages for their non-model host.

Antisense transcription and its roles in adaption to environmental stress in E. coli

It has been reported that a highly varying proportion (1% ∼ 93%) of genes in various prokaryotes have antisense RNA (asRNA) transcription. However, the extent of the pervasiveness of asRNA transcription in the well-studied E. coli K12 strain has thus far been an issue of debate. Furthermore, very little is known about the expression patterns and functions of asRNAs under various conditions. To fill these gaps, we determined the transcriptomes and proteomes of E. coli K12 at multiple time points in five culture conditions using strand-specific RNA-seq, differential RNA-seq, and quantitative mass spectrometry methods. To reduce artifacts of possible transcriptional noise, we identified asRNA using stringent criteria with biological replicate verification and transcription start sites (TSSs) information included. We identified a total of 660 asRNAs, which were generally short and largely condition-dependently transcribed. We found that the proportions of the genes which had asRNA transcription highly depended on the culture conditions and time points. We classified the transcriptional activities of the genes in six transcriptional modes according to their relative levels of asRNA to mRNA. Many genes changed their transcriptional modes at different time points of the culture conditions, and such transitions can be described in a well-defined manner. Intriguingly, the protein levels and mRNA levels of genes in the sense-only/sense-dominant mode were moderately correlated, but the same was not true for genes in the balanced/antisense-dominant mode, in which asRNAs were at a comparable or higher level to mRNAs. These observations were further validated by western blot on candidate genes, where an increase in asRNA transcription diminished gene expression in one case and enhanced it in another. These results suggest that asRNAs may directly or indirectly regulate translation by forming duplexes with cognate mRNAs. Thus, asRNAs may play an important role in the bacterium's responses to environmental changes during growth and adaption to different environments.

IMPORTANCE

The cis -antisense RNA (asRNA) is a type of understudied RNA molecules in prokaryotes, which is believed to be important in regulating gene expression. Our current understanding of asRNA is constrained by inconsistent reports about its identification and properties. These discrepancies are partially caused by a lack of sufficient samples, biological replicates, and culture conditions. This study aimed to overcome these disadvantages and identified 660 putative asRNAs using integrated information from strand-specific RNA-seq, differential RNA-seq, and mass spectrometry methods. In addition, we explored the relative expression between asRNAs and sense RNAs and investigated asRNA regulated transcriptional activity changes over different culture conditions and time points. Our work strongly suggests that asRNAs may play a crucial role in bacterium's responses to environmental changes during growth and adaption to different environments.

Molecular dynamics simulation-based trinucleotide and tetranucleotide level structural and energy characterization of the functional units of genomic DNA

Genomes of most organisms on earth are written in a universal language of life, made up of four units - adenine (A), thymine (T), guanine (G), and cytosine (C), and understanding the way they are put together has been a great challenge to date. Multiple efforts have been made to annotate this wonderfully engineered string of DNA using different methods but they lack a universal character. In this article, we have investigated the structural and energetic profiles of both prokaryotes and eukaryotes by considering two essential genomic sites, viz., the transcription start sites (TSS) and exon-intron boundaries. We have characterized these sites by mapping the structural and energy features of DNA obtained from molecular dynamics simulations, which considers all possible trinucleotide and tetranucleotide steps. For DNA, these physicochemical properties show distinct signatures at the TSS and intron-exon boundaries. Our results firmly convey the idea that DNA uses the same dialect for prokaryotes and eukaryotes and that it is worth going beyond sequence-level analyses to physicochemical space to determine the functional destiny of DNA sequences.

Transcription start site mapping and small RNA profiling of Leptospira biflexa serovar Patoc

Leptospirosis is an emerging zoonotic disease caused by bacterial species of the genus Leptospira. However, the regulatory mechanisms and pathways underlying the adaptation of pathogenic and non-pathogenic Leptospira spp. in different environmental conditions remain elusive. Leptospira biflexa is a non-pathogenic species of Leptospira that lives exclusively in a natural environment. It is an ideal model not only for exploring molecular mechanisms underlying the environmental survival of Leptospira species but also for identifying virulence factors unique to Leptospira's pathogenic species. In this study, we aim to establish the transcription start site (TSS) landscape and the small RNA (sRNA) profile of L. biflexa serovar Patoc grown to exponential and stationary phases via differential RNA-seq (dRNA-seq) and small RNA-seq (sRNA-seq) analyses, respectively. Our dRNA-seq analysis uncovered a total of 2726 TSSs, which are also used to identify other elements, e.g., promoter and untranslated regions (UTRs). Besides, our sRNA-seq analysis revealed a total of 603 sRNA candidates, comprising 16 promoter-associated sRNAs, 184 5'UTR-derived sRNAs, 230 true intergenic sRNAs, 136 5'UTR-antisense sRNAs, and 130 open reading frame (ORF)-antisense sRNAs. In summary, these findings reflect the transcriptional complexity of L. biflexa serovar Patoc under different growth conditions and help to facilitate our understanding of regulatory networks in L. biflexa. To the best of our knowledge, this is the first study reporting the TSS landscape of L. biflexa. The TSS and sRNA landscapes of L. biflexa can also be compared with its pathogenic counterparts, e.g., L. borgpetersenii and L. interrogans, to identify features contributing to their environmental survival and virulence.

Exopolysaccharide Biosynthesis in Rhizobium leguminosarum bv. trifolii Requires a Complementary Function of Two Homologous Glycosyltransferases PssG and PssI

The Pss-I region of Rhizobium leguminosarum bv. trifolii TA1 comprises more than 20 genes coding for glycosyltransferases, modifying enzymes, and polymerization/export proteins, altogether determining the biosynthesis of symbiotically relevant exopolysaccharides. In this study, the role of homologous PssG and PssI glycosyltransferases in exopolysaccharide subunit synthesis were analyzed. It was shown that the glycosyltransferase-encoding genes of the Pss-I region were part of a single large transcriptional unit with potential downstream promoters activated in specific conditions. The ΔpssG and ΔpssI mutants produced significantly lower amounts of the exopolysaccharide, while the double deletion mutant ΔpssIΔpssG produced no exopolysaccharide. Complementation of double mutation with individual genes restored exopolysaccharide synthesis, but only to the level similar to that observed for the single ΔpssI or ΔpssG mutants, indicating that PssG and PssI serve complementary functions in the process. PssG and PssI interacted with each other in vivo and in vitro. Moreover, PssI displayed an expanded in vivo interaction network comprising other GTs involved in subunit assembly and polymerization/export proteins. PssG and PssI proteins were shown to interact with the inner membrane through amphipathic helices at their C-termini, and PssG also required other proteins involved in exopolysaccharide synthesis to localize in the membrane protein fraction.

Regulation of Helicobacter pylori Urease and Acetone Carboxylase Genes by Nitric Oxide and the CrdRS Two-Component System

Helicobacter pylori colonizes the human gastric mucosa and causes various gastroduodenal diseases, including peptic ulceration and gastric cancer. Colonization requires the actions of two-component systems (TCSs) to sense and respond to changes in the host environment. In this study, we evaluated gene regulation mediated by the CrdRS TCS. Few studies have evaluated this TCS, leaving the signal(s) yet to be exhaustively determined and a need for a more complete regulon to be delineated. We performed RNA sequencing (RNA-Seq) on three isogenic H. pylori 26695 mutants: a control, a mutant with deletion of the sensory histidine kinase, ΔcrdS, and a mutant with deletion of the response regulator, ΔcrdR. Comparison of the RNA-Seq results from these mutants established a 40-gene regulon putatively controlled by the CrdRS TCS. Quantitative reverse transcriptase PCR (RT-qPCR) was used to validate 7 of 11 putative regulon members selected for analysis. We further investigated 6 confirmed CrdRS regulon genes by using phospho-incompetent H. pylori 26695 CrdR D53A and CrdS H173A mutants. These experiments further confirmed the role of CrdRS in regulation of urease, acetone carboxylase, hofD, and HP1440. Expression of these CrdRS regulon genes was also evaluated under 10 μM nitric oxide (NO) conditions. This revealed that ureA, acxA, hofD, and HP1440 expression is affected by NO in a CrdRS-dependent manner. Importantly, three of these genes (ureA, acxA, and hofD) are known to play important roles in H. pylori colonization of the stomach. IMPORTANCE The molecular strategies used by Helicobacter pylori to colonize and persist in the harsh environment of the human stomach are a critical area of study. Our study identified several genes in this gastric pathogen, including ureA, a gene encoding a protein essential to the survival of H. pylori, that are regulated via the CrdRS two-component system (TCS) in response to nitric oxide (NO). NO is a product of the innate immune system of the human host. The identification of these genes whose expression is regulated by this molecule may give insights to novel therapeutics. Two genes (ureA and acxA) determined in this study to be regulated by NO via CrdRS have been previously determined to be regulated by other TCSs, indicating that the expression of these genes may be of critical importance to H. pylori.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen B. burgdorferi

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. We employed complementary RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. By systematically mapping B. burgdorferi RNA ends at single nucleotide resolution, we delineated complex gene arrangements and operons and mapped untranslated regions (UTRs) and small RNAs (sRNAs). We experimentally tested modes of B. burgdorferi transcription termination and compared our findings to observations in E. coli , P. aeruginosa , and B. subtilis . We discovered 63% of B. burgdorferi RNA 3' ends map upstream or internal to open reading frames (ORFs), suggesting novel mechanisms of regulation. Northern analysis confirmed the presence of stable 5' derived RNAs from mRNAs encoding gene products involved in the unique infectious cycle of B. burgdorferi . We suggest these RNAs resulted from premature termination and regulatory events, including forms of cis- acting regulation. For example, we documented that the polyamine spermidine globally influences the generation of truncated mRNAs. In one case, we showed that high spermidine concentrations increased levels of RNA fragments derived from an mRNA encoding a spermidine import system, with a concomitant decrease in levels of the full- length mRNA. Collectively, our findings revealed new insight into transcription termination and uncovered an abundance of potential RNA regulators.

RNA Preparation and RNA-Seq Bioinformatics for Comparative Transcriptomics

The principal transcriptome analysis is the determination of differentially expressed genes across experimental conditions. For this, the next-generation sequencing of RNA (RNA-seq) has several advantages over other techniques, such as the capability of detecting all the transcripts in one assay over RT-qPCR, such as its higher accuracy and broader dynamic range over microarrays or the ability to detect novel transcripts, including non-coding RNA molecules, at nucleotide-level resolution over both techniques. Despite these advantages, many microbiology laboratories have not yet applied RNA-seq analyses to their investigations. The high cost of the equipment for next-generation sequencing is no longer an issue since this intermediate part of the analysis can be provided by commercial or central services. Here, we detail a protocol for the first part of the analysis, the RNA extraction and an introductory protocol to the bioinformatics analysis of the sequencing data to generate the differential expression results.