Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Molecular Biology and Physiology

A Gene Cluster That Encodes Histone Deacetylase Inhibitors Contributes to Bacterial Persistence and Antibiotic Tolerance in Burkholderia thailandensis

Sofiya N. Micheva-Viteva, Migun Shakya, Samantha H. Adikari, Cheryl D. Gleasner, Nileena Velappan, Judith R. Mourant, Patrick S. G. Chain, Elizabeth Hong-Geller
Christopher W. Marshall, Editor
Sofiya N. Micheva-Viteva
aBioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Migun Shakya
aBioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Samantha H. Adikari
aBioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cheryl D. Gleasner
aBioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nileena Velappan
aBioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Judith R. Mourant
aBioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Patrick S. G. Chain
aBioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elizabeth Hong-Geller
aBioscience Division, Los Alamos National Laboratory, Los Alamos, New Mexico, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christopher W. Marshall
Marquette University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00609-19
  • Article
  • Figures & Data
  • Info & Metrics
  • Review History
  • PDF
Loading

ABSTRACT

Persister cells are genetically identical variants in a bacterial population that have phenotypically modified their physiology to survive environmental stress. In bacterial pathogens, persisters are able to survive antibiotic treatment and reinfect patients in a frustrating cycle of chronic infection. To better define core persistence mechanisms for therapeutics development, we performed transcriptomics analyses of Burkholderia thailandensis populations enriched for persisters via three methods: flow sorting for low proton motive force, meropenem treatment, and culture aging. Although the three persister-enriched populations generally displayed divergent gene expression profiles that reflect the multimechanistic nature of stress adaptations, there were several common gene pathways activated in two or all three populations. These include polyketide and nonribosomal peptide synthesis, Clp proteases, mobile elements, enzymes involved in lipid metabolism, and ATP-binding cassette (ABC) transporter systems. In particular, identification of genes that encode polyketide synthases (PKSs) and fatty acid catabolism factors indicates that generation of secondary metabolites, natural products, and complex lipids could be part of the metabolic program that governs the persistence state. We also found that loss-of-function mutations in the PKS-encoding gene locus BTH_I2366, which plays a role in biosynthesis of histone deacetylase (HDAC) inhibitors, resulted in increased sensitivity to antibiotics targeting DNA replication. Furthermore, treatment of multiple bacterial pathogens with a fatty acid synthesis inhibitor, CP-640186, potentiated the efficacy of meropenem against the persister populations. Altogether, our results suggest that bacterial persisters may exhibit an outwardly dormant physiology but maintain active metabolic processes that are required to maintain persistence.

IMPORTANCE The discovery of antibiotics such as penicillin and streptomycin marked a historic milestone in the 1940s and heralded a new era of antimicrobial therapy as the modern standard for medical treatment. Yet, even in those early days of discovery, it was noted that a small subset of cells (∼1 in 105) survived antibiotic treatment and continued to persist, leading to recurrence of chronic infection. These persisters are phenotypic variants that have modified their physiology to survive environmental stress. In this study, we have performed three transcriptomic screens to identify persistence genes that are common between three different stressor conditions. In particular, we identified genes that function in the synthesis of secondary metabolites, small molecules, and complex lipids, which are likely required to maintain the persistence state. Targeting universal persistence genes can lead to the development of clinically relevant antipersistence therapeutics for infectious disease management.

INTRODUCTION

Many microbes can undergo phenotypic changes and enter a seemingly dormant metabolic state, termed bacterial persistence, as a survival mechanism during periods of environmental stress. In bacteria, these persisters become recalcitrant to lethal doses of antibiotic and subsequently contribute to the onset of chronic infectious disease. Although persisters do not pass antibiotic tolerance to their progeny through genetic inheritance, persistence can exacerbate antibiotic resistance by decreasing antibiotic efficacy, prolonging treatment time, and generating a reservoir of surviving cells from which antibiotic-resistant mutations can evolve (1–3).

Both stochastic and deterministic mechanisms are thought to contribute to bacterial persistence. In the former case, random fluctuations in gene transcription, translation, or protein stability lead to phenotypic variability within bacterial populations (4). Alternatively, in the latter case, phenotypic divergence could be predetermined by epigenetic memory carried in a small subpopulation of cells that slowly recover from nutrient starvation, exposure to toxins, or harsh environmental conditions. The role of stress has been implicated in the establishment of persistence, including the stringent response to starvation, DNA damage, SOS signaling, and oxidative stress. Genes that encode stress response mechanisms were activated in persisters enriched by exposure to antibiotics (5–8). Furthermore, low acidity, temperature variations, and osmotic stress were also shown to significantly increase the number of persisters in bacteria (9–13).

Multiple studies have identified a variety of genes that play a role in persistence. Disrupted equilibrium of toxin-antitoxin (TA) systems, which regulate essential cellular activities, including protein translation, DNA stability, and cell membrane integrity, leads to metabolic shutdown similar to that of the stringent response to stress (14–18). Other studies have profiled persister transcriptomes using various enrichment methods, including antibiotic lysis of susceptible bacterial populations (5, 19, 20), flow cytometry based on dormancy markers (6), and analysis of stationary-phase bacteria (8, 11, 21). However, the already highly heterogeneous nature of persisters has been further confounded by the different stressor conditions used between the experimental groups, thus generating transcriptomic data sets that were biased to the particular phenotypic assay. Currently, there is no single essential persistence gene or mechanism that can be targeted for development of antipersistence drugs.

To identify core persistence genes that are ubiquitously activated in response to various stress conditions, we performed transcriptomics analyses of Burkholderia thailandensis persisters enriched by three experimental conditions, flow sorting for low proton motive force (PMF), meropenem (Mpm) treatment, and culture aging. B. thailandensis E246, a nonpathogenic saprophyte, shares close genetic similarity with two highly infectious pathogens, Burkholderia mallei and Burkholderia pseudomallei, classified as select agents due to their high mortality rates and ability to establish chronic infections with high percentage of relapse (22). Although lacking the virulence factors found in its two closely related pathogens, B. thailandensis E264 carries highly conserved genes that regulate antibiotic resistance and stress adaptation mechanisms (23) and is used as a surrogate for B. pseudomallei in infection studies in vitro. We had previously shown that B. thailandensis harbors a high level of persisters within logarithmically expanding populations (24). Furthermore, Burkholderia species encode various stress adaptation mechanisms that can drive the persistence phenotype, including efflux pumps, biofilm formation, quorum sensing systems, and various catabolite repression/activation systems known to induce the stringent response (25–28).

In this study, we identified genes encoding polyketide synthases (PKSs) and fatty acid catabolism factors that are activated in persisters enriched via various stress conditions, suggesting that generation of secondary metabolites, natural products, and complex lipids can contribute to resiliency to antibiotics. These B. thailandensis PKS systems have been previously characterized as factories that generate a novel class of histone deacetylase (HDAC) inhibitors (29). Given the role of HDAC inhibitors as suppressors of the host immune response to microbial infections (30), our study has uncovered a link between persistence and the establishment of chronic infections. Identification of universal persistence genes or mechanisms that leads to the development of clinically relevant therapeutics could greatly improve infectious disease outcome.

RESULTS

Phenotypic screens for enrichment of B. thailandensis persister populations.To cast a wide net for discovery of genes and gene networks involved in bacterial persistence, we applied multiple methods to enrich for B. thailandensis persisters, including flow sorting for low proton motive force (PMF), cell lysis using the β-lactam antibiotic meropenem (Mpm), and culture aging on solid medium via nutrient starvation. First, we found that a subpopulation (∼0.5%) of mid-log-phase B. thailandensis intracellularly retained the fluorescent dye DiBAC4(3), signifying membrane depolarization (Fig. 1A). Bacteria normally maintain an electric potential gradient across their cytoplasmic membrane, which is essential for ATP synthesis. We further observed that the percentage of cells that retained DiBAC4(3) increased upon treatment with the amino acid analog serine hydroxamate (SHX; 10%), an inhibitor of cellular metabolism through the stringent response (31), and antibiotics such as ofloxacin (Ofl; 6%) and Mpm (38%) (Fig. 1A).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Characterization of persister phenotype. (A) Analysis of proton motive force (PMF) using fluorescence-activated cell sorting (FACS). Mid-log-growth-phase B. thailandensis cultures (OD600 of 0.2 to 0.3) were incubated in LB-Miller broth containing 50 μg/ml SHX, 20 μg/ml Ofl, or 20 μg/ml Mpm for 2 h at 37°C prior to addition of 1 μg/ml DiBAC4(3). Stationary cells (STAT) were collected from 14-day-old cultures on LB-Miller nutrient agar plates and incubated for 30 min at 37°C in LB-Miller broth containing 1 μg/ml DiBAC4(3). Flow cytometry analysis was performed on the BD FACSAria instrument using a fluorescein isothiocyanate (FITC) filter (Em: 530/30) after excitation with the 488-nm laser. Mid-log-phase cells not treated with DiBAC4(3) were used to position the signal threshold gate. One representative experiment out of five independent replicates is shown. (B) Metabolic index of B. thailandensis sorted for DiBAC4(3) retention. Mid-log-phase B. thailandensis cultures were treated with 1 μg/ml DiBAC4(3) and sorted into subpopulations with low FI (FITC-A, 500 to ≤1,000 fluorescence units [FU]) and high FI (FITC-A, 103 to ≤104 FU). Ten thousand cells were collected by FACS and combined with BacTiter-Glo reagent to quantify cellular ATP. Presented are the average and standard deviation values from three independent experiments. (C) Evaluation of bacterial tolerance to cytotoxic antibiotic concentrations. (Left) Mid-log-phase B. thailandensis cultures were treated with 1 μg/ml DiBAC4(3) and sorted to subpopulations with low FI (FITC-A, 500 to ≤1,000 FU) and high FI (FITC-A, 103 to ≤104 FU). Ten thousand cells were collected by FACS and grown to early log phase (OD600 of ∼0.2) prior to addition of 10 μg/ml Mpm. Persister populations surviving 24 h of antibiotic treatment were calculated as CFU/ml recovered on LB-Miller agar plates after antibiotic treatment relative to the CFU/ml prior to antibiotic exposure. Average and standard deviation values were derived from three independent experiments. (Middle) Mid-log-phase B. thailandensis cultures were subjected to 10 μg/ml Ofl or 10 μg/ml Mpm for 24 h. Percentage of surviving persisters was determined after 24 h of antibiotic exposure. Average and standard deviation values were derived from five independent experiments. (Right) A single CFU from 1-, 7-, 14-, and 21-day-old B. thailandensis cultures was inoculated into LB-Miller broth and cultured at 37°C until OD600 reached ∼0.2, followed by addition of 10 μg/ml Ofl for 24 h. Presented are the average and standard deviation values of percent bacterial survival calculated from 10 CFU inocula originating from three independent bacterial populations.

