Update and elucidation of Plasmodium kinomes: Prioritization of kinases as potential drug targets for malaria.

Malaria is a tropical disease caused by Plasmodium spp. and transmitted by the bite of infected Anopheles mosquitoes. Protein kinases (PKs) play key roles in the life cycle of the etiological agent of malaria, turning these proteins attractive targets for antimalarial drug discovery campaigns. As part of an effort to understand parasite signaling functions, we report the results of a bioinformatics pipeline analysis of PKs of eight Plasmodium species. To date, no P. malariae and P. ovale kinome assemble has been conducted. We classified, curated and annotated predicted kinases to update P. falciparum, P. vivax, P. yoelii, P. berghei, P. chabaudi, and P. knowlesi kinomes published to date, as well as report for the first time the kinomes of P. malariae and P. ovale. Overall, from 76 to 97 PKs were identified among all Plasmodium spp. kinomes. Most of the kinases were assigned to seven of nine major kinase groups: AGC, CAMK, CMGC, CK1, STE, TKL, OTHER; and the Plasmodium-specific group FIKK. About 30% of kinases have been deeply classified into group, family and subfamily levels and only about 10% remained unclassified. Furthermore, updating and comparing the kinomes of P. vivax and P. falciparum allowed for the prioritization and selection of kinases as potential drug targets that could be explored for discovering new drugs against malaria. This integrated approach resulted in the selection of 37 protein kinases as potential targets and the identification of investigational compounds with moderate in vitro activity against asexual P. falciparum (3D7 and Dd2 strains) stages that could serve as starting points for the search of potent antimalarial leads in the future.

Evaluation of the Inflammatory Cytokines and IL-10 Network in Individuals Co-infected With Human T-Cell Lymphotropic Virus and Hepatitis C Virus (HTLV/HCV).

BACKGROUND: Co-infection between the human T-cell lymphotropic virus (HTLV) and the hepatitis C virus (HCV) can modify the natural history of HCV infection. The aim of this study was to describe the inflammatory cytokines and IL-10 network in patients co-infected with HTLV and HCV viruses in Bahia, Brazil. METHODS: Samples from 31 HTLV/HCV co-infected individuals and 27 HCV monoinfected individuals were evaluated. IFN-γ, TNF-α, IL-10, IL-8, and IL-1 cytokines were quantified by ELISA. Clinical, laboratory data were obtained from patient records. Serum levels of the cytokines were log10-transformed and data mining was performed using Z-score statistics and correlation analysis. RESULTS: Co-infected individuals presented a tendency toward higher production of INF-γ compared to the HCV monoinfected group. Regarding cytokine pairs, there was a positive correlation (P-value < 0.05) between IL-1 and IL-8 in the HTLV/HCV co-infected group and uninfected controls, and two correlations in the HCV mono-infected group IL-8 - IL10 and IL- INF-γ - IL-10 pairs. There was no significant difference between the groups for the other parameters analyzed. CONCLUSION: The results presented herein indicated that HTLV/HCV co-infection was associated with a trend in IFN-γ production while HCV-infected individuals presented a positive correlation with both inflammatory cytokines (IL-8 and IFN-γ) and the regulatory cytokine IL-10.

LPG2 Gene Duplication in Leishmania infantum: A Case for CRISPR-Cas9 Gene Editing.

On the surface of the Leishmania promastigote, phosphoglycans (PG) such as lipophosphoglycan (LPG), proteophosphoglycan (PPG), free phosphoglycan polymers (PGs), and acid phosphatases (sAP), are dominant and contribute to the invasion and survival of Leishmania within the host cell by modulating macrophage signaling and intracellular trafficking. Phosphoglycan synthesis depends on the Golgi GDP-mannose transporter encoded by the LPG2 gene. Aiming to investigate the role of PG-containing molecules in Leishmania infantum infection process, herein we describe the generation and characterization of L. infantum LPG2-deficient parasites. This gene was unexpectedly identified as duplicated in the L. infantum genome, which impaired gene targeting using the conventional homologous recombination approach. This limitation was circumvented by the use of CRISPR/Cas9 technology. Knockout parasites were selected by agglutination assays using CA7AE antibodies followed by a lectin (RCA 120). Five clones were isolated and molecularly characterized, all revealing the expected edited genome, as well as the complete absence of LPG and PG-containing molecule expression. Finally, the deletion of LPG2 was found to impair the outcome of infection in human neutrophils, as demonstrated by a pronounced reduction (~83%) in intracellular load compared to wild-type parasite infection. The results obtained herein reinforce the importance of LPG and other PGs as virulence factors in host-parasite interactions.

