Joel S. Bader

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.

DNA Pooling: a tool for large-scale association studies

DNA pooling is a practical way to reduce the cost of large-scale association studies to identify susceptibility loci for common diseases. Pooling allows allele frequencies in groups of individuals to be measured using far fewer PCR reactions and genotyping assays than are used when genotyping individuals. Here, we discuss recent developments in quantitative genotyping assays and in the design and analysis of pooling studies. Sophisticated pooling designs are being developed that can take account of hidden population stratification, confounders and inter-loci interactions, and that allow the analysis of haplotypes.

Gaining confidence in high-throughput protein interaction networks

Although genome-scale technologies have benefited from statistical measures of data quality, extracting biologically relevant pathways from high-throughput proteomics data remains a challenge. Here we develop a quantitative method for evaluating proteomics data. We present a logistic regression approach that uses statistical and topological descriptors to predict the biological relevance of protein-protein interactions obtained from high-throughput screens for yeast. Other sources of information, including mRNA expression, genetic interactions and database annotations, are subsequently used to validate the model predictions without bias or cross-pollution. Novel topological statistics show hierarchical organization of the network of high-confidence interactions: protein complex interactions extend one to two links, and genetic interactions represent an even finer scale of organization. Knowledge of the maximum number of links that indicates a significant correlation between protein pairs (correlation distance) enables the integrated analysis of proteomics data with data from genetics and gene expression. The type of analysis presented will be essential for analyzing the growing amount of genomic and proteomics data in model organisms and humans.

Analysis of the human protein interactome and comparison with yeast, worm and fly interaction datasets.

We present the first analysis of the human proteome with regard to interactions between proteins. We also compare the human interactome with the available interaction datasets from yeast (Saccharomyces cerevisiae), worm (Caenorhabditis elegans) and fly (Drosophila melanogaster). Of >70,000 binary interactions, only 42 were common to human, worm and fly, and only 16 were common to all four datasets. An additional 36 interactions were common to fly and worm but were not observed in humans, although a coimmunoprecipitation assay showed that 9 of the interactions do occur in humans. A re-examination of the connectivity of essential genes in yeast and humans indicated that the available data do not support the presumption that the number of interaction partners can accurately predict whether a gene is essential. Finally, we found that proteins encoded by genes mutated in inherited genetic disorders are likely to interact with proteins known to cause similar disorders, suggesting the existence of disease subnetworks. The human interaction map constructed from our analysis should facilitate an integrative systems biology approach to elucidating the cellular networks that contribute to health and disease states.

Molecular model for aqueous ferrous–ferric electron transfer

We present a molecular model for studying the prototypical ferric–ferrous electron transfer process in liquid water, and we discuss its structural implications. Treatment of the nonequilibrium dynamics will be the subject of future work. The elementary constituents in the model are classical water molecules, classical ferric ions (i.e., Fe3+ particles), and a quantal electron. Pair potentials and pseudopotentials describing the interactions between these constituents are presented. These interactions lead to ligand structures of the ferric and ferrous ions that are in good agreement with those observed in nature. The validity of the tight binding model is examined. With umbrella sampling, we have computed the diabatic free energy of activation for electron transfer. The number obtained, roughly 20 kcal/mol, is in reasonable accord with the aqueous ferric–ferrous transfer activation energy of about 15 to 20 kcal/mol estimated from experiment. The Marcus relation for intersecting parabolic diabatic free ene...

Quantum and classical relaxation rates from classical simulations

The time correlation function for a harmonic quantum mechanical system can be related to the time correlation function for a corresponding classical system. Although straightforward to derive and well known in other contexts, this relationship has been unappreciated in the context of vibrational relaxation, where time correlation functions obtained from classical molecular dynamics have been used to predict relaxation rates for a quantum solute in a classical solvent. This inconsistent treatment—quantum solute, classical solvent—predicts a relaxation rate which is slower than if the entire system, both solute and solvent, were treated classically. We demonstrate that if the classical time correlation functions are rescaled to account for the ratio of quantum to classical fluctuations, providing a quantum mechanical treatment for the solute and the solvent, the relaxation rates and the entire absorption spectrum are the same as for a purely classical treatment. Our conclusions are valid when the solute and...

Total Synthesis of a Functional Designer Eukaryotic Chromosome

