Gene-Gene Interactions Dictate Ciprofloxacin Resistance in Pseudomonas aeruginosa and Facilitate Prediction of Resistance Phenotype from Genome Sequence Data.

Ciprofloxacin is one of the most widely used antibiotics for treating Pseudomonas aeruginosa infections. However, P. aeruginosa acquires mutations that confer ciprofloxacin resistance, making treatment more difficult. Resistance is multifactorial, with mutations in multiple genes influencing the resistance phenotype. However, the contributions of individual mutations and mutation combinations to the amounts of ciprofloxacin that P. aeruginosa can tolerate are not well understood. Engineering P. aeruginosa strain PAO1 to contain mutations in any one of the resistance-associated genes gyrA, nfxB, rnfC, parC, and parE showed that only gyrA mutations increased the MIC for ciprofloxacin. Mutations in parC and parE increased the MIC of a gyrA mutant, making the bacteria ciprofloxacin resistant. Mutations in nfxB and rnfC increased the MIC, conferring resistance, only if both were mutated in a gyrA background. Mutations in all of gyrA, nfxB, rnfC, and parC/E further increased the MIC. These findings reveal an epistatic network of gene-gene interactions in ciprofloxacin resistance. We used this information to predict ciprofloxacin resistance/susceptibility for 274 isolates of P. aeruginosa from their genome sequences. Antibiotic susceptibility profiles were predicted correctly for 84% of the isolates. The majority of isolates for which prediction was unsuccessful were ciprofloxacin resistant, demonstrating the involvement of additional as yet unidentified genes and mutations in resistance. Our data show that gene-gene interactions can play an important role in antibiotic resistance and can be successfully incorporated into models predicting resistance phenotype.

Genes Required for Free Phage Production are Essential for Pseudomonas aeruginosa Chronic Lung Infections.

The opportunistic pathogen Pseudomonas aeruginosa causes chronic lung infection in patients with cystic fibrosis. The Liverpool Epidemic Strain LESB58 is highly resistant to antibiotics, transmissible, and associated with increased morbidity and mortality. Its genome contains 6 prophages and 5 genomic islands. We constructed a polymerase chain reaction (PCR)-based signature-tagged mutagenesis library of 9216 LESB58 mutants and screened the mutants in a rat model of chronic lung infection. A total of 162 mutants were identified as defective for in vivo maintenance, with 11 signature-tagged mutagenesis mutants having insertions in prophage and genomic island genes. Many of these mutants showed both diminished virulence and reduced phage production. Transcription profiling by quantitative PCR and RNA-Seq suggested that disruption of these prophages had a widespread trans-acting effect on the transcriptome. This study demonstrates that temperate phages play a pivotal role in the establishment of infection through modulation of bacterial host gene expression.

Comparative genomics and biological characterization of sequential Pseudomonas aeruginosa isolates from persistent airways infection.

BACKGROUND: Pseudomonas aeruginosa establishes life-long chronic airway infections in cystic fibrosis (CF) patients. As the disease progresses, P. aeruginosa pathoadaptive variants are distinguished from the initially acquired strain. However, the genetic basis and the biology of host-bacteria interactions leading to a persistent lifestyle of P. aeruginosa are not understood. As a model system to study long term and persistent CF infections, the P. aeruginosa RP73, isolated 16.9 years after the onset of airways colonization from a CF patient, was investigated. Comparisons with strains RP1, isolated at the onset of the colonization, and clonal RP45, isolated 7 years before RP73 were carried out to better characterize genomic evolution of P. aeruginosa in the context of CF pathogenicity. RESULTS: Virulence assessments in disease animal model, genome sequencing and comparative genomics analysis were performed for clinical RP73, RP45, RP1 and prototype strains. In murine model, RP73 showed lower lethality and a remarkable capability of long-term persistence in chronic airways infection when compared to other strains. Pathological analysis of murine lungs confirmed advanced chronic pulmonary disease, inflammation and mucus secretory cells hyperplasia. Genomic analysis predicted twelve genomic islands in the RP73 genome, some of which distinguished RP73 from other prototype strains and corresponded to regions of genome plasticity. Further, comparative genomic analyses with sequential RP isolates showed signatures of pathoadaptive mutations in virulence factors potentially linked to the development of chronic infections in CF. CONCLUSIONS: The genome plasticity of P. aeruginosa particularly in the RP73 strain strongly indicated that these alterations may form the genetic basis defining host-bacteria interactions leading to a persistent lifestyle in human lungs.

