Bacteroides thetaiotaomicron uses a widespread extracellular DNase to promote bile-dependent biofilm formation.

Bacteroides thetaiotaomicron is a gut symbiont that inhabits the mucus layer and adheres to and metabolizes food particles, contributing to gut physiology and maturation. Although adhesion and biofilm formation could be key features for B. thetaiotaomicron stress resistance and gut colonization, little is known about the determinants of B. thetaiotaomicron biofilm formation. We previously showed that the B. thetaiotaomicron reference strain VPI-5482 is a poor in vitro biofilm former. Here, we demonstrated that bile, a gut-relevant environmental cue, triggers the formation of biofilm in many B. thetaiotaomicron isolates and common gut Bacteroidales species. We determined that bile-dependent biofilm formation involves the production of the DNase BT3563 or its homologs, degrading extracellular DNA (eDNA) in several B. thetaiotaomicron strains. Our study therefore shows that, although biofilm matrix eDNA provides a biofilm-promoting scaffold in many studied Firmicutes and Proteobacteria, BT3563-mediated eDNA degradation is required to form B. thetaiotaomicron biofilm in the presence of bile.

The CovR regulatory network drives the evolution of Group B Streptococcus virulence.

Virulence of the neonatal pathogen Group B Streptococcus is under the control of the master regulator CovR. Inactivation of CovR is associated with large-scale transcriptome remodeling and impairs almost every step of the interaction between the pathogen and the host. However, transcriptome analyses suggested a plasticity of the CovR signaling pathway in clinical isolates leading to phenotypic heterogeneity in the bacterial population. In this study, we characterized the CovR regulatory network in a strain representative of the CC-17 hypervirulent lineage responsible of the majority of neonatal meningitis. Transcriptome and genome-wide binding analysis reveal the architecture of the CovR network characterized by the direct repression of a large array of virulence-associated genes and the extent of co-regulation at specific loci. Comparative functional analysis of the signaling network links strain-specificities to the regulation of the pan-genome, including the two specific hypervirulent adhesins and horizontally acquired genes, to mutations in CovR-regulated promoters, and to variability in CovR activation by phosphorylation. This regulatory adaptation occurs at the level of genes, promoters, and of CovR itself, and allows to globally reshape the expression of virulence genes. Overall, our results reveal the direct, coordinated, and strain-specific regulation of virulence genes by the master regulator CovR and suggest that the intra-species evolution of the signaling network is as important as the expression of specific virulence factors in the emergence of clone associated with specific diseases.

The oxidative stress response of pathogenic Leptospira is controlled by two peroxide stress regulators which putatively cooperate in controlling virulence.

Pathogenic Leptospira are the causative agents of leptospirosis, the most widespread zoonotic infectious disease. Leptospirosis is a potentially severe and life-threatening emerging disease with highest burden in sub-tropical areas and impoverished populations. Mechanisms allowing pathogenic Leptospira to survive inside a host and induce acute leptospirosis are not fully understood. The ability to resist deadly oxidants produced by the host during infection is pivotal for Leptospira virulence. We have previously shown that genes encoding defenses against oxidants in L. interrogans are repressed by PerRA (encoded by LIMLP_10155), a peroxide stress regulator of the Fur family. In this study, we describe the identification and characterization of another putative PerR-like regulator (LIMLP_05620) in L. interrogans. Protein sequence and phylogenetic analyses indicated that LIMLP_05620 displayed all the canonical PerR amino acid residues and is restricted to pathogenic Leptospira clades. We therefore named this PerR-like regulator PerRB. In L. interrogans, the PerRB regulon is distinct from that of PerRA. While a perRA mutant had a greater tolerance to peroxide, inactivating perRB led to a higher tolerance to superoxide, suggesting that these two regulators have a distinct function in the adaptation of L. interrogans to oxidative stress. The concomitant inactivation of perRA and perRB resulted in a higher tolerance to both peroxide and superoxide and, unlike the single mutants, a double perRAperRB mutant was avirulent. Interestingly, this correlated with major changes in gene and non-coding RNA expression. Notably, several virulence-associated genes (clpB, ligA/B, and lvrAB) were repressed. By obtaining a double mutant in a pathogenic Leptospira strain, our study has uncovered an interplay of two PerRs in the adaptation of Leptospira to oxidative stress with a putative role in virulence and pathogenicity, most likely through the transcriptional control of a complex regulatory network.

