Lars Barquist

Rfam 11.0: 10 years of RNA families

The Rfam database (available via the website at http://rfam.sanger.ac.uk and through our mirror at http://rfam.janelia.org) is a collection of non-coding RNA families, primarily RNAs with a conserved RNA secondary structure, including both RNA genes and mRNA cis-regulatory elements. Each family is represented by a multiple sequence alignment, predicted secondary structure and covariance model. Here we discuss updates to the database in the latest release, Rfam 11.0, including the introduction of genome-based alignments for large families, the introduction of the Rfam Biomart as well as other user interface improvements. Rfam is available under the Creative Commons Zero license.

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.

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.

Parallel independent evolution of pathogenicity within the genus Yersinia.

Significance Our past understanding of pathogen evolution has been fragmented because of tendencies to study human clinical isolates. To understand the evolutionary trends of pathogenic bacteria though, we need the context of their nonpathogenic relatives. Our unique and detailed dataset allows description of the parallel evolution of two key human pathogens: the causative agents of plague and Yersinia diarrhea. The analysis reveals an emerging pattern where few virulence-related functions are found in all pathogenic lineages, representing key “foothold” moments that mark the emergence of these pathogens. Functional gene loss and metabolic streamlining are equally complementing the evolution of Yersinia across the pathogenic spectrum. The genus Yersinia has been used as a model system to study pathogen evolution. Using whole-genome sequencing of all Yersinia species, we delineate the gene complement of the whole genus and define patterns of virulence evolution. Multiple distinct ecological specializations appear to have split pathogenic strains from environmental, nonpathogenic lineages. This split demonstrates that contrary to hypotheses that all pathogenic Yersinia species share a recent common pathogenic ancestor, they have evolved independently but followed parallel evolutionary paths in acquiring the same virulence determinants as well as becoming progressively more limited metabolically. Shared virulence determinants are limited to the virulence plasmid pYV and the attachment invasion locus ail. These acquisitions, together with genomic variations in metabolic pathways, have resulted in the parallel emergence of related pathogens displaying an increasingly specialized lifestyle with a spectrum of virulence potential, an emerging theme in the evolution of other important human pathogens.

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.

Patterns of genome evolution that have accompanied host adaptation in Salmonella

Significance Common features have been observed in the genome sequences of bacterial pathogens that infect few hosts. These “host adaptations” include the acquisition of pathogenicity islands of multiple genes involved in disease, losses of whole genes, and even single mutations that affect gene function. Within Salmonella enterica is a natural model system of four pathogens that are each other’s closest relatives, including a host-generalist, two host-specialists, and one with strong host associations. With whole-genome sequences, we aimed to improve our understanding of the number, nature, and order of these host adaptation events, shedding light on how human and animal pathogens arose in the past, and potentially allowing us to predict how emerging pathogens will evolve in the future. Many bacterial pathogens are specialized, infecting one or few hosts, and this is often associated with more acute disease presentation. Specific genomes show markers of this specialization, which often reflect a balance between gene acquisition and functional gene loss. Within Salmonella enterica subspecies enterica, a single lineage exists that includes human and animal pathogens adapted to cause infection in different hosts, including S. enterica serovar Enteritidis (multiple hosts), S. Gallinarum (birds), and S. Dublin (cattle). This provides an excellent evolutionary context in which differences between these pathogen genomes can be related to host range. Genome sequences were obtained from ∼60 isolates selected to represent the known diversity of this lineage. Examination and comparison of the clades within the phylogeny of this lineage revealed signs of host restriction as well as evolutionary events that mark a path to host generalism. We have identified the nature and order of events for both evolutionary trajectories. The impact of functional gene loss was predicted based upon position within metabolic pathways and confirmed with phenotyping assays. The structure of S. Enteritidis is more complex than previously known, as a second clade of S. Enteritidis was revealed that is distinct from those commonly seen to cause disease in humans or animals, and that is more closely related to S. Gallinarum. Isolates from this second clade were tested in a chick model of infection and exhibited a reduced colonization phenotype, which we postulate represents an intermediate stage in pathogen–host adaptation.

A High-Resolution View of Genome-Wide Pneumococcal Transformation

Transformation is an important mechanism of microbial evolution through which bacteria have been observed to rapidly adapt in response to clinical interventions; examples include facilitating vaccine evasion and the development of penicillin resistance in the major respiratory pathogen Streptococcus pneumoniae. To characterise the process in detail, the genomes of 124 S. pneumoniae isolates produced through in vitro transformation were sequenced and recombination events detected. Those recombinations importing the selected marker were independent of unselected events elsewhere in the genome, the positions of which were not significantly affected by local sequence similarity between donor and recipient or mismatch repair processes. However, both types of recombinations were sometimes mosaic, with multiple non-contiguous segments originating from the same molecule of donor DNA. The lengths of the unselected events were exponentially distributed with a mean of 2.3 kb, implying that recombinations are stochastically resolved with a fixed per base probability of 4.4×10(-4) bp(-1). This distribution of recombination sizes, coupled with an observed under representation of large insertions within transferred sequence, suggests transformation has the potential to reduce the size of bacterial genomes, and is unlikely to act as an efficient mechanism for the uptake of accessory genomic loci.

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.

The agr locus regulates virulence and colonization genes in Clostridium difficile 027.

The transcriptional regulator AgrA, a member of the LytTR family of proteins, plays a key role in controlling gene expression in some Gram-positive pathogens, including Staphylococcus aureus and Enterococcus faecalis. AgrA is encoded by the agrACDB global regulatory locus, and orthologues are found within the genome of most Clostridium difficile isolates, including the epidemic lineage 027/BI/NAP1. Comparative RNA sequencing of the wild type and otherwise isogenic agrA null mutant derivatives of C. difficile R20291 revealed a network of approximately 75 differentially regulated transcripts at late exponential growth phase, including many genes associated with flagellar assembly and function, such as the major structural subunit, FliC. Other differentially regulated genes include several involved in bis-(3′-5′)-cyclic dimeric GMP (c-di-GMP) synthesis and toxin A expression. C. difficile 027 R20291 agrA mutant derivatives were poorly flagellated and exhibited reduced levels of colonization and relapses in the murine infection model. Thus, the agr locus likely plays a contributory role in the fitness and virulence potential of C. difficile strains in the 027/BI/NAP1 lineage.

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.