Colin F. Davenport

Genome Diversity of Pseudomonas aeruginosa PAO1 Laboratory Strains

Pseudomonas aeruginosa PAO1 is the most commonly used strain for research on this ubiquitous and metabolically versatile opportunistic pathogen. Strain PAO1, a derivative of the original Australian PAO isolate, has been distributed worldwide to laboratories and strain collections. Over decades discordant phenotypes of PAO1 sublines have emerged. Taking the existing PAO1-UW genome sequence (named after the University of Washington, which led the sequencing project) as a blueprint, the genome sequences of reference strains MPAO1 and PAO1-DSM (stored at the German Collection for Microorganisms and Cell Cultures [DSMZ]) were resolved by physical mapping and deep short read sequencing-by-synthesis. MPAO1 has been the source of near-saturation libraries of transposon insertion mutants, and PAO1-DSM is identical in its SpeI-DpnI restriction map with the original isolate. The major genomic differences of MPAO1 and PAO1-DSM in comparison to PAO1-UW are the lack of a large inversion, a duplication of a mobile 12-kb prophage region carrying a distinct integrase and protein phosphatases or kinases, deletions of 3 to 1,006 bp in size, and at least 39 single-nucleotide substitutions, 17 of which affect protein sequences. The PAO1 sublines differed in their ability to cope with nutrient limitation and their virulence in an acute murine airway infection model. Subline PAO1-DSM outnumbered the two other sublines in late stationary growth phase. In conclusion, P. aeruginosa PAO1 shows an ongoing microevolution of genotype and phenotype that jeopardizes the reproducibility of research. High-throughput genome resequencing will resolve more cases and could become a proper quality control for strain collections.

Microevolution of the major common Pseudomonas aeruginosa clones C and PA14 in cystic fibrosis lungs.

Clones C and PA14 are the worldwide most abundant clonal complexes in the Pseudomonas aeruginosa population. The microevolution of clones C and PA14 was investigated in serial cystic fibrosis (CF) airway isolates collected over 20 years since the onset of colonization. Intraclonal evolution in CF lungs was resolved by genome sequencing of first, intermediate and late isolates and subsequent multimarker SNP genotyping of the whole strain panel. Mapping of sequence reads onto the P. aeruginosa PA14 reference genome unravelled an intraclonal and interclonal sequence diversity of 0.0035% and 0.68% respectively. Clone PA14 diversified into three branches in the patients lungs, and the PA14 population acquired 15 nucleotide substitutions and a large deletion during the observation period. The clone C genome remained invariant during the first 3 years in CF lungs; however, 15 years later 947 transitions and 12 transversions were detected in a clone C mutL mutant strain. Key mutations occurred in retS, RNA polymerase, multidrug transporter, virulence and denitrification genes. Late clone C and PA14 persistors in the CF lungs were compromised in growth and cytotoxicity, but their mutation frequency was normal even in mutL mutant clades.

Pseudomonas aeruginosa Genomic Structure and Diversity.

The Pseudomonas aeruginosa genome (G + C content 65–67%, size 5.5–7 Mbp) is made up of a single circular chromosome and a variable number of plasmids. Sequencing of complete genomes or blocks of the accessory genome has revealed that the genome encodes a large repertoire of transporters, transcriptional regulators, and two-component regulatory systems which reflects its metabolic diversity to utilize a broad range of nutrients. The conserved core component of the genome is largely collinear among P. aeruginosa strains and exhibits an interclonal sequence diversity of 0.5–0.7%. Only a few loci of the core genome are subject to diversifying selection. Genome diversity is mainly caused by accessory DNA elements located in 79 regions of genome plasticity that are scattered around the genome and show an anomalous usage of mono- to tetradecanucleotides. Genomic islands of the pKLC102/PAGI-2 family that integrate into tRNALys or tRNAGly genes represent hotspots of inter- and intraclonal genomic diversity. The individual islands differ in their repertoire of metabolic genes that make a large contribution to the pangenome. In order to unravel intraclonal diversity of P. aeruginosa, the genomes of two members of the PA14 clonal complex from diverse habitats and geographic origin were compared. The genome sequences differed by less than 0.01% from each other. One hundred ninety-eight of the 231 single nucleotide substitutions (SNPs) were non-randomly distributed in the genome. Non-synonymous SNPs were mainly found in an integrated Pf1-like phage and in genes involved in transcriptional regulation, membrane and extracellular constituents, transport, and secretion. In summary, P. aeruginosa is endowed with a highly conserved core genome of low sequence diversity and a highly variable accessory genome that communicates with other pseudomonads and genera via horizontal gene transfer.

Interclonal gradient of virulence in the Pseudomonas aeruginosa pangenome from disease and environment