Role of the major determinant of polar flagellation FlhG in the endoflagella-containing spirochete Leptospira.

Spirochetes can be distinguished from other bacteria by their spiral-shaped morphology and subpolar periplasmic flagella. This study focused on FlhF and FlhG, which control the spatial and numerical regulation of flagella in many exoflagellated bacteria, in the spirochete Leptospira. In contrast to flhF which seems to be essential in Leptospira, we demonstrated that flhG- mutants in both the saprophyte L. biflexa and the pathogen L. interrogans were less motile than the wild-type strains in gel-like environments but not hyperflagellated as reported previously in other bacteria. Cryo-electron tomography revealed that the distance between the flagellar basal body and the tip of the cell decreased significantly in the flhG- mutant in comparison to wild-type and complemented strains. Additionally, comparative transcriptome analyses of L. biflexa flhG- and wild-type strains showed that FlhG acts as a negative regulator of transcription of some flagellar genes. We found that the L. interrogans flhG- mutant was attenuated for virulence in the hamster model. Cross-species complementation also showed that flhG is not interchangeable between species. Our results indicate that FlhF and FlhG in Leptospira contribute to governing cell motility but our data support the hypothesis that FlhF and FlhG function differently in each bacterial species, including among spirochetes.

The transcriptional response of pathogenic Leptospira to peroxide reveals new defenses against infection-related oxidative stress.

Pathogenic Leptospira spp. are the causative agents of the waterborne zoonotic disease leptospirosis. Leptospira are challenged by numerous adverse conditions, including deadly reactive oxygen species (ROS), when infecting their hosts. Withstanding ROS produced by the host innate immunity is an important strategy evolved by pathogenic Leptospira for persisting in and colonizing hosts. In L. interrogans, genes encoding defenses against ROS are repressed by the peroxide stress regulator, PerR. In this study, RNA sequencing was performed to characterize both the L. interrogans response to low and high concentrations of hydrogen peroxide and the PerR regulon. We showed that Leptospira solicit three main peroxidase machineries (catalase, cytochrome C peroxidase and peroxiredoxin) and heme to detoxify oxidants produced during peroxide stress. In addition, canonical molecular chaperones of the heat shock response and DNA repair proteins from the SOS response were required for Leptospira recovering from oxidative damage. Identification of the PerR regulon upon exposure to H2O2 allowed to define the contribution of this regulator in the oxidative stress response. This study has revealed a PerR-independent regulatory network involving other transcriptional regulators, two-component systems and sigma factors as well as non-coding RNAs that putatively orchestrate, in concert with PerR, the oxidative stress response. We have shown that PerR-regulated genes encoding a TonB-dependent transporter and a two-component system (VicKR) are involved in Leptospira tolerance to superoxide. This could represent the first defense mechanism against superoxide in L. interrogans, a bacterium lacking canonical superoxide dismutase. Our findings provide an insight into the mechanisms required by pathogenic Leptospira to overcome oxidative damage during infection-related conditions. This will participate in framing future hypothesis-driven studies to identify and decipher novel virulence mechanisms in this life-threatening pathogen.

Mosquito-bacteria interactions during larval development trigger metabolic changes with carry-over effects on adult fitness.

In animals with distinct life stages such as holometabolous insects, adult phenotypic variation is often shaped by the environment of immature stages, including their interactions with microbes colonizing larval habitats. Such carry-over effects were previously observed for several adult traits of the mosquito Aedes aegypti after larval exposure to different bacteria, but the mechanistic underpinnings are unknown. Here, we investigated the molecular changes triggered by gnotobiotic larval exposure to different bacteria in Ae. aegypti. We initially screened a panel of 16 bacterial isolates from natural mosquito breeding sites to determine their ability to influence adult life-history traits. We subsequently focused on four bacterial isolates (belonging to Flavobacterium, Lysobacter, Paenibacillus, and Enterobacteriaceae) with significant carry-over effects on adult survival and found that they were associated with distinct transcriptomic profiles throughout mosquito development. Moreover, we detected carry-over effects at the level of gene expression for the Flavobacterium and Paenibacillus isolates. The most prominent transcriptomic changes in gnotobiotic larvae reflected a profound remodelling of lipid metabolism, which translated into phenotypic differences in lipid storage and starvation resistance at the adult stage. Together, our findings indicate that larval exposure to environmental bacteria trigger substantial physiological changes that impact adult fitness, uncovering a possible mechanism underlying carry-over effects of mosquito-bacteria interactions during larval development.

