The primary transcriptome of the major human pathogen Helicobacter pylori

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.

Recent Advances in Genomics-Based Approaches for the Development of Intracellular Bacterial Pathogen Vaccines

Infectious diseases continue to be a leading cause of morbidity and mortality worldwide. The majority of infectious diseases are caused by intracellular pathogenic bacteria (IPB). Historically, conventional vaccination drives have helped control the pathogenesis of intracellular bacteria and the emergence of antimicrobial resistance, saving millions of lives. However, in light of various limitations, many diseases that involve IPB still do not have adequate vaccines. In response to increasing demand for novel vaccine development strategies, a new area of vaccine research emerged following the advent of genomics technology, which changed the paradigm of vaccine development by utilizing the complete genomic data of microorganisms against them. It became possible to identify genes related to disease virulence, genetic patterns linked to disease virulence, as well as the genetic components that supported immunity and favorable vaccine responses. Complete genomic databases, and advancements in transcriptomics, metabolomics, structural genomics, proteomics, immunomics, pan-genomics, synthetic genomics, and population biology have allowed researchers to identify potential vaccine candidates and predict their effects in patients. New vaccines have been created against diseases for which previously there were no vaccines available, and existing vaccines have been improved. This review highlights the key issues and explores the evolution of vaccines. The increasing volume of IPB genomic data, and their application in novel genome-based techniques for vaccine development, were also examined, along with their characteristics, and the opportunities and obstacles involved. Critically, the application of genomics technology has helped researchers rapidly select and evaluate candidate antigens. Novel vaccines capable of addressing the limitations associated with conventional vaccines have been developed and pressing healthcare issues are being addressed.

Dickkopf-related protein 1/cytoskeleton-associated protein 4 signaling activation by Helicobacter pylori-induced activator protein-1 promotes gastric tumorigenesis via the PI3K/AKT/mTOR pathway

BACKGROUND

Gastric cancer (GC) is a common malignant tumor with high incidence and mortality rates globally, especially in East Asian countries. Helicobacter pylori (H. pylori) infection is a significant and independent risk factor for GC. However, its underlying mechanism of action is not fully understood. Dickkopf-related protein (DKK) 1 is a Wnt signaling antagonist, and cytoskeleton-associated protein (CKAP) 4 is a newly identified DKK1 receptor. Recent studies found that the binding of DKK1 to CAKP4 mediated the procancer signaling of DKK1 inde-pendent of Wnt signaling. We hypothesize that H. pylori-induced activation of DKK1/CKAP4 signaling contributes to the initiation and progression of GC.

AIM

To investigate the interaction of H. pylori infection, DKK1 and CAKP4 in GC, as well as the underlying molecular mechanisms.

METHODS

RNA sequencing was used to identify differentially expressed genes (DEGs) between H. pylori-infected and uninfected primary GC cells. Gain- and loss-of-function experiments were performed to verify the H. pylori-induced upregulation of activator protein-1 (AP-1) in GC cells. A dual-luciferase reporter assay and co-immunoprecipitation were used to determine the binding of AP-1 to the DKK1 promoter and DKK1 to CKAP4. Western blotting and immunohistochemistry detected the expression of DKK1, CKAP4, and phos-phatidylinositol 3-kinase (PI3K) pathway-related proteins in GC cells and tissues. Functional experiments and tumorigenicity in nude mice detected malignant behavior of GC cells in vitro and in vivo.

RESULTS

We identified 32 DEGs between primary GC cells with and without H. pylori infection, including JUN, fos-like antigen-1 (FOSL1), and DKK1, and confirmed that the three proteins and CKAP4 were highly expressed in H. pylori-infected GC cells, H. pylori-infected gerbil gastric tissues, and human GC tissues. JUN and FOSL1 form AP-1 to transcriptionally activate DKK1 expression by binding to the DKK1 promoter. Activated DKK1 bound to CKAP4, but not the most common Wnt coreceptor low-density lipoprotein receptor-related protein 5/6, to promote GC cell growth, colony formation, migration, invasion, and xenograft tumor growth in nude mice. All these effects were driven by activation of the PI3K/AKT/mammalian target of rapamycin (mTOR) pathway. Targeting the PI3K signaling pathway by LY294002 inhibited DKK1-mediated CKAP4/PI3K signaling activity and the malignant behavior of GC cells.

CONCLUSION

H. pylori induces JUN and FOSL1 expression to form AP-1, which transcriptionally activates DKK1. Binding of DKK1 to KAKP4 contributes to gastric tumorigenesis via the PI3K/AKT/mTOR pathway.

High-Resolution Small RNAs Landscape Provides Insights into Alkane Adaptation in the Marine Alkane-Degrader Alcanivorax dieselolei B-5

Alkanes are widespread in the ocean, and Alcanivorax is one of the most ubiquitous alkane-degrading bacteria in the marine ecosystem. Small RNAs (sRNAs) are usually at the heart of regulatory pathways, but sRNA-mediated alkane metabolic adaptability still remains largely unknown due to the difficulties of identification. Here, differential RNA sequencing (dRNA-seq) modified with a size selection (~50-nt to 500-nt) strategy was used to generate high-resolution sRNAs profiling in the model species Alcanivorax dieselolei B-5 under alkane (n-hexadecane) and non-alkane (acetate) conditions. As a result, we identified 549 sRNA candidates at single-nucleotide resolution of 5'-ends, 63.4% of which are with transcription start sites (TSSs), and 36.6% of which are with processing sites (PSSs) at the 5'-ends. These sRNAs originate from almost any location in the genome, regardless of intragenic (65.8%), antisense (20.6%) and intergenic (6.2%) regions, and RNase E may function in the maturation of sRNAs. Most sRNAs locally distribute across the 15 reference genomes of Alcanivorax, and only 7.5% of sRNAs are broadly conserved in this genus. Expression responses to the alkane of several core conserved sRNAs, including 6S RNA, M1 RNA and tmRNA, indicate that they may participate in alkane metabolisms and result in more actively global transcription, RNA processing and stresses mitigation. Two novel CsrA-related sRNAs are identified, which may be involved in the translational activation of alkane metabolism-related genes by sequestering the global repressor CsrA. The relationships of sRNAs with the characterized genes of alkane sensing (ompS), chemotaxis (mcp, cheR, cheW2), transporting (ompT1, ompT2, ompT3) and hydroxylation (alkB1, alkB2, almA) were created based on the genome-wide predicted sRNA-mRNA interactions. Overall, the sRNA landscape lays the ground for uncovering cryptic regulations in critical marine bacterium, among which both the core and species-specific sRNAs are implicated in the alkane adaptive metabolisms.

Developing New Tools to Fight Human Pathogens: A Journey through the Advances in RNA Technologies

A long scientific journey has led to prominent technological advances in the RNA field, and several new types of molecules have been discovered, from non-coding RNAs (ncRNAs) to riboswitches, small interfering RNAs (siRNAs) and CRISPR systems. Such findings, together with the recognition of the advantages of RNA in terms of its functional performance, have attracted the attention of synthetic biologists to create potent RNA-based tools for biotechnological and medical applications. In this review, we have gathered the knowledge on the connection between RNA metabolism and pathogenesis in Gram-positive and Gram-negative bacteria. We further discuss how RNA techniques have contributed to the building of this knowledge and the development of new tools in synthetic biology for the diagnosis and treatment of diseases caused by pathogenic microorganisms. Infectious diseases are still a world-leading cause of death and morbidity, and RNA-based therapeutics have arisen as an alternative way to achieve success. There are still obstacles to overcome in its application, but much progress has been made in a fast and effective manner, paving the way for the solid establishment of RNA-based therapies in the future.

