ORIGINAL RESEARCH article

Front. Microbiol., 11 October 2017

Sec. Food Microbiology

Volume 8 - 2017 | https://doi.org/10.3389/fmicb.2017.01910

Home Alone: Elimination of All but One Alternative Sigma Factor in Listeria monocytogenes Allows Prediction of New Roles for σB

  • Department of Food Science, Cornell University, Ithaca, NY, United States

Article metrics

View details

35

Citations

4,5k

Views

1,4k

Downloads

Abstract

Among Listeria monocytogenes' four alternative σ factors, σB controls the largest regulon. As σB-dependent transcription of some genes may be masked by overlaps among regulons, and as some σB-dependent genes are expressed only under very specific conditions, we hypothesized that the σB regulon is not yet fully defined. To further extend our understanding of the σB regulon, we used RNA-seq to identify σB-dependent genes in an L. monocytogenes strain that expresses σB following rhamnose induction, and in which genes encoding the other alternative sigma factors have been deleted. Analysis of RNA-seq data with multiple bioinformatics approaches, including a sliding window method that detects differentially transcribed 5′ untranslated regions (UTRs), identified 105 σB-dependent transcription units (TUs) comprising 201 genes preceded by σB-dependent promoters. Of these 105 TUs, 7 TUs comprising 15 genes had not been identified previously as σB-dependent. An additional 23 genes not reported previously as σB-dependent were identified in 9 previously recognized σB-dependent TUs. Overall, 38 of these 201 genes had not been identified previously as members of the L. monocytogenes σB regulon. These newly identified σB-dependent genes encode proteins annotated as being involved in transcriptional regulation, oxidative and osmotic stress response, and in metabolism of energy, carbon and nucleotides. In total, 18 putative σB-dependent promoters were newly identified. Interestingly, a number of genes previously identified as σB-dependent did not show significant evidence for σB-dependent transcription in our experiments. Based on promoter analyses, a number of these genes showed evidence for co-regulation by σB and other transcriptional factors, suggesting that some σB-dependent genes require additional transcriptional regulators along with σB for transcription. Over-expression of a single alternative sigma factor in the absence of all other alternative sigma factors allowed us to: (i) identify new σB-dependent functions in L. monocytogenes, such as regulation of genes involved in 1,2-propanediol utilization (LMRG_00594-LMRG_00611) and biosynthesis of pyrimidine nucleotides (LMRG_00978-LMRG_00985); and (ii) identify new σB-dependent genes involved in stress response and pathogenesis functions. These data further support that σB not only regulates stress response functions, but also plays a broad role in L. monocytogenes homeostasis and resilience.

Introduction

Listeria monocytogenes is a Gram-positive foodborne pathogen that causes the serious invasive disease listeriosis, predominantly in susceptible populations such as the immunocompromised, pregnant women, and adults over 65 years old (Goulet et al., 2012). In the US, L. monocytogenes causes around 1,600 human listeriosis cases resulting in ~260 deaths annually (Scallan et al., 2011). L. monocytogenes' ability to rapidly respond to changing environmental conditions enables it to survive under a wide range of circumstances such as those that may be encountered during food processing as well as in animal or human hosts (Gray et al., 2006). Alternative sigma factors represent one key regulatory mechanism that allows bacteria to adjust rapidly to different environments. Differential association between alternative sigma factors and core RNA polymerase allows the RNA polymerase to recognize specific promoter sequences and initiate transcription of targeted genes under specific conditions. L. monocytogenes has up to four alternative sigma factors (σB, σC, σH, and σL) in addition to the housekeeping sigma factor σA. The four alternative sigma factors regulate transcription of genes important for virulence and for response to various stress and growth conditions (Chaturongakul et al., 2011). To date, the general stress response regulator σB is the most extensively studied alternative sigma factor in L. monocytogenes; σB has been shown to control a regulon of more than 180 genes (Raengpradub et al., 2008). Specifically, σB plays important roles in virulence and stress response, including transition to stationary phase and resistance to acid, osmotic, arsenate, oxidative, and cold stresses (O'byrne and Karatzas, 2008; Mujahid et al., 2013b).

Regulons controlled by other L. monocytogenes alternative sigma factors are less well defined. σC, which is present only in lineage II strains of L. monocytogenes, appears to have a small regulon (<10 genes) that contributes to thermal resistance (Zhang et al., 2005). σL, which is a member of the σ54 family, has a regulon of >70 genes involved in carbon and amino acid metabolism and cold stress (Raimann et al., 2009). σH appears to have the second largest regulon among the alternative sigma factors (>150 genes) (Chaturongakul et al., 2011) and has been reported to contribute to growth in minimal media, alkaline stress response, virulence (Rea et al., 2004), and competence (Liu et al., 2016; Medrano Romero and Morikawa, 2016). Previous studies in L. monocytogenes have identified considerable overlaps among the regulons controlled by these four alternative sigma factors (Chaturongakul et al., 2011; Mujahid et al., 2013a). Redundant regulation of a gene by multiple sigma factors may make it difficult to classify a gene as a member of a specific regulon following deletion of a single alternative sigma factor. Hence the commonly used approach of comparing transcript levels between a wild type (WT) strain and a corresponding isogenic mutant bearing a deletion of a targeted alternative sigma factor gene may not provide a complete picture of the regulon for a given alternative sigma factor as the approach may fail to identify all genes that are co-regulated by multiple transcriptional regulators. Furthermore, many of the previous transcriptional studies of the σB regulon (Chatterjee et al., 2006; Ollinger et al., 2008, 2009; Raengpradub et al., 2008; Toledo-Arana et al., 2009; Oliver et al., 2010; Chaturongakul et al., 2011; Ribeiro et al., 2014) have used microarray technology, which is heavily affected by extrinsic noise and is limited to annotated open reading frames (ORFs). In order to provide a more complete definition of the L. monocytogenes σB regulon, we performed RNA sequencing (RNA-seq) using an L. monocytogenes strain that expresses σB under an inducible promoter as the only active alternative sigma factor, thereby removing potential redundancies in regulation by σC, σH, and σL.

Materials and methods

Bacterial strains, mutant construction

The quadruple alternative σ factor mutant (ΔBCHL; FSL C3-0135) of L. monocytogenes strain 10403S (Mujahid et al., 2013a), which was used here as the background strain, was modified further to overexpress sigB from a rhamnose inducible promoter. The strain was constructed using methods analogous to those previously reported for construction of a strain overexpressing σH (Liu et al., 2016). Briefly, the sigB ORF was amplified from L. monocytogenes 10403S by PCR and cloned downstream of the rhamnose inducible promoter Prha in plasmid pLF1 (Fieseler et al., 2012); this construct was transformed into E. coli strain SM10 λpir to allow for subsequent conjugation of the plasmid into L. monocytogenes 10403S ΔBCHL, followed by chromosomal integration of the Prha-sigB construct in an arginine tRNA gene locus (yielding strain 10403S::ΔBCHL Prha-sigB; FSL B2-0425) (Lauer et al., 2002). A control strain (ΔBCHL-Prha; FSL B2-0429) was constructed by introducing plasmid pLF1 into the chromosome of L. monocytogenes 10403S ΔBCHL through conjugation and chromosomal integration (see Supplementary Table 1 for strains, plasmids and primers used).

Growth conditions and rhamnose induction

Strains were streaked from frozen Brain Heart Infusion (BHI) glycerol stocks onto a BHI agar plate, followed by incubation at 37°C for 24 h. A single colony was subsequently inoculated into 5 ml of BHI broth in a 10 ml tube and incubated at 37°C for 18 h with shaking (220 rpm). After incubation, 50 μl BHI culture was inoculated into fresh 5 ml BHI broth in a 10 ml tube and grown to OD600 0.4–0.5 at 37°C with shaking. Induction of sigB transcription was performed by adding 250 μl of 1 M rhamnose stock solution to 5 ml OD600 0.4–0.5 bacterial cultures (for a final concentration of 50 mM rhamnose), followed by incubation at 37°C for an additional 30 min. Induction with rhamnose was performed for both 10403S::ΔBCHL Prha-sigB and ΔBCHL-Prha. qRT-PCR using the SYBR Green Master Mix Reagent (Life Technologies) and the ABI Prism 7000 Sequence Detection System (Applied Biosystems, Foster City, CA) was used to determine 50 mM rhamnose as the optimal rhamnose concentration for sigB induction (Supplementary Figure 1). For these experiments, transcript levels were determined for sigB and the housekeeping gene rpoB in strain 10403S::ΔBCHL Prha-sigB. Expression level differences were determined by the ΔΔCt method using the housekeeping gene rpoB as reference gene (Livak and Schmittgen, 2001).

RNA isolation and sequencing

To define the σB regulon, strains 10403S::ΔBCHL Prha-sigB and ΔBCHL-Prha were grown to log phase, followed by induction with 50 mM rhamnose for 30 min, as previously described (Liu et al., 2016). RNA isolation and cDNA library construction were also performed as previously described (Liu et al., 2016). All experiments were performed in three biological replicates. Indexed RNA-seq libraries were quantified by digital PCR and sequencing was carried out on a Hiseq 2500 (single-end, 150-bp per read) at the Cornell Core Facility for RNA-sequencing. RNA-seq data are accessible through GEO Series accession number GSE94284.

RNA-seq alignment, coverage and differential expression analysis

RNA-seq data analysis was performed as previously described (Liu et al., 2016). Briefly, sequence reads were aligned to a 10403S genome using the BWA-mem algorithm (Li and Durbin, 2010) and the data for coverage per base on sense and antisense strands were obtained separately using samtools (Li et al., 2009). Differential expression of genes between the two strains (ΔBCHL::Prha and ΔBCHL::Prha-sigB) was initially assessed using the Bayseq package for R version 2.2.0 (Hardcastle and Kelly, 2010). Genes were considered differentially expressed (DE) if the FDR (False Discovery Rate) was <0.05 and the Fold Change (FC) was more than 2.0 or <0.5. Genes with significantly higher transcript levels in the ΔBCHL-Prha-sigB strain were identified as upregulated by σB (FC > 2.0) while genes with significantly higher transcript levels in ΔBCHL-Prha were identified as downregulated (FC < 0.5). A sliding window method described previously (Liu et al., 2016) was also applied to identify significant differential expression of fragments along the chromosome; this method provides more sensitive identification of σB-dependent genes, particularly if genes are preceded by multiple promoters (as described in Raimann et al., 2009 and Liu et al., 2016). Briefly, the 10403S genome was divided into windows of 51 nt (window size) with 25 nt overlap (sliding window) and the RNA-Seq coverage was obtained for each window. Bayseq was then used, as described above, to identify windows with significant differential expression (FDR < 0.05 and FC > 2.0) between the ΔBCHL::Prha-sigB and the ΔBCHL::Prha strains. Overlapping windows were considered as one fragment, and fragments mapped to the same transcription unit (TU) were considered to comprise one TU. TUs are defined as all of the genes and/or noncoding RNAs (ncRNAs) that are included in a transcript initiated from a specific promoter. These analyses used an L. monocytogenes 10430S genome annotation where TUs were designated based on an RNA-seq analysis of the 10403S WT strain (Oliver et al., 2009); a file with this genome annotation in Genbank format is available from the authors upon request.