A large decrease of cytosolic triosephosphate isomerase in transgenic potato roots affects the distribution of carbon in primary metabolism.

Triosephosphate isomerase (TPI, EC 5.3.1.1) catalyzes the interconversion of dihydroxyacetone-P and glyceraldehyde 3-P in the glycolytic pathway. A constitutively expressed antisense construct for cytosolic TPI was introduced into potato (Solanum tuberosum) using Agrobacterium rhizogenes to examine the metabolic effects of a reduction in cytosolic TPI in roots. We obtained a population of transgenic root clones displaying ~36 to 100 % of the TPI activity found in control clones carrying an empty binary vector. Ion exchange chromatography and immunoblot analysis showed that the antisense strategy significantly decreased the cytosolic TPI isoform, while levels of plastidial TPI activity remained apparently unaffected. Transgenic roots were characterized with respect to the activity of glycolytic enzymes, their metabolite contents and carbon fluxes. Metabolite profiling of sugars, organic acids, amino acids and lipids showed elevated levels of sucrose, glucose, fructose, fumarate, isocitrate, 4-aminobutyrate, alanine, glycine, aromatic amino acids and saturated long chain fatty acids in roots containing the lowest TPI activity. Labelings with (14)C-glucose, (14)C-sucrose and (14)C-acetate indicated that a reduction of cytosolic TPI activity in roots increased carbon metabolism through the pentose phosphate pathway, O(2) uptake and catabolism of sucrose to CO(2), and capacity for lipid synthesis. These results demonstrate that a large reduction of cytosolic TPI alters the distribution of carbon in plant primary metabolism.

Draft Genome Sequences of Two Clostridium botulinum Group II Strains Carrying Phage-Like Plasmids.

Clostridium botulinum is responsible for botulism, a potentially lethal foodborne intoxication. Here, we report the draft genome sequences of C. botulinum group II strains 202F (serotype F) and Hazen (serotype E). The genomes share many similarities, including multiple mobile genetic elements.

Transmission, adaptation and geographical spread of the Pseudomonas aeruginosa Liverpool epidemic strain.

The Liverpool epidemic strain (LES) is an important transmissible clonal lineage of Pseudomonas aeruginosa that chronically infects the lungs of people with cystic fibrosis (CF). Previous studies have focused on the genomics of the LES in a limited number of isolates, mostly from one CF centre in the UK, and from studies highlighting identification of the LES in Canada. Here we significantly extend the current LES genome database by genome sequencing 91 isolates from multiple CF centres across the UK, and we describe the comparative genomics of this large collection of LES isolates from the UK and Canada. Phylogenetic analysis revealed that the 145 LES genomes analysed formed a distinct clonal lineage when compared with the wider P. aeruginosa population. Notably, the isolates formed two clades: one associated with isolates from Canada, and the other associated with UK isolates. Further analysis of the UK LES isolates revealed clustering by clinic geography. Where isolates clustered closely together, the association was often supported by clinical data linking isolates or patients. When compared with the earliest known isolate, LESB58 (from 1988), many UK LES isolates shared common loss-of-function mutations, such as in genes gltR and fleR. Other loss-of-function mutations identified in previous studies as common adaptations during CF chronic lung infections were also identified in multiple LES isolates. Analysis of the LES accessory genome (including genomic islands and prophages) revealed variations in the carriage of large genomic regions, with some evidence for shared genomic island/prophage complement according to clinic location. Our study reveals divergence and adaptation during the spread of the LES, within the UK and between continents.

