Comparative genome sequencing of Escherichia coli allows observation of bacterial evolution on a laboratory timescale.

We applied whole-genome resequencing of Escherichia coli to monitor the acquisition and fixation of mutations that conveyed a selective growth advantage during adaptation to a glycerol-based growth medium. We identified 13 different de novo mutations in five different E. coli strains and monitored their fixation over a 44-d period of adaptation. We obtained proof that the observed spontaneous mutations were responsible for improved fitness by creating single, double and triple site-directed mutants that had growth rates matching those of the evolved strains. The success of this new genome-scale approach indicates that real-time evolution studies will now be practical in a wide variety of contexts.

Detection of hepatitis C virus transmission by use of DNA mass spectrometry.

The molecular detection of transmission of rapidly mutating pathogens such as hepatitis C virus (HCV) is commonly achieved by assessing the genetic relatedness of strains among infected patients. We describe the development of a novel mass spectrometry (MS)-based approach to identify HCV transmission. MS was used to detect products of base-specific cleavage of RNA molecules obtained from HCV polymerase chain reaction fragments. The MS-peak profiles were found to reflect variation in the HCV genomic sequence and the intrahost composition of the HCV population. Serum specimens originating from 60 case patients from 14 epidemiologically confirmed outbreaks and 25 unrelated controls were tested. Neighbor-joining trees constructed using MS-peak profile-based Hamming distances showed 100% accuracy, and linkage networks constructed using a threshold established from the Hamming distances between epidemiologically unrelated cases showed 100% sensitivity and 99.93% specificity in transmission detection. This MS-based approach is rapid, robust, reproducible, cost-effective, and applicable to investigating transmissions of other pathogens.

Molecular complexity of successive bacterial epidemics deconvoluted by comparative pathogenomics.

Understanding the fine-structure molecular architecture of bacterial epidemics has been a long-sought goal of infectious disease research. We used short-read-length DNA sequencing coupled with mass spectroscopy analysis of SNPs to study the molecular pathogenomics of three successive epidemics of invasive infections involving 344 serotype M3 group A Streptococcus in Ontario, Canada. Sequencing the genome of 95 strains from the three epidemics, coupled with analysis of 280 biallelic SNPs in all 344 strains, revealed an unexpectedly complex population structure composed of a dynamic mixture of distinct clonally related complexes. We discovered that each epidemic is dominated by micro- and macrobursts of multiple emergent clones, some with distinct strain genotype-patient phenotype relationships. On average, strains were differentiated from one another by only 49 SNPs and 11 insertion-deletion events (indels) in the core genome. Ten percent of SNPs are strain specific; that is, each strain has a unique genome sequence. We identified nonrandom temporal-spatial patterns of strain distribution within and between the epidemic peaks. The extensive full-genome data permitted us to identify genes with significantly increased rates of nonsynonymous (amino acid-altering) nucleotide polymorphisms, thereby providing clues about selective forces operative in the host. Comparative expression microarray analysis revealed that closely related strains differentiated by seemingly modest genetic changes can have significantly divergent transcriptomes. We conclude that enhanced understanding of bacterial epidemics requires a deep-sequencing, geographically centric, comparative pathogenomics strategy.

Replacing reverse line blot hybridization spoligotyping of the Mycobacterium tuberculosis complex.

Spoligotyping is a tool for the molecular characterization/typing of Mycobacterium tuberculosis complex (MTBC) strains based on target sequences (spacers) in the direct repeat (DR) region (14). The standard spoligotyping assay involves the hybridization of amplified sample DNA to nylon membrane-immobilized oligonucleotides whose sequences are representative of 43 spacer regions. Variations in the number of spacers as a result of deletions of adjacent blocks of repetitive units allow the differentiation of clinical isolates. In the present study, we developed a new multiplexed primer extension-based spoligotyping assay using automated matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) that improves the classical reverse line blot hybridization assay with respect to reproducibility, throughput, process flow, ease of use, and data analysis. Validation of the MALDI-TOF MS-based spoligotyping assay with two sample sets with a total of 326 samples resulted in 96.6% concordance (315/326) when the full spoligotype patterns were compared with the results of standard spoligotyping and 99.9% concordance when the results were compared with those of individual primer extension assays. Ten strains (including two Mycobacterium canettii strains) showed discordant results with one or two spacer differences from the membrane-based spoligotyping result. Most discordant samples were identified to be the result of ambiguities in the interpretation of weak hybridization signals in the reverse line blot assay and sequence variations in the spacer regions. We established a new automated primer extension assay and successfully validated it for use for the routine typing of MTBC strains in the research and public health laboratory environments. The present multiplex levels of up to 30 are extendable and allow the additional incorporation of controls and antibiotic resistance markers.

