Deep sequencing uncovers numerous small RNAs on all four replicons of the plant pathogen Agrobacterium tumefaciens

Methods for Bioinformatic Prediction of Genuine sRNAs from Outer Membrane Vesicles

The molecular pathogenesis of Gram-negative bacteria remains a complex and incompletely understood phenomenon. Various factors are believed to contribute to the pathogenicity of these bacteria. One key mechanism utilized by Gram-negative bacteria is the production of outer membrane vesicles (OMVs), which are small spherical particles derived from the bacterial outer membrane. These OMVs are crucial in delivering virulence factors to the host, facilitating host-pathogen interactions. Within these OMVs, small regulatory RNAs (sRNAs) have been identified as important players in modulating the host immune response. One of the main challenges in studying OMVs and their cargo of sRNAs is the difficulty in isolating and purifying sufficient quantities of OMVs, as well as accurately predicting genuine sRNAs computationally. In this chapter, we present protocols aimed at overcoming these obstacles.

The LysR-type transcription factor LsrB regulates beta-lactam resistance in Agrobacterium tumefaciens

Agrobacterium tumefaciens is a plant pathogen, broadly known as the causal agent of the crown gall disease. The soil bacterium is naturally resistant to beta-lactam antibiotics by utilizing the inducible beta-lactamase AmpC. Our picture on the condition-dependent regulation of ampC expression is incomplete. A known regulator is AmpR controlling the transcription of ampC in response to unrecycled muropeptides as a signal for cell wall stress. In our study, we uncovered the global transcriptional regulator LsrB as a critical player acting upstream of AmpR. Deletion of lsrB led to severe ampicillin and penicillin sensitivity, which could be restored to wild-type levels by lsrB complementation. By transcriptome profiling via RNA-Seq and qRT-PCR and by electrophoretic mobility shift assays, we show that ampD coding for an anhydroamidase involved in peptidoglycan recycling is under direct negative control by LsrB. Controlling AmpD levels by the LysR-type regulator in turn impacts the cytoplasmic concentration of cell wall degradation products and thereby the AmpR-mediated regulation of ampC. Our results substantially expand the existing model of inducible beta-lactam resistance in A. tumefaciens by establishing LsrB as higher-level transcriptional regulator.

Virulence and Ecology of Agrobacteria in the Context of Evolutionary Genomics

Among plant-associated bacteria, agrobacteria occupy a special place. These bacteria are feared in the field as agricultural pathogens. They cause abnormal growth deformations and significant economic damage to a broad range of plant species. However, these bacteria are revered in the laboratory as models and tools. They are studied to discover and understand basic biological phenomena and used in fundamental plant research and biotechnology. Agrobacterial pathogenicity and capability for transformation are one and the same and rely on functions encoded largely on their oncogenic plasmids. Here, we synthesize a substantial body of elegant work that elucidated agrobacterial virulence mechanisms and described their ecology. We review findings in the context of the natural diversity that has been recently unveiled for agrobacteria and emphasize their genomics and plasmids. We also identify areas of research that can capitalize on recent findings to further transform our understanding of agrobacterial virulence and ecology.

The Role of Iron in Phytopathogenic Microbe-Plant Interactions: Insights into Virulence and Host Immune Response

Iron is an essential element required for the growth and survival of nearly all forms of life. It serves as a catalytic component in multiple enzymatic reactions, such as photosynthesis, respiration, and DNA replication. However, the excessive accumulation of iron can result in cellular toxicity due to the production of reactive oxygen species (ROS) through the Fenton reaction. Therefore, to maintain iron homeostasis, organisms have developed a complex regulatory network at the molecular level. Besides catalyzing cellular redox reactions, iron also regulates virulence-associated functions in several microbial pathogens. Hosts and pathogens have evolved sophisticated strategies to compete against each other over iron resources. Although the role of iron in microbial pathogenesis in animals has been extensively studied, mechanistic insights into phytopathogenic microbe-plant associations remain poorly understood. Recent intensive research has provided intriguing insights into the role of iron in several plant-pathogen interactions. This review aims to describe the recent advances in understanding the role of iron in the lifestyle and virulence of phytopathogenic microbes, focusing on bacteria and host immune responses.

Transcriptome architecture of the three main lineages of agrobacteria

Agrobacteria are a diverse, polyphyletic group of prokaryotes with multipartite genomes capable of transferring DNA into the genomes of host plants, making them an essential tool in plant biotechnology. Despite their utility in plant transformation, genome-wide transcriptional regulation is not well understood across the three main lineages of agrobacteria. Transcription start sites (TSSs) are a necessary component of gene expression and regulation. In this study, we used differential RNA-seq and a TSS identification algorithm optimized on manually annotated TSS, then validated with existing TSS to identify thousands of TSS with nucleotide resolution for representatives of each lineage. We extend upon the 356 TSSs previously reported in Agrobacterium fabrum C58 by identifying 1,916 TSSs. In addition, we completed genomes and phenotyping of Rhizobium rhizogenes C16/80 and Allorhizobium vitis T60/94, identifying 2,650 and 2,432 TSSs, respectively. Parameter optimization was crucial for an accurate, high-resolution view of genome and transcriptional dynamics, highlighting the importance of algorithm optimization in genome-wide TSS identification and genomics at large. The optimized algorithm reduced the number of TSSs identified internal and antisense to the coding sequence on average by 90.5% and 91.9%, respectively. Comparison of TSS conservation between orthologs of the three lineages revealed differences in cell cycle regulation of ctrA as well as divergence of transcriptional regulation of chemotaxis-related genes when grown in conditions that simulate the plant environment. These results provide a framework to elucidate the mechanistic basis and evolution of pathology across the three main lineages of agrobacteria. IMPORTANCE Transcription start sites (TSSs) are fundamental for understanding gene expression and regulation. Agrobacteria, a group of prokaryotes with the ability to transfer DNA into the genomes of host plants, are widely used in plant biotechnology. However, the genome-wide transcriptional regulation of agrobacteria is not well understood, especially in less-studied lineages. Differential RNA-seq and an optimized algorithm enabled identification of thousands of TSSs with nucleotide resolution for representatives of each lineage. The results of this study provide a framework for elucidating the mechanistic basis and evolution of pathology across the three main lineages of agrobacteria. The optimized algorithm also highlights the importance of parameter optimization in genome-wide TSS identification and genomics at large.

Do small RNAs unlock the below ground microbiome-plant interaction mystery?

Over the past few decades, regulatory RNAs, such as small RNAs (sRNAs), have received increasing attention in the context of host-microbe interactions due to their diverse roles in controlling various biological processes in eukaryotes. In addition, studies have identified an increasing number of sRNAs with novel functions across a wide range of bacteria. What is not well understood is why cells regulate gene expression through post-transcriptional mechanisms rather than at the initiation of transcription. The finding of a multitude of sRNAs and their identified associated targets has allowed further investigation into the role of sRNAs in mediating gene regulation. These foundational data allow for further development of hypotheses concerning how a precise control of gene activity is accomplished through the combination of transcriptional and post-transcriptional regulation. Recently, sRNAs have been reported to participate in interkingdom communication and signalling where sRNAs originating from one kingdom are able to target or control gene expression in another kingdom. For example, small RNAs of fungal pathogens that silence plant genes and vice-versa plant sRNAs that mediate bacterial gene expression. However, there is currently a lack of evidence regarding sRNA-based inter-kingdom signalling across more than two interacting organisms. A habitat that provides an excellent opportunity to investigate interconnectivity is the plant rhizosphere, a multifaceted ecosystem where plants and associated soil microbes are known to interact. In this paper, we discuss how the interconnectivity of bacteria, fungi, and plants within the rhizosphere may be mediated by bacterial sRNAs with a particular focus on disease suppressive and non-suppressive soils. We discuss the potential roles sRNAs may play in the below-ground world and identify potential areas of future research, particularly in reference to the regulation of plant immunity genes by bacterial and fungal communities in disease-suppressive and non-disease-suppressive soils.