Candidate genes and adaptive radiation: insights from transcriptional adaptation to the limnetic niche among coregonine fishes (Coregonus spp., Salmonidae).

In the past 40 years, there has been increasing acceptance that variation in levels of gene expression represents a major source of evolutionary novelty. Gene expression divergence is therefore likely to be involved in the emergence of incipient species, namely, in a context of adaptive radiation. In the lake whitefish species complex (Coregonus clupeaformis), previous microarray experiments have led to the identification of candidate genes potentially implicated in the parallel evolution of the limnetic dwarf lake whitefish, which is highly distinct from the benthic normal lake whitefish in life history, morphology, metabolism, and behavior, and yet diverged from it only approximately 15,000 years before present. The aim of the present study was to address transcriptional divergence for six candidate genes among lake whitefish and European whitefish (Coregonus lavaretus) species pairs, as well as lake cisco (Coregonus artedi) and vendace (Coregonus albula). The main goal was to test the hypothesis that parallel phenotypic adaptation toward the use of the limnetic niche in coregonine fishes is accompanied by parallelism in candidate gene transcription as measured by quantitative real-time polymerase chain reaction. Results obtained for three candidate genes, whereby parallelism in expression was observed across all whitefish species pairs, provide strong support for the hypothesis that divergent natural selection plays an important role in the adaptive radiation of whitefish species. However, this parallelism in expression did not extend to cisco and vendace, thereby infirming transcriptional convergence between limnetic whitefish species and their limnetic congeners for these genes. As recently proposed (Lynch 2007a. The evolution of genetic networks by non-adaptive processes. Nat Rev Genet. 8:803-813), these results may suggest that convergent phenotypic evolution can result from nonadaptive shaping of genome architecture in independently evolved coregonine lineages.

The transcriptomics of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis spp., Salmonidae) divergence as revealed by next-generation sequencing.

Gene expression divergence is one of the mechanisms thought to be involved in the emergence of incipient species. Next-generation sequencing has become an extremely valuable tool for the study of this process by allowing whole transcriptome sequencing, or RNA-Seq. We have conducted a 454 GS-FLX pyrosequencing experiment to refine our understanding of adaptive divergence between dwarf and normal lake whitefish species (Coregonus clupeaformis spp.). The objectives were to: (i) investigate transcriptomic divergence as measured by liver RNA-Seq; (ii) test the correlation between divergence in expression and sequence polymorphism; and (iii) investigate the extent of allelic imbalance. We also compared the results of RNA-seq with those of a previous microarray study performed on the same fish. Following de novo assembly, results showed that normal whitefish overexpressed more contigs associated with protein synthesis while dwarf fish overexpressed more contigs related to energy metabolism, immunity and DNA replication and repair. Moreover, 63 SNPs showed significant allelic imbalance, and this phenomenon prevailed in the recently diverged dwarf whitefish. Results also showed an absence of correlation between gene expression divergence as measured by RNA-Seq and either polymorphism rate or sequence divergence between normal and dwarf whitefish. This study reiterates an important role for gene expression divergence, and provides evidence for allele-specific expression divergence as well as evolutionary decoupling of regulatory and coding sequences in the adaptive divergence of normal and dwarf whitefish. It also demonstrates how next-generation sequencing can lead to a more comprehensive understanding of transcriptomic divergence in a young species pair.

The Salmonella enterica Plasmidome as a Reservoir of Antibiotic Resistance.

