Genomic insights into a Pseudomonas amygdali isolate from Hibiscus rosa-sinensis.

The taxonomy of Pseudomonas has been extensively studied, yet the determination of species is currently difficult because of recent taxonomic changes and the lack of complete genomic sequence data. We isolated a bacterium causing a leaf spot disease on hibiscus (Hibiscus rosa-sinensis). Whole genome sequencing revealed similarity to Pseudomonas amygdali pv. tabaci and pv. lachrymans. The genome of this isolate (referred to as P. amygdali 35-1) shared 4987 genes with P. amygdali pv. hibisci, but possessed 204 unique genes and contained gene clusters encoding putative secondary metabolites and copper resistance determinants. We predicted this isolate's type III secretion effector (T3SE) repertoire and identified 64 putative T3SEs, some of which are present in other P. amygdali pv. hibisci strains. Assays showed that the isolate was resistant to copper at a concentration of 1.6 mM. This study provides an improved understanding of the genomic relatedness and diversity of the P. amygdali species.

Complete Genome Sequence Resource for Xanthomonas hortorum Isolated from Greek Oregano.

In spring 2019, necrotic leaf spots were detected on Greek oregano (Origanum vulgare var. hirtum) plants in a commercial greenhouse operation. An isolate was recovered from the diseased plants. Partial 16S ribosomal RNA sequencing and multilocus sequence analysis revealed that the isolate was a Xanthomonas sp. but proved insufficient to identify the species with certainty. Therefore, whole-genome sequencing using both Nanopore and Illumina technologies was performed. Here, we report the complete and annotated genome sequence of Xanthomonas hortorum strain 108, which was originally isolated from Greek oregano in Long Island, NY, U.S.A.

Identification of Indole-3-Acetic Acid-Regulated Genes in Pseudomonas syringae pv. tomato Strain DC3000.

The auxin indole-3-acetic acid (IAA) is a plant hormone that not only regulates plant growth and development but also plays important roles in plant-microbe interactions. We previously reported that IAA alters expression of several virulence-related genes in the plant pathogen Pseudomonas syringae pv. tomato strain DC3000 (PtoDC3000). To learn more about the impact of IAA on regulation of PtoDC3000 gene expression, we performed a global transcriptomic analysis of bacteria grown in culture, in the presence or absence of exogenous IAA. We observed that IAA repressed expression of genes involved in the type III secretion (T3S) system and motility and promoted expression of several known and putative transcriptional regulators. Several of these regulators are orthologs of factors known to regulate stress responses and accordingly expression of several stress response-related genes was also upregulated by IAA. Similar trends in expression for several genes were also observed by quantitative reverse transcription PCR. Using an Arabidopsis thaliana auxin receptor mutant that accumulates elevated auxin, we found that many of the P. syringae genes regulated by IAA in vitro were also regulated by auxin in planta. Collectively the data indicate that IAA modulates many aspects of PtoDC3000 biology, presumably to promote both virulence and survival under stressful conditions, including those encountered in or on plant leaves. IMPORTANCE Indole-3-acetic acid (IAA), a form of the plant hormone auxin, is used by many plant-associated bacteria as a cue to sense the plant environment. Previously, we showed that IAA can promote disease in interactions between the plant pathogen Pseudomonas syringae strain PtoDC000 and one of its hosts, Arabidopsis thaliana. However, the mechanisms by which IAA impacts the biology of PtoDC3000 and promotes disease are not well understood. Here, we demonstrate that IAA is a signal molecule that regulates gene expression in PtoDC3000. The presence of exogenous IAA affects expression of over 700 genes in the bacteria, including genes involved in type III secretion and genes involved in stress response. This work offers insight into the roles of auxin-promoting pathogenesis.

Complete Genome Sequence Resource for the Necrotrophic Plant-Pathogenic Bacterium Dickeya dianthicola 67-19 Isolated From New Guinea Impatiens.

New Guinea impatiens (NGI, Impatiens hawkeri) are popular bedding plants that can be affected by a number of pathogens. Using 16S rDNA sequencing and genus-specific PCR, we identified the first Dickeya dianthicola strain isolated from NGI presented with blackleg symptoms, herein designated as D. dianthicola 67-19. Here, we report a high-quality complete and annotated genome sequence of D. dianthicola 67-19. The 4,851,809 bp genome was assembled with Nanopore reads and polished with Illumina reads, yielding 422× and 105× coverage, respectively. This closed genome provides a resource for future research on comparative genomics and biology of D. dianthicola, which could translate to improved detection and disease management.

