Omri Wurtzel

Mapping the Burkholderia cenocepacia niche response via high-throughput sequencing

Determining how an organism responds to its environment by altering gene expression is key to understanding its ecology. Here, we used RNA-seq to comprehensively and quantitatively assess the transcriptional response of the bacterial opportunistic cystic fibrosis (CF) pathogen and endemic soil dweller, Burkholderia cenocepacia, in conditions mimicking these 2 environments. By sequencing 762 million bases of cDNA from 2 closely related B. cenocepacia strains (one isolated from a CF patient and one from soil), we identified a number of potential virulence factors expressed under CF-like conditions, whereas genes whose protein products are involved in nitrogen scavenging and 2-component sensing were among those induced under soil-like conditions. Interestingly, 13 new putative noncoding RNAs were discovered using this technique, 12 of which are preferentially induced in the soil environment, suggesting that ncRNAs play an important role in survival in the soil. In addition, we detected a surprisingly large number of regulatory differences between the 2 strains, which may represent specific adaptations to the niches from which each strain was isolated, despite their high degree of DNA sequence similarity. Compared with the CF strain, the soil strain shows a stronger global gene expression response to its environment, which is consistent with the need for a more dynamic reaction to the heterogeneous conditions of soil.

Genomic island variability facilitates Prochlorococcus-virus coexistence

Prochlorococcus cyanobacteria are extremely abundant in the oceans, as are the viruses that infect them. How hosts and viruses coexist in nature remains unclear, although the presence of both susceptible and resistant cells may allow this coexistence. Combined whole-genome sequencing and PCR screening technology now enables us to investigate the effect of resistance on genome evolution and the genomic mechanisms behind the long-term coexistence of Prochlorococcus and their viruses. Here we present a genome analysis of 77 substrains selected for resistance to ten viruses, revealing mutations primarily in non-conserved, horizontally transferred genes that localize to a single hypervariable genomic island. Mutations affected viral attachment to the cell surface and imposed a fitness cost to the host, manifested by significantly lower growth rates or a previously unknown mechanism of more rapid infection by other viruses. The mutant genes are generally uncommon in nature yet some carry polymorphisms matching those found experimentally. These data are empirical evidence indicating that viral-attachment genes are preferentially located in genomic islands and that viruses are a selective pressure enhancing the diversity of both island genes and island gene content. This diversity emerges as a genomic mechanism that reduces the effective host population size for infection by a given virus, thus facilitating long-term coexistence between viruses and their hosts in nature.

Comparative transcriptomics of pathogenic and non‐pathogenic Listeria species

Listeria monocytogenes is a human, food‐borne pathogen. Genomic comparisons between L. monocytogenes and Listeria innocua, a closely related non‐pathogenic species, were pivotal in the identification of protein‐coding genes essential for virulence. However, no comprehensive comparison has focused on the non‐coding genome. We used strand‐specific cDNA sequencing to produce genome‐wide transcription start site maps for both organisms, and developed a publicly available integrative browser to visualize and analyze both transcriptomes in different growth conditions and genetic backgrounds. Our data revealed conservation across most transcripts, but significant divergence between the species in a subset of non‐coding RNAs. In L. monocytogenes, we identified 113 small RNAs (33 novel) and 70 antisense RNAs (53 novel), significantly increasing the repertoire of ncRNAs in this species. Remarkably, we identified a class of long antisense transcripts (lasRNAs) that overlap one gene while also serving as the 5′ UTR of the adjacent divergent gene. Experimental evidence suggests that lasRNAs transcription inhibits expression of one operon while activating the expression of another. Such a lasRNA/operon structure, that we named ‘excludon’, might represent a novel form of regulation in bacteria.

Validation of two ribosomal RNA removal methods for microbial metatranscriptomics

The predominance of rRNAs in the transcriptome is a major technical challenge in sequence-based analysis of cDNAs from microbial isolates and communities. Several approaches have been applied to deplete rRNAs from (meta)transcriptomes, but no systematic investigation of potential biases introduced by any of these approaches has been reported. Here we validated the effectiveness and fidelity of the two most commonly used approaches, subtractive hybridization and exonuclease digestion, as well as combinations of these treatments, on two synthetic five-microorganism metatranscriptomes using massively parallel sequencing. We found that the effectiveness of rRNA removal was a function of community composition and RNA integrity for these treatments. Subtractive hybridization alone introduced the least bias in relative transcript abundance, whereas exonuclease and in particular combined treatments greatly compromised mRNA abundance fidelity. Illumina sequencing itself also can compromise quantitative data analysis by introducing a G+C bias between runs.

The Single-Nucleotide Resolution Transcriptome of Pseudomonas aeruginosa Grown in Body Temperature

One of the hallmarks of opportunistic pathogens is their ability to adjust and respond to a wide range of environmental and host-associated conditions. The human pathogen Pseudomonas aeruginosa has an ability to thrive in a variety of hosts and cause a range of acute and chronic infections in individuals with impaired host defenses or cystic fibrosis. Here we report an in-depth transcriptional profiling of this organism when grown at host-related temperatures. Using RNA-seq of samples from P. aeruginosa grown at 28°C and 37°C we detected genes preferentially expressed at the body temperature of mammalian hosts, suggesting that they play a role during infection. These temperature-induced genes included the type III secretion system (T3SS) genes and effectors, as well as the genes responsible for phenazines biosynthesis. Using genome-wide transcription start site (TSS) mapping by RNA-seq we were able to accurately define the promoters and cis-acting RNA elements of many genes, and uncovered new genes and previously unrecognized non-coding RNAs directly controlled by the LasR quorum sensing regulator. Overall we identified 165 small RNAs and over 380 cis-antisense RNAs, some of which predicted to perform regulatory functions, and found that non-coding RNAs are preferentially localized in pathogenicity islands and horizontally transferred regions. Our work identifies regulatory features of P. aeruginosa genes whose products play a role in environmental adaption during infection and provides a reference transcriptional landscape for this pathogen.

