A generic, cost-effective, and scalable cell lineage analysis platform.

Advances in single-cell genomics enable commensurate improvements in methods for uncovering lineage relations among individual cells. Current sequencing-based methods for cell lineage analysis depend on low-resolution bulk analysis or rely on extensive single-cell sequencing, which is not scalable and could be biased by functional dependencies. Here we show an integrated biochemical-computational platform for generic single-cell lineage analysis that is retrospective, cost-effective, and scalable. It consists of a biochemical-computational pipeline that inputs individual cells, produces targeted single-cell sequencing data, and uses it to generate a lineage tree of the input cells. We validated the platform by applying it to cells sampled from an ex vivo grown tree and analyzed its feasibility landscape by computer simulations. We conclude that the platform may serve as a generic tool for lineage analysis and thus pave the way toward large-scale human cell lineage discovery.

Cell lineage analysis of the mammalian female germline.

Fundamental aspects of embryonic and post-natal development, including maintenance of the mammalian female germline, are largely unknown. Here we employ a retrospective, phylogenetic-based method for reconstructing cell lineage trees utilizing somatic mutations accumulated in microsatellites, to study female germline dynamics in mice. Reconstructed cell lineage trees can be used to estimate lineage relationships between different cell types, as well as cell depth (number of cell divisions since the zygote). We show that, in the reconstructed mouse cell lineage trees, oocytes form clusters that are separate from hematopoietic and mesenchymal stem cells, both in young and old mice, indicating that these populations belong to distinct lineages. Furthermore, while cumulus cells sampled from different ovarian follicles are distinctly clustered on the reconstructed trees, oocytes from the left and right ovaries are not, suggesting a mixing of their progenitor pools. We also observed an increase in oocyte depth with mouse age, which can be explained either by depth-guided selection of oocytes for ovulation or by post-natal renewal. Overall, our study sheds light on substantial novel aspects of female germline preservation and development.

Cell lineage analysis of acute leukemia relapse uncovers the role of replication-rate heterogeneity and microsatellite instability.

Human cancers display substantial intratumoral genetic heterogeneity, which facilitates tumor survival under changing microenvironmental conditions. Tumor substructure and its effect on disease progression and relapse are incompletely understood. In the present study, a high-throughput method that uses neutral somatic mutations accumulated in individual cells to reconstruct cell lineage trees was applied to hundreds of cells of human acute leukemia harvested from multiple patients at diagnosis and at relapse. The reconstructed cell lineage trees of patients with acute myeloid leukemia showed that leukemia cells at relapse were shallow (divide rarely) compared with cells at diagnosis and were closely related to their stem cell subpopulation, implying that in these instances relapse might have originated from rarely dividing stem cells. In contrast, among patients with acute lymphoid leukemia, no differences in cell depth were observed between diagnosis and relapse. In one case of chronic myeloid leukemia, at blast crisis, most of the cells at relapse were mismatch-repair deficient. In almost all leukemia cases, > 1 lineage was observed at relapse, indicating that diverse mechanisms can promote relapse in the same patient. In conclusion, diverse relapse mechanisms can be observed by systematic reconstruction of cell lineage trees of patients with leukemia.

Comparing algorithms that reconstruct cell lineage trees utilizing information on microsatellite mutations.

Organism cells proliferate and die to build, maintain, renew and repair it. The cellular history of an organism up to any point in time can be captured by a cell lineage tree in which vertices represent all organism cells, past and present, and directed edges represent progeny relations among them. The root represents the fertilized egg, and the leaves represent extant and dead cells. Somatic mutations accumulated during cell division endow each organism cell with a genomic signature that is unique with a very high probability. Distances between such genomic signatures can be used to reconstruct an organism's cell lineage tree. Cell populations possess unique features that are absent or rare in organism populations (e.g., the presence of stem cells and a small number of generations since the zygote) and do not undergo sexual reproduction, hence the reconstruction of cell lineage trees calls for careful examination and adaptation of the standard tools of population genetics. Our lab developed a method for reconstructing cell lineage trees by examining only mutations in highly variable microsatellite loci (MS, also called short tandem repeats, STR). In this study we use experimental data on somatic mutations in MS of individual cells in human and mice in order to validate and quantify the utility of known lineage tree reconstruction algorithms in this context. We employed extensive measurements of somatic mutations in individual cells which were isolated from healthy and diseased tissues of mice and humans. The validation was done by analyzing the ability to infer known and clear biological scenarios. In general, we found that if the biological scenario is simple, almost all algorithms tested can infer it. Another somewhat surprising conclusion is that the best algorithm among those tested is Neighbor Joining where the distance measure used is normalized absolute distance. We include our full dataset in Tables S1, S2, S3, S4, S5 to enable further analysis of this data by others.

