High-resolution transcriptome maps reveal strain-specific regulatory features of multiple Campylobacter jejuni isolates

Interplay of two small RNAs fine-tunes hierarchical flagella gene expression in Campylobacter jejuni

Like for many bacteria, flagella are crucial for Campylobacter jejuni motility and virulence. Biogenesis of the flagellar machinery requires hierarchical transcription of early, middle (RpoN-dependent), and late (FliA-dependent) genes. However, little is known about post-transcriptional regulation of flagellar biogenesis by small RNAs (sRNAs). Here, we characterized two sRNAs with opposing effects on C. jejuni filament assembly and motility. We demonstrate that CJnc230 sRNA (FlmE), encoded downstream of the flagellar hook protein, is processed from the RpoN-dependent flgE mRNA by RNase III, RNase Y, and PNPase. We identify mRNAs encoding a flagella-interaction regulator and the anti-sigma factor FlgM as direct targets of CJnc230 repression. CJnc230 overexpression upregulates late genes, including the flagellin flaA, culminating in longer flagella and increased motility. In contrast, overexpression of the FliA-dependent sRNA CJnc170 (FlmR) reduces flagellar length and motility. Overall, our study demonstrates how the interplay of two sRNAs post-transcriptionally fine-tunes flagellar biogenesis through balancing of the hierarchically-expressed components.

Dissecting components of the Campylobacter jejuni fetMP-fetABCDEF gene cluster under iron limitation

IMPORTANCE

Campylobacter jejuni is a bacterium that is prevalent in the ceca of farmed poultry such as chickens. Consumption of ill-prepared poultry is thus the most common route by which C. jejuni infects the human gut to cause a typically self-limiting but severe gastrointestinal illness that can be fatal to very young, old, or immunocompromised people. The lack of a vaccine and an increasing resistance to current antibiotics highlight a need to better understand the mechanisms that make C. jejuni a successful human pathogen. This study focused on the functional components of one such mechanism-a molecular system that helps C. jejuni thrive despite the restriction on growth-available iron by the human body, which typically defends against pathogens. In providing a deeper understanding of how this system functions, this study contributes toward the goal of reducing the enormous global socioeconomic burden caused by C. jejuni.

Discovering Novel Bacterial Small RNA by RNA-seq Analysis Toolkit ANNOgesic

ANNOgesic is an RNA-seq analysis pipeline that can detect sRNAs and many other genomic features in bacteria and archaea. In addition to listing sRNA candidates, ANNOgesic also generates various formats of data files for visual examination and downstream experimental design. Based on validations from previous studies, the sRNA predictions are accurate and reliable. In this chapter, we outline the sRNA detection algorithm, important parameters used, step-by-step execution, and data interpretation with a B. pertussis study as an example. Following those procedures, novel sRNA can be revealed by ANNOgesic.

Evaluation of Genetic Content of the CRISPR Locus in Listeria monocytogenes Isolated From Clinical, Food, Seafood and Animal Samples in Iran

CRISPR arrays, which are organized to fight against non-self DNA elements, have shown sequence diversity that could be useful in evolution and typing studies. In this study, 55 samples of L. monocytogenes isolated from different sources were evaluated for CRISPR sequence polymorphism. The CRISPR loci were identified using CRISPR databases. A single PCR assay was designed to amplify the target CRISPRs using an appropriate universal primer. Sequencing results were analyzed using CRISPR databases and BLASTn, and the CRISPR locus was present in all the strains. Three hundred repeats including 55 terminal repeats were identified. Four types of consensuses direct repeat (DR) with different lengths and sequences were characterized. Sixty repeat variants were observed which possessed different polymorphisms. Two hundred and fifty spacers were identified from which 35 consensus sequences were determined, indicating the high polymorphism of the CRISPR spacers. The identified spacers showed similarities to listeria phage sequences, other bacterial phage sequences, plasmid sequences and bacterial sequences. In order to control the bacterial outbreaks, a robust and precise system of subtyping is required. High levels of polymorphism in the CRISPR loci in this study might be related to the origin and time of the samples' isolation. However, it is essential to assess, on a case-by-case basis, the characteristics of any given CRISPR locus before its use as an epidemiological marker. In conclusion, the results of this study showed that the use of sequence content of CRISPR area could provide new and valuable information on the evolution and typing of the L. monocytogenes bacterium.

Direct interaction of small non-coding RNAs CjNC140 and CjNC110 optimizes expression of key pathogenic phenotypes of Campylobacter jejuni

Small non-coding RNAs (sRNAs) are important players in modulating gene expression in bacterial pathogens, but their functions are largely undetermined in Campylobacter jejuni, an important cause of foodborne gastroenteritis in humans. In this study, we elucidated the functions of sRNA CjNC140 and its interaction with CjNC110, a previously characterized sRNA involved in the regulation of several virulence phenotypes of C. jejuni. Inactivation of CjNC140 increased motility, autoagglutination, L-methionine concentration, autoinducer-2 production, hydrogen peroxide resistance, and early chicken colonization, indicating a primarily inhibitory role of CjNC140 for these phenotypes. Apart from motility, all these effects directly contrasted the previously demonstrated positive regulation by CjNC110, suggesting that CjNC110 and CjNC140 operate in an opposite manner to modulate physiologic processes in C. jejuni. RNAseq and northern blotting further demonstrated that expression of CjNC140 increased in the absence of CjNC110, while expression of CjNC110 decreased in the absence of CjNC140, suggesting a possibility of their direct interaction. Indeed, electrophoretic mobility shift assay demonstrated a direct binding between the two sRNAs via GA- (CjNC110) and CU- (CjNC140) rich stem-loops. Additionally, RNAseq and follow-up experiments identified that CjNC140 positively regulates p19, which encodes a key iron uptake transporter in Campylobacter. Furthermore, computational analysis revealed both CjNC140 and CjNC110 are highly conserved in C. jejuni, and the predicted secondary structures support CjNC140 as a functional homolog of the iron regulatory sRNA, RyhB. These findings establish CjNC140 and CjNC110 as a key checks-and- balances mechanism in maintaining homeostasis of gene expression and optimizing phenotypes critical for C. jejuni pathobiology. IMPORTANCE Gene regulation is critical to all aspects of pathogenesis of bacterial disease, and small non-coding RNAs (sRNAs) represent a new frontier in gene regulation of bacteria. In Campylobacter jejuni, the role of sRNAs remains largely unexplored. Here, we investigate the role of two highly conserved sRNAs, CjNC110 and CjNC140, and demonstrate that CjNC140 displays a primarily inhibitory role in contrast to a primarily activating role for CjNC110 for several key virulence-associated phenotypes. Our results also revealed that the sRNA regulatory pathway is intertwined with the iron uptake system, another virulence mechanism critical for in vivo colonization. These findings open a new direction for understanding C. jejuni pathobiology and identify potential targets for intervention for this major foodborne pathogen.

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.

