Stringent response governs the oxidative stress resistance and virulence of Francisella tularensis

Abstract

Francisella tularensis is a Gram-negative bacterium responsible for causing tularemia in the northern hemisphere. F. tularensis has long been developed as a biological weapon due to its ability to cause severe illness upon inhalation of as few as ten organisms and, based on its potential to be used as a bioterror agent is now classified as a Tier 1 Category A select agent by the CDC. The stringent response facilitates bacterial survival under nutritionally challenging starvation conditions. The hallmark of stringent response is the accumulation of the effector molecules ppGpp and (p)ppGpp known as stress alarmones. The relA and spoT gene products generate alarmones in several Gram-negative bacterial pathogens. RelA is a ribosome-associated ppGpp synthetase that gets activated under amino acid starvation conditions whereas, SpoT is a bifunctional enzyme with both ppGpp synthetase and ppGpp hydrolase activities. Francisella encodes a monofunctional RelA and a bifunctional SpoT enzyme. Previous studies have demonstrated that stringent response under nutritional stresses increases expression of virulence-associated genes encoded on Francisella Pathogenicity Island. This study investigated how stringent response governs the oxidative stress response of F. tularensis. We demonstrate that RelA/SpoT-mediated ppGpp production alters global gene transcriptional profile of F. tularensis in the presence of oxidative stress. The lack of stringent response in relA/spoT gene deletion mutants of F. tularensis makes bacteria more susceptible to oxidants, attenuates survival in macrophages, and virulence in mice. This work is an important step forward towards understanding the complex regulatory network underlying the oxidative stress response of F. tularensis.

About the Authors:

Zhuo Ma

Contributed equally to this work with: Zhuo Ma, Kayla King

Roles Data curation, Formal analysis, Investigation, Methodology, Validation

Affiliation: Department of Basic and Clinical Sciences, Albany College of Pharmacy and Health Sciences, Albany, New York, United States of America

Kayla King

Contributed equally to this work with: Zhuo Ma, Kayla King

Roles Data curation, Investigation, Methodology

Affiliation: Department of Basic and Clinical Sciences, Albany College of Pharmacy and Health Sciences, Albany, New York, United States of America

Maha Alqahtani

Roles Formal analysis, Investigation, Methodology

Affiliation: Department of Microbiology and Immunology, New York Medical College, Valhalla, New York, United States of America

Madeline Worden

Roles Investigation, Methodology

Affiliation: Department of Basic and Clinical Sciences, Albany College of Pharmacy and Health Sciences, Albany, New York, United States of America

Parthasarathy Muthuraman

Roles Investigation, Methodology, Validation

Affiliation: Department of Basic and Clinical Sciences, Albany College of Pharmacy and Health Sciences, Albany, New York, United States of America

Christopher L. Cioffi

Roles Formal analysis, Investigation, Methodology

Affiliation: Department of Basic and Clinical Sciences, Albany College of Pharmacy and Health Sciences, Albany, New York, United States of America

Chandra Shekhar Bakshi

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (MM); [email protected] (CSB)

Affiliation: Department of Microbiology and Immunology, New York Medical College, Valhalla, New York, United States of America

Meenakshi Malik

Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing

* E-mail: [email protected] (MM); [email protected] (CSB)

Affiliation: Department of Basic and Clinical Sciences, Albany College of Pharmacy and Health Sciences, Albany, New York, United States of America

ORCID logo http://orcid.org/0000-0003-0051-4485

Introduction

Francisella tularensis is a Gram-negative bacterium responsible for causing tularemia in the northern hemisphere. F. tularensis has long been developed as a biological weapon due to its ability to cause severe illness upon inhalation of as few as ten organisms and based on its potential to be used as a bioterror agent is now classified as a Tier 1 Category A select agent by the CDC [1–3]. The virulent strains are classified under F. tularensis subsp. tularensis (type A), and F. tularensis subsp. holarctica (type B), whereas avirulent strains belong to F. novicida or F. philomiragia [4]. The virulent SchuS4 strain belongs to F. tularensis subsp. tularensis, while the live vaccine strain (LVS) is derived from F. tularensis subsp. holarctica. Francisella is a facultative intracellular pathogen and can replicate in a variety of cell types; however, macrophages are the primary sites of replication [5,6]. The clinical presentation of tularemia depends on the route, dose, and infecting strain of F. tularensis. The ulceroglandular, oculoglandular, or the typhoidal forms of tularemia are not fatal. However, pneumonic tularemia is a highly acute and fatal form of the disease.

The unique intracellular life cycle of Francisella exposes the bacteria to oxidative stress conditions upon its entry, brief residence in the phagosomes, and escape from phagosomes into the cytosol where replication takes place [5]. F. tularensis genome encodes conventional antioxidant enzymes such as Fe- and CuZn-containing superoxide dismutases (SodB and SodC, respectively), catalase (KatG), and alkyl hydroperoxide reductase (AhpC) homologs to counter the oxidative stress generated at these distinct intracellular locations [7–10]. Unlike other bacterial pathogens, SoxR, an oxidative stress response regulator that regulates the expression of SodB, SodC, and SodA (manganese-containing Sod) is absent in F. tularensis. Instead, another oxidative stress regulator OxyR is present, which regulates the expression of KatG and AhpC in F. tularensis [11]. In F. tularensis, SodB, KatG, and AhpC are induced in response to oxidative stress and are secreted in abundance in the extracellular milieu and into the cytosol of F. tularensis infected macrophages [12]. Both the SodB and KatG are secreted by the Type I Secretion System of F. tularensis [13]. Expression of these primary antioxidant genes starts immediately upon phagocytosis of F. tularensis SchuS4 and remains significantly upregulated during phagosomal and cytosolic phases suggesting that F. tularensis experiences oxidative stress at both of these intracellular locations [14].

The stringent response facilitates bacterial survival under nutritionally challenging starvation conditions. The hallmark of stringent response is the accumulation of the effector molecules ppGpp and (p)ppGpp (guanosine-5'-diphosphate-3'-diphosphate and guanosine-5'-triphosphate-3'-diphosphate) known as stress alarmones [15]. The relA and spoT gene products generate alarmones in several Gram-negative bacterial pathogens. RelA is a ribosome-associated ppGpp synthetase that gets activated under amino acid starvation conditions whereas, SpoT is a bifunctional enzyme with both ppGpp synthetase and ppGpp hydrolase activities. During amino acid starvation, the presence of uncharged tRNA in the acceptor site of the ribosomes sends a signal to the RelA protein associated with ribosomes to catalyze the phosphorylation of GTP in conjunction with ATP as a donor to generate (p)ppGpp. The cytosolic SpoT is required for the basal synthesis of ppGpp during bacterial growth, (p)ppGpp degradation, and elevated synthesis of ppGpp under several stress conditions, including fatty acid and carbon starvation [16]. The ppGpp cause global reprogramming of cellular and metabolic function by binding to the β’-subunit of the RNA polymerase to activate or repress several genes or by interacting directly with proteins to promote adaptation, survival, and transmission in adverse growth conditions.

Francisella encodes a monofunctional relA and a bifunctional spoT gene. A single relA gene deletion mutant of F. novicida grows better than the wild type bacteria, but exhibits attenuated virulence in a mouse model of tularemia [17]. F. tularensis utilizes relA/spoT mediated ppGpp production to promote stable physical interactions between the components of the transcription machinery to activate the expression of virulence-associated genes encoded on Francisella Pathogenicity Island (FPI). It has also been shown that SpoT rather than RelA-dependent production of ppGpp is essential for the expression of virulence genes in F. tularensis [18]. In addition to relA/spoT; migR, trmE, and chpA genes of F. tularensis through unknown mechanisms are also involved in ppGpp production, and regulation of FPI encoded virulence-associated genes [19]. A global gene expression profile of F. tularensis SchuS4 under nutrient limitation conditions induced by serine hydroxamate showed ppGpp-dependent upregulation of genes involved in virulence, metabolism, and stress responses associated with downregulation of genes required for transport and cell division [20]. These studies demonstrate that stringent response under environmental and nutritional stresses increase FPI gene expression. However, the contribution of the stringent response in gene regulation and management of the oxidative stress response when Francisella is experiencing oxidative stress conditions is not known.

In this study, we investigated how stringent response governs the oxidative stress response of F. tularensis. Our results provide a link between the stringent and oxidative stress responses. We demonstrate that relA/spoT-mediated ppGpp production alters the global gene transcriptional profile of F. tularensis in the presence of oxidative stress. The lack of stringent response in relA/spoT gene deletion mutants of F. tularensis makes bacteria more susceptible to oxidants, attenuates survival in macrophages, and virulence in mice. This study provides evidence that the stringent response in Francisella contributes to oxidative stress resistance by regulating the production of several antioxidant enzymes, and enhances our understanding of the intracellular survival mechanisms of F. tularensis.

