Deep transcriptome profiling of clinicalKlebsiella pneumoniaeisolates reveals strain and sequence type-specific adaptation

Transcriptomic Insights into Metabolism-Dependent Biosynthesis of Bacterial Nanocellulose

Bacterial nanocellulose (BNC) is an attractive green-synthesized biomaterial for biomedical applications and various other applications. However, effective engineering of BNC production has been limited by our poor knowledge of the related metabolic processes. In contrast to the traditional perception that genome critically determines biosynthesis behaviors, here we discover that the glucose metabolism could also drastically affect the BNC synthesis in Gluconacetobacter hansenii. The transcriptomic profiles of two model BNC-producing strains, G. hansenii ATCC 53582 and ATCC 23769, which have highly similar genomes but drastically different BNC yields, were compared. The results show that their BNC synthesis capacities were highly related to metabolic activities such as ATP synthesis, ion transport protein assembly, and carbohydrate metabolic processes, confirming an important role of metabolism-related transcriptomes in governing the BNC yield. Our findings provide insights into the microbial biosynthesis behaviors from a transcriptome perspective, potentially guiding cellular engineering for biomaterial synthesis.

The global RNA-RNA interactome of Klebsiella pneumoniae unveils a small RNA regulator of cell division

The ubiquitous RNA chaperone Hfq is involved in the regulation of key biological processes in many species across the bacterial kingdom. In the opportunistic human pathogen Klebsiella pneumoniae, deletion of the hfq gene affects the global transcriptome, virulence, and stress resistance; however, the ligands of the major RNA-binding protein in this species have remained elusive. In this study, we have combined transcriptomic, co-immunoprecipitation, and global RNA interactome analyses to compile an inventory of conserved and species-specific RNAs bound by Hfq and to monitor Hfq-mediated RNA-RNA interactions. In addition to dozens of RNA-RNA pairs, our study revealed an Hfq-dependent small regulatory RNA (sRNA), DinR, which is processed from the 3' terminal portion of dinI mRNA. Transcription of dinI is controlled by the master regulator of the SOS response, LexA. As DinR accumulates in K. pneumoniae in response to DNA damage, the sRNA represses translation of the ftsZ transcript by occupation of the ribosome binding site. Ectopic overexpression of DinR causes depletion of ftsZ mRNA and inhibition of cell division, while deletion of dinR antagonizes cell elongation in the presence of DNA damage. Collectively, our work highlights the important role of RNA-based gene regulation in K. pneumoniae and uncovers the central role of DinR in LexA-controlled division inhibition during the SOS response.

Genome-scale metabolic modeling reveals increased reliance on valine catabolism in clinical isolates of Klebsiella pneumoniae

Infections due to carbapenem-resistant Enterobacteriaceae have recently emerged as one of the most urgent threats to hospitalized patients within the United States and Europe. By far the most common etiological agent of these infections is Klebsiella pneumoniae, frequently manifesting in hospital-acquired pneumonia with a mortality rate of ~50% even with antimicrobial intervention. We performed transcriptomic analysis of data collected previously from in vitro characterization of both laboratory and clinical isolates which revealed shifts in expression of multiple master metabolic regulators across isolate types. Metabolism has been previously shown to be an effective target for antibacterial therapy, and genome-scale metabolic network reconstructions (GENREs) have provided a powerful means to accelerate identification of potential targets in silico. Combining these techniques with the transcriptome meta-analysis, we generated context-specific models of metabolism utilizing a well-curated GENRE of K. pneumoniae (iYL1228) to identify novel therapeutic targets. Functional metabolic analyses revealed that both composition and metabolic activity of clinical isolate-associated context-specific models significantly differs from laboratory isolate-associated models of the bacterium. Additionally, we identified increased catabolism of L-valine in clinical isolate-specific growth simulations. These findings warrant future studies for potential efficacy of valine transaminase inhibition as a target against K. pneumoniae infection.

A Novel Pathway of Chlorimuron-Ethyl Biodegradation by Chenggangzhangella methanolivorans Strain CHL1 and Its Molecular Mechanisms