Colon stem cell and crypt dynamics exposed by cell lineage reconstruction.

Stem cell dynamics in vivo are often being studied by lineage tracing methods. Our laboratory has previously developed a retrospective method for reconstructing cell lineage trees from somatic mutations accumulated in microsatellites. This method was applied here to explore different aspects of stem cell dynamics in the mouse colon without the use of stem cell markers. We first demonstrated the reliability of our method for the study of stem cells by confirming previously established facts, and then we addressed open questions. Our findings confirmed that colon crypts are monoclonal and that, throughout adulthood, the process of monoclonal conversion plays a major role in the maintenance of crypts. The absence of immortal strand mechanism in crypts stem cells was validated by the age-dependent accumulation of microsatellite mutations. In addition, we confirmed the positive correlation between physical and lineage proximity of crypts, by showing that the colon is separated into small domains that share a common ancestor. We gained new data demonstrating that colon epithelium is clustered separately from hematopoietic and other cell types, indicating that the colon is constituted of few progenitors and ruling out significant renewal of colonic epithelium from hematopoietic cells during adulthood. Overall, our study demonstrates the reliability of cell lineage reconstruction for the study of stem cell dynamics, and it further addresses open questions in colon stem cells. In addition, this method can be applied to study stem cell dynamics in other systems.

Detection of multiple disease indicators by an autonomous biomolecular computer.

The promise of biomolecular computers is their ability to interact with naturally occurring biomolecules, enabling in the future the development of context-dependent programmable drugs. Here we show a context-sensing mechanism of a biomolecular automaton that can simultaneously sense different types of molecules, allowing future integration of biomedical knowledge on a broad range of molecular disease symptoms in the decision of a biomolecular computer to release a drug molecule.

Comparison of seven single cell whole genome amplification commercial kits using targeted sequencing.

Advances in whole genome amplification (WGA) techniques enable understanding of the genomic sequence at a single cell level. Demand for single cell dedicated WGA kits (scWGA) has led to the development of several commercial kit. To this point, no robust comparison of all available kits was performed. Here, we benchmark an economical assay, comparing all commercially available scWGA kits. Our comparison is based on targeted sequencing of thousands of genomic loci, including highly mutable regions, from a large cohort of human single cells. Using this approach we have demonstrated the superiority of Ampli1 in genome coverage and of RepliG in reduced error rate. In summary, we show that no single kit is optimal across all categories, highlighting the need for a dedicated kit selection in accordance with experimental requirements.

Retrospective cell lineage reconstruction in humans by using short tandem repeats.

Cell lineage analysis aims to uncover the developmental history of an organism back to its cell of origin. Recently, novel in vivo methods utilizing genome editing enabled important insights into the cell lineages of animals. In contrast, human cell lineage remains restricted to retrospective approaches, which still lack resolution and cost-efficient solutions. Here, we demonstrate a scalable platform based on short tandem repeats targeted by duplex molecular inversion probes. With this human cell lineage tracing method, we accurately reproduced a known lineage of DU145 cells and reconstructed lineages of healthy and metastatic single cells from a melanoma patient who matched the anatomical reference while adding further refinements. This platform allowed us to faithfully recapitulate lineages of developmental tissue formation in healthy cells. In summary, our lineage discovery platform can profile informative somatic mutations efficiently and provides solid lineage reconstructions even in challenging low-mutation-rate healthy single cells.

An autonomous molecular computer for logical control of gene expression.

