Genome-wide identification of regulatory RNAs in the human pathogen Clostridium difficile

Identification of sRNAs expressed by the human pathogen Neisseria gonorrhoeae under disparate growth conditions

In the last several years, bacterial gene regulation via small RNAs (sRNAs) has been recognized as an important mechanism controlling expression of essential proteins that are critical to bacterial growth and metabolism. Technologies such as RNA-seq are rapidly expanding the field of sRNAs and are enabling a global view of the “sRNAome” of several bacterial species. While numerous sRNAs have been identified in a variety of both Gram-negative and Gram-positive bacteria, only a very small number have been fully characterized in the human pathogen Neisseria gonorrhoeae, the etiological agent of the STD gonorrhea. Here we present the first analysis of N. gonorrhoeae specifically focused on the identification of sRNAs through RNA-seq analysis of the organism cultured under different in vitro growth conditions. Using a new computational program, Rockhopper, to analyze prokaryotic RNA-seq data obtained from N. gonorrhoeae we identified several putative sRNAs and confirmed their expression and size through Northern blot analysis. In addition, RNA was collected from four different growth conditions (iron replete and deplete, as well as with and without co-culture with human endocervical cells). Many of the putative sRNAs identified shoed varying expression levels relative to the different growth conditions examine or were detected only under certain conditions but not others. Comparisons of identified sRNAs with the regulatory pattern of putative mRNA targets revealed possible functional roles for these sRNAs. These studies are the first to carry out a global analysis of N. gonorrhoeae specifically focused on sRNAs and show that RNA-mediated regulation may be an important mechanism of gene control in this human pathogen.

Abundant and diverse clustered regularly interspaced short palindromic repeat spacers in Clostridium difficile strains and prophages target multiple phage types within this pathogen

Clostridium difficile is an important human-pathogenic bacterium causing antibiotic-associated nosocomial infections worldwide. Mobile genetic elements and bacteriophages have helped shape C. difficile genome evolution. In many bacteria, phage infection may be controlled by a form of bacterial immunity called the clustered regularly interspaced short palindromic repeats/CRISPR-associated (CRISPR/Cas) system. This uses acquired short nucleotide sequences (spacers) to target homologous sequences (protospacers) in phage genomes. C. difficile carries multiple CRISPR arrays, and in this paper we examine the relationships between the host- and phage-carried elements of the system. We detected multiple matches between spacers and regions in 31 C. difficile phage and prophage genomes. A subset of the spacers was located in prophage-carried CRISPR arrays. The CRISPR spacer profiles generated suggest that related phages would have similar host ranges. Furthermore, we show that C. difficile strains of the same ribotype could either have similar or divergent CRISPR contents. Both synonymous and nonsynonymous mutations in the protospacer sequences were identified, as well as differences in the protospacer adjacent motif (PAM), which could explain how phages escape this system. This paper illustrates how the distribution and diversity of CRISPR spacers in C. difficile, and its prophages, could modulate phage predation for this pathogen and impact upon its evolution and pathogenicity.

IMPORTANCE

Clostridium difficile is a significant bacterial human pathogen which undergoes continual genome evolution, resulting in the emergence of new virulent strains. Phages are major facilitators of genome evolution in other bacterial species, and we use sequence analysis-based approaches in order to examine whether the CRISPR/Cas system could control these interactions across divergent C. difficile strains. The presence of spacer sequences in prophages that are homologous to phage genomes raises an extra level of complexity in this predator-prey microbial system. Our results demonstrate that the impact of phage infection in this system is widespread and that the CRISPR/Cas system is likely to be an important aspect of the evolutionary dynamics in C. difficile.

Primary transcriptome map of the hyperthermophilic archaeon Thermococcus kodakarensis

Background

Prokaryotes have relatively small genomes, densely-packed with protein-encoding sequences. RNA sequencing has, however, revealed surprisingly complex transcriptomes and here we report the transcripts present in the model hyperthermophilic Archaeon, Thermococcus kodakarensis, under different physiological conditions.

Results