Motility control through an anti-activation mechanism in Agrobacterium tumefaciens

Many bacteria can migrate from a free-living, planktonic state to an attached, biofilm existence. One factor regulating this transition in the facultative plant pathogen Agrobacterium tumefaciens is the ExoR-ChvG-ChvI system. Periplasmic ExoR regulates the activity of the ChvG-ChvI two-component system in response to environmental stress, most notably low pH. ChvI impacts hundreds of genes, including those required for type VI secretion, virulence, biofilm formation, and flagellar motility. Previous studies revealed that activated ChvG-ChvI represses expression of most of class II and class III flagellar biogenesis genes, but not the master motility regulator genes visN, visR, and rem. In this study, we characterized the integration of the ExoR-ChvG-ChvI and VisNR-Rem pathways. We isolated motile suppressors of the non-motile ΔexoR mutant and thereby identified the previously unannotated mirA gene encoding a 76 amino acid protein. We report that the MirA protein interacts directly with the Rem DNA-binding domain, sequestering Rem and preventing motility gene activation. The ChvG-ChvI pathway activates mirA expression and elevated mirA is sufficient to block motility. This study reveals how the ExoR-ChvG-ChvI pathway prevents flagellar motility in A. tumefaciens. MirA is also conserved among other members of the Rhizobiales suggesting similar mechanisms of motility regulation.

Agrobacterium fabrum C58 involved nitrate reductase NapA and antisense RNA NorR to denitrify

Agrobacterium fabrum C58 is a plant-associated bacterium that is able to denitrify under anoxic conditions. The cluster of denitrification genes harbored by this strain has been well characterized. It includes nir and nor operons encoding nitrite and nitric oxide reductases, respectively. However, the reductase involved in nitrate reduction has not yet been studied and little information is available on denitrification regulators in A. fabrum C58. In this study, we aimed to (i) characterize the nitrate reductase, (ii) determine its role in A. fabrum C58 fitness and root colonization and (ii) reveal the contribution of small RNA on denitrification regulation. By constructing a mutant strain defective for napA, we demonstrated that the reduction of nitrate to nitrite was catalyzed by the periplasmic nitrate reductase, NapA. We evidenced a positive role of NapA in A. fabrum C58 fitness and suggested that A. fabrum C58 is able to use components exuded by plant roots to respire anaerobically. Here, we showed that NorR small RNA increased the level of norCBQ mRNA and a decrease of NorR is correlated with a decrease in N2O emission. Together, our results underscore the importance of understanding the denitrification pathway at the strain level in order to develop strategies to mitigate N2O production at the microbial community level.

A LysR-type transcriptional regulator controls the expression of numerous small RNAs in Agrobacterium tumefaciens

Small RNAs (sRNAs) are universal posttranscriptional regulators of gene expression and hundreds of sRNAs are frequently found in each and every bacterium. In order to coordinate cellular processes in response to ambient conditions, many sRNAs are differentially expressed. Here, we asked how these small regulators are regulated using Agrobacterium tumefaciens as a model system. Among the best-studied sRNAs in this plant pathogen are AbcR1 regulating numerous ABC transporters and PmaR, a regulator of peptidoglycan biosynthesis, motility, and ampicillin resistance. We report that the LysR-type regulator VtlR (also known as LsrB) controls expression of AbcR1 and PmaR. A vtlR/lsrB deletion strain showed growth defects, was sensitive to antibiotics and severely compromised in plant tumor formation. Transcriptome profiling by RNA-sequencing revealed more than 1,200 genes with altered expression in the mutant. Consistent with the function of VtlR/LsrB as regulator of AbcR1, many ABC transporter genes were affected. Interestingly, the transcription factor did not only control the expression of AbcR1 and PmaR. In the mutant, 102 sRNA genes were significantly up- or downregulated. Thus, our study uncovered VtlR/LsrB as the master regulator of numerous sRNAs. Thereby, the transcriptional regulator harnesses the regulatory power of sRNAs to orchestrate the expression of distinct sub-regulons.

Arginine-Rich Small Proteins with a Domain of Unknown Function, DUF1127, Play a Role in Phosphate and Carbon Metabolism of Agrobacterium tumefaciens

In any given organism, approximately one-third of all proteins have a yet-unknown function. A widely distributed domain of unknown function is DUF1127. Approximately 17,000 proteins with such an arginine-rich domain are found in 4,000 bacteria. Most of them are single-domain proteins, and a large fraction qualifies as small proteins with fewer than 50 amino acids. We systematically identified and characterized the seven DUF1127 members of the plant pathogen Agrobacterium tumefaciens They all give rise to authentic proteins and are differentially expressed as shown at the RNA and protein levels. The seven proteins fall into two subclasses on the basis of their length, sequence, and reciprocal regulation by the LysR-type transcription factor LsrB. The absence of all three short DUF1127 proteins caused a striking phenotype in later growth phases and increased cell aggregation and biofilm formation. Protein profiling and transcriptome sequencing (RNA-seq) analysis of the wild type and triple mutant revealed a large number of differentially regulated genes in late exponential and stationary growth. The most affected genes are involved in phosphate uptake, glycine/serine homeostasis, and nitrate respiration. The results suggest a redundant function of the small DUF1127 paralogs in nutrient acquisition and central carbon metabolism of A. tumefaciens They may be required for diauxic switching between carbon sources when sugar from the medium is depleted. We end by discussing how DUF1127 might confer such a global impact on cell physiology and gene expression.IMPORTANCE Despite being prevalent in numerous ecologically and clinically relevant bacterial species, the biological role of proteins with a domain of unknown function, DUF1127, is unclear. Experimental models are needed to approach their elusive function. We used the phytopathogen Agrobacterium tumefaciens, a natural genetic engineer that causes crown gall disease, and focused on its three small DUF1127 proteins. They have redundant and pervasive roles in nutrient acquisition, cellular metabolism, and biofilm formation. The study shows that small proteins have important previously missed biological functions. How small basic proteins can have such a broad impact is a fascinating prospect of future research.

Interplay between Regulatory RNAs and Signal Transduction Systems during Bacterial Infection

The ability of pathogenic bacteria to stably infect the host depends on their capacity to respond and adapt to the host environment and on the efficiency of their defensive mechanisms. Bacterial envelope provides a physical barrier protecting against environmental threats. It also constitutes an important sensory interface where numerous sensing systems are located. Signal transduction systems include Two-Component Systems (TCSs) and alternative sigma factors. These systems are able to sense and respond to the ever-changing environment inside the host, altering the bacterial transcriptome to mitigate the impact of the stress. The regulatory networks associated with signal transduction systems comprise small regulatory RNAs (sRNAs) that can be directly involved in the expression of virulence factors. The aim of this review is to describe the importance of TCS- and alternative sigma factor-associated sRNAs in human pathogens during infection. The currently available genome-wide approaches for studies of TCS-regulated sRNAs will be discussed. The differences in the signal transduction mediated by TCSs between bacteria and higher eukaryotes and the specificity of regulatory RNAs for their targets make them appealing targets for discovery of new strategies to fight against multi-resistant bacteria.

Whole-genome Identification of Transcriptional Start Sites by Differential RNA-seq in Bacteria

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

