TSS-EMOTE, a refined protocol for a more complete and less biased global mapping of transcription start sites in bacterial pathogens

Staphylococcal sRNA IsrR downregulates methylthiotransferase MiaB under iron-deficient conditions

Staphylococcus aureus is a major contributor to bacterial-associated mortality, owing to its exceptional adaptability across diverse environments. Iron is vital to most organisms but can be toxic in excess. To manage its intracellular iron, S. aureus, like many pathogens, employs intricate systems. We have recently identified IsrR as a key regulatory RNA induced during iron starvation. Its role is to reduce the synthesis of non-essential iron-containing proteins under iron-depleted conditions. In this study, we unveil IsrR's regulatory action on MiaB, an enzyme responsible for methylthio group addition to specific sites on transfer RNAs (tRNAs). We use predictive tools and reporter fusion assays to demonstrate IsrR's binding to the Shine-Dalgarno sequence of miaB RNA, thereby impeding its translation. The effectiveness of IsrR hinges on the integrity of a specific C-rich region. As MiaB is non-essential and has iron-sulfur clusters, IsrR induction spares iron by downregulating miaB. This may improve S. aureus fitness and aid in navigating the host's nutritional immune defenses.IMPORTANCEIn many biotopes, including those found within an infected host, bacteria confront the challenge of iron deficiency. They employ various strategies to adapt to this scarcity of nutrients, one of which involves regulating iron-containing proteins through the action of small regulatory RNAs. Our study shows how IsrR, a small RNA from S. aureus, prevents the production of MiaB, a tRNA-modifying enzyme containing iron-sulfur clusters. With this illustration, we propose a new substrate for an iron-sparing small RNA, which, when downregulated, should reduce the need for iron and save it to essential functions.

Staphylococcal aconitase expression during iron deficiency is controlled by an sRNA-driven feedforward loop and moonlighting activity

Pathogenic bacteria employ complex systems to cope with metal ion shortage conditions and propagate in the host. IsrR is a regulatory RNA (sRNA) whose activity is decisive for optimum Staphylococcus aureus fitness upon iron starvation and for full virulence. IsrR down-regulates several genes encoding iron-containing enzymes to spare iron for essential processes. Here, we report that IsrR regulates the tricarboxylic acid (TCA) cycle by controlling aconitase (CitB), an iron-sulfur cluster-containing enzyme, and its transcriptional regulator, CcpE. This IsrR-dependent dual-regulatory mechanism provides an RNA-driven feedforward loop, underscoring the tight control required to prevent aconitase expression. Beyond its canonical enzymatic role, aconitase becomes an RNA-binding protein with regulatory activity in iron-deprived conditions, a feature that is conserved in S. aureus. Aconitase not only negatively regulates its own expression, but also impacts the enzymes involved in both its substrate supply and product utilization. This moonlighting activity concurrently upregulates pyruvate carboxylase expression, allowing it to compensate for the TCA cycle deficiency associated with iron scarcity. These results highlight the cascade of complex posttranscriptional regulations controlling S. aureus central metabolism in response to iron deficiency.

Critical factors for precise and efficient RNA cleavage by RNase Y in Staphylococcus aureus

Cellular processes require precise and specific gene regulation, in which continuous mRNA degradation is a major element. The mRNA degradation mechanisms should be able to degrade a wide range of different RNA substrates with high efficiency, but should at the same time be limited, to avoid killing the cell by elimination of all cellular RNA. RNase Y is a major endoribonuclease found in most Firmicutes, including Bacillus subtilis and Staphylococcus aureus. However, the molecular interactions that direct RNase Y to cleave the correct RNA molecules at the correct position remain unknown. In this work we have identified transcripts that are homologs in S. aureus and B. subtilis, and are RNase Y targets in both bacteria. Two such transcript pairs were used as models to show a functional overlap between the S. aureus and the B. subtilis RNase Y, which highlighted the importance of the nucleotide sequence of the RNA molecule itself in the RNase Y targeting process. Cleavage efficiency is driven by the primary nucleotide sequence immediately downstream of the cleavage site and base-pairing in a secondary structure a few nucleotides downstream. Cleavage positioning is roughly localised by the downstream secondary structure and fine-tuned by the nucleotide immediately upstream of the cleavage. The identified elements were sufficient for RNase Y-dependent cleavage, since the sequence elements from one of the model transcripts were able to convert an exogenous non-target transcript into a target for RNase Y.

The Staphylococcus aureus CamS lipoprotein is a repressor of toxin production that shapes host-pathogen interaction