Comparison between previously identified σB-dependent genes and genes identified in this study

To allow us to define newly identified σB-dependent genes, we compiled data on the L. monocytogenes σB regulon from 12 independent studies that used genome-wide transcriptomics and proteomics approaches to compare different L. monocytogenes WT strains and corresponding isogenic sigB deletion mutants (Wemekamp-Kamphuis et al., 2004; Chatterjee et al., 2006; Abram et al., 2008a,b; Ollinger et al., 2008, 2009; Orsi et al., 2008; Raengpradub et al., 2008; Oliver et al., 2009, 2010; Toledo-Arana et al., 2009; Chaturongakul et al., 2011; Ribeiro et al., 2014). Combined, these studies identified a total of 902 different genes as σB-dependent; only genes that were identified here, but not among these previously reported 902 genes were classified as “newly identified σB-dependent genes”. We also defined a “high confidence” core σB regulon for L. monocytogenes lineage II (which includes strain 10403S, used here), which contained 184 genes that have been reported as σB-dependent in at least three independent experiments among the 12 studies detailed above. This “high confidence” core σB regulon was used to define the previously identified σB-dependent genes that were not found to be differentially expressed (NDE) in our study reported here. The datasets supporting the conclusions of this article are available in the GEO repository (Private link for reviewers: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=gbktssowflaften&acc=GSE94284).

Promoter search and consensus

The transcription start sites (TSS) and 5′ untranslated regions (UTRs) of the differentially expressed TUs were manually annotated based on visual inspection of the locations where RNA-seq coverage data abruptly increased. Subsequently, the 5-50 nt regions upstream of the TSSs were manually scanned for σB-dependent promoter consensus sequences (GTTT-N12−16-GGGTAT). A consensus sequence logo of the newly identified σB-dependent promoters was generated using the WebLogo generator (Crooks et al., 2004).

5′ rapid amplification of cDNA ends (5′ RACE) analysis

Six newly identified σB-dependent promoters were selected to map promoter regions with the 5′ RACE system (ThermoFisher); these experiments were performed in three biological replicates. Briefly, isolated RNA was reverse transcribed into cDNA with gene-specific primers and cDNA was tailed with dCTP by terminal transferase. The products were then amplified with a nested gene-specific primer and an abridged anchored primer in a touchdown PCR with DreamTaq PCR master mix (ThermoFisher). PCR products were visually analyzed with agarose gel electrophoresis.

Go term enrichment analysis

Gene Ontology (GO) enrichment analysis to identify GO terms enriched among all 201 genes identified in the 105 σB-dependent TUs was performed using the GOseq 1.24.0 package for R (Young et al., 2010) as previously described (Tang et al., 2015). All GO terms containing more than five genes were analyzed for enrichment. FDR correction for multiple testing was applied and only GO terms with FDR < 0.05 were considered significant.

Comparative genomics of newly identified σB-dependent promoters

Newly identified σB-dependent promoters were searched against a database containing 27 finished genomes representing L. monocytogenes lineage I (n = 10), L. monocytogenes lineage II (n = 9), L. monocytogenes lineage III (n = 4), Listeria innocua (n = 1), Listeria seeligeri (n = 1), Listeria ivanovii (n = 1), and Listeria welshimeri (n = 1) using BLAST [32]. Matches with coverage >70% and identity >60% were considered significant and the promoters were, therefore, considered to be present in the respective genome.

Results

Combined bioinformatics analyses identify a number of new σB-dependent genes that showed significant differential transcript levels in the respective ORF

Using RNA-seq data from an L. monocytogenes strain with deletions in genes encoding all four alternative sigma factors (ΔBCHL) as well as from a corresponding isogenic strain where sigB was re-introduced under rhamnose regulation (ΔBCHL-Prha-sigB), we employed two separate bioinformatics approaches to identify σB-dependent genes. These bioinformatics approaches included (i) Bayseq analysis that assessed differential expression using RNA-seq coverage for the full length of all annotated genes (“the traditional Bayseq approach”) and (ii) Bayseq analysis that assessed differential expression using RNA-seq coverage for sliding windows that covered the complete L. monocytogenes chromosome (“the sliding window approach”).

With the traditional Bayseq approach, we initially identified 141 genes, including 7 ncRNAs, that showed significantly higher transcripts levels (FDR < 0.05; FC > 2.0) in the L. monocytogenes strain overexpressing σB(ΔBCHL-Prha-sigB), as compared to the ΔBCHL strain that does not express any alternative sigma factors (Table 1). Among the 141 genes with higher transcript levels, 117 were preceded by upstream σB-dependent promoters. When the sliding window approach was applied, we identified 299 fragments that showed significantly higher transcripts levels (FDR < 0.05; FC > 2.0) in the L. monocytogenes strain overexpressing σB. These fragments represented 177 genes, including 7 ncRNAs; 124 of these genes (including all 7 ncRNAs) were preceded by σB-dependent promoters (Supplementary Table 2). Together, these two analyses identified 193 genes as upregulated by σB; 133 of these genes (including 7 noncoding RNAs) were preceded by σB-dependent promoters (Table 1). These 133 genes are located in 93 TUs that each includes a σB-dependent promoter. Overall, 141 of the 193 upregulated genes identified here had been reported as σB-dependent in previous microarray, RNA-seq, or proteomics studies in L. monocytogenes (Wemekamp-Kamphuis et al., 2004; Chatterjee et al., 2006; Abram et al., 2008a,b; Ollinger et al., 2008, 2009; Raengpradub et al., 2008; Oliver et al., 2009, 2010; Toledo-Arana et al., 2009; Chaturongakul et al., 2011; Ribeiro et al., 2014). Among the 52 genes newly identified, with these approaches, as upregulated by σB, 16 genes (representing 10 TUs) were preceded by a putative σB-dependent promoter, while 36 did not include an upstream σB-dependent promoter. In the subsequent sections, we will only focus on those genes that showed both evidence for higher transcript levels in the presence of σB and have an upstream σB-dependent promoter.

Table 1

CategoryaBayseqSliding windowCombinedbNew
(1) No. of genes identified as σB-dependent (FC > 2.0; FDR < 0.05)14117719352
    (2) No. of genes in (1) that have upstream σB promoters11712413316
    (3) σB-dependent TUs represented by genes in (2)86c87c93c4
(4) σB-dependent TUs that were identified by only 5UTR DE fragment and contain upstream σBpromotersNAd12123
    (5) Genes found in the TUs in (4)NA16164
(6) Total number of TUs with both σBpromoters and DE fragment in ORF or 5UTReNANA105 TUsc7 TUs (represent 15 genes)
    (7) Total gene number of σB regulon members located in the TUs in (6)NANA201 genes38 genes (represent 16 TUs)

Summary of traditional Bayseq and sliding window analyses results for upregulated genes and TUs.

a

Indented categories indicates groups that represent a subset of the main, non-indented category in bold.

b

The category “combined” represents the number of genes or TUs that were identified by either the BaySeq or the sliding window approach or both.

c

LMRG_02094 has two σB-dependent promoters directly upstream of the gene; this gene was counted as two TUs identified by BaySeq and as two TUs identified by the sliding window analysis. This gene was also counted as two TUs in the “combined” column.

d

NA means “Not Applicable”.

e

Category 6 represents the combination of categories (3) and (4).

With the traditional Bayseq approach, we also identified 18 genes that showed significantly lower transcripts levels (FDR < 0.05; FC < 0.5; see Supplementary Table 3) in the L. monocytogenes strain overexpressing σB(ΔBCHL-Prha-sigB), as compared to the ΔBCHL strain. Among these 18 downregulated genes, 9 had been reported as downregulated in previous studies (Wemekamp-Kamphuis et al., 2004; Chatterjee et al., 2006; Abram et al., 2008a,b; Ollinger et al., 2008, 2009; Raengpradub et al., 2008; Oliver et al., 2009, 2010; Toledo-Arana et al., 2009; Chaturongakul et al., 2011; Ribeiro et al., 2014). Most of the nine genes newly identified as downregulated by σB encode transport proteins, including (i) LMRG_01249 and LMRG_01248, which encode subunits of a PTS galactitol transporter; (ii) LMRG_01581, which encodes ArpJ, an amino acid ABC transporter permease; and (iii) LMRG_01595, which encodes an MFS transporter permease. Other downregulated genes include LMRG_01332, which encodes Listeria adhesion protein (Lap), an alcohol acetaldehyde dehydrogenase involved in pathogenesis (Jagadeesan et al., 2010); LMRG_00198, encoding a phosphoglycerate mutase; LMRG_01596, encoding a shikimate 5-dehydrogenase, and LMRG_01597 encoding an NADH oxidase. Interestingly, one of the nine newly identified σB-downregulated genes was LMRG_01250, which encodes a putative transcriptional regulator.

Sliding window analyses of 5′UTR regions identified 12 additional TUs with σB-dependent transcription and upstream σB-dependent promoters in the absence of significant differential transcript levels in the respective ORF

The sliding window approach used here not only identified differentially regulated genes via identification of differentially transcribed fragments that mapped within an ORF (described above), but also identified differentially transcribed fragments that mapped into different 5′UTR, even if the associated downstream genes did not show evidence for differential transcription. This approach thus identified additional TUs (which may include one or multiple genes) as σB-dependent. For example, in the intergenic region upstream of srtA we identified a fragment with an FC of 3.77 located upstream of the srtA σA-dependent promoter (Figure 1). Visual inspection of the region upstream of this fragment identified a σB-dependent promoter. Hence, using the sliding window approach, we identified srtA as a new σB-dependent TU. Overall, this approach identified 12 additional TUs (Supplementary Table 4) with σB-dependent transcription and upstream σB-dependent promoters. For each of these TUs, the associated genes located in the TU did not show significant evidence for higher transcript levels in the strain overexpressing σB.

Figure 1

Figure 1

Gene co-transcribed from both σB-dependent and a σA-dependent promoters. RNA-seq coverage is shown for the ΔBCHL Prha-sigB strain (red) and ΔBCHL-Prha strain (green). Transcription start sites (TSS) are indicated by arrows labeled with the regulating sigma factor. The dotted line below the schematic represents the region found to be differentially expressed (FDR < 0.05) with the sliding window approach. The highest window fold change (FC) is shown under the dotted line.

Overall analyses reveal a total of 105 σB-dependent TUs that cover 201 genes, including a number of genes not previously identified as part of the σB regulon