Dissecting components of the Campylobacter jejuni fetMP-fetABCDEF gene cluster in iron scavenging

Campylobacter jejuni is a leading cause of bacterial gastroenteritis worldwide. Acute infection can be antecedent to highly debilitating long-term sequelae. Expression of iron acquisition systems is vital for C. jejuni to survive the low iron availability within the human gut. The C. jejuni fetMP-fetABCDEF gene cluster is known to be upregulated during human infection and under iron limitation. While FetM and FetP have been functionally linked to iron transport in prior work, here we assess the contribution by each of the downstream genes ( fetABCDEF ) to C. jejuni growth during both iron-depleted and iron-replete conditions. Significant growth impairment was observed upon disruption of fetA , fetB, fetC , and fetD , suggesting a role in iron acquisition for each encoded protein. FetA expression was modulated by iron-availability but not dependent on the presence of FetB, FetC, FetD, FetE or FetF. Functions of the putative thioredoxins FetE and FetF were redundant in iron scavenging, requiring a double deletion (Δ fetEF ) to exhibit a growth defect. C. jejuni FetE was expressed and the structure solved to 1.50 Å, revealing structural similarity to thiol-disulfide oxidases. Functional characterization in biochemical assays showed that FetE reduced insulin at a slower rate than E. coli Trx and that together, FetEF promoted substrate oxidation in cell extracts, suggesting that FetE (and presumably FetF) are oxidoreductases that can mediate oxidation in vivo . This study advances our understanding of the contributions by the fetMP-fetABCDEF gene cluster to virulence at a genetic and functional level, providing foundational knowledge towards mitigating C. jejuni -related morbidity and mortality.

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.

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.

Multiple Campylobacter jejuni proteins affecting the peptidoglycan structure and the degree of helical cell curvature

Campylobacter jejuni is a Gram-negative helical bacterium. Its helical morphology, maintained by the peptidoglycan (PG) layer, plays a key role in its transmission in the environment, colonization, and pathogenic properties. The previously characterized PG hydrolases Pgp1 and Pgp2 are important for generating C. jejuni helical morphology, with deletion mutants being rod-shaped and showing alterations in their PG muropeptide profiles in comparison to the wild type. Homology searches and bioinformatics were used to identify additional gene products involved in C. jejuni morphogenesis: the putative bactofilin 1104 and the M23 peptidase domain-containing proteins 0166, 1105, and 1228. Deletions in the corresponding genes resulted in varying curved rod morphologies with changes in their PG muropeptide profiles. All changes in the mutants complemented except 1104. Overexpression of 1104 and 1105 also resulted in changes in the morphology and in the muropeptide profiles, suggesting that the dose of these two gene products influences these characteristics. The related helical ε-Proteobacterium Helicobacter pylori has characterized homologs of C. jejuni 1104, 1105, and 1228 proteins, yet deletion of the homologous genes in H. pylori had differing effects on H. pylori PG muropeptide profiles and/or morphology compared to the C. jejuni deletion mutants. It is therefore apparent that even related organisms with similar morphologies and homologous proteins can have diverse PG biosynthetic pathways, highlighting the importance of studying PG biosynthesis in related organisms.

TSSNote-CyaPromBERT: Development of an integrated platform for highly accurate promoter prediction and visualization of Synechococcus sp. and Synechocystis sp. through a state-of-the-art natural language processing model BERT

Since the introduction of the first transformer model with a unique self-attention mechanism, natural language processing (NLP) models have attained state-of-the-art (SOTA) performance on various tasks. As DNA is the blueprint of life, it can be viewed as an unusual language, with its characteristic lexicon and grammar. Therefore, NLP models may provide insights into the meaning of the sequential structure of DNA. In the current study, we employed and compared the performance of popular SOTA NLP models (i.e., XLNET, BERT, and a variant DNABERT trained on the human genome) to predict and analyze the promoters in freshwater cyanobacterium Synechocystis sp. PCC 6803 and the fastest growing cyanobacterium Synechococcus elongatus sp. UTEX 2973. These freshwater cyanobacteria are promising hosts for phototrophically producing value-added compounds from CO₂. Through a custom pipeline, promoters and non-promoters from Synechococcus elongatus sp. UTEX 2973 were used to train the model. The trained model achieved an AUROC score of 0.97 and F1 score of 0.92. During cross-validation with promoters from Synechocystis sp. PCC 6803, the model achieved an AUROC score of 0.96 and F1 score of 0.91. To increase accessibility, we developed an integrated platform (TSSNote-CyaPromBERT) to facilitate large dataset extraction, model training, and promoter prediction from public dRNA-seq datasets. Furthermore, various visualization tools have been incorporated to address the "black box" issue of deep learning and feature analysis. The learning transfer ability of large language models may help identify and analyze promoter regions for newly isolated strains with similar lineages.

Influence of Staphylococcus aureus Strain Background on Sa3int Phage Life Cycle Switches

Staphylococcus aureus asymptomatically colonizes the nasal cavity of mammals, but it is also a leading cause of life-threatening infections. Most human nasal isolates carry Sa3 phages, which integrate into the bacterial hlb gene encoding a sphingomyelinase. The virulence factor-encoding genes carried by the Sa3-phages are highly human-specific, and most animal strains are Sa3 negative. Thus, both insertion and excision of the prophage could potentially confer a fitness advantage to S. aureus. Here, we analyzed the phage life cycle of two Sa3 phages, Φ13 and ΦN315, in different phage-cured S. aureus strains. Based on phage transfer experiments, strains could be classified into low (8325-4, SH1000, and USA300c) and high (MW2c and Newman-c) transfer strains. High-transfer strains promoted the replication of phages, whereas phage adsorption, integration, excision, or recA transcription was not significantly different between strains. RNASeq analyses of replication-deficient lysogens revealed no strain-specific differences in the CI/Mor regulatory switch. However, lytic genes were significantly upregulated in the high transfer strain MW2c Φ13 compared to strain 8325-4 Φ13. By transcriptional start site prediction, new promoter regions within the lytic modules were identified, which are likely targeted by specific host factors. Such host-phage interaction probably accounts for the strain-specific differences in phage replication and transfer frequency. Thus, the genetic makeup of the host strains may determine the rate of phage mobilization, a feature that might impact the speed at which certain strains can achieve host adaptation.

Development and validation of a CRISPR interference system for gene regulation in Campylobacter jejuni

BACKGROUND

Campylobacter spp. are the leading cause of bacterial food-borne illness in humans worldwide, with Campylobacter jejuni responsible for 80% of these infections. There is an urgent need to understand fundamental C. jejuni biology for the development of new strategies to prevent and treat infections. The range of molecular tools available to regulate gene expression in C. jejuni is limited, which in turn constrains our ability to interrogate the function of essential and conditionally essential genes. We have addressed this by developing and utilising a CRISPR-based interference system known as CRISPRi in C. jejuni to control gene expression. To achieve this, a catalytically inactive ("dead") cas9 and sgRNA backbone from the Streptococcus pyogenes CRISPRi system was combined with C. jejuni-derived promoters of predetermined expression activities to develop a CRISPRi-based repression tool in C. jejuni strains M1Cam and 81-176.

RESULTS

The CRISPRi tool was validated through successful repression of the arylsulphatase-encoding gene astA using a range of sgRNA target sequences spanning the astA gene. The tool was also applied to target astA in an M1Cam CRISPR-Cas9 deletion strain, which showed that the presence of an endogenous CRISPR-Cas9 system did not affect the activity of the CRISPRi-based repression tool. The tool was further validated against the hippicurase-encoding gene hipO. Following this, the flagella genes flgR, flaA, flaB and both flaA and flaB were targeted for CRISPRi-based repression, which resulted in varying levels of motility reduction and flagella phenotypes as determined by phenotypical assays and transmission electron microscopy (TEM).

CONCLUSIONS

This is the first report of a CRISPRi-based tool in C. jejuni, which will provide a valuable resource to the Campylobacter community.

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.

Mapping the Complex Transcriptional Landscape of the Phytopathogenic Bacterium Dickeya dadantii

Dickeya dadantii is a phytopathogenic bacterium that causes soft rot in a wide range of plant hosts worldwide and a model organism for studying virulence gene regulation. The present study provides a comprehensive and annotated transcriptomic map of D. dadantii obtained by a computational method combining five independent transcriptomic data sets: (i) paired-end RNA sequencing (RNA-seq) data for a precise reconstruction of the RNA landscape; (ii) DNA microarray data providing transcriptional responses to a broad variety of environmental conditions; (iii) long-read Nanopore native RNA-seq data for isoform-level transcriptome validation and determination of transcription termination sites; (iv) differential RNA sequencing (dRNA-seq) data for the precise mapping of transcription start sites; (v) in planta DNA microarray data for a comparison of gene expression profiles between in vitro experiments and the early stages of plant infection. Our results show that transcription units sometimes coincide with predicted operons but are generally longer, most of them comprising internal promoters and terminators that generate alternative transcripts of variable gene composition. We characterize the occurrence of transcriptional read-through at terminators, which might play a basal regulation role and explain the extent of transcription beyond the scale of operons. We finally highlight the presence of noncontiguous operons and excludons in the D. dadantii genome, novel genomic arrangements that might contribute to the basal coordination of transcription. The highlighted transcriptional organization may allow D. dadantii to finely adjust its gene expression program for a rapid adaptation to fast-changing environments.

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.

RNA Sequencing Reveals Widespread Transcription of Natural Antisense RNAs in Entamoeba Species

Entamoeba is a genus of Amoebozoa that includes the intestine-colonizing pathogenic species Entamoeba histolytica. To understand the basis of gene regulation in E. histolytica from an evolutionary perspective, we have profiled the transcriptomes of its closely related species E. dispar, E. moshkovskii and E. invadens. Genome-wide identification of transcription start sites (TSS) and polyadenylation sites (PAS) revealed the similarities and differences of their gene regulatory sequences. In particular, we found the widespread initiation of antisense transcription from within the gene coding sequences is a common feature among all Entamoeba species. Interestingly, we observed the enrichment of antisense transcription in genes involved in several processes that are common to species infecting the human intestine, e.g., the metabolism of phospholipids. These results suggest a potentially conserved and compact gene regulatory system in Entamoeba.

Application of Bacteriophages to Limit Campylobacter in Poultry Production

Campylobacter is a major foodborne pathogen with over a million United States cases a year and is typically acquired through the consumption of poultry products. The common occurrence of Campylobacter as a member of the poultry gastrointestinal tract microbial community remains a challenge for optimizing intervention strategies. Simultaneously, increasing demand for antibiotic-free products has led to the development of several alternative control measures both at the farm and in processing operations. Bacteriophages administered to reduce foodborne pathogens are one of the alternatives that have received renewed interest. Campylobacter phages have been isolated from both conventionally and organically raised poultry. Isolated and cultivated Campylobacter bacteriophages have been used as an intervention in live birds to target colonized Campylobacter in the gastrointestinal tract. Application of Campylobacter phages to poultry carcasses has also been explored as a strategy to reduce Campylobacter levels during poultry processing. This review will focus on the biology and ecology of Campylobacter bacteriophages in poultry production followed by discussion on current and potential applications as an intervention strategy to reduce Campylobacter occurrence in poultry production.

A chromosomal loop anchor mediates bacterial genome organization