Complete Genome Sequence Resource for the Necrotrophic Plant-Pathogenic Bacterium Pectobacterium carotovorum WPP14.

Pectobacterium spp. are a major cause of loss in vegetable and ornamental plant production. One of these species, Pectobacterium carotovorum, can cause soft rot disease on many plants, particularly potato. These diseases lead to significant economic loss and pose food security threats by reducing crop yields in the field, in transit, and during storage. The Gram-negative enterobacterium P. carotovorum WPP14 is a particularly virulent strain for which there is no available closed genome, limiting the molecular research for this important pathogen. Here, we report a high-quality complete and annotated genome sequence of P. carotovorum WPP14. The 4,892,225-bp genome was assembled with Nanopore reads and polished with Illumina reads, yielding 394× and 164× coverage, respectively. This closed genome provides a resource for research on improved detection and biology of P. carotovorum, which could translate into improved disease management.

Antibacterial activity and mode of action of potassium tetraborate tetrahydrate against soft-rot bacterial plant pathogens.

Bacterial soft rot caused by the bacteria Dickeya and Pectobacterium is a destructive disease of vegetables, as well as ornamental plants. Several management options exist to help control these pathogens. Because of the limited success of these approaches, there is a need for the development of alternative methods to reduce losses. In this study, we evaluated the effect of potassium tetraborate tetrahydrate (PTB) on the growth of six Dickeya and Pectobacterium spp. Disc diffusion assays showed that Dickeya spp. and Pectobacterium spp. differ in their sensitivity to PTB. Spontaneous PTB-resistant mutants of Pectobacterium were identified and further investigation of the mechanism of PTB resistance was conducted by full genome sequencing. Point mutations in genes cpdB and supK were found in a single Pectobacterium atrosepticum PTB-resistant mutant. Additionally, point mutations in genes prfB (synonym supK) and prmC were found in two independent Pectobacterium brasiliense PTB-resistant mutants. prfB and prmC encode peptide chain release factor 2 and its methyltransferase, respectively. We propose the disruption of translation activity due to PTB leads to Pectobacterium growth inhibition. The P. atrosepticum PTB-resistant mutant showed altered swimming motility. Disease severity was reduced for P. atrosepticum-inoculated potato stems sprayed with PTB. We discuss the potential risk of selecting for bacterial resistance to this chemical.

Prevention of Surface-Associated Calcium Phosphate by the Pseudomonas syringae Two-Component System CvsSR.

CvsSR is a Ca2+-induced two-component system (TCS) in the plant pathogen Pseudomonas syringae pv. tomato DC3000. Here, we discovered that CvsSR is induced by Fe3+, Zn2+, and Cd2+ However, only supplementation of Ca2+ to medium resulted in rugose, opaque colonies in ΔcvsS and ΔcvsR strains. This phenotype corresponded to formation of calcium phosphate precipitation on the surface of ΔcvsS and ΔcvsR colonies. CvsSR regulated swarming motility in P. syringae pv. tomato in a Ca2+-dependent manner, but swarming behavior was not influenced by Fe3+, Zn2+, or Cd2+ We hypothesized that reduced swarming displayed by ΔcvsS and ΔcvsR strains was due to precipitation of calcium phosphate on the surface of ΔcvsS and ΔcvsR cells grown on agar medium supplemented with Ca2+ By reducing the initial pH or adding glucose to the medium, calcium precipitation was inhibited, and swarming was restored to ΔcvsS and ΔcvsR strains, suggesting that calcium precipitation influences swarming ability. Constitutive expression of a CvsSR-regulated carbonic anhydrase and a CvsSR-regulated putative sulfate major facilitator superfamily transporter in ΔcvsS and ΔcvsR strains inhibited formation of calcium precipitates and restored the ability of ΔcvsS and ΔcvsR bacteria to swarm. Lastly, we found that glucose inhibited Ca2+-based induction of CvsSR. Hence, CvsSR is a key regulator that controls calcium precipitation on the surface of bacterial cells.IMPORTANCE Bacteria are capable of precipitating and dissolving minerals. We previously reported the characterization of the two-component system CvsSR in the plant-pathogenic bacterium Pseudomonas syringae CvsSR responds to the presence of calcium and is important for causing disease. Here, we show that CvsSR controls the ability of the bacterium to prevent calcium phosphate precipitation on the surface of cells. We also identified a carbonic anhydrase and transporter that modulate formation of surface-associated calcium precipitates. Furthermore, our results demonstrate that the ability of the bacterium to swarm is controlled by the formation and dissolution of calcium precipitates on the surface of cells. Our study describes new mechanisms for microbially induced mineralization and provides insights into the role of mineral deposits on bacterial physiology. The discoveries may lead to new technological and environmental applications.