Insights into the Orchestration of Gene Transcription Regulators in Helicobacter pylori

Bacterial pathogens employ a general strategy to overcome host defenses by coordinating the virulence gene expression using dedicated regulatory systems that could raise intricate networks. During the last twenty years, many studies of Helicobacter pylori, a human pathogen responsible for various stomach diseases, have mainly focused on elucidating the mechanisms and functions of virulence factors. In parallel, numerous studies have focused on the molecular mechanisms that regulate gene transcription to attempt to understand the physiological changes of the bacterium during infection and adaptation to the environmental conditions it encounters. The number of regulatory proteins deduced from the genome sequence analyses responsible for the correct orchestration of gene transcription appears limited to 14 regulators and three sigma factors. Furthermore, evidence is accumulating for new and complex circuits regulating gene transcription and H. pylori virulence. Here, we focus on the molecular mechanisms used by H. pylori to control gene transcription as a function of the principal environmental changes.

Role of nickel-regulated small RNA in modulation of Helicobacter pylori virulence factors

Helicobacter pylori (H. pylori) is a Gram-negative bacterium that infects about half of the world's population. H. pylori infection prevails by several mechanisms of adaptation of the bacteria and by its virulence factors including the cytotoxin associated antigen A (CagA). CagA is an oncoprotein that is the protagonist of gastric carcinogenesis associated with prolonged H. pylori infection. In this sense, small regulatory RNAs (sRNAs) are important macromolecules capable of inhibiting and activating gene expression. This function allows sRNAs to act in adjusting to unstable environmental conditions and in responding to cellular stresses in bacterial infections. Recent discoveries have shown that nickel-regulated small RNA (NikS) is a post-transcriptional regulator of virulence properties of H. pylori, including the oncoprotein CagA. Notably, high concentrations of nickel cause the reduction of NikS expression and consequently this increases the levels of CagA. In addition, NikS expression appears to be lower in clinical isolates from patients with gastric cancer when compared to patients without. With that in mind, this minireview approaches, in an accessible way, the most important and current aspects about the role of NikS in the control of virulence factors of H. pylori and the potential clinical repercussions of this modulation.

Profiling the intragenic toxicity determinants of toxin-antitoxin systems: revisiting hok/Sok regulation

Type I toxin-antitoxin systems (T1TAs) are extremely potent bacterial killing systems difficult to characterize using classical approaches. To assess the killing capability of type I toxins and to identify mutations suppressing the toxin expression or activity, we previously developed the FASTBAC-Seq (Functional AnalysiS of Toxin-Antitoxin Systems in BACteria by Deep Sequencing) method in Helicobacter pylori. This method combines a life and death selection with deep sequencing. Here, we adapted and improved our method to investigate T1TAs in the model organism Escherichia coli. As a proof of concept, we revisited the regulation of the plasmidic hok/Sok T1TA system. We revealed the death-inducing phenotype of the Hok toxin when it is expressed from the chromosome in the absence of the antitoxin and recovered previously described intragenic toxicity determinants of this system. We identified nucleotides that are essential for the transcription, translation or activity of Hok. We also discovered single-nucleotide substitutions leading to structural changes affecting either the translation or the stability of the hok mRNA. Overall, we provide the community with an easy-to-use approach to widely characterize TA systems from diverse types and bacteria.

A Tripartite Efflux System Affects Flagellum Stability in Helicobacter pylori

Helicobacter pylori uses a cluster of polar, sheathed flagella for swimming motility. A search for homologs of H. pylori proteins that were conserved in Helicobacter species that possess flagellar sheaths but were underrepresented in Helicobacter species with unsheathed flagella identified several candidate proteins. Four of the identified proteins are predicted to form part of a tripartite efflux system that includes two transmembrane domains of an ABC transporter (HP1487 and HP1486), a periplasmic membrane fusion protein (HP1488), and a TolC-like outer membrane efflux protein (HP1489). Deleting hp1486/hp1487 and hp1489 homologs in H. pylori B128 resulted in reductions in motility and the number of flagella per cell. Cryo-electron tomography studies of intact motors of the Δhp1489 and Δhp1486/hp1487 mutants revealed many of the cells contained a potential flagellum disassembly product consisting of decorated L and P rings, which has been reported in other bacteria. Aberrant motors lacking specific components, including a cage-like structure that surrounds the motor, were also observed in the Δhp1489 mutant. These findings suggest a role for the H. pylori HP1486-HP1489 tripartite efflux system in flagellum stability. Three independent variants of the Δhp1486/hp1487 mutant with enhanced motility were isolated. All three motile variants had the same frameshift mutation in fliL, suggesting a role for FliL in flagellum disassembly.

SAPPHIRE.CNN: Implementation of dRNA-seq-driven, species-specific promoter prediction using convolutional neural networks

Data availability is a consistent bottleneck for the development of bacterial species-specific promoter prediction software. In this work we leverage genome-wide promoter datasets generated with dRNA-seq in the Gram-negative bacteria Pseudomonas aeruginosa and Salmonella enterica for promoter prediction. Convolutional neural networks are presented as an optimal architecture for model training and are further modified and tailored for promoter prediction. The resulting predictors reach high binary accuracies (95% and 94.9%) on test sets and outperform each other when predicting promoters in their associated species. SAPPHIRE.CNN is available online and can also be downloaded to run locally. Our results indicate a dependency of binary promoter classification on an organism’s GC content and a decreased performance of our classifiers on genera they were not trained for, further supporting the need for dedicated, species-specific promoter classification tools.

RNA-seq Sample Preparation Kits Strongly Affect Transcriptome Profiles of a Gas-Fermenting Bacterium

Transcriptome analysis via RNA sequencing (RNA-seq) has become a standard technique employed across various biological fields of study. The rapid adoption of the RNA-seq approach has been mediated, in part, by the development of different commercial RNA-seq library preparation kits compatible with standard next-generation sequencing (NGS) platforms. Generally, the essential steps of library preparation, such as rRNA depletion and first-strand cDNA synthesis, are tailored to a specific group of organisms (e.g., eukaryotes versus prokaryotes) or genomic GC content. Therefore, the selection of appropriate commercial products is of crucial importance to capture the transcriptome of interest as closely to the native state as possible without introduction of technical bias. However, researchers rarely have the resources and time to test various commercial RNA-seq kits for their samples. This work reports a side-by-side comparison of RNA-seq data from Clostridium autoethanogenum obtained using three commercial rRNA removal and strand-specific library construction products of NuGEN Technologies, Qiagen, and Zymo Research and assesses their performance relative to published data. While all three vendors advertise their products as suitable for prokaryotes, we found significant differences in their performance regarding rRNA removal, strand specificity, and most importantly, transcript abundance distribution profiles. Notably, RNA-seq data obtained with Qiagen products were most similar to published data and delivered the best results in terms of library strandedness and transcript abundance distribution range. Our results highlight the importance of finding appropriate organism-specific workflows and library preparation products for RNA-seq studies. IMPORTANCE RNA-seq is a powerful technique for transcriptome profiling while involving elaborate sample processing before library sequencing. We show that RNA-seq library preparation kits can strongly affect the outcome of an RNA-seq experiment. Although library preparation benefits from the availability of various commercial kits, choosing appropriate products for the specific samples can be challenging for new users or for users working with unconventional organisms. Evaluating the performance of different commercial products requires significant financial and time investments infeasible for most researchers. Therefore, users are often guided in their choice of kits by published data involving similar input samples. We conclude that important consideration should be given to selecting sample processing workflows for any given organism.