The population genomics of Pseudomonas aeruginosa was analysed by genome sequencing of representative strains of the 15 most frequent clonal complexes in the P. aeruginosa population and of the five most common clones from the environment of which so far no isolate from a human infection has been detected. Gene annotation identified 5892-7187 open reading frame (ORFs; median 6381 ORFs) in the 20 6.4-7.4 Mbp large genomes. The P. aeruginosa pangenome consists of a conserved core of at least 4000 genes, a combinatorial accessory genome of a further 10 000 genes and 30 000 or more rare genes that are present in only a few strains or clonal complexes. Whole genome comparisons of single nucleotide polymorphism synteny indicated unrestricted gene flow between clonal complexes by recombination. Using standardized acute lettuce, Galleria mellonella and murine airway infection models the full spectrum of possible host responses to P. aeruginosa was observed with the 20 strains ranging from unimpaired health following infection to 100% lethality. Genome comparisons indicate that the differential genetic repertoire of clones maintains a habitat-independent gradient of virulence in the P. aeruginosa population.

The SeqWord Genome Browser: an online tool for the identification and visualization of atypical regions of bacterial genomes through oligonucleotide usage

BackgroundData mining in large DNA sequences is a major challenge in microbial genomics and bioinformatics. Oligonucleotide usage (OU) patterns provide a wealth of information for large scale sequence analysis and visualization. The purpose of this research was to make OU statistical analysis available as a novel web-based tool for functional genomics and annotation. The tool is also available as a downloadable package.ResultsThe SeqWord Genome Browser (SWGB) was developed to visualize the natural compositional variation of DNA sequences. The applet is also used for identification of divergent genomic regions both in annotated sequences of bacterial chromosomes, plasmids, phages and viruses, and in raw DNA sequences prior to annotation by comparing local and global OU patterns. The applet allows fast and reliable identification of clusters of horizontally transferred genomic islands, large multi-domain genes and genes for ribosomal RNA. Within the majority of genomic fragments (also termed genomic core sequence), regions enriched with housekeeping genes, ribosomal proteins and the regions rich in pseudogenes or genetic vestiges may be contrasted.ConclusionThe SWGB applet presents a range of comprehensive OU statistical parameters calculated for a range of bacterial species, plasmids and phages. It is available on the Internet at http://www.bi.up.ac.za/SeqWord/mhhapplet.php.

Intraclonal diversity of the Pseudomonas aeruginosa cystic fibrosis airway isolates TBCF10839 and TBCF121838: distinct signatures of transcriptome, proteome, metabolome, adherence and pathogenicity despite an almost identical genome sequence

Microevolution of closely related Pseudomonas aeruginosa was compared in the clone TB strains TBCF10839 and TBCF121838 which had been isolated from two unrelated individuals with cystic fibrosis who had acquired clone TB during a local outbreak. Compared with the strain PAO1 reference sequence the two clone TB genomes shared 23 155 nucleotide exchanges, 32 out-of-frame indels in the coding region and another repertoire of replacement and genomic islands such as PAGI-1, PAGI-2, PAGI-5, LESGI-1 and LES-prophage 4. Only TBCF121838 carried a genomic island known from Ralstonia pickettii. Six of the seven strain-specific sequence variations in the core genome were detected in genes affecting motility, biofilm formation or virulence, i.e. non-synonymous nucleotide substitutions in mexS, PA3729, PA5017, mifR, a frameshift mutation in pilF (TBCF121838) and an intragenic deletion in pilQ (TBCF10839). Despite their almost identical genome sequence the two strains differed strongly from each other in transcriptome and metabolome profiles, mucin adherence and phagocytosis assays. TBCF121838 was susceptible to killing by neutrophils, but TBCF10839 could grow in leucocytes. Microevolution in P. aeruginosa apparently can generate novel complex traits by few or even single mutations provided that predisposing mutational events had occurred before in the clonal lineage.

Reversible Silencing of Cytomegalovirus Genomes by Type I Interferon Governs Virus Latency

Herpesviruses establish a lifelong latent infection posing the risk for virus reactivation and disease. In cytomegalovirus infection, expression of the major immediate early (IE) genes is a critical checkpoint, driving the lytic replication cycle upon primary infection or reactivation from latency. While it is known that type I interferon (IFN) limits lytic CMV replication, its role in latency and reactivation has not been explored. In the model of mouse CMV infection, we show here that IFNβ blocks mouse CMV replication at the level of IE transcription in IFN-responding endothelial cells and fibroblasts. The IFN-mediated inhibition of IE genes was entirely reversible, arguing that the IFN-effect may be consistent with viral latency. Importantly, the response to IFNβ is stochastic, and MCMV IE transcription and replication were repressed only in IFN-responsive cells, while the IFN-unresponsive cells remained permissive for lytic MCMV infection. IFN blocked the viral lytic replication cycle by upregulating the nuclear domain 10 (ND10) components, PML, Sp100 and Daxx, and their knockdown by shRNA rescued viral replication in the presence of IFNβ. Finally, IFNβ prevented MCMV reactivation from endothelial cells derived from latently infected mice, validating our results in a biologically relevant setting. Therefore, our data do not only define for the first time the molecular mechanism of IFN-mediated control of CMV infection, but also indicate that the reversible inhibition of the virus lytic cycle by IFNβ is consistent with the establishment of CMV latency.