Ca2+-Induced Two-Component System CvsSR Regulates the Type III Secretion System and the Extracytoplasmic Function Sigma Factor AlgU in Pseudomonas syringae pv. tomato DC3000.

Two-component systems (TCSs) of bacteria regulate many different aspects of the bacterial life cycle, including pathogenesis. Most TCSs remain uncharacterized, with no information about the signal(s) or regulatory targets and/or role in bacterial pathogenesis. Here, we characterized a TCS in the plant-pathogenic bacterium Pseudomonas syringae pv. tomato DC3000 composed of the histidine kinase CvsS and the response regulator CvsR. CvsSR is necessary for virulence of P. syringae pv. tomato DC3000, since ΔcvsS and ΔcvsR strains produced fewer symptoms than the wild type (WT) and demonstrated reduced growth on multiple hosts. We discovered that expression of cvsSR is induced by Ca2+ concentrations found in leaf apoplastic fluid. Thus, Ca2+ can be added to the list of signals that promote pathogenesis of P. syringae pv. tomato DC3000 during host colonization. Through chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) and global transcriptome analysis (RNA-seq), we discerned the CvsR regulon. CvsR directly activated expression of the type III secretion system regulators, hrpR and hrpS, that regulate P. syringae pv. tomato DC3000 virulence in a type III secretion system-dependent manner. CvsR also indirectly repressed transcription of the extracytoplasmic sigma factor algU and production of alginate. Phenotypic analysis determined that CvsSR inversely regulated biofilm formation, swarming motility, and cellulose production in a Ca2+-dependent manner. Overall, our results show that CvsSR is a key regulatory hub critical for interaction with host plants.IMPORTANCE Pathogenic bacteria must be able to react and respond to the surrounding environment, make use of available resources, and avert or counter host immune responses. Often, these abilities rely on two-component systems (TCSs) composed of interacting proteins that modulate gene expression. We identified a TCS in the plant-pathogenic bacterium Pseudomonas syringae that responds to the presence of calcium, which is an important signal during the plant defense response. We showed that when P. syringae is grown in the presence of calcium, this TCS regulates expression of factors contributing to disease. Overall, our results provide a better understanding of how bacterial pathogens respond to plant signals and control systems necessary for eliciting disease.

Virulence of Pseudomonas syringae pv. tomato DC3000 Is Influenced by the Catabolite Repression Control Protein Crc.

Pseudomonas syringae infects diverse plant species and is widely used as a model system in the study of effector function and the molecular basis of plant diseases. Although the relationship between bacterial metabolism, nutrient acquisition, and virulence has attracted increasing attention in bacterial pathology, it is largely unexplored in P. syringae. The Crc (catabolite repression control) protein is a putative RNA-binding protein that regulates carbon metabolism as well as a number of other factors in the pseudomonads. Here, we show that deletion of crc increased bacterial swarming motility and biofilm formation. The crc mutant showed reduced growth and symptoms in Arabidopsis and tomato when compared with the wild-type strain. We have evidence that the crc mutant shows delayed hypersensitive response (HR) when infiltrated into Nicotiana benthamiana and tobacco. Interestingly, the crc mutant was more susceptible to hydrogen peroxide, suggesting that, in planta, the mutant may be sensitive to reactive oxygen species generated during pathogen-associated molecular pattern-triggered immunity (PTI). Indeed, HR was further delayed when PTI-induced tissues were challenged with the crc mutant. The crc mutant did not elicit an altered PTI response in plants compared with the wild-type strain. We conclude that Crc plays an important role in growth and survival during infection.

Disruption of the carA gene in Pseudomonas syringae results in reduced fitness and alters motility.