Proteins involved in the biosynthesis of lipophosphoglycan in Leishmania: a comparative genomic and evolutionary analysis.

BACKGROUND: Leishmania spp. are digenetic parasites capable of infecting humans and causing a range of diseases collectively known as leishmaniasis. The main mechanisms involved in the development and permanence of this pathology are linked to evasion of the immune response. Crosstalk between the immune system and particularities of each pathogenic species is associated with diverse disease manifestations. Lipophosphoglycan (LPG), one of the most important molecules present on the surface of Leishmania parasites, is divided into four regions with high molecular variability. Although LPG plays an important role in host-pathogen and vector-parasite interactions, the distribution and phylogenetic relatedness of the genes responsible for its synthesis remain poorly explored. The recent availability of full genomes and transcriptomes of Leishmania parasites offers an opportunity to leverage insight on how LPG-related genes are distributed and expressed by these pathogens. RESULTS: Using a phylogenomics-based framework, we identified a catalog of genes involved in LPG biosynthesis across 22 species of Leishmania from the subgenera Viannia and Leishmania, as well as 5 non-Leishmania trypanosomatids. The evolutionary relationships of these genes across species were also evaluated. Nine genes related to the production of the glycosylphosphatidylinositol (GPI)-anchor were highly conserved among compared species, whereas 22 genes related to the synthesis of the repeat unit presented variable conservation. Extensive gain/loss events were verified, particularly in genes SCG1-4 and SCA1-2. These genes act, respectively, on the synthesis of the side chain attached to phosphoglycans and in the transfer of arabinose residues. Phylogenetic analyses disclosed evolutionary patterns reflective of differences in host specialization, geographic origin and disease manifestation. CONCLUSIONS: The multiple gene gain/loss events identified by genomic data mining help to explain some of the observed intra- and interspecies variation in LPG structure. Collectively, our results provide a comprehensive catalog that details how LPG-related genes evolved in the Leishmania parasite specialization process.

Leveraging User-Friendly Network Approaches to Extract Knowledge From High-Throughput Omics Datasets.

Recent technological advances for the acquisition of multi-omics data have allowed an unprecedented understanding of the complex intricacies of biological systems. In parallel, a myriad of computational analysis techniques and bioinformatics tools have been developed, with many efforts directed towards the creation and interpretation of networks from this data. In this review, we begin by examining key network concepts and terminology. Then, computational tools that allow for their construction and analysis from high-throughput omics datasets are presented. We focus on the study of functional relationships such as co-expression, protein-protein interactions, and regulatory interactions that are particularly amenable to modeling using the framework of networks. We envisage that many potential users of these analytical strategies may not be completely literate in programming languages and code adaptation, and for this reason, emphasis is given to tools' user-friendliness, including plugins for the widely adopted Cytoscape software, an open-source, cross-platform tool for network analysis, visualization, and data integration.

Meta-Analysis of HTLV-1-Infected Patients Identifies CD40LG and GBP2 as Markers of ATLL and HAM/TSP Clinical Status: Two Genes Beat as One.

Human T-lymphotropic virus 1 (HTLV-1) was the first recognized human retrovirus. Infection can lead to two main symptomatologies: adult T-cell lymphoma/leukemia (ATLL) and HTLV-1 associated myelopathy/tropical spastic paraparesis (HAM/TSP). Each manifestation is associated with distinct characteristics, as ATLL presents as a leukemia-like disease, while HAM/TSP presents as severe inflammation in the central nervous system, leading to paraparesis. Previous studies have identified molecules associated with disease development, e.g., the downregulation of Foxp3 in Treg cells was associated with increased risk of HAM/TSP. In addition, elevated levels of CXCL10, CXCL9, and Neopterin in cerebrospinal fluid also present increased risk. However, these molecules were only associated with specific patient groups or viral strains. Furthermore, the majority of studies did not jointly compare all clinical manifestations, and robust analysis entails the inclusion of both ATLL and HAM/TSP. The low numbers of samples also pose difficulties in conducting gene expression analysis to identify specific molecular relationships. To address these limitations and increase the power of manifestation-specific gene associations, meta-analysis was performed using publicly available gene expression data. The application of supervised learning techniques identified alterations in two genes observed to act in tandem as potential biomarkers: GBP2 was associated with HAM/TSP, and CD40LG with ATLL. Together, both molecules demonstrated high sample-classification accuracy (AUC values: 0.88 and 1.0, respectively). Next, other genes with expression correlated to these genes were identified, and we attempted to relate the enriched pathways identified with the characteristic of each clinical manifestation. The present findings contribute to knowledge surrounding viral progression and suggest a potentially powerful new tool for the molecular classification of HTLV-associated diseases.