Lipoproteins of the opportunistic pathogen Staphylococcus aureus play a crucial role in various cellular processes and host interactions. Consisting of a protein and a lipid moiety, they support nutrient acquisition and anchor the protein to the bacterial membrane. Recently, we identified several processed and secreted small linear peptides that derive from the secretion signal sequence of S. aureus lipoproteins. Here, we show, for the first time, that the protein moiety of the S. aureus lipoprotein CamS has a biological role that is distinct from its associated linear peptide staph-cAM373. The small peptide was shown to be involved in interspecies horizontal gene transfer, the primary mechanism for the dissemination of antibiotic resistance among bacteria. We provide evidence that the CamS protein moiety is a potent repressor of cytotoxins, such as α-toxin and leukocidins. The CamS-mediated suppression of toxin transcription was reflected by altered disease severity in in vivo infection models involving skin and soft tissue, as well as bloodstream infections. Collectively, we have uncovered the role of the protein moiety of the staphylococcal lipoprotein CamS as a previously uncharacterized repressor of S. aureus toxin production, which consequently regulates virulence and disease outcomes. Notably, the camS gene is conserved in S. aureus, and we also demonstrated the muted transcriptional response of cytotoxins in 2 different S. aureus lineages. Our findings provide the first evidence of distinct biological functions of the protein moiety and its associated linear peptide for a specific lipoprotein. Therefore, lipoproteins in S. aureus consist of 3 functional components: a lipid moiety, a protein moiety, and a small linear peptide, with putative different biological roles that might not only determine the outcome of host-pathogen interactions but also drive the acquisition of antibiotic resistance determinants.

The two-component system WalKR provides an essential link between cell wall homeostasis and DNA replication in Staphylococcus aureus

IMPORTANCE

The opportunistic human pathogen Staphylococcus aureus uses an array of protein sensing systems called two-component systems (TCS) to sense environmental signals and adapt its physiology in response by regulating different genes. This sensory network is key to S. aureus versatility and success as a pathogen. Here, we reveal for the first time the full extent of the regulatory network of WalKR, the only staphylococcal TCS that is indispensable for survival under laboratory conditions. We found that WalKR is a master regulator of cell growth, coordinating the expression of genes from multiple, fundamental S. aureus cellular processes, including those involved in maintaining cell wall metabolism, protein biosynthesis, nucleotide metabolism, and the initiation of DNA replication.

Genome-wide transcription start site mapping in the facultative intracellular pathogen Brucella melitensis by Capping-seq

Brucella melitensis (B. melitensis) is an important facultative intracellular bacterium that causes global zoonotic diseases. Continuous intracellular survival and replication are the main obstruction responsible for the accessibility of prevention and treatment of brucellosis. Bacteria respond to complex environment by regulating gene expression. Many regulatory factors function at loci where RNA polymerase initiates messenger RNA synthesis. However, limited gene annotation is a current obstacle for the research on expression regulation in bacteria. To improve annotation and explore potential functional sites, we proposed a novel genome-wide method called Capping-seq for transcription start site (TSS) mapping in B. melitensis. This technique combines capture of capped primary transcripts with Single Molecule Real-Time (SMRT) sequencing technology. We identified 2,369 TSSs at single nucleotide resolution by Capping-seq. TSSs analysis of Brucella transcripts showed a preference of purine on the TSS positions. Our results revealed that -35 and -10 elements of promoter contained consensus sequences of TTGNNN and TATNNN, respectively. The 5' ends analysis showed that 57% genes are associated with more than one TSS and 47% genes contain long leader regions, suggested potential complex regulation at the 5' ends of genes in B. melitensis. Moreover, we identified 52 leaderless genes that are mainly involved in the metabolic processes. Overall, Capping-seq technology provides a unique solution for TSS determination in prokaryotes. Our findings develop a systematic insight into the primary transcriptome characterization of B. melitensis. This study represents a critical basis for investigating gene regulation and pathogenesis of Brucella.

RNase III-CLASH of multi-drug resistant Staphylococcus aureus reveals a regulatory mRNA 3'UTR required for intermediate vancomycin resistance

Treatment of methicillin-resistant Staphylococcus aureus infections is dependent on the efficacy of last-line antibiotics including vancomycin. Treatment failure is commonly linked to isolates with intermediate vancomycin resistance (termed VISA). These isolates have accumulated point mutations that collectively reduce vancomycin sensitivity, often by thickening the cell wall. Changes in regulatory small RNA expression have been correlated with antibiotic stress in VISA isolates however the functions of most RNA regulators is unknown. Here we capture RNA-RNA interactions associated with RNase III using CLASH. RNase III-CLASH uncovers hundreds of novel RNA-RNA interactions in vivo allowing functional characterisation of many sRNAs for the first time. Surprisingly, many mRNA-mRNA interactions are recovered and we find that an mRNA encoding a long 3' untranslated region (UTR) (termed vigR 3'UTR) functions as a regulatory 'hub' within the RNA-RNA interaction network. We demonstrate that the vigR 3'UTR promotes expression of folD and the cell wall lytic transglycosylase isaA through direct mRNA-mRNA base-pairing. Deletion of the vigR 3'UTR re-sensitised VISA to glycopeptide treatment and both isaA and vigR 3'UTR deletions impact cell wall thickness. Our results demonstrate the utility of RNase III-CLASH and indicate that S. aureus uses mRNA-mRNA interactions to co-ordinate gene expression more widely than previously appreciated.