Chlorimuron-ethyl is a widely used herbicide in agriculture. However, uncontrolled chlorimuron-ethyl application causes serious environmental problems. Chlorimuron-ethyl can be effectively degraded by microbes, but the underlying molecular mechanisms are not fully understood. In this study, we identified the possible pathways and key genes involved in chlorimuron-ethyl degradation by the Chenggangzhangella methanolivorans strain CHL1, a Methylocystaceae strain with the ability to degrade sulfonylurea herbicides. Using a metabolomics method, eight intermediate degradation products were identified, and three pathways, including a novel pyrimidine-ring-opening pathway, were found to be involved in chlorimuron-ethyl degradation by strain CHL1. Transcriptome sequencing indicated that three genes (atzF, atzD, and cysJ) are involved in chlorimuron-ethyl degradation by strain CHL1. The gene knock-out and complementation techniques allowed for the functions of the three genes to be identified, and the enzymes involved in the different steps of chlorimuron-ethyl degradation pathways were preliminary predicted. The results reveal a previously unreported pathway and the key genes of chlorimuron-ethyl degradation by strain CHL1, which have implications for attempts to enrich the biodegradation mechanism of sulfonylurea herbicides and to construct engineered bacteria in order to remove sulfonylurea herbicide residues from environmental media.

Prevalence of the fimbrial operon mrkABCD, mrkA expression, biofilm formation and effect of biocides on biofilm formation in carbapenemase-producing Klebsiella pneumoniae isolates belonging or not to high-risk clones

BACKGROUND

The role of mrkA adhesin expression, biofilm production, biofilm viability and biocides in the biofilms of carbapenemase-producing K. pneumoniae isolates was investigated.

METHODS

Seventeen isolates representing different sequence types and carbapenemases were investigated. mrkA expression was determined by real-time RT-PCR. Biofilm production (25°C and 37°C, with and without humidity) was determined by the crystal violet assay. The effect of isopropanol, povidone-iodine, sodium hypochlorite, chlorhexidine digluconate, benzalkonium chloride, ethanol and triclosan on biofilms was determined. The effect of povidone-iodine on biofilm biomass and thickness was also determined by Confocal Laser Scanning Microscopy (CLSM).

RESULTS

mrkA expression ranged 28.2-1.3 (high or intermediate-level; 64% of HR clones) and 21.5-1.3 (50% of non-HR clones). At 25°C biofilm formation was observed in 41% of isolates (absence of humidity) and 35% of isolates (presence of humidity), whereas at 37°C it was observed in 76% of isolates, with and without humidity. At 25°C biofilm producers were more frequently observed in HR clones (45% with humidity and 55% without humidity) than non-HR clones (17% with and without humidity). Biofilm viability from day 21 was higher at 25°C than 37°C. The greatest decrease in biofilm formation was observed with povidone iodine (29% decrease), which also decrease biofilm thickness.

CONCLUSIONS

Biofilm formation in carbapenemase-producing K. pneumoniae is related to mrkA expression. Biofilm formation is affected by temperature (37°C>25°C) whereas humidity has little effect. Biofilm viability is affected by temperature (25°C>37°C). At 25°C, HR clones are more frequently biofilm producers than non-HR clones. Povidone-iodine can decrease biofilm production and biofilm thickness.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

Antibiotic resistance in Pseudomonas aeruginosa - Mechanisms, epidemiology and evolution

