Person:

Livny, Jonathan

Background: Regulation of bacterial gene expression by small RNAs (sRNAs) have proved to be important for many biological processes. Francisella tularensis is a highly pathogenic Gram-negative bacterium that causes the disease tularaemia in humans and animals. Relatively little is known about the regulatory networks existing in this organism that allows it to survive in a wide array of environments and no sRNA regulators have been identified so far. Results: We have used a combination of experimental assays and in silico prediction to identify sRNAs in F. tularensis strain LVS. Using a cDNA cloning and sequencing approach we have shown that F. tularensis expresses homologues of several sRNAs that are well-conserved among diverse bacteria. We have also discovered two abundant putative sRNAs that share no sequence similarity or conserved genomic context with any previously annotated regulatory transcripts. Deletion of either of these two loci led to significant changes in the expression of several mRNAs that likely include the cognate target(s) of these sRNAs. Deletion of these sRNAs did not, however, significantly alter F. tularensis growth under various stress conditions in vitro, its replication in murine cells, or its ability to induce disease in a mouse model of F. tularensis infection. We also conducted a genome-wide in silico search for intergenic loci that suggests F. tularensis encodes several other sRNAs in addition to the sRNAs found in our experimental screen. Conclusion: Our findings suggest that F. tularensis encodes a significant number of non-coding regulatory RNAs, including members of well conserved families of structural and housekeeping RNAs and other poorly conserved transcripts that may have evolved more recently to help F. tularensis deal with the unique and diverse set of environments with which it must contend.

The coupling of high-density transposon mutagenesis to high-throughput DNA sequencing (transposon-insertion sequencing) enables simultaneous and genome-wide assessment of the contributions of individual loci to bacterial growth and survival. We have refined analysis of transposon-insertion sequencing data by normalizing for the effect of DNA replication on sequencing output and using a hidden Markov model (HMM)-based filter to exploit heretofore unappreciated information inherent in all transposon-insertion sequencing data sets. The HMM can smooth variations in read abundance and thereby reduce the effects of read noise, as well as permit fine scale mapping that is independent of genomic annotation and enable classification of loci into several functional categories (e.g. essential, domain essential or ‘sick’). We generated a high-resolution map of genomic loci (encompassing both intra- and intergenic sequences) that are required or beneficial for in vitro growth of the cholera pathogen, Vibrio cholerae. This work uncovered new metabolic and physiologic requirements for V. cholerae survival, and by combining transposon-insertion sequencing and transcriptomic data sets, we also identified several novel noncoding RNA species that contribute to V. cholerae growth. Our findings suggest that HMM-based approaches will enhance extraction of biological meaning from transposon-insertion sequencing genomic data.

Vibrio parahaemolyticus is the leading worldwide cause of seafood-associated gastroenteritis, yet little is known regarding its intraintestinal gene expression or physiology. To date, in vivo analyses have focused on identification and characterization of virulence factors—e.g. a crucial Type III secretion system (T3SS2)—rather than genome-wide analyses of in vivo biology. Here, we used RNA-Seq to profile V. parahaemolyticus gene expression in infected infant rabbits, which mimic human infection. Comparative transcriptomic analysis of V. parahaemolyticus isolated from rabbit intestines and from several laboratory conditions enabled identification of mRNAs and sRNAs induced during infection and of regulatory factors that likely control them. More than 12% of annotated V. parahaemolyticus genes are differentially expressed in the intestine, including the genes of T3SS2, which are likely induced by bile-mediated activation of the transcription factor VtrB. Our analyses also suggest that V. parahaemolyticus has access to glucose or other preferred carbon sources in vivo, but that iron is inconsistently available. The V. parahaemolyticus transcriptional response to in vivo growth is far more widespread than and largely distinct from that of V. cholerae, likely due to the distinct ways in which these diarrheal pathogens interact with and modulate the environment in the small intestine.

Background: High-throughput sequencing of cDNA libraries (RNA-Seq) has proven to be a highly effective approach for studying bacterial transcriptomes. A central challenge in designing RNA-Seq-based experiments is estimating a priori the number of reads per sample needed to detect and quantify thousands of individual transcripts with a large dynamic range of abundance. Results: We have conducted a systematic examination of how changes in the number of RNA-Seq reads per sample influences both profiling of a single bacterial transcriptome and the comparison of gene expression among samples. Our findings suggest that the number of reads typically produced in a single lane of the Illumina HiSeq sequencer far exceeds the number needed to saturate the annotated transcriptomes of diverse bacteria growing in monoculture. Moreover, as sequencing depth increases, so too does the detection of cDNAs that likely correspond to spurious transcripts or genomic DNA contamination. Finally, even when dozens of barcoded individual cDNA libraries are sequenced in a single lane, the vast majority of transcripts in each sample can be detected and numerous genes differentially expressed between samples can be identified. Conclusions: Our analysis provides a guide for the many researchers seeking to determine the appropriate sequencing depth for RNA-Seq-based studies of diverse bacterial species.