Niche-specific genome degradation and convergent evolution shaping Staphylococcus aureus adaptation during severe infections

During severe infections, Staphylococcus aureus moves from its colonising sites to blood and tissues and is exposed to new selective pressures, thus, potentially driving adaptive evolution. Previous studies have shown the key role of the agr locus in S. aureus pathoadaptation; however, a more comprehensive characterisation of genetic signatures of bacterial adaptation may enable prediction of clinical outcomes and reveal new targets for treatment and prevention of these infections. Here, we measured adaptation using within-host evolution analysis of 2590 S. aureus genomes from 396 independent episodes of infection. By capturing a comprehensive repertoire of single nucleotide and structural genome variations, we found evidence of a distinctive evolutionary pattern within the infecting populations compared to colonising bacteria. These invasive strains had up to 20-fold enrichments for genome degradation signatures and displayed significantly convergent mutations in a distinctive set of genes, linked to antibiotic response and pathogenesis. In addition to agr-mediated adaptation, we identified non-canonical, genome-wide significant loci including sucA-sucB and stp1. The prevalence of adaptive changes increased with infection extent, emphasising the clinical significance of these signatures. These findings provide a high-resolution picture of the molecular changes when S. aureus transitions from colonisation to severe infection and may inform correlation of infection outcomes with adaptation signatures.

6S RNA-Dependent Susceptibility to RNA Polymerase Inhibitors

Bacterial small RNAs (sRNAs) contribute to a variety of regulatory mechanisms that modulate a wide range of pathways, including metabolism, virulence, and antibiotic resistance. We investigated the involvement of sRNAs in rifampicin resistance in the opportunistic pathogen Staphylococcus aureus. Using a competition assay with an sRNA mutant library, we identified 6S RNA as being required for protection against low concentrations of rifampicin, an RNA polymerase (RNAP) inhibitor. This effect applied to rifabutin and fidaxomicin, two other RNAP-targeting antibiotics. 6S RNA is highly conserved in bacteria, and its absence in two other major pathogens, Salmonella enterica and Clostridioides difficile, also impaired susceptibility to RNAP inhibitors. In S. aureus, 6S RNA is produced from an autonomous gene and accumulates in stationary phase. In contrast to what was reported for Escherichia coli, S. aureus 6S RNA does not appear to play a critical role in the transition from exponential to stationary phase but affects σ^B-regulated expression in prolonged stationary phase. Nevertheless, its protective effect against rifampicin is independent of alternative sigma factor σ^B activity. Our results suggest that 6S RNA helps maintain RNAP-σ^A integrity in S. aureus, which could in turn help bacteria withstand low concentrations of RNAP inhibitors.

In silico analysis of promoter regions to identify regulatory elements in TetR family transcriptional regulatory genes of Mycobacterium colombiense CECT 3035

Background

Mycobacterium colombiense is an acid-fast, non-motile, rod-shaped mycobacterium confirmed to cause respiratory disease and disseminated infection in immune-compromised patients, and lymphadenopathy in immune-competent children. It has virulence mechanisms that allow them to adapt, survive, replicate, and produce diseases in the host. To tackle the diseases caused by M. colombiense, understanding of the regulation mechanisms of its genes is important. This paper, therefore, analyzes transcription start sites, promoter regions, motifs, transcription factors, and CpG islands in TetR family transcriptional regulatory (TFTR) genes of M. colombiense CECT 3035 using neural network promoter prediction, MEME, TOMTOM algorithms, and evolutionary analysis with the help of MEGA-X.

Results

The analysis of 22 protein coding TFTR genes of M. colombiense CECT 3035 showed that 86.36% and 13.64% of the gene sequences had one and two TSSs, respectively. Using MEME, we identified five motifs (MTF1, MTF2, MTF3, MTF4, and MTF5) and MTF1 was revealed as the common promoter motif for 100% TFTR genes of M. colombiense CECT 3035 which may serve as binding site for transcription factors that shared a minimum homology of 95.45%. MTF1 was compared to the registered prokaryotic motifs and found to match with 15 of them. MTF1 serves as the binding site mainly for AraC, LexA, and Bacterial histone-like protein families. Other protein families such as MATP, RR, σ-70 factor, TetR, LytTR, LuxR, and NAP also appear to be the binding candidates for MTF1. These families are known to have functions in virulence mechanisms, metabolism, quorum sensing, cell division, and antibiotic resistance. Furthermore, it was found that TFTR genes of M. colombiense CECT 3035 have many CpG islands with several fragments in their CpG islands. Molecular evolutionary genetic analysis showed close relationship among the genes.