BACKGROUND: Pseudomonas syringae infects diverse plant species and is widely used in the study of effector function and the molecular basis of disease. Although the relationship between bacterial metabolism, nutrient acquisition and virulence has attracted increasing attention in bacterial pathology, there is limited knowledge regarding these studies in Pseudomonas syringae. The aim of this study was to investigate the function of the carA gene and the small RNA P32, and characterize the regulation of these transcripts. RESULTS: Disruption of the carA gene (ΔcarA) which encodes the predicted small chain of carbamoylphosphate synthetase, resulted in arginine and pyrimidine auxotrophy in Pseudomonas syringae pv. tomato DC3000. Complementation with the wild type carA gene was able to restore growth to wild-type levels in minimal medium. Deletion of the small RNA P32, which resides immediately upstream of carA, did not result in arginine or pyrimidine auxotrophy. The expression of carA was influenced by the concentrations of both arginine and uracil in the medium. When tested for pathogenicity, ΔcarA showed reduced fitness in tomato as well as Arabidopsis when compared to the wild-type strain. In contrast, mutation of the region encoding P32 had minimal effect in planta. ΔcarA also exhibited reduced motility and increased biofilm formation, whereas disruption of P32 had no impact on motility or biofilm formation. CONCLUSIONS: Our data show that carA plays an important role in providing arginine and uracil for growth of the bacteria and also influences other factors that are potentially important for growth and survival during infection. Although we find that the small RNA P32 and carA are co-transcribed, P32 does not play a role in the phenotypes that carA is required for, such as motility, cell attachment, and virulence. Additionally, our data suggests that pyrimidines may be limited in the apoplastic space of the plant host tomato.

Analysis of the small RNA P16/RgsA in the plant pathogen Pseudomonas syringae pv. tomato strain DC3000.

Bacteria contain small non-coding RNAs (ncRNAs) that are responsible for altering transcription, translation or mRNA stability. ncRNAs are important because they regulate virulence factors and susceptibility to various stresses. Here, the regulation of a recently described ncRNA of Pseudomonas syringae pv. tomato DC3000, P16, was investigated. We determined that RpoS regulates the expression of P16. We found that deletion of P16 results in increased sensitivity to hydrogen peroxide compared to the wild-type strain, suggesting that P16 plays a role in the bacteria's susceptibility to oxidative stress. Additionally the P16 mutant displayed enhanced resistance to heat stress. Our findings provide new information on the regulation and role of this ncRNA in P. syringae.

CrcZ and CrcX regulate carbon source utilization in Pseudomonas syringae pathovar tomato strain DC3000.

Small non-coding RNAs (ncRNAs) are important components of many regulatory pathways in bacteria and play key roles in regulating factors important for virulence. Carbon catabolite repression control is modulated by small RNAs (crcZ or crcZ and crcY) in Pseudomonas aeruginosa and Pseudomonas putida. In this study, we demonstrate that expression of crcZ and crcX (formerly designated psr1 and psr2, respectively) is dependent upon RpoN together with the two-component system CbrAB, and is influenced by the carbon source present in the medium in the model plant pathogen Pseudomonas syringae pv tomato DC3000. The distribution of the members of the Crc ncRNA family was also determined by screening available genomic sequences of the Pseudomonads. Interestingly, variable numbers of the Crc family members exist in Pseudomonas genomes. The ncRNAs are comprised of three main subfamilies, named CrcZ, CrcX and CrcY. Most importantly the CrcX subfamily appears to be unique to all P. syringae strains sequenced to date.

Pseudomonas syringae Type III Secretion Protein HrpP Manipulates Plant Immunity To Promote Infection.

The bacterial plant pathogen Pseudomonas syringae deploys a type III secretion system (T3SS) to deliver effector proteins into plant cells to facilitate infection, for which many effectors have been characterized for their interactions. However, few T3SS Hrp (hypersensitive response and pathogenicity) proteins from the T3SS secretion apparatus have been studied for their direct interactions with plants. Here, we show that the P. syringae pv. tomato DC3000 T3SS protein HrpP induces host cell death, suppresses pattern-triggered immunity (PTI), and restores the effector translocation ability of the hrpP mutant. The hrpP-transgenic Arabidopsis lines exhibited decreased PTI responses to flg22 and elf18 and enhanced disease susceptibility to P. syringae pv. tomato DC3000. Transcriptome analysis reveals that HrpP sensing activates salicylic acid (SA) signaling while suppressing jasmonic acid (JA) signaling, which correlates with increased SA accumulation and decreased JA biosynthesis. Both yeast two-hybrid and bimolecular fluorescence complementation assays show that HrpP interacts with mitogen-activated protein kinase kinase 2 (MKK2) on the plant membrane and in the nucleus. The HrpP truncation HrpP1-119, rather than HrpP1-101, retains the ability to interact with MKK2 and suppress PTI in plants. In contrast, HrpP1-101 continues to cause cell death and electrolyte leakage. MKK2 silencing compromises SA signaling but has no effect on cell death caused by HrpP. Overall, our work highlights that the P. syringae T3SS protein HrpP facilitates effector translocation and manipulates plant immunity to facilitate bacterial infection. IMPORTANCE The T3SS is required for the virulence of many Gram-negative bacterial pathogens of plants and animals. This study focuses on the sensing and function of the T3SS protein HrpP during plant interactions. Our findings show that HrpP and its N-terminal truncation HrpP1-119 can interact with MKK2, promote effector translocation, and manipulate plant immunity to facilitate bacterial infection, highlighting the P. syringae T3SS component involved in the fine-tuning of plant immunity.