Sequencing cDNA libraries, generated from RNA isolated from cells under different growth and metabolic conditions has identified >2,700 sites of transcription initiation, established a genome-wide map of transcripts, and consensus sequences for transcription initiation and post-transcription regulatory elements. The primary transcription start sites (TSS) upstream of 1,254 annotated genes, plus 644 primary TSS and their promoters within genes, are identified. Most mRNAs have a 5'-untranslated region (5'-UTR) 10 to 50 nt long (median = 16 nt), but ~20% have 5'-UTRs from 50 to 300 nt long and ~14% are leaderless. Approximately 50% of mRNAs contain a consensus ribosome binding sequence. The results identify TSS for 1,018 antisense transcripts, most with sequences complementary to either the 5'- or 3'-region of a sense mRNA, and confirm the presence of transcripts from all three CRISPR loci, the RNase P and 7S RNAs, all tRNAs and rRNAs and 69 predicted snoRNAs. Two putative riboswitch RNAs were present in growing but not in stationary phase cells. The procedure used is designed to identify TSS but, assuming that the number of cDNA reads correlates with transcript abundance, the results also provide a semi-quantitative documentation of the differences in T. kodakarensis genome expression under different growth conditions and confirm previous observations of substrate-dependent specific gene expression. Many previously unanticipated small RNAs have been identified, some with relative low GC contents (≤50%) and sequences that do not fold readily into base-paired secondary structures, contrary to the classical expectations for non-coding RNAs in a hyperthermophile.

Conclusion

The results identify >2,700 TSS, including almost all of the primary sites of transcription initiation upstream of annotated genes, plus many secondary sites, sites within genes and sites resulting in antisense transcripts. The T. kodakarensis genome is small (~2.1 Mbp) and tightly packed with protein-encoding genes, but the transcriptomes established also contain many non-coding RNAs and predict extensive RNA-based regulation in this model Archaeon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-684) contains supplementary material, which is available to authorized users.

Small RNA functions in carbon metabolism and virulence of enteric pathogens

Enteric pathogens often cycle between virulent and saprophytic lifestyles. To endure these frequent changes in nutrient availability and composition bacteria possess an arsenal of regulatory and metabolic genes allowing rapid adaptation and high flexibility. While numerous proteins have been characterized with regard to metabolic control in pathogenic bacteria, small non-coding RNAs have emerged as additional regulators of metabolism. Recent advances in sequencing technology have vastly increased the number of candidate regulatory RNAs and several of them have been found to act at the interface of bacterial metabolism and virulence factor expression. Importantly, studying these riboregulators has not only provided insight into their metabolic control functions but also revealed new mechanisms of post-transcriptional gene control. This review will focus on the recent advances in this area of host-microbe interaction and discuss how regulatory small RNAs may help coordinate metabolism and virulence of enteric pathogens.

Antisense transcription regulates the expression of the enterohemorrhagic Escherichia coli virulence regulatory gene ler in response to the intracellular iron concentration

Enteric pathogens, such as enterohemorrhagic E. coli (EHEC) O157:H7, encounter varying concentrations of iron during their life cycle. In the gastrointestinal tract, the amount of available free iron is limited because of absorption by host factors. EHEC and other enteric pathogens have developed sophisticated iron-responsive systems to utilize limited iron resources, and these systems are primarily regulated by the Fur repressor protein. The iron concentration could be a signal that controls gene expression in the intestines. In this study, we explored the role of iron in LEE (locus for enterocyte effacement) virulence gene expression in EHEC. In contrast to the expression of Fur-regulated genes, the expression of LEE genes was greatly reduced in fur mutants irrespective of the iron concentration. The expression of the ler gene, the LEE-encoded master regulator, was affected at a post-transcription step by fur mutation. Further analysis showed that the loss of Fur affected the translation of the ler gene by increasing the intracellular concentration of free iron, and the transcription of the antisense strand was necessary for regulation. The results indicate that LEE gene expression is closely linked to the control of intracellular free iron homeostasis.

Bacteriophage behavioral ecology: How phages alter their bacterial host's habits