A simple method to elute cell-free HIV from dried blood spots improves their usefulness for monitoring therapy.

BACKGROUND: Dried blood spots (DBS) improve access to HIV viral load (VL) testing, but yield increased VL measurements compared to the plasma reference method because of cell-associated viral nucleic acid. In clinical settings, DBS methods may falsely categorize many patients as failing therapy. OBJECTIVES: Description of a simple method, free virus elution (FVE), to preferentially elute plasma-associated virus from DBS samples with phosphate-buffered saline, and an initial HIV VL performance comparison with standard DBS elution methods. STUDY DESIGN: The mechanism of action of FVE was studied with model DBS samples containing purified virus or washed HIV-containing cells, and with a DNA-specific HIV PCR. Using clinical samples, VL results from the new FVE method were compared to results from a dried fluid spot procedure (DFSP) protocol, which uses a guanidinium-based elution method, using plasma VL as the reference method. RESULTS: Model system experiments suggest that the method efficiently separates virus from cell-associated HIV, with a wide tolerance for incubation time and temperature. In 196 clinical samples, FVE reduced VL over-quantification from DBS, and improved DBS clinical concordance with plasma from 67% to 95%. CONCLUSIONS: A simple elution in PBS significantly reduced the over-quantification of HIV VL in DBS. Additional studies are needed to validate the method in fingerstick-collected specimens and to further understand the compartmentalization of HIV DNA and RNA in DBS specimens.

Combined genomarkers approach to Salmonella characterization reveals that nucleotide sequence differences in the phase 1 flagellin gene fliC are markers for variation within serotypes.

The characterization and tracking of pathogenic micro-organisms in the clinical laboratory and public health environment demand schemes that are easy to standardize and use, are automated and high-throughput, and provide portable data. A combined genomarkers approach for Salmonella enterica based on comparative sequence analysis by mass spectrometry has been developed. The scheme targets genes encoding synthesis and assembly of antigens, metabolic pathway enzymes, virulence factors and fluoroquinolone resistance, covering the essential sequences that distinguish between and identify variation within serotypes. This study demonstrated how this single method could replace the combination of methods currently required to determine serotypes, subtypes, antibiotic resistance profiles and the genomic relatedness of Salmonella isolates. The results revealed genomic variation within seven serotypes previously unreported. This variation can be detected by using nucleotide sequence differences in the Salmonella flagellin gene fliC as markers that are not detected by traditional serotyping methods.

Automated comparative sequence analysis by base-specific cleavage and mass spectrometry for nucleic acid-based microbial typing.

Traditional microbial typing technologies for the characterization of pathogenic microorganisms and monitoring of their global spread are often difficult to standardize and poorly portable, and they lack sufficient ease of use, throughput, and automation. To overcome these problems, we introduce the use of comparative sequencing by MALDI-TOF MS for automated high-throughput microbial DNA sequence analysis. Data derived from the public multilocus sequence typing (MLST) database (http://pubmlst.org/neisseria) established a reference set of expected peak patterns. A model pathogen, Neisseria meningitidis, was used to validate the technology and explore its applicability as an alternative to dideoxy sequencing. One hundred N. meningitidis samples were typed by comparing MALDI-TOF MS fingerprints of the standard MLST loci to reference sequences available in the public MLST database. Identification results can be obtained in 2 working days. Results were in concordance with classical dideoxy sequencing with 98% correct automatic identification. Sequence types (STs) of 89 samples were represented in the database, seven samples revealed new STs, including three new alleles, and four samples contained mixed populations of multiple STs. The approach shows interlaboratory reproducibility and allows for the exchange of mass spectrometric fingerprints to study the geographic spread of epidemic N. meningitidis strains or other microbes of clinical importance.

High-throughput mutation detection underlying adaptive evolution of Escherichia coli-K12.

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) analysis of base-specific cleavage products is an efficient, highly accurate tool for the detection of single base sequence variations. We describe the first application of this comparative sequencing strategy for automated high-throughput mutation detection in microbial genomes. The method was applied to identify DNA sequence changes that occurred in Escherichia coli K-12 MG1655 during laboratory adaptive evolution to new optimal growth phenotypes. Experiments were based on a genome-scale in silico model of E. coli metabolism and growth. This model computes several phenotypic functions and predicts optimal growth rates. To identify mutations underlying a 40-d adaptive laboratory evolution on glycerol, we resequenced 4.4% of the E. coli-K12 MG1655 genome in several clones picked at the end of the evolutionary process. The 1.54-Mb screen was completed in 13.5 h. This resequencing study is the largest reported by MALDI-TOF mass spectrometry to date. Ten mutations in 40 clones and three deviations from the reference sequence were detected. Mutations were predominantly found within the glycerol kinase gene. Functional characterization of the most prominent mutation shows its metabolic impact on the process of adaptive evolution. All sequence changes were independently confirmed by genotyping and Sanger-sequencing. We demonstrate that comparative sequencing by base-specific cleavage and MALDI-TOF mass spectrometry is an automated, fast, and highly accurate alternative to capillary sequencing.