A central role for the transcriptional regulator VtlR in small RNA-mediated gene regulation in Agrobacterium tumefaciens

LysR-type transcriptional regulators (LTTRs) are the most common type of transcriptional regulators in prokaryotes and function by altering gene expression in response to environmental stimuli. In the class Alphaproteobacteria, a conserved LTTR named VtlR is critical to the establishment of host-microbe interactions. In the mammalian pathogen Brucella abortus, VtlR is required for full virulence in a mouse model of infection, and VtlR activates the expression of abcR2, which encodes a small regulatory RNA (sRNA). In the plant symbiont Sinorhizobium meliloti, the ortholog of VtlR, named LsrB, is involved in the symbiosis of the bacterium with alfalfa. Agrobacterium tumefaciens is a close relative of both B. abortus and S. meliloti, and this bacterium is the causative agent of crown gall disease in plants. In the present study, we demonstrate that VtlR is involved in the ability of A. tumefaciens to grow appropriately in artificial medium, and an A. tumefaciens vtlR deletion strain is defective in motility, biofilm formation, and tumorigenesis of potato discs. RNA-sequencing analyses revealed that more than 250 genes are dysregulated in the ∆vtlR strain, and importantly, VtlR directly controls the expression of three sRNAs in A. tumefaciens. Taken together, these data support a model in which VtlR indirectly regulates hundreds of genes via manipulation of sRNA pathways in A. tumefaciens, and moreover, while the VtlR/LsrB protein is present and structurally conserved in many members of the Alphaproteobacteria, the VtlR/LsrB regulatory circuitry has diverged in order to accommodate the unique environmental niche of each organism.

Structure of a Core Promoter in Bifidobacterium longum NCC2705

Bacterial promoters consist of core sequence motifs termed -35 and -10 boxes. The consensus motifs are TTGACA and TATAAT, respectively, which were identified from leading investigations on Escherichia coli However, the consensus sequences are not likely to fit genetically divergent bacteria. The sigma factor of the genus Bifidobacterium has a characteristic polar domain in the N terminus, suggesting the possibility of specific promoter recognition. We reevaluated the structure of Bifidobacterium longum NCC2705 promoters and compared them to other bacteria. Transcriptional start sites (TSSs) of the B. longum NCC2705 strain were identified using transcriptome sequencing (RNA-Seq) analysis to extract promoter regions. Conserved motifs of a bifidobacterial promoter were determined using regions upstream of TSSs and a hidden Markov model. As a result, consensus motifs of the -35 and -10 boxes were TTGTGC and TACAAT, respectively. To assess each base of both motifs, we constructed 37 plasmids based on pKO403-TPCTcon, including the hup promoter connected with a chloramphenicol acetyltransferase as a reporter gene. This reporter assay showed two optimal motifs of the -35 and -10 boxes, namely, TTGNNN and TANNNT, respectively. We further analyzed spacer lengths between the -35 and -10 boxes via a bioinformatics approach. The spacer lengths predominant in bacteria have been generally reported to be approximately 17 bp. In contrast, the predominant spacer lengths in the genus Bifidobacterium and related species were 11 bp, in addition to 17 bp. A reporter assay to assess the spacer lengths indicated that the 11-bp spacer length produced unusually high activity.IMPORTANCE The structures of sigma factors vary among bacterial strains, indicating that recognition rules may also vary. Therefore, we investigated the promoter structure of Bifidobacterium longum NCC2705 using a bioinformatics approach and wet analyses. The most frequent and optimal motifs were similar to other bacterial consensus motifs. The optimal spacer length between the two boxes was reported to be 17 bp. It is widely applied to a bioinformatics approach for other bacteria. Unexpectedly, conserved spacer lengths were 11 bp as well as 17 bp in the genus Bifidobacterium Moreover, the sigma factor of the genus Bifidobacterium has a characteristic domain in the N terminus which may contribute to the additional functions. Hence, it would be valuable to reevaluate the promoter in other organisms.

Transcription attenuation-derived small RNA rnTrpL regulates tryptophan biosynthesis gene expression in trans

Ribosome-mediated transcription attenuation is a basic posttranscriptional regulation mechanism in bacteria. Liberated attenuator RNAs arising in this process are generally considered nonfunctional. In Sinorhizobium meliloti, the tryptophan (Trp) biosynthesis genes are organized into three operons, trpE(G), ppiD-trpDC-moaC-moeA, and trpFBA-accD-folC, of which only the first one, trpE(G), contains a short ORF (trpL) in the 5′-UTR and is regulated by transcription attenuation. Under conditions of Trp sufficiency, transcription is terminated between trpL and trpE(G), and a small attenuator RNA, rnTrpL, is produced. Here, we show that rnTrpL base-pairs with trpD and destabilizes the polycistronic trpDC mRNA, indicating rnTrpL-mediated downregulation of the trpDC operon in trans. Although all three trp operons are regulated in response to Trp availability, only in the two operons trpE(G) and trpDC the Trp-mediated regulation is controlled by rnTrpL. Together, our data show that the trp attenuator coordinates trpE(G) and trpDC expression posttranscriptionally by two fundamentally different mechanisms: ribosome-mediated transcription attenuation in cis and base-pairing in trans. Also, we present evidence that rnTrpL-mediated regulation of trpDC genes expression in trans is conserved in Agrobacterium and Bradyrhizobium, suggesting that the small attenuator RNAs may have additional conserved functions in the control of bacterial gene expression.

baerhunter: an R package for the discovery and analysis of expressed non-coding regions in bacterial RNA-seq data

SUMMARY

Standard bioinformatics pipelines for the analysis of bacterial transcriptomic data commonly ignore non-coding but functional elements e.g. small RNAs, long antisense RNAs or untranslated regions (UTRs) of mRNA transcripts. The root of this problem is the use of incomplete genome annotation files. Here, we present baerhunter, a coverage-based method implemented in R, that automates the discovery of expressed non-coding RNAs and UTRs from RNA-seq reads mapped to a reference genome. The core algorithm is part of a pipeline that facilitates downstream analysis of both coding and non-coding features. The method is simple, easy to extend and customize and, in limited tests with simulated and real data, compares favourably against the currently most popular alternative.

AVAILABILITY AND IMPLEMENTATION

The baerhunter R package is available from: https://github.com/irilenia/baerhunter.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Functional Characterization of a Global Virulence Regulator Hfq and Identification of Hfq-Dependent sRNAs in the Plant Pathogen Pantoea ananatis

To successfully infect plant hosts, the collective regulation of virulence factors in a bacterial pathogen is crucial. Hfq is an RNA chaperone protein that facilitates the small RNA (sRNA) regulation of global gene expression at the post-transcriptional level. In this study, the functional role of Hfq in a broad host range phytopathogen Pantoea ananatis was determined. Inactivation of the hfq gene in P. ananatis LMG 2665^T resulted in the loss of pathogenicity and motility. In addition, there was a significant reduction of quorum sensing signal molecule acyl-homoserine lactone (AHL) production and biofilm formation. Differential sRNA expression analysis between the hfq mutant and wild-type strains of P. ananatis revealed 276 sRNAs affected in their abundance by the loss of hfq at low (OD₆₀₀ = 0.2) and high cell (OD₆₀₀ = 0.6) densities. Further analysis identified 25 Hfq-dependent sRNAs, all showing a predicted Rho-independent terminator of transcription and mapping within intergenic regions of the P. ananatis genome. These included known sRNAs such as ArcZ, FnrS, GlmZ, RprA, RyeB, RyhB, RyhB2, Spot42, and SsrA, and 16 novel P. ananatis sRNAs. The current study demonstrated that Hfq is an important component of the collective regulation of virulence factors and sets a foundation for understanding Hfq-sRNA mediated regulation in the phytopathogen P. ananatis.