Nucleoprotein complexes play an integral role in genome organization of both eukaryotes and prokaryotes. Apart from their role in locally structuring and compacting DNA, several complexes are known to influence global organization by mediating long-range anchored chromosomal loop formation leading to spatial segregation of large sections of DNA. Such megabase-range interactions are ubiquitous in eukaryotes, but have not been demonstrated in prokaryotes. Here, using a genome-wide sedimentation-based approach, we found that a transcription factor, Rok, forms large nucleoprotein complexes in the bacterium Bacillus subtilis. Using chromosome conformation capture and live-imaging of DNA loci, we show that these complexes robustly interact with each other over large distances. Importantly, these Rok-dependent long-range interactions lead to anchored chromosomal loop formation, thereby spatially isolating large sections of DNA, as previously observed for insulator proteins in eukaryotes.

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.

The Campylobacter jejuni Response Regulator and Cyclic-Di-GMP Binding CbrR Is a Novel Regulator of Flagellar Motility

A leading cause of bacterial gastroenteritis, Campylobacter jejuni is also associated with broad sequelae, including extragastrointestinal conditions such as reactive arthritis and Guillain-Barré Syndrome (GBS). CbrR is a C. jejuni response regulator that is annotated as a diguanylate cyclase (DGC), an enzyme that catalyzes the synthesis of c-di-GMP, a universal bacterial second messenger, from GTP. In C. jejuni DRH212, we constructed an unmarked deletion mutant, cbrR^-, and complemented mutant, cbrR⁺. Motility assays indicated a hyper-motile phenotype associated with cbrR^-, whereas motility was deficient in cbrR⁺. The overexpression of CbrR in cbrR⁺ was accompanied by a reduction in expression of FlaA, the major flagellin. Biofilm assays and scanning electron microscopy demonstrated similarities between DRH212 and cbrR^-; however, cbrR⁺ was unable to form significant biofilms. Transmission electron microscopy showed similar cell morphology between the three strains; however, cbrR⁺ cells lacked flagella. Differential radial capillary action of ligand assays (DRaCALA) showed that CbrR binds GTP and c-di-GMP. Liquid chromatography tandem mass spectrometry detected low levels of c-di-GMP in C. jejuni and in E. coli expressing CbrR. CbrR is therefore a negative regulator of FlaA expression and motility, a critical virulence factor in C. jejuni pathogenesis.

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.

The tracrRNA in CRISPR Biology and Technologies

CRISPR-Cas adaptive immune systems in bacteria and archaea utilize short CRISPR RNAs (crRNAs) to guide sequence-specific recognition and clearance of foreign genetic material. Multiple crRNAs are stored together in a compact format called a CRISPR array that is transcribed and processed into the individual crRNAs. While the exact processing mechanisms vary widely, some CRISPR-Cas systems, including those encoding the Cas9 nuclease, rely on a trans-activating crRNA (tracrRNA). The tracrRNA was discovered in 2011 and was quickly co-opted to create single-guide RNAs as core components of CRISPR-Cas9 technologies. Since then, further studies have uncovered processes extending beyond the traditional role of tracrRNA in crRNA biogenesis, revealed Cas nucleases besides Cas9 that are dependent on tracrRNAs, and established new applications based on tracrRNA engineering. In this review, we describe the biology of the tracrRNA and how its ongoing characterization has garnered new insights into prokaryotic immune defense and enabled key technological advances.

Promotech: a general tool for bacterial promoter recognition

Promoters are genomic regions where the transcription machinery binds to initiate the transcription of specific genes. Computational tools for identifying bacterial promoters have been around for decades. However, most of these tools were designed to recognize promoters in one or few bacterial species. Here, we present Promotech, a machine-learning-based method for promoter recognition in a wide range of bacterial species. We compare Promotech's performance with the performance of five other promoter prediction methods. Promotech outperforms these other programs in terms of area under the precision-recall curve (AUPRC) or precision at the same level of recall. Promotech is available at https://github.com/BioinformaticsLabAtMUN/PromoTech .

Digging into the lesser-known aspects of CRISPR biology

A long time has passed since regularly interspaced DNA repeats were discovered in prokaryotes. Today, those enigmatic repetitive elements termed clustered regularly interspaced short palindromic repeats (CRISPR) are acknowledged as an emblematic part of multicomponent CRISPR-Cas (CRISPR associated) systems. These systems are involved in a variety of roles in bacteria and archaea, notably, that of conferring protection against transmissible genetic elements through an adaptive immune-like response. This review summarises the present knowledge on the diversity, molecular mechanisms and biology of CRISPR-Cas. We pay special attention to the most recent findings related to the determinants and consequences of CRISPR-Cas activity. Research on the basic features of these systems illustrates how instrumental the study of prokaryotes is for understanding biology in general, ultimately providing valuable tools for diverse fields and fuelling research beyond the mainstream.

Comparative genomics provides structural and functional insights into Bacteroides RNA biology

Bacteria employ noncoding RNA molecules for a wide range of biological processes, including scaffolding large molecular complexes, catalyzing chemical reactions, defending against phages, and controlling gene expression. Secondary structures, binding partners, and molecular mechanisms have been determined for numerous small noncoding RNAs (sRNAs) in model aerobic bacteria. However, technical hurdles have largely prevented analogous analyses in the anaerobic gut microbiota. While experimental techniques are being developed to investigate the sRNAs of gut commensals, computational tools and comparative genomics can provide immediate functional insight. Here, using Bacteroides thetaiotaomicron as a representative microbiota member, we illustrate how comparative genomics improves our understanding of the RNA biology in an understudied gut bacterium. We investigate putative RNA-binding proteins and predict a Bacteroides cold-shock protein homologue to have an RNA-related function. We apply an in-silico protocol incorporating both sequence and structural analysis to determine the consensus structures and conservation of nine Bacteroides noncoding RNA families. Using structure probing, we validate and refine these predictions, and deposit them in the Rfam database. Through synteny analyses, we illustrate how genomic co-conservation can serve as a predictor of sRNA function. Altogether, this work showcases the power of RNA informatics for investigating the RNA biology of anaerobic microbiota members.

RNA landscape of the emerging cancer-associated microbe Fusobacterium nucleatum

