Sequence diversity of Pseudomonas aeruginosa: impact on population structure and genome evolution

Comparative sequencing of Pseudomonas aeruginosa genes oriC, citS, ampC, oprI, fliC, and pilA in 19 environmental and clinical isolates revealed the sequence diversity to be about 1 order of magnitude lower than in comparable housekeeping genes of Salmonella. In contrast to the low nucleotide substitution rate, the frequency of recombination among different P. aeruginosa genotypes was high, leading to the random association of alleles. The P. aeruginosa population consists of equivalent genotypes that form a net-like population structure. However, each genotype represents a cluster of closely related strains which retain their sequence signature in the conserved gene pool and carry a set of genotype-specific DNA blocks. The codon adaptation index, a quantitative measure of synonymous codon bias of genes, was found to be consistently high in the P. aeruginosa genome irrespective of the metabolic category and the abundance of the encoded gene product. Such uniformly high codon adaptation indices of 0.55 to 0.85 fit the ubiquitous lifestyle of P. aeruginosa.     

Genomic Variation among Contemporary Pseudomonas aeruginosa Isolates from Chronically Infected Cystic Fibrosis Patients

The airways of individuals with cystic fibrosis (CF) often become chronically infected with unique strains of the opportunistic pathogen Pseudomonas aeruginosa. Several lines of evidence suggest that the infecting P. aeruginosa lineage diversifies in the CF lung niche, yet so far this contemporary diversity has not been investigated at a genomic level. In this work, we sequenced the genomes of pairs of randomly selected contemporary isolates sampled from the expectorated sputum of three chronically infected adult CF patients. Each patient was infected by a distinct strain of P. aeruginosa. Single nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were identified in the DNA common to the paired isolates from different patients. The paired isolates from one patient differed due to just 1 SNP and 8 indels. The paired isolates from a second patient differed due to 54 SNPs and 38 indels. The pair of isolates from the third patient both contained a mutS mutation, which conferred a hypermutator phenotype; these isolates cumulatively differed due to 344 SNPs and 93 indels. In two of the pairs of isolates, a different accessory genome composition, specifically integrated prophage, was identified in one but not the other isolate of each pair. We conclude that contemporary isolates from a single sputum sample can differ at the SNP, indel, and accessory genome levels and that the cross-sectional genomic variation among coeval pairs of P. aeruginosa CF isolates can be comparable to the variation previously reported to differentiate between paired longitudinally sampled isolates.     

A major Pseudomonas aeruginosa clone common to patients and aquatic habitats.

The genomic relatedness of 573 Pseudomonas aeruginosa strains from environmental and clinical habitats was examined by digesting the genome with the rare-cutting enzyme SpeI. Thirty-nine strains were collected from environmental habitats mainly of aquatic origin, like rivers, lakes, or sanitary facilities. Four hundred fifty strains were collected from 76 patients with cystic fibrosis (CF) treated at four different centers, and 25 additional clinical isolates were collected from patients suffering from other diseases. Twenty-nine P. aeruginosa isolates were collected from the environment of one CF clinic. Thirty strains from culture collections were of environmental and clinic origin. A common macrorestriction fingerprint pattern was found in 13 of 46 CF patients, 5 of 29 environmental isolates from the same hospital, in a single ear infection isolate from another hospital, and 8 of 38 isolates from aquatic habitats about 300 km away from the CF clinic. The data indicate that closely related variants of one major clone (called clone C) persisted in various spatially and temporally separated habitats. Southern analysis of the clonal variants with six gene probes and two probes for genes coding for rRNA revealed almost the same hybridization patterns. With the exception of the phenotypically rapidly evolving CF isolates, the close relatedness of the strains of the clone was also shown by their identical responses in pyocin typing, phage typing, and serotyping. Besides clone C, three other P. aeruginosa clones were isolated from more than one clinical or environmental source.     

Genetic and phenotypic characterization of a Pseudomonas aeruginosa population with high frequency of genomic islands