Antibiotics are powerful drugs used in the treatment of bacterial infections. The inappropriate use of these medicines has driven the dissemination of antibiotic resistance (AR) in most bacteria. Pseudomonas aeruginosa is an opportunistic pathogen commonly involved in environmental- and difficult-to-treat hospital-acquired infections. This species is frequently resistant to several antibiotics, being in the "critical" category of the WHO's priority pathogens list for research and development of new antibiotics. In addition to a remarkable intrinsic resistance to several antibiotics, P. aeruginosa can acquire resistance through chromosomal mutations and acquisition of AR genes. P. aeruginosa has one of the largest bacterial genomes and possesses a significant assortment of genes acquired by horizontal gene transfer (HGT), which are frequently localized within integrons and mobile genetic elements (MGEs), such as transposons, insertion sequences, genomic islands, phages, plasmids and integrative and conjugative elements (ICEs). This genomic diversity results in a non-clonal population structure, punctuated by specific clones that are associated with significant morbidity and mortality worldwide, the so-called high-risk clones. Acquisition of MGEs produces a fitness cost in the host, that can be eased over time by compensatory mutations during MGE-host coevolution. Even though plasmids and ICEs are important drivers of AR, the underlying evolutionary traits that promote this dissemination are poorly understood. In this review, we provide a comprehensive description of the main strategies involved in AR in P. aeruginosa and the leading drivers of HGT in this species. The most recently developed genomic tools that allowed a better understanding of the features contributing for the success of P. aeruginosa are discussed.

Fatty Acid Methyl Ester (FAME) Profiling Identifies Carbapenemase-Producing Klebsiella pneumoniae Belonging to Clonal Complex 258

Carbapenem-resistant Klebsiella pneumoniae (CRKP) is one of the most extensively antibiotic-resistant pathogens encountered in the clinical setting today. A few studies to-date suggest that CRKP and carbapenem-susceptible K. pneumoniae (CSKP) differ from one another not only with respect to their underlying genetics, but also their transcriptomic and metabolomic fingerprints. Within this context, we characterize the fatty acid methyl ester (FAME) profiles of these pathogens in vitro. Specifically, we evaluated the FAME profiles of six Klebsiella pneumoniae carbapenemase (KPC)-producing isolates belonging to the CC258 lineage (KPC+/258+), six KPC-producing isolates belonging to non-CC258 lineages (KPC+/258−), and six non-KPC-producing isolates belonging to non-CC258 lineages (KPC−/258−). We utilized a single-step sample preparation method to simultaneously lyse bacterial cells and transesterify the lipid fraction, and identified 14 unique FAMEs using gas chromatography-mass spectrometry. The machine learning algorithm Random Forest identified four FAMEs that were highly discriminatory between CC258 and non-CC258 isolates (9(Z)-octadecenoate, 2-phenylacetate, pentadecanoate, and hexadecanoate), of which three were also significantly different in relative abundance between these two groups. These findings suggest that distinct differences exist between CC258 and non-CC258 K. pneumoniae isolates with respect to the metabolism of both fatty acids and amino acids, a hypothesis that is supported by previously-acquired transcriptomic data.

Phenotypical profile and global transcriptomic profile of Hypervirulent Klebsiella pneumoniae due to carbapenemase-encoding plasmid acquisition

Background

Plasmids play an vital role in driving the rapid global spread of antimicrobial resistance and adaptation to changing ambient conditions. It has been suggested that the presence of plasmids can pose tremendous impacts on the host physiology. However, little is known regarding the contributions of carbapenemase-encoding plasmid carriage on the physiology and pathogenicity of hypervirulent K. pneumoniae (hvKP).

Results

Here we performed a transcriptomic analysis of hvKP with or without carbapenemase-encoding plasmid p24835-NDM5. The results had shown 683 genes with differential expression (false discovery rate, ≤0.001; > 2-fold change), of which 107 were up-regulated and 576 were down-regulated. Gene groups with functions relating to carbohydrate metabolism and multidrug efflux system were increased in genes with increased expression, and those relating to capsule biosynthesis and virulence factors were increased in the genes with decreased expression. In agreement with these changes, survival rate of TfpNDM-hvKP in the presence of normal human serum decreased, and competitive index (CI values) indicated significant fitness defects in the plasmid-carrying hvKP strain when co-cultured with its plasmid-free isogenic ancestor and the ATCC control. Moreover, the p24835-NDM5-containing hvKP strain retained its high neutrophil-mediated phagocytosis and murine lethality.

Conclusion