Fusobacterium nucleatum, long known as a constituent of the oral microflora, has recently garnered renewed attention for its association with several different human cancers. The growing interest in this emerging cancer-associated bacterium contrasts with a paucity of knowledge about its basic gene expression features and physiological responses. As fusobacteria lack all established small RNA-associated proteins, post-transcriptional networks in these bacteria are also unknown. In the present study, using differential RNA-sequencing, we generate high-resolution global RNA maps for five clinically relevant fusobacterial strains-F. nucleatum subspecies nucleatum, animalis, polymorphum and vincentii, as well as F. periodonticum-for early, mid-exponential growth and early stationary phase. These data are made available in an online browser, and we use these to uncover fundamental aspects of fusobacterial gene expression architecture and a suite of non-coding RNAs. Developing a vector for functional analysis of fusobacterial genes, we discover a conserved fusobacterial oxygen-induced small RNA, FoxI, which serves as a post-transcriptional repressor of the major outer membrane porin FomA. Our findings provide a crucial step towards delineating the regulatory networks enabling F. nucleatum adaptation to different environments, which may elucidate how these bacteria colonize different compartments of the human body.

Widespread divergent transcription from bacterial and archaeal promoters is a consequence of DNA-sequence symmetry

Transcription initiates at promoters, DNA regions recognized by a DNA-dependent RNA polymerase. We previously identified horizontally acquired Escherichia coli promoters from which the direction of transcription was unclear. In the present study, we show that more than half of these promoters are bidirectional and drive divergent transcription. Using genome-scale approaches, we demonstrate that 19% of all transcription start sites detected in E. coli are associated with a bidirectional promoter. Bidirectional promoters are similarly common in diverse bacteria and archaea, and have inherent symmetry: specific bases required for transcription initiation are reciprocally co-located on opposite DNA strands. Bidirectional promoters enable co-regulation of divergent genes and are enriched in both intergenic and horizontally acquired regions. Divergent transcription is conserved among bacteria, archaea and eukaryotes, but the underlying mechanisms for bidirectionality are different.

Noncanonical crRNAs derived from host transcripts enable multiplexable RNA detection by Cas9

CRISPR-Cas systems recognize foreign genetic material using CRISPR RNAs (crRNAs). In type II systems, a trans-activating crRNA (tracrRNA) hybridizes to crRNAs to drive their processing and utilization by Cas9. While analyzing Cas9-RNA complexes from Campylobacter jejuni, we discovered tracrRNA hybridizing to cellular RNAs, leading to formation of "noncanonical" crRNAs capable of guiding DNA targeting by Cas9. Our discovery inspired the engineering of reprogrammed tracrRNAs that link the presence of any RNA of interest to DNA targeting with different Cas9 orthologs. This capability became the basis for a multiplexable diagnostic platform termed LEOPARD (leveraging engineered tracrRNAs and on-target DNAs for parallel RNA detection). LEOPARD allowed simultaneous detection of RNAs from different viruses in one test and distinguished severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its D614G (Asp614→Gly) variant with single-base resolution in patient samples.

Inactivation of the core cheVAWY chemotaxis genes disrupts chemotactic motility and organised biofilm formation in Campylobacter jejuni

Flagellar motility plays a central role in the bacterial foodborne pathogen Campylobacter jejuni, as flagellar motility is required for reaching the intestinal epithelium and subsequent colonisation or disease. Flagellar proteins also contribute strongly to biofilm formation during transmission. Chemotaxis is the process directing flagellar motility in response to attractant and repellent stimuli, but its role in biofilm formation of C. jejuni is not well understood. Here we show that inactivation of the core chemotaxis genes cheVAWY in C. jejuni strain NCTC 11168 affects both chemotactic motility and biofilm formation. Inactivation of any of the core chemotaxis genes (cheA, cheY, cheV or cheW) impaired chemotactic motility but did not affect flagellar assembly or growth. The ∆cheY mutant swam in clockwise loops, while complementation restored normal motility. Inactivation of the core chemotaxis genes interfered with the ability to form a discrete biofilm at the air-media interface, and the ∆cheY mutant displayed reduced dispersal/shedding of bacteria into the planktonic fraction. This suggests that while the chemotaxis system is not required for biofilm formation per se, it is necessary for organized biofilm formation. Hence interference with the Campylobacter chemotaxis system at any level disrupts optimal chemotactic motility and transmission modes such as biofilm formation.

Survival Strategies of Streptococcus pyogenes in Response to Phage Infection

Bacteriophages exert strong evolutionary pressure on their microbial hosts. In their lytic lifecycle, complete bacterial subpopulations are utilized as hosts for bacteriophage replication. However, during their lysogenic lifecycle, bacteriophages can integrate into the host chromosome and alter the host’s genomic make-up, possibly resulting in evolutionary important adjustments. Not surprisingly, bacteria have evolved sophisticated immune systems to protect against phage infection. Streptococcus pyogenes isolates are frequently lysogenic and their prophages have been shown to be major contributors to the virulence of this pathogen. Most S. pyogenes phage research has focused on genomic prophages in relation to virulence, but little is known about the defensive arsenal of S. pyogenes against lytic phage infection. Here, we characterized Phage A1, an S. pyogenes bacteriophage, and investigated several mechanisms that S. pyogenes utilizes to protect itself against phage predation. We show that Phage A1 belongs to the Siphoviridae family and contains a circular double-stranded DNA genome that follows a modular organization described for other streptococcal phages. After infection, the Phage A1 genome can be detected in isolated S. pyogenes survivor strains, which enables the survival of the bacterial host and Phage A1 resistance. Furthermore, we demonstrate that the type II-A CRISPR-Cas system of S. pyogenes acquires new spacers upon phage infection, which are increasingly detectable in the absence of a capsule. Lastly, we show that S. pyogenes produces membrane vesicles that bind to phages, thereby limiting the pool of phages available for infection. Altogether, this work provides novel insight into survival strategies employed by S. pyogenes to combat phage predation.

Identification of a non-coding RNA and its putative involvement in the regulation of tetanus toxin synthesis in Clostridium tetani