The RNase YbeY Is Vital for Ribosome Maturation, Stress Resistance, and Virulence of the Natural Genetic Engineer Agrobacterium tumefaciens

Riboregulation involving regulatory RNAs, RNA chaperones, and ribonucleases is fundamental for the rapid adaptation of gene expression to changing environmental conditions. The gene coding for the RNase YbeY belongs to the minimal prokaryotic genome set and has a profound impact on physiology in a wide range of bacteria. Here, we show that the Agrobacterium tumefaciens ybeY gene is not essential. Deletion of the gene in the plant pathogen reduced growth, motility, and stress tolerance. Most interestingly, YbeY is crucial for A. tumefaciens-mediated T-DNA transfer and tumor formation. Comparative proteomics by using isobaric tags for relative and absolute quantitation (iTRAQ) revealed dysregulation of 59 proteins, many of which have previously been found to be dependent on the RNA chaperone Hfq. YbeY and Hfq have opposing effects on production of these proteins. Accumulation of a 16S rRNA precursor in the ybeY mutant suggests that A. tumefaciens YbeY is involved in rRNA processing. RNA coimmunoprecipitation-sequencing (RIP-Seq) showed binding of YbeY to the region immediately upstream of the 16S rRNA. Purified YbeY is an oligomer with RNase activity. It does not physically interact with Hfq and thus plays a partially overlapping but distinct role in the riboregulatory network of the plant pathogen.IMPORTANCE Although ybeY gene belongs to the universal bacterial core genome, its biological function is incompletely understood. Here, we show that YbeY is critical for fitness and host-microbe interaction in the plant pathogen Agrobacterium tumefaciens Consistent with the reported endoribonuclease activity of YbeY, A. tumefaciens YbeY acts as a RNase involved in maturation of 16S rRNA. This report adds a worldwide plant pathogen and natural genetic engineer of plants to the growing list of bacteria that require the conserved YbeY protein for host-microbe interaction.

Prioritizing bona fide bacterial small RNAs with machine learning classifiers

Bacterial small (sRNAs) are involved in the control of several cellular processes. Hundreds of putative sRNAs have been identified in many bacterial species through RNA sequencing. The existence of putative sRNAs is usually validated by Northern blot analysis. However, the large amount of novel putative sRNAs reported in the literature makes it impractical to validate each of them in the wet lab. In this work, we applied five machine learning approaches to construct twenty models to discriminate bona fide sRNAs from random genomic sequences in five bacterial species. Sequences were represented using seven features including free energy of their predicted secondary structure, their distances to the closest predicted promoter site and Rho-independent terminator, and their distance to the closest open reading frames (ORFs). To automatically calculate these features, we developed an sRNA Characterization Pipeline (sRNACharP). All seven features used in the classification task contributed positively to the performance of the predictive models. The best performing model obtained a median precision of 100% at 10% recall and of 64% at 40% recall across all five bacterial species, and it outperformed previous published approaches on two benchmark datasets in terms of precision and recall. Our results indicate that even though there is limited sRNA sequence conservation across different bacterial species, there are intrinsic features in the genomic context of sRNAs that are conserved across taxa. We show that these features are utilized by machine learning approaches to learn a species-independent model to prioritize bona fide bacterial sRNAs.

Prediction and identification of novel sRNAs involved in Agrobacterium strains by integrated genome-wide and transcriptome-based methods

Small RNAs (sRNAs) are a class of gene regulators in bacteria, playing a central role in their response to environmental changes. Bioinformatic prediction facilitates the identification of sRNAs expressed under different conditions. We propose a novel method of prediction of sRNAs from the genome of Agrobacterium based on a positional weight matrix of conditional sigma factors. sRNAs predicted from the genome are integrated with the virulence-specific transcriptome data to identify putative sRNAs that are overexpressed during Agrobacterial virulence induction. A total of 384 sRNAs are predicted from transcriptome data analysis of Agrobacterium fabrum and 100-500 sRNAs from the genome of different Agrobacterial strains. In order to refine our study, a final set of 10 novel sRNAs with best features across different replicons targeting virulence genes were experimentally identified using semi-quantitative polymerase chain reaction. Since Ti plasmid plays a major role in virulence, out of 10 sRNAs across the replicons, 4 novel sRNAs differentially expressed under virulence induced and non-induced conditions are predicted to be present in the Ti plasmid T-DNA region flanking virulence-related genes like agrocinopine synthase, indole 3-lactate synthase, mannopine synthase and tryptophan monooxygenase. Further validation of the function of these sRNAs in conferring virulence would be relevant to explore their role in Agrobacterium-mediated plant transformation.

A Small Regulatory RNA Controls Cell Wall Biosynthesis and Antibiotic Resistance

Small regulatory RNAs play an important role in the adaptation to changing conditions. Here, we describe a differentially expressed small regulatory RNA (sRNA) that affects various cellular processes in the plant pathogen Agrobacterium tumefaciens Using a combination of bioinformatic predictions and comparative proteomics, we identified nine targets, most of which are positively regulated by the sRNA. According to these targets, we named the sRNA PmaR for peptidoglycan biosynthesis, motility, and ampicillin resistance regulator. Agrobacterium spp. are long known to be naturally resistant to high ampicillin concentrations, and we can now explain this phenotype by the positive PmaR-mediated regulation of the beta-lactamase gene ampC Structure probing revealed a spoon-like structure of the sRNA, with a single-stranded loop that is engaged in target interaction in vivo and in vitro Several riboregulators have been implicated in antibiotic resistance mechanisms, such as uptake and efflux transporters, but PmaR represents the first example of an sRNA that directly controls the expression of an antibiotic resistance gene.IMPORTANCE The alphaproteobacterium Agrobacterium tumefaciens is able to infect various eudicots causing crown gall tumor formation. Based on its unique ability of interkingdom gene transfer, Agrobacterium serves as a crucial biotechnological tool for genetic manipulation of plant cells. The presence of hundreds of putative sRNAs in this organism suggests a considerable impact of riboregulation on A. tumefaciens physiology. Here, we characterized the biological function of the sRNA PmaR that controls various processes crucial for growth, motility, and virulence. Among the genes directly targeted by PmaR is ampC coding for a beta-lactamase that confers ampicillin resistance, suggesting that the sRNA is crucial for fitness in the competitive microbial composition of the rhizosphere.

Effective removal of a range of Ti/Ri plasmids using a pBBR1-type vector having a repABC operon and a lux reporter system

Ti and Ri plasmids of pathogenic Agrobacterium strains are stably maintained by the function of a repABC operon and have been classified into four incompatibility groups, namely, incRh1, incRh2, incRh3, and incRh4. Removal of these plasmids from their bacterial cells is an important step in determining strain-specific virulence characteristics and to construct strains useful for transformation. Here, we developed two powerful tools to improve this process. We first established a reporter system to detect the presence and absence of Ti/Ri plasmids in cells by using an acetosyringone (AS)-inducible promoter of the Ti2 small RNA and luxAB from Vibrio harveyi. This system distinguished a Ti/Ri plasmid-free cell colony among plasmid-harboring cell colonies by causing the latter colonies to emit light in response to AS. We then constructed new "Ti/Ri eviction plasmids," each of which carries a repABC from one of four Ti/Ri plasmids that belonged to incRh1, incRh2, incRh3, and incRh4 groups in the suicidal plasmid pK18mobsacB and in a broad-host-range pBBR1 vector. Introduction of the new eviction plasmids into Agrobacterium cells harboring the corresponding Ti/Ri plasmids led to Ti/Ri plasmid-free cells in every incRh group. The Ti/Ri eviction was more effective by plasmids with the pBBR1 backbone than by those with the pK18mobsacB backbone. Furthermore, the highly stable cryptic plasmid pAtC58 in A. tumefaciens C58 was effectively evicted by the introduction of a pBBR1 vector containing the repABC of pAtC58. These results indicate that the set of pBBR1-repABC plasmids is a powerful tool for the removal of stable rhizobial plasmids.