4-Methylcytosine DNA modification is critical for global epigenetic regulation and virulence in the human pathogen Leptospira interrogans.

In bacteria, DNA methylation can be facilitated by 'orphan' DNA methyltransferases lacking cognate restriction endonucleases, but whether and how these enzymes control key cellular processes are poorly understood. The effects of a specific modification, 4-methylcytosine (4mC), are even less clear, as this epigenetic marker is unique to bacteria and archaea, whereas the bulk of epigenetic research is currently performed on eukaryotes. Here, we characterize a 4mC methyltransferase from the understudied pathogen Leptospira spp. Inactivating this enzyme resulted in complete abrogation of CTAG motif methylation, leading to genome-wide dysregulation of gene expression. Mutants exhibited growth defects, decreased adhesion to host cells, higher susceptibility to LPS-targeting antibiotics, and, importantly, were no longer virulent in an acute infection model. Further investigation resulted in the discovery of at least one gene, that of an ECF sigma factor, whose transcription was altered in the methylase mutant and, subsequently, by mutation of the CTAG motifs in the promoter of the gene. The genes that comprise the regulon of this sigma factor were, accordingly, dysregulated in the methylase mutant and in a strain overexpressing the sigma factor. Our results highlight the importance of 4mC in Leptospira physiology, and suggest the same of other understudied species.

Mycobacterium abscessus virulence traits unraveled by transcriptomic profiling in amoeba and macrophages.

Free-living amoebae are thought to represent an environmental niche in which amoeba-resistant bacteria may evolve towards pathogenicity. To get more insights into factors playing a role for adaptation to intracellular life, we characterized the transcriptomic activities of the emerging pathogen Mycobacterium abscessus in amoeba and murine macrophages (Mϕ) and compared them with the intra-amoebal transcriptome of the closely related, but less pathogenic Mycobacterium chelonae. Data on up-regulated genes in amoeba point to proteins that allow M. abscessus to resist environmental stress and induce defense mechanisms, as well as showing a switch from carbohydrate carbon sources to fatty acid metabolism. For eleven of the most upregulated genes in amoeba and/or Mϕ, we generated individual gene knock-out M. abscessus mutant strains, from which ten were found to be attenuated in amoeba and/or Mϕ in subsequence virulence analyses. Moreover, transfer of two of these genes into the genome of M. chelonae increased the intra-Mϕ survival of the recombinant strain. One knock-out mutant that had the gene encoding Eis N-acetyl transferase protein (MAB_4532c) deleted, was particularly strongly attenuated in Mϕ. Taken together, M. abscessus intra-amoeba and intra-Mϕ transcriptomes revealed the capacity of M. abscessus to adapt to an intracellular lifestyle, with amoeba largely contributing to the enhancement of M. abscessus intra-Mϕ survival.

Insights into amebiasis using a human 3D-intestinal model.

Entamoeba histolytica is the causative agent of amebiasis, an infectious disease targeting the intestine and the liver in humans. Two types of intestinal infection are caused by this parasite: silent infection, which occurs in the majority of cases, and invasive disease, which affects 10% of infected persons. To understand the intestinal pathogenic process, several in vitro models, such as cell cultures, human tissue explants or human intestine xenografts in mice, have been employed. Nevertheless, our knowledge on the early steps of amebic intestinal infection and the molecules involved during human-parasite interaction is scarce, in part due to limitations in the experimental settings. In the present work, we took advantage of tissue engineering approaches to build a three-dimensional (3D)-intestinal model that is able to replicate the general characteristics of the human colon. This system consists of an epithelial layer that develops tight and adherens junctions, a mucus layer and a lamina propria-like compartment made up of collagen containing macrophages and fibroblast. By means of microscopy imaging, omics assays and the evaluation of immune responses, we show a very dynamic interaction between E. histolytica and the 3D-intestinal model. Our data highlight the importance of several virulence markers occurring in patients or in experimental models, but they also demonstrate the involvement of under described molecules and regulatory factors in the amoebic invasive process.