Conclusion

We believe these findings will provide a better understanding of the regulation of TFTR genes in M. colombiense CECT 3035 involved in vital processes such as cell division, pathogenesis, and drug resistance and are likely to provide insights for drug development important to tackle the diseases caused by this mycobacterium. We believe this is the first report of in silico analyses of the transcriptional regulation of M. colombiense TFTR genes.

Expanding the Staphylococcus aureus SarA Regulon to Small RNAs

SarA, a transcriptional regulator of Staphylococcus aureus, is a major global regulatory system that coordinates the expression of target genes involved in its pathogenicity. Various studies have identified a large number of SarA target genes, but an in-depth characterization of the sarA regulon, including small regulatory RNAs (sRNAs), has not yet been done. In this study, we utilized transcriptome sequencing (RNA-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) to determine a comprehensive list of SarA-regulated targets, including both mRNAs and sRNAs. RNA-Seq analysis indicated 390 mRNAs and 51 sRNAs differentially expressed in a ΔsarA mutant, while ChIP-Seq revealed 354 mRNAs and 55 sRNA targets in the S. aureus genome. We confirmed the authenticity of several novel SarA targets by Northern blotting and electrophoretic mobility shift assays. Among them, we characterized repression of sprG2, a gene that encodes the toxin of a type I toxin-antitoxin system, indicating a multilayer lockdown of toxin expression by both SarA and its cognate antitoxin, SprF2. Finally, a novel SarA consensus DNA binding sequence was generated using the upstream promoter sequences of 15 novel SarA-regulated sRNA targets. A genome-wide scan with a deduced SarA motif enabled the discovery of new potential SarA target genes which were not identified in our RNA-Seq and ChIP-Seq analyses. The strength of this new consensus was confirmed with one predicted sRNA target. The RNA-Seq and ChIP-Seq combinatory analysis gives a snapshot of the regulation, whereas bioinformatic analysis reveals a permanent view of targets based on sequence. Altogether these experimental and in silico methodologies are effective to characterize transcriptional factor (TF) regulons and functions.

IMPORTANCEStaphylococcus aureus, a commensal and opportunist pathogen, is responsible for a large number of human and animal infections, from benign to severe. Gene expression adaptation during infection requires a complex network of regulators, including transcriptional factors (TF) and sRNAs. TF SarA influences virulence, metabolism, biofilm formation, and resistance to some antibiotics. SarA directly regulates expression of around 20 mRNAs and a few sRNAs. Here, we combined high-throughput expression screening methods combined with binding assays and bioinformatics for an in-depth investigation of the SarA regulon. This combinatory approach allowed the identification of 85 unprecedented mRNAs and sRNAs targets, with at least 14 being primary. Among novel SarA direct targets, we characterized repression of sprG2, a gene that encodes the toxin of a toxin-antitoxin system, indicating a multilayer lockdown of toxin expression by both SarA and its cognate antitoxin, SprF2.

RNase J1 and J2 Are Host-Encoded Factors for Plasmid Replication

Plasmids need to ensure their transmission to both daughter-cells when their host divides, but should at the same time avoid overtaxing their hosts by directing excessive host-resources toward production of plasmid factors. Naturally occurring plasmids have therefore evolved regulatory mechanisms to restrict their copy-number in response to the volume of the cytoplasm. In many plasmid families, copy-number control is mediated by a small plasmid-specified RNA, which is continuously produced and rapidly degraded, to ensure that its concentration is proportional to the current plasmid copy-number. We show here that pSA564 from the RepA_N-family is regulated by a small antisense RNA (RNA1), which, when over-expressed in trans, blocks plasmid replication and cures the bacterial host. The 5' untranslated region (5'UTR) of the plasmid replication initiation gene (repA) potentially forms two mutually exclusive secondary structures, ON and OFF, where the latter both sequesters the repA ribosome binding site and acts as a rho-independent transcriptional terminator. Duplex formation between RNA1 and the 5'UTR shifts the equilibrium to favor the putative OFF-structure, enabling a single small RNA to down-regulate repA expression at both transcriptional and translational levels. We further examine which sequence elements on the antisense RNA and on its 5'UTR target are needed for this regulation. Finally, we identify the host-encoded exoribonucleases RNase J1 and J2 as the enzymes responsible for rapidly degrading the replication-inhibiting section of RNA1. This region accumulates and blocks RepA expression in the absence of either RNase J1 or J2, which are therefore essential host factors for pSA564 replication in Staphylococcus aureus.