A Genome-Wide Prediction and Identification of Intergenic Small RNAs by Comparative Analysis in Mesorhizobium huakuii 7653R

In bacteria, small non-coding RNAs (sRNAs) are critical regulators of cellular adaptation to changes in metabolism, physiology, or the external environment. In the last decade, more than 2000 of sRNA families have been reported in the Rfam database and have been shown to exert various regulatory functions in bacterial transcription and translation. However, little is known about sRNAs and their functions in Mesorhizobium. Here, we predicted putative sRNAs in the intergenic regions (IGRs) of M. huakuii 7653R by genome-wide comparisons with four related Mesorhizobial strains. The expression and transcribed regions of candidate sRNAs were analyzed using a set of high-throughput RNA deep sequencing data. In all, 39 candidate sRNAs were found, with 5 located in the symbiotic megaplasmids and 34 in the chromosome of M. huakuii 7653R. Of these, 24 were annotated as functional sRNAs in the Rfam database and 15 were recognized as putative novel sRNAs. The expression of nine selected sRNAs was confirmed by Northern blotting, and most of the nine selected sRNAs were highly expressed in 28 dpi nodules and under symbiosis-mimicking conditions. For those putative novel sRNAs, functional categorizations of their target genes were performed by analyzing the enriched GO terms. In addition, MH_s15 was shown to be an abundant and conserved sRNA.

Small RNAs in Bacterial Virulence and Communication

Bacterial pathogens must endure or adapt to different environments and stresses during transmission and infection. Posttranscriptional gene expression control by regulatory RNAs, such as small RNAs and riboswitches, is now considered central to adaptation in many bacteria, including pathogens. The study of RNA-based regulation (riboregulation) in pathogenic species has provided novel insight into how these bacteria regulate virulence gene expression. It has also uncovered diverse mechanisms by which bacterial small RNAs, in general, globally control gene expression. Riboregulators as well as their targets may also prove to be alternative targets or provide new strategies for antimicrobials. In this article, we present an overview of the general mechanisms that bacteria use to regulate with RNA, focusing on examples from pathogens. In addition, we also briefly review how deep sequencing approaches have aided in opening new perspectives in small RNA identification and the study of their functions. Finally, we discuss examples of riboregulators in two model pathogens that control virulence factor expression or survival-associated phenotypes, such as stress tolerance, biofilm formation, or cell-cell communication, to illustrate how riboregulation factors into regulatory networks in bacterial pathogens.

Genome-Wide Detection of Small Regulatory RNAs in Deep-Sea Bacterium Shewanella piezotolerans WP3

Shewanella are one of the most abundant Proteobacteria in the deep-sea and are renowned for their versatile electron accepting capacities. The molecular mechanisms involved in their adaptation to diverse and extreme environments are not well understood. Small non-coding RNAs (sRNAs) are known for modulating the gene expression at transcriptional and posttranscriptional levels, subsequently playing a key role in microbial adaptation. To understand the potential roles of sRNAs in the adaptation of Shewanella toward deep-sea environments, here an in silico approach was utilized to detect the sRNAs in the genome of Shewanella piezotolerans WP3, a piezotolerant and psychrotolerant deep-sea iron reducing bacterium. After scanning 3673 sets of 5' and 3' UTRs of orthologous genes, 209 sRNA candidates were identified with high confidence in S. piezotolerans WP3. About 92% (193 out of 209) of these putative sRNAs belong to the class trans-encoded RNAs, suggesting that trans-regulatory RNAs are the dominant class of sRNAs in S. piezotolerans WP3. The remaining 16 cis-regulatory RNAs were validated through quantitative polymerase chain reaction. Five cis-sRNAs were further shown to act as cold regulated sRNAs. Our study provided additional evidence at the transcriptional level to decipher the microbial adaptation mechanisms to extreme environmental conditions.

Host-Pathogen interactions modulated by small RNAs

Biological processes such as defense mechanisms and microbial offense strategies are regulated through RNA induced interference in eukaryotes. Genetic mutations are modulated through biogenesis of small RNAs which directly impacts upon host development. Plant defense mechanisms are regulated and supported by a diversified group of small RNAs which are involved in streamlining several RNA interference pathways leading toward the initiation of pathogen gene silencing mechanisms. In the similar context, pathogens also utilize the support of small RNAs to launch their offensive attacks. Also there are strong evidences about the active involvement of these RNAs in symbiotic associations. Interestingly, small RNAs are not limited to the individuals in whom they are produced; they also show cross kingdom influences through variable interactions with other species thus leading toward the inter-organismic gene silencing. The phenomenon is understandable in the microbes which utilize these mechanisms to overcome host defense line. Understanding the mechanism of triggering host defense strategies can be a valuable step toward the generation of disease resistant host plants. We think that the cross kingdom trafficking of small RNA is an interesting insight that is needed to be explored for its vitality.

Genome-wide identification and characterization of small RNAs in Rhodobacter capsulatus and identification of small RNAs affected by loss of the response regulator CtrA

Small non-coding RNAs (sRNAs) are involved in the control of numerous cellular processes through various regulatory mechanisms, and in the past decade many studies have identified sRNAs in a multitude of bacterial species using RNA sequencing (RNA-seq). Here, we present the first genome-wide analysis of sRNA sequencing data in Rhodobacter capsulatus, a purple nonsulfur photosynthetic alphaproteobacterium. Using a recently developed bioinformatics approach, sRNA-Detect, we detected 422 putative sRNAs from R. capsulatus RNA-seq data. Based on their sequence similarity to sRNAs in a sRNA collection, consisting of published putative sRNAs from 23 additional bacterial species, and RNA databases, the sequences of 124 putative sRNAs were conserved in at least one other bacterial species; and, 19 putative sRNAs were assigned a predicted function. We bioinformatically characterized all putative sRNAs and applied machine learning approaches to calculate the probability of a nucleotide sequence to be a bona fide sRNA. The resulting quantitative model was able to correctly classify 95.2% of sequences in a validation set. We found that putative cis-targets for antisense and partially overlapping sRNAs were enriched with protein-coding genes involved in primary metabolic processes, photosynthesis, compound binding, and with genes forming part of macromolecular complexes. We performed differential expression analysis to compare the wild type strain to a mutant lacking the response regulator CtrA, an important regulator of gene expression in R. capsulatus, and identified 18 putative sRNAs with differing levels in the two strains. Finally, we validated the existence and expression patterns of four novel sRNAs by Northern blot analysis.

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.

Perspectives on the Transition From Bacterial Phytopathogen Genomics Studies to Applications Enhancing Disease Management: From Promise to Practice

The advent of genomics has advanced science into a new era, providing a plethora of "toys" for researchers in many related and disparate fields. Genomics has also spawned many new fields, including proteomics and metabolomics, furthering our ability to gain a more comprehensive view of individual organisms and of interacting organisms. Genomic information of both bacterial pathogens and their hosts has provided the critical starting point in understanding the molecular bases of how pathogens disrupt host cells to cause disease. In addition, knowledge of the complete genome sequence of the pathogen provides a potentially broad slate of targets for the development of novel virulence inhibitors that are desperately needed for disease management. Regarding plant bacterial pathogens and disease management, the potential for utilizing genomics resources in the development of durable resistance is enhanced because of developing technologies that enable targeted modification of the host. Here, we summarize the role of genomics studies in furthering efforts to manage bacterial plant diseases and highlight novel genomics-enabled strategies heading down this path.