Early biomolecular computer research focused on laboratory-scale, human-operated computers for complex computational problems. Recently, simple molecular-scale autonomous programmable computers were demonstrated allowing both input and output information to be in molecular form. Such computers, using biological molecules as input data and biologically active molecules as outputs, could produce a system for 'logical' control of biological processes. Here we describe an autonomous biomolecular computer that, at least in vitro, logically analyses the levels of messenger RNA species, and in response produces a molecule capable of affecting levels of gene expression. The computer operates at a concentration of close to a trillion computers per microlitre and consists of three programmable modules: a computation module, that is, a stochastic molecular automaton; an input module, by which specific mRNA levels or point mutations regulate software molecule concentrations, and hence automaton transition probabilities; and an output module, capable of controlled release of a short single-stranded DNA molecule. This approach might be applied in vivo to biochemical sensing, genetic engineering and even medical diagnosis and treatment. As a proof of principle we programmed the computer to identify and analyse mRNA of disease-related genes associated with models of small-cell lung cancer and prostate cancer, and to produce a single-stranded DNA molecule modelled after an anticancer drug.

Recursive construction of perfect DNA molecules from imperfect oligonucleotides.

Making faultless complex objects from potentially faulty building blocks is a fundamental challenge in computer engineering, nanotechnology and synthetic biology. Here, we show for the first time how recursion can be used to address this challenge and demonstrate a recursive procedure that constructs error-free DNA molecules and their libraries from error-prone oligonucleotides. Divide and Conquer (D&C), the quintessential recursive problem-solving technique, is applied in silico to divide the target DNA sequence into overlapping oligonucleotides short enough to be synthesized directly, albeit with errors; error-prone oligonucleotides are recursively combined in vitro, forming error-prone DNA molecules; error-free fragments of these molecules are then identified, extracted and used as new, typically longer and more accurate, inputs to another iteration of the recursive construction procedure; the entire process repeats until an error-free target molecule is formed. Our recursive construction procedure surpasses existing methods for de novo DNA synthesis in speed, precision, amenability to automation, ease of combining synthetic and natural DNA fragments, and ability to construct designer DNA libraries. It thus provides a novel and robust foundation for the design and construction of synthetic biological molecules and organisms.

Estimating cell depth from somatic mutations.

The depth of a cell of a multicellular organism is the number of cell divisions it underwent since the zygote, and knowing this basic cell property would help address fundamental problems in several areas of biology. At present, the depths of the vast majority of human and mouse cell types are unknown. Here, we show a method for estimating the depth of a cell by analyzing somatic mutations in its microsatellites, and provide to our knowledge for the first time reliable depth estimates for several cells types in mice. According to our estimates, the average depth of oocytes is 29, consistent with previous estimates. The average depth of B cells ranges from 34 to 79, linearly related to the mouse age, suggesting a rate of one cell division per day. In contrast, various types of adult stem cells underwent on average fewer cell divisions, supporting the notion that adult stem cells are relatively quiescent. Our method for depth estimation opens a window for revealing tissue turnover rates in animals, including humans, which has important implications for our knowledge of the body under physiological and pathological conditions.

Differential activation of cysteine-substitution mutants of fibroblast growth factor receptor 3 is determined by cysteine localization.

Various human skeletal disorders are thought to be caused by mutations in fibroblast growth factor receptor 3 (FGFR3). These result in chronic FGFR3 hyperactivation and inhibition of bone growth. One such disorder, thanatophoric dysplasia, the most common form of sporadic, lethal dwarfism, is associated frequently with cysteine substitutions (G370C, S371C, and Y373C) in the extracellular juxtamembrane region of the receptor. These mutations have been suggested to induce disulfide-mediated receptor dimerization and constitutive activation. An adjacent cysteine substitution (G375C) leads to a less severe form of human dwarfism, achondroplasia, suggesting that the intensity of FGFR3 activation by these cross-links may be position dependent. To test this hypothesis, we have sequentially replaced each amino acid at positions 370-375 of FGFR3 with cysteine. Expression of each of these mutant forms in 293T cells led to their spontaneous, ligand-independent dimerization and increased basal phosphorylation. Wild-type (WT) FGFR3 became dimerized and phosphorylated only on FGF stimulation. Among the mutants, only two (G370C and S371C) caused high basal phosphorylation with significantly increased constitutive levels of mitogen-activated protein kinase (MAPK) phosphorylation and c-fos transcription. This activity was probably caused by mutant homodimer pairs, because WT-mutant heterodimers were observed only in the presence, but not in the absence, of FGF1. The high spontaneous activity of the mutants in positions 370-371, unlike those in 372-375, affirms their known involvement with thanatophoric dysplasia. We conclude that the G370C and S371C mutant receptors spontaneously dimerize in the correct spatial orientation required for effective signal transduction, whereas the 372-5 mutants, like the WT receptor, may achieve this orientation only on ligand binding.