Proteomic Analysis Reveals a Predominant NFE2L2 (NRF2) Signature in Canonical Pathway and Upstream Regulator Analysis of Leishmania-Infected Macrophages.

CBA mice macrophages (MØ) control infection by Leishmania major and are susceptive to Leishmania amazonensis, suggesting that both parasite species induce distinct responses that play important roles in infection outcome. To evaluate the MØ responses to infection arising from these two Leishmania species, a proteomic study using a Multidimensional Protein Identification Technology (MudPIT) approach with liquid chromatography tandem mass spectrometry (LC-MS/MS) was carried out on CBA mice bone-marrow MØ (BMMØ). Following SEQUEST analysis, which revealed 2,838 proteins detected in BMMØ, data mining approach found six proteins significantly associated with the tested conditions. To investigate their biological significance, enrichment analysis was performed using Ingenuity Pathway Analysis (IPA). A three steps IPA approach revealed 4 Canonical Pathways (CP) and 7 Upstream Transcriptional Factors (UTFs) strongly associated with the infection process. NRF2 signatures were present in both CPs and UTFs pathways. Proteins involved in iron metabolism, such as heme oxigenase 1 (HO-1) and ferritin besides sequestosome (SQSMT1 or p62) were found in the NRF2 CPs and the NRF2 UTFs. Differences in the involvement of iron metabolism pathway in Leishmania infection was revealed by the presence of HO-1 and ferritin. Noteworty, HO-1 was strongly associated with L. amazonensis infection, while ferritin was regulated by both species. As expected, higher HO-1 and p62 expressions were validated in L. amazonensis-infected BMMØ, in addition to decreased expression of ferritin and nitric oxide production. Moreover, BMMØ incubated with L. amazonensis LPG also expressed higher levels of HO-1 in comparison to those stimulated with L. major LPG. In addition, L. amazonensis-induced uptake of holoTf was higher than that induced by L. major in BMMØ, and holoTf was also detected at higher levels in vacuoles induced by L. amazonensis. Taken together, these findings indicate that NRF2 pathway activation and increased HO-1 production, together with higher levels of holoTf uptake, may promote permissiveness to L. amazonensis infection. In this context, differences in protein signatures triggered in the host by L. amazonensis and L. major infection could drive the outcomes in distinct clinical forms of leishmaniasis.

In Search of Biomarkers for Pathogenesis and Control of Leishmaniasis by Global Analyses of Leishmania-Infected Macrophages.

Leishmaniasis is a vector-borne, neglected tropical disease with a worldwide distribution that can present in a variety of clinical forms, depending on the parasite species and host genetic background. The pathogenesis of this disease remains far from being elucidated because the involvement of a complex immune response orchestrated by host cells significantly affects the clinical outcome. Among these cells, macrophages are the main host cells, produce cytokines and chemokines, thereby triggering events that contribute to the mediation of the host immune response and, subsequently, to the establishment of infection or, alternatively, disease control. There has been relatively limited commercial interest in developing new pharmaceutical compounds to treat leishmaniasis. Moreover, advances in the understanding of the underlying biology of Leishmania spp. have not translated into the development of effective new chemotherapeutic compounds. As a result, biomarkers as surrogate disease endpoints present several potential advantages to be used in the identification of targets capable of facilitating therapeutic interventions considered to ameliorate disease outcome. More recently, large-scale genomic and proteomic analyses have allowed the identification and characterization of the pathways involved in the infection process in both parasites and the host, and these analyses have been shown to be more effective than studying individual molecules to elucidate disease pathogenesis. RNA-seq and proteomics are large-scale approaches that characterize genes or proteins in a given cell line, tissue, or organism to provide a global and more integrated view of the myriad biological processes that occur within a cell than focusing on an individual gene or protein. Bioinformatics provides us with the means to computationally analyze and integrate the large volumes of data generated by high-throughput sequencing approaches. The integration of genomic expression and proteomic data offers a rich multi-dimensional analysis, despite the inherent technical and statistical challenges. We propose that these types of global analyses facilitate the identification, among a large number of genes and proteins, those that hold potential as biomarkers. The present review focuses on large-scale studies that have identified and evaluated relevant biomarkers in macrophages in response to Leishmania infection.