Materials and methods

Ethics statement

This study was carried out in strict accordance with the recommendations and guidelines of the National Council for Research (NCR) for care and use of animals. All the animal experiments were conducted in the centralized Animal Resources Facility of New York Medical College licensed by the USDA and the NYS Department of Health, Division of Laboratories and Research, and accredited by the American Association for the Accreditation of Laboratory Care. The use of animals and protocols were approved by the Institutional Animal Care and Use Committee (IACUC) of New York Medical College (Protocol Number 69-2-0914H). Mice were administered an anesthetic cocktail consisting of ketamine (5 mg/kg) and xylazine (4 mg/kg) and underwent experimental manipulation only after they failed to exhibit a toe pinch reflex. Mice exhibiting more than 25% weight loss, anorexia, dehydration and impairment of mobility were removed from the study and euthanized by approved means. Humane endpoints were also necessary for mice which survived at the conclusion of the experiments. Mice were administered an anesthetic cocktail of ketamine and xylazine intraperitoneally and then euthanized via cervical dislocation, a method that is consistent with recommendations of the Panel on Euthanasia of the American Veterinary Medical Association. In all experimental procedures, efforts were made to minimize pain and suffering.

Bacterial strains and cultures

F. tularensis subspecies holarctica live vaccine strain (LVS; ATCC 29684; American Culture Collection, Rockville, MD) obtained from BEI Resources, Manassas, VA was used in this study. Bacterial cultures were grown on modified MH-chocolate agar plates (MMH) supplemented with 2% IsoVitaleX. For broth cultures, bacteria were grown in Muller-Hinton Broth (MHB). Working stocks of all bacterial cultures were prepared by growing to mid-log phase at 37°C with 5% CO2 in MHB, snap-frozen in liquid nitrogen, and stored at -80°C until further use. Escherichia coli DH5-α strain (Invitrogen) used for cloning was grown either on Luria-Bertani (LB) broth or LB agar plates (Invitrogen). For the selection of transformants, mutants, or transcomplemented strains, LB-broth or LB-agar media was supplemented with kanamycin (10μg/mL) or hygromycin (200μg/mL). The bacterial strains used in this study are shown in Table 1.

[Figure omitted. See PDF.]

Table 1. List of bacterial strains and plasmid vectors used in this study.

https://doi.org/10.1371/journal.pone.0224094.t001

Generation of F. tularensis LVS ΔrelA and ΔrelAΔspoT gene deletion mutants and transcomplemented strains

In-frame, single-gene deletion mutants of the relA gene (FTL_0285) (ΔrelA), and a double gene deletion mutant of relA and spoT gene (FTL_1413) (ΔrelAΔspoT) of F. tularensis were generated using a previously described allelic replacement method [10,22,23]. For confirmation of the mutants, a duplex colony PCR was performed using internal relA and/or spoT gene-specific primers along with sodB gene-specific primers as internal controls. In-frame deletions of both the relA and spoT genes were confirmed by DNA sequencing of the flanking regions. The deletion strain was designed in a way to preserve the downstream open reading frames to avoid any potential polar effects due to gene deletions.

The ΔrelA gene deletion mutant was transcomplemented with full-length relA gene cloned in the BamHI and Xhol restriction sites of E. coli-Francisella vector pMP882. The ΔrelAΔspoT was transcomplemented with a full-length spoT gene cloned in the pMP822 vector in a similar fashion. The transcomplementation vectors were electroporated into the F. tularensis relA and ΔrelAΔspoT gene deletion mutants and selected on MMH agar supplemented with 200μg/mL hygromycin. The transcomplementation was confirmed by PCR. The primer sequences used for generation and screening of mutants and transcomplementation are shown in Tables 2 and 3.

[Figure omitted. See PDF.]

Table 2. List of primers used for the generation of gene deletion mutants.

https://doi.org/10.1371/journal.pone.0224094.t002

[Figure omitted. See PDF.]

Table 3. List of primers used for the generation of complementation constructs.

https://doi.org/10.1371/journal.pone.0224094.t003

Detection of ppGpp by high-performance liquid chromatography (HPLC)

The ppGpp extraction was done as follows. Briefly, bacterial cultures of F. tularensis LVS, ΔrelA, relA/spoT, and ΔrelAΔspoT mutants were grown to OD600 of 0.5 at 37°C in 10 ml MH-broth. Then 5 ml of cultures were directly mixed with 0.5 ml of 11M formic acid followed by freezing on dry ice and kept at -80°C until use. After that, the mixtures were thawed and allowed to sit on ice for 30 minutes. 1ml aliquots of mixtures were centrifuged at 4C at maximum speed for 5 minutes. The supernatant was then filtered through 0.2μM filters and stored at -80°C until HPLC analysis.

The ppGpp was quantified by HPLC. Isocratic BEH HILIC C18-HPLC column (5 μm, 4.6×150 mm) was used to determine the retention time and absorbance spectrum of a ppGpp standard (TriLink Biosciences). The presence of degradation products—ppGpp, and most likely, pppGpp—was apparent at 2.01 min retention time. The column was run with buffer containing 0.36 M NH4H2PO4, pH 3.4, 2.5% acetonitrile at 26°C with a flow rate of 1.5 mL/minute with a 30 minutes run time. Samples were run for three separate times. Relative levels from samples were represented for comparison.

Bacterial growth curves

Wild type F. tularensis LVS, the ΔrelA mutant, the ΔrelAΔspoT mutant and the corresponding transcomplemented strains were grown on MH chocolate agar plates. The bacterial cultures were then resuspended to an OD600 of 0.05 in MHB. The cultures were either left untreated or treated with 500μM serine hydroxamate, 0.5 and 1mM hydrogen peroxide (H2O2). All bacterial cultures were incubated at 37°C with shaking (150 rpm), and the OD600 values were recorded for 28 hours.

Disc diffusion assays

For disc diffusion assays, bacterial cultures of the wild type F. tularensis LVS, the ΔrelA mutant, the ΔrelAΔspoT mutant, and the corresponding transcomplemented strains grown on MH-chocolate agar plates adjusted to an OD600 of 2.0 in MHB. Two hundred microliters of the bacterial suspensions were spread onto MH-chocolate agar plates using a sterile cotton swab. Sterile filter paper discs were placed onto the agar and were impregnated with 10μL of paraquat (0.88μg/mL), pyrogallol (40μg/mL), menadione (0.39μg/μL), tert-butyl hydroperoxide (TBH; 21.9% solution), and cumene hydroperoxide (CHP; 1.25% solution). The plates were incubated for 72 hours at 37°C and 5% CO2. The zones of inhibition around the discs were measured to determine the sensitivity to the compounds tested.

Bacterial killing assays

For bacterial killing assays, equal numbers of bacteria (1x109 CFU/mL) were diluted in MHB and exposed to H2O2 (1.0mM), paraquat (1.0mM), pyrogallol (1.0mM), or a temperature of 37 or 48°C. Both treated and untreated bacterial suspensions were allowed to incubate for 1 hour and 3 hours. The cultures were then serially diluted 10-fold in PBS and plated on MH-chocolate agar plates and incubated for 72 hours at 37°C with 5% CO2. Viable bacteria were enumerated by counting the colonies and were expressed as Log10 CFU/mL.

Macrophage cell culture assays

RAW macrophage cell line and Bone Marrow-Derived Macrophages (BMDMs) isolated from wild type C57BL/6 or phox-/- mice were infected with wild type F. tularensis LVS, the ΔrelA mutant, the ΔrelAΔspoT mutant, and the corresponding transcomplemented strains at a multiplicity of infection (MOI) of 100. The extracellular bacteria were killed by treating with gentamycin after 2 hours of infection to allow the intracellular bacteria to replicate for 4 or 24 hours. The cells were lysed at both 4- and 24-hours post-infection, diluted 10-fold and plated for enumeration of bacterial colonies.

Mouse survival studies

Mice experiments were conducted at the Animal Resource Facility at New York Medical College according to approved IACUC protocols. All in vivo survival studies were conducted in wild type C57BL/6 mice aged 6 to 8 weeks. Before inoculation, mice were anesthetized with a cocktail of Ketamine and Xylazine. The mice were then intranasally inoculated with doses ranging from 1x105 to 1x106 of the wild type F. tularensis LVS, the ΔrelA mutant, the ΔrelAΔspoT mutant in 20μL PBS. Mice were monitored up to 21 days post-infection for morbidity and mortality. Mice were weighed daily to monitor the infection. The survival results were expressed as Kaplan-Meier survival curves, and the data were analyzed using the Log-rank test.

RNA sequencing