Various genomic islands, PAPI-1, PAPI-2, PAGI-1, PAGI-2, PAGI-3, and PAGI-4, and the element pKLC102 have been characterized in different P. aeruginosa strains from diverse habitats and geographical locations. Chromosomal DNA macroarray of 100 P. aeruginosa strains isolated from 85 unrelated patients hospitalized in an intensive care unit was created to assess the occurrence of these genomic islands (GEIs). The macroarray was then hybridized with labeled probes derived from each genomic island. In addition, PFGE patterns with SpeI, frequency of virulence genes, and antimicrobial resistance patterns of the strains were studied. Our results showed that almost all P. aeruginosa strains presented up to eight virulence genes. By SpeI macrorestriction fragment analysis we were able to identify 49 restriction patterns; 35 patterns correspond to single strains and the remaining 14 to strains subgroup (a-n). Most of the strains showed variation in number or composition of GEIs and a specific antimicrobial pattern indicating that each strain was an unrelated isolate. In terms of the number of genomic islands per strain, 7 GEIs were found in 34% of the strains, 6 in 18%, 5 in 12%, 4 in 14%, 3 in 10%, 2 in 7%, and 1 in 4%; only one isolate did not present any GEI. The genomic islands PAPI-1 and PAPI-2 and the element pKLC102 were the most frequently detected. The analysis of the location of each GEI in the chromosome of two strains show that the islands PAGI-3, PAPI-1, PAPI-2 and pKLC102 are present in the insertion site previously reported, but that PAGI-2 and PAGI-4 are inserted in another chromosome place in a site not characterized yet. In conclusion our data show that P. aeruginosa strains exhibited an epidemic population structure with horizontal transfer of DNA resulting in a high frequency of GEIs.     

Comparative genome mapping of Pseudomonas aeruginosa PAO with P. aeruginosa C, which belongs to a major clone in cystic fibrosis patients and aquatic habitats.

A physical and genetic map was constructed for Pseudomonas aeruginosa C. Mainly, two-dimensional methods were used to place 47 SpeI, 8 PacI, 5 SwaI, and 4 I-CeuI sites onto the 6.5-Mb circular chromosome. A total of 21 genes, including the rrn operons and the origin of replication, were located on the physical map. Comparison of the physical and genetic map of strain C with that of the almost 600-kb-smaller genome of P. aeruginosa reference strain PAO revealed conservation of gene order between the two strains. A large-scale mosaic structure which was due to insertions of blocks of new genetic elements which had sizes of 23 to 155 kb and contained new SpeI sites was detected in the strain C chromosome. Most of these insertions were concentrated in three locations: two are congruent with the ends of the region rich in biosynthetic genes, and the third is located in the proposed region of the replication terminus. In addition, three insertions were scattered in the region rich in biosynthetic genes. The arrangement of the rrn operons around the origin of replication was conserved in C, PAO, and nine other examined independent strains.     

Association of virulence plasmid and antibiotic resistance determinants with chromosomal multilocus genotypes in Mexican Salmonella enterica serovar Typhimurium strains

Background  

Bacterial genomes are mosaic structures composed of genes present in every strain of the same species (core genome), and genes present in some but not all strains of a species (accessory genome). The aim of this study was to compare the genetic diversity of core and accessory genes of a Salmonella enterica subspecies enterica serovar Typhimurium (Typhimurium) population isolated from food-animal and human sources in four regions of Mexico. Multilocus sequence typing (MLST) and macrorestriction fingerprints by pulsed-field gel electrophoresis (PFGE) were used to address the core genetic variation, and genes involved in pathogenesis and antibiotic resistance were selected to evaluate the accessory genome.     

Single-nucleotide-polymorphism genotyping of Coxiella burnetii during a Q fever outbreak in The Netherlands

Coxiella burnetii is the etiological agent of Q fever. Currently, the Netherlands is facing the largest Q fever epidemic ever, with almost 4,000 notified human cases. Although the presence of a hypervirulent strain is hypothesized, epidemiological evidence, such as the animal reservoir(s) and genotype of the C. burnetii strain(s) involved, is still lacking. We developed a single-nucleotide-polymorphism (SNP) genotyping assay directly applicable to clinical samples. Ten discriminatory SNPs were carefully selected and detected by real-time PCR. SNP genotyping appeared to be highly suitable for discrimination of C. burnetii strains and easy to perform with clinical samples. With this new method, we show that the Dutch outbreak is caused by at least 5 different C. burnetii genotypes. SNP typing of 14 human samples from the outbreak revealed the presence of 3 dissimilar genotypes. Two genotypes were also present in livestock at 9 farms in the outbreak area. SNP analyses of bulk milk from 5 other farms, commercial cow milk, and cow colostrum revealed 2 additional genotypes that were not detected in humans. SNP genotyping data from clinical samples clearly demonstrate that at least 5 different C. burnetii genotypes are involved in the Dutch outbreak.     

