Systematic interrogation of CRISPR antimicrobials in Klebsiella pneumoniae reveals nuclease-, guide- and strain-dependent features influencing antimicrobial activity

CRISPR-Cas systems can be utilized as programmable-spectrum antimicrobials to combat bacterial infections. However, how CRISPR nucleases perform as antimicrobials across target sites and strains remains poorly explored. Here, we address this knowledge gap by systematically interrogating the use of CRISPR antimicrobials using multidrug-resistant and hypervirulent strains of Klebsiella pneumoniae as models. Comparing different Cas nucleases, DNA-targeting nucleases outperformed RNA-targeting nucleases based on the tested targets. Focusing on AsCas12a that exhibited robust targeting across different strains, we found that the elucidated modes of escape varied widely, restraining opportunities to enhance killing. We also encountered individual guide RNAs yielding different extents of clearance across strains, which were linked to an interplay between improper gRNA folding and strain-specific DNA repair and survival. To explore features that could improve targeting across strains, we performed a genome-wide screen in different K. pneumoniae strains that yielded guide design rules and trained an algorithm for predicting guide efficiency. Finally, we showed that Cas12a antimicrobials can be exploited to eliminate K. pneumoniae when encoded in phagemids delivered by T7-like phages. Altogether, our results highlight the importance of evaluating antimicrobial activity of CRISPR antimicrobials across relevant strains and define critical parameters for efficient CRISPR-based targeting.

Expanding the flexibility of base editing for high-throughput genetic screens in bacteria

Genome-wide screens have become powerful tools for elucidating genotype-to-phenotype relationships in bacteria. Of the varying techniques to achieve knockout and knockdown, CRISPR base editors are emerging as promising options. However, the limited number of available, efficient target sites hampers their use for high-throughput screening. Here, we make multiple advances to enable flexible base editing as part of high-throughput genetic screening in bacteria. We first co-opt the Streptococcus canis Cas9 that exhibits more flexible protospacer-adjacent motif recognition than the traditional Streptococcus pyogenes Cas9. We then expand beyond introducing premature stop codons by mutating start codons. Next, we derive guide design rules by applying machine learning to an essentiality screen conducted in Escherichia coli. Finally, we rescue poorly edited sites by combining base editing with Cas9-induced cleavage of unedited cells, thereby enriching for intended edits. The efficiency of this dual system was validated through a conditional essentiality screen based on growth in minimal media. Overall, expanding the scope of genome-wide knockout screens with base editors could further facilitate the investigation of new gene functions and interactions in bacteria.

A method to correct for local alterations in DNA copy number that bias functional genomics assays applied to antibiotic-treated bacteria

Functional genomics techniques, such as transposon insertion sequencing and RNA-sequencing, are key to studying relative differences in bacterial mutant fitness or gene expression under selective conditions. However, certain stress conditions, mutations, or antibiotics can directly interfere with DNA synthesis, resulting in systematic changes in local DNA copy numbers along the chromosome. This can lead to artifacts in sequencing-based functional genomics data when comparing antibiotic treatment to an unstressed control. Further, relative differences in gene-wise read counts may result from alterations in chromosomal replication dynamics, rather than selection or direct gene regulation. We term this artifact "chromosomal location bias" and implement a principled statistical approach to correct it by calculating local normalization factors along the chromosome. These normalization factors are then directly incorporated into statistical analyses using standard RNA-sequencing analysis methods without modifying the read counts themselves, preserving important information about the mean-variance relationship in the data. We illustrate the utility of this approach by generating and analyzing a ciprofloxacin-treated transposon insertion sequencing data set in Escherichia coli as a case study. We show that ciprofloxacin treatment generates chromosomal location bias in the resulting data, and we further demonstrate that failing to correct for this bias leads to false predictions of mutant drug sensitivity as measured by minimum inhibitory concentrations. We have developed an R package and user-friendly graphical Shiny application, ChromoCorrect, that detects and corrects for chromosomal bias in read count data, enabling the application of functional genomics technologies to the study of antibiotic stress.IMPORTANCEAltered gene dosage due to changes in DNA replication has been observed under a variety of stresses with a variety of experimental techniques. However, the implications of changes in gene dosage for sequencing-based functional genomics assays are rarely considered. We present a statistically principled approach to correcting for the effect of changes in gene dosage, enabling testing for differences in the fitness effects or regulation of individual genes in the presence of confounding differences in DNA copy number. We show that failing to correct for these effects can lead to incorrect predictions of resistance phenotype when applying functional genomics assays to investigate antibiotic stress, and we provide a user-friendly application to detect and correct for changes in DNA copy number.

Improved RNA stability estimation through Bayesian modeling reveals most Salmonella transcripts have subminute half-lives

RNA decay is a crucial mechanism for regulating gene expression in response to environmental stresses. In bacteria, RNA-binding proteins (RBPs) are known to be involved in posttranscriptional regulation, but their global impact on RNA half-lives has not been extensively studied. To shed light on the role of the major RBPs ProQ and CspC/E in maintaining RNA stability, we performed RNA sequencing of Salmonella enterica over a time course following treatment with the transcription initiation inhibitor rifampicin (RIF-seq) in the presence and absence of these RBPs. We developed a hierarchical Bayesian model that corrects for confounding factors in rifampicin RNA stability assays and enables us to identify differentially decaying transcripts transcriptome-wide. Our analysis revealed that the median RNA half-life in Salmonella in early stationary phase is less than 1 min, a third of previous estimates. We found that over half of the 500 most long-lived transcripts are bound by at least one major RBP, suggesting a general role for RBPs in shaping the transcriptome. Integrating differential stability estimates with cross-linking and immunoprecipitation followed by RNA sequencing (CLIP-seq) revealed that approximately 30% of transcripts with ProQ binding sites and more than 40% with CspC/E binding sites in coding or 3' untranslated regions decay differentially in the absence of the respective RBP. Analysis of differentially destabilized transcripts identified a role for ProQ in the oxidative stress response. Our findings provide insights into posttranscriptional regulation by ProQ and CspC/E, and the importance of RBPs in regulating gene expression.

An expanded transcriptome atlas for Bacteroides thetaiotaomicron reveals a small RNA that modulates tetracycline sensitivity