Essential gene analysis in Acinetobacter baumannii by high-density transposon mutagenesis and CRISPR interference

Acinetobacter baumannii is a poorly understood bacterium capable of life-threatening infections in hospitals. Few antibiotics remain effective against this highly resistant pathogen. Developing rationally-designed antimicrobials that can target A. baumannii requires improved knowledge of the proteins that carry out essential processes allowing growth of the organism. Unfortunately, studying essential genes has been challenging using traditional techniques, which usually require time-consuming recombination-based genetic manipulations. Here, we performed saturating mutagenesis with dual transposon systems to identify essential genes in A. baumannii and we developed a CRISPR-interference (CRISPRi) system for facile analysis of these genes. We show that the CRISPRi system enables efficient transcriptional silencing in A. baumannii Using these tools, we confirmed the essentiality of the novel cell division protein AdvA and discovered a previously uncharacterized AraC-family transcription factor (ACX60_RS03245) that is necessary for growth. In addition, we show that capsule biosynthesis is a conditionally essential process, with mutations in late-acting steps causing toxicity in strain ATCC 17978 that can be bypassed by blocking early-acting steps or activating the BfmRS stress response. These results open new avenues for analysis of essential pathways in A. baumannii ImportanceNew approaches are urgently needed to control A. baumannii, one of the most drug resistant pathogens known. To facilitate the development of novel targets that allow inhibition of the pathogen, we performed a large-scale identification of genes whose products the bacterium needs for growth. We also developed a CRISPR-based gene knockdown tool that operates efficiently in A. baumannii, allowing rapid analysis of these essential genes. We used these methods to define multiple processes vital to the bacterium, including a previously uncharacterized gene-regulatory factor and export of a protective polymeric capsule. These tools will enhance our ability to investigate processes critical for the essential biology of this challenging hospital-acquired pathogen.

A CRISPR interference platform for selective downregulation of gene expression in Borrelia burgdorferi

The spirochete Borrelia burgdorferi causes Lyme disease, an increasingly prevalent infection. While previous studies have provided important insight into B. burgdorferi biology, many aspects, including basic cellular processes, remain underexplored. To help speed up the discovery process, we adapted a CRISPR interference (CRISPRi) platform for use in B. burgdorferi For efficiency and flexibility of use, we generated various CRISPRi template constructs that produce different basal and induced levels of dcas9 and carry different antibiotic resistance markers. We characterized the effectiveness of our CRISPRi platform by targeting the motility and cell morphogenesis genes flaB, mreB, rodA, and ftsI, whose native expression levels span two orders of magnitude. For all four genes, we obtained gene repression efficiencies of at least 95%. We showed by darkfield microscopy and cryo-electron tomography that flagellin (FlaB) depletion reduced the length and number of periplasmic flagella, which impaired cellular motility and resulted in cell straightening. Depletion of FtsI caused cell filamentation, implicating this protein in cell division in B. burgdorferi Finally, localized cell bulging in MreB- and RodA-depleted cells matched the locations of new peptidoglycan insertion specific to spirochetes of the Borrelia genus. These results therefore implicate MreB and RodA in the particular mode of cell wall elongation of these bacteria. Collectively, our results demonstrate the efficiency and ease of use of our B. burgdorferi CRISPRi platform, which should facilitate future genetic studies of this important pathogen.IMPORTANCE Gene function studies are facilitated by the availability of rapid and easy-to-use genetic tools. Homologous recombination-based methods traditionally used to genetically investigate gene function remain cumbersome to perform in B. burgdorferi, as they often are relatively inefficient. In comparison, our CRISPRi platform offers an easy and fast method to implement as it only requires a single plasmid transformation step and IPTG addition to obtain potent (>95%) downregulation of gene expression. To facilitate studies of various genes in wild-type and genetically modified strains, we provide over 30 CRISPRi plasmids that produce distinct levels of dcas9 expression and carry different antibiotic resistance markers. Our CRISPRi platform represents a useful and efficient complement to traditional genetic and chemical methods to study gene function in B. burgdorferi.

Construction of prokaryotic strand-specific primary-transcripts saturated RNASeq library by controlled heat magnesium-dependent mRNA degradation

The main limiting factors for RNA-Seq analysis are quality and quantity of the isolated mRNA. In prokaryotes, the proportion of messenger RNA to total RNA is rather low. Therefore, the main strategy of library preparation for sequencing is mRNA enrichment. Ribosomal and transfer RNAs, both monophosphorylated at the 5'-ends, are the major fractions of total RNA, while the bulk of primary transcripts is triphosphorylated at the 5'-teminus. Due to its low molecular weight, transfer RNA could be easily removed by a quick precipitation in LiCl solution. Ribosomal RNA may be degraded enzymatically by 5'-end terminal exonuclease XRN-1. These steps allow enriching samples in mRNA during the first stages of RNA-Seq library preparation. The desired level of fragmentation of enriched mRNA necessary for the 2^nd generation sequencing can be controlled by the duration of incubation at elevated temperatures in the presence of Mg²⁺-ions. Here, we describe a simple protocol for construction of the primary prokaryotic mRNA-saturated library without long depletion procedures.