Advances in computational analysis of metagenome sequences.

Second-generation sequencing technologies are revolutionizing the study of metagenomes. Whole-genome shotgun sequencing of metagenomic DNA may become an attractive alternative to the current widely used ribosomal RNA gene studies. Large data sets of short sequence reads are mapped onto a custom microbial reference sequence. If a bacterial pangenome of completely sequenced genomes is taken as a reference, the output consists of the distribution of bacterial taxa in and bacterial gene contents of the metagenome. The relative abundance of functional categories and of individual pathways and fitness traits encoded by the metagenomic gene pool provides insight into habitat-specific features of the microbial community. Polymorphic sites in sequence reads may resolve the number and abundance of individual clonal complexes of dominant species in the polymicrobial community. These SNPs and de novo mutations may be exploited to trace the spatiotemporal spread of clones and the emergence of novel traits such as fitness or resistance determinants. In conclusion, massively parallel sequencing of metagenomic DNA allows deep insights into the composition and the genetic repertoire of polymicrobial communities.

Genometa - A Fast and Accurate Classifier for Short Metagenomic Shotgun Reads

Summary Metagenomic studies use high-throughput sequence data to investigate microbial communities in situ. However, considerable challenges remain in the analysis of these data, particularly with regard to speed and reliable analysis of microbial species as opposed to higher level taxa such as phyla. We here present Genometa, a computationally undemanding graphical user interface program that enables identification of bacterial species and gene content from datasets generated by inexpensive high-throughput short read sequencing technologies. Our approach was first verified on two simulated metagenomic short read datasets, detecting 100% and 94% of the bacterial species included with few false positives or false negatives. Subsequent comparative benchmarking analysis against three popular metagenomic algorithms on an Illumina human gut dataset revealed Genometa to attribute the most reads to bacteria at species level (i.e. including all strains of that species) and demonstrate similar or better accuracy than the other programs. Lastly, speed was demonstrated to be many times that of BLAST due to the use of modern short read aligners. Our method is highly accurate if bacteria in the sample are represented by genomes in the reference sequence but cannot find species absent from the reference. This method is one of the most user-friendly and resource efficient approaches and is thus feasible for rapidly analysing millions of short reads on a personal computer. Availability The Genometa program, a step by step tutorial and Java source code are freely available from http://genomics1.mh-hannover.de/genometa/ and on http://code.google.com/p/genometa/. This program has been tested on Ubuntu Linux and Windows XP/7.

Intraclonal genome diversity of Pseudomonas aeruginosa clones CHA and TB

BackgroundAdaptation of Pseudomonas aeruginosa to different living conditions is accompanied by microevolution resulting in genomic diversity between strains of the same clonal lineage. In order to detect the impact of colonized habitats on P. aeruginosa microevolution we determined the genomic diversity between the highly virulent cystic fibrosis (CF) isolate CHA and two temporally and geographically unrelated clonal variants. The outcome was compared with the intraclonal genome diversity between three more closely related isolates of another clonal complex.ResultsThe three clone CHA isolates differed in their core genome in several dozen strain specific nucleotide exchanges and small deletions from each other. Loss of function mutations and non-conservative amino acid replacements affected several habitat- and lifestyle-associated traits, for example, the key regulator GacS of the switch between acute and chronic disease phenotypes was disrupted in strain CHA. Intraclonal genome diversity manifested in an individual composition of the respective accessory genome whereby the highest number of accessory DNA elements was observed for isolate PT22 from a polluted aquatic habitat. Little intraclonal diversity was observed between three spatiotemporally related outbreak isolates of clone TB. Although phenotypically different, only a few individual SNPs and deletions were detected in the clone TB isolates. Their accessory genome mainly differed in prophage-like DNA elements taken up by one of the strains.ConclusionsThe higher geographical and temporal distance of the clone CHA isolates was associated with an increased intraclonal genome diversity compared to the more closely related clone TB isolates derived from a common source demonstrating the impact of habitat adaptation on the microevolution of P. aeruginosa. However, even short-term habitat differentiation can cause major phenotypic diversification driven by single genomic variation events and uptake of phage DNA.