A systematic review on omics data (metagenomics, metatranscriptomics, and metabolomics) in the role of microbiome in gallbladder disease

Microbiotas are the range of microorganisms (mainly bacteria and fungi) colonizing multicellular, macroscopic organisms. They are crucial for several metabolic functions affecting the health of the host. However, difficulties hamper the investigation of microbiota composition in cultivating microorganisms in standard growth media. For this reason, our knowledge of microbiota can benefit from the analysis of microbial macromolecules (DNA, transcripts, proteins, or by-products) present in various samples collected from the host. Various omics technologies are used to obtain different data. Metagenomics provides a taxonomical profile of the sample. It can also be used to obtain potential functional information. At the same time, metatranscriptomics can characterize members of a microbiome responsible for specific functions and elucidate genes that drive the microbiotas relationship with its host. Thus, while microbiota refers to microorganisms living in a determined environment (taxonomy of microorganisms identified), microbiome refers to the microorganisms and their genes living in a determined environment and, of course, metagenomics focuses on the genes and collective functions of identified microorganisms. Metabolomics completes this framework by determining the metabolite fluxes and the products released into the environment. The gallbladder is a sac localized under the liver in the human body and is difficult to access for bile and tissue sampling. It concentrates the bile produced in the hepatocytes, which drains into bile canaliculi. Bile promotes fat digestion and is released from the gallbladder into the upper small intestine in response to food. Considered sterile originally, recent data indicate that bile microbiota is associated with the biliary tract’s inflammation and carcinogenesis. The sample size is relevant for omic studies of rare diseases, such as gallbladder carcinoma. Although in its infancy, the study of the biliary microbiota has begun taking advantage of several omics strategies, mainly based on metagenomics, metabolomics, and mouse models. Here, we show that omics analyses from the literature may provide a more comprehensive image of the biliary microbiota. We review studies performed in this environmental niche and focus on network-based approaches for integrative studies.

Differentially Expressed Genes of Pseudomonas aeruginosa Isolates from Eyes with Keratitis and Healthy Conjunctival Sacs

Background

Pseudomonas aeruginosa (P. aeruginosa) is the second-most common commensal bacterium in healthy conjunctival sacs. When the corneal epithelial barrier is damaged, P. aeruginosa in a healthy conjunctival sac can cause infectious keratitis, which can result in the loss of vision. This study was designed to investigate the differentially expressed genes (DEGs) of P. aeruginosa isolates from eyes with keratitis and from healthy conjunctival sacs to predict their functions and pathways through Illumina high-throughput RNA sequencing (RNA-seq).

Methods

P. aeruginosa isolates from keratitis and healthy conjunctival sacs were obtained. The transcriptome profile of P. aeruginosa was characterized by a high throughput RNA-seq strategy using the Illumina HiSeq 2500 platform. The DEGs were analyzed with DESeq and validated through quantitative real-time polymerase chain reaction (PCR) and with experimental mice. GO enrichment and the KEGG pathway were also analyzed.

Results

In genome-wide transcriptional analysis, 557 genes (332 upregulated and 225 downregulated) were found to be differentially expressed (fold change ≥ 2, p ≤ 0.05) in the strains from keratitis. GO enrichment analysis suggested that DEGs tended to be associated with cellular and metabolic processes. KEGG pathway analysis revealed the DEGs were typically associated with the pathways of the bacterial secretion system and pyoverdine metabolism. Eleven DEGs were validated using quantitative reverse-transcription PCR and verified with experimental mice. The results were consistent with those obtained in RNA-seq.

Conclusion

The DEGs related to pilin, T2SS, T3SS, and pyoverdine metabolisms were significantly altered in the strains from keratitis. The findings may be helpful for further investigations on genes or pathways related to the pathogenesis of and therapeutic targets for P. aeruginosa keratitis.

Genome-wide transcription start site mapping in the facultative intracellular pathogen Brucella melitensis by Capping-seq

Brucella melitensis (B. melitensis) is an important facultative intracellular bacterium that causes global zoonotic diseases. Continuous intracellular survival and replication are the main obstruction responsible for the accessibility of prevention and treatment of brucellosis. Bacteria respond to complex environment by regulating gene expression. Many regulatory factors function at loci where RNA polymerase initiates messenger RNA synthesis. However, limited gene annotation is a current obstacle for the research on expression regulation in bacteria. To improve annotation and explore potential functional sites, we proposed a novel genome-wide method called Capping-seq for transcription start site (TSS) mapping in B. melitensis. This technique combines capture of capped primary transcripts with Single Molecule Real-Time (SMRT) sequencing technology. We identified 2,369 TSSs at single nucleotide resolution by Capping-seq. TSSs analysis of Brucella transcripts showed a preference of purine on the TSS positions. Our results revealed that -35 and -10 elements of promoter contained consensus sequences of TTGNNN and TATNNN, respectively. The 5' ends analysis showed that 57% genes are associated with more than one TSS and 47% genes contain long leader regions, suggested potential complex regulation at the 5' ends of genes in B. melitensis. Moreover, we identified 52 leaderless genes that are mainly involved in the metabolic processes. Overall, Capping-seq technology provides a unique solution for TSS determination in prokaryotes. Our findings develop a systematic insight into the primary transcriptome characterization of B. melitensis. This study represents a critical basis for investigating gene regulation and pathogenesis of Brucella.

Positive Selection of Mutations in the Helicobacter pylori katA 5' Untranslated Region in a Mongolian Gerbil Model of Gastric Disease

To evaluate potential effects of gastric inflammation on Helicobacter pylori diversification and evolution within the stomach, we experimentally infected Mongolian gerbils with an H. pylori strain in which Cag type IV secretion system (T4SS) activity is controlled by a TetR/tetO system. Gerbils infected with H. pylori under conditions in which Cag T4SS activity was derepressed had significantly higher levels of gastric inflammation than gerbils infected under conditions with repressed Cag T4SS activity. Mutations in the 5' untranslated region (UTR) of katA (encoding catalase) were detected in strains cultured from 8 of the 17 gerbils infected with Cag T4SS-active H. pylori and none of the strains from 17 gerbils infected with Cag T4SS-inactive H. pylori. Catalase enzymatic activity, steady-state katA transcript levels, and katA transcript stability were increased in strains with these single nucleotide polymorphisms (SNPs) compared to strains in which these SNPs were absent. Moreover, strains harboring these SNPs exhibited increased resistance to bactericidal effects of hydrogen peroxide, compared to control strains. Experimental introduction of the SNPs into the wild-type katA 5' UTR resulted in increased katA transcript stability, increased katA steady-state levels, and increased catalase enzymatic activity. Based on site-directed mutagenesis and modeling of RNA structure, increased katA transcript levels were correlated with higher predicted thermal stability of the katA 5' UTR secondary structure. These data suggest that high levels of gastric inflammation positively select for H. pylori strains producing increased levels of catalase, which may confer survival advantages to the bacteria in an inflammatory gastric environment.

Inverse folding based pre-training for the reliable identification of intrinsic transcription terminators

It is well-established that neural networks can predict or identify structural motifs of non-coding RNAs (ncRNAs). Yet, the neural network based identification of RNA structural motifs is limited by the availability of training data that are often insufficient for learning features of specific ncRNA families or structural motifs. Aiming to reliably identify intrinsic transcription terminators in bacteria, we introduce a novel pre-training approach that uses inverse folding to generate training data for predicting or identifying a specific family or structural motif of ncRNA. We assess the ability of neural networks to identify secondary structure by systematic in silico mutagenesis experiments. In a study to identify intrinsic transcription terminators as functionally well-understood RNA structural motifs, our inverse folding based pre-training approach significantly boosts the performance of neural network topologies, which outperform previous approaches to identify intrinsic transcription terminators. Inverse-folding based pre-training provides a simple, yet highly effective way to integrate the well-established thermodynamic energy model into deep neural networks for identifying ncRNA families or motifs. The pre-training technique is broadly applicable to a range of network topologies as well as different types of ncRNA families and motifs.