The emergence of multidrug-resistant bacterial strains worldwide has become a serious problem for public health over recent decades. The increase in antimicrobial resistance has been expanding via plasmids as mobile genetic elements encoding antimicrobial resistance (AMR) genes that are transferred vertically and horizontally. This study focuses on Salmonella enterica, one of the leading foodborne pathogens in industrialized countries. S. enterica is known to carry several plasmids involved not only in virulence but also in AMR. In the current paper, we present an integrated strategy to detect plasmid scaffolds in whole genome sequencing (WGS) assemblies. We developed a two-step procedure to predict plasmids based on i) the presence of essential elements for plasmid replication and mobility, as well as ii) sequence similarity to a reference plasmid. Next, to confirm the accuracy of the prediction in 1750 S. enterica short-read sequencing data, we combined Oxford Nanopore MinION long-read sequencing with Illumina MiSeq short-read sequencing in hybrid assemblies for 84 isolates to evaluate the proportion of plasmid that has been detected. At least one scaffold with an origin of replication (ORI) was predicted in 61.3% of the Salmonella isolates tested. The results indicated that IncFII and IncI1 ORIs were distributed in many S. enterica serotypes and were the most prevalent AMR genes carrier, whereas IncHI2A/IncHI2 and IncA/C2 were more serotype restricted but bore several AMR genes. Comparison between hybrid and short-read assemblies revealed that 81.1% of plasmids were found in the short-read sequencing using our pipeline. Through this process, we established that plasmids are prevalent in S. enterica and we also substantially expand the AMR genes in the resistome of this species.

The phenomics and expression quantitative trait locus mapping of brain transcriptomes regulating adaptive divergence in lake whitefish species pairs (Coregonus sp.).

We used microarrays and a previously established linkage map to localize the genetic determinants of brain gene expression for a backcross family of lake whitefish species pairs (Coregonus sp.). Our goals were to elucidate the genomic distribution and sex specificity of brain expression QTL (eQTL) and to determine the extent to which genes controlling transcriptional variation may underlie adaptive divergence in the recently evolved dwarf (limnetic) and normal (benthic) whitefish. We observed a sex bias in transcriptional genetic architecture, with more eQTL observed in males, as well as divergence in genome location of eQTL between the sexes. Hotspots of nonrandom aggregations of up to 32 eQTL in one location were observed. We identified candidate genes for species pair divergence involved with energetic metabolism, protein synthesis, and neural development on the basis of colocalization of eQTL for these genes with eight previously identified adaptive phenotypic QTL and four previously identified outlier loci from a genome scan in natural populations. Eighty-eight percent of eQTL-phenotypic QTL colocalization involved growth rate and condition factor QTL, two traits central to adaptive divergence between whitefish species pairs. Hotspots colocalized with phenotypic QTL in several cases, revealing possible locations where master regulatory genes, such as a zinc-finger protein in one case, control gene expression directly related to adaptive phenotypic divergence. We observed little evidence of colocalization of brain eQTL with behavioral QTL, which provides insight into the genes identified by behavioral QTL studies. These results extend to the transcriptome level previous work illustrating that selection has shaped recent parallel divergence between dwarf and normal lake whitefish species pairs and that metabolic, more than morphological, differences appear to play a key role in this divergence.

Major Release of 161 Whole-Genome Sequences from the International Pseudomonas Consortium Database.

Pseudomonas aeruginosa is an environmental bacterium and opportunistic pathogen. Here, we present draft genome sequences for 161 isolates from diverse clinical and environmental sources. This set of genome sequences complements other major public data releases from the International Pseudomonas Consortium Database.

Genomics of antibiotic-resistance prediction in Pseudomonas aeruginosa.