We next isolated exponentially dividing cells that exhibited high and low fluorescence intensity (FI) upon incubation with DiBAC4(3) and found that FI was inversely correlated with metabolic index, as measured by the total ATP pool (Fig. 1B). Upon expansion of sorted cells in nutrient-rich medium, we observed a >100-fold increase in the survival rate of cells that originated from populations sorted for high FI compared to low FI after 24-h treatment with twice the minimal bactericidal concentration (2× MBC) of Ofl (Fig. 1C, left panel). All together, these results suggest that persisters exhibit low PMF compared to wild-type cells and that we can leverage this property to sort persisters that display enhanced DiBAC4(3) uptake via increased fluorescence. A decrease in PMF has been previously linked to cell dormancy and multidrug tolerance phenotype in bacterial persistence (32–35).

Our second assay was based on enrichment for Burkholderia persisters through lysis of the actively dividing cells with Mpm, a broad-spectrum carbapenem antibiotic that inhibits cell wall synthesis. In a previous study, we had reported varied potencies of antibiotics against B. thailandensis persisters, suggesting that there are different mechanisms that lead to persister establishment (24). Treatment with Mpm has led to improved outcomes in melioidosis patients diagnosed with severe sepsis (36). We found that Mpm is less effective at killing Burkholderia persisters than the DNA gyrase inhibitor Ofl (Fig. 1C, middle panel). Furthermore, treatment with the MIC of Mpm significantly increased the population of B. thailandensis cells with low PMF (and high FI) within the span of one replication cycle, 2 h of exposure (Fig. 1A). For our transcriptome studies, we collected RNA from bacterial populations that survived 2× MIC of Mpm after more than 2 cycles of replication, with ∼80% enrichment of persister populations (see Fig. S1A in the supplemental material).

FIG S1

(A) Isolation of B. thailandensis populations with low PMF. Mid-log-phase cultures (OD600 of 0.2 to 0.3) of B. thailandensis were incubated in 1 μg/ml DiBAC4(3) for 30 min prior to cell sorting by FACS on the BD FACSAria using an FITC filter (Em: 530/30). Hierarchical gating was performed, first, on the cell distribution plot of the log scale of side scatter (SSC-A) versus forward scatter (FSC), followed by gating of log-scale fluorescence intensity at the 530/30 emission filter versus FSC cell distribution. Bacteria treated with the metabolic inhibitor SHX served as a positive control to set the gate for collection of persisters with high FI levels due to DiBAC4(3) dye retention. The negative control was B. thailandensis not treated with DiBAC4(3). The right panel of FACS histograms presents the percentage of cells with low PMF at 4 h of treatment with 2× MBC of Mpm, which was the point of transcriptome collection for next-generation sequencing (NGS). (B) Hierarchical clustering of genome-wide expression profiles (rows) and experimental conditions (columns) based on normalized counts per million reads (CPM). Control groups (C) clustered on the right side of the plot, and persister samples (P) clustered together on the left side of the plot. One of the three experimental samples for low PMF (PMF_006P) clustered with the respective control group (PMF_005C), indicating potential sample cross-contamination in this particular experiment. Additionally, one of the three experimental samples exposed to Mpm at 1× MBC (mpm_018P) experienced a stronger bacteriostatic effect, as evidenced by its clustering with the static populations of 21-day-old cultures. Download FIG S1, TIF file, 0.2 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

Finally, we subjected cells to nutrient starvation by long-term incubation on LB-Miller agar plates at 4°C, as a model for population aging and biofilm-induced resistance to antibiotics that could naturally occur in nutrient-limited soil microcosms. In this case, population aging is comparable to long-term stationary phase, bearing the complexity of respiratory and metabolic adaptations that are quite different from aging resulting from asymmetrical division. We observed a gradual increase in the survival rates of B. thailandensis persisters isolated from 7-, 14-, and 21-day-old cultures on agar plates (Fig. 1C, right panel). As a measure of tolerance, we used two doses of MBC of Ofl that normally killed ≥99.9% of the bacterial populations isolated from a 1-day culture on LB agar. At least ∼60% of the cells collected from 14-day-old cultures had low PMF (Fig. 1A, B.th. STAT). In summary, our transcriptome screens were based on (i) physical separation of persisters via cell sorting (low PMF), (ii) ∼80% persister enrichment following Mpm treatment (Mpm), and (iii) aged bacterial populations harboring ∼60% of cells with the persister phenotype (aged).

Genome-wide transcriptome analysis of persisters.We performed RNA Mi-Seq transcription profiling of B. thailandensis gene expression patterns (i) in DiBAC4(3)-treated mid-log-phase cultures sorted by high FI (low PMF), (ii) in cells lysed with 2× MBC of meropenem for 4 h (Mpm treated), and (iii) following 21 days of aging on nutrient agar at 8°C (aged). We analyzed relative changes in gene expression between persister-enriched and the following controls processed in parallel: (i) DiBAC4(3)-treated cells sorted for low FI, (ii) untreated exponential cultures, and (iii) cells following 1 day of growth on nutrient agar.

Hierarchical clustering analysis of the global transcriptomes revealed good data reproducibility between experimental replicates, where control samples from all tested conditions clustered together and showed distinct gene expression patterns from the persister-enriched samples (Fig. S1B). Applying principal-component analysis (PCA), we observed that the transcriptome of the aged population showed the most distinct differential gene expression (Fig. 2A). The persisters enriched via cell sorting for low PMF or exposure to Mpm displayed relatively convergent global gene expression profiles to their respective control groups. These results suggest that the mechanism governing resilience in the aged populations differs significantly from that of bacterial membrane depolarization and Mpm tolerance. We also note that the low-PMF controls (red dots in Fig. 2) clustered with the controls from the other phenotypic assays on the left side of the plot, signifying that DiBAC4(3) dye itself did not contribute to the persister phenotype through interference with bacterial membrane function.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Comparative analysis of genome-wide transcriptomics data from three persistence phenotypic screens. (A) Principal-component analysis of RNA-seq data was performed using DESeq2 as part of PiReT analyses on the EDGE bioinformatics platform. The control and persister populations for the three experimental methods are indicated in the legend. (B) Comparison of the number of differentially expressed genes (log2-fold) in persister populations enriched via three methods (low PMF, Mpm, and aged). The P value cutoffs applied for differentially regulated genes were ≤0.001 for the populations with low PMF, and P values of ≤0.01 were used for Mpm-treated and aged populations featuring higher variability. Colored ovals in the legend show five functional groups of interest with gene representatives differentially expressed in persister-enriched populations.

DESeq statistical analysis revealed that persisters with low PMF had the highest number (300) of differentially expressed genes (log2-fold), with a P value of ≤0.001 and low probability of false discoveries (false-discovery rate [FDR] < 0.01) (Fig. 2B and Table S1). Our selection criteria for differentially expressed genes was a log2 (xp/xc) value of ≥1 for gene overexpression in the persister cells (xp) versus the control population (xc) and a log2 (xp/xc) value of ≤−1 for gene repression. In comparison, the number of differentially expressed genes in aged (Table S2) and Mpm-treated (Table S3) populations was lower (41 and 146, respectively), even at less stringent statistical parameters, FDR < 0.5 and P value ≤ 0.01 (Fig. 2B and Tables S2 and S3). Furthermore, the low-PMF cells expressed a much smaller number of downregulated genes (24) than the upregulated gene pool (276), compared to the differential gene expression profiles of the aged and Mpm-treated persisters and previous genome-wide transcriptomics studies of persistence (5, 19).

TABLE S1

Differentially expressed genes in the low-PMF screen. Download Table S1, XLSX file, 0.04 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

TABLE S2

Differentially expressed genes in the aged screen. Download Table S2, XLSX file, 0.02 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

TABLE S3

Differentially expressed gene in the Mpm screen. Download Table S3, XLSX file, 0.01 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

Genes encoding polyketide synthases (PKSs) and nonribosomal peptide synthetases (NRPSs) were upregulated in persisters from all three screens (Fig. 2B and Fig. 3A). PKSs polymerize simple fatty acids into a large variety of secondary metabolite products called polyketides. Due to the complexity of polyketide synthesis, PKS genes are organized into operon systems (Fig. S2A). There are 11 gene clusters with predicted PKS and NRPS functions in the Burkholderia genome (29). We found PKS genes from two operons, including BTH_I2366 and BTH_II1677, to be upregulated in low-PMF and Mpm-treated persisters, whereas aged persisters showed upregulated PKS genes located only on chromosome I (Fig. 3A). PKS operons often include NRPS genes that work in cooperation with PKSs to synthesize polyketide-peptide hybrid macromolecules. An NRPS gene, BTH_I1956, located in the PKS operon on chromosome I was upregulated in both low-PMF and aged persisters (Fig. 3A).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Functional pathways differentially regulated in persisters enriched via different experimental methods. Heatmaps show the fold change (FC) in log2 scale of genes in selected pathways between persisters and control groups. Gene names are labeled using locus tags. FC values were applied in the hierarchical clustering of genes (rows) and experimental conditions (columns).

FIG S2