欧洲黑杨基因资源材性关联基因的SNP分析

以115个欧洲黑杨(Populus nigra L.)无性系为材料, 利用TaqMan技术分析了欧洲黑杨基因资源参与木质素和纤维素合成的酶(4CL、PAL和CesA2)的单核苷酸多态性, 并对分型的SNPs与木材材性性状(物理性状：基本密度、纤维长、纤维宽、微纤丝角; 化学性状: 木质素含量、纤维素含量、a 纤维素含量等)进行了相关分析。结果如下: (1)在对4CL、PAL和CesA2等3个基因进行检测时, 共获得27个SNPs标记, 对其中转换(A-G, C-T)有17个位点, 颠换(A-C, G-C, G-T, A-T等)有10个位点; (2)对其中的3个SNPs进行了分型, 分别记作SNP1、SNP2和SNP3; (2)对已经分型SNPs位点与材性性状进行方差分析, 结果显示, 3个SNPs中只有SNP1与4年生欧洲黑杨综纤维素含量显著相关, 表现为负效应, 贡献率为11.11%; (3)对欧洲黑杨4CL基因的SNP1标记的不同基因型所对应的材性性状进行方差分析, 结果显示基因型为CC和CT的欧洲黑杨相对于基因型为TT的欧洲黑杨有较高的纤维素含量。     

Imputation of single-nucleotide polymorphisms in inbred mice using local phylogeny

We present full-genome genotype imputations for 100 classical laboratory mouse strains, using a novel method. Using genotypes at 549,683 SNP loci obtained with the Mouse Diversity Array, we partitioned the genome of 100 mouse strains into 40,647 intervals that exhibit no evidence of historical recombination. For each of these intervals we inferred a local phylogenetic tree. We combined these data with 12 million loci with sequence variations recently discovered by whole-genome sequencing in a common subset of 12 classical laboratory strains. For each phylogenetic tree we identified strains sharing a leaf node with one or more of the sequenced strains. We then imputed high- and medium-confidence genotypes for each of 88 nonsequenced genomes. Among inbred strains, we imputed 92% of SNPs genome-wide, with 71% in high-confidence regions. Our method produced 977 million new genotypes with an estimated per-SNP error rate of 0.083% in high-confidence regions and 0.37% genome-wide. Our analysis identified which of the 88 nonsequenced strains would be the most informative for improving full-genome imputation, as well as which additional strain sequences will reveal more new genetic variants. Imputed sequences and quality scores can be downloaded and visualized online.     

Tracing of Salmonella spp. in two pork slaughter and cutting plants using serotyping and macrorestriction genotyping

AIMS: The origin of Salmonella contamination of pork products is not well established. In order to further this knowledge, the transmission of Salmonella spp. from live pigs to pork cuts was investigated in two pork slaughter and cutting plants. METHODS AND RESULTS: Salmonella spp. were isolated from both pork (pigs, carcasses, cuts) and the environment before and during slaughterhouse activities. Eight serotypes were identified. XbaI and SpeI macrorestriction distinguished 20 genotypes of Salmonella Typhimurium and 16 genotypes of Salmonella Derby. A major cluster of Salmonella Typhimurium genotypes was common to both plants and all pig-related genotypes, while a predominant pig-related Salmonella Derby genotype was common to both plants. CONCLUSION: None of the Salmonella strains persisted for long periods in the pork-processing environments. SIGNIFICANCE AND IMPACT OF THE STUDY: This work shows that contaminated live pigs, because of bacterial spread due to the process and ineffective cleaning procedures, are involved in Salmonella contamination.     

The physical map of the chromosome of a serogroup A strain of Neisseria meningitidis shows complex rearrangements relative to the chromosomes of the two mapped strains of the closely related species N. gonorrhoeae.