Bacteriophages have an essential gene kit that enables their invasion, replication, and production. In addition to this "core" genome, they can carry "accessory" genes that dramatically impact bacterial biology, and presumably boost their own success. The content of phage genomes continue to surprise us by revealing new ways that viruses impact bacterial biology. The genome of a Clostridium difficile myovirus, phiCDHM1, contains homologs of three bacterial accessory gene regulator (agr) genes. The agr system is a type of quorum sensing (QS), via which the phage may modify C. difficile interactions with its environment. Although their mechanism of action is unknown, mutants in bacterial versions of these genes impact sporulation and virulence. To explore how phage QS genes may influence C. difficile biology, we examine the main categories of bacterial behavior that phages have been shown to influence and discuss how interactions via QS could influence behavior at a wider level.

Pleiotropic role of the RNA chaperone protein Hfq in the human pathogen Clostridium difficile

Clostridium difficile is an emergent human pathogen and the most common cause of nosocomial diarrhea. Our recent data strongly suggest the importance of RNA-based mechanisms for the control of gene expression in C. difficile. In an effort to understand the function of the RNA chaperone protein Hfq, we constructed and characterized an Hfq-depleted strain in C. difficile. Hfq depletion led to a growth defect, morphological changes, an increased sensitivity to stresses, and a better ability to sporulate and to form biofilms. The transcriptome analysis revealed pleiotropic effects of Hfq depletion on gene expression in C. difficile, including genes encoding proteins involved in sporulation, stress response, metabolic pathways, cell wall-associated proteins, transporters, and transcriptional regulators and genes of unknown function. Remarkably, a great number of genes of the regulon dependent on sporulation-specific sigma factor, SigK, were upregulated in the Hfq-depleted strain. The altered accumulation of several sRNAs and interaction of Hfq with selected sRNAs suggest potential involvement of Hfq in these regulatory RNA functions. Altogether, these results suggest the pleiotropic role of Hfq protein in C. difficile physiology, including processes important for the C. difficile infection cycle, and expand our knowledge of Hfq-dependent regulation in Gram-positive bacteria.

A novel secreted metalloprotease (CD2830) from Clostridium difficile cleaves specific proline sequences in LPXTG cell surface proteins

Bacterial secreted proteins constitute a biologically important subset of proteins involved in key processes related to infection such as adhesion, colonization, and dissemination. Bacterial extracellular proteases, in particular, have attracted considerable attention, as they have been shown to be indispensable for bacterial virulence. Here, we analyzed the extracellular subproteome of Clostridium difficile and identified a hypothetical protein, CD2830, as a novel secreted metalloprotease. Following the identification of a CD2830 cleavage site in human HSP90β, a series of synthetic peptide substrates was used to identify the favorable CD2830 cleavage motif. This motif was characterized by a high prevalence of proline residues. Intriguingly, CD2830 has a preference for cleaving Pro-Pro bonds, unique among all hitherto described proteases. Strikingly, within the C. difficile proteome two putative adhesion molecules, CD2831 and CD3246, were identified that contain multiple CD2830 cleavage sites (13 in total). We subsequently found that CD2830 efficiently cleaves CD2831 between two prolines at all predicted cleavage sites. Moreover, native CD2830, secreted by live cells, cleaves endogenous CD2831 and CD3246. These findings highlight CD2830 as a highly specific endoproteinase with a preference for proline residues surrounding the scissile bond. Moreover, the efficient cleavage of two putative surface adhesion proteins points to a possible role of CD2830 in the regulation of C. difficile adhesion.

When a virus is not a parasite: the beneficial effects of prophages on bacterial fitness

Most organisms on the planet have viruses that infect them. Viral infection may lead to cell death, or to a symbiotic relationship where the genomes of both virus and host replicate together. In the symbiotic state, both virus and cell potentially experience increased fitness as a result of the other. The viruses that infect bacteria, called bacteriophages (or phages), well exemplify the symbiotic relationships that can develop between viruses and their host. In this review, we will discuss the many ways that prophages, which are phage genomes integrated into the genomes of their hosts, influence bacterial behavior and virulence.

To acquire or resist: the complex biological effects of CRISPR-Cas systems

Prokaryotic CRISPR-Cas (clustered regularly interspaced short palindromic repeat-CRISPR associated) systems provide a sophisticated adaptive immune system that offers protection against foreign DNA. These systems are widely distributed in prokaryotes and exert an important influence on bacterial behavior and evolution. However, interpreting the biological effects of a CRISPR-Cas system within a given species can be complicated because the outcome of rejecting foreign DNA does not always provide a fitness advantage, as foreign DNA uptake is sometimes beneficial. To address these issues, here we review data pertaining to the potential in vivo costs and benefits of CRISPR-Cas systems, novel functions for these systems, and how they may be inactivated.

