Strain-specific parallel evolution drives short-term diversification during Pseudomonas aeruginosa biofilm formation.

Generation of genetic diversity is a prerequisite for bacterial evolution and adaptation. Short-term diversification and selection within populations is, however, largely uncharacterised, as existing studies typically focus on fixed substitutions. Here, we use whole-genome deep-sequencing to capture the spectrum of mutations arising during biofilm development for two Pseudomonas aeruginosa strains. This approach identified single nucleotide variants with frequencies from 0.5% to 98.0% and showed that the clinical strain 18A exhibits greater genetic diversification than the type strain PA01, despite its lower per base mutation rate. Mutations were found to be strain specific: the mucoid strain 18A experienced mutations in alginate production genes and a c-di-GMP regulator gene; while PA01 acquired mutations in PilT and PilY1, possibly in response to a rapid expansion of a lytic Pf4 bacteriophage, which may use type IV pili for infection. The Pf4 population diversified with an evolutionary rate of 2.43 × 10(-3) substitutions per site per day, which is comparable to single-stranded RNA viruses. Extensive within-strain parallel evolution, often involving identical nucleotides, was also observed indicating that mutation supply is not limiting, which was contrasted by an almost complete lack of noncoding and synonymous mutations. Taken together, these results suggest that the majority of the P. aeruginosa genome is constrained by negative selection, with strong positive selection acting on an accessory subset of genes that facilitate adaptation to the biofilm lifecycle. Long-term bacterial evolution is known to proceed via few, nonsynonymous, positively selected mutations, and here we show that similar dynamics govern short-term, within-population bacterial diversification.

GemSIM: general, error-model based simulator of next-generation sequencing data.

BACKGROUND: GemSIM, or General Error-Model based SIMulator, is a next-generation sequencing simulator capable of generating single or paired-end reads for any sequencing technology compatible with the generic formats SAM and FASTQ (including Illumina and Roche/454). GemSIM creates and uses empirically derived, sequence-context based error models to realistically emulate individual sequencing runs and/or technologies. Empirical fragment length and quality score distributions are also used. Reads may be drawn from one or more genomes or haplotype sets, facilitating simulation of deep sequencing, metagenomic, and resequencing projects. RESULTS: We demonstrate GemSIM's value by deriving error models from two different Illumina sequencing runs and one Roche/454 run, and comparing and contrasting the resulting error profiles of each run. Overall error rates varied dramatically, both between individual Illumina runs, between the first and second reads in each pair, and between datasets from Illumina and Roche/454 technologies. Indels were markedly more frequent in Roche/454 than Illumina and both technologies suffered from an increase in error rates near the end of each read.The effects of these different profiles on low-frequency SNP-calling accuracy were investigated by analysing simulated sequencing data for a mixture of bacterial haplotypes. In general, SNP-calling using VarScan was only accurate for SNPs with frequency > 3%, independent of which error model was used to simulate the data. Variation between error profiles interacted strongly with VarScan's 'minumum average quality' parameter, resulting in different optimal settings for different sequencing runs. CONCLUSIONS: Next-generation sequencing has unprecedented potential for assessing genetic diversity, however analysis is complicated as error profiles can vary noticeably even between different runs of the same technology. Simulation with GemSIM can help overcome this problem, by providing insights into the error profiles of individual sequencing runs and allowing researchers to assess the effects of these errors on downstream data analysis.