An integrative, multi-omics approach towards the prioritization of Klebsiella pneumoniae drug targets.

Klebsiella pneumoniae (Kp) is a globally disseminated opportunistic pathogen that can cause life-threatening infections. It has been found as the culprit of many infection outbreaks in hospital environments, being particularly aggressive towards newborns and adults under intensive care. Many Kp strains produce extended-spectrum ß-lactamases, enzymes that promote resistance against antibiotics used to fight these infections. The presence of other resistance determinants leading to multidrug-resistance also limit therapeutic options, and the use of 'last-resort' drugs, such as polymyxins, is not uncommon. The global emergence and spread of resistant strains underline the need for novel antimicrobials against Kp and related bacterial pathogens. To tackle this great challenge, we generated multiple layers of 'omics' data related to Kp and prioritized proteins that could serve as attractive targets for antimicrobial development. Genomics, transcriptomics, structuromic and metabolic information were integrated in order to prioritize candidate targets, and this data compendium is freely available as a web server. Twenty-nine proteins with desirable characteristics from a drug development perspective were shortlisted, which participate in important processes such as lipid synthesis, cofactor production, and core metabolism. Collectively, our results point towards novel targets for the control of Kp and related bacterial pathogens.

Follow up of a robust meta-signature to identify Zika virus infection in Aedes aegypti: another brick in the wall.

The mosquito Aedes aegypti is the main vector of several arthropod-borne diseases that have global impacts. In a previous meta-analysis, our group identified a vector gene set containing 110 genes strongly associated with infections of dengue, West Nile and yellow fever viruses. Of these 110 genes, four genes allowed a highly accurate classification of infected status. More recently, a new study of Ae. aegypti infected with Zika virus (ZIKV) was published, providing new data to investigate whether this "infection" gene set is also altered during a ZIKV infection. Our hypothesis is that the infection-associated signature may also serve as a proxy to classify the ZIKV infection in the vector. Raw data associated with the NCBI/BioProject were downloaded and re-analysed. A total of 18 paired-end replicates corresponding to three ZIKV-infected samples and three controls were included in this study. The nMDS technique with a logistic regression was used to obtain the probabilities of belonging to a given class. Thus, to compare both gene sets, we used the area under the curve and performed a comparison using the bootstrap method. Our meta-signature was able to separate the infected mosquitoes from the controls with good predictive power to classify the Zika-infected mosquitoes.

Integrated Analysis Reveals That miR-193b, miR-671, and TREM-1 Correlate With a Good Response to Treatment of Human Localized Cutaneous Leishmaniasis Caused by Leishmania braziliensis.

Localized cutaneous leishmaniasis (LCL) is a chronic disease characterized by ulcerated skin lesion(s) and uncontrolled inflammation. The mechanisms underlying the pathogenesis of LCL are not completely understood, and little is known about posttranscriptional regulation during LCL. MicroRNAs (miRNAs) are non-coding small RNAs that regulate gene expression and can be implicated in the pathogenesis of LCL. We investigated the involvement of miRNAs and their targets genes in human LCL using publicly available transcriptome data sets followed by ex vivo validation. Initial analysis highlighted that miRNA expression is altered during LCL, as patients clustered separately from controls. Joint analysis identified eight high confidence miRNAs that had altered expression (-1.5 ≤ fold change ≥ 1.5; p < 0.05) between cutaneous ulcers and uninfected skin. We found that the expression of miR-193b and miR-671 are greatly associated with their target genes, CD40 and TNFR, indicating the important role of these miRNAs in the expression of genes related to the inflammatory response observed in LCL. In addition, network analysis revealed that miR-193b, miR-671, and TREM1 correlate only in patients who show faster wound healing (up to 59 days) and not in patients who require longer cure times (more than 60 days). Given that these miRNAs are associated with control of inflammation and healing time, our findings reveal that they might influence the pathogenesis and prognosis of LCL.

Comparative Genomic Analysis of a Clinical Isolate of Klebsiella quasipneumoniae subsp. similipneumoniae, a KPC-2 and OKP-B-6 Beta-Lactamases Producer Harboring Two Drug-Resistance Plasmids from Southeast Brazil.

The aim of this study was to unravel the genetic determinants responsible for multidrug (including carbapenems) resistance and virulence in a clinical isolate of Klebsiella quasipneumoniae subsp. similipneumoniae by whole-genome sequencing and comparative analyses. Eighty-three clinical isolates initially identified as carbapenem-resistant K. pneumoniae were collected from nosocomial infections in southeast Brazil. After RAPD screening, the KPC-142 isolate, showing the most divergent DNA pattern, was selected for complete genome sequencing in an Illumina HiSeq 2500 instrument. Reads were assembled into scaffolds, gaps between scaffolds were resolved by in silico gap filling and extensive bioinformatics analyses were performed, using multiple comparative analysis tools and databases. Genome sequencing allowed to correct the classification of the KPC-142 isolate as K. quasipneumoniae subsp. similipneumoniae. To the best of our knowledge this is the first complete genome reported to date of a clinical isolate of this subspecies harboring both class A beta-lactamases KPC-2 and OKP-B-6 from South America. KPC-142 has one 5.2 Mbp chromosome (57.8% G+C) and two plasmids: 190 Kbp pKQPS142a (50.7% G+C) and 11 Kbp pKQPS142b (57.3% G+C). The 3 Kbp region in pKQPS142b containing the blaKPC-2 was found highly similar to that of pKp13d of K. pneumoniae Kp13 isolated in Southern Brazil in 2009, suggesting the horizontal transfer of this resistance gene between different species of Klebsiella. KPC-142 additionally harbors an integrative conjugative element ICEPm1 that could be involved in the mobilization of pKQPS142b and determinants of resistance to other classes of antimicrobials, including aminoglycoside and silver. We present the completely assembled genome sequence of a clinical isolate of K. quasipneumoniae subsp. similipneumoniae, a KPC-2 and OKP-B-6 beta-lactamases producer and discuss the most relevant genomic features of this important resistant pathogen in comparison to several strains belonging to K. quasipneumoniae subsp. similipneumoniae (phylogroup II-B), K. quasipneumoniae subsp. quasipneumoniae (phylogroup II-A), K. pneumoniae (phylogroup I), and K. variicola (phylogroup III). Our study contributes to the description of the characteristics of a novel K. quasipneumoniae subsp. similipneumoniae strain circulating in South America that currently represent a serious potential risk for nosocomial settings.

Follow up of a robust meta-signature to identify Zika virus infection in Aedes aegypti: another brick in the wall

The mosquito Aedes aegypti is the main vector of several arthropod-borne diseases that have global impacts. In a previous meta-analysis, our group identified a vector gene set containing 110 genes strongly associated with infections of dengue, West Nile and yellow fever viruses. Of these 110 genes, four genes allowed a highly accurate classification of infected status. More recently, a new study of Ae. aegypti infected with Zika virus (ZIKV) was published, providing new data to investigate whether this "infection" gene set is also altered during a ZIKV infection. Our hypothesis is that the infection-associated signature may also serve as a proxy to classify the ZIKV infection in the vector. Raw data associated with the NCBI/BioProject were downloaded and re-analysed. A total of 18 paired-end replicates corresponding to three ZIKV-infected samples and three controls were included in this study. The nMDS technique with a logistic regression was used to obtain the probabilities of belonging to a given class. Thus, to compare both gene sets, we used the area under the curve and performed a comparison using the bootstrap method. Our meta-signature was able to separate the infected mosquitoes from the controls with good predictive power to classify the Zika-infected mosquitoes.

Meta-Analysis of Aedes aegypti Expression Datasets: Comparing Virus Infection and Blood-Fed Transcriptomes to Identify Markers of Virus Presence.

The mosquito Aedes aegypti (L.) is vector of several arboviruses including dengue, yellow fever, chikungunya, and more recently zika. Previous transcriptomic studies have been performed to elucidate altered pathways in response to viral infection. However, the intrinsic coupling between alimentation and infection were unappreciated in these studies. Feeding is required for the initial mosquito contact with the virus and these events are highly dependent. Addressing this relationship, we reinterrogated datasets of virus-infected mosquitoes with two different diet schemes (fed and unfed mosquitoes), evaluating the metabolic cross-talk during both processes. We constructed coexpression networks with the differentially expressed genes of these comparison: virus-infected versus blood-fed mosquitoes and virus-infected versus unfed mosquitoes. Our analysis identified one module with 110 genes that correlated with infection status (representing ~0.7% of the A. aegypti genome). Furthermore, we performed a machine-learning approach and summarized the infection status using only four genes (AAEL012128, AAEL014210, AAEL002477, and AAEL005350). While three of the four genes were annotated as hypothetical proteins, AAEL012128 gene is a membrane amino acid transporter correlated with viral envelope binding. This gene alone is able to discriminate all infected samples and thus should have a key role to discriminate viral infection in the A. aegypti mosquito. Moreover, validation using external datasets found this gene as differentially expressed in four transcriptomic experiments. Therefore, these genes may serve as a proxy of viral infection in the mosquito and the others 106 identified genes provides a framework to future studies.

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 (Kp13pH), magnesium deprivation (Kp13Mg), high concentrations of calcium (Kp13Ca) and iron (Kp13Fe), and a control condition with PB (Kp13PolB). 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.

Analysis of the ergosterol biosynthesis pathway cloning, molecular characterization and phylogeny of lanosterol 14 α-demethylase (ERG11) gene of Moniliophthora perniciosa.

The phytopathogenic fungus Moniliophthora perniciosa (Stahel) Aime & Philips-Mora, causal agent of witches' broom disease of cocoa, causes countless damage to cocoa production in Brazil. Molecular studies have attempted to identify genes that play important roles in fungal survival and virulence. In this study, sequences deposited in the M. perniciosa Genome Sequencing Project database were analyzed to identify potential biological targets. For the first time, the ergosterol biosynthetic pathway in M. perniciosa was studied and the lanosterol 14α-demethylase gene (ERG11) that encodes the main enzyme of this pathway and is a target for fungicides was cloned, characterized molecularly and its phylogeny analyzed. ERG11 genomic DNA and cDNA were characterized and sequence analysis of the ERG11 protein identified highly conserved domains typical of this enzyme, such as SRS1, SRS4, EXXR and the heme-binding region (HBR). Comparison of the protein sequences and phylogenetic analysis revealed that the M. perniciosa enzyme was most closely related to that of Coprinopsis cinerea.

Comparative analysis of the complete genome of KPC-2-producing Klebsiella pneumoniae Kp13 reveals remarkable genome plasticity and a wide repertoire of virulence and resistance mechanisms.

BACKGROUND: Klebsiella pneumoniae is an important opportunistic pathogen associated with nosocomial and community-acquired infections. A wide repertoire of virulence and antimicrobial resistance genes is present in K. pneumoniae genomes, which can constitute extra challenges in the treatment of infections caused by some strains. K. pneumoniae Kp13 is a multidrug-resistant strain responsible for causing a large nosocomial outbreak in a teaching hospital located in Southern Brazil. Kp13 produces K. pneumoniae carbapenemase (KPC-2) but is unrelated to isolates belonging to ST 258 and ST 11, the main clusters associated with the worldwide dissemination of KPC-producing K. pneumoniae. In this report, we perform a genomic comparison between Kp13 and each of the following three K. pneumoniae genomes: MGH 78578, NTUH-K2044 and 342. RESULTS: We have completely determined the genome of K. pneumoniae Kp13, which comprises one chromosome (5.3 Mbp) and six plasmids (0.43 Mbp). Several virulence and resistance determinants were identified in strain Kp13. Specifically, we detected genes coding for six beta-lactamases (SHV-12, OXA-9, TEM-1, CTX-M-2, SHV-110 and KPC-2), eight adhesin-related gene clusters, including regions coding for types 1 (fim) and 3 (mrk) fimbrial adhesins. The rmtG plasmidial 16S rRNA methyltransferase gene was also detected, as well as efflux pumps belonging to five different families. Mutations upstream the OmpK35 porin-encoding gene were evidenced, possibly affecting its expression. SNPs analysis relative to the compared strains revealed 141 mutations falling within CDSs related to drug resistance which could also influence the Kp13 lifestyle. Finally, the genetic apparatus for synthesis of the yersiniabactin siderophore was identified within a plasticity region. Chromosomal architectural analysis allowed for the detection of 13 regions of difference in Kp13 relative to the compared strains. CONCLUSIONS: Our results indicate that the plasticity occurring at many hierarchical levels (from whole genomic segments to individual nucleotide bases) may play a role on the lifestyle of K. pneumoniae Kp13 and underlie the importance of whole-genome sequencing to study bacterial pathogens. The general chromosomal structure was somewhat conserved among the compared bacteria, and recombination events with consequent gain/loss of genomic segments appears to be driving the evolution of these strains.