Transcriptome-Wide Mapping of 5-methylcytidine RNA Modifications in Bacteria, Archaea, and Yeast Reveals m5C within Archaeal mRNAs

The presence of 5-methylcytidine (m5C) in tRNA and rRNA molecules of a wide variety of organisms was first observed more than 40 years ago. However, detection of this modification was limited to specific, abundant, RNA species, due to the usage of low-throughput methods. To obtain a high resolution, systematic, and comprehensive transcriptome-wide overview of m5C across the three domains of life, we used bisulfite treatment on total RNA from both gram positive (B. subtilis) and gram negative (E. coli) bacteria, an archaeon (S. solfataricus) and a eukaryote (S. cerevisiae), followed by massively parallel sequencing. We were able to recover most previously documented m5C sites on rRNA in the four organisms, and identified several novel sites in yeast and archaeal rRNAs. Our analyses also allowed quantification of methylated m5C positions in 64 tRNAs in yeast and archaea, revealing stoichiometric differences between the methylation patterns of these organisms. Molecules of tRNAs in which m5C was absent were also discovered. Intriguingly, we detected m5C sites within archaeal mRNAs, and identified a consensus motif of AUCGANGU that directs methylation in S. solfataricus. Our results, which were validated using m5C-specific RNA immunoprecipitation, provide the first evidence for mRNA modifications in archaea, suggesting that this mode of post-transcriptional regulation extends beyond the eukaryotic domain.

The excludon: a new concept in bacterial antisense RNA-mediated gene regulation

In recent years, non-coding RNAs have emerged as key regulators of gene expression. Among these RNAs, the antisense RNAs (asRNAs) are particularly abundant, but in most cases the function and mechanism of action for a particular asRNA remains elusive. Here, we highlight a recently discovered paradigm termed the excludon, which defines a genomic locus encoding an unusually long asRNA that spans divergent genes or operons with related or opposing functions. Because these asRNAs can inhibit the expression of one operon while functioning as an mRNA for the adjacent operon, they act as fine-tuning regulatory switches in bacteria.

A vast collection of microbial genes that are toxic to bacteria

In the process of clone-based genome sequencing, initial assemblies frequently contain cloning gaps that can be resolved using cloning-independent methods, but the reason for their occurrence is largely unknown. By analyzing 9,328,693 sequencing clones from 393 microbial genomes, we systematically mapped more than 15,000 genes residing in cloning gaps and experimentally showed that their expression products are toxic to the Escherichia coli host. A subset of these toxic sequences was further evaluated through a series of functional assays exploring the mechanisms of their toxicity. Among these genes, our assays revealed novel toxins and restriction enzymes, and new classes of small, non-coding toxic RNAs that reproducibly inhibit E. coli growth. Further analyses also revealed abundant, short, toxic DNA fragments that were predicted to suppress E. coli growth by interacting with the replication initiator DnaA. Our results show that cloning gaps, once considered the result of technical problems, actually serve as a rich source for the discovery of biotechnologically valuable functions, and suggest new modes of antimicrobial interventions.

Differential Translation Tunes Uneven Production of Operon-Encoded Proteins

Clustering of functionally related genes in operons allows for coregulated gene expression in prokaryotes. This is advantageous when equal amounts of gene products are required. Production of protein complexes with an uneven stoichiometry, however, requires tuning mechanisms to generate subunits in appropriate relative quantities. Using comparative genomic analysis, we show that differential translation is a key determinant of modulated expression of genes clustered in operons and that codon bias generally is the best in silico indicator of unequal protein production. Variable ribosome density profiles of polycistronic transcripts correlate strongly with differential translation patterns. In addition, we provide experimental evidence that de novo initiation of translation can occur at intercistronic sites, allowing for differential translation of any gene irrespective of its position on a polycistronic messenger. Thus, modulation of translation efficiency appears to be a universal mode of control in bacteria and archaea that allows for differential production of operon-encoded proteins.

Mutation detection with next-generation resequencing through a mediator genome.

The affordability of next generation sequencing (NGS) is transforming the field of mutation analysis in bacteria. The genetic basis for phenotype alteration can be identified directly by sequencing the entire genome of the mutant and comparing it to the wild-type (WT) genome, thus identifying acquired mutations. A major limitation for this approach is the need for an a-priori sequenced reference genome for the WT organism, as the short reads of most current NGS approaches usually prohibit de-novo genome assembly. To overcome this limitation we propose a general framework that utilizes the genome of relative organisms as mediators for comparing WT and mutant bacteria. Under this framework, both mutant and WT genomes are sequenced with NGS, and the short sequencing reads are mapped to the mediator genome. Variations between the mutant and the mediator that recur in the WT are ignored, thus pinpointing the differences between the mutant and the WT. To validate this approach we sequenced the genome of Bdellovibrio bacteriovorus 109J, an obligatory bacterial predator, and its prey-independent mutant, and compared both to the mediator species Bdellovibrio bacteriovorus HD100. Although the mutant and the mediator sequences differed in more than 28,000 nucleotide positions, our approach enabled pinpointing the single causative mutation. Experimental validation in 53 additional mutants further established the implicated gene. Our approach extends the applicability of NGS-based mutant analyses beyond the domain of available reference genomes.