Intrinsic transcriptional terminators are essential regulators in determining the 3’ end of transcripts in bacteria. The underlying mechanism involves RNA secondary structure, where nucleotides fold into a specific hairpin motif. Identifying terminator sequences in bacterial genomes has conventionally been approached with well-established energy models for structural motifs. However, the folding mechanism of transcription terminators is understood only partially, limiting the success of energy-model based identification.

Neural networks have been proposed to overcome these limitations. However, their adoption for predicting and identifying RNA secondary structure has been a double edged sword: Neural networks promise to learn features that are not represented by the energy models, while they are black boxes that lack explicit modeling assumptions and may fail to account for features that are well understandable based on decades-old energy models.

Here, we introduce a pre-training approach for neural networks that uses energy-model based inverse folding of structural motifs. As we demonstrate, this approach “brings back the energy model” to identify transcriptional terminators and overcomes the limitations of previous energy-model based predictions. Our approach works for diverse types of neural networks, and is suitable for the identification of structural motifs of many other RNA molecules beyond transcriptional terminators.

Posttranscriptional Regulation in Response to Different Environmental Stresses in Campylobacter jejuni

The survival strategies that Campylobacter jejuni (C. jejuni) employ throughout its transmission and infection life cycles remain largely elusive. Specifically, there is a lack of understanding about the posttranscriptional regulation of stress adaptations resulting from small noncoding RNAs (sRNAs). Published C. jejuni sRNAs have been discovered in specific conditions but with limited insights into their biological activities. Many more sRNAs are yet to be discovered as they may be condition-dependent. Here, we have generated transcriptomic data from 21 host- and transmission-relevant conditions. The data uncovered transcription start sites, expression patterns and posttranscriptional regulation during various stress conditions. This data set helped predict a list of putative sRNAs. We further explored the sRNAs' biological functions by integrating differential gene expression analysis, coexpression analysis, and genome-wide sRNA target prediction. The results showed that the C. jejuni gene expression was influenced primarily by nutrient deprivation and food storage conditions. Further exploration revealed a putative sRNA (CjSA21) that targeted tlp1 to 4 under food processing conditions. tlp1 to 4 are transcripts that encode methyl-accepting chemotaxis proteins (MCPs), which are responsible for chemosensing. These results suggested CjSA21 inhibits chemotaxis and promotes survival under food processing conditions. This study presents the broader research community with a comprehensive data set and highlights a novel sRNA as a potential chemotaxis inhibitor. IMPORTANCE The foodborne pathogen C. jejuni is a significant challenge for the global health care system. It is crucial to investigate C. jejuni posttranscriptional regulation by small RNAs (sRNAs) in order to understand how it adapts to different stress conditions. However, limited data are available for investigating sRNA activity under stress. In this study, we generate gene expression data of C. jejuni under 21 stress conditions. Our data analysis indicates that one of the novel sRNAs mediates the adaptation to food processing conditions. Results from our work shed light on the posttranscriptional regulation of C. jejuni and identify an sRNA associated with food safety.

RNase III-CLASH of multi-drug resistant Staphylococcus aureus reveals a regulatory mRNA 3'UTR required for intermediate vancomycin resistance

Treatment of methicillin-resistant Staphylococcus aureus infections is dependent on the efficacy of last-line antibiotics including vancomycin. Treatment failure is commonly linked to isolates with intermediate vancomycin resistance (termed VISA). These isolates have accumulated point mutations that collectively reduce vancomycin sensitivity, often by thickening the cell wall. Changes in regulatory small RNA expression have been correlated with antibiotic stress in VISA isolates however the functions of most RNA regulators is unknown. Here we capture RNA-RNA interactions associated with RNase III using CLASH. RNase III-CLASH uncovers hundreds of novel RNA-RNA interactions in vivo allowing functional characterisation of many sRNAs for the first time. Surprisingly, many mRNA-mRNA interactions are recovered and we find that an mRNA encoding a long 3' untranslated region (UTR) (termed vigR 3'UTR) functions as a regulatory 'hub' within the RNA-RNA interaction network. We demonstrate that the vigR 3'UTR promotes expression of folD and the cell wall lytic transglycosylase isaA through direct mRNA-mRNA base-pairing. Deletion of the vigR 3'UTR re-sensitised VISA to glycopeptide treatment and both isaA and vigR 3'UTR deletions impact cell wall thickness. Our results demonstrate the utility of RNase III-CLASH and indicate that S. aureus uses mRNA-mRNA interactions to co-ordinate gene expression more widely than previously appreciated.

Transcriptome profiles of genes related to growth and virulence potential in Vibrio alginolyticus treated with modified clay

Vibrio alginolyticus is a globally distributed opportunistic pathogen that causes different degrees of disease in various marine organisms, such as fish, shrimp and shellfish. At present, vibriosis caused by V. alginolyticus has a wide epidemic range and causes frequent outbreaks, resulting in substantial losses in aquaculture. According to previous studies, modified clay (MC) could effectively flocculate and reduce the density of Vibrio in water, but it is still unknown whether MC inhibits growth and how it affects virulence in bottom flocs. Here, we studied the response mechanism of V. alginolyticus in flocs treated with MC at the transcriptome level and verified the transcriptomic data combined with relevant physiological experiments and reverse transcription quantitative real-time PCR (RT-qPCR) for the first time. It was found that the morphology of Vibrio in the MC flocs changed, the membrane function was damaged, the antioxidant system was activated, and the material and energy metabolism also changed. In addition, MC could inhibit the expression of virulence factors of V. alginolyticus; for example, flagella, pilus, siderophores, quorum sensing, and other related genes were significantly downregulated. In general, MC effectively inhibited the growth of Vibrio and reduced its virulence potential in flocs, which could provide theoretical support for a new model of healthy aquaculture.

Early and late genome-wide gastric epithelial transcriptome response during infection with the human carcinogen Helicobacterpylori

Infection of the stomach by Helicobacter pylori is a major risk factor for the development of gastric cancer. Colonization of the gastric epithelium leads to the activation of multiple disease-related signaling pathways. Serine protease HtrA represents an important secreted virulence factor that mediates cleavage of cellular junctions. However, its potential role in nuclear responses is unknown. Here, we performed a genome-wide RNA-seq analysis of polarized gastric epithelial cells infected by wild-type (wt) and ΔhtrA mutant bacteria. Fluorescence microscopy showed that H. pylori wt, but not ΔhtrA bacteria, preferably localized at cellular junctions. Our results pinpointed early (2 h) and late (6 h) transcriptional responses, with most differentially expressed genes at 6 h post infection. The transcriptomes revealed HtrA-dependent targeting of genes associated with inflammation and apoptosis (e.g. IL8, ZFP36, TNF). Accordingly, infection with the ΔhtrA mutant induced increased apoptosis rates in host cells, which was associated with reduced H. pylori CagA expression. In contrast, transcription of various carcinogenesis-associated genes (e.g. DKK1, DOCK8) was affected by H. pylori independent of HtrA. These findings suggest that H. pylori disturbs previously unknown molecular pathways in an HtrA-dependent and HtrA-independent manner, and provide valuable new insights of this significant pathogen in humans and thus potential targets for better controlling the risk of malignant transformation.