Macromolecular crowding links ribosomal protein gene dosage to growth rate in Vibrio cholerae.

BACKGROUND: In fast-growing bacteria, the genomic location of ribosomal protein (RP) genes is biased towards the replication origin (oriC). This trait allows optimizing their expression during exponential phase since oriC neighboring regions are in higher dose due to multifork replication. Relocation of s10-spc-α locus (S10), which codes for most of the RP, to ectopic genomic positions shows that its relative distance to the oriC correlates to a reduction on its dosage, its expression, and bacterial growth rate. However, a mechanism linking S10 dosage to cell physiology has still not been determined. RESULTS: We hypothesized that S10 dosage perturbations impact protein synthesis capacity. Strikingly, we observed that in Vibrio cholerae, protein production capacity was independent of S10 position. Deep sequencing revealed that S10 relocation altered chromosomal replication dynamics and genome-wide transcription. Such changes increased as a function of oriC-S10 distance. Since RP constitutes a large proportion of cell mass, lower S10 dosage could lead to changes in macromolecular crowding, impacting cell physiology. Accordingly, cytoplasm fluidity was higher in mutants where S10 is most distant from oriC. In hyperosmotic conditions, when crowding differences are minimized, the growth rate and replication dynamics were highly alleviated in these strains. CONCLUSIONS: The genomic location of RP genes ensures its optimal dosage. However, besides of its essential function in translation, their genomic position sustains an optimal macromolecular crowding essential for maximizing growth. Hence, this could be another mechanism coordinating DNA replication to bacterial growth.

Enzymatic synthesis of random sequences of RNA and RNA analogues by DNA polymerase theta mutants for the generation of aptamer libraries.

Nucleic acid aptamers, especially RNA, exhibit valuable advantages compared to protein therapeutics in terms of size, affinity and specificity. However, the synthesis of libraries of large random RNAs is still difficult and expensive. The engineering of polymerases able to directly generate these libraries has the potential to replace the chemical synthesis approach. Here, we start with a DNA polymerase that already displays a significant template-free nucleotidyltransferase activity, human DNA polymerase theta, and we mutate it based on the knowledge of its three-dimensional structure as well as previous mutational studies on members of the same polA family. One mutant exhibited a high tolerance towards ribonucleotides (NTPs) and displayed an efficient ribonucleotidyltransferase activity that resulted in the assembly of long RNA polymers. HPLC analysis and RNA sequencing of the products were used to quantify the incorporation of the four NTPs as a function of initial NTP concentrations and established the randomness of each generated nucleic acid sequence. The same mutant revealed a propensity to accept other modified nucleotides and to extend them in long fragments. Hence, this mutant can deliver random natural and modified RNA polymers libraries ready to use for SELEX, with custom lengths and balanced or unbalanced ratios.

The Legionella pneumophila genome evolved to accommodate multiple regulatory mechanisms controlled by the CsrA-system.

The carbon storage regulator protein CsrA regulates cellular processes post-transcriptionally by binding to target-RNAs altering translation efficiency and/or their stability. Here we identified and analyzed the direct targets of CsrA in the human pathogen Legionella pneumophila. Genome wide transcriptome, proteome and RNA co-immunoprecipitation followed by deep sequencing of a wild type and a csrA mutant strain identified 479 RNAs with potential CsrA interaction sites located in the untranslated and/or coding regions of mRNAs or of known non-coding sRNAs. Further analyses revealed that CsrA exhibits a dual regulatory role in virulence as it affects the expression of the regulators FleQ, LqsR, LetE and RpoS but it also directly regulates the timely expression of over 40 Dot/Icm substrates. CsrA controls its own expression and the stringent response through a regulatory feedback loop as evidenced by its binding to RelA-mRNA and links it to quorum sensing and motility. CsrA is a central player in the carbon, amino acid, fatty acid metabolism and energy transfer and directly affects the biosynthesis of cofactors, vitamins and secondary metabolites. We describe the first L. pneumophila riboswitch, a thiamine pyrophosphate riboswitch whose regulatory impact is fine-tuned by CsrA, and identified a unique regulatory mode of CsrA, the active stabilization of RNA anti-terminator conformations inside a coding sequence preventing Rho-dependent termination of the gap operon through transcriptional polarity effects. This allows L. pneumophila to regulate the pentose phosphate pathway and the glycolysis combined or individually although they share genes in a single operon. Thus the L. pneumophila genome has evolved to acclimate at least five different modes of regulation by CsrA giving it a truly unique position in its life cycle.