Dynamic insights on transcription initiation and RNA processing during bacterial adaptation

Transcription initiation and RNA processing govern gene expression and enable bacterial adaptation by reshaping the RNA landscape. The aim of this study was to simultaneously observe these two fundamental processes in a transcriptome responding to an environmental signal. A controlled σ^E system in E. coli was coupled to our previously described tagRNA-seq method to yield process kinetics information. Changes in transcription initiation frequencies (TIF) and RNA processing frequencies (PF) were followed using 5' RNA tags. Changes in TIF showed a binary increased/decreased pattern that alternated between transcriptionally activated and repressed promoters, providing the bacterial population with transcriptional oscillation. PF variation fell into three categories of cleavage activity: (i) constant and independent of RNA levels, (ii) increased once RNA has accumulated, and (iii) positively correlated to changes in TIF. This work provides a comprehensive and dynamic view of major events leading to transcriptomic reshaping during bacterial adaptation. It unveils an interplay between transcription initiation and the activity of specific RNA cleavage sites. This study utilized a well-known genetic system to analyze fundamental processes and can serve as a blueprint for comprehensive studies that exploit the RNA metabolism to decipher and understand bacterial gene expression control.

The Staphylococcus aureus ArlRS two-component system regulates virulence factor expression through MgrA

The Gram-positive bacterium, Staphylococcus aureus, is a versatile pathogen that can sense and adapt to a wide variety of environments within the human host, in part through its 16 two-component regulatory systems. The ArlRS two-component system has been shown to affect many cellular processes in S. aureus, including autolysis, biofilm formation, capsule synthesis and virulence. Yet the molecular details of this regulation remained largely unknown. We used RNA sequencing to identify the ArlRS regulon, and found 70% overlap with that of the global regulator MgrA. These genes included cell wall-anchored adhesins (ebh, sdrD), polysaccharide and capsule synthesis genes, cell wall remodeling genes (lytN, ddh), the urease operon, genes involved in metal transport (feoA, mntH, sirA), anaerobic metabolism genes (adhE, pflA, nrdDG) and a large number of virulence factors (lukSF, lukAB, nuc, gehB, norB, chs, scn and esxA). We show that ArlR directly activates expression of mgrA and identify a probable ArlR-binding site (TTTTCTCAT-N4 -TTTTAATAA). A highly similar sequence is also found in the spx P2 promoter, which was recently shown to be regulated by ArlRS. We also demonstrate that ArlS has kinase activity toward ArlR in vitro, although it has slower kinetics than other similar histidine kinases.

APERO: a genome-wide approach for identifying bacterial small RNAs from RNA-Seq data

Small non-coding RNAs (sRNAs) regulate numerous cellular processes in all domains of life. Several approaches have been developed to identify them from RNA-seq data, which are efficient for eukaryotic sRNAs but remain inaccurate for the longer and highly structured bacterial sRNAs. We present APERO, a new algorithm to detect small transcripts from paired-end bacterial RNA-seq data. In contrast to previous approaches that start from the read coverage distribution, APERO analyzes boundaries of individual sequenced fragments to infer the 5′ and 3′ ends of all transcripts. Since sRNAs are about the same size as individual fragments (50–350 nucleotides), this algorithm provides a significantly higher accuracy and robustness, e.g., with respect to spontaneous internal breaking sites. To demonstrate this improvement, we develop a comparative assessment on datasets from Escherichia coli and Salmonella enterica, based on experimentally validated sRNAs. We also identify the small transcript repertoire of Dickeya dadantii including putative intergenic RNAs, 5′ UTR or 3′ UTR-derived RNA products and antisense RNAs. Comparisons to annotations as well as RACE-PCR experimental data confirm the precision of the detected transcripts. Altogether, APERO outperforms all existing methods in terms of sRNA detection and boundary precision, which is crucial for comprehensive genome annotations. It is freely available as an open source R package on https://github.com/Simon-Leonard/APERO

Revisiting the regulation of the capsular polysaccharide biosynthesis gene cluster in Staphylococcus aureus