Semiquantitative analysis of clinical heat stress in Clostridium difficile strain 630 using a GeLC/MS workflow with emPAI quantitation

Clostridium difficile is considered to be the most frequent cause of infectious bacterial diarrhoea in hospitals worldwide yet its adaptive ability remains relatively uncharacterised. Here, we used GeLC/MS and the exponentially modified protein abundance index (emPAI) calculation to determine proteomic changes in response to a clinically relevant heat stress. Reproducibility between both biological and technical replicates was good, and a 37°C proteome of 224 proteins was complemented by a 41°C proteome of 202 proteins at a 1% false discovery rate. Overall, 236 C. difficile proteins were identified and functionally categorised, of which 178 were available for comparative purposes. A total of 65 proteins (37%) were modulated by 1.5-fold or more at 41°C compared to 37°C and we noted changes in the majority of proteins associated with amino acid metabolism, including upregulation of the reductive branch of the leucine fermentation pathway. Motility was reduced at 41°C as evidenced by a 2.7 fold decrease in the flagellar filament protein, FliC, and a global increase in proteins associated with detoxification and adaptation to atypical conditions was observed, concomitant with decreases in proteins mediating transcriptional elongation and the initiation of protein synthesis. Trigger factor was down regulated by almost 5-fold. We propose that under heat stress, titration of the GroESL and dnaJK/grpE chaperones by misfolded proteins will, in the absence of trigger factor, prevent nascent chains from emerging efficiently from the ribosome causing translational stalling and also an increase in secretion. The current work has thus allowed development of a heat stress model for the key cellular processes of protein folding and export.

Genome-wide analysis of small RNAs expressed by Yersinia pestis identifies a regulator of the Yop-Ysc type III secretion system

Small noncoding RNA (sRNA) molecules are integral components of the regulatory machinery for many bacterial species and are known to posttranscriptionally regulate metabolic and stress-response pathways, quorum sensing, virulence factors, and more. The Yop-Ysc type III secretion system (T3SS) is a critical virulence component for the pathogenic Yersinia species, and the regulation of this system is tightly controlled at each step from transcription to translocation of effectors into host cells. The contribution of sRNAs to the regulation of the T3SS in Yersinia has been largely unstudied, however. Previously, our lab identified a role for the sRNA chaperone protein Hfq in the regulation of components of the T3SS in the gastrointestinal pathogen Yersinia pseudotuberculosis. Here we present data demonstrating a similar requirement for Hfq in the closely related species Yersinia pestis. Through deep sequencing analysis of the Y. pestis sRNA-ome, we found 63 previously unidentified putative sRNAs in this species. We identified a Yersinia-specific sRNA, Ysr141, carried by the T3SS plasmid pCD1 that is required for the production of multiple T3SS proteins. In addition, we show that Ysr141 targets an untranslated region upstream of yopJ to posttranscriptionally activate the synthesis of the YopJ protein. Furthermore, Ysr141 may be an unstable and/or processed sRNA, which could contribute to its function in the regulation of the T3SS. The discovery of an sRNA that influences the synthesis of the T3SS adds an additional layer of regulation to this tightly controlled virulence determinant of Y. pestis.

Small regulatory RNAs from low-GC Gram-positive bacteria

Small regulatory RNAs (sRNAs) that act by base-pairing were first discovered in so-called accessory DNA elements—plasmids, phages, and transposons—where they control replication, maintenance, and transposition. Since 2001, a huge body of work has been performed to predict and identify sRNAs in a multitude of bacterial genomes. The majority of chromosome-encoded sRNAs have been investigated in E. coli and other Gram-negative bacteria. However, during the past five years an increasing number of sRNAs were found in Gram-positive bacteria. Here, we outline our current knowledge on chromosome-encoded sRNAs from low-GC Gram-positive species that act by base-pairing, i.e., an antisense mechanism. We will focus on sRNAs with known targets and defined regulatory mechanisms with special emphasis on Bacillus subtilis.