A library of programmable DNAzymes that operate in a cellular environment.

DNAzymes were used as inhibitory agents in a variety of experimental disease settings, such as cancer, viral infections and even HIV. Drugs that become active only upon the presence of preprogrammed abnormal environmental conditions may enable selective molecular therapy by targeting abnormal cells without injuring normal cells. Here we show a novel programmable DNAzyme library composed of variety of Boolean logic gates, including YES, AND, NOT, OR, NAND, ANDNOT, XOR, NOR and 3-input-AND gate, that uses both miRNAs and mRNAs as inputs. Each gate is based on the c-jun cleaving Dz13 DNAzyme and active only in the presence of specific input combinations. The library is modular, supports arbitrary inputs and outputs, cascadable, highly specific and robust. We demonstrate the library's potential diagnostic abilities on miRNA and mRNA combinations in cell lysate and its ability to operate in a cellular environment by using beacon-like c-jun mimicking substrate in living mammalian cells.

Recursive construction and error correction of DNA molecules and libraries from synthetic and natural DNA.

Making error-free, custom DNA assemblies from potentially faulty building blocks is a fundamental challenge in synthetic biology. Here, we show how recursion can be used to address this challenge using a recursive procedure that constructs error-free DNA molecules and their libraries from error-prone synthetic oligonucleotides and naturally existing DNA. Specifically, we describe how divide and conquer (D&C), the quintessential recursive problem-solving technique, is applied in silico to divide target DNA sequences into overlapping, albeit error prone, oligonucleotides, and how recursive construction is applied in vitro to combine them to form error-prone DNA molecules. To correct DNA sequence errors, error-free fragments of these molecules are then identified, extracted, and used as new, typically longer and more accurate, inputs to another iteration of the recursive construction procedure; the entire process repeats until an error-free target molecule is formed. The method allows combining synthetic and natural DNA fragments into error-free designer DNA libraries, thus providing a foundation for the design and construction of complex synthetic DNA assemblies.

Muscle-bound primordial stem cells give rise to myofiber-associated myogenic and non-myogenic progenitors.

Myofiber cultures give rise to myogenic as well as to non-myogenic cells. Whether these myofiber-associated non-myogenic cells develop from resident stem cells that possess mesenchymal plasticity or from other stem cells such as mesenchymal stem cells (MSCs) remain unsolved. To address this question, we applied a method for reconstructing cell lineage trees from somatic mutations to MSCs and myogenic and non-myogenic cells from individual myofibers that were cultured at clonal density.Our analyses show that (i) in addition to myogenic progenitors, myofibers also harbor non-myogenic progenitors of a distinct, yet close, lineage; (ii) myofiber-associated non-myogenic and myogenic cells share the same muscle-bound primordial stem cells of a lineage distinct from bone marrow MSCs; (iii) these muscle-bound primordial stem-cells first part to individual muscles and then differentiate into myogenic and non-myogenic stem cells.

Reconstruction of cell lineage trees in mice.