Ms1 RNA Interacts With the RNA Polymerase Core in Streptomyces coelicolor and Was Identified in Majority of Actinobacteria Using a Linguistic Gene Synteny Search

Bacteria employ small non-coding RNAs (sRNAs) to regulate gene expression. Ms1 is an sRNA that binds to the RNA polymerase (RNAP) core and affects the intracellular level of this essential enzyme. Ms1 is structurally related to 6S RNA that binds to a different form of RNAP, the holoenzyme bearing the primary sigma factor. 6S RNAs are widespread in the bacterial kingdom except for the industrially and medicinally important Actinobacteria. While Ms1 RNA was identified in Mycobacterium, it is not clear whether Ms1 RNA is present also in other Actinobacteria species. Here, using a computational search based on secondary structure similarities combined with a linguistic gene synteny approach, we identified Ms1 RNA in Streptomyces. In S. coelicolor, Ms1 RNA overlaps with the previously annotated scr3559 sRNA with an unknown function. We experimentally confirmed that Ms1 RNA/scr3559 associates with the RNAP core without the primary sigma factor HrdB in vivo. Subsequently, we applied the computational approach to other Actinobacteria and identified Ms1 RNA candidates in 824 Actinobacteria species, revealing Ms1 RNA as a widespread class of RNAP binding sRNAs, and demonstrating the ability of our multifactorial computational approach to identify weakly conserved sRNAs in evolutionarily distant genomes.

Antimicrobial Resistance Profile by Metagenomic and Metatranscriptomic Approach in Clinical Practice: Opportunity and Challenge

The burden of bacterial resistance to antibiotics affects several key sectors in the world, including healthcare, the government, and the economic sector. Resistant bacterial infection is associated with prolonged hospital stays, direct costs, and costs due to loss of productivity, which will cause policy makers to adjust their policies. Current widely performed procedures for the identification of antibiotic-resistant bacteria rely on culture-based methodology. However, some resistance determinants, such as free-floating DNA of resistance genes, are outside the bacterial genome, which could be potentially transferred under antibiotic exposure. Metagenomic and metatranscriptomic approaches to profiling antibiotic resistance offer several advantages to overcome the limitations of the culture-based approach. These methodologies enhance the probability of detecting resistance determinant genes inside and outside the bacterial genome and novel resistance genes yet pose inherent challenges in availability, validity, expert usability, and cost. Despite these challenges, such molecular-based and bioinformatics technologies offer an exquisite advantage in improving clinicians’ diagnoses and the management of resistant infectious diseases in humans. This review provides a comprehensive overview of next-generation sequencing technologies, metagenomics, and metatranscriptomics in assessing antimicrobial resistance profiles.

Genome-Wide TSS Distribution in Three Related Clostridia with Normalized Capp-Switch Sequencing

Transcription initiation is a tightly regulated process that is crucial for many aspects of prokaryotic physiology. High-throughput transcription start site (TSS) mapping can shed light on global and local regulation of transcription initiation, which in turn may help us understand and predict microbial behavior. In this study, we used Capp-Switch sequencing to determine the TSS positions in the genomes of three model solventogenic clostridia: Clostridium acetobutylicum ATCC 824, C. beijerinckii DSM 6423, and C. beijerinckii NCIMB 8052. We first refined the approach by implementing a normalization pipeline accounting for gene expression, yielding a total of 12,114 mapped TSSs across the species. We further compared the distributions of these sites in the three strains. Results indicated similar distribution patterns at the genome scale, but also some sharp differences, such as for the butyryl-CoA synthesis operon, particularly when comparing C. acetobutylicum to the C. beijerinckii strains. Lastly, we found that promoter structure is generally poorly conserved between C. acetobutylicum and C. beijerinckii. A few conserved promoters across species are discussed, showing interesting examples of how TSS determination and comparison can improve our understanding of gene expression regulation at the transcript level. IMPORTANCE Solventogenic clostridia have been employed in industry for more than a century, initially being used in the acetone-butanol-ethanol (ABE) fermentation process for acetone and butanol production. Interest in these bacteria has recently increased in the context of green chemistry and sustainable development. However, our current understanding of their genomes and physiology limits their optimal use as industrial solvent production platforms. The gene regulatory mechanisms of solventogenesis are still only partly understood, impeding efforts to increase rates and yields. Genome-wide mapping of transcription start sites (TSSs) for three model solventogenic Clostridium strains is an important step toward understanding mechanisms of gene regulation in these industrially important bacteria.

Small RNAs beyond Model Organisms: Have We Only Scratched the Surface?

Small RNAs (sRNAs) are essential regulators in the adaptation of bacteria to environmental changes and act by binding targeted mRNAs through base complementarity. Approximately 550 distinct families of sRNAs have been identified since their initial characterization in the 1980s, accelerated by the emergence of RNA-sequencing. Small RNAs are found in a wide range of bacterial phyla, but they are more prominent in highly researched model organisms compared to the rest of the sequenced bacteria. Indeed, Escherichia coli and Salmonella enterica contain the highest number of sRNAs, with 98 and 118, respectively, with Enterobacteriaceae encoding 145 distinct sRNAs, while other bacteria families have only seven sRNAs on average. Although the past years brought major advances in research on sRNAs, we have perhaps only scratched the surface, even more so considering RNA annotations trail behind gene annotations. A distinctive trend can be observed for genes, whereby their number increases with genome size, but this is not observable for RNAs, although they would be expected to follow the same trend. In this perspective, we aimed at establishing a more accurate representation of the occurrence of sRNAs in bacteria, emphasizing the potential for novel sRNA discoveries.

The Small RNA CyaR Activates Translation of the Outer Membrane Haem Receptor chuA in Enterohemorrhagic Escherichia coli

To sense the transition from environment to host, bacteria use a range of environmental cues to control expression of virulence genes. Iron is tightly sequestered in host tissues and in the human pathogen enterohemorrhagic Escherichia coli (EHEC) iron-limitation induces transcription of the outer membrane haem transporter encoded by chuAS. ChuA expression is post-transcriptionally activated at 37°C by a FourU RNA thermometer ensuring that the haem receptor is only expressed under low iron, high temperature conditions that indicate the host. Here we demonstrate that expression of chuA is also independently regulated by the cAMP-responsive small RNA (sRNA) CyaR and transcriptional terminator Rho. These results indicate that chuAS expression is regulated at the transcription initiation, transcript elongation, and translational level. We speculate that additional sensing of the gluconeogenic environment allows further precision in determining when EHEC is at the gastrointestinal epithelium of the host. With previous studies, it appears that the chuAS transcript is controlled by eight regulatory inputs that control expression through six different transcriptional and post-transcriptional mechanisms. The results highlight the ability of regulatory sRNAs to integrate multiple environmental signals into a layered hierarchy of signal input.

An overview of gene regulation in bacteria by small RNAs derived from mRNA 3' ends

Over the past two decades, small noncoding RNAs (sRNAs) that regulate mRNAs by short base pairing have gone from a curiosity to a major class of post-transcriptional regulators in bacteria. They are integral to many stress responses and regulatory circuits, affecting almost all aspects of bacterial life. Following pioneering sRNA searches in the early 2000s, the field quickly focused on conserved sRNA genes in the intergenic regions of bacterial chromosomes. Yet, it soon emerged that there might be another rich source of bacterial sRNAs—processed 3′ end fragments of mRNAs. Several such 3′ end-derived sRNAs have now been characterized, often revealing unexpected, conserved functions in diverse cellular processes. Here, we review our current knowledge of these 3′ end-derived sRNAs—their biogenesis through ribonucleases, their molecular mechanisms, their interactions with RNA-binding proteins such as Hfq or ProQ and their functional scope, which ranges from acting as specialized regulators of single metabolic genes to constituting entire noncoding arms in global stress responses. Recent global RNA interactome studies suggest that the importance of functional 3′ end-derived sRNAs has been vastly underestimated and that this type of cross-regulation between genes at the mRNA level is more pervasive in bacteria than currently appreciated.