Designer Chromosome One of the ultimate aims of synthetic biology is to build designer organisms from the ground up. Rapid advances in DNA synthesis has allowed the assembly of complete bacterial genomes. Eukaryotic organisms, with their generally much larger and more complex genomes, present an additional challenge to synthetic biologists. Annaluru et al. (p. 55, published online 27 March) designed a synthetic eukaryotic chromosome based on yeast chromosome III. The designer chromosome, shorn of destabilizing transfer RNA genes and transposons, is ∼14% smaller than its wild-type template and is fully functional with every gene tagged for easy removal. A synthetic version of yeast chromosome III with every gene tagged can substitute for the original. Rapid advances in DNA synthesis techniques have made it possible to engineer viruses, biochemical pathways and assemble bacterial genomes. Here, we report the synthesis of a functional 272,871–base pair designer eukaryotic chromosome, synIII, which is based on the 316,617–base pair native Saccharomyces cerevisiae chromosome III. Changes to synIII include TAG/TAA stop-codon replacements, deletion of subtelomeric regions, introns, transfer RNAs, transposons, and silent mating loci as well as insertion of loxPsym sites to enable genome scrambling. SynIII is functional in S. cerevisiae. Scrambling of the chromosome in a heterozygous diploid reveals a large increase in a-mater derivatives resulting from loss of the MATα allele on synIII. The complete design and synthesis of synIII establishes S. cerevisiae as the basis for designer eukaryotic genome biology.

The minimum information required for reporting a molecular interaction experiment (MIMIx)

A wealth of molecular interaction data is available in the literature, ranging from large-scale datasets to a single interaction confirmed by several different techniques. These data are all too often reported either as free text or in tables of variable format, and are often missing key pieces of information essential for a full understanding of the experiment. Here we propose MIMIx, the minimum information required for reporting a molecular interaction experiment. Adherence to these reporting guidelines will result in publications of increased clarity and usefulness to the scientific community and will support the rapid, systematic capture of molecular interaction data in public databases, thereby improving access to valuable interaction data.

Synthetic chromosome arms function in yeast and generate phenotypic diversity by design

Recent advances in DNA synthesis technology have enabled the construction of novel genetic pathways and genomic elements, furthering our understanding of system-level phenomena. The ability to synthesize large segments of DNA allows the engineering of pathways and genomes according to arbitrary sets of design principles. Here we describe a synthetic yeast genome project, Sc2.0, and the first partially synthetic eukaryotic chromosomes, Saccharomyces cerevisiae chromosome synIXR, and semi-synVIL. We defined three design principles for a synthetic genome as follows: first, it should result in a (near) wild-type phenotype and fitness; second, it should lack destabilizing elements such as tRNA genes or transposons; and third, it should have genetic flexibility to facilitate future studies. The synthetic genome features several systemic modifications complying with the design principles, including an inducible evolution system, SCRaMbLE (synthetic chromosome rearrangement and modification by loxP-mediated evolution). We show the utility of SCRaMbLE as a novel method of combinatorial mutagenesis, capable of generating complex genotypes and a broad variety of phenotypes. When complete, the fully synthetic genome will allow massive restructuring of the yeast genome, and may open the door to a new type of combinatorial genetics based entirely on variations in gene content and copy number.

A systems biology analysis of the Drosophila phagosome.

Phagocytes have a critical function in remodelling tissues during embryogenesis and thereafter are central effectors of immune defence. During phagocytosis, particles are internalized into ‘phagosomes’, organelles from which immune processes such as microbial destruction and antigen presentation are initiated. Certain pathogens have evolved mechanisms to evade the immune system and persist undetected within phagocytes, and it is therefore evident that a detailed knowledge of this process is essential to an understanding of many aspects of innate and adaptive immunity. However, despite the crucial role of phagosomes in immunity, their components and organization are not fully defined. Here we present a systems biology analysis of phagosomes isolated from cells derived from the genetically tractable model organism Drosophila melanogaster and address the complex dynamic interactions between proteins within this organelle and their involvement in particle engulfment. Proteomic analysis identified 617 proteins potentially associated with Drosophila phagosomes; these were organized by protein–protein interactions to generate the ‘phagosome interactome’, a detailed protein–protein interaction network of this subcellular compartment. These networks predicted both the architecture of the phagosome and putative biomodules. The contribution of each protein and complex to bacterial internalization was tested by RNA-mediated interference and identified known components of the phagocytic machinery. In addition, the prediction and validation of regulators of phagocytosis such as the ‘exocyst’, a macromolecular complex required for exocytosis but not previously implicated in phagocytosis, validates this strategy. In generating this ‘systems-based model’, we show the power of applying this approach to the study of complex cellular processes and organelles and expect that this detailed model of the phagosome will provide a new framework for studying host–pathogen interactions and innate immunity.