Rapid whole-genome mutational profiling using next-generation sequencing technologies.

Forward genetic mutational studies, adaptive evolution, and phenotypic screening are powerful tools for creating new variant organisms with desirable traits. However, mutations generated in the process cannot be easily identified with traditional genetic tools. We show that new high-throughput, massively parallel sequencing technologies can completely and accurately characterize a mutant genome relative to a previously sequenced parental (reference) strain. We studied a mutant strain of Pichia stipitis, a yeast capable of converting xylose to ethanol. This unusually efficient mutant strain was developed through repeated rounds of chemical mutagenesis, strain selection, transformation, and genetic manipulation over a period of seven years. We resequenced this strain on three different sequencing platforms. Surprisingly, we found fewer than a dozen mutations in open reading frames. All three sequencing technologies were able to identify each single nucleotide mutation given at least 10-15-fold nominal sequence coverage. Our results show that detecting mutations in evolved and engineered organisms is rapid and cost-effective at the whole-genome level using new sequencing technologies. Identification of specific mutations in strains with altered phenotypes will add insight into specific gene functions and guide further metabolic engineering efforts.

Whole-genome sequencing and variant discovery in C. elegans.

Massively parallel sequencing instruments enable rapid and inexpensive DNA sequence data production. Because these instruments are new, their data require characterization with respect to accuracy and utility. To address this, we sequenced a Caernohabditis elegans N2 Bristol strain isolate using the Solexa Sequence Analyzer, and compared the reads to the reference genome to characterize the data and to evaluate coverage and representation. Massively parallel sequencing facilitates strain-to-reference comparison for genome-wide sequence variant discovery. Owing to the short-read-length sequences produced, we developed a revised approach to determine the regions of the genome to which short reads could be uniquely mapped. We then aligned Solexa reads from C. elegans strain CB4858 to the reference, and screened for single-nucleotide polymorphisms (SNPs) and small indels. This study demonstrates the utility of massively parallel short read sequencing for whole genome resequencing and for accurate discovery of genome-wide polymorphisms.

Pyrobayes: an improved base caller for SNP discovery in pyrosequences.

Previously reported applications of the 454 Life Sciences pyrosequencing technology have relied on deep sequence coverage for accurate polymorphism discovery because of frequent insertion and deletion sequence errors. Here we report a new base calling program, Pyrobayes, for pyrosequencing reads. Pyrobayes permits accurate single-nucleotide polymorphism (SNP) calling in resequencing applications, even in shallow read coverage, primarily because it produces more confident base calls than the native base calling program.

Three dimensional transport lattice model for describing action potentials in axons stimulated by external electrodes.

Conditions that stimulate action potentials in one or more nerves is of widespread interest. Axon and nerve models are usually based on two dimensional pre-specified lumped equivalents that assume where currents will flow. In contrast, here we illustrate creation of three dimensional (3D) system models with a transport lattice of interconnected local models for external and internal electrolyte and axon membrane. The transport lattice solves Laplace's equation in the extracellular medium and is coupled to the Hodgkin-Huxley model at local membrane sites. These space-filling models incorporate the geometric scale, which allows explicit representation of confined axons and external electrodes. The present results demonstrate feasibility of the basic approach. These models are spatially coarse and approximate, but can be straightforwardly improved. The transport lattice system models are modular and multiscale (spatial scales ranging from the membrane thickness of 5 nm to the axon segment length of 2 cm).

Model of a confined spherical cell in uniform and heterogeneous applied electric fields.

Cells exposed to electric fields are often confined to a small volume within a solid tissue or within or near a device. Here we report on an approach to describing the frequency and time domain electrical responses of a spatially confined spherical cell by using a transport lattice system model. Two cases are considered: (1) a uniform applied field created by parallel plane electrodes, and (2) a heterogeneous applied field created by a planar electrode and a sharp microelectrode. Here fixed conductivities and dielectric permittivities of the extra- and intracellular media and of the membrane are used to create local transport models that are interconnected to create the system model. Consistent with traditional analytical solutions for spherical cells in an electrolyte of infinite extent, in the frequency domain the field amplification, G(m) (f) is large at low frequencies, f<1 MHz. G(m) (f) gradually decreases above 1 MHz and reaches a lower plateau at about 300 MHz, with the cell becoming almost "electrically invisible". In the time domain the application of a field pulse can result in altered localized transmembrane voltage changes due to a single microelectrode. The transport lattice approach provides modular, multiscale modeling capability that here ranges from cell membranes (5 nm scale) to the cell confinement volume ( approximately 40 microm scale).

Cylindrical cell membranes in uniform applied electric fields: validation of a transport lattice method.

The frequency and time domain transmembrane voltage responses of a cylindrical cell in an external electric field are calculated using a transport lattice, which allows solution of a variety of biologically relevant transport problems with complex cell geometry and field interactions. Here we demonstrate the method for a cylindrical membrane geometry and compare results with known analytical solutions. Results of transport lattice simulations on a Cartesian lattice are found to have discrepancies with the analytical solutions due to the limited volume of the system model and approximations for the local membrane model on the Cartesian lattice. Better agreement is attained when using a triangular mesh to represent the geometry rather than a Cartesian lattice. The transport lattice method can be readily extended to more sophisticated cell, organelle, and tissue configurations. Local membrane models within a system lattice can also include nonlinear responses such as electroporation and ion-channel gating.

Transport lattice models of heat transport in skin with spatially heterogeneous, temperature-dependent perfusion.

BACKGROUND: Investigation of bioheat transfer problems requires the evaluation of temporal and spatial distributions of temperature. This class of problems has been traditionally addressed using the Pennes bioheat equation. Transport of heat by conduction, and by temperature-dependent, spatially heterogeneous blood perfusion is modeled here using a transport lattice approach. METHODS: We represent heat transport processes by using a lattice that represents the Pennes bioheat equation in perfused tissues, and diffusion in nonperfused regions. The three layer skin model has a nonperfused viable epidermis, and deeper regions of dermis and subcutaneous tissue with perfusion that is constant or temperature-dependent. Two cases are considered: (1) surface contact heating and (2) spatially distributed heating. The model is relevant to the prediction of the transient and steady state temperature rise for different methods of power deposition within the skin. Accumulated thermal damage is estimated by using an Arrhenius type rate equation at locations where viable tissue temperature exceeds 42 degrees C. Prediction of spatial temperature distributions is also illustrated with a two-dimensional model of skin created from a histological image. RESULTS: The transport lattice approach was validated by comparison with an analytical solution for a slab with homogeneous thermal properties and spatially distributed uniform sink held at constant temperatures at the ends. For typical transcutaneous blood gas sensing conditions the estimated damage is small, even with prolonged skin contact to a 45 degrees C surface. Spatial heterogeneity in skin thermal properties leads to a non-uniform temperature distribution during a 10 GHz electromagnetic field exposure. A realistic two-dimensional model of the skin shows that tissue heterogeneity does not lead to a significant local temperature increase when heated by a hot wire tip. CONCLUSIONS: The heat transport system model of the skin was solved by exploiting the mathematical analogy between local thermal models and local electrical (charge transport) models, thereby allowing robust, circuit simulation software to obtain solutions to Kirchhoff's laws for the system model. Transport lattices allow systematic introduction of realistic geometry and spatially heterogeneous heat transport mechanisms. Local representations for both simple, passive functions and more complex local models can be easily and intuitively included into the system model of a tissue.