AI-2 Induces Urease Expression Through Downregulation of Orphan Response Regulator HP1021 in Helicobacter pylori

Helicobacter pylori causes gastric infections in more than half of the world's population. The bacterium's survival in the stomach is mediated by the abundant production of urease to enable acid acclimation. In this study, our transcriptomic analysis demonstrated that the expression of urease structural proteins, UreA and UreB, is induced by the autoinducer AI-2 in H. pylori. We also found that the orphan response regulator HP1021 is downregulated by AI-2, resulting in the induction of urease expression. HP1021 represses the expression of urease by directly binding to the promoter region of ureAB, ranging from −47 to +3 with respect to the transcriptional start site. The study findings suggest that quorum sensing via AI-2 enhances acid acclimation when bacterial density increases, and might enable bacterial dispersal to other sites when entering gastric acid.

Synthetic 3'-UTR valves for optimal metabolic flux control in Escherichia coli

As the design of genetic circuitry for synthetic biology becomes more sophisticated, diverse regulatory bioparts are required. Despite their importance, well-characterized 3′-untranslated region (3′-UTR) bioparts are limited. Thus, transcript 3′-ends require further investigation to understand the underlying regulatory role and applications of the 3′-UTR. Here, we revisited the use of Term-Seq in the Escherichia coli strain K-12 MG1655 to enhance our understanding of 3′-UTR regulatory functions and to provide a diverse collection of tunable 3′-UTR bioparts with a wide termination strength range. Comprehensive analysis of 1,629 transcript 3′-end positions revealed multiple 3′-termini classes generated through transcription termination and RNA processing. The examination of individual Rho-independent terminators revealed a reduction in downstream gene expression over a wide range, which led to the design of novel synthetic metabolic valves that control metabolic fluxes in branched pathways. These synthetic metabolic valves determine the optimal balance of heterologous pathways for maximum target biochemical productivity. The regulatory strategy using 3′-UTR bioparts is advantageous over promoter- or 5′-UTR-based transcriptional control as it modulates gene expression at transcription levels without trans-acting element requirements (e.g. transcription factors). Our results provide a foundational platform for 3′-UTR engineering in synthetic biology applications.

Spacer prioritization in CRISPR-Cas9 immunity is enabled by the leader RNA

CRISPR-Cas systems store fragments of foreign DNA called spacers as immunological recordings used to combat future infections. Of the many spacers stored in a CRISPR array, the newest spacers are known to be prioritized for immune defense. However, the underlying mechanism remains unclear. Here we show that the leader region upstream of CRISPR arrays in CRISPR-Cas9 systems enhances CRISPR RNA (crRNA) processing from the newest spacer, prioritizing defense against the matching invader. Using the CRISPR-Cas9 system from Streptococcus pyogenes as a model, we found that the transcribed leader interacts with the conserved repeats bordering the newest spacer. The resulting interaction promotes tracrRNA hybridization with the second repeat, accelerating crRNA processing. Accordingly, disrupting this structure reduces the abundance of the associated crRNA and immune defense against targeted plasmids and bacteriophages. Beyond the S. pyogenes system, bioinformatics analyses revealed that leader-repeat structures appear across CRISPR-Cas9 systems. CRISPR-Cas systems thus possess an RNA-based mechanism to prioritize defense against the most recently encountered invaders.

Identification of Antimotilins, Novel Inhibitors of Helicobacter pylori Flagellar Motility That Inhibit Stomach Colonization in a Mouse Model

New treatment options against the widespread cancerogenic gastric pathogen Helicobacter pylori are urgently needed. We describe a novel screening procedure for inhibitors of H. pylori flagellar biosynthesis. The assay is based on a flaA flagellin gene-luciferase reporter fusion in H. pylori and was amenable to multi-well screening formats with an excellent Z factor. We screened various compound libraries to identify virulence blockers ("antimotilins") that inhibit H. pylori motility or the flagellar type III secretion apparatus. We identified compounds that either inhibit both motility and the bacterial viability, or the flagellar system only, without negatively affecting bacterial growth. Novel anti-virulence compounds which suppressed flagellar biosynthesis in H. pylori were active on pure H. pylori cultures in vitro and partially suppressed motility directly, reduced flagellin transcript and flagellin protein amounts. We performed a proof-of-principle treatment study in a mouse model of chronic H. pylori infection and demonstrated a significant effect on H. pylori colonization for one antimotilin termed Active2 even as a monotherapy. The diversity of the intestinal microbiota was not significantly affected by Active2. In conclusion, the novel antimotilins active against motility and flagellar assembly bear promise to complement commonly used antibiotic-based combination therapies for treating and eradicating H. pylori infections. IMPORTANCE Helicobacter pylori is one of the most prevalent bacterial pathogens, inflicting hundreds of thousands of peptic ulcers and gastric cancers to patients every year. Antibacterial treatment of H. pylori is complicated due to the need of combining multiple antibiotics, entailing serious side effects and increasing selection for antibiotic resistance. Here, we aimed to explore novel nonantibiotic approaches to H. pylori treatment. We selected an antimotility approach since flagellar motility is essential for H. pylori colonization. We developed a screening system for inhibitors of H. pylori motility and flagellar assembly, and identified numerous novel antibacterial and anti-motility compounds (antimotilins). Selected compounds were further characterized, and one was evaluated in a preclinical therapy study in mice. The antimotilin compound showed a good efficacy to reduce bacterial colonization in the model, such that the antimotilin approach bears promise to be further developed into a therapy against H. pylori infection in humans.

Is Laryngeal Squamous Cell Carcinoma Related to Helicobacter pylori?

Background

Laryngeal carcinoma is a primary malignant tumor originating from the laryngeal mucosa. In recent years, an increasing number of studies have confirmed that Helicobacter pylori may play a role in the occurrence and development of laryngeal cancer. We conducted a systematic review and meta-analysis to identify and emphasize the relationship between laryngeal cancer and Helicobacter pylori infection.

Methods

We actively searched for systematic reviews of PubMed, Embase, Web of Science, and Cochrane libraries to select studies that met the recommended guidelines. A total of 1756 patients with laryngeal cancer were included in this study to assess the association of Helicobacter pylori in the larynx with laryngeal cancer. A subassessment of the risk of bias for each study that met the inclusion criteria was carried out. To illustrate the reasons for heterogeneity, we performed a subgroup analysis to determine the type of study, the quality of the article, the diagnostic method, and the impact of exposure factors.

Results

The meta-analysis included a total of 17 case-control studies on the association between Helicobacter pylori in the larynx and laryngeal cancer. Our meta-analysis showed that Helicobacter pylori infection in the larynx significantly increased the risk of malignant tumors in the larynx (odds ratio, 2.96; 95% CI, 1.83-4.78; P<0.00001; I^{2 =} 86%). They still existed when we controlled for patients’ exposure to smoking factors (odds ratio, 3.86; 95% CI, 1.89–7.88).

Conclusion

This systematic review and meta-analysis showed an association between Helicobacter pylori and laryngeal cancer. These findings are consistent with the understanding that chronic inflammatory tissue insult may lead to the development of malignancy. Controlling risk factors will help us identify patients with an increased risk of laryngeal cancer.