Posttranscriptional regulation of oral bacterial adaptive responses

Within the past 10 years, it has become increasingly evident that posttranscriptional regulation is among the most important mechanisms used by bacteria to modulate gene expression in response to environmental perturbations. Posttranscriptional mechanisms provide a much faster response and lower energy burden compared to most transcription regulatory pathways and they have the unique advantage that they can override existing transcriptional responses once the environment changes. Because of this, virulence factor gene expression is particularly suited for posttranscriptional control, and not surprisingly, an abundance of recent evidence indicates that posttranscriptional regulators are the predominant virulence regulators of human pathogens. Typically, this involves global riboregulators that primarily serve as modulators of virulence gene translation initiation and/or mRNA stability. Surprisingly little has been reported about posttranscriptional regulatory pathways in oral bacteria, but recent results suggest that oral species are equally dependent upon posttranscriptional control of their adaptive genetic responses. In this report, we discuss the major themes in RNA-based regulation of gene expression and review the available literature related to the most commonly studied oral bacterial species.

The Clostridium small RNome that responds to stress: the paradigm and importance of toxic metabolite stress in C. acetobutylicum

Background

Small non-coding RNAs (sRNA) are emerging as major components of the cell’s regulatory network, several possessing their own regulons. A few sRNAs have been reported as being involved in general or toxic-metabolite stress, mostly in Gram^- prokaryotes, but hardly any in Gram⁺ prokaryotes. Significantly, the role of sRNAs in the stress response remains poorly understood at the genome-scale level. It was previously shown that toxic-metabolite stress is one of the most comprehensive and encompassing stress responses in the cell, engaging both the general stress (or heat-shock protein, HSP) response as well as specialized metabolic programs.

Results

Using RNA deep sequencing (RNA-seq) we examined the sRNome of C. acetobutylicum in response to the native but toxic metabolites, butanol and butyrate. 7.5% of the RNA-seq reads mapped to genome outside annotated ORFs, thus demonstrating the richness and importance of the small RNome. We used comparative expression analysis of 113 sRNAs we had previously computationally predicted, and of annotated mRNAs to set metrics for reliably identifying sRNAs from RNA-seq data, thus discovering 46 additional sRNAs. Under metabolite stress, these 159 sRNAs displayed distinct expression patterns, a select number of which was verified by Northern analysis. We identified stress-related expression of sRNAs affecting transcriptional (6S, S-box & solB) and translational (tmRNA & SRP-RNA) processes, and 65 likely targets of the RNA chaperone Hfq.

Conclusions

Our results support an important role for sRNAs for understanding the complexity of the regulatory network that underlies the stress response in Clostridium organisms, whether related to normophysiology, pathogenesis or biotechnological applications.

Genome-wide mapping of transcriptional start sites defines an extensive leaderless transcriptome in Mycobacterium tuberculosis

Deciphering physiological changes that mediate transition of Mycobacterium tuberculosis between replicating and nonreplicating states is essential to understanding how the pathogen can persist in an individual host for decades. We have combined RNA sequencing (RNA-seq) of 5′ triphosphate-enriched libraries with regular RNA-seq to characterize the architecture and expression of M. tuberculosis promoters. We identified over 4,000 transcriptional start sites (TSSs). Strikingly, for 26% of the genes with a primary TSS, the site of transcriptional initiation overlapped with the annotated start codon, generating leaderless transcripts lacking a 5′ UTR and, hence, the Shine-Dalgarno sequence commonly used to initiate ribosomal engagement in eubacteria. Genes encoding proteins with active growth functions were markedly depleted from the leaderless transcriptome, and there was a significant increase in the overall representation of leaderless mRNAs in a starvation model of growth arrest. The high percentage of leaderless genes may have particular importance in the physiology of nonreplicating M. tuberculosis.

•

A resource for the identification of in vitro active promoters in M. tuberculosis

•

A quarter of all genes in M. tuberculosis are expressed as leaderless mRNAs

•

Leaderless mRNAs are differentially associated with toxin-antitoxin modules

•

Abundance of leaderless mRNAs increases during starvation-induced growth arrest