Activation of PKS/NRPS operon on B. thailandensis chromosome I upon Mpm treatment contributes to increased antibiotic tolerance. (A) RNA-seq data sets were analyzed using the JBrowse tools for data visualization (https://jbrowse.org/). The operon encompassing BTH_I2366 is activated in bacterial populations exposed to 1× MBC of Mpm. Reads mapping to the BTH_I2363–I2367 megaenzyme system involved in the synthesis of natural products were ∼5- to 10-fold higher in antibiotic-stressed cells than in the untreated control sample. PKS, polyketide synthetase; NRPS, nonribosomal peptide synthase; DAAS, dihydroaeruginoic acid synthetase. (B) BTH_I2366 loss-of-function phenotype correlated with decreased resistance to various antibiotics including Ofl. The left panel shows a significant decrease in the bacterial counts (CFU) of BTH_I2366 mutant strain BT08684 compared to wild-type B. thailandensis E264 strain upon plating of limiting dilutions of cells exposed to antibiotic treatment. Sequential passaging of the BT08684 cells on LB-Miller agar plates resulted in spontaneous loss of the Tn23 insertion sequence with the BTH_I2366 locus, and bacterial cells reverted to a wild-type phenotype with recovered resilience to antibiotic treatment. The right panel shows PCR products validating Tn23 insertion 1 kb downstream of the BTH_I2366 start codon. Genomic DNA isolated from wild-type B. thailandensis E264 and BT086864 mutant cells was used as the template in PCR with two primer pairs, ORFfw0.5 and LacZrev, to amplify a 1-kb region between 0.5 kb and 1.5 kb downstream of the start codon (lane 2). In the BTH08684 mutant strain, insertion of Tn23 1 kb downstream of the BTH_I2366 start codon generated an ∼600-bp PCR product (lane 3) and disrupted the formation of the 1-kb PCR product (lane 4 versus lane 2). Upon reverse mutation, the 600-bp product was mostly lost (lane 5) and the 1-kb product was restored (lane 6). Note that the abundant PCR products in lanes 1 and 5 at ∼3 kb resulted from nonspecific primer annealing. The 1-kb Plus DNA ladder (lane M; New England Biolabs) was used as the standard. Download FIG S2, TIF file, 1.7 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

Four genes that encode various classes of Clp proteases (BTH_II1046, BTH_II1895, BTH_II0140, and BTH_II1476) were primarily upregulated in persisters with low PMF (Fig. 3B). The ATP-binding subunit of type B Clp protease (BTH_II1895) was expressed at higher levels in all three screens. Clp proteases have been previously linked to persistence through proteolytic regulation of toxin-antitoxin systems (5, 37) and participate in bacterial stress response mechanisms via the degradation of misfolded proteins (38, 39) and activation of outer membrane biosynthesis enzymes (40).

The three screens yielded other upregulated genes that likely play a role in regulating bacterial persistence. We observed overexpression of mobile elements in bacteriophage operon genes from primarily the low-PMF screen (Fig. 3C and Table S1). Bacteriophage genes have previously been associated with biofilm formation and antibiotic tolerance in Pseudomonas aeruginosa (41, 42). Genes regulating lipid biosynthesis were identified in all three screens (Fig. 2B, Fig. 3D, and Tables S1 to S3). We detected upregulation of multiple gene loci involved in fatty acid catabolism in persisters, including acyl coenzyme A (acyl-CoA) dehydrogenase (ACAD), which catalyzes fatty acid β-oxidation (BTH_I0060 and BTH_I3298) ,and long-chain acyl-CoA ligase, which catalyzes the breakdown of complex fatty acids (BTH_II2096 and BTH_I2529) (43–45). ATP-binding cassette (ABC) transporter systems for amino acids, such as the arginine/ornithine transporter (BTH_I2383), were upregulated in persisters from the low-PMF and aged screens (Fig. 2B, Fig. 3E, and Tables S1 and S2).

We also noted other interesting gene expression patterns among the three screens. N-acetylmuramoyl-l-alanine amidases (NAMLAAs), which function in bacterial cell wall remodeling, were upregulated in cells with low PMF (BTH_I0722) and Mpm treatment (BTH_II0537) but remained silenced in aged cells (Tables S1 to S3). Finally, the gene encoding MarR, multiple antibiotic resistance transcription factor (BTH_I2391), was upregulated in aged cultures but not in cells exposed to Mpm (Tables S2 and S3).

Validation of RNA-seq differential gene expression using RT-qPCR.To validate the RNA sequencing (RNA-seq) results, we performed targeted RT-qPCR using culture aging (for 7, 14, and 21 days) and antibiotic treatment (Mpm, Ofl, and ceftazidime (Cef) at 1× MBC) as stress factors for the enrichment of persister cells (Table 1). We observed consistent upregulation of the PKS gene BTH_I2366 in response to all tested aging and antibiotic stress conditions. NRPS (BTH_I1956) expression was upregulated in response to the three antibiotics but not in the aged cultures. Upregulation of PKS and NRPS transcripts was the only common gene expression pattern seen by RT-qPCR across bacteria treated with each of the three antibiotics. Genes encoding ClpP (BTH_II1046) and ClpA/B (BTH_II0140) were upregulated in persisters from both aging and antibiotic stress cultures. However, we noted that Ofl treatment inhibited clp expression compared to the β-lactam antibiotics Mpm and Cef. We observed a similar divergent effect of Ofl on the cell wall-remodeling enzyme NAMLAA (BTH_II0537), which was downregulated in Burkholderia treated with Ofl but upregulated upon Mpm and Cef treatment. The gene encoding HK97 phage major capsid protein (BTH_II1047) was upregulated in the aged cultures but not in response to antibiotics. The RT-qPCR results also confirmed the RNA-seq data showing that MarR (BTH_I2391) was upregulated in the aged (14- and 21-day) cultures. Additionally, MarR transcript levels were increased in response to Cef but remained the same with Mpm and Ofl treatment, indicating divergent cell responses to different antibiotics and mechanisms of action.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

RT-qPCR validation of select genes identified by RNA-seq to be differentially expressed in persistersa

Polyketide synthase gene function augments bacterial resilience to antibiotics targeting DNA replication.Upon identification of the PKS gene BTH_I2366 as a general persistence biomarker in B. thailandensis (Fig. 3A), we further investigated its role in bacterial survival. We tested the antibiotic sensitivity of three transposon mutants (BT02489, BT06658, and BT08684), generated by random insertion of Tn23 transposon sequences within the BTH_I2366 gene locus (46). Treatment with quinolone antibiotics Ofl and ciprofloxacin (Cipro), which inhibit DNA replication, led to a significant decrease in the number of persisters in the mutants compared to the wild type (WT), indicating that transposon mutagenesis of BTH_I2366 increased the susceptibility of persister cells to antibiotic killing (Fig. 4A). Among the mutant strains, BT08684 exhibited the highest variability in persister cell survival upon antibiotic treatment. This observation prompted the study of possible reverse BT08684 mutations, resulting from spontaneous Tn23 excisions. PCR analysis of genomic DNA from BT08684 cell populations that exhibit wild-type sensitivity to Ofl displayed loss of the transposon insertion in the BTH_I2366 gene (Fig. S2B), thus validating the role of PKS gene function in antibiotic tolerance.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

PKS gene activity is a determinant of B. thailandensis persistence. (A) PKS BTH_I2366 loss of function increased B. thailandensis sensitivity to antibiotics that target DNA replication. Mid-log-growth bacterial cultures (OD600 of ∼0.2 to 0.3) were treated with 1× MBC of Mpm (10 μg/ml), Ofl (10 μg/ml), and Cipro (5 μg/ml) for 24 h. The percentages of persister populations were determined as the ratio of antibiotic-surviving CFU/ml versus the cell number (CFU/ml) prior to antibiotic exposure. Presented are data from four independent antibiotic treatment experiments. The significance of the registered sensitivity to antibiotic treatment was evaluated using one-way ANOVA, where # indicates a P value of <0.002. (B) The PKS-encoding BTH_I2366 gene is responsible for membrane depolarization upon Mpm treatment. Wild-type B. thailandensis and mutant strains were grown to mid-log phase (OD600 of 0.2 to 0.3) and were treated with 10 μg/ml Mpm for 3 h. Bacterial cultures were incubated in 1 μg/ml DiBAC4(3) for 30 min prior to flow cytometry analysis on the BD FACSAria using an FITC filter (Em: 530/30). Gate P2 (top panel) represents the population of cells with low PMF that retain the fluorescent dye, shown in green. The total analyzed bacterial population (P1) is shown in red. Data are representative of five independent experiments.

Complementation of the BTH_I2366 mutant background is challenging due to the broad-spectrum antibiotic resistance of the mutant strains, both natural (resulting from β-lactamase activity or efflux pump-mediated exclusion of aminoglycosides) and acquired (Tn23-associated trimethoprim resistance [Tpr]), that limited our ability to use selective pressure for transgene delivery. Thus, we generated a B. thailandensis strain that expresses green fluorescent protein (GFP) with a Tpr selection marker incorporated 25 bp downstream of the glucosamine-6-phosphate synthetase (GlmS) stop codon using the mini-Tn7 system. This B. thailandensis strain, designated B. thailandensis TpR-Tn7, served as a persistence control to evaluate the impact of transposon-acquired antibiotic resistance from a mobile element. We did not observe a significant difference in the survival rates of B. thailandensis wild type (WT) and the B. thailandensis TpR-Tn7 strain (Tn7-G) when treated with quinolone antibiotics (Fig. 4A).

Interestingly, treatment with Mpm led to a different outcome from that with the antibiotics targeting DNA replication (Ofl and Cipro). First, we observed no difference between the survival rates of persisters in wild-type B. thailandensis and in BTH_I2366 mutants treated with Mpm (Fig. 4A). Second, there was a marked decrease in the number of DiBAC4(3)-stained cells (∼1 to 2%) for BTH_I2366 mutants BT02489 and BT06658 compared to wild type (20%) in response to Mpm, indicating that BTH_I2366 mutagenesis counteracted the membrane depolarization induced by an antibiotic that targets cell wall synthesis (Fig. 4B). Treatment with Ofl slightly increased the population of persisters with low PMF, at ∼6 to 8% compared to 40% with Mpm exposure (data not shown). However, we did not observe any significant difference between DiBAC4(3) staining in the wild type and in the PKS mutants exposed to quinolone antibiotics. Collectively, these data indicate that PKS function contributes broadly to the establishment of a persistence phenotype in B. thailandensis exposed to antibiotics.

Natural products with HDAC inhibitor properties regulate DNA synthesis in B. thailandensis.A previous study had identified that BTH_I2366 is part of a megaenzyme PKS/NRPS gene cluster complex that produces natural products in the form of acyldepsipeptides with histone deacetylase (HDAC) inhibitor properties, termed burkholdacs (29). In that study, a galactose-regulated promoter (Ptac) was engineered to drive PKS/NRPS overexpression and the resultant natural products were synthetically modified via dithiothreitol reduction to increase their potency as HDAC inhibitors. Here, we observed that the conditioned medium of wild-type cells treated with Mpm for at least 12 h exhibited ∼70% HDAC inhibition, indicating release of secondary metabolites with HDAC-inhibitory properties via Mpm-induced cell lysis (Fig. 5A). HDAC inhibition in conditioned medium was reduced to ∼30 to 40% in Mpm-treated BTH_I2366 mutants, indicating that BTH_I2366 function is required for efficient production of metabolites with HDAC-inhibitory properties in response to antibiotic stress.

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

The PKS operon encodes secondary metabolites with HDAC-inhibitory properties that regulate DNA synthesis in B. thailandensis. (A) Treatment with Mpm stimulated accumulation of HDAC inhibitors in wild-type B. thailandensis cultures. Mpm at 10 μg/ml was added to cultures of wild-type B. thailandensis E264 and BTH_I2366 mutants grown to an OD600 of ∼0.2. Crude conditioned medium from overnight (18-h) cultures grown without agitation in a 96-well plate was analyzed for inhibition of HDAC activity. Trichostatin A (TSA; 1 μM) was used as an inhibitor control. The percent inhibition of HDAC activity (%) was calculated relative to samples with no drug or conditioned medium. Shown are the average and standard deviation values from three independent experiments. (B) Loss of PKS gene function results in accumulation of DNA in stationary-phase B. thailandensis cultures. Flow cytometry analysis was performed on mid-log-phase (LOG) and stationary-phase (STAT) cultures to determine the DNA content in wild-type B. thailandensis and BTH_I2366 mutant cells. The membrane-permeant Vybrant DyeCycle green stain was added to bacterial cultures 30 min prior to analysis on a BD FACSAria instrument. DNA content was measured as mean fluorescence intensity (MFI) units. Shown are data from 10,000 cells analyzed in one representative experiment out of three independent replicates.

Given the key role of HDAC inhibitors in stimulating DNA replication in both eukaryotes (47) and bacteria (48), we measured DNA content in both the wild type and BTH_I2366 mutants by flow cytometry. DNA content was decreased by ∼50% in the BTH_I2366 mutants compared to wild-type cells in the mid-log phase, consistent with reduced expression of HDAC inhibitors in the mutants (Fig. 5B). Interestingly, DNA content was comparatively increased in the BTH_I2366 mutants in stationary phase, suggesting that secondary metabolites with HDAC-inhibitory properties may negatively regulate DNA accumulation in stationary-phase bacteria. We found no difference in the rate of cell growth between the wild-type and mutant strains in nutrient-rich medium (data not shown).

The fatty acid synthesis inhibitor CP-640186 enhanced antibiotic activity against persister populations of various Gram-negative bacteria.Polyketide synthase function is closely linked to fatty acid synthesis, and the two pathways share common components, including acetyl-CoA and malonyl-CoA (49). We tested the effect of CP-640186, an inhibitor of mammalian acetyl-CoA carboxylase (mACC1/2), which catalyzes the synthesis of malonyl-CoA, on the survival of B. thailandensis persisters. Since the compound exhibited no antibacterial activity on its own (data not shown), we combined CP-640186 with Mpm treatment and observed an ∼10-fold decline in persister survival in wild-type B. thailandensis (Fig. 6A). Interestingly, we registered no additive effect of CP-640186 on Ofl potency against the persister population.

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Acetyl-CoA carboxylase inhibits persister formation in Gram-negative bacteria. (A) Combination of Mpm with acetyl-CoA carboxylase inhibitor CP-640186 reduced persisters in B. thailandensis. Mid-log-phase cultures (OD600 of 0.2 to 0.3) of B. thailandensis were treated with 10 μg/ml Mpm, 10 μg/ml Ofl, or a combination of antibiotic and 50 μg/ml CP-640186 for 24 h. The percentage of persisters surviving antibiotic treatment was calculated as previously described. Shown are average and standard deviation values derived from four independent experiments. (B) CP-640186 potentiated Mpm activity against persisters in three Gram-negative bacterial pathogens. Mpm at 2× MBC alone or in combination with 50 μg/ml CP-640186 was added to each bacterial species (OD600 of ∼0.2). Persister populations were evaluated after 24 h of antibiotic exposure as described above. Shown are average and standard deviation values derived from four independent experiments.

We further tested the effect of CP-640186 on Mpm potency against persisters from other Gram-negative pathogens, including Pseudomonas aeruginosa, Salmonella enterica serovar Typhimurium, and an attenuated strain of Yersinia pestis A1122. Upon addition of CP-640186 to Mpm treatment, we observed ∼3-, 4-, and 8-fold decreases in the percentage of P. aeruginosa, S. Typhimurium, and Y. pestis persister cells, respectively, compared to antibiotic monotherapy (Fig. 6B). Altogether, our results indicate that fatty acid metabolism is a viable target for development of antimicrobials to combat bacterial persistence.

DISCUSSION

Bacterial persistence is a transient survival mechanism that mediates noninherited antibiotic resistance and contributes to extended latency in chronic disease (2). Persisters comprise a small fraction (0.0001% to 1%) of isogenic bacterial populations and exhibit highly dynamic phenotypic changes in response to a variety of environmental stress conditions (52–54). Our functional genomics studies indicate that nutrient starvation and antibiotic toxicity stimulate gene clusters that drive production of secondary metabolites with HDAC-inhibitory properties. These natural products contribute to the establishment of a persistence phenotype manifested in bacterial membrane depolarization, low metabolic index, and bacterial resilience toward antibiotics targeting DNA and protein synthesis. A statistical comparison of the global transcriptome profiles showed that persisters enriched by the three methods generated very different gene expression patterns, underscoring the multimechanistic and complex nature of bacterial persistence and antibiotic tolerance that can shape cellular phenotypic behaviors. However, there were also pathways activated in common.

We found upregulation of multiple gene loci that encode PKSs, in particular BTH_I2366, in B. thailandensis persisters enriched by all three phenotypic screens. PKSs form large modules of megaenzymes that synthesize a variety of natural products and complex lipids needed for cell growth and metabolism (55). Previous studies have characterized B. thailandensis natural products produced by the gene cluster that encompasses BTH_I2366 as histone deacetylase inhibitors (termed burkholdacs) (29). One possible role for the burkholdacs in bacterial cells is regulation of posttranscriptional modification (PTM) of proteins upon onset of stressful conditions to allow for quick adaptations. PTM through lysine acetylation is an evolutionarily conserved strategy for regulation of enzyme activity, DNA-protein interactions, and protein stability (48, 56, 57). Although the role of burkholdacs in global bacterial metabolism remains to be determined, our study underscores their contribution to bacterial resilience to antibiotics.

It remains an open question whether pathogenic Burkholderia species synthesize HDAC inhibitors in response to stress, including antibiotic exposure, to counteract the immune defenses of the eukaryotic host. Treatment of B. thailandensis with the antibiotic trimethoprim has been previously shown to induce synthesis of a variety of secondary metabolites (many of them PKS and NRPS products) that act as quorum-sensing molecules or antibiotics (58). Burkholderia species, including the highly pathogenic select agent B. pseudomallei, are known for their large genomes that contain multiple pathogenic islands and large numbers of natural product biosynthetic gene clusters (59). For example, gene cluster BTH_II0204–207/BPSS0130–133, which is highly conserved between B. thailandensis E264 and B. pseudomallei K96243, was found to produce chemical analogs of betulinan and terferol, which are potent inhibitors of eukaryotic phosphodiesterase (PDE) (60). As a catalyst for cAMP and cGMP hydrolysis, PDE has been implicated in host inflammatory responses and represents a therapeutic target for inflammatory disease (63, 64). PKS-encoding gene islands are not specific to the Burkholderia genus and have been previously associated with pathogenicity in the gut microflora. A genotoxin encoded by the pks pathogenicity island in Escherichia coli was found to induce DNA double-strand breaks in host mammalian cells, resulting in higher mutation rates (61, 62).

In addition to PKSs, we also identified genes that regulate lipid catabolism in persisters. PKSs are analogous in architecture to, and share sequence homology with, type I fatty acid synthases, such as FAS-I, responsible for de novo fatty acid synthesis (49, 65). The two enzyme assembly lines share homologous domains (acyl carrier protein [ACP]) that carry the growing polyketide or fatty acid chains and utilize acyl-CoA and malonyl-CoA as precursors (49). We observed that Mpm treatment in combination with the inhibitor of malonyl-CoA synthesis (CP-640186) resulted in a decrease of persisters in multiple pathogens, including B. thailandensis, P. aeruginosa, S. Typhimurium, and attenuated Y. pestis. Our findings suggest that targeting acetyl-CoA carboxylase enzymatic function could be a viable therapeutic strategy for the elimination of persister populations in Gram-negative bacteria.

Fatty acids are a major component of the bacterial cell membrane and an important source of metabolic energy. Intermediates of the saturated fatty acid biosynthetic pathway provide short-chain 3-hydroxy fatty acids that are needed for synthesis of biotin, the lipid A component of Gram-negative bacteria outer membranes, and acyl-homoserine lactones (AHLs) of quorum-sensing molecules (66–68). The by-products of fatty acid degradation can fuel pathways for amino acid synthesis and energy production or serve as building blocks for a variety of secondary metabolites. We detected upregulation of multiple gene loci involved in fatty acid catabolism in persisters enriched by at least two of the three phenotypic screens, including acyl-CoA dehydrogenase (ACAD), which catalyzes fatty acid β-oxidation (BTH_I0060 and BTH_I3298), and long-chain acyl-CoA ligase, which catalyzes the breakdown of complex fatty acids (BTH_II2096 and BTH_I2529) (43–45). Other groups have reported activation of bifunctional enoyl-CoA hydratase (BTH_II0417) upon exposure of B. thailandensis DW503 to 10× MIC of the fluoroquinolone antibiotic ciprofloxacin and the DNA-damaging chemical mitomycin C (8, 69). Enoyl-CoA hydratase is an essential enzyme in fatty acid β-oxidation, which catalyzes the complete oxidation of long-chain fatty acids to provide an energy source and building blocks for primary and secondary metabolites.

Our transcriptomics studies found that ∼10% of genes differentially upregulated in persister populations encoded ABC transporters (Fig. 3E; see also Tables S1 to S3 in the supplemental material). ABC transporter systems are essential for the import and export of peptides, proteins, polysaccharides, and many other molecules to maintain metabolic homeostasis in the cell (70). All ABC transporter systems share a highly conserved ATP-hydrolyzing domain that couples ATP binding/hydrolysis to the translocation of essential nutrients or expulsion of antibiotics from the periplasm, thus contributing to drug resistance (76). Various gene loci encoding branched-chain amino acid ABC transporters were upregulated in both the aged and low-PMF persister populations (Fig. 3E and Tables S1 and 2). Persisters with low PMF had 38 upregulated ABC transporter genes, including sugar and amino acid transporters (Table S1). Branched-chain amino acid ABC transporters and peptide ABC transporters have previously been reported to be significantly upregulated in stationary-phase Listeria monocytogenes and B. thailandensis DW503 treated with 10× MIC of ciprofloxacin (8, 77). Additionally, other genes encoding ABC transporters were upregulated in B. thailandensis with low PMF (BTH_II1288 and BTH_I0880, Table S1) and in aged populations of L. monocytogenes (77).

Several independent transcriptomics studies have previously implicated bacteriophage and mobile element activation as universal stress response mechanisms that may contribute to antibiotic tolerance. For example, ciprofloxacin treatment of B. thailandensis at 10× MIC was found to stimulate bacteriophage encoding genomic islands and SOS response genes, and sublethal concentrations of β-lactam antibiotics activated phage development in various bacterial species, including Burkholderia cepacia (8, 78). Interestingly, we observed upregulation of phage operons in B. thailandensis persisters enriched by low PMF and aging but not in response to Mpm (Fig. 3C and Tables S1 to S3). These results demonstrate that species-specific gene regulatory mechanisms might drive divergent cellular responses to various antibiotic classes. Activation of bacteriophage operons in response to environmental stress can shape microbial evolution and accelerate genetic adaptations leading to antibiotic resistance (79).

Collectively, our results support a model in which bacterial persisters exhibit an outwardly dormant physiology but sustain active metabolic processes that are required to maintain persistence. We found that antibiotic treatment activated a Burkholderia gene cluster (BTH_I2363–2369) that produces natural products with known HDAC-inhibitory effects on host immune responses. Our study also sheds light on the role of fatty acid catabolism in antibiotic tolerance, identifying PKS systems as biomarkers of persistence. The modular architecture of the PKS megaenzymes provides a highly adaptable platform for the generation of polyketides with diverse structures and biological activities. Further studies will be needed to understand the mechanistic link between the antibiotic-induced production of natural products with HDAC-inhibitory properties and the establishment of persistent infections through suppression of host immune defenses. We expect that a more comprehensive understanding of the persistence state will enhance our ability to develop broad-spectrum, high-efficiency inhibitors to target persister populations.

MATERIALS AND METHODS

Bacterial strains, growth conditions, and reagents.B. thailandensis E264, P. aeruginosa Migula (ATCC 25347), S. Typhimurium (Loeffler) Chi4064 (ATCC 53648), and attenuated Y. pestis A1122 bacterial cultures were maintained on LB-Miller agar plates at room temperature. To test antibiotic sensitivity and metabolic activity, cultures initiated from single colonies were grown overnight in 5 ml LB-Miller broth (ThermoFisher Scientific) at 37°C with agitation at 120 rpm. To make mid-log cultures, overnight stationary cultures were diluted 500-fold in fresh LB-Miller broth and incubated at 37°C with agitation at 120 rpm for 4 to 5 h until culture density reached an optical density at 600 nm (OD600) of 0.2 to 0.3. Antibiotics were acquired from Sigma-Aldrich (St. Louis, MO): M2574 (Mpm), C0690500 (Cef), 1478108 (Ofl), 17850 (Cipro), and PZ0362 (CP-640186). The minimal bactericidal concentration (MBC) of antibiotics was determined for each new batch prior to performing persister assays according to standard CLSI protocols for MIC. MIC was determined as the lowest antibiotic concentration at which the OD600 remained <0.5 during 24 h of bacterial culture with agitation. MBC was calculated as the antibiotic concentration at which ∼90% loss of population viability was achieved.

Metabolic index determination.We measured the total ATP pool in bacterial populations using the BacTiter-Glo microbial cell viability assay (Promega, Madison. WI). We added equal volumes (100 μl) of BacTiter-Glo reagent and bacterial culture (100 μl) each and measured luminescence in white 96-well plates (Corning) using the BioTek Synergy HTX microplate reader. The cell number for the corresponding bacterial populations was derived from the OD600 measurements correlated with CFU/ml. Cell number was derived from a standard curve reflecting the linear correlation between OD600 (ranging from 0.1 to 0.5) and the CFU/ml. Metabolic index was calculated as the ratio between the ATP pool enumerated as relative light units (RLU) and cell number.

Enumeration of persister populations.We measured persister populations as a percentage of surviving cells after 24 h of antibiotic treatment (at 1× MBC or higher dose) relative to the number of cells (CFU/ml) prior to antibiotic exposure. Limiting dilutions (log10) of bacterial cultures were plated before and after antibiotic treatment on antibiotic-free LB-Miller agar plates and incubated at 37°C for 24 h (P. aeruginosa and S. Typhimurium) or 48 h (B. thailandensis) or up to 5 days (Y. pestis A1122). In parallel, bacterial cultures were plated on LB-Miller agar plates containing antibiotic at 1× MIC to evaluate the occurrence of genetic antibiotic-resistant mutants.

Enrichment of persister populations and RNA extraction.We applied three different methods to enrich for B. thailandensis persisters for RNA-seq transcriptomic analysis.

(i) Low PMF. We treated mid-log-phase B. thailandensis cultures (OD600 of 0.2 to 0.3) with 1 μg/ml DiBAC4(3) for 30 min, washed bacteria twice in dye-free LB-Miller broth, and resuspended bacteria in fresh LB broth to a density of OD600 of ∼0.1. Persisters were sorted using a blue laser [excitation (Ex)/emission (Em), 488/530(30) nm, respectively] with the BD FACSAria instrument (Becton Dickinson, NJ) (see Fig. S1A in the supplemental material). Total RNA was isolated from 106 cells sorted for either high or low FI. Collection of 106 cells with high FI was a lengthy process (about 4 h) given the low frequency of this population (0.5 to 1%). To restrict changes in the global transcriptome during sorting, cells were stained, sorted, washed, and resuspended in QIAzol lysis reagent (Qiagen, Germantown, MD) in batches within 15-min intervals.

(ii) Mpm-treated. We exposed B. thailandensis to 2× MBC of Mpm at the mid-log phase (OD600 of ∼0.3) for 4 h based on our previous observation that the metabolic index increased rapidly within 1 to 3 h of antibiotic exposure followed by a sudden drop (24). Furthermore, we had detected ∼80% absorption of the fluorescent dye DiBAC4(3) at this time point, which was indicative of high enrichment for persisters (Fig. S1A). Untreated B. thailandensis populations were used for the control samples. RNA was isolated from three independent bacterial cultures per experiment.

(iii) Aged. RNA was isolated from B. thailandensis cultured on LB-Miller agar plates for 21 days at 8 to 10°C. In addition to temperature stress and nutrient limitation, the low temperature allowed for the retention of moisture within the agar medium, which was essential for cell viability. Approximately 50 CFU was scraped from the agar cultures, resuspended in 1 ml phosphate-buffered saline (PBS), and incubated with 1,000 U lysozyme for 15 min at 37°C prior to lysing the cells in QIAzol reagent. Total RNA was isolated using the miRNeasy minikit (Qiagen) (24). Total RNA for the control samples was harvested from 1-day-old cultures and processed in parallel to the aged populations. Three independent starter cultures were used for sample preparation.

RNA-seq procedure and analysis. We utilized the RNA 6000 Pico kit (Agilent, catalog no. 5067-1513) to evaluate quality of RNA isolated from the bacterial cultures. Total RNA was depleted of rRNA using the Ribo-Zero Gram-negative bacterium kit (Illumina, catalog no. MRZGN126), and the concentration was determined with the Qubit RNA high-sensitivity (HS) assay (ThermoFisher Scientific, catalog no. Q32855). cDNA libraries were generated with the ScriptSeq v2 library preparation kit (Illumina, catalog no. SSV21106), titrated with the Qubit double-stranded DNA (dsDNA) HS assay (ThermoFisher Scientific, catalog no. Q32854), and sequenced on the Illumina Mi-Seq platform, generating paired-end 151-bp reads (25 to 40 million reads per sample). We used the high-sensitivity DNA kit (Agilent, catalog no. 5067-4626) to determine the average size of the library and normalized the input of sequencing libraries to 5 pmol. For RNA-seq analysis, we applied PiReT (https://github.com/mshakya/piret) implementation on the EDGE Bioinformatics platform (71). Briefly, PiReT performed quality control on raw RNA-seq data using FaQC v2.09 (72) and then mapped quality-controlled reads to the reference genome B. thailandensis E264 (GCA_000012365.1) using HISAT2 v2.1.0 (73). We counted the number of reads that mapped to different features including genes and coding sequences (CDS) using featureCounts v1.6.2 (74). We applied count data to calculate fragments per kilobase per million (FPKM) and found differentially expressed genes using R package DESeq2 v1.18.1 (75). We applied hierarchical clustering analysis as a quality control measure of the transcriptomic data and observed distinct differentiation between control groups and persisters from the three assays (Fig. S1B). Heatmaps of fold change (FC) were constructed using R package heatmap v1.0.10 wherein both rows and columns were clustered based on FC values.

Gene transcription validation via RT-qPCR.The Power SYBR green RNA-to-Ct one-step kit (ThermoFisher Scientific) was applied to perform qPCR analysis of bacterial transcripts on the ABI 7500 Fast System platform (Applied Biosystems). Total RNA was isolated from bacterial cultures, RNA concentration was determined using an ND1000 NanoDrop spectrophotometer (Thermo Scientific), and RNA quality was evaluated on a denaturing 1% agarose gel. Each RT-qPCR mixture (25 μl) was run in duplicate with transcript-specific forward and reverse primers (Table S4) at a final concentration of 600 pmol. Temperature melting curves of the PCR products showed a single amplicon for each primer pair. Purified PCR products were sequence verified for specificity. MinE transcripts served as an internal sample control since their relative abundance was comparable to the average transcript levels of the genes investigated in this study. rRNA (rrsD/16S) was used as a sample quality control and a reference gene transcript for evaluation of MinE fluctuations between experimental conditions. The RNA level of each gene of interest was calculated as a ratio relative to the gene transcript level in the control sample (log2), applying the ΔΔCT threshold cycle (CT) method.

TABLE S4

Forward and reverse primers for qRT-PCR and PCR validation of Tn23 insertion into BTH_I2366 locus. Download Table S4, DOCX file, 0.01 MB.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

Histone deacetylase assay.We applied the histone deacetylase (HDAC) activity assay (fluorometric) kit (Abcam, Cambridge, MA). The enzymatic reaction mixtures (50 μl) were prepared according to the manufacturer’s two-step protocol for quantification of HDAC activity. Trichostatin A (TSA; 1 μM) was used as the inhibitor control in the assay. Mid-log-phase B. thailandensis cultures were treated with 10 μg/ml Mpm and incubated for 12 to 24 h at 37°C in a 96-well plate (150 μl per well). Crude conditioned medium (30 μl) was then assayed for HDAC-inhibitory properties. The enzymatic reactions were conducted in duplicate at room temperature (RT) for 30 min prior to the addition of stop solution and developer. Fluorescence intensity was measured using a microplate reader (Biotec) at Ex/Em 350(30)/440(20) nm. One hundred percent HDAC activity was defined as fluorescence intensity (F100) in the absence of drug or B. thailandensis crude conditioned medium. Zero percent activity was defined as fluorescence intensity for each data set without HDAC nuclear extract (F0). Percent HDAC inhibition was determined according to the equation % inhibition = [1 − (FI − F0)/(F100 − F0)], where FI is the fluorescence intensity in the presence of TSA or B. thailandensis conditioned medium.

Analysis of in vivo DNA content in bacterial cells.Vybrant DyeCycle green stain (ThermoFisher Scientific) was used to analyze DNA content in logarithmically expanding (LOG) and stationary-phase (STAT) B. thailandensis cultures. DyeCycle green stain was added at a final concentration of 10 μM to 108 cells resuspended in 1 ml LB-Miller broth. Bacterial cells were incubated in the stain-containing medium for 30 min at 37°C and were analyzed by flow cytometry (BD FACSAria II) using an argon laser. Unstained cells were used as negative controls to position the gates for analysis of fluorescence intensity. The cell distribution plots of log scale of side scatter (SSC-A) versus forward scatter (FSC) were used to gate the total population of analyzed cells, followed by gating of log-scale fluorescence intensity at the 530/30 emission filter versus FSC cell distribution. Mean fluorescence intensity (MFI) and standard deviation data were obtained through BD FACSDiva software (Becton, Dickinson, NJ). MFI was used as a measure of DNA content in live bacterial cells.

Generation of B. thailandensis TpR-Tn7 strain.Insertion of the mini-Tn7T-Tp-GFP sequence into B. thailandensis E264 cells was facilitated via triparental mating of B. thailandensis E264 with E. coli BW29427 strains harboring pUC18T-mini-Tn7T-Tp-gfpmut3 and pTNS2, Addgene plasmids 65033 and #64968, respectively (50). E. coli BW29427 is a diaminopimelic acid (DAP) auxotrophic mutant that provides convenient and efficient conjugation-mediated plasmid transfer followed by counterselection of the donor strain on DAP-deficient medium. The pTNS2 and pUC18T-mini-Tn7T-Tp-gfpmut3 plasmids were independently transfected in E. coli BW29427 cells by heat shock and were maintained on ampicillin (Amp) or Tp medium, respectively, supplemented with 0.3 mM DAP. The pUC18T-mini-Tn7T-Tp-gfpmut3 plasmid was mobilized into B. thailandensis E264 cells via triparental mating with the helper pTNS2 (carrying the transposase gene) at a 1:1:1 ratio of bacterial cells (106 in 100 μl each), collected on a 20-mm cellulose filter, and incubated on LB-Amp-DAP overnight at 37°C. Single clones of B. thailandensis transformants were isolated on LB agar supplemented with 100 μg/ml Tp, and the insertion of Tn7T-Tp-GFP sequences was confirmed via colony PCR with the primer pair glmS1-DN (5′-GTTCGTCGTCCACTGGGATCA) and Tn7L (PTn7L, 5′-ATTAGCTTACGACGCTACACCC) as previously described (51).

Statistical analysis.Student’s t test and one-way analysis of variance (ANOVA) (under R version 3.4.2) were applied to determine the statistical significance between target and control samples. The Wilcoxon rank test was used to confirm the normal distributions for multiple-sample comparisons with a family error of <0.05.

Data availability.All transcriptomics data have been submitted to GEO under accession number GSE127838 and NCBI tracking system no. 19795101.

ACKNOWLEDGMENTS

This work was funded by a grant to E.H.-G. from the Defense Threat Reduction Agency (DTRA) of the U.S. Department of Defense (DoD). The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

We thank Colin Manoil (University of Washington, Seattle) for his gift of the three BTH_I2366 transposon mutants: BT02489, BT06658, and BT08684. We also thank Herbert Schweizer (University of Florida) for the pUC18T-mini-Tn7T-Tp-gfpmut3 and pTNS2 plasmids and William W. Metcalf (University of Illinois, Urbana-Champaign) for E. coli strain BW29427.

FOOTNOTES

    • Received October 2, 2019.
    • Accepted January 3, 2020.
  • This article followed an open peer review process. The review history can be read here.

This is a work of the U.S. Government and is not subject to copyright protection in the United States. Foreign copyrights may apply.

REFERENCES

  1. 1.↵
    1. Levin BR,
    2. Rozen DE
    . 2006. Non-inherited antibiotic resistance. Nat Rev Microbiol 4:556–562. doi:10.1038/nrmicro1445.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Cohen NR,
    2. Lobritz MA,
    3. Collins JJ
    . 2013. Microbial Persistence and the road to drug resistance. Cell Host Microbe 13:632–642. doi:10.1016/j.chom.2013.05.009.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Mulcahy LR,
    2. Burns JL,
    3. Lory S,
    4. Lewis K
    . 2010. Emergence of Pseudomonas aeruginosa strains producing high levels of persister cells in patients with cystic fibrosis. J Bacteriol 192:6191–6199. doi:10.1128/JB.01651-09.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Lidstrom ME,
    2. Konopka MC
    . 2010. The role of physiological heterogeneity in microbial population behavior. Nat Chem Biol 6:705–712. doi:10.1038/nchembio.436.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Keren I,
    2. Shah D,
    3. Spoering A,
    4. Kaldalu N,
    5. Lewis K
    . 2004. Specialized persister cells and the mechanism of multidrug tolerance in Escherichia coli. J Bacteriol 186:8172–8180. doi:10.1128/JB.186.24.8172-8180.2004.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Shah D,
    2. Zhang Z,
    3. Khodursky A,
    4. Kaldalu N,
    5. Kurg K,
    6. Lewis K
    . 2006. Persisters: a distinct physiological state of E. coli. BMC Microbiol 6:53. doi:10.1186/1471-2180-6-53.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Orman MA,
    2. Brynildsen MP
    . 2013. Establishment of a method to rapidly assay bacterial persister metabolism. Antimicrob Agents Chemother 57:4398–4409. doi:10.1128/AAC.00372-13.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Ulrich RL,
    2. DeShazer D,
    3. Kenny TA,
    4. Ulrich MP,
    5. Moravusova A,
    6. Opperman T,
    7. Bavari S,
    8. Bowlin TL,
    9. Moir DT,
    10. Panchal RG
    . 2013. Characterization of the Burkholderia thailandensis SOS response by using whole-transcriptome shotgun sequencing. Appl Environ Microbiol 79:5830–5843. doi:10.1128/AEM.00538-13.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Leung V,
    2. Lévesque CM
    . 2012. A stress-inducible quorum-sensing peptide mediates the formation of persister cells with noninherited multidrug tolerance. J Bacteriol 194:2265–2274. doi:10.1128/JB.06707-11.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Murakami K,
    2. Ono T,
    3. Viducic D,
    4. Kayama S,
    5. Mori M,
    6. Hirota K,
    7. Nemoto K,
    8. Miyake Y
    . 2005. Role for rpoS gene of Pseudomonas aeruginosa in antibiotic tolerance. FEMS Microbiol Lett 242:161–167. doi:10.1016/j.femsle.2004.11.005.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Leszczynska D,
    2. Matuszewska E,
    3. Kuczynska-Wisnik D,
    4. Furmanek-Blaszk B,
    5. Laskowska E
    . 2013. The formation of persister cells in stationary-phase cultures of Escherichia coli is associated with the aggregation of endogenous proteins. PLoS One 8:e54737. doi:10.1371/journal.pone.0054737.
    OpenUrlCrossRefPubMed
  12. 12.↵
    1. Javid B,
    2. Sorrentino F,
    3. Toosky M,
    4. Zheng W,
    5. Pinkham JT,
    6. Jain N,
    7. Pan M,
    8. Deighan P,
    9. Rubin EJ
    . 2014. Mycobacterial mistranslation is necessary and sufficient for rifampicin phenotypic resistance. Proc Natl Acad Sci U S A 111:1132–1137. doi:10.1073/pnas.1317580111.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Mordukhova EA,
    2. Pan J-G
    . 2014. Stabilization of homoserine-O-succinyltransferase (MetA) decreases the frequency of persisters in Escherichia coli under stressful conditions. PLoS One 9:e110504. doi:10.1371/journal.pone.0110504.
    OpenUrlCrossRef
  14. 14.↵
    1. Van Melderen L,
    2. Saavedra De Bast M
    . 2009. Bacterial toxin-antitoxin systems: more than selfish entities? PLoS Genet 5:e1000437. doi:10.1371/journal.pgen.1000437.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Gerdes K,
    2. Maisonneuve E
    . 2012. Bacterial persistence and toxin-antitoxin loci. Annu Rev Microbiol 66:103–123. doi:10.1146/annurev-micro-092611-150159.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Wakamoto Y,
    2. Dhar N,
    3. Chait R,
    4. Schneider K,
    5. Signorino-Gelo F,
    6. Leibler S,
    7. McKinney JD
    . 2013. Dynamic persistence of antibiotic-stressed mycobacteria. Science 339:91–95. doi:10.1126/science.1229858.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Loris R,
    2. Garcia-Pino A
    . 2014. Disorder- and dynamics-based regulatory mechanisms in toxin-antitoxin modules. Chem Rev 114:6933–6947. doi:10.1021/cr400656f.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Conlon BP,
    2. Rowe SE,
    3. Gandt AB,
    4. Nuxoll AS,
    5. Donegan NP,
    6. Zalis EA,
    7. Clair G,
    8. Adkins JN,
    9. Cheung AL,
    10. Lewis K
    . 2016. Persister formation in Staphylococcus aureus is associated with ATP depletion. Nat Microbiol 1:16051. doi:10.1038/nmicrobiol.2016.51.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Keren I,
    2. Minami S,
    3. Rubin E,
    4. Lewis K
    . 2011. Characterization and transcriptome analysis of Mycobacterium tuberculosis persisters. mBio 2:e00100-11. doi:10.1128/mBio.00100-11.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Torrey HL,
    2. Keren I,
    3. Via LE,
    4. Lee JS,
    5. Lewis K
    . 2016. High persister mutants in Mycobacterium tuberculosis. PLoS One 11:e0155127. doi:10.1371/journal.pone.0155127.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Wood DN,
    2. Chaussee MA,
    3. Chaussee MS,
    4. Buttaro BA
    . 2005. Persistence of Streptococcus pyogenes in stationary-phase cultures. J Bacteriol 187:3319–3328. doi:10.1128/JB.187.10.3319-3328.2005.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Limmathurotsakul D,
    2. Golding N,
    3. Dance DA,
    4. Messina JP,
    5. Pigott DM,
    6. Moyes CL,
    7. Rolim DB,
    8. Bertherat E,
    9. Day NP,
    10. Peacock SJ,
    11. Hay SI
    . 2016. Predicted global distribution of Burkholderia pseudomallei and burden of melioidosis. Nat Microbiol 1:15008. doi:10.1038/nmicrobiol.2015.8.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Kim HS,
    2. Schell MA,
    3. Yu Y,
    4. Ulrich RL,
    5. Sarria SH,
    6. Nierman WC,
    7. DeShazer D
    . 2005. Bacterial genome adaptation to niches: divergence of the potential virulence genes in three Burkholderia species of different survival strategies. BMC Genomics 6:174. doi:10.1186/1471-2164-6-174.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Micheva-Viteva SN,
    2. Ross BN,
    3. Gao J,
    4. Adikari S,
    5. Zhang P,
    6. Mourant JR,
    7. Wu TH,
    8. Werner JH,
    9. Torres AG,
    10. Hong-Geller E
    . 2018. Increased mortality in mice following immunoprophylaxis therapy with high dosage of nicotinamide in Burkholderia persistent infections. Infect Immun 87:e00592-18. doi:10.1128/IAI.00592-18.
    OpenUrlCrossRef
  25. 25.↵
    1. Rhodes KA,
    2. Schweizer HP
    . 2016. Antibiotic resistance in Burkholderia species. Drug Resist Updat 28:82–90. doi:10.1016/j.drup.2016.07.003.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Fazli M,
    2. Almblad H,
    3. Rybtke ML,
    4. Givskov M,
    5. Eberl L,
    6. Tolker-Nielsen T
    . 2014. Regulation of biofilm formation in Pseudomonas and Burkholderia species. Environ Microbiol 16:1961–1981. doi:10.1111/1462-2920.12448.
    OpenUrlCrossRef
  27. 27.↵
    1. Eberl L
    . 2006. Quorum sensing in the genus Burkholderia. Int J Med Microbiol 296:103–110. doi:10.1016/j.ijmm.2006.01.035.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    1. Gierok P,
    2. Kohler C,
    3. Steinmetz I,
    4. Lalk M
    . 2016. Burkholderia pseudomallei colony morphotypes show a synchronized metabolic pattern after acute infection. PLoS Negl Trop Dis 10:e0004483. doi:10.1371/journal.pntd.0004483.
    OpenUrlCrossRef
  29. 29.↵
    1. Biggins JB,
    2. Gleber CD,
    3. Brady SF
    . 2011. Acyldepsipeptide HDAC inhibitor production induced in Burkholderia thailandensis. Org Lett 13:1536–1539. doi:10.1021/ol200225v.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Roger T,
    2. Lugrin J,
    3. Le Roy D,
    4. Goy G,
    5. Mombelli M,
    6. Koessler T,
    7. Ding XC,
    8. Chanson A-L,
    9. Reymond MK,
    10. Miconnet I,
    11. Schrenzel J,
    12. François P,
    13. Calandra T
    . 2011. Histone deacetylase inhibitors impair innate immune responses to Toll-like receptor agonists and to infection. Blood 117:1205–1217. doi:10.1182/blood-2010-05-284711.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Kuroda A,
    2. Murphy H,
    3. Cashel M,
    4. Kornberg A
    . 1997. Guanosine tetra- and pentaphosphate promote accumulation of inorganic polyphosphate in Escherichia coli. J Biol Chem 272:21240–21243. doi:10.1074/jbc.272.34.21240.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Meylan S,
    2. Porter CBM,
    3. Yang JH,
    4. Belenky P,
    5. Gutierrez A,
    6. Lobritz MA,
    7. Park J,
    8. Kim SH,
    9. Moskowitz SM,
    10. Collins JJ
    . 2017. Carbon sources tune antibiotic susceptibility in Pseudomonas aeruginosa via tricarboxylic acid cycle control. Cell Chem Biol 24:195–206. doi:10.1016/j.chembiol.2016.12.015.
    OpenUrlCrossRef
  33. 33.↵
    1. Harms A,
    2. Maisonneuve E,
    3. Gerdes K
    . 2016. Mechanisms of bacterial persistence during stress and antibiotic exposure. Science 354:aaf4268. doi:10.1126/science.aaf4268.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Liu Z,
    2. Hu K,
    3. Tang Y,
    4. Li H,
    5. Tang H,
    6. Hu X,
    7. Ma X,
    8. Liu Z
    . 2018. SmpB down-regulates proton-motive force for the persister tolerance to aminoglycosides in Aeromonas veronii. Biochem Biophys Res Commun 507:407–413. doi:10.1016/j.bbrc.2018.11.052.
    OpenUrlCrossRef
  35. 35.↵
    1. Allison KR,
    2. Brynildsen MP,
    3. Collins JJ
    . 2011. Metabolite-enabled eradication of bacterial persisters by aminoglycosides. Nature 473:216–220. doi:10.1038/nature10069.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Cheng AC,
    2. Fisher DA,
    3. Anstey NM,
    4. Stephens DP,
    5. Jacups SP,
    6. Currie BJ
    . 2004. Outcomes of patients with melioidosis treated with meropenem. Antimicrob Agents Chemother 48:1763–1765. doi:10.1128/aac.48.5.1763-1765.2004.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Springer MT,
    2. Singh VK,
    3. Cheung AL,
    4. Donegan NP,
    5. Chamberlain NR
    . 2016. Effect of clpP and clpC deletion on persister cell number in Staphylococcus aureus. J Med Microbiol 65:848–857. doi:10.1099/jmm.0.000304.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Thomsen LE,
    2. Olsen JE,
    3. Foster JW,
    4. Ingmer H
    . 2002. ClpP is involved in the stress response and degradation of misfolded proteins in Salmonella enterica serovar Typhimurium. Microbiology 148:2727–2733. doi:10.1099/00221287-148-9-2727.
    OpenUrlCrossRefPubMedWeb of Science
  39. 39.↵
    1. Krüger E,
    2. Witt E,
    3. Ohlmeier S,
    4. Hanschke R,
    5. Hecker M
    . 2000. The Clp proteases of Bacillus subtilis are directly involved in degradation of misfolded proteins. J Bacteriol 182:3259–3265. doi:10.1128/jb.182.11.3259-3265.2000.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Culp E,
    2. Wright GD
    . 2017. Bacterial proteases, untapped antimicrobial drug targets. J Antibiot (Tokyo) 70:366–377. doi:10.1038/ja.2016.138.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Petrova OE,
    2. Schurr JR,
    3. Schurr MJ,
    4. Sauer K
    . 2011. The novel Pseudomonas aeruginosa two-component regulator BfmR controls bacteriophage-mediated lysis and DNA release during biofilm development through PhdA. Mol Microbiol 81:767–783. doi:10.1111/j.1365-2958.2011.07733.x.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Secor PR,
    2. Sweere JM,
    3. Michaels LA,
    4. Malkovskiy AV,
    5. Lazzareschi D,
    6. Katznelson E,
    7. Rajadas J,
    8. Birnbaum ME,
    9. Arrigoni A,
    10. Braun KR,
    11. Evanko SP,
    12. Stevens DA,
    13. Kaminsky W,
    14. Singh PK,
    15. Parks WC,
    16. Bollyky PL
    . 2015. Filamentous bacteriophage promote biofilm assembly and function. Cell Host Microbe 18:549–559. doi:10.1016/j.chom.2015.10.013.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Fujita Y,
    2. Matsuoka H,
    3. Hirooka K
    . 2007. Regulation of fatty acid metabolism in bacteria. Mol Microbiol 66:829–839. doi:10.1111/j.1365-2958.2007.05947.x.
    OpenUrlCrossRefPubMedWeb of Science
  44. 44.↵
    1. DiRusso CC,
    2. Heimert TL,
    3. Metzger AK
    . 1992. Characterization of FadR, a global transcriptional regulator of fatty acid metabolism in Escherichia coli. Interaction with the fadB promoter is prevented by long chain fatty acyl coenzyme A. J Biol Chem 267:8685–8691.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Farewell A,
    2. Diez AA,
    3. DiRusso CC,
    4. Nyström T
    . 1996. Role of the Escherichia coli FadR regulator in stasis survival and growth phase-dependent expression of the uspA, fad, and fab genes. J Bacteriol 178:6443–6450. doi:10.1128/jb.178.22.6443-6450.1996.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Gallagher LA,
    2. Ramage E,
    3. Patrapuvich R,
    4. Weiss E,
    5. Brittnacher M,
    6. Manoil C
    . 2013. Sequence-defined transposon mutant library of Burkholderia thailandensis. mBio 4:e00604-13. doi:10.1128/mBio.00604-13.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Bar-Ziv R,
    2. Voichek Y,
    3. Barkai N
    . 2016. Chromatin dynamics during DNA replication. Genome Res 26:1245–1256. doi:10.1101/gr.201244.115.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Zhang Q,
    2. Zhou A,
    3. Li S,
    4. Ni J,
    5. Tao J,
    6. Lu J,
    7. Wan B,
    8. Li S,
    9. Zhang J,
    10. Zhao S,
    11. Zhao G-P,
    12. Shao F,
    13. Yao Y-F
    . 2016. Reversible lysine acetylation is involved in DNA replication initiation by regulating activities of initiator DnaA in Escherichia coli. Sci Rep 6:30837. doi:10.1038/srep30837.
    OpenUrlCrossRef
  49. 49.↵
    1. Cronan JE,
    2. Thomas J
    . 2009. Bacterial fatty acid synthesis and its relationships with polyketide synthetic pathways. Methods Enzymol 459:395–433. doi:10.1016/S0076-6879(09)04617-5.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Choi K-H,
    2. Schweizer HP
    . 2006. Mini-Tn7 insertion in bacteria with single attTn7 sites: example Pseudomonas aeruginosa. Nat Protoc 1:153–161. doi:10.1038/nprot.2006.24.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Choi K-H,
    2. Mima T,
    3. Casart Y,
    4. Rholl D,
    5. Kumar A,
    6. Beacham IR,
    7. Schweizer HP
    . 2008. Genetic tools for select-agent-compliant manipulation of Burkholderia pseudomallei. Appl Environ Microbiol 74:1064–1075. doi:10.1128/AEM.02430-07.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Jõers A,
    2. Kaldalu N,
    3. Tenson T
    . 2010. The frequency of persisters in Escherichia coli reflects the kinetics of awakening from dormancy. J Bacteriol 192:3379–3384. doi:10.1128/JB.00056-10.
    OpenUrlAbstract/FREE Full Text
  53. 53.↵
    1. Johnson PJT,
    2. Levin BR
    . 2013. Pharmacodynamics, population dynamics, and the evolution of persistence in Staphylococcus aureus. PLoS Genet 9:e1003123. doi:10.1371/journal.pgen.1003123.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Van den Bergh B,
    2. Fauvart M,
    3. Michiels J
    . 2017. Formation, physiology, ecology, evolution and clinical importance of bacterial persisters. FEMS Microbiol Rev 41:219–251. doi:10.1093/femsre/fux001.
    OpenUrlCrossRef
  55. 55.↵
    1. Smith JL,
    2. Sherman DH
    . 2008. An enzyme assembly line. Science 321:1304–1305. doi:10.1126/science.1163785.
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    1. Bernal V,
    2. Castaño-Cerezo S,
    3. Gallego-Jara J,
    4. Écija-Conesa A,
    5. de Diego T,
    6. Iborra JL,
    7. Cánovas M
    . 2014. Regulation of bacterial physiology by lysine acetylation of proteins. N Biotechnol 31:586–595. doi:10.1016/j.nbt.2014.03.002.
    OpenUrlCrossRef
  57. 57.↵
    1. Kim G-W,
    2. Yang X-J
    . 2011. Comprehensive lysine acetylomes emerging from bacteria to humans. Trends Biochem Sci 36:211–220. doi:10.1016/j.tibs.2010.10.001.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Okada BK,
    2. Wu Y,
    3. Mao D,
    4. Bushin LB,
    5. Seyedsayamdost MR
    . 2016. Mapping the trimethoprim-induced secondary metabolome of Burkholderia thailandensis. ACS Chem Biol 11:2124–2130. doi:10.1021/acschembio.6b00447.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. Liu X,
    2. Cheng Y-Q
    . 2014. Genome-guided discovery of diverse natural products from Burkholderia sp. J Ind Microbiol Biotechnol 41:275–284. doi:10.1007/s10295-013-1376-1.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Biggins JB,
    2. Liu X,
    3. Feng Z,
    4. Brady SF
    . 2011. Metabolites from the induced expression of cryptic single operons found in the genome of Burkholderia pseudomallei. J Am Chem Soc 133:1638–1641. doi:10.1021/ja1087369.
    OpenUrlCrossRefPubMedWeb of Science
  61. 61.↵
    1. Nougayrède J-P,
    2. Homburg S,
    3. Taieb F,
    4. Boury M,
    5. Brzuszkiewicz E,
    6. Gottschalk G,
    7. Buchrieser C,
    8. Hacker J,
    9. Dobrindt U,
    10. Oswald E
    . 2006. Escherichia coli induces DNA double-strand breaks in eukaryotic cells. Science 313:848–851. doi:10.1126/science.1127059.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Cuevas-Ramos G,
    2. Petit CR,
    3. Marcq I,
    4. Boury M,
    5. Oswald E,
    6. Nougayrède J-P
    . 2010. Escherichia coli induces DNA damage in vivo and triggers genomic instability in mammalian cells. Proc Natl Acad Sci U S A 107:11537–11542. doi:10.1073/pnas.1001261107.
    OpenUrlAbstract/FREE Full Text
  63. 63.↵
    1. Aoki M,
    2. Kobayashi M,
    3. Ishikawa J,
    4. Saita Y,
    5. Terai Y,
    6. Takayama K,
    7. Miyata K,
    8. Yamada T
    . 2000. A novel phosphodiesterase type 4 inhibitor, YM976 (4–(3-chlorophenyl)-1,7-diethylpyrido[2,3-d]pyrimidin-2(1H)-one), with little emetogenic activity. J Pharmacol Exp Ther 295:255–260.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    1. Banner KH,
    2. Press NJ
    . 2009. Dual PDE3/4 inhibitors as therapeutic agents for chronic obstructive pulmonary disease. Br J Pharmacol 157:892–906. doi:10.1111/j.1476-5381.2009.00170.x.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Maier T,
    2. Leibundgut M,
    3. Ban N
    . 2008. The crystal structure of a mammalian fatty acid synthase. Science 321:1315–1322. doi:10.1126/science.1161269.
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    1. Raetz CRH,
    2. Whitfield C
    . 2002. Lipopolysaccharide endotoxins. Annu Rev Biochem 71:635–700. doi:10.1146/annurev.biochem.71.110601.135414.
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.↵
    1. Whitehead NA,
    2. Barnard AML,
    3. Slater H,
    4. Simpson NJL,
    5. Salmond G
    . 2001. Quorum-sensing in Gram-negative bacteria. FEMS Microbiol Rev 25:365–404. doi:10.1111/j.1574-6976.2001.tb00583.x.
    OpenUrlCrossRefPubMedWeb of Science
  68. 68.↵
    1. Salaemae W,
    2. Booker GW,
    3. Polyak SW
    . 2016. The role of biotin in bacterial physiology and virulence: a novel antibiotic target for Mycobacterium tuberculosis. Microbiol Spectr 4(2):VMBF-0008-2015. doi:10.1128/microbiolspec.VMBF-0008-2015.
    OpenUrlCrossRef
  69. 69.↵
    1. Agnihotri G,
    2. Liu H-W
    . 2003. Enoyl-CoA hydratase: reaction, mechanism, and inhibition. Bioorg Med Chem 11:9–20. doi:10.1016/s0968-0896(02)00333-4.
    OpenUrlCrossRefPubMed
  70. 70.↵
    1. Young J,
    2. Holland IB
    . 1999. ABC transporters: bacterial exporters-revisited five years on. Biochim Biophys Acta Biomembr 1461:177–200. doi:10.1016/S0005-2736(99)00158-3.
    OpenUrlCrossRefPubMedWeb of Science
  71. 71.↵
    1. Li P-E,
    2. Lo C-C,
    3. Anderson J,
    4. Davenport KW,
    5. Bishop-Lilly KA,
    6. Xu Y,
    7. Ahmed S,
    8. Feng S,
    9. Mokashi VP,
    10. Chain PSG
    . 2017. Enabling the democratization of the genomics revolution with a fully integrated web-based bioinformatics platform. Nucleic Acids Res 45:67–80. doi:10.1093/nar/gkw1027.
    OpenUrlCrossRefPubMed
  72. 72.↵
    1. Lo C-C,
    2. Chain P
    . 2014. Rapid evaluation and quality control of next generation sequencing data with FaQCs. BMC Bioinformatics 15:366. doi:10.1186/s12859-014-0366-2.
    OpenUrlCrossRef
  73. 73.↵
    1. Kim D,
    2. Langmead B,
    3. Salzberg SL
    . 2015. HISAT: a fast spliced aligner with low memory requirements. Nat Methods 12:357–360. doi:10.1038/nmeth.3317.
    OpenUrlCrossRefPubMed
  74. 74.↵
    1. Liao Y,
    2. Smyth GK,
    3. Shi W
    . 2014. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30:923–930. doi:10.1093/bioinformatics/btt656.
    OpenUrlCrossRefPubMedWeb of Science
  75. 75.↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    . 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15:550. doi:10.1186/s13059-014-0550-8.
    OpenUrlCrossRefPubMed
  76. 76.↵
    1. Nishino K,
    2. Yamada J,
    3. Hirakawa H,
    4. Hirata T,
    5. Yamaguchi A
    . 2003. Roles of TolC-dependent multidrug transporters of Escherichia coli in resistance to beta-lactams. Antimicrob Agents Chemother 47:3030–3033. doi:10.1128/aac.47.9.3030-3033.2003.
    OpenUrlAbstract/FREE Full Text
  77. 77.↵
    1. Wu S,
    2. Yu P-L,
    3. Wheeler D,
    4. Flint S
    . 2018. Transcriptomic study on persistence and survival of Listeria monocytogenes following lethal treatment with nisin. J Glob Antimicrob Resist 15:25–31. doi:10.1016/j.jgar.2018.06.003.
    OpenUrlCrossRef
  78. 78.↵
    1. Kamal F,
    2. Dennis JJ
    . 2015. Burkholderia cepacia complex phage-antibiotic synergy (PAS): antibiotics stimulate lytic phage activity. Appl Environ Microbiol 81:1132–1138. doi:10.1128/AEM.02850-14.
    OpenUrlAbstract/FREE Full Text
  79. 79.↵
    1. Olszak T,
    2. Latka A,
    3. Rozniowski B,
    4. Valvano MA,
    5. Drulis-Kawa Z
    . 2017. Phage life cycles behind bacterial biodiversity. Curr Med Chem 24:3987–4001. doi:10.2174/0929867324666170413100136.
    OpenUrlCrossRef
View Abstract
PreviousNext
Back to top
Download PDF
Citation Tools
A Gene Cluster That Encodes Histone Deacetylase Inhibitors Contributes to Bacterial Persistence and Antibiotic Tolerance in Burkholderia thailandensis
Sofiya N. Micheva-Viteva, Migun Shakya, Samantha H. Adikari, Cheryl D. Gleasner, Nileena Velappan, Judith R. Mourant, Patrick S. G. Chain, Elizabeth Hong-Geller
mSystems Feb 2020, 5 (1) e00609-19; DOI: 10.1128/mSystems.00609-19

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
A Gene Cluster That Encodes Histone Deacetylase Inhibitors Contributes to Bacterial Persistence and Antibiotic Tolerance in Burkholderia thailandensis
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
A Gene Cluster That Encodes Histone Deacetylase Inhibitors Contributes to Bacterial Persistence and Antibiotic Tolerance in Burkholderia thailandensis
Sofiya N. Micheva-Viteva, Migun Shakya, Samantha H. Adikari, Cheryl D. Gleasner, Nileena Velappan, Judith R. Mourant, Patrick S. G. Chain, Elizabeth Hong-Geller
mSystems Feb 2020, 5 (1) e00609-19; DOI: 10.1128/mSystems.00609-19
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • Review History
  • PDF

KEYWORDS

Burkholderia
antibiotic resistance
bacterial persistence
stress
transcriptomics

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077