Capsular polysaccharide (CP) biosynthesis in Staphylococcus aureus is tightly controlled resulting in a heterogeneous phenotype within a population and CP being mainly detectable in nongrowing cells. Expression of the corresponding biosynthesis gene cluster is driven by one promoter element (P_cap ). Here, we demonstrate that P_cap contains a main SigB-dependent promoter. The SigB consensus motif overlaps with a previously described inverted repeat (IR) that is crucial for cap expression. The essentiality of the IR is derived from this region acting as a SigB binding site rather than as an operator site for the proposed cap activators RbsR and MsaB. Furthermore, P_cap contains an extensive upstream region harboring a weak SigA-dependent promoter and binding sites for cap repressors such as SaeR, CodY and Rot. Heterogeneous CP synthesis is determined by SigB activity and repressor binding to the upstream region. SigB dependency and regulation by the upstream repressors are also sufficient to explain the temporal gene expression pattern at the transcriptional level. However, CP synthesis remains growth phase-dependent even when transcription is rendered constitutive, suggesting additional posttranscriptional regulatory circuits. Thus, the interference of multiple repressors with SigB-dependent promoter activity as well as post-transcriptional mechanisms ensure the appropriate regulation of CP synthesis.

Architecture of divergent flagellar promoters controlled by CtrA in Rhodobacter sphaeroides

BACKGROUND

Rhodobacter sphaeroides has two sets of flagellar genes, fla1 and fla2, that are responsible for the synthesis of two different flagellar structures. The expression of the fla2 genes is under control of CtrA. In several α-proteobacteria CtrA is also required for the expression of the flagellar genes, but the architecture of CtrA-dependent promoters has only been studied in detail in Caulobacter crescentus. In many cases the expression of fla genes originates from divergent promoters located a few base pairs apart, suggesting a particular arrangement of the cis-acting sites.

RESULTS

Here we characterized several control regions of the R. sphaeroides fla2 genes and analyzed in detail two regions containing the divergent promoters flgB2p-fliI2p, and fliL2p-fliF2p. Binding sites for CtrA of these promoters were identified in silico and tested by site directed mutagenesis. We conclude that each one of these promoter regions has a particular arrangement, either a single CtrA binding site for activation of fliL2p and fliF2p, or two independent sites for activation of flgB2p and fliI2p. ChIP experiments confirmed that CtrA binds to the control region containing the flgB2 and fliI2 promoters, supporting the notion that CtrA directly controls the expression of the fla2 genes. The flgB and fliI genes are syntenic and show a short intercistronic region in closely related bacterial species. We analyzed these regions and found that the arrangement of the CtrA binding sites varies considerably.

CONCLUSIONS

The results in this work reveal the arrangement of the fla2 divergent promoters showing that CtrA promotes transcriptional activation using more than a single architecture.

Mapping 5'-Ends and Their Phosphorylation State With EMOTE, TSS-EMOTE, and nEMOTE

Knowing the exact 5'-end of a specific RNA molecule has long been of great interest, since this information can reveal transcription start sites, which subsequently pinpoints promoter regions and potential binding sites for transcriptional regulators. Mapping of the 5'-end can also identify endoribonucleolytic cleavage sites, for the maturation of ribosomal RNA or events which initiate degradation of all or parts of the RNA molecule. If the 5' nucleotide of an RNA is triphosphorylated, then the nucleotide represents the original, unaltered 5'-end of the RNA. Similarly, if an endoribonuclease cleaves the RNA with a mechanism that leaves a monophosphate on the 5'-end downstream cleavage product, then it is useful to search exclusively for monophosphorylated 5'-ends to identify the cleavage sites of the endoribonuclease. This chapter details three variations of the EMOTE (exact mapping of transcriptome ends) protocol, which can identify RNA 5'-ends on a transcriptome-wide scale, and simultaneously provide information about the phosphorylation state of the detected ends. EMOTE, TSS-EMOTE (transcription start specific EMOTE), and nEMOTE (nonphosphorylated EMOTE) reveals mono-, tri-, and nonphosphorylated 5'-ends, respectively, and the protocols can be performed individually or in parallel on the same RNA sample.

Post-transcriptional control of virulence gene expression in Staphylococcus aureus

Opportunistic pathogens have to be ready to change life-style whenever the occasion arises, and therefore need to keep tight control over the expression of their virulence factors. Doubly so for commensal bacteria, such as Staphylococcus aureus, which should avoid harming their hosts when they are in a state of peaceful co-existence. S. aureus carries very few sigma factors to help define the transcriptional programs, but instead uses a plethora of small RNA molecules and RNA-RNA interactions to regulate gene expression post-transcriptionally. The endoribonucleases RNase III and RNase Y contribute to this regulatory diversity, and provide a link to RNA-decay and intra-cellular spatiotemporal control of expression. In this review we describe some of these post-transcriptional mechanisms as well as some of the novel transcriptomic approaches that have been used to find and to study them.

Genome-scale analysis of Methicillin-resistant Staphylococcus aureus USA300 reveals a tradeoff between pathogenesis and drug resistance