Antibiotic resistance is a worldwide health issue spreading quickly among human and animal pathogens, as well as environmental bacteria. Misuse of antibiotics has an impact on the selection of resistant bacteria, thus contributing to an increase in the occurrence of resistant genotypes that emerge via spontaneous mutation or are acquired by horizontal gene transfer. There is a specific and urgent need not only to detect antimicrobial resistance but also to predict antibiotic resistance in silico. We now have the capability to sequence hundreds of bacterial genomes per week, including assembly and annotation. Novel and forthcoming bioinformatics tools can predict the resistome and the mobilome with a level of sophistication not previously possible. Coupled with bacterial strain collections and databases containing strain metadata, prediction of antibiotic resistance and the potential for virulence are moving rapidly toward a novel approach in molecular epidemiology. Here, we present a model system in antibiotic-resistance prediction, along with its promises and limitations. As it is commonly multidrug resistant, Pseudomonas aeruginosa causes infections that are often difficult to eradicate. We review novel approaches for genotype prediction of antibiotic resistance. We discuss the generation of microbial sequence data for real-time patient management and the prediction of antimicrobial resistance.

The Pseudomonas aeruginosa Pan-Genome Provides New Insights on Its Population Structure, Horizontal Gene Transfer, and Pathogenicity.

The huge increase in the availability of bacterial genomes led us to a point in which we can investigate and query pan-genomes, for example, the full set of genes of a given bacterial species or clade. Here, we used a data set of 1,311 high-quality genomes from the human pathogen Pseudomonas aeruginosa, 619 of which were newly sequenced, to show that a pan-genomic approach can greatly refine the population structure of bacterial species, provide new insights to define species boundaries, and generate hypotheses on the evolution of pathogenicity. The 665-gene P. aeruginosa core genome presented here, which constitutes only 1% of the entire pan-genome, is the first to be in the same order of magnitude as the minimal bacterial genome and represents a conservative estimate of the actual core genome. Moreover, the phylogeny based on this core genome provides strong evidence for a five-group population structure that includes two previously undescribed groups of isolates. Comparative genomics focusing on antimicrobial resistance and virulence genes showed that variation among isolates was partly linked to this population structure. Finally, we hypothesized that horizontal gene transfer had an important role in this respect, and found a total of 3,010 putative complete and fragmented plasmids, 5% and 12% of which contained resistance or virulence genes, respectively. This work provides data and strategies to study the evolutionary trajectories of resistance and virulence in P. aeruginosa.

Reservoirs of resistance: polymyxin resistance in veterinary-associated companion animal isolates of Pseudomonas aeruginosa.

BACKGROUND: Pseudomonas aeruginosa is an opportunistic pathogen and a major cause of infections. Widespread resistance in human infections are increasing the use of last resort antimicrobials such as polymyxins. However, these have been used for decades in veterinary medicine. Companion animals are an understudied source of antimicrobial resistant P. aeruginosa isolates. This study evaluated the susceptibility of P. aeruginosa veterinary isolates to polymyxins to determine whether the veterinary niche represents a potential reservoir of resistance genes for pathogenic bacteria in both animals and humans. METHODS AND RESULTS: Clinical P. aeruginosa isolates (n=24) from UK companion animals were compared for antimicrobial susceptibility to a panel of human-associated isolates (n=37). Minimum inhibitory concentration (MIC) values for polymyxin B and colistin in the companion animals was significantly higher than in human isolates (P=0.033 and P=0.013, respectively). Genotyping revealed that the veterinary isolates were spread throughout the P. aeruginosa population, with shared array types from human infections such as keratitis and respiratory infections, suggesting the potential for zoonotic transmission. Whole genome sequencing revealed mutations in genes associated with polymyxin resistance and other antimicrobial resistance-related genes. CONCLUSION: The high levels of resistance to polymyxin shown here, along with genetic similarities between some human and animal isolates, together suggest a need for sustained surveillance of this veterinary niche as a potential reservoir for resistant, clinically relevant bacteria in both animals and humans.

Genomic characterisation of an international Pseudomonas aeruginosa reference panel indicates that the two major groups draw upon distinct mobile gene pools.