Association testing by DNA pooling: an effective initial screen.

With an ever-increasing resource of validated single-nucleotide polymorphisms (SNPs), the limiting factors in genome-wide association analysis have become genotyping capacity and the availability of DNA. We provide a proof of concept of the use of pooled DNA as a means of efficiently screening SNPs and prioritizing them for further study. This approach reduces the final number of SNPs that undergo full, sample-by-sample genotyping as well as the quantity of DNA used overall. We have examined 15 SNPs in the cholesteryl ester transfer protein (CETP) gene, a gene previously demonstrated to be associated with serum high-density lipoprotein cholesterol levels. The SNPs were amplified in two pools of DNA derived from groups of individuals with extremely high and extremely low serum high-density lipoprotein cholesterol levels, respectively. P values <0.05 were obtained for 14 SNPs, supporting the described association. Genotyping of the individual samples showed that the average margin of error in frequency estimate was approximately 4% when pools were used. These findings clearly demonstrate the potential of pooling techniques and their associated technologies as an initial screen in the search for genetic associations.

Novel mass spectrometry-based tool for genotypic identification of mycobacteria.

Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) after base-specific cleavage of PCR amplified and in vitro-transcribed 16S rRNA gene (rDNA) was used for the identification of mycobacteria. Full-length 16S rDNA reference sequences of 12 type strains of Mycobacterium spp. frequently isolated from clinical specimens were determined by PCR, cloning, and sequencing. For MALDI-TOF MS-based comparative sequence analysis, mycobacterial 16S rDNA signature sequences ( approximately 500 bp) of the 12 type strains and 24 clinical isolates were PCR amplified using RNA promoter-tagged forward primers. T7 RNA polymerase-mediated transcription of forward strands in the presence of 5-methyl ribo-CTP maximized mass differences of fragments generated by base-specific cleavage. In vitro transcripts were subsequently treated with RNase T1, resulting in G-specific cleavage. Sample analysis by MALDI-TOF MS showed a specific mass signal pattern for each of the 12 type strains, allowing unambiguous identification. All 24 clinical isolates were identified unequivocally by comparing their detected mass signal pattern to the reference sequence-derived in silico pattern of the type strains and to the in silico mass patterns of published 16S rDNA sequences. A 16S rDNA microheterogeneity of the Mycobacterium xenopi type strain (DSM 43995) was detected by MALDI-TOF MS and later confirmed by Sanger dideoxy sequencing. In conclusion, analysis of 16S rDNA amplicons by MS after base-specific cleavage of RNA transcripts allowed fast and reliable identification of the Mycobacterium tuberculosis complex and ubiquitous mycobacteria (mycobacteria other than tuberculosis). The technology delivers an open platform for high-throughput microbial identification on the basis of any specific genotypic marker region.

High-throughput MALDI-TOF discovery of genomic sequence polymorphisms.

We describe a comparative sequencing strategy that is based on matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) analyses of complete base-specific cleavage reactions of a target sequence. The target is converted to a DNA/RNA mosaic structure after PCR amplification using in vitro transcription. Cleavage with defined specificity is achieved by ribonucleases. The set of cleavage products is subjected to mass spectrometry without prior fractionation. The presented resequencing assay is particularly useful for single-nucleotide polymorphism (SNP) discovery. The combination of mass spectra from four complementary cleavage reactions detects approximately 98% of all possible homozygous and heterozygous SNPs in target sequences with a length of up to 500 bases. In general, both the identity and location of the sequence variation are determined. This was exemplified by the discovery of SNPs in the human gene coding for the cholesteryl ester transfer protein using a panel of 96 genomic DNAs.