Overnight cultures of wild type F. tularensis LVS and the ΔrelAΔspoT mutant were adjusted to an OD600 of 0.2 and were grown for 3 hours at 37°C with shaking in 10 ml MH-broth in the absence or presence 1mM H2O2. The bacterial cells were pelleted, and total RNA was purified using the Purelink RNA Mini Kit (Ambion). The contaminating DNA from RNA preparations was removed using TURBO DNA-free kits (Invitrogen). The experiments were repeated three times, and cumulative data were analyzed for differential expression of genes. The RNA samples were submitted to the Genomics Core Facility at New York Medical College for RNA sequencing. EDGE-pro software was used to determine expression levels from RNA-seq data. The comparisons were made for differential expression of genes between untreated F. tularensis LVS and the ΔrelAΔspoT mutant (-H2O2) or upon treatment with H2O2 (+H2O2). DESeq2 software was used for pairwise detection and quantification of differential gene expression. A cutoff for multiple testing corrected P-value (adjusted P-value) of <0.05 was used to determine statistical significance and filter differentially expressed genes.

Transcriptional analysis of the target genes

The RNA from wild type F. tularensis LVS, the ΔrelA, ΔrelAΔspoT mutants and the corresponding complemented strains were isolated as described above. cDNA was synthesized using the iScript cDNA Synthesis Kit (Ambion). Quantitative real-time PCR (qPCR) was performed using iQ SYBR Green Supermix (BioRad) to determine the transcriptional levels of genes of interest. The expression of tul4 gene was used as an internal control [18]. Relative levels are represented as fold change and were calculated as follows: 2−ΔΔCT = 2-[ΔCT (mutant) − ΔCT (WT)], where ΔCT = CT target gene − CT internal control. The primer sequences used for qRT-PCR are shown in Table 4.

[Figure omitted. See PDF.]

Table 4. List of primers used for transcriptional analysis by qRT-PCR.

https://doi.org/10.1371/journal.pone.0224094.t004

Statistical analysis

All results were expressed as means ± S.E.M. or S.D. One-Way ANOVA followed by Tukey-Kramer Multiple Comparison tests and Student’s t-test were used for statistical analysis of the data. Survival data were analyzed using the Log-rank test and graphed using Kaplan-Meier survival curves. A P value of less than 0.05 was considered significant.

Results

Generation and characterization of ΔrelA and ΔrelAΔspoT mutants of F. tularensis

The relA gene (FTL_0285) encoding GTP Pyrophosphokinase in Francisella is 647 base pairs in length and is transcribed as a single transcriptional unit; while the spoT gene (FTL_1413) encoding a guanosine-3’ 5’-bis (diphosphate) 3’-phosphohydrolase/(p)ppGpp synthase is 704 bp in length is transcribed as an operon along with capA, capB and capC genes (Fig 1A). To investigate the role of stringent response in oxidative stress resistance of F. tularensis LVS, we constructed unmarked, in-frame relA single gene deletion (ΔrelA) and relA spoT (ΔrelAΔspoT) double gene deletion mutants. The gene deletions in single and double mutants were confirmed by PCR (Fig 1B), followed by DNA sequencing. The transcomplemented strains were constructed by providing full copies of the relA gene in the ΔrelA mutant, and spoT gene in the ΔrelAΔspoT mutant. To determine the ability of the wild type F. tularensis LVS, ΔrelA and the ΔrelAΔspoT mutants to produce ppGpp, the bacterial strains were grown to late exponential phase, and ppGpp production was determined by HPLC. The production of ppGpp was dropped to 45 and 28% in the ΔrelA and the ΔrelAΔspoT mutants respectively, as compared to that observed for the wild type F. tularensis LVS, indicating that ΔrelA and more so the ΔrelAΔspoT mutant is deficient in ppGpp production (Fig 1C and 1D). The growth characteristics of the wild type F. tularensis LVS, the ΔrelA, the ΔrelAΔspoT double mutant and the transcomplemented strains were determined by growing in MH-broth in the absence or presence of 500μM of serine hydroxamate, a serine homolog responsible for inducing amino acid starvation-like conditions. The ΔrelA mutant did not exhibit any growth defect and grew similarly to the wild type F. tularensis LVS when grown in the MH-broth. The ΔrelAΔspoT mutant showed retarded growth and entered the stationary phase after 16 hours of growth. Transcomplementation with spoT gene restored the growth of the ΔrelAΔspoT mutant as well as prevented its early entry into the stationary phase (Fig 1E). Growth of all Francisella strains was slightly reduced in the presence of serine hydroxamate; however, the ΔrelAΔspoT mutant failed to grow in its presence. Transcomplementation of the ΔrelAΔspoT with the spoT gene restored its growth in the presence of serine hydroxamate (Fig 1F). Sensitivities of the wild type F. tularensis LVS, the ΔrelA, the ΔrelAΔspoT mutants and the transcomplemented strains were also tested against streptomycin, nitrofurantoin, and tetracycline by disc diffusion assays. It was observed that the ΔrelAΔspoT mutant exhibited enhanced sensitivities towards all the three antibiotics tested as compared to the wild type F. tularensis LVS or the ΔrelA mutant. Transcomplementation of the ΔrelAΔspoT mutant restored the sensitivities similar to those observed for the wild type or the ΔrelA mutant (Fig 1G). Collectively, these results demonstrate that induction of RelA/SpoT-mediated stringent response associated with the production of ppGpp is required for growth of F. tularensis under normal growth conditions, amino acid starvation, as well as resistance towards antibiotics.

[Figure omitted. See PDF.]

Fig 1. Generation and characterization of ΔrelA and ΔrelAΔspoT mutants of F. tularensis.

(A) The genomic organization of the relA and the spoT genes of the F. tularensis LVS. (B) Multiplex colony PCR using relA and spoT gene-specific primers and sodB gene primers as internal controls. Amplification of the sodB gene confirmed the presence of the DNA template in the reaction, whereas the absence of the relA gene product in the ΔrelA and ΔrelAΔspoT mutants confirmed the deletion of the relA gene and that of spoT confirmed spoT gene deletion in the ΔrelAΔspoT mutant. (C and D) Determination of ppGpp production and its quantitation in the F. tularensis (Ft) LVS, the ΔrelA and the ΔrelAΔspoT mutants by HPLC. (E and F) Growth curves of the Ft LVS, the ΔrelA, the ΔrelAΔspoT mutants, and the transcomplemented strains in the presence of the indicated concentrations of serine hydroxamate. (G) The susceptibility of the Ft LVS, the ΔrelA and the ΔrelAΔspoT mutants and the transcomplemented strains to indicated antibiotics was tested by disk diffusion assay. The plates were incubated for 48–72 hours, and the zones of inhibition around the discs were measured. The results shown in C, D, E, F, and G are representative of three independent experiments with identical results. The data in G are represented as Mean±S.D. (n = 3 biological replicates) and were analyzed using one-way ANOVA (*P<0.05; **P<0.01).

https://doi.org/10.1371/journal.pone.0224094.g001

The transcriptional profile of the ΔrelAΔspoT mutant of F. tularensis

RNA sequencing was used to obtain the transcriptional profiles of the wild type F. tularensis LVS and the ΔrelAΔspoT mutant with or without the treatment with H2O2. Since the ΔrelAΔspoT mutant grows slower than the wild type F. tularensis LVS in the presence of H2O2, 3-hour treatment with H2O2 was selected. At this time point post-treatment, the viability of the mutant strain is not affected significantly. An adjusted P value of <0.05 was used as a cut-off to verify differentially expressed genes in the ΔrelAΔspoT mutant. A total of 318 genes were differentially expressed in untreated ΔrelAΔspoT mutant as compared to the wild type F. tularensis LVS. Of these 163 genes (51%) were downregulated while 155 genes (49%) were upregulated in the ΔrelAΔspoT mutant. The majority of the downregulated genes were involved in metabolism (n = 63), hypothetical proteins (n = 30), Francisella Pathogenicity Island (FPI) genes (n = 30) as well as transport, replication, transcription, translation and stress response genes. The upregulated genes mostly belonged to the hypothetical proteins (n = 65) metabolism (n = 47) and translation (n = 12) categories (Fig 2). Treatment of the ΔrelAΔspoT mutant with H2O2 resulted in differential expression of a total of 855 genes in the ΔrelAΔspoT mutant as compared to the wild type F. tularensis LVS. Of the 450 (53%) downregulated genes, majority of the genes belonged to metabolism (n = 155), hypothetical proteins (n = 95), others (n = 49) and FPI (n = 32) categories. The genes that were upregulated following exposure to H2O2 also belonged to the similar categories except for the FPI genes (Fig 2). The overall gene transcription profiles of the ΔrelAΔspoT mutant as compared to the wild type F. tularensis in the absence or the presence of oxidative stress-induced upon exposure to H2O2 are also shown in S1 and S2 Tables. Collectively these results demonstrate that RelA and SpoT not only control the expression of several genes of F. tularensis under normal growth conditions but also regulate genes both positively and negatively when the bacteria are exposed to oxidative stress.