Bacterial small RNAs (sRNAs) are short transcripts that typically do not encode proteins and often act as regulators of gene expression through a variety of mechanisms. Regulatory sRNAs have been identified in many species, including Mycobacterium tuberculosis, the causative agent of tuberculosis. Here, we use a computational algorithm to predict sRNA candidates in the mycobacterial species M. smegmatis and M. bovis BCG and confirmed the expression of many sRNAs using Northern blotting. Thus, we have identified 17 and 23 novel sRNAs in M. smegmatis and M. bovis BCG, respectively. We have also applied a high-throughput technique (Deep-RACE) to map the 5′ and 3′ ends of many of these sRNAs and identified potential regulators of sRNAs by analysis of existing ChIP-seq datasets. The sRNAs identified in this work likely contribute to the unique biology of mycobacteria.

The ability to simultaneously characterize the bacterial and host expression programs during infection would facilitate a comprehensive understanding of pathogen-host interactions. While RNA-Seq has greatly advanced our ability to study the transcriptomes of prokaryote and eukaryotes separately, limitations in existing protocols for generating and analyzing RNA-Seq data have hindered simultaneous profiling of host and bacterial pathogen transcripts from the same sample. Here we provide a detailed protocol for simultaneous analysis of host and bacterial transcripts by RNA-Seq. Importantly, this protocol details the steps required for efficient host and bacteria lysis, barcoding of samples, technical advances in sample preparation for low yield sample inputs and a computational pipeline to analyze both mammalian and microbial reads from mixed hostpathogen RNA-Seq data. Sample preparation takes 3 d from cultured cells to pooled libraries. Data analysis takes an additional day. Compared with previous methods, the protocol detailed here provides a sensitive, facile and generalizable approach, suitable for large-scale studies, which will enable the field to obtain in-depth analysis of hostpathogen interactions in infection models.

Regulatory 5′ untranslated regions (r5′UTRs) of mRNAs such as riboswitches modulate the expression of genes involved in varied biological processes in both bacteria and eukaryotes. New high-throughput sequencing technologies could provide powerful tools for discovery of novel r5′UTRs, but the size and complexity of the datasets generated by these technologies makes it difficult to differentiate r5′UTRs from the multitude of other types of RNAs detected. Here, we developed and implemented a bioinformatic approach to identify putative r5′UTRs from within large datasets of RNAs recently identified by pyrosequencing of the Vibrio cholerae small transcriptome. This screen yielded only ∼1% of all non-overlapping RNAs along with 75% of previously annotated r5′UTRs and 69 candidate V. cholerae r5′UTRs. These candidates include several putative functional homologues of diverse r5′UTRs characterized in other species as well as numerous candidates upstream of genes involved in pathways not known to be regulated by r5′UTRs, such as fatty acid oxidation and peptidoglycan catabolism. Two of these novel r5′UTRs were experimentally validated using a GFP reporter-based approach. Our findings suggest that the number and diversity of pathways regulated by r5′UTRs has been underestimated and that deep sequencing-based transcriptomics will be extremely valuable in the search for novel r5′UTRs.

Tuberculosis (TB) is a major global health problem, infecting millions of people each year. The causative agent of TB, Mycobacterium tuberculosis, is one of the world’s most ancient and successful pathogens. However, until recently, no work on small regulatory RNAs had been performed in this organism. Regulatory RNAs are found in all three domains of life, and have already been shown to regulate virulence in well-known pathogens, such as Staphylococcus aureus and Vibrio cholera. Here we report the discovery of 34 novel small RNAs (sRNAs) in the TB-complex M. bovis BCG, using a combination of experimental and computational approaches. Putative homologues of many of these sRNAs were also identified in M. tuberculosis and/or M. smegmatis. Those sRNAs that are also expressed in the non-pathogenic M. smegmatis could be functioning to regulate conserved cellular functions. In contrast, those sRNAs identified specifically in M. tuberculosis could be functioning in mediation of virulence, thus rendering them potential targets for novel antimycobacterials. Various features and regulatory aspects of some of these sRNAs are discussed.