A physical map of the chromosome of N. meningitidis Z2491 (serogroup A, subgroup IV-1) has been constructed. Z2491 DNA was digested with NheI, SpeI, SgfI, PacI, BglII, or PmeI, resulting in a limited number of fragments that were resolved by contour-clamped homogeneous electric field (CHEF) electrophoresis. The estimated genome size for this strain was 2,226 kb. To construct the map, probes corresponding to single-copy genes or sequences were used on Southern blots of chromosomal DNA digested with the different mapping enzymes and subjected to CHEF electrophoresis. By determining which fragments from different digests hybridized to each specific probe, it was possible to walk back and forth between digests to form a circular macrorestriction map. The intervals between mapped restriction sites range from 10 to 143 kb in size. A total of 117 markers have been placed on the map; 75 represent identified genes, with the remaining markers defined by anonymous cloned fragments of neisserial DNA. Comparison of the arrangement of genetic loci in Z2491 with that in gonococcal strain FA1090, for which a physical map was previously constructed, revealed complex genomic rearrangements between the two strains. Although gene order is generally conserved over much of the chromosome, a region of approximately 500 kb shows translocation and/or inversion of multiple blocks of markers between the two strains. Even within the relatively conserved portions of the maps, several genetic markers are in different positions in Z2491 and FA1090.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

SNP/RD typing of Mycobacterium tuberculosis Beijing strains reveals local and worldwide disseminated clonal complexes

The Beijing strain is one of the most successful genotypes of Mycobacterium tuberculosis worldwide and appears to be highly homogenous according to existing genotyping methods. To type Beijing strains reliably we developed a robust typing scheme using single nucleotide polymorphisms (SNPs) and regions of difference (RDs) derived from whole-genome sequencing data of eight Beijing strains. SNP/RD typing of 259 M. tuberculosis isolates originating from 45 countries worldwide discriminated 27 clonal complexes within the Beijing genotype family. A total of 16 Beijing clonal complexes contained more than one isolate of known origin, of which two clonal complexes were strongly associated with South African origin. The remaining 14 clonal complexes encompassed isolates from different countries. Even highly resolved clonal complexes comprised isolates from distinct geographical sites. Our results suggest that Beijing strains spread globally on multiple occasions and that the tuberculosis epidemic caused by the Beijing genotype is at least partially driven by modern migration patterns. The SNPs and RDs presented in this study will facilitate future molecular epidemiological and phylogenetic studies on Beijing strains.     

The dominant Australian community-acquired methicillin-resistant Staphylococcus aureus clone ST93-IV [2B] is highly virulent and genetically distinct

Community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) USA300 has spread rapidly across North America, and CA-MRSA is also increasing in Australia. However, the dominant Australian CA-MRSA strain, ST93-IV [2B] appears distantly related to USA300 despite strikingly similar clinical and epidemiological profiles. Here, we compared the virulence of a recent Australian ST93 isolate (JKD6159) to other MRSA, including USA300, and found that JKD6159 was the most virulent in a mouse skin infection model. We fully sequenced the genome of JKD6159 and confirmed that JKD6159 is a distinct clone with 7616 single nucleotide polymorphisms (SNPs) distinguishing this strain from all other S. aureus genomes. Despite its high virulence there were surprisingly few virulence determinants. However, genes encoding α-hemolysin, Panton-Valentine leukocidin (PVL) and α-type phenol soluble modulins were present. Genome comparisons revealed 32 additional CDS in JKD6159 but none appeared to encode new virulence factors, suggesting that this clone's enhanced pathogenicity could lie within subtler genome changes, such as SNPs within regulatory genes. To investigate the role of accessory genome elements in CA-MRSA epidemiology, we next sequenced three additional Australian non-ST93 CA-MRSA strains and compared them with JKD6159, 19 completed S. aureus genomes and 59 additional S. aureus genomes for which unassembled genome sequence data was publicly available (82 genomes in total). These comparisons showed that despite its distinctive genotype, JKD6159 and other CA-MRSA clones (including USA300) share a conserved repertoire of three notable accessory elements (SSCmecIV, PVL prophage, and pMW2). This study demonstrates that the genetically distinct ST93 CA-MRSA from Australia is highly virulent. Our comparisons of geographically and genetically diverse CA-MRSA genomes suggest that apparent convergent evolution in CA-MRSA may be better explained by the rapid dissemination of a highly conserved accessory genome from a common source.     