The cell lineage tree of a multicellular organism represents its history of cell divisions from the very first cell, the zygote. A new method for high-resolution reconstruction of parts of such cell lineage trees was recently developed based on phylogenetic analysis of somatic mutations accumulated during normal development of an organism. In this study we apply this method in mice to reconstruct the lineage trees of distinct cell types. We address for the first time basic questions in developmental biology of higher organisms, namely what is the correlation between the lineage relation among cells and their (1) function, (2) physical proximity and (3) anatomical proximity. We analyzed B-cells, kidney-, mesenchymal- and hematopoietic-stem cells, as well as satellite cells, which are adult skeletal muscle stem cells isolated from their niche on the muscle fibers (myofibers) from various skeletal muscles. Our results demonstrate that all analyzed cell types are intermingled in the lineage tree, indicating that none of these cell types are single exclusive clones. We also show a significant correlation between the physical proximity of satellite cells within muscles and their lineage. Furthermore, we show that satellite cells obtained from a single myofiber are significantly clustered in the lineage tree, reflecting their common developmental origin. Lineage analysis based on somatic mutations enables performing high resolution reconstruction of lineage trees in mice and humans, which can provide fundamental insights to many aspects of their development and tissue maintenance.

Stochastic computing with biomolecular automata.

Stochastic computing has a broad range of applications, yet electronic computers realize its basic step, stochastic choice between alternative computation paths, in a cumbersome way. Biomolecular computers use a different computational paradigm and hence afford novel designs. We constructed a stochastic molecular automaton in which stochastic choice is realized by means of competition between alternative biochemical pathways, and choice probabilities are programmed by the relative molar concentrations of the software molecules coding for the alternatives. Programmable and autonomous stochastic molecular automata have been shown to perform direct analysis of disease-related molecular indicators in vitro and may have the potential to provide in situ medical diagnosis and cure.

DNA molecule provides a computing machine with both data and fuel.

The unique properties of DNA make it a fundamental building block in the fields of supramolecular chemistry, nanotechnology, nano-circuits, molecular switches, molecular devices, and molecular computing. In our recently introduced autonomous molecular automaton, DNA molecules serve as input, output, and software, and the hardware consists of DNA restriction and ligation enzymes using ATP as fuel. In addition to information, DNA stores energy, available on hybridization of complementary strands or hydrolysis of its phosphodiester backbone. Here we show that a single DNA molecule can provide both the input data and all of the necessary fuel for a molecular automaton. Each computational step of the automaton consists of a reversible software molecule input molecule hybridization followed by an irreversible software-directed cleavage of the input molecule, which drives the computation forward by increasing entropy and releasing heat. The cleavage uses a hitherto unknown capability of the restriction enzyme FokI, which serves as the hardware, to operate on a noncovalent software input hybrid. In the previous automaton, software input ligation consumed one software molecule and two ATP molecules per step. As ligation is not performed in this automaton, a fixed amount of software and hardware molecules can, in principle, process any input molecule of any length without external energy supply. Our experiments demonstrate 3 x 10(12) automata per microl performing 6.6 x 10(10) transitions per second per microl with transition fidelity of 99.9%, dissipating about 5 x 10(-9) W microl as heat at ambient temperature.

Human combinatorial Fab library yielding specific and functional antibodies against the human fibroblast growth factor receptor 3.

The human combinatorial antibody library Fab 1 (HuCAL-Fab 1) was generated by transferring the heavy and light chain variable regions from the previously constructed single-chain Fv library (Knappik, A., Ge, L., Honegger, A., Pack, P., Fischer, M., Wellnhofer, G., Hoess, A., Wölle, J., Plückthun, A., and Virnekäs, B. (2000) J. Mol. Biol. 296, 57-86), diversified in both complementarity-determining regions 3 into a novel Fab display vector, yielding 2.1 x 10(10) different antibody fragments. The modularity has been retained in the Fab display and screening plasmids, ensuring rapid conversion into various antibody formats as well as antibody optimization using prebuilt maturation cassettes. HuCAL-Fab 1 was challenged against the human fibroblast growth factor receptor 3, a potential therapeutic antibody target, against which, to the best of our knowledge, no functional antibodies could be generated so far. A unique screening mode was designed utilizing recombinant functional proteins and cell lines differentially expressing fibroblast growth factor receptor isoforms diversified in expression and receptor dependence. Specific Fab fragments with subnanomolar affinities were isolated by selection without any maturation steps as determined by fluorescence flow cytometry. Some of the selected Fab fragments completely inhibit target-mediated cell proliferation, rendering them the first monoclonal antibodies against fibroblast growth factor receptors having significant function blocking activity. This study validates HuCAL-Fab 1 as a valuable source for the generation of target-specific antibodies for therapeutic applications.