Clostridium tetani produces the tetanus toxin (TeNT), one of the most powerful bacterial toxins known to humankind and responsible for tetanus. The regulation of toxin expression is complex and involves the alternative sigma factor TetR as well as other regulators. Here, a transcriptional analysis of the TeNT-encoding large plasmid of C. tetani identified a putative non-coding small RNA (sRNA), located in close vicinity of the 3′ untranslated region of the tent gene. A northern blot experiment could identify a respective sRNA with a size of approx. 140 nucleotides. Sequence analysis showed that the sRNA contains a 14-nucleotide region that is complementary to a 5′ located region of tent. In order to investigate the function of the sRNA, we applied a RNA interference approach targeting the sRNA in two C. tetani wild-type strains; the constructed antisense C. tetani strains showed an approx. threefold increase in both extracellular and total TeNT production compared to the respective wild-type strains. In addition, recombinant C. tetani strains were constructed that contained tent-locus harboring plasmids with and without the sRNA. However, the introduction of the tent-locus without the sRNA in a C. tetani strain lacking the wild-type TeNT-encoding large plasmid resulted in a lower TeNT production compared to the same strain with recombinant tent-locus with the sRNA. This suggests that the expression or the effect of the sRNA is modulated by the C. tetani genetic background, notably that of the wild-type TeNT-encoding large plasmid. In addition, some recombinant strains exhibited modulated growth patterns, characterized by premature bacterial cell lysis. Taken together, our data indicate that the sRNA acts as a negative regulator of TeNT synthesis, with a possible impact on the growth of C. tetani. We hypothesize that the role of this sRNA is to limit toxin levels in the exponential growth phase in order to prevent premature bacterial lysis.

Cross-species RNA-seq for deciphering host-microbe interactions

The human body is constantly exposed to microorganisms, which entails manifold interactions between human cells and diverse commensal or pathogenic bacteria. The cellular states of the interacting cells are decisive for the outcome of these encounters such as whether bacterial virulence programmes and host defence or tolerance mechanisms are induced. This Review summarizes how next-generation RNA sequencing (RNA-seq) has become a primary technology to study host-microbe interactions with high resolution, improving our understanding of the physiological consequences and the mechanisms at play. We illustrate how the discriminatory power and sensitivity of RNA-seq helps to dissect increasingly complex cellular interactions in time and space down to the single-cell level. We also outline how future transcriptomics may answer currently open questions in host-microbe interactions and inform treatment schemes for microbial disorders.

Coupled Transcriptomics for Differential Expression Analysis and Determination of Transcription Start Sites: Design and Bioinformatics

The term coupled transcriptomics is coined to describe a design of an RNA-seq experiment intended for both differential expression analysis and genome-wide determination of the transcription start sites (TSS). The minimal requirements for the first analysis are two experimental conditions with at least two biological replicates enabling statistical tests. The second analysis involves the bioinformatics comparison of the data generated from a control RNA-seq library with another library enriched in primary transcripts using Terminator^™ 5'-phosphate-dependent exonuclease, in an experiment denominated differential RNA-seq (dRNA-seq). Usually, dRNA-seq is carried out with specific protocols for library construction, different of those used for common differential expression analysis. Our experimental design allows to use the same data for both analyses, reducing the number of libraries to be generated and sequenced. This is a guide for designing a coupled transcriptomics experiment and for the subsequent bioinformatics procedures. The proposed methods can be applied to the detection and study of small RNA genes.

Natural Competence and Horizontal Gene Transfer in Campylobacter

Thermophilic Campylobacter, in particular Campylobacter jejuni, C. coli and C. lari are the main relevant Campylobacter species for human infections. Due to their high capacity of genetic exchange by horizontal gene transfer (HGT), rapid adaptation to changing environmental and host conditions contribute to successful spreading and persistence of these foodborne pathogens. However, extensive HGT can exert dangerous side effects for the bacterium, such as the incorporation of gene fragments leading to disturbed gene functions. Here we discuss mechanisms of HGT, notably natural transformation, conjugation and bacteriophage transduction and limiting regulatory strategies of gene transfer. In particular, we summarize the current knowledge on how the DNA macromolecule is exchanged between single cells. Mechanisms to stimulate and to limit HGT obviously coevolved and maintained an optimal balance. Chromosomal rearrangements and incorporation of harmful mutations are risk factors for survival and can result in drastic loss of fitness. In Campylobacter, the restricted recognition and preferential uptake of free DNA from relatives are mediated by a short methylated DNA pattern and not by a classical DNA uptake sequence as found in other bacteria. A class two CRISPR-Cas system is present but also other DNases and restriction-modification systems appear to be important for Campylobacter genome integrity. Several lytic and integrated bacteriophages have been identified, which contribute to genome diversity. Furthermore, we focus on the impact of gene transfer on the spread of antibiotic resistance genes (resistome) and persistence factors. We discuss remaining open questions in the HGT field, supposed to be answered in the future by current technologies like whole-genome sequencing and single-cell approaches.

Comparative Methylome Analysis of Campylobacter jejuni Strain YH002 Reveals a Putative Novel Motif and Diverse Epigenetic Regulations of Virulence Genes

Campylobacter jejuni is a major cause of foodborne gastroenteritis worldwide inflicting palpable socioeconomic costs. The ability of this pathogen to successfully infect its hosts is determined not only by the presence of specific virulence genes but also by the pathogen’s capacity to appropriately regulate those virulence genes. Therefore, DNA methylation can play a critical role in both aspects of this process because it serves as both a means to protect the integrity of the cellular DNA from invasion and as a mechanism to control transcriptional regulation within the cell. In the present study we report the comparative methylome data of C. jejuni YH002, a multidrug resistant strain isolated from retail beef liver. Investigation into the methylome identified a putative novel motif (CGCGA) of a type II restriction-modification (RM) system. Comparison of methylomes of the strain to well-studied C. jejuni strains highlighted non-uniform methylation patterns among the strains though the existence of the typical type I and type IV RM systems were also observed. Additional investigations into the existence of DNA methylation sites within gene promoters, which may ultimately result in altered levels of transcription, revealed several virulence genes putatively regulated using this mode of action. Of those identified, a flagella gene (flhB), a RNA polymerase sigma factor (rpoN), a capsular polysaccharide export protein (kpsD), and a multidrug efflux pump were highly notable.

Single-Nucleotide RNA Maps for the Two Major Nosocomial Pathogens Enterococcus faecalis and Enterococcus faecium