Differential stress transcriptome landscape of historic and recently emerged hypervirulent strains of Clostridium difficile strains determined using RNA-seq

C. difficile is the most common cause of nosocomial diarrhea in North America and Europe. Genomes of individual strains of C. difficile are highly divergent. To determine how divergent strains respond to environmental changes, the transcriptomes of two historic and two recently isolated hypervirulent strains were analyzed following nutrient shift and osmotic shock. Illumina based RNA-seq was used to sequence these transcriptomes. Our results reveal that although C. difficile strains contain a large number of shared and strain specific genes, the majority of the differentially expressed genes were core genes. We also detected a number of transcriptionally active regions that were not part of the primary genome annotation. Some of these are likely to be small regulatory RNAs.

Genome-wide analysis of cell type-specific gene transcription during spore formation in Clostridium difficile

Clostridium difficile, a Gram positive, anaerobic, spore-forming bacterium is an emergent pathogen and the most common cause of nosocomial diarrhea. Although transmission of C. difficile is mediated by contamination of the gut by spores, the regulatory cascade controlling spore formation remains poorly characterized. During Bacillus subtilis sporulation, a cascade of four sigma factors, σ(F) and σ(G) in the forespore and σ(E) and σ(K) in the mother cell governs compartment-specific gene expression. In this work, we combined genome wide transcriptional analyses and promoter mapping to define the C. difficile σ(F), σ(E), σ(G) and σ(K) regulons. We identified about 225 genes under the control of these sigma factors: 25 in the σ(F) regulon, 97 σ(E)-dependent genes, 50 σ(G)-governed genes and 56 genes under σ(K) control. A significant fraction of genes in each regulon is of unknown function but new candidates for spore coat proteins could be proposed as being synthesized under σ(E) or σ(K) control and detected in a previously published spore proteome. SpoIIID of C. difficile also plays a pivotal role in the mother cell line of expression repressing the transcription of many members of the σ(E) regulon and activating sigK expression. Global analysis of developmental gene expression under the control of these sigma factors revealed deviations from the B. subtilis model regarding the communication between mother cell and forespore in C. difficile. We showed that the expression of the σ(E) regulon in the mother cell was not strictly under the control of σ(F) despite the fact that the forespore product SpoIIR was required for the processing of pro-σ(E). In addition, the σ(K) regulon was not controlled by σ(G) in C. difficile in agreement with the lack of pro-σ(K) processing. This work is one key step to obtain new insights about the diversity and evolution of the sporulation process among Firmicutes.

The second messenger cyclic Di-GMP regulates Clostridium difficile toxin production by controlling expression of sigD

The Gram-positive obligate anaerobe Clostridium difficile causes potentially fatal intestinal diseases. How this organism regulates virulence gene expression is poorly understood. In many bacterial species, the second messenger cyclic di-GMP (c-di-GMP) negatively regulates flagellar motility and, in some cases, virulence. c-di-GMP was previously shown to repress motility of C. difficile. Recent evidence indicates that flagellar gene expression is tightly linked with expression of the genes encoding the two C. difficile toxins TcdA and TcdB, which are key virulence factors for this pathogen. Here, the effect of c-di-GMP on expression of the toxin genes tcdA and tcdB was determined, and the mechanism connecting flagellar and toxin gene expressions was examined. In C. difficile, increasing c-di-GMP levels reduced the expression levels of tcdA and tcdB, as well as that of tcdR, which encodes an alternative sigma factor that activates tcdA and tcdB expression. We hypothesized that the C. difficile orthologue of the flagellar alternative sigma factor SigD (FliA; σ(28)) mediates regulation of toxin gene expression in response to c-di-GMP. Indeed, ectopic expression of sigD in C. difficile resulted in increased expression levels of tcdR, tcdA, and tcdB. Furthermore, sigD expression enhanced toxin production and increased the cytopathic effect of C. difficile on cultured fibroblasts. Finally, evidence is provided that SigD directly activates tcdR expression and that SigD cannot activate tcdA or tcdB expression independent of TcdR. Taken together, these data suggest that SigD positively regulates toxin genes in C. difficile and that c-di-GMP can inhibit both motility and toxin production via SigD, making this signaling molecule a key virulence gene regulator in C. difficile.