Isolation and characterization of active promoters from Gluconacetobacter diazotrophicus strain PAL5 using a promoter-trapping plasmid

Gluconacetobacter diazotrophicus is a nitrogen-fixing, endophytic bacterium that has the potential to promote plant growth and increase yield. Genetically modified strains might get more benefits to host plants, including through expression of useful proteins, such as Cry toxins from B. thuringiensis, or enzymes involved in phytohormone production, proteins with antagonistic activity for phytopathogens, or that improve nutrient utilization by the plant. For that, expression systems for G. diazotrophicus are needed, which requires active promoters fused to foreign (or innate) genes. This article describes the construction of a G. diazotrophicus PAL5 promoter library using a promoter-less lacZ-bearing vector, and the identification of six active promoters through β-galactosidase activity assays, sequencing and localization in the bacterial genome. The characterized promoters, which are located on distinct regions of the bacterial genome and encoding either sense or antisense transcripts, present variable expression strengths and might be used in the future for expressing useful proteins.

Discovery and profiling of small RNAs responsive to stress conditions in the plant pathogen Pectobacterium atrosepticum

Background

Small RNAs (sRNAs) have emerged as important regulatory molecules and have been studied in several bacteria. However, to date, there have been no whole-transcriptome studies on sRNAs in any of the Soft Rot Enterobacteriaceae (SRE) group of pathogens. Although the main ecological niches for these pathogens are plants, a significant part of their life cycle is undertaken outside their host within adverse soil environment. However, the mechanisms of SRE adaptation to this harsh nutrient-deficient environment are poorly understood.

Results

In the study reported herein, by using strand-specific RNA-seq analysis and in silico sRNA predictions, we describe the sRNA pool of Pectobacterium atrosepticum and reveal numerous sRNA candidates, including those that are induced during starvation-activated stress responses. Consequently, strand-specific RNA-seq enabled detection of 137 sRNAs and sRNA candidates under starvation conditions; 25 of these sRNAs were predicted for this bacterium in silico. Functional annotations were computationally assigned to 68 sRNAs. The expression of sRNAs in P. atrosepticum was compared under growth-promoting and starvation conditions: 68 sRNAs were differentially expressed with 47 sRNAs up-regulated under nutrient-deficient conditions. Conservation analysis using BLAST showed that most of the identified sRNAs are conserved within the SRE. Subsequently, we identified 9 novel sRNAs within the P. atrosepticum genome.

Conclusions

Since many of the identified sRNAs are starvation-induced, the results of our study suggests that sRNAs play key roles in bacterial adaptive response. Finally, this work provides a basis for future experimental characterization and validation of sRNAs in plant pathogens.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2376-0) contains supplementary material, which is available to authorized users.

Thermodynamic matchers for the construction of the cuckoo RNA family

RNA family models describe classes of functionally related, non-coding RNAs based on sequence and structure conservation. The most important method for modeling RNA families is the use of covariance models, which are stochastic models that serve in the discovery of yet unknown, homologous RNAs. However, the performance of covariance models in finding remote homologs is poor for RNA families with high sequence conservation, while for families with high structure but low sequence conservation, these models are difficult to built in the first place. A complementary approach to RNA family modeling involves the use of thermodynamic matchers. Thermodynamic matchers are RNA folding programs, based on the established thermodynamic model, but tailored to a specific structural motif. As thermodynamic matchers focus on structure and folding energy, they unfold their potential in discovering homologs, when high structure conservation is paired with low sequence conservation. In contrast to covariance models, construction of thermodynamic matchers does not require an input alignment, but requires human design decisions and experimentation, and hence, model construction is more laborious. Here we report a case study on an RNA family that was constructed by means of thermodynamic matchers. It starts from a set of known but structurally different members of the same RNA family. The consensus secondary structure of this family consists of 2 to 4 adjacent hairpins. Each hairpin loop carries the same motif, CCUCCUCCC, while the stems show high variability in their nucleotide content. The present study describes (1) a novel approach for the integration of the structurally varying family into a single RNA family model by means of the thermodynamic matcher methodology, and (2) provides the results of homology searches that were conducted with this model in a wide spectrum of bacterial species.

Expression of Agrobacterium Homolog Genes Encoding T-complex Recruiting Protein under Virulence Induction Conditions

The proteins encoded by three Agrobacterial genes, atu5117, atu4860, and atu4856, are highly homologous to each other in amino acid sequence. All three proteins can bind to VirD2 and are named VBP1, VBP2, and VBP3 (VirD2-binding protein), respectively. VBP is involved in T-DNA transfer by recruiting the T-complex from the cytosol to the polar transport apparatus T4SS (type IVsecretion system) and is defined as the “T-complex recruiting protein.” However, it remains unknown how these three homologous genes co-exist in a relatively small prokaryotic genome. To understand whether these three homologous genes are expressed differentially under virulence induction conditions, we examined the effects of virulence induction conditions, including various pH values, temperatures and acetosyringone (AS, an effective virulence inducer to Agrobacterium tumefaciens) concentrations, on the expression of the three VBP-encoding genes. Our data showed that vbp1 (atu5117) and vbp3 (atu4856) maintained constant expression under the tested induction conditions, whereas the expression of vbp2 (atu4860) was affected by the conditions. Culture conditions favorable to the expression of vbp2 differed from the reported induction conditions for other virulence proteins. In particular, the pH value was a crucial factor for the expression of vbp2. In addition, the deletion of vbp1 affected the expression of vbp2. Taken together, these results suggest that the mechanisms regulating the expression of these three homologous genes are different from the virulence induction mechanism and that VBP homologs are presumably involved in other biological processes in addition to T-complex recruitment.

Small RNAs--the secret agents in the plant-pathogen interactions

Eukaryotic regulatory small RNAs (sRNAs) that induce RNA interference (RNAi) are involved in a plethora of biological processes, including host immunity and pathogen virulence. In plants, diverse classes of sRNAs contribute to the regulation of host innate immunity. These immune-regulatory sRNAs operate through distinct RNAi pathways that trigger transcriptional or post-transcriptional gene silencing. Similarly, many pathogen-derived sRNAs also regulate pathogen virulence. Remarkably, the influence of regulatory sRNAs is not limited to the individual organism in which they are generated. It can sometimes extend to interacting species from even different kingdoms. There they trigger gene silencing in the interacting organism, a phenomenon called cross-kingdom RNAi. This is exhibited in advanced pathogens and parasites that produce sRNAs to suppress host immunity. Conversely, in host-induced gene silencing (HIGS), diverse plants are engineered to trigger RNAi against pathogens and pests to confer host resistance. Cross-kingdom RNAi opens up a vastly unexplored area of research on mobile sRNAs in the battlefield between hosts and pathogens.

Small RNA Deep-Sequencing Analyses Reveal a New Regulator of Virulence in Agrobacterium fabrum C58