Genome-Wide Identification of Genes Important for Growth of Dickeya dadantii and Dickeya dianthicola in Potato (Solanum tuberosum) Tubers.

Dickeya species are causal agents of soft rot diseases in many economically important crops, including soft rot disease of potato (Solanum tuberosum). Using random barcode transposon-site sequencing (RB-TnSeq), we generated genome-wide mutant fitness profiles of Dickeya dadantii 3937, Dickeya dianthicola ME23, and Dickeya dianthicola 67-19 isolates collected after passage through several in vitro and in vivo conditions. Though all three strains are pathogenic on potato, D. dadantii 3937 is a well-characterized model while D. dianthicola strains ME23 and 67-19 are recent isolates. Strain ME23 specifically was identified as a representative strain from a 2014 outbreak on potato. This study generated comparable gene fitness measurements across ecologically relevant conditions for both model and non-model strains. Tubers from the potato cultivars "Atlantic," "Dark Red Norland," and "Upstate Abundance" provided highly similar conditions for bacterial growth. Using the homolog detection software PyParanoid, we matched fitness values for orthologous genes in the three bacterial strains. Direct comparison of fitness among the strains highlighted shared and variable traits important for growth. Bacterial growth in minimal medium required many metabolic traits that were also essential for competitive growth in planta, such as amino acid, carbohydrate, and nucleotide biosynthesis. Growth in tubers specifically required the pectin degradation gene kduD. Disruption in three putative DNA-binding proteins had strain-specific effects on competitive fitness in tubers. Though the Soft Rot Pectobacteriaceae can cause disease with little host specificity, it remains to be seen the extent to which strain-level variation impacts virulence.

Characterization of the Fur regulon in Pseudomonas syringae pv. tomato DC3000.

The plant pathogen Pseudomonas syringae pv. tomato DC3000 (DC3000) is found in a wide variety of environments and must monitor and respond to various environmental signals such as the availability of iron, an essential element for bacterial growth. An important regulator of iron homeostasis is Fur (ferric uptake regulator), and here we present the first study of the Fur regulon in DC3000. Using chromatin immunoprecipitation followed by massively parallel sequencing (ChIP-seq), 312 chromosomal regions were highly enriched by coimmunoprecipitation with a C-terminally tagged Fur protein. Integration of these data with previous microarray and global transcriptome analyses allowed us to expand the putative DC3000 Fur regulon to include genes both repressed and activated in the presence of bioavailable iron. Using nonradioactive DNase I footprinting, we confirmed Fur binding in 41 regions, including upstream of 11 iron-repressed genes and the iron-activated genes encoding two bacterioferritins (PSPTO_0653 and PSPTO_4160), a ParA protein (PSPTO_0855), and a two-component system (TCS) (PSPTO_3382 to PSPTO_3380).

Transcriptome analysis of Pseudomonas syringae identifies new genes, noncoding RNAs, and antisense activity.

To fully understand how bacteria respond to their environment, it is essential to assess genome-wide transcriptional activity. New high-throughput sequencing technologies make it possible to query the transcriptome of an organism in an efficient unbiased manner. We applied a strand-specific method to sequence bacterial transcripts using Illumina's high-throughput sequencing technology. The resulting sequences were used to construct genome-wide transcriptional profiles. Novel bioinformatics analyses were developed and used in combination with proteomics data for the qualitative classification of transcriptional activity in defined regions. As expected, most transcriptional activity was consistent with predictions from the genome annotation. Importantly, we identified and confirmed transcriptional activity in areas of the genome inconsistent with the annotation and in unannotated regions. Further analyses revealed potential RpoN-dependent promoter sequences upstream of several noncoding RNAs (ncRNAs), suggesting a role for these ncRNAs in RpoN-dependent phenotypes. We were also able to validate a number of transcriptional start sites, many of which were consistent with predicted promoter motifs. Overall, our approach provides an efficient way to survey global transcriptional activity in bacteria and enables rapid discovery of specific areas in the genome that merit further investigation.