Enterococcus faecalis and faecium are two major representative clinical strains of the Enterococcus genus and are sadly notorious to be part of the top agents responsible for nosocomial infections. Despite their critical implication in worldwide public healthcare, essential and available resources such as deep transcriptome annotations remain poor, which also limits our understanding of post-transcriptional control small regulatory RNA (sRNA) functions in these bacteria. Here, using the dRNA-seq technique in combination with ANNOgesic analysis, we successfully mapped and annotated transcription start sites (TSS) of both E. faecalis V583 and E. faecium AUS0004 at single nucleotide resolution. Analyzing bacteria in late exponential phase, we capture ~40% (E. faecalis) and 43% (E. faecium) of the annotated protein-coding genes, determine 5′ and 3′ UTR (untranslated region) length, and detect instances of leaderless mRNAs. The transcriptome maps revealed sRNA candidates in both bacteria, some found in previous studies and new ones. Expression of candidate sRNAs is being confirmed under biologically relevant environmental conditions. This comprehensive global TSS mapping atlas provides a valuable resource for RNA biology and gene expression analysis in the Enterococci. It can be accessed online at www.helmholtz-hiri.de/en/datasets/enterococcus through an instance of the genomic viewer JBrowse.

Genomic determinants for initiation and length of natural antisense transcripts in Entamoeba histolytica

Natural antisense transcripts (NAT) have been reported in prokaryotes and eukaryotes. While the functions of most reported NATs remain unknown, their potentials in regulating the transcription of their counterparts have been speculated. Entamoeba histolytica, which is a unicellular eukaryotic parasite, has a compact protein-coding genome with very short intronic and intergenic regions. The regulatory mechanisms of gene expression in this compact genome are under-described. In this study, by genome-wide mapping of RNA-Seq data in the genome of E. histolytica, we show that a substantial fraction of its protein-coding genes (28%) has significant transcription on their opposite strand (i.e. NAT). Intriguingly, we found the location of transcription start sites or polyadenylation sites of NAT are determined by the specific motifs encoded on the opposite strand of the gene coding sequences, thereby providing a compact regulatory system for gene transcription. Moreover, we demonstrated that NATs are globally up-regulated under various environmental conditions including temperature stress and pathogenicity. While NATs do not appear to be consequences of spurious transcription, they may play a role in regulating gene expression in E. histolytica, a hypothesis which needs to be tested.

Transcription termination and antitermination of bacterial CRISPR arrays

A hallmark of CRISPR-Cas immunity systems is the CRISPR array, a genomic locus consisting of short, repeated sequences ('repeats') interspersed with short, variable sequences ('spacers'). CRISPR arrays are transcribed and processed into individual CRISPR RNAs that each include a single spacer, and direct Cas proteins to complementary sequences in invading nucleic acid. Most bacterial CRISPR array transcripts are unusually long for untranslated RNA, suggesting the existence of mechanisms to prevent premature transcription termination by Rho, a conserved bacterial transcription termination factor that rapidly terminates untranslated RNA. We show that Rho can prematurely terminate transcription of bacterial CRISPR arrays, and we identify a widespread antitermination mechanism that antagonizes Rho to facilitate complete transcription of CRISPR arrays. Thus, our data highlight the importance of transcription termination and antitermination in the evolution of bacterial CRISPR-Cas systems.

RNA-based thermoregulation of a Campylobacter jejuni zinc resistance determinant

RNA thermometers (RNATs) trigger bacterial virulence factor expression in response to the temperature shift on entering a warm-blooded host. At lower temperatures these secondary structures sequester ribosome-binding sites (RBSs) to prevent translation initiation, whereas at elevated temperatures they "melt" allowing translation. Campylobacter jejuni is the leading bacterial cause of human gastroenteritis worldwide yet little is known about how it interacts with the host including host induced gene regulation. Here we demonstrate that an RNAT regulates a C. jejuni gene, Cj1163c or czcD, encoding a member of the Cation Diffusion Facilitator family. The czcD upstream untranslated region contains a predicted stem loop within the mRNA that sequesters the RBS to inhibit translation at temperatures below 37°C. Mutations that disrupt or enhance predicted secondary structure have significant and predictable effects on temperature regulation. We also show that in an RNAT independent manner, CzcD expression is induced by Zn(II). Mutants lacking czcD are hypersensitive to Zn(II) and also over-accumulate Zn(II) relative to wild-type, all consistent with CzcD functioning as a Zn(II) exporter. Importantly, we demonstrate that C. jejuni Zn(II)-tolerance at 32°C, a temperature at which the RNAT limits CzcD production, is increased by RNAT disruption. Finally we show that czcD inactivation attenuates larval killing in a Galleria infection model and that at 32°C disrupting RNAT secondary structure to allow CzcD production can enhance killing. We hypothesise that CzcD regulation by metals and temperature provides a mechanism for C. jejuni to overcome innate immune system-mediated Zn(II) toxicity in warm-blooded animal hosts.

Interplay between Regulatory RNAs and Signal Transduction Systems during Bacterial Infection

The ability of pathogenic bacteria to stably infect the host depends on their capacity to respond and adapt to the host environment and on the efficiency of their defensive mechanisms. Bacterial envelope provides a physical barrier protecting against environmental threats. It also constitutes an important sensory interface where numerous sensing systems are located. Signal transduction systems include Two-Component Systems (TCSs) and alternative sigma factors. These systems are able to sense and respond to the ever-changing environment inside the host, altering the bacterial transcriptome to mitigate the impact of the stress. The regulatory networks associated with signal transduction systems comprise small regulatory RNAs (sRNAs) that can be directly involved in the expression of virulence factors. The aim of this review is to describe the importance of TCS- and alternative sigma factor-associated sRNAs in human pathogens during infection. The currently available genome-wide approaches for studies of TCS-regulated sRNAs will be discussed. The differences in the signal transduction mediated by TCSs between bacteria and higher eukaryotes and the specificity of regulatory RNAs for their targets make them appealing targets for discovery of new strategies to fight against multi-resistant bacteria.

Native Processing of Single Guide RNA Transcripts to Create Catalytic Cas9/Single Guide RNA Complexes in Planta