[Figure omitted. See PDF.]

Fig 2. The transcriptional profile of the ΔrelAΔspoT mutant of F. tularensis.

A total number of differentially expressed genes in the ΔrelAΔspoT mutant compared to the wild-type F. tularensis LVS.

https://doi.org/10.1371/journal.pone.0224094.g002

Expression of the FPI genes in ΔrelAΔspoT mutant of F. tularensis

We further analyzed the differential expression of the FPI genes in the ΔrelAΔspoT mutant under normal as well as oxidative stress conditions created by exposing the bacteria to H2O2. The expression of iglA, B, C, D, pdpE, C, IglI, H, G, F, E, pdpB/icmF and pdpA genes encoded on FPI was significantly downregulated in untreated ΔrelAΔspoT mutant as compared to the wild type F. tularensis LVS. However, the expression of FPI encoded genes iglJ, dotU and vgrG in untreated ΔrelAΔspoT mutant did not alter significantly from that observed in untreated F. tularensis LVS (Fig 3A). The expression of all FPI genes including iglJ, dotU, and vgrG genes was significantly downregulated when the gene expression profile of the ΔrelAΔspoT mutant was compared with F. tularensis LVS upon exposure to H2O2 (Fig 3B). We also confirmed the expression profiles of select FPI genes in the wild type F. tularensis LVS, the ΔrelA, and the ΔrelAΔspoT mutants by quantitative real-time PCR (qRT-PCR). It was observed that deletion of the relA gene affected the expression of iglA as well as the pdpA gene. However, these genes were significantly downregulated in the ΔrelAΔspoT mutant as compared to the wild type F. tularensis LVS and the ΔrelA mutant both in untreated, as well as H2O2 treated bacteria. Transcomplementation of the ΔrelA mutant completely restored the expression levels of the FPI genes, while the expression of these genes was only partially restored in the ΔrelAΔspoT mutant transcomplemented with the spoT gene (Fig 3C and 3D). These results indicate that stringent response induced by RelA/SpoT positively regulates the expression of virulence-associated genes encoded on the FPI under normal as well as oxidative stress conditions.

[Figure omitted. See PDF.]

Fig 3. Expression of the FPI genes in ΔrelAΔspoT mutant of F. tularensis.

RNAseq analysis comparing the differential expression of genes encoded on FPI in untreated (A) and H2O2 treated (B) ΔrelAΔspoT mutant compared to that of wild type F. tularensis (Ft) LVS. The data are represented as Log2 fold change and are cumulative of three independent experiments. The brackets on the right represent operon and genes in red fonts indicate significantly downregulated expression in the ΔrelAΔspoT mutant compared to that of wild type Ft LVS upon treatment with H2O2. (C and D) Quantitative reverse transcriptase PCR (qRT-PCR) was performed to evaluate the transcription of select genes. The amount of target gene amplification was normalized to a tul4 gene internal control. The relative mRNA levels are presented as Mean±SD (n = 3 biological replicates). A and B, DESeq2 software was used for pairwise detection and quantification of differential gene expression. A cutoff for multiple testing corrected P-value (adjusted P-value) of <0.05 was used to determine statistical significance and filter differentially expressed genes. Data shown in C and D were analyzed using ANOVA. (*P<0.05, **P<0.01).

https://doi.org/10.1371/journal.pone.0224094.g003

Expression of stress and heat shock proteins in ΔrelAΔspoT mutant of F. tularensis

Induction of stress-related proteins is a hallmark of the RelA-SpoT-dependent stringent response. We studied the expression profiles of genes involved in degradation and disaggregation accumulated or misfolded proteins, and DNA damage repair mechanisms and heat shock in the ΔrelAΔspoT mutant in the absence or presence of oxidative stress and compared those with the wild type F. tularensis LVS. All genes examined except clpB, hsp90 and hsp40 were downregulated in untreated ΔrelAΔspoT mutant as compared to F. tularensis LVS. However, none of the up- or downregulated genes except the hslU gene achieved statistical significance (Fig 4A). On the contrary, all stress-related genes examined were found to be significantly downregulated in the ΔrelAΔspoT mutant as compared to the wild type F. tularensis LVS when exposed to the oxidative stress (Fig 4B). The differential expression of select genes was also confirmed by qRT-PCR. The expression profile of the representative stress response genes was similar to that observed by RNAseq. However, all these genes were significantly downregulated in untreated as well as H2O2 treated ΔrelAΔspoT mutant as compared to the F. tularensis LVS. Transcomplementation of the ΔrelA mutant restored the wild type phenotype, while only a partial restoration of the gene expression was observed in the transcomplemented ΔrelAΔspoT mutant (Fig 4C and 4D).

[Figure omitted. See PDF.]

Fig 4. Expression of stress and heat shock proteins in ΔrelAΔspoT mutant of F. tularensis.

RNAseq analysis comparing the differential expression of important stress response genes in untreated (A) and H2O2 treated (B) ΔrelAΔspoT mutant compared to that of wild type Ft LVS. The data are represented as Log2 fold change and are cumulative of three independent experiments. Genes in red fonts indicate significantly downregulated expression in the ΔrelAΔspoT mutant compared to that of wild type Ft LVS. A cutoff for multiple testing corrected P-value (adjusted P-value) of <0.05 was used to determine statistical significance and filter differentially expressed genes. (C and D) Quantitative reverse transcriptase PCR (qRT-PCR) was performed to evaluate the transcription of select genes. The amount of target gene amplification was normalized to a tul4 internal control. The relative mRNA levels are presented as Mean±SD (n = 3 biological replicates). (E) A cell viability assay was performed by growing the indicated bacterial strains at 37 and 48°C for 1 and 3 hours. The cultures were diluted 10-fold and spotted on MH-chocolate agar plates. The data in C, D and E were analyzed by ANOVA, and the P values were recorded. *P<0.05; **P<0.01; ***P<0.001.

https://doi.org/10.1371/journal.pone.0224094.g004

Since RelA-SpoT-dependent induction of heat shock proteins is associated with survival at higher temperatures, we also examined the effect of the loss of relA and relA/spoT on bacterial viability when exposed to a higher temperature of 48°C. It was observed that the viability of the ΔrelA mutant was unaffected and remained similar to the wild type F. tularensis when exposed to a temperature of 48°C for 1 hour and 5–7 fold fewer bacteria were recovered after 3 hours of exposure. On the other hand, nearly 10-100-fold less viable bacteria were recovered when the ΔrelAΔspoT mutant was exposed to a temperature of 48°C for 1 and 3 hours, respectively. Transcomplementation of the ΔrelAΔspoT mutant restored the wild type phenotype (Fig 4E). Collectively, these results indicate that RelA-SpoT-mediated stringent response induces protective mechanisms by altering the expression of genes involved in degradation and disaggregation accumulated or misfolded proteins, and DNA damage repair mechanisms under the conditions of oxidative stress. Also, the RelA-SpoT-mediated stringent response facilitates bacterial survival at higher temperatures.

Expression of MglA-dependent genes and transcriptional regulators in ΔrelAΔspoT mutant of F. tularensis

MglA is a major transcriptional regulator of F. tularensis [24]. It has been reported the ppGpp produced by RelA and SpoT plays an important role in MglA-dependent gene regulation of F. tularensis [18]. We next investigated if the transcription of the MglA-dependent genes is altered in untreated or H2O2 treated ΔrelAΔspoT mutant. It was observed that all major MglA-regulated genes were significantly downregulated in untreated ΔrelAΔspoT mutant as compared to untreated F. tularensis LVS. However, the downregulation observed for FTL_0129, FTL_0130, FTL_1546, FTL_1790, and FTL_1876 genes in the ΔrelAΔspoT mutant did not achieve statistical significance. Upon exposure to oxidative stress, all MglA regulated genes examined including FTL_0129, FTL_0130, FTL_1546, FTL_1790, and FTL_1876 genes were significantly downregulated in the ΔrelAΔspoT mutant (Fig 5A and 5B). These results were also confirmed by qRT-PCR for the highly downregulated MglA-dependent gene FTL_1219 (Fig 5C).

[Figure omitted. See PDF.]

Fig 5. Expression of MglA-dependent genes in ΔrelAΔspoT mutant of F. tularensis.

RNAseq analysis comparing the differential expression of MglA regulated genes in untreated (A) and H2O2 treated (B) ΔrelAΔspoT mutant compared to that of wild type Ft LVS. The data are represented as Log2 fold change and are cumulative of three independent experiments. Genes in red fonts indicate significantly downregulated expression in the ΔrelAΔspoT mutant compared to that of wild type Ft LVS. A cutoff for multiple testing corrected P-value (adjusted P-value) of <0.05 was used to determine statistical significance and filter differentially expressed genes. (C) qRT-PCR was performed to evaluate the transcription of a select gene. The amount of target gene amplification was normalized to a tul4 internal control. The relative mRNA levels are presented as Mean±SD (n = 3 biological replicates). The data in C were analyzed by ANOVA, and the P values were recorded. **P<0.01.