2SNP: scalable phasing based on 2-SNP haplotypes

2SNP software package implements a new very fast scalable algorithm for haplotype inference based on genotype statistics collected only for pairs of SNPs. This software can be used for comparatively accurate phasing of large number of long genome sequences, e.g. obtained from DNA arrays. As an input 2SNP takes genotype matrix and outputs the corresponding haplotype matrix. On datasets across 79 regions from HapMap 2SNP is several orders of magnitude faster than GERBIL and PHASE while matching them in quality measured by the number of correctly phased genotypes, single-site and switching errors. For example, 2SNP requires 41 s on Pentium 4 2 Ghz processor to phase 30 genotypes with 1381 SNPs (ENm010.7p15:2 data from HapMap) versus GERBIL and PHASE requiring more than a week and admitting no less errors than 2SNP.     

The development of a genome wide SNP set for the Barnacle goose Branta leucopsis

Migratory birds are of particular interest for population genetics because of the high connectivity between habitats and populations. A high degree of connectivity requires using many genetic markers to achieve the required statistical power, and a genome wide SNP set can fit this purpose. Here we present the development of a genome wide SNP set for the Barnacle Goose Branta leucopsis, a model species for the study of bird migration. We used the genome of a different waterfowl species, Mallard Anas platyrhynchos, as a reference to align Barnacle Goose second generation sequence reads from an RRL library and detected 2188 SNPs genome wide. Furthermore, we used chimeric flanking sequences, merged from both Mallard and Barnacle Goose DNA sequence information, to create primers for validation by genotyping. Validation with a 384 SNP genotyping set resulted in 374 (97%) successfully typed SNPs in the assay, of which 358 (96%) were polymorphic. Additionally, we validated our SNPs on relatively old (30 years) museum samples, which resulted in a success rate of at least 80%. This shows that museum samples could be used in standard SNP genotyping assays. Our study also shows that the genome of a related species can be used as reference to detect genome wide SNPs in birds, because genomes of birds are highly conserved. This is illustrated by the use of chimeric flanking sequences, which showed that the incorporation of flanking nucleotides from Mallard into Barnacle Goose sequences lead to equal genotyping performance when compared to flanking sequences solely composed of Barnacle Goose sequence.     

A combined physical and genetic map of Pseudomonas aeruginosa PAO

A combined physical and genetic map of Pseudomonas aeruginosa PAO was constructed by pulsed-field gel electrophoresis and Southern hybridization using cosmid clones from a genomic library carrying known genes. A total of 37 SpeI restriction fragments have been mapped on the 5862 kb genome, and fragment contiguity demonstrated by hybridization with clones from a SpeI junction fragment library and fragments obtained by partial SpeI digestion, both derived from the P. aeruginosa PAO chromosome.     

Flagellin gene sequence variation in the genus Pseudomonas

Flagellin gene (fliC) sequences from 18 strains of Pseudomonas sensu stricto representing 8 different species, and 9 representative fliC sequences from other members of the gamma sub-division of proteobacteria, were compared. Analysis was performed on N-terminal, C-terminal and whole fliC sequences. The fliC analyses confirmed the inferred relationship between P. mendocina, P. oleovorans and P. aeruginosa based on 16S rRNA sequence comparisons. In addition, the analyses indicated that P. putida PRS2000 was closely related to P. fluorescens SBW25 and P. fluorescens NCIMB 9046T, but suggested that P. putida PaW8 and P. putida PRS2000 were more closely related to other Pseudomonas spp. than they were to each other. There were a number of inconsistencies in inferred evolutionary relationships between strains, depending on the analysis performed. In particular, whole flagellin gene comparisons often differed from those obtained using N- and C-terminal sequences. However, there were also inconsistencies between the terminal region analyses, suggesting that phylogenetic relationships inferred on the basis of fliC sequence should be treated with caution. Although the central domain of fliC is highly variable between Pseudomonas strains, there was evidence of sequence similarities between the central domains of different Pseudomonas fliC sequences. This indicates the possibility of recombination in the central domain of fliC genes within Pseudomonas species, and between these genes and those from other bacteria.