Novel ways of regulating Ti plasmid functions were investigated by studying small RNAs (sRNAs) that are known to act as posttranscriptional regulators in plant pathogenic bacteria. sRNA-seq analyses of Agrobacterium fabrum C58 allowed us to identify 1,108 small transcripts expressed in several growth conditions that could be sRNAs. A quarter of them were confirmed by bioinformatics or by biological experiments. Antisense RNAs represent 24% of the candidates and they are over-represented on the pTi (with 62% of pTi sRNAs), suggesting differences in the regulatory mechanisms between the essential and accessory replicons. Moreover, a large number of these pTi antisense RNAs are transcribed opposite to those genes involved in virulence. Others are 5'- and 3'-untranslated region RNAs and trans-encoded RNAs. We have validated, by rapid amplification of cDNA ends polymerase chain reaction, the transcription of 14 trans-encoded RNAs, among which RNA1111 is expressed from the pTiC58. Its deletion decreased the aggressiveness of A. fabrum C58 on tomatoes, tobaccos, and kalanchoe, suggesting that this sRNA activates virulence. The identification of its putative target mRNAs (6b gene, virC2, virD3, and traA) suggests that this sRNA may coordinate two of the major pTi functions, the infection of plants and its dissemination among bacteria.

Genome-wide analyses in bacteria show small-RNA enrichment for long and conserved intergenic regions

Interest in finding small RNAs (sRNAs) in bacteria has significantly increased in recent years due to their regulatory functions. Development of high-throughput methods and more sophisticated computational algorithms has allowed rapid identification of sRNA candidates in different species. However, given their various sizes (50 to 500 nucleotides [nt]) and their potential genomic locations in the 5' and 3' untranslated regions as well as in intergenic regions, identification and validation of true sRNAs have been challenging. In addition, the evolution of bacterial sRNAs across different species continues to be puzzling, given that they can exert similar functions with various sequences and structures. In this study, we analyzed the enrichment patterns of sRNAs in 13 well-annotated bacterial species using existing transcriptome and experimental data. All intergenic regions were analyzed by WU-BLAST to examine conservation levels relative to species within or outside their genus. In total, more than 900 validated bacterial sRNAs and 23,000 intergenic regions were analyzed. The results indicate that sRNAs are enriched in intergenic regions, which are longer and more conserved than the average intergenic regions in the corresponding bacterial genome. We also found that sRNA-coding regions have different conservation levels relative to their flanking regions. This work provides a way to analyze how noncoding RNAs are distributed in bacterial genomes and also shows conserved features of intergenic regions that encode sRNAs. These results also provide insight into the functions of regions surrounding sRNAs and into optimization of RNA search algorithms.

Non-coding RNA regulation in pathogenic bacteria located inside eukaryotic cells

Intracellular bacterial pathogens have evolved distinct lifestyles inside eukaryotic cells. Some pathogens coexist with the infected cell in an obligate intracellular state, whereas others transit between the extracellular and intracellular environment. Adaptation to these intracellular lifestyles is regulated in both space and time. Non-coding small RNAs (sRNAs) are post-transcriptional regulatory molecules that fine-tune important processes in bacterial physiology including cell envelope architecture, intermediate metabolism, bacterial communication, biofilm formation, and virulence. Recent studies have shown production of defined sRNA species by intracellular bacteria located inside eukaryotic cells. The molecules targeted by these sRNAs and their expression dynamics along the intracellular infection cycle remain, however, poorly characterized. Technical difficulties linked to the isolation of “intact” intracellular bacteria from infected host cells might explain why sRNA regulation in these specialized pathogens is still a largely unexplored field. Transition from the extracellular to the intracellular lifestyle provides an ideal scenario in which regulatory sRNAs are intended to participate; so much work must be done in this direction. This review focuses on sRNAs expressed by intracellular bacterial pathogens during the infection of eukaryotic cells, strategies used with these pathogens to identify sRNAs required for virulence, and the experimental technical challenges associated to this type of studies. We also discuss varied techniques for their potential application to study RNA regulation in intracellular bacterial infections.

Profound impact of Hfq on nutrient acquisition, metabolism and motility in the plant pathogen Agrobacterium tumefaciens

As matchmaker between mRNA and sRNA interactions, the RNA chaperone Hfq plays a key role in riboregulation of many bacteria. Often, the global influence of Hfq on the transcriptome is reflected by substantially altered proteomes and pleiotropic phenotypes in hfq mutants. Using quantitative proteomics and co-immunoprecipitation combined with RNA-sequencing (RIP-seq) of Hfq-bound RNAs, we demonstrate the pervasive role of Hfq in nutrient acquisition, metabolism and motility of the plant pathogen Agrobacterium tumefaciens. 136 of 2544 proteins identified by iTRAQ (isobaric tags for relative and absolute quantitation) were affected in the absence of Hfq. Most of them were associated with ABC transporters, general metabolism and motility. RIP-seq of chromosomally encoded Hfq^3xFlag revealed 1697 mRNAs and 209 non-coding RNAs (ncRNAs) associated with Hfq. 56 ncRNAs were previously undescribed. Interestingly, 55% of the Hfq-bound ncRNAs were encoded antisense (as) to a protein-coding sequence suggesting that A. tumefaciens Hfq plays an important role in asRNA-target interactions. The exclusive enrichment of 296 mRNAs and 31 ncRNAs under virulence conditions further indicates a role for post-transcriptional regulation in A. tumefaciens-mediated plant infection. On the basis of the iTRAQ and RIP-seq data, we assembled a comprehensive model of the Hfq core regulon in A. tumefaciens.

Differential RNA-seq: the approach behind and the biological insight gained

RNA-sequencing has revolutionized the quantitative and qualitative analysis of transcriptomes in both prokaryotes and eukaryotes. It provides a generic approach for gene expression profiling, annotation of transcript boundaries and operons, as well as identifying novel transcripts including small noncoding RNA molecules and antisense RNAs. We recently developed a differential RNA-seq (dRNA-seq) method which in addition to the above, yields information as to whether a given RNA is a primary or processed transcript. Originally applied to describe the primary transcriptome of the gastric pathogen Helicobacter pylori, dRNA-seq has since provided global maps of transcriptional start sites in diverse species, informed new biology in the CRISPR-Cas9 system, advanced to a tool for comparative transcriptomics, and inspired simultaneous RNA-seq of pathogen and host.

Riboregulation in plant-associated α-proteobacteria

The symbiotic α-rhizobia Sinorhizobium meliloti, Bradyrhizobium japonicum, Rhizobium etli and the related plant pathogen Agrobacterium tumefaciens are important model organisms for studying plant-microbe interactions. These metabolically versatile soil bacteria are characterized by complex lifestyles and large genomes. Here we summarize the recent knowledge on their small non-coding RNAs (sRNAs) including conservation, function, and interaction of the sRNAs with the RNA chaperone Hfq. In each of these organisms, an inventory of hundreds of cis- and trans-encoded sRNAs with regulatory potential was uncovered by high-throughput approaches and used for the construction of 39 sRNA family models. Genome-wide analyses of hfq mutants and co-immunoprecipitation with tagged Hfq revealed a major impact of the RNA chaperone on the physiology of plant-associated α-proteobacteria including symbiosis and virulence. Highly conserved members of the SmelC411 family are the AbcR sRNAs, which predominantly regulate ABC transport systems. AbcR1 of A. tumefaciens controls the uptake of the plant-generated signaling molecule GABA and is a central regulator of nutrient uptake systems. It has similar functions in S. meliloti and the human pathogen Brucella abortus. As RNA degradation is an important process in RNA-based gene regulation, a short overview on ribonucleases in plant-associated α-proteobacteria concludes this review.

Transcriptome sequencing revealed the transcriptional organization at ribosome-mediated attenuation sites in Corynebacterium glutamicum and identified a novel attenuator involved in aromatic amino acid biosynthesis