Construction of an rsmX co-variance model and identification of five rsmX non-coding RNAs in Pseudomonas syringae pv. tomato DC3000.

Non-coding RNAs (ncRNAs) are important components of many regulatory pathways and have key roles in regulating diverse functions. In the Pseudomonads, the two-component system, GacA/S, directly regulates at least two well-characterized ncRNAs, RsmZ and RsmY, which act by sequestration of translation repressor proteins to control expression of various exoproducts. Pseudomonas fluorescens CHA0 possesses a third ncRNA, RsmX, which also participates in this regulatory pathway. In this study we confirmed expression of five rsmX ncRNAs in Pseudomonas syringae pv. tomato DC3000, and determined the distribution of the members of the rsmX ncRNA family by screening available genomic sequences of the Pseudomonads. Variable numbers of the rsmX family exist in Pseudomonas genomes, with up to five paralogs in Pseudomonas syringae strains. In Pseudomonas syringae pv. tomato DC3000, the rsmX genes are 112 to 120 nucleotides in size and are predicted by structural analysis to contain multiple exposed GGA motifs, which is consistent with structural features of the Rsm ncRNAs. We also found that these rsmX ncRNA genes share a conserved upstream region suggesting that their expression is dependent upon the global response regulator, GacA.

Complete Genome Sequence of the Necrotrophic Plant-Pathogenic Bacterium Pectobacterium brasiliense 1692.

We report here the complete and annotated genome sequence of Pectobacterium brasiliense 1692, a Gram-negative enterobacterium that can cause soft rot disease in many plant hosts.

Complete Genome Sequence of a Gram-Positive Bacterium, Leifsonia sp. Strain PS1209, a Potato Endophyte.

We report the complete and annotated genome sequence of a Gram-positive bacterium, Leifsonia sp. strain PS1209, a potato endophyte that was isolated from apparently healthy tubers of potato cultivar NY166. The circular genome is 4,091,164 bp long, with a GC content of 69.08%, containing 3,926 genes.

Deep RNA sequencing of L. monocytogenes reveals overlapping and extensive stationary phase and sigma B-dependent transcriptomes, including multiple highly transcribed noncoding RNAs.

BACKGROUND: Identification of specific genes and gene expression patterns important for bacterial survival, transmission and pathogenesis is critically needed to enable development of more effective pathogen control strategies. The stationary phase stress response transcriptome, including many sigmaB-dependent genes, was defined for the human bacterial pathogen Listeria monocytogenes using RNA sequencing (RNA-Seq) with the Illumina Genome Analyzer. Specifically, bacterial transcriptomes were compared between stationary phase cells of L. monocytogenes 10403S and an otherwise isogenic DeltasigB mutant, which does not express the alternative sigma factor sigmaB, a major regulator of genes contributing to stress response, including stresses encountered upon entry into stationary phase. RESULTS: Overall, 83% of all L. monocytogenes genes were transcribed in stationary phase cells; 42% of currently annotated L. monocytogenes genes showed medium to high transcript levels under these conditions. A total of 96 genes had significantly higher transcript levels in 10403S than in DeltasigB, indicating sigmaB-dependent transcription of these genes. RNA-Seq analyses indicate that a total of 67 noncoding RNA molecules (ncRNAs) are transcribed in stationary phase L. monocytogenes, including 7 previously unrecognized putative ncRNAs. Application of a dynamically trained Hidden Markov Model, in combination with RNA-Seq data, identified 65 putative sigmaB promoters upstream of 82 of the 96 sigmaB-dependent genes and upstream of the one sigmaB-dependent ncRNA. The RNA-Seq data also enabled annotation of putative operons as well as visualization of 5'- and 3'-UTR regions. CONCLUSIONS: The results from these studies provide powerful evidence that RNA-Seq data combined with appropriate bioinformatics tools allow quantitative characterization of prokaryotic transcriptomes, thus providing exciting new strategies for exploring transcriptional regulatory networks in bacteria.