These data indicate that hvKP responds to carbapenemase-encoding plasmid by altering the expression of genes involved in carbohydrate metabolism, antibiotic resistance, capsule biosynthesis and virulence expression. Apart from antibiotic resistance selective advantages, carbapenemase-encoding plasmid carriage may also lead to virulence change or adaption to specific habitats in hvKP strain.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5705-2) contains supplementary material, which is available to authorized users.

Integrated analysis of population genomics, transcriptomics and virulence provides novel insights into Streptococcus pyogenes pathogenesis

Streptococcus pyogenes causes 700 million human infections annually worldwide, yet, despite a century of intensive effort, there is no licensed vaccine against this bacterium. Although a number of large-scale genomic studies of bacterial pathogens have been published, the relationships among the genome, transcriptome, and virulence in large bacterial populations remain poorly understood. We sequenced the genomes of 2,101 emm28 S. pyogenes invasive strains, from which we selected 492 phylogenetically diverse strains for transcriptome analysis and 50 strains for virulence assessment. Data integration provided a novel understanding of the virulence mechanisms of this model organism. Genome-wide association study, expression quantitative trait loci analysis, machine learning, and isogenic mutant strains identified and confirmed a one-nucleotide indel in an intergenic region that significantly alters global transcript profiles and ultimately virulence. The integrative strategy that we used is generally applicable to any microbe and may lead to new therapeutics for many human pathogens.

Detection of high-risk carbapenem-resistant Klebsiella pneumoniae and Enterobacter cloacae isolates using volatile molecular profiles

Infections caused by carbapenem-resistant Enterobacteriaceae (CRE) are alarming in the clinical setting, as CRE isolates often exhibit resistance to most clinically-available antibiotics. Klebsiella pneumoniae carbapenemase (KPC) is the most common carbapenemase carried by CRE in North America and Europe, frequently detected in isolates of K. pneumoniae, Escherichia coli, and Enterobacter cloacae. Notably, KPC-expressing strains often arise from clonal lineages, with sequence type 258 (ST258) representing the dominant lineage in K. pneumoniae, ST131 in E. coli, and ST78 and ST171 in E. cloacae. Prior studies have demonstrated that carbapenem-resistant K. pneumoniae differs from carbapenem-susceptible K. pneumoniae at both the transcriptomic and soluble metabolomic levels. In the present study, we sought to determine whether carbapenem-resistant and carbapenem-susceptible isolates of K. pneumoniae, E. coli, and E. cloacae produce distinct volatile metabolic profiles. We were able to identify a volatile metabolic fingerprint that could discriminate between CRE and non-CRE with an area under the receiver operating characteristic curve (AUROC) as high as 0.912. Species-specific AUROCs were as high as 0.988 for K. pneumoniae and 1.000 for E. cloacae. Paradoxically, curing of KPC-expressing plasmids from a subset of K. pneumoniae isolates further accentuated the metabolic differences observed between ST258 and non-ST258.

Whole Genome Sequence Analysis of CTX-M-15 Producing Klebsiella Isolates Allowed Dissecting a Polyclonal Outbreak Scenario