The Gram-positive bacterium Corynebacterium glutamicum belongs to the order Corynebacteriales and is used as a producer of amino acids at industrial scales. Due to its economic importance, gene expression and particularly the regulation of amino acid biosynthesis has been investigated extensively. Applying the high-resolution technique of transcriptome sequencing (RNA-seq), recently a vast amount of data has been generated that was used to comprehensively analyze the C. glutamicum transcriptome. By analyzing RNA-seq data from a small RNA cDNA library of C. glutamicum, short transcripts in the known transcriptional attenuators sites of the trp operon, the ilvBNC operon and the leuA gene were verified. Furthermore, whole transcriptome RNA-seq data were used to elucidate the transcriptional organization of these three amino acid biosynthesis operons. In addition, we discovered and analyzed the novel attenuator aroR, located upstream of the aroF gene (cg1129). The DAHP synthase encoded by aroF catalyzes the first step in aromatic amino acid synthesis. The AroR leader peptide contains the amino acid sequence motif F-Y-F, indicating a regulatory effect by phenylalanine and tyrosine. Analysis by real-time RT-PCR suggests that the attenuator regulates the transcription of aroF in dependence of the cellular amount of tRNA loaded with phenylalanine when comparing a phenylalanine-auxotrophic C. glutamicum mutant fed with limiting and excess amounts of a phenylalanine-containing dipeptide. Additionally, the very interesting finding was made that all analyzed attenuators are leaderless transcripts.

Genome-wide identification of Hfq-regulated small RNAs in the fire blight pathogen Erwinia amylovora discovered small RNAs with virulence regulatory function

BACKGROUND

Erwinia amylovora is a phytopathogenic bacterium and causal agent of fire blight disease in apples and pears. Although many virulence factors have been characterized, the coordination of expression of these virulence factors in E. amylovora is still not clear. Regulatory small RNAs (sRNAs) are important post-transcriptional regulatory components in bacteria. A large number of sRNAs require the RNA chaperone Hfq for both stability and functional activation. In E. amylovora, Hfq was identified as a major regulator of virulence and various virulence traits. However, information is still lacking about Hfq-dependent sRNAs on a genome scale, including the virulence regulatory functions of these sRNAs in E. amylovora.

RESULTS

Using both an RNA-seq analysis and a Rho-independent terminator search, 40 candidate Hfq-dependent sRNAs were identified in E. amylovora. The expression and sizes of 12 sRNAs and the sequence boundaries of seven sRNAs were confirmed by Northern blot and 5' RACE assay respectively. Sequence conservation analysis identified sRNAs conserved only in the Erwinia genus as well as E. amylovora species-specific sRNAs. In addition, a dynamic re-patterning of expression of Hfq-dependent sRNAs was observed at 6 and 12 hours after induction in Hrp-inducing minimal medium. Furthermore, sRNAs that control virulence traits were characterized, among which ArcZ positively controls the type III secretion system (T3SS), amylovoran exopolysaccahride production, biofilm formation, and motility, and negatively modulates attachment while RmaA (Hrs6) and OmrAB both negatively regulate amylovoran production and positively regulate motility.

CONCLUSIONS

This study has significantly enhanced our understanding of the Hfq-dependent sRNAs in E. amylovora at the genome level. The identification of multiple virulence-regulating sRNAs also suggests that post-transcriptional regulation by sRNAs may play a role in the deployment of virulence factors needed during varying stages of pathogenesis during host invasion by E. amylovora.

Two separate modules of the conserved regulatory RNA AbcR1 address multiple target mRNAs in and outside of the translation initiation region

The small RNA AbcR1 regulates the expression of ABC transporters in the plant pathogen Agrobacterium tumefaciens, the plant symbiont Sinorhizobium meliloti, and the human pathogen Brucella abortus. A combination of proteomic and bioinformatic approaches suggested dozens of AbcR1 targets in A. tumefaciens. Several of these newly discovered targets are involved in the uptake of amino acids, their derivatives, and sugars. Among the latter is the periplasmic sugar-binding protein ChvE, a component of the virulence signal transduction system. We examined 16 targets and their interaction with AbcR1 in close detail. In addition to the previously described mRNA interaction site of AbcR1 (M1), the CopraRNA program predicted a second functional module (M2) as target-binding site. Both M1 and M2 contain single-stranded anti-SD motifs. Using mutated AbcR1 variants, we systematically tested by band shift experiments, which sRNA region is responsible for mRNA binding and gene regulation. On the target site, we find that AbcR1 interacts with some mRNAs in the translation initiation region and with others far into their coding sequence. Our data show that AbcR1 is a versatile master regulator of nutrient uptake systems in A. tumefaciens and related bacteria.

Membrane lipids in Agrobacterium tumefaciens: biosynthetic pathways and importance for pathogenesis

Many cellular processes critically depend on the membrane composition. In this review, we focus on the biosynthesis and physiological roles of membrane lipids in the plant pathogen Agrobacterium tumefaciens. The major components of A. tumefaciens membranes are the phospholipids (PLs), phosphatidylethanolamine (PE), phosphatidylglycerol, phosphatidylcholine (PC) and cardiolipin, and ornithine lipids (OLs). Under phosphate-limited conditions, the membrane composition shifts to phosphate-free lipids like glycolipids, OLs and a betaine lipid. Remarkably, PC and OLs have opposing effects on virulence of A. tumefaciens. OL-lacking A. tumefaciens mutants form tumors on the host plant earlier than the wild type suggesting a reduced host defense response in the absence of OLs. In contrast, A. tumefaciens is compromised in tumor formation in the absence of PC. In general, PC is a rare component of bacterial membranes but amount to ~22% of all PLs in A. tumefaciens. PC biosynthesis occurs via two pathways. The phospholipid N-methyltransferase PmtA methylates PE via the intermediates monomethyl-PE and dimethyl-PE to PC. In the second pathway, the membrane-integral enzyme PC synthase (Pcs) condenses choline with CDP-diacylglycerol to PC. Apart from the virulence defect, PC-deficient A. tumefaciens pmtA and pcs double mutants show reduced motility, enhanced biofilm formation and increased sensitivity towards detergent and thermal stress. In summary, there is cumulative evidence that the membrane lipid composition of A. tumefaciens is critical for agrobacterial physiology and tumor formation.

Small non-coding RNAs in plant-pathogenic Xanthomonas spp

The genus Xanthomonas comprises a large group of plant-pathogenic bacteria. The infection and bacterial multiplication in the plant tissue depends on the type III secretion system and other virulence determinants. Recent studies revealed that bacterial virulence is also controlled at the post-transcriptional level by small non-coding RNAs (sRNAs). In this review, we highlight our current knowledge about sRNAs and RNA-binding proteins in Xanthomonas species.

Riboregulators and the role of Hfq in photosynthetic bacteria

Anoxygenic and oxygenic bacteria directly convert solar energy into biomass using photosynthesis. The formation and composition of photosynthetic complexes has to be tightly controlled in response to environmental conditions, as exposure to sunlight can be harmful due to the generation of reactive oxygen species and the damaging effects of UV irradiation. Therefore, photosynthetic bacteria are exposed to a particular set of regulatory challenges in addition to those that also affect other bacteria, requiring sophisticated regulatory systems. Indeed, hundreds of potential regulatory RNAs have been identified in photosynthetic model bacteria as well as antisense RNAs (asRNAs) of up to several kb in length that protect certain mRNAs from degradation. The trans-acting small non-coding RNAs (sRNAs), PcrZ and PsrR1, control pigment and photosystem biogenesis in Rhodobacter sphaeroides and cyanobacteria, respectively. The asRNAs IsrR and As1_flv4 act as negative regulators and the asRNAs PsbA2R and PsbA3R as positive effectors of photosynthesis gene expression in Synechocystis 6803.