Genome-scale analysis of genetic regulatory elements in Streptomyces avermitilis MA-4680 using transcript boundary information

Background

The gram-positive bacterium, Streptomyces avermitilis, holds industrial importance as the producer of avermectin, a widely used anthelmintic agent, and a heterologous expression host of secondary metabolite-biosynthetic gene clusters. Despite its industrial importance, S. avermitilis’ genome organization and regulation of gene expression remain poorly understood. In this study, four different types of Next-Generation Sequencing techniques, including dRNA-Seq, Term-Seq, RNA-Seq and ribosome profiling, were applied to S. avermitilis to determine transcription units of S. avermitilis at a genome-wide level and elucidate regulatory elements for transcriptional and translational control of individual transcription units.

Result

By applying dRNA-Seq and Term-Seq to S. avermitilis MA-4680, a total of 2361 transcription start sites and 2017 transcript 3′-end positions were identified, respectively, leading to determination of 1601 transcription units encoded in S. avermitilis’ genome. Cataloguing the transcription units and integrated analysis of multiple high-throughput data types revealed the presence of diverse regulatory elements for gene expression, such as promoters, 5′-UTRs, terminators, 3′-UTRs and riboswitches. The conserved promoter motifs were identified from 2361 transcription start sites as 5′-TANNNT and 5′-BTGACN for the − 10 and − 35 elements, respectively. The − 35 element and spacer lengths between − 10 and − 35 elements were critical for transcriptional regulation of functionally distinct genes, suggesting the involvement of unique sigma factors. In addition, regulatory sequences recognized by antibiotic regulatory proteins were identified from the transcription start site information. Analysis of the 3′-end of RNA transcript revealed that stem structure formation is a major determinant for transcription termination of most transcription units.

Conclusions

The transcription unit architecture elucidated from the transcripts’ boundary information provides insights for unique genetic regulatory mechanisms of S. avermitilis. Our findings will elevate S. avermitilis’ potential as a production host for a diverse set of secondary metabolites.

A Practical Guide to Small Protein Discovery and Characterization Using Mass Spectrometry

Small proteins of up to ∼50 amino acids are an abundant class of biomolecules across all domains of life. Yet due to the challenges inherent in their size, they are often missed in genome annotations, and are difficult to identify and characterize using standard experimental approaches. Consequently, we still know few small proteins even in well-studied prokaryotic model organisms. Mass spectrometry (MS) has great potential for the discovery, validation, and functional characterization of small proteins. However, standard MS approaches are poorly suited to the identification of both known and novel small proteins due to limitations at each step of a typical proteomics workflow, i.e., sample preparation, protease digestion, liquid chromatography, MS data acquisition, and data analysis. Here, we outline the major MS-based workflows and bioinformatic pipelines used for small protein discovery and validation. Special emphasis is placed on highlighting the adjustments required to improve detection and data quality for small proteins. We discuss both the unbiased detection of small proteins and the targeted analysis of small proteins of interest. Finally, we provide guidelines to prioritize novel small proteins, and an outlook on methods with particular potential to further improve comprehensive discovery and characterization of small proteins.

Regulation of Heterogenous LexA Expression in Staphylococcus aureus by an Antisense RNA Originating from Transcriptional Read-Through upon Natural Mispairings in the sbrB Intrinsic Terminator

Bacterial genomes are pervasively transcribed, generating a wide variety of antisense RNAs (asRNAs). Many of them originate from transcriptional read-through events (TREs) during the transcription termination process. Previous transcriptome analyses revealed that the lexA gene from Staphylococcus aureus, which encodes the main SOS response regulator, is affected by the presence of an asRNA. Here, we show that the lexA antisense RNA (lexA-asRNA) is generated by a TRE on the intrinsic terminator (TT_sbrB) of the sbrB gene, which is located downstream of lexA, in the opposite strand. Transcriptional read-through occurs by a natural mutation that destabilizes the TT_sbrB structure and modifies the efficiency of the intrinsic terminator. Restoring the mispairing mutation in the hairpin of TT_sbrB prevented lexA-asRNA transcription. The level of lexA-asRNA directly correlated with cellular stress since the expressions of sbrB and lexA-asRNA depend on the stress transcription factor SigB. Comparative analyses revealed strain-specific nucleotide polymorphisms within TT_sbrB, suggesting that this TT could be prone to accumulating natural mutations. A genome-wide analysis of TREs suggested that mispairings in TT hairpins might provide wider transcriptional connections with downstream genes and, ultimately, transcriptomic variability among S. aureus strains.

Metal-Organic-Framework-Based Hydrogen-Release Platform for Multieffective Helicobacter Pylori Targeting Therapy and Intestinal Flora Protective Capabilities

Helicobacter pylori (H. pylori) infection is the leading cause of chronic gastritis, peptic ulcer, and gastric cancer. Antibiotics, as traditional method for eliminating H. pylori, have no targeting effect, which causes serious bacterial resistance and gut dysbacteriosis. Moreover, antibiotics can hardly address hyperactive inflammatory response or damaged gastric mucosal barrier caused by H. pylori infection. Here, a pH-responsive metal-organic framework hydrogen-generation nanoparticle (Pd(H) @ ZIF-8) is reported, which is encapsulated with ascorbate palmitate (AP) hydrogel. Both in vitro and in vivo experiments demonstrate that the outer AP hydrogel can target and adhere to the inflammatory site through electrostatic interactions, and is then hydrolyzed by matrix metalloproteinase (MMP) enriching in inflammatory sites. The released Pd(H) @ ZIF-8 nanoparticles are further decomposed by gastric acid to generate zinc ions (Zn²⁺ ) and hydrogen, thus effectively killing H. pylori, alleviating inflammation and restoring impaired gastric mucosa simultaneously. Unexpectedly, this metal-organic framework hydrogen-generation platform (Pd(H) @ ZIF-8 @ AP) also has an effect toward avoiding the imbalance of intestinal flora, which thus provides a more precise, effective, and healthy strategy for the treatment of H. pylori infection.

Identification of RNA Binding Partners of CRISPR-Cas Proteins in Prokaryotes Using RIP-Seq

CRISPR-Cas systems consist of a complex ribonucleoprotein (RNP) machinery encoded in prokaryotic genomes to confer adaptive immunity against foreign mobile genetic elements. Of these, especially the class 2, Type II CRISPR-Cas9 RNA-guided systems with single protein effector modules have recently received much attention for their application as programmable DNA scissors that can be used for genome editing in eukaryotes. While many studies have concentrated their efforts on improving RNA-mediated DNA targeting with these Type II systems, little is known about the factors that modulate processing or binding of the CRISPR RNA (crRNA) guides and the trans-activating tracrRNA to the nuclease protein Cas9, and whether Cas9 can also potentially interact with other endogenous RNAs encoded within the host genome. Here, we describe RIP-seq as a method to globally identify the direct RNA binding partners of CRISPR-Cas RNPs using the Cas9 nuclease as an example. RIP-seq combines co-immunoprecipitation (coIP) of an epitope-tagged Cas9 followed by isolation and deep sequencing analysis of its co-purified bound RNAs. This method can not only be used to study interactions of Cas9 with its known interaction partners, crRNAs and tracrRNA in native systems, but also to reveal potential additional RNA substrates of Cas9. For example, in RIP-seq analysis of Cas9 from the foodborne pathogen Campylobacter jejuni (CjeCas9), we recently identified several endogenous RNAs bound to CjeCas9 RNP in a crRNA-dependent manner, leading to the discovery of PAM-independent RNA cleavage activity of CjeCas9 as well as non-canonical crRNAs. RIP-seq can be easily adapted to any other effector RNP of choice from other CRISPR-Cas systems, allowing for the identification of target RNAs. Deciphering novel RNA-protein interactions for CRISPR-Cas proteins within host bacterial genomes will lead to a better understanding of the molecular mechanisms and functions of these systems and enable us to use the in vivo identified interaction rules as design principles for nucleic acid-targeting applications, fitted to each nuclease of interest.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