Pseudomonas aeruginosa is an important opportunistic pathogen, especially in the context of infections of cystic fibrosis (CF). In order to facilitate coordinated study of this pathogen, an international reference panel of P. aeruginosa isolates was assembled. Here we report the genome sequencing and analysis of 33 of these isolates and 7 reference genomes to further characterise this panel. Core genome single nucleotide variant phylogeny demonstrated that the panel strains are widely distributed amongst the P. aeruginosa population. Common loss-of-function mutations reported as adaptive during CF (such as in mucA and mexA) were identified amongst isolates from chronic respiratory infections. From the 40 strains analysed, 37 unique resistomes were predicted, based on the Resistance Gene Identifier method using the Comprehensive Antibiotic Resistance Database. Notably, hierarchical clustering and phylogenetic reconstructions based on the presence/absence of genomic islands (GIs), prophages and other regions of genome plasticity (RGPs) supported the subdivision of P. aeruginosa into two main groups. This is the largest, most diverse analysis of GIs and associated RGPs to date, and the results suggest that, at least at the largest clade grouping level (group 1 vs group 2), each group may be drawing upon distinct mobile gene pools.

Salmonella enterica Prophage Sequence Profiles Reflect Genome Diversity and Can Be Used for High Discrimination Subtyping.

Non-typhoidal Salmonella is a leading cause of foodborne illness worldwide. Prompt and accurate identification of the sources of Salmonella responsible for disease outbreaks is crucial to minimize infections and eliminate ongoing sources of contamination. Current subtyping tools including single nucleotide polymorphism (SNP) typing may be inadequate, in some instances, to provide the required discrimination among epidemiologically unrelated Salmonella strains. Prophage genes represent the majority of the accessory genes in bacteria genomes and have potential to be used as high discrimination markers in Salmonella. In this study, the prophage sequence diversity in different Salmonella serovars and genetically related strains was investigated. Using whole genome sequences of 1,760 isolates of S. enterica representing 151 Salmonella serovars and 66 closely related bacteria, prophage sequences were identified from assembled contigs using PHASTER. We detected 154 different prophages in S. enterica genomes. Prophage sequences were highly variable among S. enterica serovars with a median ± interquartile range (IQR) of 5 ± 3 prophage regions per genome. While some prophage sequences were highly conserved among the strains of specific serovars, few regions were lineage specific. Therefore, strains belonging to each serovar could be clustered separately based on their prophage content. Analysis of S. Enteritidis isolates from seven outbreaks generated distinct prophage profiles for each outbreak. Taken altogether, the diversity of the prophage sequences correlates with genome diversity. Prophage repertoires provide an additional marker for differentiating S. enterica subtypes during foodborne outbreaks.

AnCo3, a New Member of the Emerging Family of Phage-Like Plasmids.

A phage-like plasmid isolated from a clinical isolate of Salmonella enterica serovar Derby has strong nucleotide sequence identity to the phage-like plasmids pSTM_phi isolated from Salmonella enterica serovar Typhimurium L495, AnCo1 and AnCo2 from Escherichia coli 243 and Escherichia coli 244, and the virulent Salmonella-specific SSU5 bacteriophage.

Whole-Genome Sequencing of Lactobacillus Species from Two Commercial Probiotic Products.

Eight Lactobacillus strains, each intrinsically resistant to an antibiotic, were isolated from two commercial probiotic products. Whole-genome sequencing identified two efflux transporters, a multidrug and extrusion protein (MATE) efflux transporter, and LmrCD, which may contribute to their intrinsic antibiotic resistance and may therefore facilitate their survival in the intestinal microbiota following antibiotic therapy.

Complete Genome Sequences of Two Phage-Like Plasmids Carrying the CTX-M-15 Extended-Spectrum ß-Lactamase Gene.

Two similar phage-like plasmids carrying CTX-M-15 resistance cassettes were identified from two environmental Escherichia coli isolates. They demonstrate strong nucleotide sequence identity to the phage-like plasmid pECOH89 and Salmonella bacteriophage SSU5.