https://doi.org/10.1371/journal.pone.0224094.g005

We also investigated if the loss of stringent response mediated by RelA/SpoT affects the expression of other transcription regulators in the absence or the presence of oxidative stress. It was observed that the expression of the majority of the transcriptional regulators was downregulated in untreated ΔrelAΔspoT mutant. However, only fevR/pigR and cphA achieved statistical significance (Fig 6A). On the other hand, expression of all the transcriptional regulator genes examined except migR was significantly downregulated in the ΔrelAΔspoT mutant as compared to the wild type F. tularensis LVS upon exposure to the oxidative stress (Fig 6B). Expression of several regulators including pmrA, FTL_0632, FTL_ 1126, and mglA remained unaltered in untreated ΔrelAΔspoT mutant (Fig 6A and 6C) however, was significantly upregulated in H2O2 treated ΔrelAΔspoT mutant (Fig 6B and 6C). These results were also confirmed by qRT-PCR (Fig 6D and 6E). Taken together, these results demonstrate that RelA-SpoT mediates the differential expression of several transcriptional regulators under normal growth as well as the conditions of oxidative stress.

[Figure omitted. See PDF.]

Fig 6. Expression of transcriptional regulators in ΔrelAΔspoT mutant of F. tularensis.

(A, B and C) RNAseq analysis comparing the differential expression of transcriptional regulator genes in untreated and H2O2 treated ΔrelAΔspoT mutant compared to that of wild type Ft LVS. The data are represented as Log2 fold change and are cumulative of three independent experiments. Genes in red fonts indicate significantly downregulated expression in the ΔrelAΔspoT mutant compared to that of wild type Ft LVS. Genes in green fonts indicate significantly upregulated expression in the ΔrelAΔspoT mutant compared to that of wild type Ft LVS. A cutoff for multiple testing corrected P-value (adjusted P-value) of <0.05 was used to determine statistical significance and filter differentially expressed genes. (D, E) qRT-PCR was performed to evaluate the transcription of select genes. The amount of target gene amplification was normalized to a tul4 internal control. The relative mRNA levels are presented as Mean±SD (n = 3 biological replicates). The data in D and E were analyzed by ANOVA, and the P values were recorded. *P<0.05; **P<0.01.

https://doi.org/10.1371/journal.pone.0224094.g006

Expression of genes involved in antioxidant defense mechanisms in ΔrelAΔspoT mutant of F. tularensis

We next investigated the effect of loss of stringent response in the expression of genes encoding antioxidant enzymes of F. tularensis. Expression of methionine sulfoxide reductase B (msrB), thioredoxin (FTL_0611) glutaredoxin 2 and 1, glutathione peroxidase, katG, and sodB were downregulated in untreated ΔrelAΔspoT mutant. However, except thioredoxin and hypothetical protein genes (FTL_1224 and 1225) expression of none of the other antioxidant enzyme genes achieved statistical significance as compared to those observed for untreated F. tularensis LVS (Fig 7A). The expression of methionine sulfoxide reductase B (msrB), thioredoxin (FTL_0611) glutaredoxin 2 and 1, glutathione peroxidase, FTL_1224, and FTL_1225 was significantly downregulated in the ΔrelAΔspoT mutant when exposed to oxidative stress. The expression of katG and sodB remained unaltered while those of short chain dehydrogenase, sodC and ahpC were significantly upregulated in the ΔrelAΔspoT mutant as compared to the F. tularensis LVS upon exposure to oxidative stress (Fig 7B).

[Figure omitted. See PDF.]

Fig 7. Expression of genes involved in antioxidant defense mechanisms in ΔrelAΔspoT mutant of F. tularensis.

RNAseq analysis comparing the differential expression of antioxidant enzyme genes in untreated (A) and H2O2 (B) treated ΔrelAΔspoT mutant compared to that of wild type Ft LVS. The data are represented as Log2 fold change and are cumulative of three independent experiments. Genes in red fonts indicate significantly downregulated expression in the ΔrelAΔspoT mutant compared to that of wild type Ft LVS. Genes in green fonts indicate significantly upregulated expression in the ΔrelAΔspoT mutant compared to that of wild type Ft LVS. A cutoff for multiple testing corrected P-value (adjusted P-value) of <0.05 was used to determine statistical significance and filter differentially expressed genes. (C) qRT-PCR was performed to evaluate the transcription of select genes. The amount of target gene amplification was normalized to a tul4 internal control. The relative mRNA levels are presented as Mean±SD (n = 3 biological replicates). The data were analyzed by ANOVA, and the P values were recorded. *P<0.05. (D) The western blots of the lysates of the indicated Francisella strains in the absence or presence of H2O2 probed with anti-KatG antibodies, that were stripped and re-probed with antibodies against SodB.

https://doi.org/10.1371/journal.pone.0224094.g007

The expression levels of sodB, katG, thioredoxin 1, and glutaredoxin 2 were also confirmed by qRT-PCR, and the results similar to those obtained from the RNAseq were observed (Fig 7C). We also tested the protein levels of two primary antioxidant enzymes KatG and SodB in bacterial lysates from the wild type F. tularensis LVS, the ΔrelA, ΔrelAΔspoT mutant and the transcomplemented strains by western blot analysis using anti-KatG and anti-SodB antibodies. Diminished protein levels of both KatG and SodB were observed in the ΔrelA mutant and more so in the ΔrelAΔspoT mutant, and these levels were marginally improved in the transcomplemented strains (Fig 7D). However, the levels of both KatG and SodB in the ΔrelA mutant and in the ΔrelAΔspoT mutant were similar to those observed for F. tularensis LVS when treated with H2O2, further confirming the results observed with RNA sequencing. Collectively, these results indicate that the loss of stringent response is associated with differential expression of genes involved in antioxidant defense mechanisms of F. tularensis.

ΔrelAΔspoT mutant of F. tularensis LVS is sensitive to oxidative stress

We next investigated the sensitivity of the ΔrelAΔspoT mutant by performing growth curve analysis, bacterial killing, and disc diffusion assays in the presence of oxidants. As observed earlier, the ΔrelAΔspoT mutant grew slowly and entered the stationary phase earlier than the wild type F. tularensis LVS or the ΔrelA mutant. Transcomplementation of the ΔrelAΔspoT mutant restored its growth (Fig 8A). However, the ΔrelAΔspoT mutant failed to grow when the growth curves were generated in the presence of 0.5 and 1.0 mM H2O2 and again, transcomplementation of the ΔrelAΔspoT mutant restored its growth similar to that observed for the ΔrelA mutant (Fig 8B). We confirmed these findings by performing bacterial killing assays by exposing wild type F. tularensis, the ΔrelA and the ΔrelAΔspoT mutants and the transcomplemented strains to superoxide generating compounds paraquat and pyrogallol for 1 and 3 hours and counting the colonies for bacterial viability. The viability of both the ΔrelA and the ΔrelAΔspoT mutant was similarly affected after 1 hour of treatment. However, the viability of the ΔrelAΔspoT mutant was further reduced after 3 hours of exposure. Nearly 100- and 10-fold less viable ΔrelAΔspoT mutant bacteria treated with paraquat and pyrogallol, respectively, were recovered as compared to those observed for the wild type F. tularensis LVS and the ΔrelA mutant (Fig 8C and 8D). The ΔrelAΔspoT mutant also exhibited enhanced sensitivity as compared to the wild type F. tularensis LVS or the ΔrelA mutant towards diamide and superoxide-generating compounds menadione and organic peroxide TBH when tested by disc diffusion assays (Fig 8E, 8F and 8G). Collectively, these results demonstrate that the loss of relA and spoT is associated with enhanced susceptibility of the ΔrelAΔspoT mutant towards oxidative stress. These results also indicate that the stringent response induced by the RelA and SpoT is also linked with the oxidative stress resistance of F. tularensis.

[Figure omitted. See PDF.]

Fig 8. ΔrelAΔspoT mutant of F. tularensis LVS is sensitive to oxidative stress.

F. tularensis (Ft) LVS, the ΔrelA and ΔrelAΔspoT mutants and the transcomplemented strains were grown in (A) MH-broth or (B) MH-broth containing 500 μM and 1mM of H2O2. The cultures were grown for 28 hours, and OD600 readings were recorded every 4 hours. The bacterial cultures were exposed to (C) 1mM Paraquat or (D) 1mM Pyrogallol for 1 and 3 hours (n = 3 biological replicates). The cultures were diluted 10-fold and plated on MH-chocolate agar plates for bacterial enumeration. The results are expressed as Log10 CFU/ml. The data are representative of three independent experiments conducted with identical results. Disc diffusion assay with diamide (E) and superoxide-generating compounds, menadione (F) and organic peroxide tert-butyl hydroperoxide (TBH) (G). The results are expressed as a zone of inhibitions around the discs impregnated with compounds. The data in C, D, E, F, and G were analyzed by ANOVA, and the P values were recorded. *P<0.05; **P<0.01.