Overall, our analyses identified 105 σB-dependent TUs (i.e., TUs upregulated in the presence of σB) that are preceded by σB-dependent promoters; this includes (i) 93 TUs that showed significant differential transcript levels in one or multiple ORFs that are part of a TU and (ii) 12 TUs that showed significant differential transcript levels in a fragment located in the 5′UTR. Seven of these 105 TUs had not been previously identified as σB-dependent (meaning none of the 15 genes in these seven TUs had been identified as σB-dependent in previous studies (Wemekamp-Kamphuis et al., 2004; Chatterjee et al., 2006; Abram et al., 2008a,b; Ollinger et al., 2008, 2009; Raengpradub et al., 2008; Oliver et al., 2009, 2010; Toledo-Arana et al., 2009; Chaturongakul et al., 2011; Ribeiro et al., 2014). Overall, the 105 σB-dependent TUs include 201 different genes; 38 of these 201 genes had not been previously identified as members of the L. monocytogenes σB regulon (Tables 1, 2). While additional genes were newly identified as σB-dependent through our bioinformatics approaches (see Table 1), not all of these genes met the dual criteria of having both an upstream σB-dependent promoter and differential upregulation in the ΔBCHL-Prha-sigB strain.

Table 2

New σB regulon memberLocus tag for strain EGDeTUGene product
OXIDATIVE STRESS RESPONSE
LMRG_00891lmo1439LMRG_00891Manganese superoxide dismutase
LMRG_02644lmo0222LMRG_02644-LMRG_02643Heat shock protein 33
LMRG_02443lmo0014qoxABCDA quinol oxidase
OSMOTIC STRESS RESPONSE
LMRG_02775lmo1701LMRG_02775-LMRG_02778An hypothetical protein
LMRG_02115 LMRG_02116lmo1015 lmo1016LMRG_02114-LMRG_02116Permease protein OpuAB, glycine betaine-binding protein Op of glycine betaine ABC transport system,
PYRIMIDINE BIOSYNTHESIS
LMRG_00978-LMRG_00985lmo1831- lmo1838LMRG_00978-LMRG_00985Enzymes in uridine-5′-phosphate biosynthesis
1,2-PROPANEDIOL UTILIZATION
LMRG_00596-00599, 00601,00603-00606, 00608-00611lmo1153-1156,1160-1163, 1165-1168LMRG_00594-LMRG_00611Enzymes that degrade carbon compound 1,2-propanediol
METABOLISM OF CARBON
LMRG_01960 LMRG_01962lmo2736 lmo2734LMRG_01960-LMRG_01963A glycerate kinase and a glycosyl hydrolase
TRANSCRIPTIONAL REGULATION
LMRG_01913lmo2784LMRG_01913Transcriptional antiterminator of lichenan operon, BglG family
LMRG_00050 LMRG_00051lmo0359 lmo0360LMRG_00050-LMRG_00051A fructose-biphosphate hydrolase and a DeoR family putative transcriptional regulator
NO CLEARLY ANNOTATED FUNCTIONS
LMRG_01659lmo2173LMRG_01658-LMRG_01659σL -dependent activator
LMRG_01784lmo2464LMRG_01784-LMRG_01786Transcriptional regulator of tetR Family
LMRG_01676lmo2156LMRG_01676A hypothetic protein
LMRG_01199_asN/ALMRG_01199_asAntisense RNA
LMRG_00529lmo1067LMRG_00529GTP-binding protein TypA/BipA
LMRG_02699lmo2569LMRG_02695-LMRG_02699Periplasmic oligopeptide-binding protein oppA

Newly identified σB regulon members grouped by function.

Newly identified σB regulon members include genes involved in stress response, metabolism of carbon and nucleotides, and virulence

The 38 genes newly identified as members of the σB regulon (Table 2) encode proteins involved in oxidative and osmotic stress response, metabolism of energy, carbon and nucleotides, transcriptional regulation, and other functions of L. monocytogenes. Newly identified σB-dependent genes involved in oxidative stress response include LMRG_00891, LMRG_02644, and LMRG_02443, which is the second gene in the qoxABCD operon. Our data newly identified σB-dependent promoters upstream of all of these operons. We also newly identified LMRG_02775 as a member of the σB regulon. While the other members of the LMRG_02775-LMRG_02778 operon had previously been identified as σB-dependent (Oliver et al., 2010; Chaturongakul et al., 2011), we newly identified a σB-dependent promoter upstream of this operon, suggesting direct σB-dependent transcription of this 4-gene operon, which encodes four proteins, including one OsmC/Ohr family protein. While the functions of the genes in this operon have not yet been elucidated, OsmC proteins are induced by ethanol and osmotic stresses and Ohr proteins are induced by organic peroxide and are involved in organic hydroperoxide detoxification in Gram-negative bacteria (Atichartpongkul et al., 2001). In other intracellular pathogens such as Mycobacterium tuberculosis, OsmC proteins have been reported to protect against organic hydroperoxide stress (Saikolappan et al., 2011), suggesting this operon may contribute to oxidative or osmotic stress response in L. monocytogenes. We also newly identified two members of the LMRG_02114-LMRG_02116 operon as part of the σB regulon and newly identified a σB-dependent promoter upstream of this operon. This operon encodes the OpuAA, OpuAB, and OpuAC subunits of the glycine betaine ABC transport system. Expression of the opuA operon is under osmotic control in B. subtilis, and a B. subtilis mutant strain lacking the OpuA transport system showed a considerably decreased ability to uptake the osmoprotectant glycine betaine (Kempf and Bremer, 1995).

Our study also identified a number of new members of the σB regulon involved in energy, carbon and nucleotide metabolism. Specifically, we newly identified LMRG_00978-LMRG_00985 as a σB-dependent TU; genes in this TU encode a set of enzymes involved in uridine-5′-phosphate biosynthesis, a pyrimidine ribonucleotides de novo biosynthesis pathway (Figure 2). Importantly, our data also indicate that all genes in the LMRG_00594-LMRG_00611 operon are part of the σB regulon; this operon encodes all enzymes involved in utilization of carbon compound 1,2-propanediol, which is an important carbon source during infection for bacterial pathogens such as S. enterica (Bobik et al., 1999) and was reported to be degraded by L. innocua (Xue et al., 2008). In addition to providing direct RNA-seq evidence for significant differential transcript levels for 6 genes in this operon, which had not been previously identified as σB-dependent, we also identified a σB-dependent promoter upstream of this operon; in previous studies, only 5 genes in this 18-gene operon had been identified as σB-dependent (Wemekamp-Kamphuis et al., 2004; Chatterjee et al., 2006; Abram et al., 2008a,b; Ollinger et al., 2008, 2009; Raengpradub et al., 2008; Oliver et al., 2009, 2010; Toledo-Arana et al., 2009; Chaturongakul et al., 2011; Ribeiro et al., 2014). Our data also newly identified a σB-dependent promoter upstream of the three-gene operon LMRG_01960-LMRG_01962 with LMRG_01960 and LMRG_01962 newly identified as σB-dependent. We also newly identified a σB-dependent promoter upstream of LMRG_01432-LMRG_01431, suggesting direct σB regulation of this operon.

Figure 2

Figure 2

Pathway of uridine-5′-phosphate biosynthesis. Genes encoding the enzymes involved in this pathway are labeled in purple. All genes in the TU LMRG_00985-LMRG_00978 except LMRG_00981 are shown in this pathway.

We also identified new members of the σB regulon that encode transcriptional regulators, further supporting that σB is part of a complex regulatory network in L. monocytogenes. Specifically, we newly identified LMRG_01913, which encodes a lichenan operon transcriptional antiterminator as a σB-dependent TU. In addition, we newly identified LMRG_00050-LMRG_00051 as a σB-dependent TU, in which LMRG_00051 encodes a putative transcriptional regulator. LMRG_01659, encoding a σ54-dependent activator (Arous et al., 2004) and LMRG_01784, encoding a TetR family transcriptional regulator, were also newly identified as members of the σB regulon. While the other genes in TUs LMRG_01658-LMRG_01659 and LMRG_01784-LMRG_01786 had previously been reported as σB-dependent (Raengpradub et al., 2008; Oliver et al., 2009, 2010; Ollinger et al., 2009; Toledo-Arana et al., 2009; Chaturongakul et al., 2011; Ribeiro et al., 2014), we newly identified σB-dependent promoters upstream of both operons.

Our data also identified a σB-dependent promoter upstream of a TU including a single gene encoding sortase A, srtA. SrtA is involved in proteolysis and virulence in L. monocytogenes (Garandeau et al., 2002). Additional newly identified σB-dependent TUs (where none of the genes in a TU had previously been identified as σB-dependent) included LMRG_00529 and two TUs (LMRG_01676 and LMRG_01199_as) with no clearly annotated functions.

Gene ontology analyses identify a role for σB in regulating genes encoding functions involved in transport, homeostasis, pathogenesis, and nucleotide metabolism

Gene ontology (GO) enrichment analysis performed for the σB regulon (i.e., all 201 genes that were part of the 105 σB-dependent TUs identified here) found 18 GO terms to be overrepresented (see Table 3), including GO terms related to transport, homeostasis, metabolic and catabolic processes. In addition, GO terms related to “pyrimidine nucleotide biosynthetic and metabolic processes”; “glycolysis, metabolic and catabolic processes of propanediol, diol and alcohol,” were also over-represented, suggesting a role for σB in regulating general metabolic functions.

Table 3

GO IDDE genes in this GO term (LMRG_)GO term
GO:000606600050, 00236,00272, 00594, 00600-00607, 00611, 01432, 01627, 01789-01793, 01962, 02219, 02632, 02695Alcohol metabolic process
GO:000609100050, 00981, 01627, 01790-01793, 02442-02444, 02632Generation of precursor metabolites and energy
GO:000609600050, 01627, 01790-01793, 02632Glycolytic process
GO:000622000978-00985, 01076Pyrimidine nucleotide metabolic process
GO:000622100978-00985, 01076Pyrimidine nucleotide biosynthetic process
GO:000905600050, 00221, 00600-00603, 00605, 00607, 00611, 01030, 01627, 01790-01793, 02219, 02304, 02472, 02632Catabolic process
GO:001541800873, 00874, 00877-00880, 02114-02116Quaternary-ammonium-compound-transporting ATPase activity
GO:001569500873, 00874, 02114-02116Organic cation transport
GO:001569700873, 00874, 02114-02116Quaternary ammonium group transport
GO:003010400873, 00874, 02114-02116Water homeostasis
GO:003431100594, 00600-00603, 00605, 00607, 00611Diol metabolic process
GO:003431300600-00603, 00605, 00607, 00611Diol catabolic process
GO:004424800050, 00221, 00600-00603, 00605, 00607, 00611, 01030, 01627, 01790-01793, 02219, 02304, 02472, 02632Cellular catabolic process
GO:004427500050, 00221, 01030, 01627, 01790-01793, 02219, 02632Cellular carbohydrate catabolic process
GO:004616400050, 00600-00607, 00611, 01627, 01790-01793, 02219, 02632Alcohol catabolic process
GO:004887800873, 00874, 02041, 02114-02116Chemical homeostasis
GO:005114300594, 00600-00603, 00605, 00607, 00611Propanediol metabolic process
GO:005114400600-00603, 00605, 00607, 00611Propanediol catabolic process

Enriched GO terms for 201 genes in 105 σB-dependent TUsa.

a

σB regulon members are defined as genes in the σB-dependent TUs that are preceded by σB-dependent promoters and showed significant differential transcript levels in one or more fragments located within the TU.

Conservation of newly identified σB-dependent promoters in L. monocytogenes and Listeria spp.

Overall, we newly identified 18 putative σB-dependent promoters here (see Figure 3 for a consensus sequence). 5′ RACE performed on six of these promoters confirmed five of these promoters as σB-dependent (Figure 4). For pduA, where 5′RACE did not detect a clear σB-dependent transcript, a putative σA-dependent promoter was found 14 nt upstream of the putative σB-dependent promoter; transcripts from this σA-dependent promoter may have masked the σB-dependent transcript. Across all 18 newly identified putative σB-dependent promoters, the first three nucleotides of the −10 signal sequence (GGG) are generally conserved among both these newly identified putative σB-dependent promoters and the previously published σB-dependent promoter sequences (Oliver et al., 2009), the last three nucleotides appear to be less conserved among the new promoters identified here, as compared to the previously reported σB-dependent promoter sequences (Figure 3). Through comparative genomics, we assessed that all but one of the 18 σB-dependent promoters newly identified in L. monocytogenes 10403S were present in the 23 additional L. monocytogenes genomes analyzed here; the σB-dependent promoter upstream of the anti-sense RNA (asRNA), which overlaps LMRG_00359 in the opposite strand, was not found in the four lineage III genomes analyzed (Table 4). Overall, −35 and −10 promoter regions were highly conserved among different strains and lineages, even though some promoters show lineage- or strain-specific sequence features (Table 4). Further analysis showed that only 10 of the 18 σB-dependent promoters newly identified in L. monocytogenes were found in all 5 Listeria species analyzed (i.e., L. monocytogenes, L. innocua, L. ivanovii, L. welshimeri, and L. seeligeri) (Table 4). While these comparisons suggest that we identified new σB-dependent functions that are largely conserved in L. monocytogenes and to a lesser extent in other Listeria sensu strictu species, future comparative genomics analyses that utilize additional closed genomes as they become available will help to further define lineage-, strain-, and species-specific members of the SigB regulon.

Figure 3

Figure 3

Sequence logo for newly identified putative σB-dependent promoters. (A) Sequence logo for the 18 σB-dependent promoters associated with newly identified σB-dependent genes. (B) Sequence logo for σB-dependent promoters published in previous study (Oliver et al., 2009).

Figure 4

Figure 4

5′RACE PCR confirmation of six putative σB-dependent promoters newly identified here by RNA-seq. The six putative σB-dependent promoters selected for confirmation are located upstream of qoxA, srtA, pduA, LMRG_01431, LMRG_02778, and LMRG_01913. LMRG_01948 was used as a positive control, as it was reported to have a strong σB-dependent promoter (Abram et al., 2008b). The image shows the 5′RACE PCR products in a 3% agarose gel; the two lanes for each target represent reaction performed with RNA isolated from the ΔBCHL Prha-sigB strain, which expresses σB (+) and the ΔBCHL-Prha strain, which does not express σB (−). The σB-dependent transcript bands are shown in rectangles. Multiple 5′RACE bands were expected for these genes as they were all initially characterized by having a putative σB-dependent promoter as well as an additional σA-dependent promoter. Five of six promoters selected for confirmation displayed clear σB-dependent transcript bands, supporting the existence of σB-dependent promoters. 5′RACE PCR for pduA did not yield a band that confirmed the putative σB-dependent promoter that was identified by RNA-seq; a putative σA-dependent promoter was found 14 nt upstream of the putative σB-dependent promoter; transcripts from this σA-dependent promoter may have masked the σB-dependent transcript. The expected sizes of the σB-dependent transcript bands are around 250 bp for qoxA, 170 bp for srtA, 165 bp for pduA, 275 bp for LMRG_01431, 185 bp for LMRG_02778, 275 bp for LMRG_01913, and 210 bp for LMRG_01948. The image was adjusted for contrast, intensity levels and saturation using Photoshop; manipulations were performed on the whole picture and no specific bands were enhanced or modified. The results shown here are representative of three biological replicates.

Table 4

PromotersL. monocytogenesL. innocuaL. welshimeriL. ivanoviiL. seeligeri
Lineage ILineage IILineage III
LMRG_00986
−35 signalGCACGCACGCACGCACGCACACACGCAC
−10 signalGGCTCTGGCTCTGGCTCTGGCTCTGGCTCTGGCTCGGGCTCT
LMRG_02643
−35 signalGAATGAATGAATGAATGAATGAATGAAT
−10 signalGGGCTTGGGCTTGGGCTTGGGCTTGGGCTTGGGTTAGGGTTA
gbuA
−35 signalGCTTGCTTaGCTTaGCTTGCTTGCTTGCTT
−10 signalGGGCCGGGGCCGGGGCCGGGGCCGGGGCCGGGGCCGGGGCCG
LMRG_00529
−35 signalCGCTGGCTCGCTCGCTCGCTCGCCCGCT
−10 signalGGTAAGGGTAAGGGTAAGGGTAAGGGTAAGGGTAAGGGTAAG
pduA
−35 signalGTAAGTAAGTAAGTAAGTAAGTAAGTAA
−10 signalAGGAGGAGGAGGAGGAGGAGGAGGAGGAGGAGGAGGAGGAGG
LMRG_02813
−35 signalGTACGTACGTACGTACNFbNFNF
−10 signalGGGTATGGGTATGGGTATGGGTATNFNFNF
LMRG_01676
−35 signalGATTGATTAATTGATATATTGATTGATT
−10 signalGGGATAGGGATAGGGATAAGGATAAGGAAAAGGATAAGGATA
qoxA
−35 signalGTTTcGTTTGTTTGTTTGTTTGTTTGTTT
−10 signalGGGAATGGGAATGGGAATGGGAATGGGAATGGGAATGGGAAT
LMRG_00050
−35 signalGATAGATAGATAGATAGATANFNF
−10 signalTGGTTCTGGTTCTGGTTCTGGTTCTGGTCCNFNF
LMRG_00271
−35 signalGTTTGTTTGTTTGTTTGTTTNFNF
−10 signalGGGGAAGGGGAAGGGGAAGGGGAAGGGGAANFNF
asRNA(LMRG_00359)
−35 signalGTTTGTTTNFGTTTNFNFNF
−10 signalGGGCTAGGGCTANFGAGCTANFNFNF
srtA
−35 signalGATTGATTGATTGATTGATTGATTGATT
−10 signalGGGAAAGGGAAAGGGAAAGGGAAAGGGAAAGGGAAAGGGAAA
LMRG_02094
−35 signalGGTTGGTTGGTTGGTTGGTTGATTGGTT
−10 signalGGGTAAGGGTAAGGGTAAGGGTAAGGGTAAGGGTAAGGGAAA
LMRG_01431
−35 signalGTTAGTTAGTTAGTTAGTTANFNF
−10 signalGGGTAAGGGTAAGGGTAAGGGTAAGGGTAANFNF
LMRG_02778
−35 signalGCTTGCTTGCTTGCTTGCTTGCTTGCTT
−10 signalGGGTATGGGTATGGGTATGGGTATGAGTATGCGTATGCGTAT
asRNA (LMRG_01199)
−35 signalGATTGATTGATTGTTTNFATTTNF
−10 signalTGGAATTGGAATTGGAATTGGAGTNFTGGAATNF
LMRG_01963
−35 signalGTTTGTTTGTTTNFNFNFNF
−10 signalAGCAATGGCAATGGCAATNFNFNFNF
LMRG_01913
−35 signalGTTTGTTTGTTTGTTTGTTTGTTTNF
−10 signalGG(G/A)CTAdGGGCTA(T/G)AGCTTeGGGCGAGGGCGAGAATTANF

σB-dependent −35 and −10 promoter regions for the newly identified σB-dependent operons.

a

With the exception of strains 10403S and SLCC5850 (both lineage II strains), all the lineage II and III strains present a deletion of a thymine between the −35 and −10 regions.

b

NF, Not found.

c

Strain SLCC2482 (lineage I) presents an insertion of a thymine between the −35 and −10 region.

d

Strain ATCC has the sequence GGACTA at the −10 region while all other Lineage I strains have the sequence GGGCTA.

e

Strain SLCC2376 has the sequence GAGCTT at the −10 region while all other lineage III strains have the sequence TAGCTT.

Comparison with previously reported σB-dependent regulon members identifies a number of σB-dependent genes that are also regulated by σA

In order to compare the σB-dependent genes identified here to previously reported σB-dependent genes, we used data from previous microarray, RNA-seq and proteomic studies in different strains and under different stress conditions (Wemekamp-Kamphuis et al., 2004; Chatterjee et al., 2006; Abram et al., 2008a,b; Ollinger et al., 2008, 2009; Raengpradub et al., 2008; Oliver et al., 2009, 2010; Toledo-Arana et al., 2009; Chaturongakul et al., 2011; Ribeiro et al., 2014), to define a core σB regulon for L. monocytogenes lineage II, which includes the strain we used here (10403S). Overall, 184 genes have been reported as σB-dependent at least three times in prior publications (representing high confidence for each as a member of the σB regulon). Our data showed that 110 of the 193 genes that were identified as significantly upregulated by σB in our study also were included in this “high confidence regulon” of 184 genes. Therefore, overexpression of σB in ΔBCHL failed to identify 74 genes previously reported as upregulated by σB (Supplementary Table 5).

Promoter searches upstream of the 5′ UTR of these 74 genes identified upstream σB-dependent promoters for 35 of these genes; these 35 genes represent 28 TUs with corresponding σB-dependent promoters. Twenty-three of these 28 TUs had both σA- and σB-dependent promoters (Table 5). Therefore, it is possible that expression from the σA-dependent promoter masked expression from the σB-dependent promoter in these genes, not allowing their identification as σB-dependent in our study. While no clear σA-dependent promoters were identified for the other five TUs, RNA-seq coverage patterns in the quadruple mutant, which only expresses σA, suggested that four of these TUs may be co-regulated by σA. These four TUs were LMRG_00261-LMRG_00262, LMRG_00335-LMRG_00338, LMRG_01794, and htrA (LMRG_02579) whose σB-dependent transcription generates a long (473 nt) 5'UTR, which may have a regulatory function.

Table 5

TUGenesOther genes in the same operonPromoter identified (from upstream to downstream)Genes or fragmentsa identified as differentially expressed in this studyComments
1LMRG_00100, LMRG_00101σB, σA, σANo
2LMRG_00262LMRG_00261σBLMRG_00261σA-dependent expression pattern
3LMRG_00334σB, σANo
4LMRG_00336LMRG_00335, LMRG_00337, LMRG_00338,σBLMRG_00335, LMRG_00337,σA-dependent expression pattern
5LMRG_00530σA, σBFrag 43108, Frag 43109Two DE fragments identified between the promoters
6LMRG_00672LMRG_00671σA, σBNo
7LMRG_00745σA, σBNoHost factor-I protein
8LMRG_00884LMRG_00885σB, σANo
9LMRG_00906LMRG_00910, LMRG_00909, LMRG_00908, LMRG_00907σA, σH, σBNoRpoD (σA)
10LMRG_01076, LMRG_01077LMRG_01078, LMRG_01079, LMRG_01080σA, σBLMRG_01078, LMRG_01079, LMRG_01080
11LMRG_01140σA, σBNoAlpha-acetotactate decarboxylase
12LMRG_01284σB, σAFrag 86921 to 86928DE fragments identified in 5'UTR
13LMRG_01361LMRG_01360σA, σA, σBLMRG_01360
14LMRG_01432LMRG_01431σB, σALMRG_01431Glycerol kinase
15LMRG_01737σB, σAFragments 101803 to 101809DE fragments identified in 5'UTR
16LMRG_01794σBNoσA-dependent transcription initiating from LMRG_01793 (eno) is running through LMRG_01794 despite the presence of a Rho-independent transcription terminator in between
17LMRG_02000, LMRG_02001, LMRG_02002σA, σBNo
18LMRG_02028σA, σBsrtA, DE fragments identified in 5'UTR
19LMRG_02055, LMRG_02056LMRG_02057σB, σAfrag 38333nagA and nagB operon. Two NagR attenuatormotives were identified suggesting co-regulation with NagR
20LMRG_02215, LMRG_02216LMRG_02217σA, σBLMRG_02217
21LMRG_02219LMRG_02218σB, σALMRG_02218rpiB
22LMRG_02317, LMRG_02320LMRG_02311, LMRG_02312, LMRG_02313, LMRG_02314, LMRG_02315, LMRG_02316, LMRG_02318, LMRG_02319σA, σBrsbV and rsbX from sigB operon
23LMRG_02579LMRG_02580σBNohtrA, σA-dependent expression pattern; σB-dependent transcription generates a long (473 nt) 5′UTR, which may have regulatory function
24LMRG_02624σB, σANohly, has a PrfA box, a σA-dependent promoter, and a putative σB-dependent promoter suggested in previous studies
25LMRG_02632σB, σAFrags 8396 to 8402ldh
26LMRG_02731σA, σBNo
27LMRG_02772σA, σBFrags 68894 to 68895

Co-regulation of σA and σB among the NDE σB-dependent genes.

a

Each fragment is composed of one or more 51 nt overlapping windows.

Some σB-dependent TUs are also regulated by other transcriptional regulators

Three TUs previously reported as σB-dependent, but not identified as such in our study reported here, were found to be co-regulated by σB-dependent promoters and other alternative sigma factors and transcriptional factors. LMRG_02294, encoding a hypothetical protein, is co-transcribed with LMRG_02293, a gene encoding a putative amidohydrolase. A σB-dependent−10 region and a σL-dependent promoter in close proximity have been identified for this TU, suggesting co-regulation of transcription (Supplementary Table 6). As the strains used in this study lack an active σL, overexpression of σB itself might not be sufficient to significantly upregulate the TU. LMRG_02294 only showed a non-significant upregulation in the strain expressing σB.

LMRG_00906 (rpoD, encoding σA), which is part of the LMRG_00911-LMRG_00906 TU, is co-transcribed from an upstream σB-dependent promoter, a σH-dependent promoter upstream of LMRG_00908 and a σA-dependent promoter upstream of LMRG_00911. Strong expression from the σA-dependent promoter as well as the absence of σH may explain why no significant differential expression of rpoD was observed in our study. LMRG_02624 (hly) has a PrfA box, a σA-dependent promoter, and a putative σB-dependent promoter suggested in previous studies (Orsi et al., 2008; Tsai et al., 2011). hly was only reported as σB-dependent with active expression of PrfA (Ollinger et al., 2009; Toledo-Arana et al., 2009; Ribeiro et al., 2014), which was not shown to be significantly differentially expressed in our study.

Discussion

A number of studies have analyzed L. monocytogenes alternative-sigma-factor-dependent expression profiles at the transcriptional and protein levels based on comparisons between WT and single deletion mutant strains. However, regulon redundancy and overlaps in regulation by multiple sigma factors is common, as has been shown in Listeria (Chaturongakul et al., 2011), B. subtilis (Mascher et al., 2007), and in other bacteria (Nuss et al., 2010). We thus used an experimental design that enabled us to study the functions and genes regulated by σB in the absence of all other alternative sigma factors in order to further explore the L. monocytogenes σB regulon and σB-dependent gene regulation. Over-expression of a single alternative sigma factor in the absence of all other alternative sigma factors, the specific strategy used here, has previously been shown to eliminate redundant regulation by multiple alternative sigma factors, thus allowing improved insights into both the function of a single sigma factor and coregulation by different sigma factors (Mujahid et al., 2013a; Wang et al., 2014; Liu et al., 2016). As detailed above, this approach allowed us to (i) identify new σB-dependent functions in L. monocytogenes, such as regulation of genes involved in 1,2-propanediol utilization and biosynthesis of pyrimidine nucleotides; (ii) identify new σB-dependent genes involved in stress response and virulence functions; and (iii) further define a role for σB in L. monocytogenes homeostasis.

L. monocytogenes σB directly upregulates >100 TUs and >200 genes

Our approach shows that the σB regulon includes >100 TUs and >200 genes that are directly up-regulated by σB. As RNA-seq has been well documented, in both L. monocytogenes (Oliver et al., 2009) and other organisms (Vivancos et al., 2010), to provide quantitative data that are well correlated with qPCR data, qPCR confirmation of genes newly identified as σB-dependent was not deemed necessary; rather σB-dependent transcriptional start sites were confirmed for selected genes via 5′ RACE. By comparison, previous studies have reported the σB regulon in 10403S as including >110 genes and >80 promoters (Raengpradub et al., 2008; Abram et al., 2008a,b; Oliver et al., 2009; Ollinger et al., 2009; Chaturongakul et al., 2011). Combined results from previous microarray, RNA-seq and proteomic studies in different L. monocytogenes strains under different stress conditions show a large σB pan-regulon of up to 902 genes (Wemekamp-Kamphuis et al., 2004; Chatterjee et al., 2006; Abram et al., 2008a,b; Ollinger et al., 2008, 2009; Raengpradub et al., 2008; Oliver et al., 2009, 2010; Toledo-Arana et al., 2009; Chaturongakul et al., 2011; Ribeiro et al., 2014) and a core regulon in lineage II strains of 184 genes, as supported by three or more previous experiments. Overall, our data further confirm that σB is the alternative sigma factor that regulates the largest regulon; other regulators that directly regulate large gene sets in L. monocytogenes include σH, which directly regulates more than 60 genes (Chaturongakul et al., 2011; Liu et al., 2016).

Newly identified σB-dependent promoters reveal a novel role for σB in oxidative stress, 1,2-propanediol utilization and biosynthesis of pyrimidine nucleotides

While L. monocytogenes σB clearly has been shown to be important for survival under many stress conditions, such as exposure to heat, bile, osmotic and acid stress (Ferreira et al., 2001, 2003; Sue et al., 2004; Seifart Gomes et al., 2011), the role of σB in oxidative stress response is less clear. Previous studies have reported results ranging from hypersensitivity (Ferreira et al., 2001; Oliver et al., 2010) to hyperresistance (Moorhead and Dykes, 2003; Boura et al., 2016) to oxidative stress for σB null mutants. For example, a previous study reported that viability of a 10403S ΔsigB mutant was 100-fold lower than the WT under oxidative stress induced by cumene hydroperoxide in stationary phase (Ferreira et al., 2001). However, a recent study using L. monocytogenes 10403S and EGD-e showed that deletion of sigB led to hyperresistance to H2O2 in stationary phase cells grown under aerobic conditions (Boura et al., 2016). Further, while studies in B. cereus showed that loss of σB or σB regulon members showed impaired glucose-starvation-induced resistance to H2O2 (Engelmann and Hecker, 1996; Antelmann et al., 1997; Zuber, 2009), the ΔsigB mutant showed either a similar degree of resistance (Engelmann and Hecker, 1996) or hyperresistance to H2O2 induced oxidative stress (Van Schaik et al., 2005). Interestingly, we found both hslO and the qoxABCD operon to be transcribed from a σB-dependent promoter, suggesting a direct role of σB in resistance to oxidative stress. Among heat shock proteins, HslO has a unique chaperone activity because it is redox regulated and protects both thermally and oxidatively damaged proteins from irreversible aggregation (Kim et al., 2001). The importance of qoxABCD in oxidative stress response is also supported by observations in B. subtilis where a qoxABCD mutant shows reduced growth under aerobic conditions (Winstedt and Von Wachenfeldt, 2000). Upregulation of qox genes in L. monocytogenes in the presence of glycerol further indicates their involvement in adaption to an aerobic environment (Joseph et al., 2008). Our results suggest that the operon qoxABCD should be a focus for future studies on the molecular mechanisms behind oxidative stress resistance. As the exact role of σB in oxidative stress response in L. monocytogenes has not been elucidated to date, identification of these genes as regulated by σB-dependent promoters may provide fresh insight into the role of σB in oxidative stress response.

Interestingly, our data also found the propanediol utilization (pdu) operon to be transcribed from an upstream σB-dependent promoter, revealing a new role for σB in the utilization of 1,2-propanediol. 1,2-propanediol is produced during the catabolism of rhamnose and fucose, two sugars that are abundant in mammalian glycoconjugates (Staib and Fuchs, 2014). In addition, 1,2-propanediol is also used in various food products as a direct food additive (Cameron et al., 1998). The metabolism of 1,2-propanediol has been well studied in Salmonella; this organism is not only able to use 1,2-propanediol as a carbon source, but the pdu operon has also been shown to be important for virulence in Salmonella (Bobik et al., 1999). So far, studies in L. monocytogenes showed that transcription of some of the pdu genes increases under certain conditions, including lack of glucose (Nilsson et al., 2005), during intracellular growth (Joseph et al., 2006), and during growth on the surface of cold smoked salmon slices (Tang et al., 2015). In L. innocua, the pduD gene is required for 1,2-propanediol metabolism and 17 genes within the pduA-to-pduF gene cluster are induced when cells are grown in a medium containing 1,2-propanediol (Xue et al., 2008). The studies detailed above suggest the importance of 1,2-propanediol utilization in Listeria during growth in foods and in the host environment; our study reveals evidence for a regulatory role for σB in 1,2-propanediol utilization. This finding further illustrates important regulatory roles for σB in adaptation and transition to food and host environments.

Another σB-dependent promoter newly identified here is upstream of the LMRG_00978 to LMRG_00986 operon, which encodes proteins involved in pyrimidine ribonucleotide biosynthesis. Previous studies showed that genes involved in the de novo synthesis of purines (purA, purQ, lmo1771) and pyrimidines (pyrE) are required for intracellular proliferation of L. monocytogenes, suggesting that these bases and nucleotides are not provided by the host cell, but must be synthesized by the bacterium (Schauer et al., 2010). Our study supports σB as one of the key regulators for nucleotide metabolism, which is essential for L. monocytogenes to succeed in the intracellular environment.

Newly identified σB-dependent promoters reveal additional mechanisms for σB-dependent roles in osmotic stress, glycerol utilization, and virulence

In addition to identification of new σB-dependent functions, as detailed above, our data also newly identified σB-dependent genes involved in functions previously linked to σB regulation, including oxidative stress, osmotic stress, glycerol utilization, and virulence. The response to osmotic stress is one of the first reported specific functions for σB in L. monocytogenes (Becker et al., 1998). Specific genes involved in osmotic stress response, such as opuCA, have a confirmed σB-dependent promoter (Kazmierczak et al., 2003) and are activated by osmotic stress in a σB-dependent manner (Sue et al., 2004). The osmotic activation of σB in L. monocytogenes is rapid, transient, and proportional to the magnitude of the osmotic stress applied (Utratna et al., 2011). Here, two operons with a role in osmotic stress response were newly identified as being directly regulated by σB, including (i) LMRG_02114-LMRG_02116 and (ii) LMRG_02775-LMRG_02778, which includes one gene encoding an OsmC/Ohr family protein. Previous studies have found that one of two genes encoding OsmC/Ohr family proteins is a member of the σB regulon in B. subtilis (Kempf and Bremer, 1995; Volker et al., 1998), suggesting a general role for σB-dependent regulation of OsmC/Ohr family proteins. These osmotic stress response proteins, newly identified to be σB-dependent, broaden our knowledge of how σB fine-tunes gene regulation to support chemical and water homeostasis in L. monocytogenes. We hypothesize that these σB-dependent genes could play an important role under high osmotic stress conditions such as those encountered by L. monocytogenes during gastrointestinal passage (Sue et al., 2004).

Glycerol is used as an alternative carbon source by intracellular L. monocytogenes (Eylert et al., 2008). Previous publications also suggested the importance of the ability to metabolize a series of carbon sources for successful infection by intracellular human pathogens (Olive and Sassetti, 2016) such as M. tuberculosis (de Carvalho et al., 2010). In L. monocytogenes, the mutant defective in the uptake and metabolism of glycerol showed impaired intracellular growth (Eylert et al., 2008). In a previous proteomics study, σB was found to regulate the LMRG_02000-LMRG_02002 operon, which encodes subunits of dihydroxyacetone kinase, an enzyme involved in glycerol metabolism in L. monocytogenes (Abram et al., 2008b). In addition, a sigB mutant had a diminished ability to use glycerol as a sole carbon source (Abram et al., 2008b). In our study, the LMRG_01432-LMRG_01431 operon, which encodes a glycerol kinase and a glycerol uptake facilitator protein, was found to have a σB-dependent promoter, further supporting an important role for σB in regulation of glycerol metabolism. As growth in glycerol also upregulates PrfA activity (Joseph et al., 2008; Stoll et al., 2008), this newly revealed role for σB in glycerol metabolism further supports a link between carbohydrate metabolism and pathogenesis in L. monocytogenes.

While σB previously has been recognized as a transcriptional regulator of a number of virulence genes, including prfA, inlA, inlC2, and inlD (Mcgann et al., 2007; Guldimann et al., 2016), our study revealed that another virulence gene, srtA, which is involved in proteolysis and processing of internalin proteins (Garandeau et al., 2002), is also directly regulated by σB. SrtA is required for the cell wall anchoring of InlA and, presumably, for the anchoring of other LPXTG-containing proteins that are involved in listerial infections (Bierne et al., 2002). A recent study also reported that inhibition of sortase A by chalcone could prevent L. monocytogenes infection (Li et al., 2016). In S. aureus, a number of sortase substrate proteins were previously observed to have higher or lower expression in a the ΔsigB mutant (Hempel et al., 2010). Our finding confirms and extends the role σB plays in pathogenesis, especially in the interaction between pathogen and host surface molecules during host cell invasion.

Co-regulation of genes by σB and σA as well as other transcriptional regulators supports importance of σB in resilience and homeostasis

While L. monocytogenes σB has previously been well defined as a regulator of stress response and virulence functions, our data reported here further support a broader role for σB in resilience and homeostasis. The enriched GO terms related to general transport and homeostasis confirmed contributions of σB in coordinating complex networks responsive to changing environmental conditions.

Importantly, we also found a number of σB-dependent TUs and genes preceded by both σB and σA-dependent promoters (Supplementary Table 7). Key genes in this category include prfA, rsbV, qoxABCD (aerobic respiration), mogR (repressor of genes involved in flagella), phoU (phosphate transport regulation), cggR (central glycolytic genes regulator), ltrC (response to cold) and others. Genes preceded by both a promoter regulated by a housekeeping sigma factor and a promoter regulated by an alternative sigma factor have been identified in other bacteria. For example, in B. subtilis, clpC is preceded by two promoters (σA and σB) (Kruger et al., 1996) and expression of the ureABC operon, which encodes urease, is dependent on σA and σH (Wray et al., 1997). Additionally, the Extracytoplasmic Function (ECF) sigma factors of B. subtilis are known to co-regulate promoter regions with an existing σA promoter (Eiamphungporn and Helmann, 2008; Kingston et al., 2011). In E. coli, extensive overlaps between promoters of the primary sigma factor σ70 and alternative sigma factors such as σ32 and others also suggest coregulation of gene expression by multiple sigma factors under various growth conditions (Wade et al., 2006).

Seventy-four genes identified in previous studies as σB-dependent were not differentially expressed as σB-dependent in this study. A considerable number of these genes were co-regulated by more than one sigma factor, predominantly by σB and σA. A likely explanation regarding why some previously identified σB-dependent genes showed no differential expression in this study may reflect the growth conditions under which the experiments were conducted: the cells were in logarithmic growth with no additional imposed stress conditions. Previous studies used stress conditions such as stationary phase and salt stress to induce σB expression. In contrast, our experiment used artificial induction of σB expression by addition of rhamnose to log-phase cells. The transcriptional activity of the σA housekeeping sigma factor has been found to be diminished under stress conditions (Sharma and Chatterji, 2010; Delumeau et al., 2011). Fast-growing cells depend on σA to regulate expression of housekeeping genes such as those involved in protein biosynthesis, DNA replication and structural proteins. Therefore, under log-phase, σA is usually highly active to support the reproductive needs of the cells (Delumeau et al., 2011). Under the conditions used in this study, genes co-transcribed by both σA and σB under log-phase might have very little expression originating from the sometimes weaker σB promoter if the stronger σA promoter is located upstream or very close to the σB promoter.

Conclusions

In this study, we used RNA-seq to explore the role of σB in L. monocytogenes by overexpressing σB in a strain where genes encoding all other alternative sigma factors had been deleted. Combined with prior data revealing important roles for σB in pathogenesis and stress response, identification of new putative σB-dependent promoters upstream of a number of genes indicates a broader regulon for this alternative sigma factor, which also appears to contribute to cellular homeostasis.

Transcriptomic approaches such as RNA-seq and microarrays as well as proteomic approaches are powerful tools to explore the regulation of sigma factors, however, they have limited abilities to distinguish indirectly regulated genes from directly regulated genes. Further experiments with approaches such as ChIP-seq may allow a better definition of the direct regulons of sigma factors.

Statements

Author contributions

YL performed the RNA-seq experiments and initial data analyses and was a major contributor in writing the manuscript; RO performed sliding window analysis of RNA-Seq data and promoter identification; KB, MW, and VG co-wrote the manuscript and conceived the study. All authors read and approved the final manuscript.

Funding

This project was supported by a grant from the National Institute of Allergy and Infectious Diseases (NIAID) of the National Institutes of Health (NIH) (2 RO1 AI052151-05A1).

Acknowledgments

We thank Dr. Martin J. Loessner at ETH Zurich for the kind gift of plasmid pLF1.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2017.01910/full#supplementary-material

Supplementary Figure 1

qRT-PCR results of rhamnose induction of sigB and the σB-dependent gene LMRG_01602. Expression of σB under rhamnose induction was tested by qRT-PCR. LMRG_01602 was used as the targeted σB-dependent gene and rpoB as the reference gene; error bars show the standard deviation for the 2 biological replicates; Y axis represents the fold change of the gene expression level. Fold change is calculated as −ΔCt = CtrpoB − Cttarget gene. The σB-dependent gene referred in second graph is LMRG_01602.

Supplementary Table 1

Plasmid, primers and strains used in this study.

Supplementary Table 2

Genes and TUs upregulated by σB identified by traditional BaySeq and sliding window analyses.

Supplementary Table 3

Genes downregulated by σB identified by traditional BaySeq.

Supplementary Table 4

σB-dependent TUs identified by sliding window analyses, but not by the traditional BaySeq approach.

Supplementary Table 5

List of NDE σB-dependent genes.

Supplementary Table 6

5′UTRs of the NDE σB-dependent genes.

Supplementary Table 7

TUs with σB-dependent promoters.

    Abbreviations

  • asRNA

    anti-sense RNA

  • BHI

    Brain Heart Infusion

  • BWA

    Burrows-Wheeler Aligner

  • DE

    differentially expressed

  • ECF

    Extracytoplasmic Function

  • FC

    fold change

  • FDR

    false discovery rate

  • GEO

    Gene Expression Omnibus

  • GO

    Gene Ontology

  • ncRNA

    noncoding RNA

  • NDE

    not differentially expressed

  • NEB

    New England Biolabs

  • NIAID

    National Institute of Allergy and Infectious Diseases

  • NIH

    National Institutes of Health

  • NRC

    normalized RNA-seq coverage

  • ORF

    open reading frame

  • 5′ RACE

    5′ Rapid amplification of cDNA ends

  • RNA-seq

    RNA sequencing

  • TU

    transcription unit

  • UTR

    untranslated region

  • WT

    wild type

  • TSS

    transcription start sites.

References

  • 1

    AbramF.StarrE.KaratzasK. A.Matlawska-WasowskaK.BoydA.WiedmannM.et al. (2008a). Identification of components of the Sigma B regulon in Listeria monocytogenes that contribute to acid and salt tolerance. Appl. Environ. Microbiol.74, 68486858. 10.1128/AEM.00442-08

  • 2

    AbramF.Wan-LinS.WiedmannM.BoorK. J.CooteP.BottingC.et al. (2008b). Proteomic analyses of a Listeria monocytogenes mutant lacking σB identify new components of the σB regulon and highlight a role for σB in the utilization of glycerol. Appl. Environ. Microbiol.74, 594604. 10.1128/AEM.01921-07

  • 3

    AntelmannH.EngelmannS.SchmidR.SorokinA.LapidusA.HeckerM. (1997). Expression of a stress- and starvation-induced dps/pexB-homologous gene is controlled by the alternative sigma factor SigmaB in Bacillus subtilis. J. Bacteriol.179, 72517256. 10.1128/jb.179.23.7251-7256.1997

  • 4

    ArousS.BuchrieserC.FolioP.GlaserP.NamaneA.HebraudM.et al. (2004). Global analysis of gene expression in an rpoN mutant of Listeria monocytogenes. Microbiology150, 15811590. 10.1099/mic.0.26860-0

  • 5

    AtichartpongkulS.LoprasertS.VattanaviboonP.WhangsukW.HelmannJ. D.MongkolsukS. (2001). Bacterial Ohr and OsmC paralogues define two protein families with distinct functions and patterns of expression. Microbiology147, 17751782. 10.1099/00221287-147-7-1775

  • 6

    BeckerL. A.CetinM. S.HutkinsR. W.BensonA. K. (1998). Identification of the gene encoding the alternative sigma factor SigmaB from Listeria monocytogenes and its role in osmotolerance. J. Bacteriol.180, 45474554.

  • 7

    BierneH.MazmanianS. K.TrostM.PucciarelliM. G.LiuG.DehouxP.et al. (2002). Inactivation of the srtA gene in Listeria monocytogenes inhibits anchoring of surface proteins and affects virulence. Mol. Microbiol.43, 869881. 10.1046/j.1365-2958.2002.02798.x

  • 8

    BobikT. A.HavemannG. D.BuschR. J.WilliamsD. S.AldrichH. C. (1999). The propanediol utilization (pdu) operon of Salmonella enterica serovar Typhimurium LT2 includes genes necessary for formation of polyhedral organelles involved in coenzyme B12-dependent 1, 2-propanediol degradation. J. Bacteriol.181, 59675975.

  • 9

    BouraM.KeatingC.RoyetK.PaudyalR.O'donoghueB.O'byrneC. P.et al. (2016). The presence of SigB in Listeria monocytogenes strains EGD-e and 10403S leads to hypersensitivity to hydrogen peroxide in stationary phase under aerobic conditions. Appl. Environ. Microbiol.82, 45844591. 10.1128/AEM.00709-16

  • 10

    CameronD. C.AltarasN. E.HoffmanM. L.ShawA. J. (1998). Metabolic engineering of propanediol pathways. Biotechnol. Prog.14, 116125. 10.1021/bp9701325

  • 11

    ChatterjeeS. S.HossainH.OttenS.KuenneC.KuchminaK.MachataS.et al. (2006). Intracellular gene expression profile of Listeria monocytogenes. Infect. Immun.74, 13231338. 10.1128/IAI.74.2.1323-1338.2006

  • 12

    ChaturongakulS.RaengpradubS.PalmerM. E.BergholzT. M.OrsiR. H.HuY.et al. (2011). Transcriptomic and phenotypic analyses identify coregulated, overlapping regulons among PrfA, CtsR, HrcA, and the alternative sigma factors σB, σC, σH, and σL in Listeria monocytogenes. Appl. Environ. Microbiol.77, 187200. 10.1128/AEM.00952-10

  • 13

    CrooksG. E.HonG.ChandoniaJ. M.BrennerS. E. (2004). WebLogo: a sequence logo generator. Genome Res.14, 11881190. 10.1101/gr.849004

  • 14

    de CarvalhoL. P.FischerS. M.MarreroJ.NathanC.EhrtS.RheeK. Y. (2010). Metabolomics of Mycobacterium tuberculosis reveals compartmentalized co-catabolism of carbon substrates. Chem. Biol.17, 11221131. 10.1016/j.chembiol.2010.08.009

  • 15

    DelumeauO.LecointeF.MuntelJ.GuillotA.GuedonE.MonnetV.et al. (2011). The dynamic protein partnership of RNA polymerase in Bacillus subtilis. Proteomics11, 29923001. 10.1002/pmic.201000790

  • 16

    EiamphungpornW.HelmannJ. D. (2008). The Bacillus subtilis σM regulon and its contribution to cell envelope stress responses. Mol. Microbiol.67, 830848. 10.1111/j.1365-2958.2007.06090.x

  • 17

    EngelmannS.HeckerM. (1996). Impaired oxidative stress resistance of Bacillus subtilis sigB mutants and the role of katA and katE. FEMS Microbiol. Lett.145, 6369. 10.1111/j.1574-6968.1996.tb08557.x

  • 18

    EylertE.ScharJ.MertinsS.StollR.BacherA.GoebelW.et al. (2008). Carbon metabolism of Listeria monocytogenes growing inside macrophages. Mol. Microbiol.69, 10081017. 10.1111/j.1365-2958.2008.06337.x

  • 19

    FerreiraA.O'byrneC. P.BoorK. J. (2001). Role of σB in heat, ethanol, acid, and oxidative stress resistance and during carbon starvation in Listeria monocytogenes. Appl. Environ. Microbiol.67, 44544457. 10.1128/AEM.67.10.4454-4457.2001

  • 20

    FerreiraA.SueD.O'byrneC. P.BoorK. J. (2003). Role of Listeria monocytogenes σB in survival of lethal acidic conditions and in the acquired acid tolerance response. Appl. Environ. Microbiol.69, 26922698. 10.1128/AEM.69.5.2692-2698.2003

  • 21

    FieselerL.SchmitterS.TeiserskasJ.LoessnerM. J. (2012). Rhamnose-inducible gene expression in Listeria monocytogenes. PLoS ONE7:e43444. 10.1371/journal.pone.0043444

  • 22

    GarandeauC.Reglier-PoupetH.DubailI.BerettiJ. L.BercheP.CharbitA. (2002). The sortase SrtA of Listeria monocytogenes is involved in processing of internalin and in virulence. Infect. Immun.70, 13821390. 10.1128/IAI.70.3.1382-1390.2002

  • 23

    GouletV.HebertM.HedbergC.LaurentE.VaillantV.De ValkH.et al. (2012). Incidence of listeriosis and related mortality among groups at risk of acquiring listeriosis. Clin. Infect. Dis.54, 652660. 10.1093/cid/cir902

  • 24

    GrayM. J.FreitagN. E.BoorK. J. (2006). How the bacterial pathogen Listeria monocytogenes mediates the switch from environmental Dr. Jekyll to pathogenic Mr. Hyde. Infect. Immun.74, 25052512. 10.1128/IAI.74.5.2505-2512.2006

  • 25

    GuldimannC.BoorK. J.WiedmannM.Guariglia-OropezaV. (2016). Resilience in the face of uncertainty: sigma B fine-tunes gene expression to support homeostasis in Gram-positive bacteria. Appl. Environ. Microbiol.82, 44564469. 10.1128/AEM.00714-16

  • 26

    HardcastleT. J.KellyK. A. (2010). BaySeq: empirical Bayesian methods for identifying differential expression in sequence count data. BMC Bioinformatics11:422. 10.1186/1471-2105-11-422

  • 27

    HempelK.Pane-FarreJ.OttoA.SieversS.HeckerM.BecherD. (2010). Quantitative cell surface proteome profiling for SigB-dependent protein expression in the human pathogen Staphylococcus aureus via biotinylation approach. J. Proteome Res.9, 15791590. 10.1021/pr901143a

  • 28

    JagadeesanB.KooO. K.KimK. P.BurkholderK. M.MishraK. K.AroonnualA.et al. (2010). LAP, an alcohol acetaldehyde dehydrogenase enzyme in Listeria, promotes bacterial adhesion to enterocyte-like Caco-2 cells only in pathogenic species. Microbiology156, 27822795. 10.1099/mic.0.036509-0

  • 29

    JosephB.MertinsS.StollR.ScharJ.UmeshaK. R.LuoQ.et al. (2008). Glycerol metabolism and PrfA activity in Listeria monocytogenes. J. Bacteriol.190, 54125430. 10.1128/JB.00259-08

  • 30

    JosephB.PrzybillaK.StuhlerC.SchauerK.SlaghuisJ.FuchsT. M.et al. (2006). Identification of Listeria monocytogenes genes contributing to intracellular replication by expression profiling and mutant screening. J. Bacteriol.188, 556568. 10.1128/JB.188.2.556-568.2006

  • 31

    KazmierczakM. J.MithoeS. C.BoorK. J.WiedmannM. (2003). Listeria monocytogenes σB regulates stress response and virulence functions. J. Bacteriol.185, 57225734. 10.1128/JB.185.19.5722-5734.2003

  • 32

    KempfB.BremerE. (1995). OpuA, an osmotically regulated binding protein-dependent transport system for the osmoprotectant glycine betaine in Bacillus subtilis. J. Biol. Chem.270, 1670116713. 10.1074/jbc.270.28.16701

  • 33

    KimS. J.JeongD. G.ChiS. W.LeeJ. S.RyuS. E. (2001). Crystal structure of proteolytic fragments of the redox-sensitive Hsp33 with constitutive chaperone activity. Nat. Struct. Biol.8, 459466. 10.1038/87639

  • 34

    KingstonA. W.SubramanianC.RockC. O.HelmannJ. D. (2011). A σW-dependent stress response in Bacillus subtilis that reduces membrane fluidity. Mol. Microbiol.81, 6979. 10.1111/j.1365-2958.2011.07679.x

  • 35

    KrugerE.MsadekT.HeckerM. (1996). Alternate promoters direct stress-induced transcription of the Bacillus subtilis clpC operon. Mol. Microbiol.20, 713723. 10.1111/j.1365-2958.1996.tb02511.x

  • 36

    LauerP.ChowM. Y. N.LoessnerM. J.PortnoyD. A.CalendarR. (2002). Construction, characterization, and use of two Listeria monocytogenes site-specific phage integration vectors. J. Bacteriol.184, 41774186. 10.1128/JB.184.15.4177-4186.2002

  • 37

    LiH.ChenY.ZhangB.NiuX.SongM.LuoZ.et al. (2016). Inhibition of sortase A by chalcone prevents Listeria monocytogenes infection. Biochem. Pharmacol.106, 1929. 10.1016/j.bcp.2016.01.018

  • 38

    LiH.DurbinR. (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics26, 589595. 10.1093/bioinformatics/btp698

  • 39

    LiH.HandsakerB.WysokerA.FennellT.RuanJ.HomerN.et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics25, 20782079. 10.1093/bioinformatics/btp352

  • 40

    LiuY.OrsiR. H.BoorK. J.WiedmannM.Guariglia-OropezaV. (2016). An advanced bioinformatics approach for analyzing RNA-seq data reveals sigma H-dependent regulation of competence genes in Listeria monocytogenes. BMC Genomics17:115. 10.1186/s12864-016-2432-9

  • 41

    LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the method. Methods25, 402408. 10.1006/meth.2001.1262

  • 42

    MascherT.HachmannA. B.HelmannJ. D. (2007). Regulatory overlap and functional redundancy among Bacillus subtilis extracytoplasmic function sigma factors. J. Bacteriol.189, 69196927. 10.1128/JB.00904-07

  • 43

    McgannP.WiedmannM.BoorK. J. (2007). The alternative sigma factor σB and the virulence gene regulator PrfA both regulate transcription of Listeria monocytogenes internalins. Appl. Environ. Microbiol.73, 29192930. 10.1128/AEM.02664-06

  • 44

    Medrano RomeroV.MorikawaK. (2016). Listeria monocytogenes SigmaH contributes to expression of competence genes and intracellular growth. J. Bacteriol.198, 12071217. 10.1128/JB.00718-15

  • 45

    MoorheadS. M.DykesG. A. (2003). The role of the sigB gene in the general stress response of Listeria monocytogenes varies between a strain of serotype 1/2a and a strain of serotype 4c. Curr. Microbiol.46, 461466. 10.1007/s00284-002-3867-6

  • 46

    MujahidS.OrsiR. H.BoorK. J.WiedmannM. (2013a). Protein level identification of the Listeria monocytogenes Sigma H, Sigma L, and Sigma C regulons. BMC Microbiol.13:156. 10.1186/1471-2180-13-156

  • 47

    MujahidS.OrsiR. H.VangayP.BoorK. J.WiedmannM. (2013b). Refinement of the Listeria monocytogenes σB regulon through quantitative proteomic analysis. Microbiology159, 11091119. 10.1099/mic.0.066001-0

  • 48

    NilssonL.HansenT. B.GarridoP.BuchrieserC.GlaserP.KnochelS.et al. (2005). Growth inhibition of Listeria monocytogenes by a nonbacteriocinogenic Carnobacterium piscicola. J. Appl. Microbiol.98, 172183. 10.1111/j.1365-2672.2004.02438.x

  • 49

    NussA. M.GlaeserJ.BerghoffB. A.KlugG. (2010). Overlapping alternative sigma factor regulons in the response to singlet oxygen in Rhodobacter sphaeroides. J. Bacteriol.192, 26132623. 10.1128/JB.01605-09

  • 50

    O'byrneC. P.KaratzasK. A. (2008). The role of Sigma B (σB) in the stress adaptations of Listeria monocytogenes: overlaps between stress adaptation and virulence. Adv. Appl. Microbiol.65, 115140. 10.1016/S0065-2164(08)00605-9

  • 51

    OliveA. J.SassettiC. M. (2016). Metabolic crosstalk between host and pathogen: sensing, adapting and competing. Nat. Rev. Microbiol.14, 221234. 10.1038/nrmicro.2016.12

  • 52

    OliverH. F.OrsiR. H.PonnalaL.KeichU.WangW.SunQ.et al. (2009). Deep RNA sequencing of L. monocytogenes reveals overlapping and extensive stationary phase and Sigma B-dependent transcriptomes, including multiple highly transcribed noncoding RNAs. BMC Genomics10:641. 10.1186/1471-2164-10-641

  • 53

    OliverH. F.OrsiR. H.WiedmannM.BoorK. J. (2010). Listeria monocytogenes σB has a small core regulon and a conserved role in virulence but makes differential contributions to stress tolerance across a diverse collection of strains. Appl. Environ. Microbiol.76, 42164232. 10.1128/AEM.00031-10

  • 54

    OllingerJ.BowenB.WiedmannM.BoorK. J.BergholzT. M. (2009). Listeria monocytogenes σB modulates PrfA-mediated virulence factor expression. Infect. Immun.77, 21132124. 10.1128/IAI.01205-08

  • 55

    OllingerJ.WiedmannM.BoorK. J. (2008). Sigma(B)- and PrfA-dependent transcription of genes previously classified as putative constituents of the Listeria monocytogenes PrfA regulon. Foodborne Pathog. Dis.5, 281293. 10.1089/fpd.2008.0079

  • 56

    OrsiR. H.MaronS. B.NightingaleK. K.JeromeM.TaborH.WiedmannM. (2008). Lineage specific recombination and positive selection in coding and intragenic regions contributed to evolution of the main Listeria monocytogenes virulence gene cluster. Infection Genetics and Evolution8, 566576. 10.1016/j.meegid.2008.04.006

  • 57

    RaengpradubS.WiedmannM.BoorK. J. (2008). Comparative analysis of the σB-dependent stress responses in Listeria monocytogenes and Listeria innocua strains exposed to selected stress conditions. Appl. Environ. Microbiol.74, 158171. 10.1128/AEM.00951-07

  • 58

    RaimannE.SchmidB.StephanR.TasaraT. (2009). The alternative sigma factor σL of L. monocytogenes promotes growth under diverse environmental stresses. Foodborne Pathog. Dis.6, 583591. 10.1089/fpd.2008.0248

  • 59

    ReaR. B.GahanC. G.HillC. (2004). Disruption of putative regulatory loci in Listeria monocytogenes demonstrates a significant role for Fur and PerR in virulence. Infect. Immun.72, 717727. 10.1128/IAI.72.2.717-727.2004

  • 60

    RibeiroV. B.MujahidS.OrsiR. H.BergholzT. M.WiedmannM.BoorK. J.et al. (2014). Contributions of σB and PrfA to Listeria monocytogenes salt stress under food relevant conditions. Int. J. Food Microbiol.177, 98108. 10.1016/j.ijfoodmicro.2014.02.018

  • 61

    SaikolappanS.DasK.SasindranS. J.JagannathC.DhandayuthapaniS. (2011). OsmC proteins of Mycobacterium tuberculosis and Mycobacterium smegmatis protect against organic hydroperoxide stress. Tuberculosis (Edinb)91(Suppl. 1), S119S127. 10.1016/j.tube.2011.10.021

  • 62

    ScallanE.HoekstraR. M.AnguloF. J.TauxeR. V.WiddowsonM.-A.RoyS. L.et al. (2011). Foodborne illness acquired in the United States—major pathogens. Emerg. Infect. Dis.17, 715. 10.3201/eid1701.P11101

  • 63

    SchauerK.GeginatG.LiangC.GoebelW.DandekarT.FuchsT. M. (2010). Deciphering the intracellular metabolism of Listeria monocytogenes by mutant screening and modelling. BMC Genomics11:573. 10.1186/1471-2164-11-573

  • 64

    Seifart GomesC.IzarB.PazanF.MohamedW.MraheilM. A.MukherjeeK.et al. (2011). Universal stress proteins are important for oxidative and acid stress resistance and growth of Listeria monocytogenes EGD-e in vitro and in vivo. PLoS ONE6:e24965. 10.1371/journal.pone.0024965

  • 65

    SharmaU. K.ChatterjiD. (2010). Transcriptional switching in Escherichia coli during stress and starvation by modulation of sigma activity. FEMS Microbiol. Rev.34, 646657. 10.1111/j.1574-6976.2010.00223.x

  • 66

    StaibL.FuchsT. M. (2014). From food to cell: nutrient exploitation strategies of enteropathogens. Microbiology160, 10201039. 10.1099/mic.0.078105-0

  • 67

    StollR.MertinsS.JosephB.Müller-AltrockS.GoebelW. (2008). Modulation of PrfA activity in Listeria monocytogenes upon growth in different culture media. Microbiology154, 38563876. 10.1099/mic.0.2008/018283-0

  • 68

    SueD.FinkD.WiedmannM.BoorK. J. (2004). σB-dependent gene induction and expression in Listeria monocytogenes during osmotic and acid stress conditions simulating the intestinal environment. Microbiology150, 38433855. 10.1099/mic.0.27257-0

  • 69

    TangS.OrsiR. H.Den BakkerH. C.WiedmannM.BoorK. J.BergholzT. M. (2015). Transcriptomic analysis of Listeria monocytogenes adaptation to growth on vacuum-packed cold smoked salmon. Appl. Environ. Microbiol. 81, 68126824. 10.1128/AEM.01752-15

  • 70

    Toledo-AranaA.DussurgetO.NikitasG.SestoN.Guet-RevilletH.BalestrinoD.et al. (2009). The Listeria transcriptional landscape from saprophytism to virulence. Nature459, 950956. 10.1038/nature08080

  • 71

    TsaiY. H. L.MaronS. B.McgannP.NightingaleK. K.WiedmannM.OrsiR. H. (2011). Recombination and positive selection contributed to the evolution of Listeria monocytogenes lineages III and IV, two distinct and well supported uncommon L. monocytogenes lineages. Infect. Genet. Evol.11, 18811890. 10.1016/j.meegid.2011.08.001

  • 72

    UtratnaM.ShawI.StarrE.O'byrneC. P. (2011). Rapid, transient, and proportional activation of σB in response to osmotic stress in Listeria monocytogenes. Appl. Environ. Microbiol.77, 78417845. 10.1128/AEM.05732-11

  • 73

    Van SchaikW.ZwieteringM. H.De VosW. M.AbeeT. (2005). Deletion of the sigB Gene in Bacillus cereus ATCC 14579 leads to hydrogen peroxide hyperresistance. Appl. Environ. Microbiol.71, 64276430. 10.1128/AEM.71.10.6427-6430.2005

  • 74

    VivancosA. P.GuellM.DohmJ. C.SerranoL.HimmelbauerH. (2010). Strand-specific deep sequencing of the transcriptome. Genome Res.20, 989999. 10.1101/gr.094318.109

  • 75

    VolkerU.AndersenK. K.AntelmannH.DevineK. M.HeckerM. (1998). One of two osmC homologs in Bacillus subtilis is part of the SigmaB-dependent general stress regulon. J. Bacteriol.180, 42124218.

  • 76

    WadeJ. T.Castro RoaD.GraingerD. C.HurdD.BusbyS. J.StruhlK.et al. (2006). Extensive functional overlap between sigma factors in Escherichia coli. Nat. Struct. Mol. Biol.13, 806814. 10.1038/nsmb1130

  • 77

    WangS. Y.OrsiR. H.TangS. L.ZhangW.WiedmannM.BoorK. J. (2014). Phosphotransferase system-dependent extracellular growth of Listeria monocytogenes is regulated by alternative sigma factors σL and σH. Appl. Environ. Microbiol.80, 76737682. 10.1128/AEM.02530-14

  • 78

    Wemekamp-KamphuisH. H.WoutersJ. A.De LeeuwP. P.HainT.ChakrabortyT.AbeeT. (2004). Identification of sigma factor Sigma B-controlled genes and their impact on acid stress, high hydrostatic pressure, and freeze survival in Listeria monocytogenes EGD-e. Appl. Environ. Microbiol.70, 34573466. 10.1128/AEM.70.6.3457-3466.2004

  • 79

    WinstedtL.Von WachenfeldtC. (2000). Terminal oxidases of Bacillus subtilis strain 168: one quinol oxidase, cytochrome aa3 or cytochrome bd, is required for aerobic growth. J. Bacteriol.182, 65576564. 10.1128/JB.182.23.6557-6564.2000

  • 80

    WrayL. V.FersonA. E.FisherS. H. (1997). Expression of the Bacillus subtilis ureABC operon is controlled by multiple regulatory factors including CodY, GlnR, TnrA, and Spo0H. J. Bacteriol.179, 54945501. 10.1128/jb.179.17.5494-5501.1997

  • 81

    XueJ.MurrietaC. M.RuleD. C.MillerK. W. (2008). Exogenous or L-rhamnose-derived 1,2-propanediol is metabolized via a PduD-dependent pathway in Listeria innocua. Appl. Environ. Microbiol.74, 70737079. 10.1128/AEM.01074-08

  • 82

    YoungM. D.WakefieldM. J.SmythG. K.OshlackA. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol.11:R14. 10.1186/gb-2010-11-2-r14

  • 83

    ZhangC.NietfeldtJ.ZhangM.BensonA. K. (2005). Functional consequences of genome evolution in Listeria monocytogenes: the lmo0423 and lmo0422 genes encode sigmaC and LstR, a lineage II-specific heat shock system. J. Bacteriol.187, 72437253. 10.1128/JB.187.21.7243-7253.2005

  • 84

    ZuberP. (2009). Management of oxidative stress in Bacillus. Annu. Rev. Microbiol.63, 575597. 10.1146/annurev.micro.091208.073241

Summary

Keywords

RNA-seq, Listeria monocytogenes, sigma B, overlapping regulons, promoter

Citation

Liu Y, Orsi RH, Boor KJ, Wiedmann M and Guariglia-Oropeza V (2017) Home Alone: Elimination of All but One Alternative Sigma Factor in Listeria monocytogenes Allows Prediction of New Roles for σB. Front. Microbiol. 8:1910. doi: 10.3389/fmicb.2017.01910

Received

24 June 2017

Accepted

19 September 2017

Published

11 October 2017

Volume

8 - 2017

Edited by

Avelino Alvarez-Ordóñez, Universidad de León, Spain

Reviewed by

Weili Liang, Chinese Center for Disease Control and Prevention, China; Kathrin Rychli, Veterinärmedizinische Universität Wien, Austria

Updates

Copyright

*Correspondence: Veronica Guariglia-Oropeza

This article was submitted to Food Microbiology, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics