Pathogenomic analysis and characterization of Pasteurella multocida strains recovered from human infections

Pasteurella multocida is an upper respiratory tract commensal in several mammal and bird species but can also cause severe disease in humans and in production animals such as poultry, cattle, and pigs. In this study, we performed whole-genome sequencing of P. multocida isolates recovered from a range of human infections, from the mouths of cats, and from wounds on dogs. Together with publicly available P. multocida genome sequences, we performed phylogenetic and comparative genomic analyses. While isolates from cats and dogs were spread across the phylogenetic tree, human infections were caused almost exclusively by subsp. septica strains. Most of the human isolates were capsule type A and LPS type L1 and L3; however, some strains lacked a capsule biosynthesis locus, and some strains contained a novel LPS outer-core locus, distinct from the eight LPS loci that can currently be identified using an LPS multiplex PCR. In addition, the P. multocida strains isolated from human infections contained novel mobile genetic elements. We compiled a curated database of known P. multocida virulence factor and antibiotic resistance genes (PastyVRDB) allowing for detailed characterization of isolates. The majority of human P. multocida isolates encoded a reduced range of iron receptors and contained only one filamentous hemagglutinin gene. Finally, gene-trait analysis identified a putative L-fucose uptake and utilization pathway that was over-represented in subsp. septica strains and may represent a novel host predilection mechanism in this subspecies. Together, these analyses have identified pathogenic mechanisms likely important for P. multocida zoonotic infections.IMPORTANCEPasteurella multocida can cause serious infections in humans, including skin and wound infections, pneumonia, peritonitis, meningitis, and bacteraemia. Cats and dogs are known vectors of human pasteurellosis, transmitting P. multocida via bite wounds or contact with animal saliva. The mechanisms that underpin P. multocida human predilection and pathogenesis are poorly understood. With increasing identification of antibiotic-resistant P. multocida strains, understanding these mechanisms is vital for developing novel treatments and control strategies to combat P. multocida human infection. Here, we show that a narrow range of P. multocida strains cause disease in humans, while cats and dogs, common vectors for zoonotic infections, can harbor a wide range of P. multocida strains. We also present a curated P. multocida-specific database, allowing quick and detailed characterization of newly sequenced P. multocida isolates.

The evolutionary mechanism of non-carbapenemase carbapenem-resistant phenotypes in Klebsiella spp

Antibiotic resistance is driven by selection, but the degree to which a bacterial strain's evolutionary history shapes the mechanism and strength of resistance remains an open question. Here, we reconstruct the genetic and evolutionary mechanisms of carbapenem resistance in a clinical isolate of Klebsiella quasipneumoniae. A combination of short- and long-read sequencing, machine learning, and genetic and enzymatic analyses established that this carbapenem-resistant strain carries no carbapenemase-encoding genes. Genetic reconstruction of the resistance phenotype confirmed that two distinct genetic loci are necessary in order for the strain to acquire carbapenem resistance. Experimental evolution of the carbapenem-resistant strains in growth conditions without the antibiotic revealed that both loci confer a significant cost and are readily lost by de novo mutations resulting in the rapid evolution of a carbapenem-sensitive phenotype. To explain how carbapenem resistance evolves via multiple, low-fitness single-locus intermediates, we hypothesised that one of these loci had previously conferred adaptation to another antibiotic. Fitness assays in a range of drug concentrations show how selection in the antibiotic ceftazidime can select for one gene (bla_DHA-1) potentiating the evolution of carbapenem resistance by a single mutation in a second gene (ompK36). These results show how a patient's treatment history might shape the evolution of antibiotic resistance and could explain the genetic basis of carbapenem-resistance found in many enteric-pathogens.

A Structurally Characterized Staphylococcus aureus Evolutionary Escape Route from Treatment with the Antibiotic Linezolid

Methicillin-resistant Staphylococcus aureus (MRSA) is a bacterial pathogen that presents great health concerns. Treatment requires the use of last-line antibiotics, such as members of the oxazolidinone family, of which linezolid is the first member to see regular use in the clinic. Here, we report a short time scale selection experiment in which strains of MRSA were subjected to linezolid treatment. Clonal isolates which had evolved a linezolid-resistant phenotype were characterized by whole-genome sequencing. Linezolid-resistant mutants were identified which had accumulated mutations in the ribosomal protein uL3. Multiple clones which had two mutations in uL3 exhibited resistance to linezolid, 2-fold higher than the clinical breakpoint. Ribosomes from this strain were isolated and subjected to single-particle cryo-electron microscopic analysis and compared to the ribosomes from the parent strain. We found that the mutations in uL3 lead to a rearrangement of a loop that makes contact with Helix 90, propagating a structural change over 15 Å away. This distal change swings nucleotide U2504 into the binding site of the antibiotic, causing linezolid resistance. IMPORTANCE Antibiotic resistance poses a critical problem to human health and decreases the utility of these lifesaving drugs. Of particular concern is the "superbug" methicillin-resistant Staphylococcus aureus (MRSA), for which treatment of infection requires the use of last-line antibiotics, including linezolid. In this paper, we characterize the atomic rearrangements which the ribosome, the target of linezolid, undergoes during its evolutionary journey toward becoming drug resistant. Using cryo-electron microscopy, we describe a particular molecular mechanism which MRSA uses to become resistant to linezolid.

Forensic genomics of a novel Klebsiella quasipneumoniae type from a neonatal intensive care unit in China reveals patterns of colonization, evolution and epidemiology

During March 2017, a neonatal patient with severe diarrhoea subsequently developed septicaemia and died, with Klebsiella isolated as the causative microorganism. In keeping with infection control protocols, the coincident illness of an attending staff member and three other neonates with Klebsiella infection triggered an outbreak response, leading to microbiological assessment of isolates collected from the staff member and all 21 co-housed neonates. Multilocus sequence typing and genomic sequencing identified that the isolates from the 21 neonates were of a new Klebsiella sequence type, ST2727, and taxonomically belonged to K. quasipneumoniae subsp. similipneumoniae (formerly referred to as KpIIB). Genomic characterization showed that the isolated ST2727 strains had diverged from other K. quasipneumoniae subsp. similipneumoniae strains at least 90 years ago, whereas the neonatal samples were highly similar with a genomic divergence of 3.6 months. There was no relationship to the Klebsiella isolate from the staff member. This demonstrates that no transmission occurred from staff to patient or between patients. Rather, the data suggest that ST2727 colonized each neonate from a common hospital source. Sequence-based analysis of the genomes revealed several genes for antimicrobial resistance and some virulence features, but suggest that ST2727 is neither extremely-drug resistant nor hypervirulent. Our results highlight the clinical significance and genomic properties of ST2727 and urge genome-based measures be implemented for diagnostics and surveillance within hospital environments. Additionally, the present study demonstrates the need to scale the power of genomic analysis in retrospective studies where relatively few samples are available.

First description of antimicrobial resistance in carbapenem-susceptible Klebsiella pneumoniae after imipenem treatment, driven by outer membrane remodeling

BACKGROUND

The emergence of carbapenem-resistant Klebsiella pneumoniae (CRKP) poses a looming threat to human health. Although there are numerous studies regarding porin alteration in association with the production of ESBLs and/or AmpC β-lactamase, a systematic study on the treatment-emergence of porins alteration in antibiotic resistance does not yet exist. The aim of this study was to investigate the underlying mechanism of resistance of K. pneumoniae during carbapenem treatment.

RESULTS

Here, we report three strains (FK-2624, FK-2723 and FK-2820) isolated from one patient before and after imipenem treatment during hospitalization. Antibiotic susceptibility testing indicated that that the first isolate, FK-2624, was susceptible to almost all tested antimicrobials, being resistant only to fosfomycin. The subsequent isolates FK-2723 and FK-2820 were multidrug resistant (MDR). After imipenem therapy, FK-2820 was found to be carbapenem-resistant. PCR and Genome Sequencing analysis indicated that oqxA, and fosA5, were identified in all three strains. In addition, FK-2624 also harbored bla_SHV-187 and bla_TEM-116. The bla_SHV-187 and bla_TEM-116 genes were not detected in FK-2723 and FK-2820. bla_DHA-1, qnrB4, aac (6')-IIc, and bla_SHV-12, EreA2, CatA2, SulI, and tetD, were identified in both FK-2723 and FK-2820. Moreover, the genes bla_DHA-1, qnrB4, aac (6')-IIc were co-harbored on a plasmid. Of the virulence factors found in this study, ybtA, ICEKp6, mrkD, entB, iroN, rmpA2-6, wzi16 and capsular serotype K57 were found in the three isolates. The results of pairwise comparisons, multi-locus sequencing typing (MLST) and pulsed-field gel electrophoresis (PFGE) revealed high homology among the isolates. Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) results showed that isolate FK-2820 lacked OmpK36, with genome sequence data validating that there was a premature stop codon in the ompK36 gene and real-time RT-PCR suggesting high turnover of the ompK36 non-sense transcript in FK-2820, with the steady-state mRNA level 0.007 relative to the initial isolate.

CONCLUSION

This study in China highlight that the alteration of outer membrane porins due to the 14-day use of imipenem play a potential role in leading to clinical presentation of carbapenem-resistance. This is the first description of increased resistance developing from a carbapenem-susceptible K. pneumoniae with imipenem treatment driven by outer membrane remodeling.

A genome-wide scan for correlated mutations detects macromolecular and chromatin interactions in Arabidopsis thaliana

The concept of exploiting correlated mutations has been introduced and applied successfully to identify interactions within and between biological macromolecules. Its rationale lies in the preservation of physical interactions via compensatory mutations. With the massive increase of available sequence information, approaches based on correlated mutations have regained considerable attention. We analyzed a set of 10 707 430 single nucleotide polymorphisms detected in 1135 accessions of the plant Arabidopsis thaliana. To measure their covariance and to reveal the global genome-wide sequence correlation structure of the Arabidopsis genome, the adjusted mutual information has been estimated for each possible pair of polymorphic sites. We developed a series of filtering steps to account for genetic linkage and lineage relations between Arabidopsis accessions, as well as transitive covariance as possible confounding factors. We show that upon appropriate filtering, correlated mutations prove indeed informative with regard to molecular interactions, and furthermore, appear to reflect on chromosomal interactions. Our study demonstrates that the concept of correlated mutations can also be applied successfully to within-species sequence variation and establishes a promising approach to help unravel the complex molecular interactions in A. thaliana and other species with broad sequence information.

A genetic map of cassava (Manihot esculenta Crantz) with integrated physical mapping of immunity-related genes

BACKGROUND

Cassava, Manihot esculenta Crantz, is one of the most important crops world-wide representing the staple security for more than one billion of people. The development of dense genetic and physical maps, as the basis for implementing genetic and molecular approaches to accelerate the rate of genetic gains in breeding program represents a significant challenge. A reference genome sequence for cassava has been made recently available and community efforts are underway for improving its quality. Cassava is threatened by several pathogens, but the mechanisms of defense are far from being understood. Besides, there has been a lack of information about the number of genes related to immunity as well as their distribution and genomic organization in the cassava genome.

RESULTS

A high dense genetic map of cassava containing 2,141 SNPs has been constructed. Eighteen linkage groups were resolved with an overall size of 2,571 cM and an average distance of 1.26 cM between markers. More than half of mapped SNPs (57.4%) are located in coding sequences. Physical mapping of scaffolds of cassava whole genome sequence draft using the mapped markers as anchors resulted in the orientation of 687 scaffolds covering 45.6% of the genome. One hundred eighty nine new scaffolds are anchored to the genetic cassava map leading to an extension of the present cassava physical map with 30.7 Mb. Comparative analysis using anchor markers showed strong co-linearity to previously reported cassava genetic and physical maps. In silico based searching for conserved domains allowed the annotation of a repertory of 1,061 cassava genes coding for immunity-related proteins (IRPs). Based on physical map of the corresponding sequencing scaffolds, unambiguous genetic localization was possible for 569 IRPs.

CONCLUSIONS

This is the first study reported so far of an integrated high density genetic map using SNPs with integrated genetic and physical localization of newly annotated immunity related genes in cassava. These data build a solid basis for future studies to map and associate markers with single loci or quantitative trait loci for agronomical important traits. The enrichment of the physical map with novel scaffolds is in line with the efforts of the cassava genome sequencing consortium.