https://doi.org/10.1371/journal.pone.0224094.g008

ΔrelAΔspoT mutant of F. tularensis LVS is attenuated for intramacrophage growth and virulence in mice, and its replication is partially restored in phox-/- BMDMs

We next performed cell culture assays using the Raw 264.7 cell line and BMDMs to investigate if the downregulated expression of the FPI genes is associated with attenuated intramacrophage survival. Except for the ΔrelAΔspoT mutant, equal numbers of bacteria were taken up by the infected Raw cells or BMDMs after 4 hours post-infection. Nearly 5-fold fewer bacteria were taken up by the macrophages infected with the ΔrelAΔspoT mutant at 4-hour time-point. At 24 hours post-infection, the ΔrelA mutant was found to be only partially attenuated for intramacrophage growth. However, nearly 1000-fold less ΔrelAΔspoT mutant bacteria were recovered at 24 hours post-infection as compared to the wild type F. tularensis indicating an attenuation of the intramacrophage growth. Transcomplementation of both the ΔrelA and ΔrelAΔspoT mutants restored the intramacrophage replication (Fig 9A and 9B).

[Figure omitted. See PDF.]

Fig 9. ΔrelAΔspoT mutant of F. tularensis LVS is attenuated for intramacrophage growth and virulence in mice, and its replication is partially restored in phox-/- BMDMs.

RAW macrophage cell line (A) and Bone Marrow-Derived Macrophages (BMDMs) (B) or the BMDMs derived from phox-/- macrophages (C) were infected with wild type F. tularensis (Ft) LVS, the indicated mutants and the transcomplemented strains at a multiplicity of infection (MOI) of 100. The cells were lysed 4 and 24 hours post-infection, diluted 10-fold, and plated for enumeration of bacterial colonies. The data are represented as Mean±SD (n = 3 biological replicates) and are representative of 4 independent experiments. The data were analyzed using ANOVA and P values were determined. **P<0.01; ***P<0.001. (D and E) C57BL/6 mice (n = 5/group) were infected intranasally with 1x105 CFUs of F. tularensis LVS and 1x105 or 1x106 CFUs of the ΔrelA and the ΔrelAΔspoT mutants. A higher challenge dose (1x106 CFUs) of F. tularensis LVS was not used and hence not shown. The mice were observed for mortality and morbidity for 21 days. The body weights of the infected mice were recorded daily to monitor the progress of infection. The data in D were analyzed using the Log-rank test.

https://doi.org/10.1371/journal.pone.0224094.g009

Our preceding results demonstrated that the ΔrelAΔspoT mutant is sensitive to oxidants and attenuated for intramacrophage growth. We next investigated if the replication-deficient phenotype of the ΔrelAΔspoT mutant is restored by infecting phox-/- macrophages which are incapable of generating reactive oxygen species. It was observed that nearly 100-fold higher ΔrelAΔspoT mutant bacteria were recovered from the phox-/-macrophages than those recovered from the wild type macrophages, indicating that the mutant bacteria can replicate in the absence of oxidative stress (Fig 9C). These results demonstrate that the stringent response mediated by RelA/SpoT of F. tularensis is required to overcome the oxidative stress in macrophages to establish its intracellular replicative niche.

We also investigated if the loss of relA and relA/spoT is associated with the attenuation of virulence in mice. It was observed that 100% of C57BL/6 mice infected intranasally either with 1×105 or 1×106 CFUs of the ΔrelA or the ΔrelAΔspoT mutant survived the infection. Mice infected with the ΔrelA mutant exhibited initial loss of body weight from days 4–8 post-infection and then regained their initial body weight. On the other hand, mice infected with the ΔrelAΔspoT mutant did not show any bodyweight loss. All the control mice infected with 1×105 CFU of F. tularensis LVS succumbed to infection by day ten post-infection and experienced severe body weight loss (Fig 8D and 8E). Collectively, these results demonstrate that RelA and SpoT mediated stringent response contributes to virulence in mice.

Discussion

The synthesis of ppGpp by RelA in Gram-negative bacteria on the stalled ribosome diminishes the synthesis of stable RNA and ribosomal proteins, reprogram transcription machinery by inducing the expression of alternative sigma factors and enhance their association with core RNA polymerases. This response described as stringent response alters the half-life of RNA polymerase–promoter complexes to up- or downregulate the expression of several genes [16]. The SpoT is a bifunctional synthetase and hydrolase that regulates the intracellular concentrations of (p)ppGpp. The SpoT-dependent accumulation of ppGpp occurs under carbon, iron, and fatty acid starvation conditions [15]. The SpoT in-turn, is regulated by ribosome-associated proteases and acyl carrier proteins [25,26]. In Francisella, although not mechanistically demonstrated, it has been proposed that SpoT is regulated by migR, trmE, and chpA genes [19]. These genes, through unidentified mechanisms, are involved in SpoT-mediated ppGpp production, and regulation of FPI encoded virulence-associated genes. An additional protein mediates the promoter interaction that affects the gene transcription in E. coli is known as DksA [27]. The DksA homolog is absent in F. tularensis.

A very limited number of studies conducted to date on the stringent response of Francisella have primarily been focused on understanding the regulation of virulence gene expression encoded on the FPI. The ΔrelA mutant of F. novicida fails to produce ppGpp under the conditions of amino acid starvation but grows better than the wild type strain at different stages of growth, including the stationary phase [17]. Moreover, the ΔrelA mutant of F. novicida exhibits reduced viability at higher temperatures and attenuated virulence in mice. The expression of virulence genes requires the association of MglA and SspA proteins with RNA polymerase to form a complex which acts in conjunction with a DNA binding protein known as PigR. Additionally, RelA/SpoT-mediated production of ppGpp is required for stable interactions of PigR with MglA-SspA-RNA polymerase complex for the expression of virulence genes [18]. Similarly, MigR, TrmE, and CphA proteins of F. tularensis are required for accumulation ppGpp and virulence gene expression [19]. Transcriptome analysis of F. tularensis SchuS4 under starvation-like conditions induced artificially by exposing Francisella to serine hydroxamate has shown upregulation of FPI encoded genes involved in virulence, stress response, and metabolism and downregulation of genes involved in transport and cell division [20]. Collectively, this handful of studies have demonstrated the association of nutritional stress with virulence gene expression in Francisella.

Herein, we report a critical role of the stringent response in oxidative stress resistance of F. tularensis. The unique intracellular lifestyle of Francisella in addition to the nutritional stress also exposes bacteria to oxidative stress. However, how stringent response governs the oxidative stress response of F. tularensis is not known. To address this question, we adopted a genetic approach and generated a single gene deletion mutant ΔrelA and a ΔrelAΔspoT double gene mutant. Unlike previous studies that mostly used amino acid starvation as a means to induce stringent response, in this study, we investigated the role of stringent response in a nutritionally rich environment and when the bacteria are exposed to oxidative stress. We observed that ppGpp production was reduced drastically in both the ΔrelA and the ΔrelAΔspoT mutants as compared to wild type F. tularensis LVS even when the bacteria were still in the exponential phase of growth and not exposed to any stress, indicating that expression of ppGpp is also required under homeostatic growth conditions. Our initial characterization revealed that the ΔrelAΔspoT but not the ΔrelA mutant has a growth defect, enter stationary phase early, and exhibit high resistance towards streptomycin, nitrofurantoin, and tetracycline. These phenotypic attributes, although have been reported for the ΔrelAΔspoT mutants of several other Gram-negative bacteria but have not been reported for the ΔrelAΔspoT mutant of F. tularensis LVS.

Our results show that a large number of genes are regulated by RelA/SpoT-dependent ppGpp even when bacteria are growing in a rich environment and not exposed to any nutritional stress. The majority of the genes that are regulated under the homeostatic growth conditions included those involved in metabolism, transport, and FPI genes. A number of genes encoding hypothetical proteins were also both up-and down-regulated in the ΔrelAΔspoT mutant. These results indicate that stringent response plays an essential role in the growth and survival of F. tularensis even in a nutritionally rich environment. When the bacteria were exposed to oxidative stress, twice the numbers of genes belonging to similar categories were differentially expressed in the ΔrelAΔspoT mutant indicating that stringent response regulates the oxidative stress response experienced by the bacteria in a nutritionally rich environment.