Expansion of the SOS regulon of Vibrio cholerae through extensive transcriptome analysis and experimental validation.

BACKGROUND: The SOS response is an almost ubiquitous response of cells to genotoxic stresses. The full complement of genes in the SOS regulon for Vibrio species has only been addressed through bioinformatic analyses predicting LexA binding box consensus and in vitro validation. Here, we perform whole transcriptome sequencing from Vibrio cholerae treated with mitomycin C as an SOS inducer to characterize the SOS regulon and other pathways affected by this treatment. RESULTS: Comprehensive transcriptional profiling allowed us to define the full landscape of promoters and transcripts active in V. cholerae. We performed extensive transcription start site (TSS) mapping as well as detection/quantification of the coding and non-coding RNA (ncRNA) repertoire in strain N16961. To improve TSS detection, we developed a new technique to treat RNA extracted from cells grown in various conditions. This allowed for identification of 3078 TSSs with an average 5'UTR of 116 nucleotides, and peak distribution between 16 and 64 nucleotides; as well as 629 ncRNAs. Mitomycin C treatment induced transcription of 737 genes and 28 ncRNAs at least 2 fold, while it repressed 231 genes and 17 ncRNAs. Data analysis revealed that in addition to the core genes known to integrate the SOS regulon, several metabolic pathways were induced. This study allowed for expansion of the Vibrio SOS regulon, as twelve genes (ubiEJB, tatABC, smpA, cep, VC0091, VC1190, VC1369-1370) were found to be co-induced with their adjacent canonical SOS regulon gene(s), through transcriptional read-through. Characterization of UV and mitomycin C susceptibility for mutants of these newly identified SOS regulon genes and other highly induced genes and ncRNAs confirmed their role in DNA damage rescue and protection. CONCLUSIONS: We show that genotoxic stress induces a pervasive transcriptional response, affecting almost 20% of the V. cholerae genes. We also demonstrate that the SOS regulon is larger than previously known, and its syntenic organization is conserved among Vibrio species. Furthermore, this specific co-localization is found in other γ-proteobacteria for genes recN-smpA and rmuC-tatABC, suggesting SOS regulon conservation in this phylum. Finally, we comment on the limitations of widespread NGS approaches for identification of all RNA species in bacteria.

Common Cell Shape Evolution of Two Nasopharyngeal Pathogens.

Respiratory infectious diseases are the third cause of worldwide death. The nasopharynx is the portal of entry and the ecological niche of many microorganisms, of which some are pathogenic to humans, such as Neisseria meningitidis and Moraxella catarrhalis. These microbes possess several surface structures that interact with the actors of the innate immune system. In our attempt to understand the past evolution of these bacteria and their adaption to the nasopharynx, we first studied differences in cell wall structure, one of the strongest immune-modulators. We were able to show that a modification of peptidoglycan (PG) composition (increased proportion of pentapeptides) and a cell shape change from rod to cocci had been selected for along the past evolution of N. meningitidis. Using genomic comparison across species, we correlated the emergence of the new cell shape (cocci) with the deletion, from the genome of N. meningitidis ancestor, of only one gene: yacF. Moreover, the reconstruction of this genetic deletion in a bacterium harboring the ancestral version of the locus together with the analysis of the PG structure, suggest that this gene is coordinating the transition from cell elongation to cell division. Accompanying the loss of yacF, the elongation machinery was also lost by several of the descendants leading to the change in the PG structure observed in N. meningitidis. Finally, the same evolution was observed for the ancestor of M. catarrhalis. This suggests a strong selection of these genetic events during the colonization of the nasopharynx. This selection may have been forced by the requirement of evolving permissive interaction with the immune system, the need to reduce the cellular surface exposed to immune attacks without reducing the intracellular storage capacity, or the necessity to better compete for adhesion to target cells.

