Integrating Sequencing Technologies in Personal Genomics: Optimal Low Cost Reconstruction of Structural Variants

The goal of human genome re-sequencing is obtaining an accurate assembly of an individual's genome. Recently, there has been great excitement in the development of many technologies for this (e.g. medium and short read sequencing from companies such as 454 and SOLiD, and high-density oligo-arrays from Affymetrix and NimbelGen), with even more expected to appear. The costs and sensitivities of these technologies differ considerably from each other. As an important goal of personal genomics is to reduce the cost of re-sequencing to an affordable point, it is worthwhile to consider optimally integrating technologies. Here, we build a simulation toolbox that will help us optimally combine different technologies for genome re-sequencing, especially in reconstructing large structural variants (SVs). SV reconstruction is considered the most challenging step in human genome re-sequencing. (It is sometimes even harder than de novo assembly of small genomes because of the duplications and repetitive sequences in the human genome.) To this end, we formulate canonical problems that are representative of issues in reconstruction and are of small enough scale to be computationally tractable and simulatable. Using semi-realistic simulations, we show how we can combine different technologies to optimally solve the assembly at low cost. With mapability maps, our simulations efficiently handle the inhomogeneous repeat-containing structure of the human genome and the computational complexity of practical assembly algorithms. They quantitatively show how combining different read lengths is more cost-effective than using one length, how an optimal mixed sequencing strategy for reconstructing large novel SVs usually also gives accurate detection of SNPs/indels, how paired-end reads can improve reconstruction efficiency, and how adding in arrays is more efficient than just sequencing for disentangling some complex SVs. Our strategy should facilitate the sequencing of human genomes at maximum accuracy and low cost

Effective Rheology of Bubbles Moving in a Capillary Tube

We calculate the average volumetric flux versus pressure drop of bubbles moving in a single capillary tube with varying diameter, finding a square-root relation from mapping the flow equations onto that of a driven overdamped pendulum. The calculation is based on a derivation of the equation of motion of a bubble train from considering the capillary forces and the entropy production associated with the viscous flow. We also calculate the configurational probability of the positions of the bubbles.Comment: 4 pages, 1 figur

Rise and Demise of Bioinformatics? Promise and Progress

The field of bioinformatics and computational biology has gone through a number of transformations during the past 15 years, establishing itself as a key component of new biology. This spectacular growth has been challenged by a number of disruptive changes in science and technology. Despite the apparent fatigue of the linguistic use of the term itself, bioinformatics has grown perhaps to a point beyond recognition. We explore both historical aspects and future trends and argue that as the field expands, key questions remain unanswered and acquire new meaning while at the same time the range of applications is widening to cover an ever increasing number of biological disciplines. These trends appear to be pointing to a redefinition of certain objectives, milestones, and possibly the field itself

A Large Gene Network in Immature Erythroid Cells Is Controlled by the Myeloid and B Cell Transcriptional Regulator PU.1

PU.1 is a hematopoietic transcription factor that is required for the development of myeloid and B cells. PU.1 is also expressed in erythroid progenitors, where it blocks erythroid differentiation by binding to and inhibiting the main erythroid promoting factor, GATA-1. However, other mechanisms by which PU.1 affects the fate of erythroid progenitors have not been thoroughly explored. Here, we used ChIP-Seq analysis for PU.1 and gene expression profiling in erythroid cells to show that PU.1 regulates an extensive network of genes that constitute major pathways for controlling growth and survival of immature erythroid cells. By analyzing fetal liver erythroid progenitors from mice with low PU.1 expression, we also show that the earliest erythroid committed cells are dramatically reduced in vivo. Furthermore, we find that PU.1 also regulates many of the same genes and pathways in other blood cells, leading us to propose that PU.1 is a multifaceted factor with overlapping, as well as distinct, functions in several hematopoietic lineages

Kinin B1 Receptor Enhances the Oxidative Stress in a Rat Model of Insulin Resistance: Outcome in Hypertension, Allodynia and Metabolic Complications

BACKGROUND: Kinin B(1) receptor (B(1)R) is induced by the oxidative stress in models of diabetes mellitus. This study aims at determining whether B(1)R activation could perpetuate the oxidative stress which leads to diabetic complications. METHODS AND FINDINGS: Young Sprague-Dawley rats were fed with 10% D-Glucose or tap water (controls) for 8-12 weeks. A selective B(1)R antagonist (SSR240612) was administered acutely (3-30 mg/kg) or daily for a period of 7 days (10 mg/kg) and the impact was measured on systolic blood pressure, allodynia, protein and/or mRNA B(1)R expression, aortic superoxide anion (O(2)(*-)) production and expression of superoxide dismutase (MnSOD) and catalase. SSR240612 reduced dose-dependently (3-30 mg/kg) high blood pressure in 12-week glucose-fed rats, but had no effect in controls. Eight-week glucose-fed rats exhibited insulin resistance (HOMA index), hypertension, tactile and cold allodynia and significant increases of plasma levels of glucose and insulin. This was associated with higher aortic levels of O(2)(*-), NADPH oxidase activity, MnSOD and catalase expression. All these abnormalities including B(1)R overexpression (spinal cord, aorta, liver and gastrocnemius muscle) were normalized by the prolonged treatment with SSR240612. The production of O(2)(*-) in the aorta of glucose-fed rats was also measured in the presence and absence of inhibitors (10-100 microM) of NADPH oxidase (apocynin), xanthine oxidase (allopurinol) or nitric oxide synthase (L-NAME) with and without Sar[D-Phe(8)]des-Arg(9)-BK (20 microM; B(1)R agonist). Data show that the greater aortic O(2)(*-) production induced by the B(1)R agonist was blocked only by apocynin. CONCLUSIONS: Activation of kinin B(1)R increased O(2)(*-) through the activation of NADPH oxidase in the vasculature. Prolonged blockade of B(1)R restored cardiovascular, sensory and metabolic abnormalities by reducing oxidative stress and B(1)R gene expression in this model