The present CRISPR/Cas9 gene editing dogma for single guide RNA (sgRNA) delivery is based on the premise that 5'-and 3'-nucleotide overhangs negate Cas9/sgRNA catalytic activity in vivo. This has led to engineering strategies designed to either avoid or remove extraneous nucleotides at the 5' and 3' termini of sgRNAs. Previously, we used a Tobacco mosaic virus viral vector to express both GFP and a sgRNA from a single virus-derived mRNA in Nicotiana benthamiana This vector yielded high levels of GFP and catalytically active sgRNAs. Here, in an effort to understand the biochemical interactions of this result, we used in vitro assays to demonstrate that nucleotide overhangs 5', but not 3', proximal to the sgRNA do in fact inactivate Cas9 catalytic activity at the specified target site. Next we showed that in planta sgRNAs bound to Cas9 are devoid of the expected 5' overhangs transcribed by the virus. Furthermore, when a plant nuclear promoter was used for expression of the GFP-sgRNA fusion transcript, it also produced indels when delivered with Cas9. These results reveal that 5' auto-processing of progenitor sgRNAs occurs natively in plants. Toward a possible mechanism for the perceived auto-processing, we found, using in vitro-generated RNAs and those isolated from plants, that the 5' to 3' exoribonuclease XRN1 can degrade elongated progenitor sgRNAs, whereas the mature sgRNA end products are resistant. Comparisons with other studies suggest that sgRNA auto-processing may be a phenomenon not unique to plants, but present in other eukaryotes as well.

An RNA-centric view on gut Bacteroidetes

Bacteria employ noncoding RNAs to maintain cellular physiology, adapt global gene expression to fluctuating environments, sense nutrients, coordinate their interaction with companion microbes and host cells, and protect themselves against bacteriophages. While bacterial RNA research has made fundamental contributions to biomedicine and biotechnology, the bulk of our knowledge of RNA biology stems from the study of a handful of aerobic model species. In comparison, RNA research is lagging in many medically relevant obligate anaerobic species, in particular the numerous commensal bacteria comprising our gut microbiota. This review presents a guide to RNA-based regulatory mechanisms in the phylum Bacteroidetes, focusing on the most abundant bacterial genus in the human gut, Bacteroides spp. This includes recent case reports on riboswitches, an mRNA leader, cis- and trans-encoded small RNAs (sRNAs) in Bacteroides spp., and a survey of CRISPR-Cas systems across Bacteroidetes. Recent work from our laboratory now suggests the existence of hundreds of noncoding RNA candidates in Bacteroides thetaiotaomicron, the emerging model organism for functional microbiota research. Based on these collective observations, we predict mechanistic and functional commonalities and differences between Bacteroides sRNAs and those of other model bacteria, and outline open questions and tools needed to boost Bacteroidetes RNA research.

Early termination of the Shiga toxin transcript generates a regulatory small RNA

Enterohemorrhagic E. coli is a significant human pathogen that can cause severe disease due to the release of Shiga toxins. The toxins are encoded within lysogenic bacteriophage and controlled by antitermination of the phage late promoter, P_R′. This promoter is always active, but terminated immediately downstream during lysogeny. A byproduct of antitermination regulation is transcription of a short RNA that is thought to be nonfunctional. Here we demonstrate that in Shiga toxin-encoding phages, this short RNA is a Hfq-binding regulatory small RNA. The small RNA represses toxin production threefold under lysogenic conditions and promotes high cell density growth. Lysogenic bacteriophages are highly abundant and our results suggest that antiterminated phage promoters may be a rich source of regulatory RNAs.

Enterohemorrhagic Escherichia coli is a significant human pathogen that causes disease ranging from hemorrhagic colitis to hemolytic uremic syndrome. The latter can lead to potentially fatal renal failure and is caused by the release of Shiga toxins that are encoded within lambdoid bacteriophages. The toxins are encoded within the late transcript of the phage and are regulated by antitermination of the P_R′ late promoter during lytic induction of the phage. During lysogeny, the late transcript is prematurely terminated at t_R′ immediately downstream of P_R′, generating a short RNA that is a byproduct of antitermination regulation. We demonstrate that this short transcript binds the small RNA chaperone Hfq, and is processed into a stable 74-nt regulatory small RNA that we have termed StxS. StxS represses expression of Shiga toxin 1 under lysogenic conditions through direct interactions with the stx1AB transcript. StxS acts in trans to activate expression of the general stress response sigma factor, RpoS, through direct interactions with an activating seed sequence within the 5′ UTR. Activation of RpoS promotes high cell density growth under nutrient-limiting conditions. Many phages utilize antitermination to regulate the lytic/lysogenic switch and our results demonstrate that short RNAs generated as a byproduct of this regulation can acquire regulatory small RNA functions that modulate host fitness.

Whole-genome Identification of Transcriptional Start Sites by Differential RNA-seq in Bacteria

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

The CARF Protein MM_0565 Affects Transcription of the Casposon-Encoded cas1-solo Gene in Methanosarcina mazei Gö1

Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) loci are found in bacterial and archaeal genomes where they provide the molecular machinery for acquisition of immunity against foreign DNA. In addition to the cas genes fundamentally required for CRISPR activity, a second class of genes is associated with the CRISPR loci, of which many have no reported function in CRISPR-mediated immunity. Here, we characterize MM_0565 associated to the type I-B CRISPR-locus of Methanosarcina mazei Gö1. We show that purified MM_0565 composed of a CRISPR-Cas Associated Rossmann Fold (CARF) and a winged helix-turn-helix domain forms a dimer in solution; in vivo, the dimeric MM_0565 is strongly stabilized under high salt stress. While direct effects on CRISPR-Cas transcription were not detected by genetic approaches, specific binding of MM_0565 to the leader region of both CRISPR-Cas systems was observed by microscale thermophoresis and electromobility shift assays. Moreover, overexpression of MM_0565 strongly induced transcription of the cas1-solo gene located in the recently reported casposon, the gene product of which shows high similarity to classical Cas1 proteins. Based on our findings, and taking the absence of the expressed CRISPR locus-encoded Cas1 protein into account, we hypothesize that MM_0565 might modulate the activity of the CRISPR systems on different levels.