Plasticity in gene expression allows bacteria to adapt to diverse environments. This is particularly relevant in the dynamic niche of the human intestinal tract; however, transcriptional networks remain largely unknown for gut-resident bacteria. Here we apply differential RNA sequencing (RNA-seq) and conventional RNA-seq to the model gut bacterium Bacteroides thetaiotaomicron to map transcriptional units and profile their expression levels across 15 in vivo-relevant growth conditions. We infer stress- and carbon source-specific transcriptional regulons and expand the annotation of small RNAs (sRNAs). Integrating this expression atlas with published transposon mutant fitness data, we predict conditionally important sRNAs. These include MasB, which downregulates tetracycline tolerance. Using MS2 affinity purification and RNA-seq, we identify a putative MasB target and assess its role in the context of the MasB-associated phenotype. These data-publicly available through the Theta-Base web browser ( http://micromix.helmholtz-hiri.de/bacteroides/ )-constitute a valuable resource for the microbiome community.

A comparative analysis of peptide-delivered antisense antibiotics employing diverse nucleotide mimics

Antisense oligomer (ASO)-based antibiotics that target mRNAs of essential bacterial genes have great potential for counteracting antimicrobial resistance and for precision microbiome editing. To date, the development of such antisense antibiotics has primarily focused on using phosphorodiamidate morpholino (PMO) and peptide nucleic acid (PNA) backbones, largely ignoring the growing number of chemical modalities that have spurred the success of ASO-based human therapy. Here, we directly compare the activities of seven chemically distinct 10mer ASOs, all designed to target the essential gene acpP upon delivery with a KFF-peptide carrier into Salmonella. Our systematic analysis of PNA, PMO, phosphorothioate-modified DNA (PTO), 2'-methylated RNA (RNA-OMe), 2'-methoxyethylated RNA (RNA-MOE), 2'-fluorinated RNA (RNA-F) and 2'-4'-locked RNA (LNA) is based on a variety of in vitro and in vivo methods to evaluate ASO uptake, target pairing and inhibition of bacterial growth. Our data show that only PNA and PMO are efficiently delivered by the KFF peptide into Salmonella to inhibit bacterial growth. Nevertheless, the strong target binding affinity and in vitro translational repression activity of LNA and RNA-MOE make them promising modalities for antisense antibiotics that will require the identification of an effective carrier.

Improved prediction of bacterial CRISPRi guide efficiency from depletion screens through mixed-effect machine learning and data integration

CRISPR interference (CRISPRi) is the leading technique to silence gene expression in bacteria; however, design rules remain poorly defined. We develop a best-in-class prediction algorithm for guide silencing efficiency by systematically investigating factors influencing guide depletion in genome-wide essentiality screens, with the surprising discovery that gene-specific features substantially impact prediction. We develop a mixed-effect random forest regression model that provides better estimates of guide efficiency. We further apply methods from explainable AI to extract interpretable design rules from the model. This study provides a blueprint for predictive models for CRISPR technologies where only indirect measurements of guide activity are available.

Ribosome profiling reveals the fine-tuned response of Escherichia coli to mild and severe acid stress

IMPORTANCE

Bacteria react very differently to survive in acidic environments, such as the human gastrointestinal tract. Escherichia coli is one of the extremely acid-resistant bacteria and has a variety of acid-defense mechanisms. Here, we provide the first genome-wide overview of the adaptations of E. coli K-12 to mild and severe acid stress at both the transcriptional and translational levels. Using ribosome profiling and RNA sequencing, we uncover novel adaptations to different degrees of acidity, including previously hidden stress-induced small proteins and novel key transcription factors for acid defense, and report mRNAs with pH-dependent differential translation efficiency. In addition, we distinguish between acid-specific adaptations and general stress response mechanisms using denoising autoencoders. This workflow represents a powerful approach that takes advantage of next-generation sequencing techniques and machine learning to systematically analyze bacterial stress responses.

Vector-borne Trypanosoma brucei parasites develop in artificial human skin and persist as skin tissue forms

Transmission of Trypanosoma brucei by tsetse flies involves the deposition of the cell cycle-arrested metacyclic life cycle stage into mammalian skin at the site of the fly's bite. We introduce an advanced human skin equivalent and use tsetse flies to naturally infect the skin with trypanosomes. We detail the chronological order of the parasites' development in the skin by single-cell RNA sequencing and find a rapid activation of metacyclic trypanosomes and differentiation to proliferative parasites. Here we show that after the establishment of a proliferative population, the parasites enter a reversible quiescent state characterized by slow replication and a strongly reduced metabolism. We term these quiescent trypanosomes skin tissue forms, a parasite population that may play an important role in maintaining the infection over long time periods and in asymptomatic infected individuals.

DksA is a conserved master regulator of stress response in Acinetobacter baumannii

Coordination of bacterial stress response mechanisms is critical for long-term survival in harsh environments for successful host infection. The general and specific stress responses of well-studied Gram-negative pathogens like Escherichia coli are controlled by alternative sigma factors, archetypically RpoS. The deadly hospital pathogen Acinetobacter baumannii is notoriously resistant to environmental stresses, yet it lacks RpoS, and the molecular mechanisms driving this incredible stress tolerance remain poorly defined. Here, using functional genomics, we identified the transcriptional regulator DksA as a master regulator for broad stress protection and virulence in A. baumannii. Transcriptomics, phenomics and in vivo animal studies revealed that DksA controls ribosomal protein expression, metabolism, mutation rates, desiccation, antibiotic resistance, and host colonization in a niche-specific manner. Phylogenetically, DksA was highly conserved and well-distributed across Gammaproteobacteria, with 96.6% containing DksA, spanning 88 families. This study lays the groundwork for understanding DksA as a major regulator of general stress response and virulence in this important pathogen.

Design and off-target prediction for antisense oligomers targeting bacterial mRNAs with the MASON web server

Antisense oligomers (ASOs), such as peptide nucleic acids (PNAs), designed to inhibit the translation of essential bacterial genes, have emerged as attractive sequence- and species-specific programmable RNA antibiotics. Yet, potential drawbacks include unwanted side effects caused by their binding to transcripts other than the intended target. To facilitate the design of PNAs with minimal off-target effects, we developed MASON (make antisense oligomers now), a web server for the design of PNAs that target bacterial mRNAs. MASON generates PNA sequences complementary to the translational start site of a bacterial gene of interest and reports critical sequence attributes and potential off-target sites. We based MASON's off-target predictions on experiments in which we treated Salmonella enterica serovar Typhimurium with a series of 10-mer PNAs derived from a PNA targeting the essential gene acpP but carrying two serial mismatches. Growth inhibition and RNA-sequencing (RNA-seq) data revealed that PNAs with terminal mismatches are still able to target acpP, suggesting wider off-target effects than anticipated. Comparison of these results to an RNA-seq data set from uropathogenic Escherichia coli (UPEC) treated with eleven different PNAs confirmed that our findings are not unique to Salmonella We believe that MASON's off-target assessment will improve the design of specific PNAs and other ASOs.

An integrated transcriptomics-functional genomics approach reveals a small RNA that modulates Bacteroides thetaiotaomicron sensitivity to tetracyclines

Gene expression plasticity allows bacteria to adapt to diverse environments, tie their metabolism to available nutrients, and cope with stress. This is particularly relevant in a niche as dynamic and hostile as the human intestinal tract, yet transcriptional networks remain largely unknown in gut Bacteroides spp. Here, we map transcriptional units and profile their expression levels in Bacteroides thetaiotaomicron over a suite of 15 defined experimental conditions that are relevant in vivo , such as variation of temperature, pH, and oxygen tension, exposure to antibiotic stress, and growth on simple carbohydrates or on host mucin-derived glycans. Thereby, we infer stress- and carbon source-specific transcriptional regulons, including conditional expression of capsular polysaccharides and polysaccharide utilization loci, and expand the annotation of small regulatory RNAs (sRNAs) in this organism. Integrating this comprehensive expression atlas with transposon mutant fitness data, we identify conditionally important sRNAs. One example is MasB, whose inactivation led to increased bacterial tolerance of tetracyclines. Using MS2 affinity purification coupled with RNA sequencing, we predict targets of this sRNA and discuss their potential role in the context of the MasB-associated phenotype. Together, this transcriptomic compendium in combination with functional sRNA genomics-publicly available through a new iteration of the 'Theta-Base' web browser (www.helmholtz-hiri.de/en/datasets/bacteroides-v2)-constitutes a valuable resource for the microbiome and sRNA research communities alike.

Grad-seq analysis of Enterococcus faecalis and Enterococcus faecium provides a global view of RNA and protein complexes in these two opportunistic pathogens

Enterococcus faecalis and Enterococcus faecium are major nosocomial pathogens. Despite their relevance to public health and their role in the development of bacterial antibiotic resistance, relatively little is known about gene regulation in these species. RNA-protein complexes serve crucial functions in all cellular processes associated with gene expression, including post-transcriptional control mediated by small regulatory RNAs (sRNAs). Here, we present a new resource for the study of enterococcal RNA biology, employing the Grad-seq technique to comprehensively predict complexes formed by RNA and proteins in E. faecalis V583 and E. faecium AUS0004. Analysis of the generated global RNA and protein sedimentation profiles led to the identification of RNA-protein complexes and putative novel sRNAs. Validating our data sets, we observe well-established cellular RNA-protein complexes such as the 6S RNA-RNA polymerase complex, suggesting that 6S RNA-mediated global control of transcription is conserved in enterococci. Focusing on the largely uncharacterized RNA-binding protein KhpB, we use the RIP-seq technique to predict that KhpB interacts with sRNAs, tRNAs, and untranslated regions of mRNAs, and might be involved in the processing of specific tRNAs. Collectively, these datasets provide departure points for in-depth studies of the cellular interactome of enterococci that should facilitate functional discovery in these and related Gram-positive species. Our data are available to the community through a user-friendly Grad-seq browser that allows interactive searches of the sedimentation profiles (https://resources.helmholtz-hiri.de/gradseqef/).

INRI-seq enables global cell-free analysis of translation initiation and off-target effects of antisense inhibitors

Ribosome profiling (Ribo-seq) is a powerful method for the transcriptome-wide assessment of protein synthesis rates and the study of translational control mechanisms. Yet, Ribo-seq also has limitations. These include difficulties with the analysis of translation-modulating molecules such as antibiotics, which are often toxic or challenging to deliver into living cells. Here, we have developed in vitro Ribo-seq (INRI-seq), a cell-free method to analyze the translational landscape of a fully customizable synthetic transcriptome. Using Escherichia coli as an example, we show how INRI-seq can be used to analyze the translation initiation sites of a transcriptome of interest. We also study the global impact of direct translation inhibition by antisense peptide nucleic acid (PNA) to analyze PNA off-target effects. Overall, INRI-seq presents a scalable, sensitive method to study translation initiation in a transcriptome-wide manner without the potentially confounding effects of extracting ribosomes from living cells.

Ushering in a new era of single-cell transcriptomics in bacteria

Transcriptome analysis of individual cells by single-cell RNA-seq (scRNA-seq) has become routine for eukaryotic tissues, even being applied to whole multicellular organisms. In contrast, developing methods to read the transcriptome of single bacterial cells has proven more challenging, despite a general perception of bacteria as much simpler than eukaryotes. Bacterial cells are harder to lyse, their RNA content is about two orders of magnitude lower than that of eukaryotic cells, and bacterial mRNAs are less stable than their eukaryotic counterparts. Most importantly, bacterial transcripts lack functional poly(A) tails, precluding simple adaptation of popular standard eukaryotic scRNA-seq protocols that come with the double advantage of specific mRNA amplification and concomitant depletion of rRNA. However, thanks to very recent breakthroughs in methodology, bacterial scRNA-seq is now feasible. This short review will discuss recently published bacterial scRNA-seq approaches (MATQ-seq, microSPLiT, and PETRI-seq) and a spatial transcriptomics approach based on multiplexed in situ hybridization (par-seqFISH). Together, these novel approaches will not only enable a new understanding of cell-to-cell variation in bacterial gene expression, they also promise a new microbiology by enabling high-resolution profiling of gene activity in complex microbial consortia such as the microbiome or pathogens as they invade, replicate, and persist in host tissue.

A target expression threshold dictates invader defense and prevents autoimmunity by CRISPR-Cas13

CRISPR-Cas systems must enact robust immunity against foreign genetic material without inducing cytotoxic autoimmunity. For type VI systems that use Cas13 nucleases and recognize RNA targets, immune activation requires extensive CRISPR RNA (crRNA) guide-target complementarity and a target-flanking motif. Here, we report a third requirement shaping the immune response: the expression of the target transcript exceeding a threshold. We found that endogenous non-essential transcripts targeted by crRNAs rarely elicited autoimmunity. Instead, autoimmune induction required over-expressing the targeted transcripts above a threshold. A genome-wide screen confirmed target expression levels as a global determinant of cytotoxic autoimmunity and revealed that this threshold shifts with each guide-target pair. This threshold further ensured defense against a lytic bacteriophage yet allowed the tolerance of a targeted beneficial gene expressed from an invading plasmid. These findings establish target expression levels as an additional criterion for immune defense by RNA-targeting CRISPR-Cas systems, preventing autoimmunity and distinguishing pathogenic and benign invaders.

Comprehensive analysis of PNA-based antisense antibiotics targeting various essential genes in uropathogenic Escherichia coli

Antisense peptide nucleic acids (PNAs) that target mRNAs of essential bacterial genes exhibit specific bactericidal effects in several microbial species, but our mechanistic understanding of PNA activity and their target gene spectrum is limited. Here, we present a systematic analysis of PNAs targeting 11 essential genes with varying expression levels in uropathogenic Escherichia coli (UPEC). We demonstrate that UPEC is susceptible to killing by peptide-conjugated PNAs, especially when targeting the widely-used essential gene acpP. Our evaluation yields three additional promising target mRNAs for effective growth inhibition, i.e.dnaB, ftsZ and rpsH. The analysis also shows that transcript abundance does not predict target vulnerability and that PNA-mediated growth inhibition is not universally associated with target mRNA depletion. Global transcriptomic analyses further reveal PNA sequence-dependent but also -independent responses, including the induction of envelope stress response pathways. Importantly, we show that 9mer PNAs are generally as effective in inhibiting bacterial growth as their 10mer counterparts. Overall, our systematic comparison of a range of PNAs targeting mRNAs of different essential genes in UPEC suggests important features for PNA design, reveals a general bacterial response to PNA conjugates and establishes the feasibility of using PNA antibacterials to combat UPEC.

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.

Global RNA profiles show target selectivity and physiological effects of peptide-delivered antisense antibiotics

Antisense peptide nucleic acids (PNAs) inhibiting mRNAs of essential genes provide a straight-forward way to repurpose our knowledge of bacterial regulatory RNAs for development of programmable species-specific antibiotics. While there is ample proof of PNA efficacy, their target selectivity and impact on bacterial physiology are poorly understood. Moreover, while antibacterial PNAs are typically designed to block mRNA translation, effects on target mRNA levels are not well-investigated. Here, we pioneer the use of global RNA-seq analysis to decipher PNA activity in a transcriptome-wide manner. We find that PNA-based antisense oligomer conjugates robustly decrease mRNA levels of the widely-used target gene, acpP, in Salmonella enterica, with limited off-target effects. Systematic analysis of several different PNA-carrier peptides attached not only shows different bactericidal efficiency, but also activation of stress pathways. In particular, KFF-, RXR- and Tat-PNA conjugates especially induce the PhoP/Q response, whereas the latter two additionally trigger several distinct pathways. We show that constitutive activation of the PhoP/Q response can lead to Tat-PNA resistance, illustrating the utility of RNA-seq for understanding PNA antibacterial activity. In sum, our study establishes an experimental framework for the design and assessment of PNA antimicrobials in the long-term quest to use these for precision editing of microbiota.

Functional analysis of colonization factor antigen I positive enterotoxigenic Escherichia coli identifies genes implicated in survival in water and host colonization

Enterotoxigenic Escherichia coli (ETEC) expressing the colonization pili CFA/I are common causes of diarrhoeal infections in humans. Here, we use a combination of transposon mutagenesis and transcriptomic analysis to identify genes and pathways that contribute to ETEC persistence in water environments and colonization of a mammalian host. ETEC persisting in water exhibit a distinct RNA expression profile from those growing in richer media. Multiple pathways were identified that contribute to water survival, including lipopolysaccharide biosynthesis and stress response regulons. The analysis also indicated that ETEC growing in vivo in mice encounter a bottleneck driving down the diversity of colonizing ETEC populations.

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.

Global identification of RsmA/N binding sites in Pseudomonas aeruginosa by in vivo UV CLIP-seq

Pseudomonas aeruginosa harbours two redundant RNA-binding proteins RsmA/RsmN (RsmA/N), which play a critical role in balancing acute and chronic infections. However, in vivo binding sites on target transcripts and the overall impact on the physiology remains unclear. In this study, we applied in vivo UV crosslinking immunoprecipitation followed by RNA-sequencing (UV CLIP-seq) to detect RsmA/N-binding sites at single-nucleotide resolution and mapped more than 500 binding sites to approximately 400 genes directly bound by RsmA/N in P. aeruginosa. This also verified the ANGGA sequence in apical loops skewed towards 5'UTRs as a consensus motif for RsmA/N binding. Genetic analysis combined with CLIP-seq results suggested previously unrecognized RsmA/N targets involved in LPS modification. Moreover, the RsmA/N-titrating RNAs RsmY/RsmZ may be positively regulated by the RsmA/N-mediated translational repression of their upstream regulators, thus providing a possible mechanistic explanation for homoeostasis of the Rsm system. Thus, our study provides a detailed view of RsmA/N-RNA interactions and a resource for further investigation of the pleiotropic effects of RsmA/N on gene expression in P. aeruginosa.

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.

A global data-driven census of Salmonella small proteins and their potential functions in bacterial virulence

Small proteins are an emerging class of gene products with diverse roles in bacterial physiology. However, a full understanding of their importance has been hampered by insufficient genome annotations and a lack of comprehensive characterization in microbes other than Escherichia coli. We have taken an integrative approach to accelerate the discovery of small proteins and their putative virulence-associated functions in Salmonella Typhimurium. We merged the annotated small proteome of Salmonella with new small proteins predicted with in silico and experimental approaches. We then exploited existing and newly generated global datasets that provide information on small open reading frame expression during infection of epithelial cells (dual RNA-seq), contribution to bacterial fitness inside macrophages (Transposon-directed insertion sequencing), and potential engagement in molecular interactions (Grad-seq). This integrative approach suggested a new role for the small protein MgrB beyond its known function in regulating PhoQ. We demonstrate a virulence and motility defect of a Salmonella ΔmgrB mutant and reveal an effect of MgrB in regulating the Salmonella transcriptome and proteome under infection-relevant conditions. Our study highlights the power of interpreting available 'omics' datasets with a focus on small proteins, and may serve as a blueprint for a data integration-based survey of small proteins in diverse bacteria.

Global discovery of bacterial RNA-binding proteins by RNase-sensitive gradient profiles reports a new FinO domain protein

RNA-binding proteins (RBPs) play important roles in bacterial gene expression and physiology but their true number and functional scope remain little understood even in model microbes. To advance global RBP discovery in bacteria, we here establish glycerol gradient sedimentation with RNase treatment and mass spectrometry (GradR). Applied to Salmonella enterica, GradR confirms many known RBPs such as CsrA, Hfq, and ProQ by their RNase-sensitive sedimentation profiles, and discovers the FopA protein as a new member of the emerging family of FinO/ProQ-like RBPs. FopA, encoded on resistance plasmid pCol1B9, primarily targets a small RNA associated with plasmid replication. The target suite of FopA dramatically differs from the related global RBP ProQ, revealing context-dependent selective RNA recognition by FinO-domain RBPs. Numerous other unexpected RNase-induced changes in gradient profiles suggest that cellular RNA helps to organize macromolecular complexes in bacteria. By enabling poly(A)-independent generic RBP discovery, GradR provides an important element in the quest to build a comprehensive catalog of microbial RBPs.

An amphipathic peptide with antibiotic activity against multidrug-resistant Gram-negative bacteria

Peptide antibiotics are an abundant and synthetically tractable source of molecular diversity, but they are often cationic and can be cytotoxic, nephrotoxic and/or ototoxic, which has limited their clinical development. Here we report structure-guided optimization of an amphipathic peptide, arenicin-3, originally isolated from the marine lugworm Arenicola marina. The peptide induces bacterial membrane permeability and ATP release, with serial passaging resulting in a mutation in mlaC, a phospholipid transport gene. Structure-based design led to AA139, an antibiotic with broad-spectrum in vitro activity against multidrug-resistant and extensively drug-resistant bacteria, including ESBL, carbapenem- and colistin-resistant clinical isolates. The antibiotic induces a 3–4 log reduction in bacterial burden in mouse models of peritonitis, pneumonia and urinary tract infection. Cytotoxicity and haemolysis of the progenitor peptide is ameliorated with AA139, and the ‘no observable adverse effect level’ (NOAEL) dose in mice is ~10-fold greater than the dose generally required for efficacy in the infection models.

Peptide antibiotics often display a very narrow therapeutic index. Here, the authors present an optimized peptide antibiotic with broad-spectrum in vitro activities, in vivo efficacy in multiple disease models against multidrug-resistant Gram-negative infections, and reduced toxicity.

A decade of advances in transposon-insertion sequencing

It has been 10 years since the introduction of modern transposon-insertion sequencing (TIS) methods, which combine genome-wide transposon mutagenesis with high-throughput sequencing to estimate the fitness contribution or essentiality of each genetic component in a bacterial genome. Four TIS variations were published in 2009: transposon sequencing (Tn-Seq), transposon-directed insertion site sequencing (TraDIS), insertion sequencing (INSeq) and high-throughput insertion tracking by deep sequencing (HITS). TIS has since become an important tool for molecular microbiologists, being one of the few genome-wide techniques that directly links phenotype to genotype and ultimately can assign gene function. In this Review, we discuss the recent applications of TIS to answer overarching biological questions. We explore emerging and multidisciplinary methods that build on TIS, with an eye towards future applications.

In this Review, several experts discuss progress in the decade since the development of transposon-based approaches for bacterial genetic screens. They describe how advances in both experimental technologies and analytical strategies are resulting in insights into diverse biological processes.

Functional analysis of Salmonella Typhi adaptation to survival in water

Contaminated water is a major risk factor associated with the transmission of Salmonella enterica serovar Typhi (S. Typhi), the aetiological agent of human typhoid. However, little is known about how this pathogen adapts to living in the aqueous environment. We used transcriptome analysis (RNA‐seq) and transposon mutagenesis (TraDIS) to characterize these adaptive changes and identify multiple genes that contribute to survival. Over half of the genes in the S. Typhi genome altered expression level within the first 24 h following transfer from broth culture to water, although relatively few did so in the first 30 min. Genes linked to central metabolism, stress associated with arrested proton motive force and respiratory chain factors changed expression levels. Additionally, motility and chemotaxis genes increased expression, consistent with a scavenging lifestyle. The viaB‐associated gene tviC encoding a glcNAc epimerase that is required for Vi polysaccharide biosynthesis was, along with several other genes, shown to contribute to survival in water. Thus, we define regulatory adaptation operating in S. Typhi that facilitates survival in water.

Rapid transcriptional responses to serum exposure are associated with sensitivity and resistance to antibody-mediated complement killing in invasive Salmonella Typhimurium ST313

Background: Salmonella Typhimurium ST313 exhibits signatures of adaptation to invasive human infection, including higher resistance to humoral immune responses than gastrointestinal isolates. Full resistance to antibody-mediated complement killing (serum resistance) among nontyphoidal Salmonellae is uncommon, but selection of highly resistant strains could compromise vaccine-induced antibody immunity. Here, we address the hypothesis that serum resistance is due to a distinct genotype or transcriptome response in S. Typhimurium ST313.

Methods: Six S. Typhimurium ST313 bloodstream isolates, three of which were antibody resistant, were studied. Genomic content (single nucleotide polymorphisms and larger chromosomal modifications) of the strains was determined by Illumina and PACBIO sequencing, and functionally characterized using RNA-seq, transposon directed insertion site sequencing (TraDIS), targeted gene deletion and transfer of selected point mutations in an attempt to identify features associated with serum resistance.

  Results: Sequence polymorphisms in genes from strains with atypical serum susceptibility when transferred from strains that were highly resistant or susceptible to a strain that exhibited intermediate susceptibility did not significantly alter serum killing phenotype. No large chromosomal modifications typified serum resistance or susceptibility. Genes required for resistance to serum identified by TraDIS and RNA-seq included those involved in exopolysaccharide synthesis, iron scavenging and metabolism. Most of the down-regulated genes were associated with membrane proteins. Resistant and susceptible strains had distinct transcriptional responses to serum, particularly related to genes responsible for polysaccharide biosynthesis. There was higher upregulation of wca locus genes, involved in the biosynthesis of colanic acid exopolysaccharide, in susceptible strains and increased expression of fepE, a regulator of very long-chain lipopolysaccharide in resistant strains.

Conclusion: Clinical isolates of S. Typhimurium ST313 exhibit distinct antibody susceptibility phenotypes that may be associated with changes in gene expression on exposure to serum.

A global genomic approach uncovers novel components for twitching motility-mediated biofilm expansion in Pseudomonas aeruginosa

Pseudomonas aeruginosa is an extremely successful pathogen able to cause both acute and chronic infections in a range of hosts, utilizing a diverse arsenal of cell-associated and secreted virulence factors. A major cell-associated virulence factor, the Type IV pilus (T4P), is required for epithelial cell adherence and mediates a form of surface translocation termed twitching motility, which is necessary to establish a mature biofilm and actively expand these biofilms. P. aeruginosa twitching motility-mediated biofilm expansion is a coordinated, multicellular behaviour, allowing cells to rapidly colonize surfaces, including implanted medical devices. Although at least 44 proteins are known to be involved in the biogenesis, assembly and regulation of the T4P, with additional regulatory components and pathways implicated, it is unclear how these components and pathways interact to control these processes. In the current study, we used a global genomics-based random-mutagenesis technique, transposon directed insertion-site sequencing (TraDIS), coupled with a physical segregation approach, to identify all genes implicated in twitching motility-mediated biofilm expansion in P. aeruginosa. Our approach allowed identification of both known and novel genes, providing new insight into the complex molecular network that regulates this process in P. aeruginosa. Additionally, our data suggest that the flagellum-associated gene products have a differential effect on twitching motility, based on whether components are intra- or extracellular. Overall the success of our TraDIS approach supports the use of this global genomic technique for investigating virulence genes in bacterial pathogens.

Morphological, genomic and transcriptomic responses of Klebsiella pneumoniae to the last-line antibiotic colistin

Colistin remains one of the few antibiotics effective against multi-drug resistant (MDR) hospital pathogens, such as Klebsiella pneumoniae. Yet resistance to this last-line drug is rapidly increasing. Characterized mechanisms of col^R in K. pneumoniae are largely due to chromosomal mutations in two-component regulators, although a plasmid-mediated col^R mechanism has recently been uncovered. However, the effects of intrinsic colistin resistance are yet to be characterized on a whole-genome level. Here, we used a genomics-based approach to understand the mechanisms of adaptive col^R acquisition in K. pneumoniae. In controlled directed-evolution experiments we observed two distinct paths to colistin resistance acquisition. Whole genome sequencing identified mutations in two colistin resistance genes: in the known col^R regulator phoQ which became fixed in the population and resulted in a single amino acid change, and unstable minority variants in the recently described two-component sensor crrB. Through RNAseq and microscopy, we reveal the broad range of effects that colistin exposure has on the cell. This study is the first to use genomics to identify a population of minority variants with mutations in a col^R gene in K. pneumoniae.

Machine learning identifies signatures of host adaptation in the bacterial pathogen Salmonella enterica

Emerging pathogens are a major threat to public health, however understanding how pathogens adapt to new niches remains a challenge. New methods are urgently required to provide functional insights into pathogens from the massive genomic data sets now being generated from routine pathogen surveillance for epidemiological purposes. Here, we measure the burden of atypical mutations in protein coding genes across independently evolved Salmonella enterica lineages, and use these as input to train a random forest classifier to identify strains associated with extraintestinal disease. Members of the species fall along a continuum, from pathovars which cause gastrointestinal infection and low mortality, associated with a broad host-range, to those that cause invasive infection and high mortality, associated with a narrowed host range. Our random forest classifier learned to perfectly discriminate long-established gastrointestinal and invasive serovars of Salmonella. Additionally, it was able to discriminate recently emerged Salmonella Enteritidis and Typhimurium lineages associated with invasive disease in immunocompromised populations in sub-Saharan Africa, and within-host adaptation to invasive infection. We dissect the architecture of the model to identify the genes that were most informative of phenotype, revealing a common theme of degradation of metabolic pathways in extraintestinal lineages. This approach accurately identifies patterns of gene degradation and diversifying selection specific to invasive serovars that have been captured by more labour-intensive investigations, but can be readily scaled to larger analyses.

Molecular phenotyping of infection-associated small non-coding RNAs

Infection is a complicated balance, with both pathogen and host struggling to tilt the result in their favour. Bacterial infection biology has relied on forward genetics for many of its advances, defining phenotype in terms of replication in model systems. However, many known virulence factors fail to produce robust phenotypes, particularly in the systems most amenable to genetic manipulation, such as cell-culture models. This has particularly been limiting for the study of the bacterial regulatory small RNAs (sRNAs) in infection. We argue that new sequencing-based technologies can work around this problem by providing a 'molecular phenotype', defined in terms of the specific transcriptional dysregulation in the infection system induced by gene deletion. We illustrate this using the example of our recent study of the PinT sRNA using dual RNA-seq, that is, simultaneous RNA sequencing of host and pathogen during infection. We additionally discuss how other high-throughput technologies, in particular genetic interaction mapping using transposon insertion sequencing, may be used to further dissect molecular phenotypes. We propose a strategy for how high-throughput technologies can be integrated in the study of non-coding regulators as well as bacterial virulence factors, enhancing our ability to rapidly generate hypotheses with regards to their function.This article is part of the themed issue 'The new bacteriology'.

The primary transcriptome of Neisseria meningitidis and its interaction with the RNA chaperone Hfq

Neisseria meningitidis is a human commensal that can also cause life-threatening meningitis and septicemia. Despite growing evidence for RNA-based regulation in meningococci, their transcriptome structure and output of regulatory small RNAs (sRNAs) are incompletely understood. Using dRNA-seq, we have mapped at single-nucleotide resolution the primary transcriptome of N. meningitidis strain 8013. Annotation of 1625 transcriptional start sites defines transcription units for most protein-coding genes but also reveals a paucity of classical σ70-type promoters, suggesting the existence of activators that compensate for the lack of −35 consensus sequences in N. meningitidis. The transcriptome maps also reveal 65 candidate sRNAs, a third of which were validated by northern blot analysis. Immunoprecipitation with the RNA chaperone Hfq drafts an unexpectedly large post-transcriptional regulatory network in this organism, comprising 23 sRNAs and hundreds of potential mRNA targets. Based on this data, using a newly developed gfp reporter system we validate an Hfq-dependent mRNA repression of the putative colonization factor PrpB by the two trans-acting sRNAs RcoF1/2. Our genome-wide RNA compendium will allow for a better understanding of meningococcal transcriptome organization and riboregulation with implications for colonization of the human nasopharynx.

Resolving host-pathogen interactions by dual RNA-seq

The transcriptome is a powerful proxy for the physiological state of a cell, healthy or diseased. As a result, transcriptome analysis has become a key tool in understanding the molecular changes that accompany bacterial infections of eukaryotic cells. Until recently, such transcriptomic studies have been technically limited to analyzing mRNA expression changes in either the bacterial pathogen or the infected eukaryotic host cell. However, the increasing sensitivity of high-throughput RNA sequencing now enables "dual RNA-seq" studies, simultaneously capturing all classes of coding and noncoding transcripts in both the pathogen and the host. In the five years since the concept of dual RNA-seq was introduced, the technique has been applied to a range of infection models. This has not only led to a better understanding of the physiological changes in pathogen and host during the course of an infection but has also revealed hidden molecular phenotypes of virulence-associated small noncoding RNAs that were not visible in standard infection assays. Here, we use the knowledge gained from these recent studies to suggest experimental and computational guidelines for the design of future dual RNA-seq studies. We conclude this review by discussing prospective applications of the technique.

A profile-based method for identifying functional divergence of orthologous genes in bacterial genomes

Motivation: Next generation sequencing technologies have provided us with a wealth of information on genetic variation, but predicting the functional significance of this variation is a difficult task. While many comparative genomics studies have focused on gene flux and large scale changes, relatively little attention has been paid to quantifying the effects of single nucleotide polymorphisms and indels on protein function, particularly in bacterial genomics.

Results: We present a hidden Markov model based approach we call delta-bitscore (DBS) for identifying orthologous proteins that have diverged at the amino acid sequence level in a way that is likely to impact biological function. We benchmark this approach with several widely used datasets and apply it to a proof-of-concept study of orthologous proteomes in an investigation of host adaptation in Salmonella enterica. We highlight the value of the method in identifying functional divergence of genes, and suggest that this tool may be a better approach than the commonly used dN/dS metric for identifying functionally significant genetic changes occurring in recently diverged organisms.

Availability and Implementation: A program implementing DBS for pairwise genome comparisons is freely available at: https://github.com/UCanCompBio/deltaBS.

Contact:nicole.wheeler@pg.canterbury.ac.nz or lars.barquist@uni-wuerzburg.de

Supplementary information:Supplementary data are available at Bioinformatics online.

The in vitro and in vivo effects of constitutive light expression on a bioluminescent strain of the mouse enteropathogen Citrobacter rodentium

Bioluminescent reporter genes, such as those from fireflies and bacteria, let researchers use light production as a non-invasive and non-destructive surrogate measure of microbial numbers in a wide variety of environments. As bioluminescence needs microbial metabolites, tagging microorganisms with luciferases means only live metabolically active cells are detected. Despite the wide use of bioluminescent reporter genes, very little is known about the impact of continuous (also called constitutive) light expression on tagged bacteria. We have previously made a bioluminescent strain of Citrobacter rodentium, a bacterium which infects laboratory mice in a similar way to how enteropathogenic Escherichia coli (EPEC) and enterohaemorrhagic E. coli (EHEC) infect humans. In this study, we compared the growth of the bioluminescent C. rodentium strain ICC180 with its non-bioluminescent parent (strain ICC169) in a wide variety of environments. To understand more about the metabolic burden of expressing light, we also compared the growth profiles of the two strains under approximately 2,000 different conditions. We found that constitutive light expression in ICC180 was near-neutral in almost every non-toxic environment tested. However, we also found that the non-bioluminescent parent strain has a competitive advantage over ICC180 during infection of adult mice, although this was not enough for ICC180 to be completely outcompeted. In conclusion, our data suggest that constitutive light expression is not metabolically costly to C. rodentium and supports the view that bioluminescent versions of microbes can be used as a substitute for their non-bioluminescent parents to study bacterial behaviour in a wide variety of environments.

Studying RNA Homology and Conservation with Infernal: From Single Sequences to RNA Families

Emerging high-throughput technologies have led to a deluge of putative non-coding RNA (ncRNA) sequences identified in a wide variety of organisms. Systematic characterization of these transcripts will be a tremendous challenge. Homology detection is critical to making maximal use of functional information gathered about ncRNAs: identifying homologous sequence allows us to transfer information gathered in one organism to another quickly and with a high degree of confidence. ncRNA presents a challenge for homology detection, as the primary sequence is often poorly conserved and de novo secondary structure prediction and search remain difficult. This unit introduces methods developed by the Rfam database for identifying "families" of homologous ncRNAs starting from single "seed" sequences, using manually curated sequence alignments to build powerful statistical models of sequence and structure conservation known as covariance models (CMs), implemented in the Infernal software package. We provide a step-by-step iterative protocol for identifying ncRNA homologs and then constructing an alignment and corresponding CM. We also work through an example for the bacterial small RNA MicA, discovering a previously unreported family of divergent MicA homologs in genus Xenorhabdus in the process. © 2016 by John Wiley & Sons, Inc.

Global RNA recognition patterns of post-transcriptional regulators Hfq and CsrA revealed by UV crosslinking in vivo

The molecular roles of many RNA‐binding proteins in bacterial post‐transcriptional gene regulation are not well understood. Approaches combining in vivo UV crosslinking with RNA deep sequencing (CLIP‐seq) have begun to revolutionize the transcriptome‐wide mapping of eukaryotic RNA‐binding protein target sites. We have applied CLIP‐seq to chart the target landscape of two major bacterial post‐transcriptional regulators, Hfq and CsrA, in the model pathogen Salmonella Typhimurium. By detecting binding sites at single‐nucleotide resolution, we identify RNA preferences and structural constraints of Hfq and CsrA during their interactions with hundreds of cellular transcripts. This reveals 3′‐located Rho‐independent terminators as a universal motif involved in Hfq–RNA interactions. Additionally, Hfq preferentially binds 5′ to sRNA‐target sites in mRNAs, and 3′ to seed sequences in sRNAs, reflecting a simple logic in how Hfq facilitates sRNA–mRNA interactions. Importantly, global knowledge of Hfq sites significantly improves sRNA‐target predictions. CsrA binds AUGGA sequences in apical loops and targets many Salmonella virulence mRNAs. Overall, our generic CLIP‐seq approach will bring new insights into post‐transcriptional gene regulation by RNA‐binding proteins in diverse bacterial species.

The TraDIS toolkit: sequencing and analysis for dense transposon mutant libraries

Summary: Transposon insertion sequencing is a high-throughput technique for assaying large libraries of otherwise isogenic transposon mutants providing insight into gene essentiality, gene function and genetic interactions. We previously developed the Transposon Directed Insertion Sequencing (TraDIS) protocol for this purpose, which utilizes shearing of genomic DNA followed by specific PCR amplification of transposon-containing fragments and Illumina sequencing. Here we describe an optimized high-yield library preparation and sequencing protocol for TraDIS experiments and a novel software pipeline for analysis of the resulting data. The Bio-Tradis analysis pipeline is implemented as an extensible Perl library which can either be used as is, or as a basis for the development of more advanced analysis tools. This article can serve as a general reference for the application of the TraDIS methodology.

Availability and implementation: The optimized sequencing protocol is included as supplementary information. The Bio-Tradis analysis pipeline is available under a GPL license at https://github.com/sanger-pathogens/Bio-Tradis

Contact:parkhill@sanger.ac.uk

Supplementary information:Supplementary data are available at Bioinformatics online.

Dual RNA-seq unveils noncoding RNA functions in host-pathogen interactions

Bacteria express many small RNAs for which the regulatory roles in pathogenesis have remained poorly understood due to a paucity of robust phenotypes in standard virulence assays. Here we use a generic 'dual RNA-seq' approach to profile RNA expression simultaneously in pathogen and host during Salmonella enterica serovar Typhimurium infection and reveal the molecular impact of bacterial riboregulators. We identify a PhoP-activated small RNA, PinT, which upon bacterial internalization temporally controls the expression of both invasion-associated effectors and virulence genes required for intracellular survival. This riboregulatory activity causes pervasive changes in coding and noncoding transcripts of the host. Interspecies correlation analysis links PinT to host cell JAK-STAT signalling, and we identify infection-specific alterations in multiple long noncoding RNAs. Our study provides a paradigm for a sensitive RNA-based analysis of intracellular bacterial pathogens and their hosts without physical separation, as well as a new discovery route for hidden functions of pathogen genes.

Accelerating Discovery and Functional Analysis of Small RNAs with New Technologies

Over the past decade, bacterial small RNAs (sRNAs) have gone from a biological curiosity to being recognized as a major class of regulatory molecules. High-throughput methods for sampling the transcriptional output of bacterial cells demonstrate that sRNAs are universal features of bacterial transcriptomes, are plentiful, and appear to vary extensively over evolutionary time. With ever more bacteria coming under study, the question becomes how can we accelerate the discovery and functional characterization of sRNAs in diverse organisms. New technologies built on high-throughput sequencing are emerging that can rapidly provide global insight into the numbers and functions of sRNAs in bacteria of interest, providing information that can shape hypotheses and guide research. In this review, we describe recent developments in transcriptomics (RNA-seq) and functional genomics that we expect to help us develop an integrated, systems-level view of sRNA biology in bacteria.

Signatures of adaptation in human invasive Salmonella Typhimurium ST313 populations from sub-Saharan Africa

Two lineages of Salmonella enterica serovar Typhimurium (S. Typhimurium) of multi-locus sequence type ST313 have been linked with the emergence of invasive Salmonella disease across sub-Saharan Africa. The expansion of these lineages has a temporal association with the HIV pandemic and antibiotic usage. We analysed the whole genome sequence of 129 ST313 isolates representative of the two lineages and found evidence of lineage-specific genome degradation, with some similarities to that observed in S. Typhi. Individual ST313 S. Typhimurium isolates exhibit a distinct metabolic signature and modified enteropathogenesis in both a murine and cattle model of colitis, compared to S. Typhimurium outside of the ST313 lineages. These data define phenotypes that distinguish ST313 isolates from other S. Typhimurium and may represent adaptation to a distinct pathogenesis and lifestyle linked to an-immuno-compromised human population.

Robust identification of noncoding RNA from transcriptomes requires phylogenetically-informed sampling

Noncoding RNAs are integral to a wide range of biological processes, including translation, gene regulation, host-pathogen interactions and environmental sensing. While genomics is now a mature field, our capacity to identify noncoding RNA elements in bacterial and archaeal genomes is hampered by the difficulty of de novo identification. The emergence of new technologies for characterizing transcriptome outputs, notably RNA-seq, are improving noncoding RNA identification and expression quantification. However, a major challenge is to robustly distinguish functional outputs from transcriptional noise. To establish whether annotation of existing transcriptome data has effectively captured all functional outputs, we analysed over 400 publicly available RNA-seq datasets spanning 37 different Archaea and Bacteria. Using comparative tools, we identify close to a thousand highly-expressed candidate noncoding RNAs. However, our analyses reveal that capacity to identify noncoding RNA outputs is strongly dependent on phylogenetic sampling. Surprisingly, and in stark contrast to protein-coding genes, the phylogenetic window for effective use of comparative methods is perversely narrow: aggregating public datasets only produced one phylogenetic cluster where these tools could be used to robustly separate unannotated noncoding RNAs from a null hypothesis of transcriptional noise. Our results show that for the full potential of transcriptomics data to be realized, a change in experimental design is paramount: effective transcriptomics requires phylogeny-aware sampling.

We have analysed more than 400 public transcriptomes, generated using RNA-seq, from almost 40 strains of Bacteria and Archaea. We discovered that the capacity to identify noncoding RNA outputs from this data is strongly dependent on phylogenetic sampling. Our results show that, for the full potential of transcriptomics data as a discovery tool to be realized, a change in experimental design is critical: effective comparative transcriptomics requires phylogeny-aware sampling. We also examined how comparative transcriptomics experiments can be used to effectively identify RNA elements. We find that, for RNA element discovery, a phylogeny-informed sampling approach is more effective than analyses of individual species. Phylogeny-informed sampling reveals a narrow ‘Goldilocks Zone’ (where species are not too similar and not too divergent) for RNA identification using clusters of related species. In stark contrast to protein-coding genes, not only is the phylogenetic window for the effective use of comparative methods for noncoding RNA identification perversely narrow, but few existing datasets sit within this Goldilocks Zone: by aggregating public datasets, we were only able to create one phylogenetic cluster where comparative tools could be used to confidently separate unannotated noncoding RNAs from transcriptional noise.

Functional genomics reveals that Clostridium difficile Spo0A coordinates sporulation, virulence and metabolism

BACKGROUND

Clostridium difficile is an anaerobic, Gram-positive bacterium that can reside as a commensal within the intestinal microbiota of healthy individuals or cause life-threatening antibiotic-associated diarrhea in immunocompromised hosts. C. difficile can also form highly resistant spores that are excreted facilitating host-to-host transmission. The C. difficile spo0A gene encodes a highly conserved transcriptional regulator of sporulation that is required for relapsing disease and transmission in mice.

RESULTS

Here we describe a genome-wide approach using a combined transcriptomic and proteomic analysis to identify Spo0A regulated genes. Our results validate Spo0A as a positive regulator of putative and novel sporulation genes as well as components of the mature spore proteome. We also show that Spo0A regulates a number of virulence-associated factors such as flagella and metabolic pathways including glucose fermentation leading to butyrate production.

CONCLUSIONS

The C. difficile spo0A gene is a global transcriptional regulator that controls diverse sporulation, virulence and metabolic phenotypes coordinating pathogen adaptation to a wide range of host interactions. Additionally, the rich breadth of functional data allowed us to significantly update the annotation of the C. difficile 630 reference genome which will facilitate basic and applied research on this emerging pathogen.

Dominant role of nucleotide substitution in the diversification of serotype 3 pneumococci over decades and during a single infection

Streptococcus pneumoniae of serotype 3 possess a mucoid capsule and cause disease associated with high mortality rates relative to other pneumococci. Phylogenetic analysis of a complete reference genome and 81 draft sequences from clonal complex 180, the predominant serotype 3 clone in much of the world, found most sampled isolates belonged to a clade affected by few diversifying recombinations. However, other isolates indicate significant genetic variation has accumulated over the clonal complex's entire history. Two closely related genomes, one from the blood and another from the cerebrospinal fluid, were obtained from a patient with meningitis. The pair differed in their behaviour in a mouse model of disease and in their susceptibility to antimicrobials, with at least some of these changes attributable to a mutation that up-regulated the patAB efflux pump. This indicates clinically important phenotypic variation can accumulate rapidly through small alterations to the genotype.

Approaches to querying bacterial genomes with transposon-insertion sequencing

In this review, we discuss transposon-insertion sequencing, variously known in the literature as TraDIS, Tn-seq, INSeq, and HITS. By monitoring a large library of single transposon-insertion mutants with high-throughput sequencing, these methods can rapidly identify genomic regions that contribute to organismal fitness under any condition assayable in the laboratory with exquisite resolution. We discuss the various protocols that have been developed and methods for analysis. We provide an overview of studies that have examined the reproducibility and accuracy of these methods, as well as studies showing the advantages offered by the high resolution and dynamic range of high-throughput sequencing over previous methods. We review a number of applications in the literature, from predicting genes essential for in vitro growth to directly assaying requirements for survival under infective conditions in vivo. We also highlight recent progress in assaying non-coding regions of the genome in addition to known coding sequences, including the combining of RNA-seq with high-throughput transposon mutagenesis.

A comparison of dense transposon insertion libraries in the Salmonella serovars Typhi and Typhimurium

Salmonella Typhi and Typhimurium diverged only ∼50 000 years ago, yet have very different host ranges and pathogenicity. Despite the availability of multiple whole-genome sequences, the genetic differences that have driven these changes in phenotype are only beginning to be understood. In this study, we use transposon-directed insertion-site sequencing to probe differences in gene requirements for competitive growth in rich media between these two closely related serovars. We identify a conserved core of 281 genes that are required for growth in both serovars, 228 of which are essential in Escherichia coli. We are able to identify active prophage elements through the requirement for their repressors. We also find distinct differences in requirements for genes involved in cell surface structure biogenesis and iron utilization. Finally, we demonstrate that transposon-directed insertion-site sequencing is not only applicable to the protein-coding content of the cell but also has sufficient resolution to generate hypotheses regarding the functions of non-coding RNAs (ncRNAs) as well. We are able to assign probable functions to a number of cis-regulatory ncRNA elements, as well as to infer likely differences in trans-acting ncRNA regulatory networks.