Challenges in defining the functional, non-coding, expressed genome of members of the Mycobacterium tuberculosis complex

A definitive transcriptome atlas for the non-coding expressed elements of the members of the Mycobacterium tuberculosis complex (MTBC) does not exist. Incomplete lists of non-coding transcripts can be obtained for some of the reference genomes (e.g., M. tuberculosis H37Rv) but to what extent these transcripts have homologues in closely related species or even strains is not clear. This has implications for the analysis of transcriptomic data; non-coding parts of the transcriptome are often ignored in the absence of formal, reliable annotation. Here, we review the state of our knowledge of non-coding RNAs in pathogenic mycobacteria, emphasizing the disparities in the information included in commonly used databases. We then proceed to review ways of combining computational solutions for predicting the non-coding transcriptome with experiments that can help refine and confirm these predictions.

T1TAdb: the database of type I toxin-antitoxin systems

Type I toxin-antitoxin (T1TA) systems constitute a large class of genetic modules with antisense RNA (asRNA)-mediated regulation of gene expression. They are widespread in bacteria and consist of an mRNA coding for a toxic protein and a noncoding asRNA that acts as an antitoxin preventing the synthesis of the toxin by directly base-pairing to its cognate mRNA. The co- and post-transcriptional regulation of T1TA systems is intimately linked to RNA sequence and structure, therefore it is essential to have an accurate annotation of the mRNA and asRNA molecules to understand this regulation. However, most T1TA systems have been identified by means of bioinformatic analyses solely based on the toxin protein sequences, and there is no central repository of information on their specific RNA features. Here we present the first database dedicated to type I TA systems, named T1TAdb. It is an open-access web database (https://d-lab.arna.cnrs.fr/t1tadb) with a collection of ∼1900 loci in ∼500 bacterial strains in which a toxin-coding sequence has been previously identified. RNA molecules were annotated with a bioinformatic procedure based on key determinants of the mRNA structure and the genetic organization of the T1TA loci. Besides RNA and protein secondary structure predictions, T1TAdb also identifies promoter, ribosome-binding, and mRNA-asRNA interaction sites. It also includes tools for comparative analysis, such as sequence similarity search and computation of structural multiple alignments, which are annotated with covariation information. To our knowledge, T1TAdb represents the largest collection of features, sequences, and structural annotations on this class of genetic modules.

Genome-wide association study of gastric cancer- and duodenal ulcer-derived Helicobacter pylori strains reveals discriminatory genetic variations and novel oncoprotein candidates

Genome-wide association studies (GWASs) can reveal genetic variations associated with a phenotype in the absence of any hypothesis of candidate genes. The problem of false-positive sites linked with the responsible site might be bypassed in bacteria with a high homologous recombination rate, such as Helicobacter pylori, which causes gastric cancer. We conducted a small-sample GWAS (125 gastric cancer cases and 115 controls) followed by prediction of gastric cancer and control (duodenal ulcer) H. pylori strains. We identified 11 single nucleotide polymorphisms (eight amino acid changes) and three DNA motifs that, combined, allowed effective disease discrimination. They were often informative of the underlying molecular mechanisms, such as electric charge alteration at the ligand-binding pocket, alteration in subunit interaction, and mode-switching of DNA methylation. We also identified three novel virulence factors/oncoprotein candidates. These results provide both defined targets for further informatic and experimental analyses to gain insights into gastric cancer pathogenesis and a basis for identifying a set of biomarkers for distinguishing these H. pylori-related diseases.

Methylation Motifs in Promoter Sequences May Contribute to the Maintenance of a Conserved m5C Methyltransferase in Helicobacter pylori

DNA methylomes of Helicobacter pylori strains are complex due to the large number of DNA methyltransferases (MTases) they possess. H. pylori J99 M.Hpy99III is a 5-methylcytosine (^m5C) MTase that converts GCGC motifs to G^m5CGC. Homologs of M.Hpy99III are found in essentially all H. pylori strains. Most of these homologs are orphan MTases that lack a cognate restriction endonuclease, and their retention in H. pylori strains suggest they have roles in gene regulation. To address this hypothesis, green fluorescent protein (GFP) reporter genes were constructed with six putative promoters that had a GCGC motif in the extended −10 region, and the expression of the reporter genes was compared in wild-type H. pylori G27 and a mutant lacking the M.Hpy99III homolog (M.HpyGIII). The expression of three of the GFP reporter genes was decreased significantly in the mutant lacking M.HpyGIII. In addition, the growth rate of the H. pylori G27 mutant lacking M.HpyGIII was reduced markedly compared to that of the wild type. These findings suggest that the methylation of the GCGC motif in many H. pylori GCGC-containing promoters is required for the robust expression of genes controlled by these promoters, which may account for the universal retention of M.Hpy99III homologs in H. pylori strains.

RNase III-mediated processing of a trans-acting bacterial sRNA and its cis-encoded antagonist

Bacterial small RNAs (sRNAs) are important post-transcriptional regulators in stress responses and virulence. They can be derived from an expanding list of genomic contexts, such as processing from parental transcripts by RNase E. The role of RNase III in sRNA biogenesis is less well understood despite its well-known roles in rRNA processing, RNA decay, and cleavage of sRNA-mRNA duplexes. Here, we show that RNase III processes a pair of cis-encoded sRNAs (CJnc190 and CJnc180) of the food-borne pathogen Campylobacter jejuni. While CJnc180 processing by RNase III requires CJnc190, RNase III processes CJnc190 independent of CJnc180 via cleavage of an intramolecular duplex. We also show that CJnc190 directly represses translation of the colonization factor PtmG by targeting a G-rich ribosome-binding site, and uncover that CJnc180 is a cis-acting antagonist of CJnc190, indirectly affecting ptmG regulation. Our study highlights a role for RNase III in sRNA biogenesis and adds cis-encoded RNAs to the expanding diversity of transcripts that can antagonize bacterial sRNAs.

Small RNAs that target G-rich sequences are generated by diverse biogenesis pathways in Epsilonproteobacteria

Bacterial small RNAs (sRNAs) are widespread post-transcriptional regulators controlling bacterial stress responses and virulence. Nevertheless, little is known about how they arise and evolve. Homologues can be difficult to identify beyond the strain level using sequence-based approaches, and similar functionalities can arise by convergent evolution. Here, we found that the virulence-associated CJnc190 sRNA of the foodborne pathogen Campylobacter jejuni resembles the RepG sRNA from the gastric pathogen Helicobacter pylori. However, while both sRNAs bind G-rich sites in their target mRNAs using a C/U-rich loop, they largely differ in their biogenesis. RepG is transcribed from a stand-alone gene and does not require processing, whereas CJnc190 is transcribed from two promoters as precursors that are processed by RNase III and also has a cis-encoded antagonist, CJnc180. By comparing CJnc190 homologues in diverse Campylobacter species, we show that RNase III-dependent processing of CJnc190 appears to be a conserved feature even outside of C. jejuni. We also demonstrate the CJnc180 antisense partner is expressed in C. coli, yet here might be derived from the 3'UTR of the upstream flagella-related gene. Our analysis of G-tract targeting sRNAs in Epsilonproteobacteria demonstrates that similar sRNAs can have markedly different biogenesis pathways.