Francisella encounters oxidative stress in phagosomes as well as in the cytosol of the phagocytic cells during its intracellular residence. These results suggest that the expression of T6SS components encoded on FPI is regulated by stringent response under normal growth conditions as well as when the bacteria are exposed to oxidative stress. Specifically, the expression of genes that constitute the membrane complex of T6SS such as pdpB, iglE and dotU and the genes that encode T6SS effector proteins such as pdpA, vgrG, iglG, and iglF is regulated by the stringent response when the bacteria sense oxidative stress conditions. Coherent with the expression profile of the FPI genes, the ΔrelAΔspoT mutant but not the ΔrelA was severely attenuated for intramacrophage growth. It was observed that both the ΔrelA and the ΔrelAΔspoT mutants were highly attenuated for virulence upon intranasal infection in C57BL/6 mice. A previous study has reported the attenuation of ΔrelAΔspoT mutant of F. tularensis LVS in BALB/c mice upon intradermal challenge [18].

Similar to the FPI genes, the general stress proteins were also differentially expressed in the ΔrelAΔspoT mutant of F. tularensis LVS and three transcriptional profiles emerged when the bacteria were exposed to oxidative stress. The expression profile of clpB, hsp90, and hsp40 genes involved in disaggregation and reactivation of the aggregated proteins, the stress response genes that serve as disaggregation or degradation machines and are transcribed as operons (clpP, clpX, and lon; hslU, hslV; grpE, dnaK, and dnaJ) and usp, sohB and starvation protein A genes involved in response to DNA damage were significantly downregulated when ΔrelAΔspoT mutant bacteria were exposed to oxidative stress. These results indicate that RelA/SpoT-dependent stringent response regulates the expression of the stress response genes in the presence of oxidative stress to facilitate bacterial survival. The differential expression of these stress response genes was also reflected in the enhanced sensitivity of the ΔrelAΔspoT mutant towards high temperature. Contrary to what has been reported for the ΔrelA mutant of F. novicida [17], the viability of the ΔrelA mutant of F. tularensis LVS is not affected by exposure to a higher temperature and remains similar to its wild type counterpart.

MglA is a master regulator of F. tularensis [24,28,29]. Several MglA-regulated genes were downregulated in the ΔrelAΔspoT mutant and remained similarly downregulated under oxidative stress conditions indicating that positive regulation of the MglA-dependent genes by RelA/SpoT is independent of the oxidative stress. In addition to MglA, in this study, it was also observed that a number of transcriptional regulators are also regulated both positively and negatively by RelA/SpoT. The positive regulation of the fevR and cphA genes is independent of the oxidative stress. FevR regulates the virulence gene expression of F. tularensis [30]; while cphA regulate FPI genes as well as the production of ppGpp [19]. On the other hand, a LySR family transcriptional regulator (FTL_1193), fur and oxyR are regulated positively by RelA/SpoT only under oxidative stress conditions. However, the transcription of migR and mglA genes is independent of RelA/SpoT-dependent stringent response. The expression of several transcriptional regulators, including pmrA, is negatively regulated by RelA/SpoT under oxidative stress conditions. Collectively, these results indicate that stringent response mediated by RelA/SpoT regulates expression of a number of transcriptional regulators under oxidative stress response and thus, serves as a master regulator to overcome oxidative stress.

Induction of stringent response is associated with increased transcription of oxidative stress response genes in Gram-positive bacterial pathogens such as Enterococcus faecalis, Streptococci and Staphylococcus aureus [31–34]. A similar association has also been demonstrated in Gram-negative pathogen Pseudomonas aeruginosa. It has been reported that stringent response mediated oxidative stress response protects P. aeruginosa against antibiotics exerting their bactericidal activities primarily by the induction of oxidative stress [35,36]. Recent reports have shown that ppGpp-mediated stringent response is required for the optimal catalase activity that provides H2O2 tolerance and protects from antibiotic-mediated killing of nutrient-starved P. aeruginosa during planktonic and biofilm growth [37,38]. These studies demonstrated a link between stringent and oxidative stress responses.

In the present study, it was observed that RelA/SpoT positively regulated a number of genes encoding antioxidant enzymes in F. tularensis. The most prominent ones being the methionine sulfoxide reductase B (msrB), thioredoxin, glutaredoxin 1 and 3 under oxidative stress conditions. In E. coli, the trx1 and grx2 genes are positively regulated by ppGpp during the stationary phase of growth [39]. On the other hand, expression of primary antioxidant enzyme genes katG and sodB are not regulated by the stringent response. Moreover, the result indicates that expression of sodC and ahpC is negatively regulated by RelA/SpoT under the conditions of oxidative stress. Superoxide dismutase B (SodB) and catalase (KatG) activities are regulated by the stringent response in Pseudomonas aeruginosa [37,38]; however, the expression profiles of both KatG and SodB under oxidative stress conditions are not known. Similar to P. aeruginosa, we observed that in addition to the transcript levels, the protein levels of primary antioxidant enzymes SodB and KatG of F. tularensis were reduced in untreated ΔrelA and ΔrelAΔspoT mutants. However, exposure to H2O2 restored the transcript as well as the protein levels of SodB and KatG, similar to those observed for the wild type F. tularensis LVS, indicating that under oxidative stress conditions, the expression of these primary antioxidant enzymes is not regulated by the stringent response. Despite the levels of primary antioxidant enzymes SodB, KatG, SodC, and AhpC remained unaltered or slightly elevated under oxidative stress, the ΔrelAΔspoT mutant exhibited enhanced sensitivity towards peroxides and superoxide-generating compounds. These observations indicate that thioredoxin and glutaredoxin systems, which are regulated by stringent response under oxidative stress conditions play an essential role in maintaining the redox-homeostasis and serves as an oxidative stress defense mechanism. Collectively, these observations indicate that stringent response of F. tularensis regulates the expression of antioxidant enzymes to prevent damage from oxidative stress.

To conclude, the unique intracellular lifestyle of Francisella in addition to the nutritional stress, also exposes bacteria to oxidative stress. This study demonstrates how stringent response governs the complex oxidative stress response of F. tularensis by regulating FPI genes, stress response genes, transcriptional regulators, and antioxidant defense mechanisms. This study provides evidence that in addition to nutritional stresses, the stringent response also plays an essential role in the oxidative stress resistance of F. tularensis.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Gene expression profile of untreated (-H2O2) ΔrelAΔspoT mutant as compared to the wild type F. tularensis LVS as determined by RNAseq analysis.

https://doi.org/10.1371/journal.pone.0224094.s001

(PDF)

S2 Table. Gene expression profile of ΔrelAΔspoT mutant as compared to the wild type F. tularensis LVS following treatment with hydrogen peroxide (+H2O2) as determined by RNAseq analysis.

https://doi.org/10.1371/journal.pone.0224094.s002

(PDF)

Citation: Ma Z, King K, Alqahtani M, Worden M, Muthuraman P, Cioffi CL, et al. (2019) Stringent response governs the oxidative stress resistance and virulence of Francisella tularensis. PLoS ONE 14(10): e0224094. https://doi.org/10.1371/journal.pone.0224094

References

1. Pechous RD, McCarthy TR, Zahrt TC. Working toward the future: insights into Francisella tularensis pathogenesis and vaccine development. Microbiol Mol Biol Rev. 2009;73: 684–711. pmid:19946137

2. Eisen L. Francisella tularensis : Biology, Pathogenicity, Epidemiology, and Biodefense. Emerg Infect Dis. 2007;13: 1973–1973.

3. Hirschmann JV. From Squirrels to Biological Weapons: The Early History of Tularemia. Am J Med Sci. 2018;356: 319–328. pmid:30146078

4. Kingry LC, Petersen JM. Comparative review of Francisella tularensis and Francisella novicida. Front Cell Infect Microbiol. 2014;4: 35. pmid:24660164

5. Santic M, Molmeret M, Klose KE, Abu Kwaik Y. Francisella tularensis travels a novel, twisted road within macrophages. Trends in Microbiology. 2006. pp. 37–44. pmid:16356719

6. Hall JD, Woolard MD, Gunn BM, Craven RR, Taft-Benz S, Frelinger JA, et al. Infected-host-cell repertoire and cellular response in the lung following inhalation of Francisella tularensis Schu S4, LVS, or U112. Infect Immun. 2008; pmid:18852251

7. Bakshi CS, Malik M, Regan K, Melendez JA, Metzger DW, Pavlov VM, et al. Superoxide dismutase B gene (sodB)-deficient mutants of Francisella tularensis demonstrate hypersensitivity to oxidative stress and attenuated virulence. J Bacteriol. Center for Immunology and Microbial Disease, Albany Medical College, 47 New Scotland Avenue, Albany, NY 12208–3479, USA; 2006;188: 6443–6448. pmid:16923916

