DNA assembly for nanopore data storage readout.

Synthetic DNA is becoming an attractive substrate for digital data storage due to its density, durability, and relevance in biological research. A major challenge in making DNA data storage a reality is that reading DNA back into data using sequencing by synthesis remains a laborious, slow and expensive process. Here, we demonstrate successful decoding of 1.67 megabytes of information stored in short fragments of synthetic DNA using a portable nanopore sequencing platform. We design and validate an assembly strategy for DNA storage that drastically increases the throughput of nanopore sequencing. Importantly, this assembly strategy is generalizable to any application that requires nanopore sequencing of small DNA amplicons.

Random access in large-scale DNA data storage.

Synthetic DNA is durable and can encode digital data with high density, making it an attractive medium for data storage. However, recovering stored data on a large-scale currently requires all the DNA in a pool to be sequenced, even if only a subset of the information needs to be extracted. Here, we encode and store 35 distinct files (over 200 MB of data), in more than 13 million DNA oligonucleotides, and show that we can recover each file individually and with no errors, using a random access approach. We design and validate a large library of primers that enable individual recovery of all files stored within the DNA. We also develop an algorithm that greatly reduces the sequencing read coverage required for error-free decoding by maximizing information from all sequence reads. These advances demonstrate a viable, large-scale system for DNA data storage and retrieval.

Can one hear the shape of a population history?

Reconstructing past population size from present day genetic data is a major goal of population genetics. Recent empirical studies infer population size history using coalescent-based models applied to a small number of individuals. Here we provide tight bounds on the amount of exact coalescence time data needed to recover the population size history of a single, panmictic population at a certain level of accuracy. In practice, coalescence times are estimated from sequence data and so our lower bounds should be taken as rather conservative.

Analysis of unstable behavior in a mathematical model for erythropoiesis.

We consider an age-structured model that describes the regulation of erythropoiesis through the negative feedback loop between erythropoietin and hemoglobin. This model is reduced to a system of two ordinary differential equations with two constant delays for which we show existence of a unique steady state. We determine all instances at which this steady state loses stability via a Hopf bifurcation through a theoretical bifurcation analysis establishing analytical expressions for the scenarios in which they arise. We show examples of supercritical Hopf bifurcations for parameter values estimated according to physiological values for humans found in the literature and present numerical simulations in agreement with the theoretical analysis. We provide a strategy for parameter estimation to match empirical measurements and predict dynamics in experimental settings, and compare existing data on hemoglobin oscillation in rabbits with predictions of our model.