Extended-spectrum β-lactamase (ESBL) producing Klebsiella pneumoniae pose an important threat of infection with increased morbidity and mortality, especially for immunocompromised patients. Here, we use the rise of multidrug-resistant K. pneumoniae in a German neurorehabilitation center from April 2015 to April 2016 to dissect the benefit of whole genome sequencing (WGS) for outbreak analyses. In total, 53 isolates were obtained from 52 patients and examined using WGS. Two independent analysis strategies (reference-based and -free) revealed the same distinct clusters of two CTX-M-15 producing K. pneumoniae clones (ST15, n = 31; ST405, n = 7) and one CTX-M-15 producing Klebsiella quasipneumoniae strain (ST414, n = 8). Additionally, we determined sequence variations associated with antimicrobial resistance phenotypes in single isolates expressing carbapenem and colistin resistance, respectively. For rapid detection of the major K. pneumoniae outbreak clone (ST15), a selective triplex PCR was deduced from WGS data of the major outbreak strain and K. pneumoniae genome data deposited in central databases. Moreover, we introduce two novel open-source applications supporting reference genome selection (refRank; https://gitlab.com/s.fuchs/refRank) and alignment-based SNP-filtering (SNPfilter; https://gitlab.com/s.fuchs/snpfilter) in NGS analyses.

Evaluating Different Virulence Traits of Klebsiella pneumoniae Using Dictyostelium discoideum and Zebrafish Larvae as Host Models

Multiresistant and invasive hypervirulent Klebsiella pneumoniae strains have become one of the most urgent bacterial pathogen threats. Recent analyses revealed a high genomic plasticity of this species, harboring a variety of mobile genetic elements associated with virulent strains, encoding proteins of unknown function whose possible role in pathogenesis have not been addressed. K. pneumoniae virulence has been studied mainly in animal models such as mice and pigs, however, practical, financial, ethical and methodological issues limit the use of mammal hosts. Consequently, the development of simple and cost-effective experimental approaches with alternative host models is needed. In this work we described the use of both, the social amoeba and professional phagocyte Dictyostelium discoideum and the fish Danio rerio (zebrafish) as surrogate host models to study K. pneumoniae virulence. We compared three K. pneumoniae clinical isolates evaluating their resistance to phagocytosis, intracellular survival, lethality, intestinal colonization, and innate immune cells recruitment. Optical transparency of both host models permitted studying the infective process in vivo, following the Klebsiella-host interactions through live-cell imaging. We demonstrated that K. pneumoniae RYC492, but not the multiresistant strains 700603 and BAA-1705, is virulent to both host models and elicits a strong immune response. Moreover, this strain showed a high resistance to phagocytosis by D. discoideum, an increased ability to form biofilms and a more prominent and irregular capsule. Besides, the strain 700603 showed the unique ability to replicate inside amoeba cells. Genomic comparison of the K. pneumoniae strains showed that the RYC492 strain has a higher overall content of virulence factors although no specific genes could be linked to its phagocytosis resistance, nor to the intracellular survival observed for the 700603 strain. Our results indicate that both zebrafish and D. discoideum are advantageous host models to study different traits of K. pneumoniae that are associated with virulence.

Klebsiella pneumoniae: a major worldwide source and shuttle for antibiotic resistance

Klebsiella pneumoniae is an important multidrug-resistant (MDR) pathogen affecting humans and a major source for hospital infections associated with high morbidity and mortality due to limited treatment options. We summarize the wide resistome of this pathogen, which encompasses plentiful chromosomal and plasmid-encoded antibiotic resistance genes (ARGs). Under antibiotic selective pressure, K. pneumoniae continuously accumulates ARGs, by de novo mutations, and via acquisition of plasmids and transferable genetic elements, leading to extremely drug resistant (XDR) strains harboring a 'super resistome'. In the last two decades, numerous high-risk (HiR) MDR and XDR K. pneumoniae sequence types have emerged showing superior ability to cause multicontinent outbreaks, and continuous global dissemination. The data highlight the complex evolution of MDR and XDR K. pneumoniae, involving transfer and spread of ARGs, and epidemic plasmids in highly disseminating successful clones. With the worldwide catastrophe of antibiotic resistance and the urgent need to identify the main pathogens that pose a threat on the future of infectious diseases, further studies are warranted to determine the epidemic traits and plasmid acquisition in K. pneumoniae. There is a need for future genomic and translational studies to decipher specific targets in HiR clones to design targeted prevention and treatment.

Exploring the Genome and Phenotype of Multi-Drug Resistant Klebsiella pneumoniae of Clinical Origin

Klebsiella pneumoniae is an important nosocomial pathogen with an extraordinary resistant phenotype due to a combination of acquired resistant-elements and efflux mechanisms. In this study a detailed molecular characterization of 11 K. pneumoniae isolates of clinical origin was carried out. Eleven clinical isolates were tested for their susceptibilities, by disk diffusion and broth microdilution and interpreted according to CLSI guidelines. Efflux activity was determined by measuring the extrusion of ethidium bromide and biofilm formation was assessed following static growth in Müeller-Hinton and minimal media M9 broths at two temperatures and time points. Template DNA from all 11 isolates was extracted and sequenced. The study collection was found to be resistant to several (extended-spectrum beta-lactam) ESBL-type compounds along with several (fluoro)quinolones (FQ). Resistance to tetracycline accounted for 55% of the study collection (n = 6) and three of the 11 isolates were resistance to carbapenems. Genotyping identified bla_CTX-M-15 (82%), bla_SHV-12 (55%), and bla_TEM-1B (45%) ESBL encoding genes and FQ resistance was associated the presence of the oqxAB operon, identified in 10 of the 11 isolates and qnrB gene in one isolate. The polymorphisms detected in the quinolone resistance-determining regions (QRDRs) were associated with isolates of the clonal group CG15. Sequence types (ST) identified were representative of previously described clonal groups including CG258 (n = 7), CG15 (n = 3), and CG147 (n = 1). Plasmid replicon type databases were queried indicating the presence of IncFII and IncFIB replicon types in the majority of the isolates (91%), followed by IncFIA (45%), and IncR (45%). Two of the 11 isolates were found positive for yersiniabactin siderophore-encoding genes. No differences in the ability to efflux ethidium bromide were identified. Biofilm formation was stronger when the isolates were grown under stressed conditions at 37°C for a period up to 96 h. These data confirm the fact that well-recognized clonal groups of K. pneumoniae of importance to human health carries a diverse repertoire of antimicrobial resistance determinants, particularly related to critically important drugs in the ESBL and FQ classes. The capacity of most isolates to form strong biofilms, when stressed under laboratory-simulated conditions, supports the risk to human health associated with nosocomial infections deriving from indwelling medical devices.

The polymyxin B-induced transcriptomic response of a clinical, multidrug-resistant Klebsiella pneumoniae involves multiple regulatory elements and intracellular targets

Background

The emergence of multidrug-resistant Klebsiella pneumoniae is a major public health concern. Many K. pneumoniae infections can only be treated when resorting to last-line drugs such as polymyxin B (PB). However, resistance to this antibiotic is also observed, although insufficient information is described on its mode of action as well as the mechanisms used by resistant bacteria to evade its effects. We aimed to study PB resistance and the influence of abiotic stresses in a clinical K. pneumoniae strain using whole transcriptome profiling.

Results

We sequenced 12 cDNA libraries of K. pneumoniae Kp13 bacteria, from two biological replicates of the original strain Kp13 (Kp13) and five derivative strains: induced high-level PB resistance in acidic pH (Kp13_pH), magnesium deprivation (Kp13_Mg), high concentrations of calcium (Kp13_Ca) and iron (Kp13_Fe), and a control condition with PB (Kp13_PolB). Our results show the involvement of multiple regulatory loci that differentially respond to each condition as well as a shared gene expression response elicited by PB treatment, and indicate the participation of two-regulatory components such as ArcA-ArcB, which could be involved in re-routing the K. pneumoniae metabolism following PB treatment. Modules of co-expressed genes could be determined, which correlated to growth in acid stress and PB exposure. We hypothesize that polymyxin B induces metabolic shifts in K. pneumoniae that could relate to surviving against the action of this antibiotic.

Conclusions

We obtained whole transcriptome data for K. pneumoniae under different environmental conditions and PB treatment. Our results supports the notion that the K. pneumoniae response to PB exposure goes beyond damaged membrane reconstruction and involves recruitment of multiple gene modules and intracellular targets.

Electronic supplementary material

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

Klebsiella pneumoniae Population Genomics and Antimicrobial-Resistant Clones

Antimicrobial-resistant Klebsiella pneumoniae (Kp) has emerged as a major global public health problem. While resistance can occur across a broad range of Kp clones, a small number have become globally distributed and commonly cause outbreaks in hospital settings. Here we describe recent comparative genomics investigations that have shed light on Kp population structure and the evolution of antimicrobial-resistant clones. These studies provide the basic framework within which genomic epidemiology and evolution can be understood, but have merely scratched the surface of what can and should be explored. We assert that further large-scale comparative and functional genomics studies are urgently needed to better understand the biology of this clinically important bacterium.