8. Bakshi CS, Mahawar M, Melendez JA, Sellati TJ, Malik M, Metzger DW, et al. Identification of Francisella tularensis Live Vaccine Strain CuZn Superoxide Dismutase as Critical for Resistance to Extracellularly Generated Reactive Oxygen Species. J Bacteriol. Center for Immunology and Microbial Disease, MC 151, Albany Medical College, 47 New Scotland Ave., Albany, NY 12208, USA; 2009;191: 6447–6456. pmid:19684141

9. Lindgren H, Shen H, Zingmark C, Golovliov I, Conlan W, Sjöstedt A. Resistance of Francisella tularensis strains against reactive nitrogen and oxygen species with special reference to the role of KatG. Infect Immun. Department of Clinical Microbiology, Clinical Bacteriology, Umea University, SE-901 85 Umea, Sweden. [email protected]; 2007;75: 1303–1309. pmid:17210667

10. Alharbi A, Rabadi SM, Alqahtani M, Marghani D, Worden M, Ma Z, et al. Role of peroxiredoxin of the AhpC/TSA family in antioxidant defense mechanisms of Francisella tularensis. PLoS One. 2019; pmid:30870480

11. Malik M, Bakshi CS, Jen Y, Rabadi SM, Catlett S V., Ma Z, et al. Elucidation of a mechanism of oxidative stress regulation in Francisella tularensis live vaccine strain. Mol Microbiol. 2016;101: 856–878. pmid:27205902

12. Clemens DL, Lee BY, Horwitz MA. Francisella tularensis enters macrophages via a novel process involving pseudopod loops. Infect Immun. Division of Infectious Diseases, Dept. of Medicine, UCLA School of Medicine, CHS 37–121, 10833 LeConte Ave., Los Angeles, CA 90095–1688, USA. [email protected]; 2005;73: 5892–5902. pmid:16113308

13. Ma Z, Banik S, Rane H, Mora VT, Rabadi SM, Doyle CR, et al. EmrA1 membrane fusion protein of F rancisella tularensis LVS is required for resistance to oxidative stress, intramacrophage survival and virulence in mice. Mol Microbiol. Department of Basic and Social Sciences, Albany College of Pharmacy and Health Sciences, Albany, New York, USA; 2014;91: 976–995. pmid:24397487

14. Wehrly TD, Chong A, Virtaneva K, Sturdevant DE, Child R, Edwards JA, et al. Intracellular biology and virulence determinants of Francisella tularensis revealed by transcriptional profiling inside macrophages. Cell Microbiol. 2009;11: 1128–1150. pmid:19388904

15. Potrykus K, Cashel M. (p)ppGpp: Still Magical? Annu Rev Microbiol. 2008; pmid:18454629

16. Cashel M, Potrykus K. Stringent Response. Brenner’s Encyclopedia of Genetics: Second Edition. 2013. https://doi.org/10.1016/B978-0-12-374984-0.01486-8

17. Dean RE, Ireland PM, Jordan JE, Titball RW, Oyston PCF. RelA regulates virulence and intracellular survival of Francisella novicida. Microbiology. 2009;155: 4104–4113. pmid:19762448

18. Charity JC, Blalock LT, Costante-Hamm MM, Kasper DL, Dove SL. Small molecule control of virulence gene expression in Francisella tularensis. PLoS Pathog. 2009;5: e1000641. pmid:19876386

19. Faron M, Fletcher JR, Rasmussen JA, Long ME, Allen L-AH, Jones BD. The Francisella tularensis migR, trmE, and cphA Genes Contribute to F. tularensis Pathogenicity Island Gene Regulation and Intracellular Growth by Modulation of the Stress Alarmone ppGpp. Infect Immun. 2013; pmid:23716606

20. Murch AL, Skipp PJ, Roach PL, Oyston PCF. Whole genome transcriptomics reveals global effects including up-regulation of Francisella pathogenicity island gene expression during active stringent response in the highly virulent Francisella tularensis subsp. tularensis SCHU S4. Microbiology. 2017;163: 1664–1679. pmid:29034854

21. Lovullo ED, Sherrill LA, Pavelka MS. Improved shuttle vectors for Francisella tularensis genetics: Research Letter. FEMS Microbiol Lett. 2009;291: 95–102. pmid:19067747

22. Alqahtani M, Ma Z, Ketkar H, Suresh RV, Malik M, Bakshi CS. Characterization of a Unique Outer Membrane Protein Required for Oxidative Stress Resistance and Virulence of Francisella tularensis. J Bacteriol. 2018; JB.00693-17. pmid:29378894

23. Malik M, Bakshi CS, Jen Y, Rabadi SM, Catlett S V., Ma Z, et al. Elucidation of a mechanism of oxidative stress regulation in Francisella tularensis live vaccine strain. Mol Microbiol. 2016;101: 856–878. pmid:27205902

24. Lauriano CM, Barker JR, Yoon S-S, Nano FE, Arulanandam BP, Hassett DJ, et al. MglA regulates transcription of virulence factors necessary for Francisella tularensis intraamoebae and intramacrophage survival. Proc Natl Acad Sci U S A. 2004;101: 4246–9. pmid:15010524

25. Battesti A, Bouveret E. Acyl carrier protein/SpoT interaction, the switch linking SpoT-dependent stress response to fatty acid metabolism. Mol Microbiol. 2006; pmid:17078815

26. Battesti A, Bouveret E. Bacteria possessing two RelA/SpoT-like proteins have evolved a specific stringent response involving The acyl carrier protein-SpoT interaction. J Bacteriol. 2009; pmid:18996989

27. Kim J-S, Liu L, Fitzsimmons LF, Wang Y, Crawford MA, Mastrogiovanni M, et al. DksA–DnaJ redox interactions provide a signal for the activation of bacterial RNA polymerase. Proc Natl Acad Sci. 2018; pmid:30429329

28. Guina T, Radulovic D, Bahrami AJ, Bolton DL, Rohmer L, Jones-Isaac KA, et al. MglA regulates Francisella tularensis subsp. novicida (Francisella novicida) response to starvation and oxidative stress. J Bacteriol. 2007;189: 6580–6586. pmid:17644593

29. Brotcke A, Weiss DS, Kim CC, Chain P, Malfatti S, Garcia E, et al. Identification of MglA-regulated genes reveals novel virulence factors in Francisella tularensis. Infect Immun. 2006;74: 6642–6655. pmid:17000729

30. Brotcke A, Monack DM. Identification of fevR, a novel regulator of virulence gene expression in Francisella novicida. Infect Immun. 2008;76: 3473–3480. pmid:18559431

31. Gaca AO, Abranches J, Kajfasz JK, Lemos JA. Global transcriptional analysis of the stringent response in Enterococcus faecalis. Microbiology. 2012;158: 1994–2004. pmid:22653948

32. Reiß S, Pané-Farré J, Fuchs S, François P, Liebeke M, Schrenzel J, et al. Global Analysis of the Staphylococcus aureus Response to Mupirocin. Antimicrob Agents Chemother. 2012;56: 787–804. pmid:22106209

33. Kazmierczak KM, Wayne KJ, Rechtsteiner A, Winkler ME. Roles of rel Spn in stringent response, global regulation and virulence of serotype 2 Streptococcus pneumoniae D39. Mol Microbiol. 2009;72: 590–611. pmid:19426208

34. Nascimento MM, Lemos JA, Abranches J, Lin VK, Burne RA. Role of RelA of Streptococcus mutans in Global Control of Gene Expression. J Bacteriol. 2008;190: 28–36. pmid:17951382

35. Nguyen D, Joshi-Datar A, Lepine F, Bauerle E, Olakanmi O, Beer K, et al. Active starvation responses mediate antibiotic tolerance in biofilms and nutrient-limited bacteria. Science (80-). 2011; pmid:22096200

36. Kohanski MA, Dwyer DJ, Hayete B, Lawrence CA, Collins JJ. A Common Mechanism of Cellular Death Induced by Bactericidal Antibiotics. Cell. 2007;130: 797–810. pmid:17803904

37. Khakimova M, Ahlgren HG, Harrison JJ, English AM, Nguyen D. The stringent response controls catalases in Pseudomonas aeruginosa and is required for hydrogen peroxide and antibiotic tolerance. J Bacteriol. 2013; pmid:23457248

38. Martins D, McKay G, Sampathkumar G, Khakimova M, English AM, Nguyen D. Superoxide dismutase activity confers (p)ppGpp-mediated antibiotic tolerance to stationary-phase Pseudomonas aeruginosa. Proc Natl Acad Sci. 2018;115: 9797–9802. pmid:30201715

39. Potamitou A, Neubauer P, Holmgren A, Vlamis-Gardikas A. Expression of Escherichia coli Glutaredoxin 2 Is Mainly Regulated by ppGpp and ς S. J Biol Chem. 2002;277: 17775–17780. pmid:11889138