GH16 and GH81 family ß-(1,3)-glucanases in Aspergillus fumigatus are essential for conidial cell wall morphogenesis.

The fungal cell wall is a rigid structure because of fibrillar and branched ß-(1,3)-glucan linked to chitin. Softening of the cell wall is an essential phenomenon during fungal morphogenesis, wherein rigid cell wall structures are cleaved by glycosylhydrolases. During the search for glycosylhydrolases acting on ß-(1,3)-glucan, we identified seven genes in the Aspergillus fumigatus genome coding for potential endo-ß-(1,3)-glucanase. ENG1 (previously characterized and named ENGL1, Mouyna et al., ), belongs to the Glycoside-Hydrolase 81 (GH81) family, while ENG2 to ENG7, to GH16 family. ENG1 and four GH16 genes (ENG2-5) were expressed in the resting conidia as well as during germination, suggesting an essential role during A. fumigatus morphogenesis. Here, we report the effect of sequential deletion of AfENG2-5 (GH16) followed by AfENG1 (GH81) deletion in the Δeng2,3,4,5 mutant. The Δeng1,2,3,4,5 mutant showed conidial defects, with linear chains of conidia unable to separate while the germination rate was not affected. These results show, for the first time in a filamentous fungus, that endo ß-(1,3)-glucanases are essential for proper conidial cell wall assembly and thus segregation of conidia during conidiation.

Repressor activity of the RpoS/σS-dependent RNA polymerase requires DNA binding.

The RpoS/σ(S) sigma subunit of RNA polymerase (RNAP) activates transcription of stationary phase genes in many Gram-negative bacteria and controls adaptive functions, including stress resistance, biofilm formation and virulence. In this study, we address an important but poorly understood aspect of σ(S)-dependent control, that of a repressor. Negative regulation by σ(S) has been proposed to result largely from competition between σ(S) and other σ factors for binding to a limited amount of core RNAP (E). To assess whether σ(S) binding to E alone results in significant downregulation of gene expression by other σ factors, we characterized an rpoS mutant of Salmonella enterica serovar Typhimurium producing a σ(S) protein proficient for Eσ(S) complex formation but deficient in promoter DNA binding. Genome expression profiling and physiological assays revealed that this mutant was defective for negative regulation, indicating that gene repression by σ(S) requires its binding to DNA. Although the mechanisms of repression by σ(S) are likely specific to individual genes and environmental conditions, the study of transcription downregulation of the succinate dehydrogenase operon suggests that σ competition at the promoter DNA level plays an important role in gene repression by Eσ(S).

RNase J depletion leads to massive changes in mRNA abundance in Helicobacter pylori.

Degradation of RNA as an intermediate message between genes and corresponding proteins is important for rapid attenuation of gene expression and maintenance of cellular homeostasis. This process is controlled by ribonucleases that have different target specificities. In the bacterial pathogen Helicobacter pylori, an exo- and endoribonuclease RNase J is essential for growth. To explore the role of RNase J in H. pylori, we identified its putative targets at a global scale using next generation RNA sequencing. We found that strong depletion for RNase J led to a massive increase in the steady-state levels of non-rRNAs. mRNAs and RNAs antisense to open reading frames were most affected with over 80% increased more than 2-fold. Non-coding RNAs expressed in the intergenic regions were much less affected by RNase J depletion. Northern blotting of selected messenger and non-coding RNAs validated these results. Globally, our data suggest that RNase J of H. pylori is a major RNase involved in degradation of most cellular RNAs.

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

