Genome sequencing in microfabricated high-density picolitre reactors.

The proliferation of large-scale DNA-sequencing projects in recent years has driven a search for alternative methods to reduce time and cost. Here we describe a scalable, highly parallel sequencing system with raw throughput significantly greater than that of state-of-the-art capillary electrophoresis instruments. The apparatus uses a novel fibre-optic slide of individual wells and is able to sequence 25 million bases, at 99% or better accuracy, in one four-hour run. To achieve an approximately 100-fold increase in throughput over current Sanger sequencing technology, we have developed an emulsion method for DNA amplification and an instrument for sequencing by synthesis using a pyrosequencing protocol optimized for solid support and picolitre-scale volumes. Here we show the utility, throughput, accuracy and robustness of this system by shotgun sequencing and de novo assembly of the Mycoplasma genitalium genome with 96% coverage at 99.96% accuracy in one run of the machine.

Assessing the feasibility of GS FLX Pyrosequencing for sequencing the Atlantic salmon genome.

BACKGROUND: With a whole genome duplication event and wealth of biological data, salmonids are excellent model organisms for studying evolutionary processes, fates of duplicated genes and genetic and physiological processes associated with complex behavioral phenotypes. It is surprising therefore, that no salmonid genome has been sequenced. Atlantic salmon (Salmo salar) is a good representative salmonid for sequencing given its importance in aquaculture and the genomic resources available. However, the size and complexity of the genome combined with the lack of a sequenced reference genome from a closely related fish makes assembly challenging. Given the cost and time limitations of Sanger sequencing as well as recent improvements to next generation sequencing technologies, we examined the feasibility of using the Genome Sequencer (GS) FLX pyrosequencing system to obtain the sequence of a salmonid genome. Eight pooled BACs belonging to a minimum tiling path covering approximately 1 Mb of the Atlantic salmon genome were sequenced by GS FLX shotgun and Long Paired End sequencing and compared with a ninth BAC sequenced by Sanger sequencing of a shotgun library. RESULTS: An initial assembly using only GS FLX shotgun sequences (average read length 248.5 bp) with approximately 30x coverage allowed gene identification, but was incomplete even when 126 Sanger-generated BAC-end sequences (approximately 0.09x coverage) were incorporated. The addition of paired end sequencing reads (additional approximately 26x coverage) produced a final assembly comprising 175 contigs assembled into four scaffolds with 171 gaps. Sanger sequencing of the ninth BAC (approximately 10.5x coverage) produced nine contigs and two scaffolds. The number of scaffolds produced by the GS FLX assembly was comparable to Sanger-generated sequencing; however, the number of gaps was much higher in the GS FLX assembly. CONCLUSION: These results represent the first use of GS FLX paired end reads for de novo sequence assembly. Our data demonstrated that this improved the GS FLX assemblies; however, with respect to de novo sequencing of complex genomes, the GS FLX technology is limited to gene mining and establishing a set of ordered sequence contigs. Currently, for a salmonid reference sequence, it appears that a substantial portion of sequencing should be done using Sanger technology.

Comprehensive resequence analysis of a 136 kb region of human chromosome 8q24 associated with prostate and colon cancers.

Recently, genome-wide association studies have identified loci across a segment of chromosome 8q24 (128,100,000-128,700,000) associated with the risk of breast, colon and prostate cancers. At least three regions of 8q24 have been independently associated with prostate cancer risk; the most centromeric of which appears to be population specific. Haplotypes in two contiguous but independent loci, marked by rs6983267 and rs1447295, have been identified in the Cancer Genetic Markers of Susceptibility project ( http://cgems.cancer.gov ), which genotyped more than 5,000 prostate cancer cases and 5,000 controls of European origin. The rs6983267 locus is also strongly associated with colorectal cancer. To ascertain a comprehensive catalog of common single-nucleotide polymorphisms (SNPs) across the two regions, we conducted a resequence analysis of 136 kb (chr8: 128,473,000-128,609,802) using the Roche/454 next-generation sequencing technology in 39 prostate cancer cases and 40 controls of European origin. We have characterized a comprehensive catalog of common (MAF > 1%) SNPs within this region, including 442 novel SNPs and have determined the pattern of linkage disequilibrium across the region. Our study has generated a detailed map of genetic variation across the region, which should be useful for choosing SNPs for fine mapping of association signals in 8q24 and investigations of the functional consequences of select common variants.

Whole-genome sequencing reveals a link between ß-lactam resistance and synthetases of the alarmone (p)ppGpp in Staphylococcus aureus.

The overwhelming majority of methicillin-resistant Staphylococcus aureus (MRSA) clinical isolates exhibit a peculiar heterogeneous resistance to ß-lactam antibiotics: in cultures of such strains, the majority of cells display only a low level of methicillin resistance--often close to the MIC breakpoint of susceptible strains. Yet, in the same cultures, subpopulations of bacteria exhibiting very high levels of resistance are also present with variable frequencies, which are characteristic of the particular MRSA lineage. The mechanism of heterogeneous resistance is not understood. We describe here an experimental system for exploring the mechanism of heterogeneous resistance. Copies of the resistance gene mecA cloned into a temperature-sensitive plasmid were introduced into the fully sequenced methicillin-susceptible clinical isolate S. aureus strain 476. Transductants of strain 476 expressed methicillin resistance in a heterogeneous fashion: the great majority of cells showed only low MIC (0.75 µg/ml) for the antibiotic, but a minority population of highly resistant bacteria (MIC >300 µg/ml) was also present with a frequency of ∼10(-4). The genetic backgrounds of the majority and minority cells were compared by whole-genome sequencing: the only differences detectable were two point mutations in relA of the highly resistant minority population of bacteria. The relA gene codes for the synthesis of (p)ppGpp, an effector of the stringent stress response. Titration of (p)ppGpp showed increased amounts of this effector in the highly resistant cells. Involvement of (p)ppGpp synthesis genes may explain some of the perplexing aspects of ß-lactam resistance in MRSA, since many environmental and genetic changes can modulate cellular levels of (p)ppGpp.