Staphylococcus aureus infection is a rising public health care threat. S. aureus is believed to have elaborate regulatory networks that orchestrate its virulence. Despite its importance, the systematic understanding of the transcriptional landscape of S. aureus is limited. Here, we describe the primary transcriptome landscape of an epidemic USA300 isolate of community-acquired methicillin-resistant S. aureus. We experimentally determined 1,861 transcription start sites with their principal promoter elements, including well-conserved -35 and -10 elements and weakly conserved -16 element and 5' untranslated regions containing AG-rich Shine-Dalgarno sequence. In addition, we identified 225 genes whose transcription was initiated from multiple transcription start sites, suggesting potential regulatory functions at transcription level. Along with the transcription unit architecture derived by integrating the primary transcriptome analysis with operon prediction, the measurement of differential gene expression revealed the regulatory framework of the virulence regulator Agr, the SarA-family transcriptional regulators, and β-lactam resistance regulators. Interestingly, we observed a complex interplay between virulence regulation, β-lactam resistance, and metabolism, suggesting a possible tradeoff between pathogenesis and drug resistance in the USA300 strain. Our results provide platform resource for the location of transcription initiation and an in-depth understanding of transcriptional regulation of pathogenesis, virulence, and antibiotic resistance in S. aureus.

Deep sequencing approaches for the analysis of prokaryotic transcriptional boundaries and dynamics

The identification of the protein-coding regions of a genome is straightforward due to the universality of start and stop codons. However, the boundaries of the transcribed regions, conditional operon structures, non-coding RNAs and the dynamics of transcription, such as pausing of elongation, are non-trivial to identify, even in the comparatively simple genomes of prokaryotes. Traditional methods for the study of these areas, such as tiling arrays, are noisy, labour-intensive and lack the resolution required for densely-packed bacterial genomes. Recently, deep sequencing has become increasingly popular for the study of the transcriptome due to its lower costs, higher accuracy and single nucleotide resolution. These methods have revolutionised our understanding of prokaryotic transcriptional dynamics. Here, we review the deep sequencing and data analysis techniques that are available for the study of transcription in prokaryotes, and discuss the bioinformatic considerations of these analyses.

If it transcribes, we can sequence it: mining the complexities of host-pathogen-environment interactions using RNA-seq

Host-pathogen interactions are exceedingly complex because they involve multiple host tissues, often occur in the context of normal microflora, and can span diverse microenvironments. Although decades of gene expression studies have provided detailed insights into infection processes, technical challenges have restricted experiments to single pathogenic species or host tissues. RNA-sequencing (RNA-seq) has revolutionized the study of gene expression because in addition to quantifying transcriptional output, it allows detection and characterization of all transcripts in a genome. Here, we review how refined approaches to RNA-seq are used to map the transcriptional networks that control host-pathogen interactions. These enhanced techniques include dRNA-seq and term-seq for the fine-scale mapping of transcriptional start and termination sites, and dual RNA-seq for simultaneous sequencing of host and bacterial pathogen transcriptomes. Dual RNA-seq experiments are currently limited to in vitro infection systems that do not fully reflect the complexities of the in vivo environment, thus a challenge is to develop in vivo model systems and experimental approaches that address the biological heterogeneity of host environments, followed by the integration of RNA-seq with other genome-scale datasets to identify the transcriptional networks that mediate host-pathogen interactions.

Cell Cycle Constraints and Environmental Control of Local DNA Hypomethylation in α-Proteobacteria

Heritable DNA methylation imprints are ubiquitous and underlie genetic variability from bacteria to humans. In microbial genomes, DNA methylation has been implicated in gene transcription, DNA replication and repair, nucleoid segregation, transposition and virulence of pathogenic strains. Despite the importance of local (hypo)methylation at specific loci, how and when these patterns are established during the cell cycle remains poorly characterized. Taking advantage of the small genomes and the synchronizability of α-proteobacteria, we discovered that conserved determinants of the cell cycle transcriptional circuitry establish specific hypomethylation patterns in the cell cycle model system Caulobacter crescentus. We used genome-wide methyl-N6-adenine (m6A-) analyses by restriction-enzyme-cleavage sequencing (REC-Seq) and single-molecule real-time (SMRT) sequencing to show that MucR, a transcriptional regulator that represses virulence and cell cycle genes in S-phase but no longer in G1-phase, occludes 5'-GANTC-3' sequence motifs that are methylated by the DNA adenine methyltransferase CcrM. Constitutive expression of CcrM or heterologous methylases in at least two different α-proteobacteria homogenizes m6A patterns even when MucR is present and affects promoter activity. Environmental stress (phosphate limitation) can override and reconfigure local hypomethylation patterns imposed by the cell cycle circuitry that dictate when and where local hypomethylation is instated.