Clostridium difficile is an emergent pathogen, and the most common cause of nosocomial diarrhea. In an effort to understand the role of small noncoding RNAs (sRNAs) in C. difficile physiology and pathogenesis, we used an in silico approach to identify 511 sRNA candidates in both intergenic and coding regions. In parallel, RNA-seq and differential 5'-end RNA-seq were used for global identification of C. difficile sRNAs and their transcriptional start sites at three different growth conditions (exponential growth phase, stationary phase, and starvation). This global experimental approach identified 251 putative regulatory sRNAs including 94 potential trans riboregulators located in intergenic regions, 91 cis-antisense RNAs, and 66 riboswitches. Expression of 35 sRNAs was confirmed by gene-specific experimental approaches. Some sRNAs, including an antisense RNA that may be involved in control of C. difficile autolytic activity, showed growth phase-dependent expression profiles. Expression of each of 16 predicted c-di-GMP-responsive riboswitches was observed, and experimental evidence for their regulatory role in coordinated control of motility and biofilm formation was obtained. Finally, we detected abundant sRNAs encoded by multiple C. difficile CRISPR loci. These RNAs may be important for C. difficile survival in bacteriophage-rich gut communities. Altogether, this first experimental genome-wide identification of C. difficile sRNAs provides a firm basis for future RNome characterization and identification of molecular mechanisms of sRNA-based regulation of gene expression in this emergent enteropathogen.

Single nucleotide resolution RNA-seq uncovers new regulatory mechanisms in the opportunistic pathogen Streptococcus agalactiae.

BACKGROUND: Streptococcus agalactiae, or Group B Streptococcus, is a leading cause of neonatal infections and an increasing cause of infections in adults with underlying diseases. In an effort to reconstruct the transcriptional networks involved in S. agalactiae physiology and pathogenesis, we performed an extensive and robust characterization of its transcriptome through a combination of differential RNA-sequencing in eight different growth conditions or genetic backgrounds and strand-specific RNA-sequencing. RESULTS: Our study identified 1,210 transcription start sites (TSSs) and 655 transcript ends as well as 39 riboswitches and cis-regulatory regions, 39 cis-antisense non-coding RNAs and 47 small RNAs potentially acting in trans. Among these putative regulatory RNAs, ten were differentially expressed in response to an acid stress and two riboswitches sensed directly or indirectly the pH modification. Strikingly, 15% of the TSSs identified were associated with the incorporation of pseudo-templated nucleotides, showing that reiterative transcription is a pervasive process in S. agalactiae. In particular, 40% of the TSSs upstream genes involved in nucleotide metabolism show reiterative transcription potentially regulating gene expression, as exemplified for pyrG and thyA encoding the CTP synthase and the thymidylate synthase respectively. CONCLUSIONS: This comprehensive map of the transcriptome at the single nucleotide resolution led to the discovery of new regulatory mechanisms in S. agalactiae. It also provides the basis for in depth analyses of transcriptional networks in S. agalactiae and of the regulatory role of reiterative transcription following variations of intra-cellular nucleotide pools.

Quantification of stochastic noise of splicing and polyadenylation in Entamoeba histolytica.

Alternative splicing and polyadenylation were observed pervasively in eukaryotic messenger RNAs. These alternative isoforms could either be consequences of physiological regulation or stochastic noise of RNA processing. To quantify the extent of stochastic noise in splicing and polyadenylation, we analyzed the alternative usage of splicing and polyadenylation sites in Entamoeba histolytica using RNA-Seq. First, we identified a large number of rarely spliced alternative junctions and then showed that the occurrence of these alternative splicing events is correlated with splicing site sequence, occurrence of constitutive splicing events and messenger RNA abundance. Our results implied the majority of these alternative splicing events are likely to be stochastic error of splicing machineries, and we estimated the corresponding error rates. Second, we observed extensive microheterogeneity of polyadenylation cleavage sites, and the extent of such microheterogeneity is correlated with the occurrence of constitutive cleavage events, suggesting most of such microheterogeneity is likely to be stochastic. Overall, we only observed a small fraction of alternative splicing and polyadenylation isoforms that are unlikely to be solely stochastic, implying the functional relevance of alternative splicing and polyadenylation in E. histolytica is limited. Lastly, we revised the gene models and annotated their 3'UTR in AmoebaDB, providing valuable resources to the community.