Evolutionary dynamics of adult stem cells: Comparison of random and immortal strand segregation mechanisms

This paper develops a point-mutation model describing the evolutionary dynamics of a population of adult stem cells. Such a model may prove useful for quantitative studies of tissue aging and the emergence of cancer. We consider two modes of chromosome segregation: (1) Random segregation, where the daughter chromosomes of a given parent chromosome segregate randomly into the stem cell and its differentiating sister cell. (2) ``Immortal DNA strand'' co-segregation, for which the stem cell retains the daughter chromosomes with the oldest parent strands. Immortal strand co-segregation is a mechanism, originally proposed by Cairns (J. Cairns, {\it Nature} {\bf 255}, 197 (1975)), by which stem cells preserve the integrity of their genomes. For random segregation, we develop an ordered strand pair formulation of the dynamics, analogous to the ordered strand pair formalism developed for quasispecies dynamics involving semiconservative replication with imperfect lesion repair (in this context, lesion repair is taken to mean repair of postreplication base-pair mismatches). Interestingly, a similar formulation is possible with immortal strand co-segregation, despite the fact that this segregation mechanism is age-dependent. From our model we are able to mathematically show that, when lesion repair is imperfect, then immortal strand co-segregation leads to better preservation of the stem cell lineage than random chromosome segregation. Furthermore, our model allows us to estimate the optimal lesion repair efficiency for preserving an adult stem cell population for a given period of time. For human stem cells, we obtain that mispaired bases still present after replication and cell division should be left untouched, to avoid potentially fixing a mutation in both DNA strands.Comment: 9 pages, 3 figure

Irreversible thermodynamics of open chemical networks I: Emergent cycles and broken conservation laws

In this and a companion paper we outline a general framework for the thermodynamic description of open chemical reaction networks, with special regard to metabolic networks regulating cellular physiology and biochemical functions. We first introduce closed networks "in a box", whose thermodynamics is subjected to strict physical constraints: the mass-action law, elementarity of processes, and detailed balance. We further digress on the role of solvents and on the seemingly unacknowledged property of network independence of free energy landscapes. We then open the system by assuming that the concentrations of certain substrate species (the chemostats) are fixed, whether because promptly regulated by the environment via contact with reservoirs, or because nearly constant in a time window. As a result, the system is driven out of equilibrium. A rich algebraic and topological structure ensues in the network of internal species: Emergent irreversible cycles are associated to nonvanishing affinities, whose symmetries are dictated by the breakage of conservation laws. These central results are resumed in the relation $a + b = s^Y$ between the number of fundamental affinities $a$, that of broken conservation laws $b$ and the number of chemostats $s^Y$. We decompose the steady state entropy production rate in terms of fundamental fluxes and affinities in the spirit of Schnakenberg's theory of network thermodynamics, paving the way for the forthcoming treatment of the linear regime, of efficiency and tight coupling, of free energy transduction and of thermodynamic constraints for network reconstruction.Comment: 18 page

Design and Implementation of an Open Source Indexing Solution for a Large Set of Radiological Reports and Images

This paper hopes to share the insights we experienced during designing, building, and running an indexing solution for a large set of radiological reports and images in a production environment for more than 3 years. Several technical challenges were encountered and solved in the course of this project. One hundred four million words in 1.8 million radiological reports from 1989 to the present were indexed and became instantaneously searchable in a user-friendly fashion; the median query duration is only 31 ms. Currently, our highly tuned index holds 332,088 unique words in four languages. The indexing system is feature-rich and language-independent and allows for making complex queries. For research and training purposes it certainly is a valuable and convenient addition to our radiology informatics toolbox. Extended use of open-source technology dramatically reduced both implementation time and cost. All software we developed related to the indexing project has been made available to the open-source community covered by an unrestricted Berkeley Software Distribution-style license

A method for sensitivity analysis to assess the effects of measurement error in multiple exposure variables using external validation data

Measurement error in self-reported dietary intakes is known to bias the association between dietary intake and a health outcome of interest such as risk of a disease. The association can be distorted further by mismeasured confounders, leading to invalid results and conclusions. It is, however, difficult to adjust for the bias in the association when there is no internal validation data

Fluorescence Fluctuation Techniques for the Investigation of Structure-Function Relationships of G-Protein-Coupled Receptors

G-protein-coupled receptors (GPCRs) are seven transmembrane receptors that form the largest superfamily of signaling proteins, and the family members function in a diverse array of metabolic pathways including cardiac function, immune response, neurotransmission, smell, taste, cell differentiation and growth, and vision. It is becoming clear that alteration in the quaternary structure of the GPCR receptor can have a profound impact on signaling capabilities. Biochemical, biophysical, physiological, x-ray crystallographic, and computational methods have been used over the last 40–50 years to study the structure and function of GPCRs. Evidence from these studies confirm that GPCRs can be organized as monomers, dimers, and higher-order oligomers. However, many times, these methods require extraction of the receptor from its native environment and high levels of expression and only provide a snapshot of information. A need arose for techniques that could measure the assembly and disassembly of receptors at few-to-single molecule resolution in their native environment at fast time scales. In the last 20 years, fluorescence fluctuation techniques have filled this need and provided new insight into the dynamics of GPCR organization in the absence and presence of ligands, agonists, and antagonists. In this book chapter, we provide a brief introduction to GPCR structure and function [Section 1]. An overview of the theoretical basis for fluorescence fluctuations techniques (FFTs) and how FFTs can be used to study the oligomeric structure of GPCRs in live and fixed cells is explained [Section 2]. We discuss the advantages and limitations of FFTs [Section 3], and finally, we summarize select case studies on GPCR structure and function revealed by FFT experiments [Section 4]

Identification of Low Allele Frequency Mosaic Mutations in Alzheimer Disease

Germline mutations ofAPP,PSEN1, andPSEN2 genes cause autosomal dominant Alzheimer disease (AD). Somatic variants of the same genes may underlie pathogenesis in sporadic AD, which is the most prevalent form of the disease. Importantly, such somatic variants may be present at very low allelic frequency, confined to the brain, and are thus very difficult or impossible to detect in blood-derived DNA. Ever-refined methodologies to identify mutations present in a fraction of the DNA of the original tissue are rapidly transforming our understanding of DNA mutation and their role in complex pathologies such as tumors. These methods stand poised to test to what extend somatic variants may play a role in AD and other neurodegenerative diseases

A phenomenological approach to the simulation of metabolism and proliferation dynamics of large tumour cell populations

A major goal of modern computational biology is to simulate the collective behaviour of large cell populations starting from the intricate web of molecular interactions occurring at the microscopic level. In this paper we describe a simplified model of cell metabolism, growth and proliferation, suitable for inclusion in a multicell simulator, now under development (Chignola R and Milotti E 2004 Physica A 338 261-6). Nutrients regulate the proliferation dynamics of tumor cells which adapt their behaviour to respond to changes in the biochemical composition of the environment. This modeling of nutrient metabolism and cell cycle at a mesoscopic scale level leads to a continuous flow of information between the two disparate spatiotemporal scales of molecular and cellular dynamics that can be simulated with modern computers and tested experimentally.Comment: 58 pages, 7 figures, 3 tables, pdf onl

Principles guiding embryo selection following genome-wide haplotyping of preimplantation embryos.

STUDY QUESTION How to select and prioritize embryos during PGD following genome-wide haplotyping? SUMMARY ANSWER In addition to genetic disease-specific information, the embryo selected for transfer is based on ranking criteria including the existence of mitotic and/or meiotic aneuploidies, but not carriership of mutations causing recessive disorders. WHAT IS KNOWN ALREADY Embryo selection for monogenic diseases has been mainly performed using targeted disease-specific assays. Recently, these targeted approaches are being complemented by generic genome-wide genetic analysis methods such as karyomapping or haplarithmisis, which are based on genomic haplotype reconstruction of cell(s) biopsied from embryos. This provides not only information about the inheritance of Mendelian disease alleles but also about numerical and structural chromosome anomalies and haplotypes genome-wide. Reflections on how to use this information in the diagnostic laboratory are lacking. STUDY DESIGN, SIZE, DURATION We present the results of the first 101 PGD cycles (373 embryos) using haplarithmisis, performed in the Centre for Human Genetics, UZ Leuven. The questions raised were addressed by a multidisciplinary team of clinical geneticist, fertility specialists and ethicists. PARTICIPANTS/MATERIALS, SETTING, METHODS Sixty-three couples enrolled in the genome-wide haplotyping-based PGD program. Families presented with either inherited genetic variants causing known disorders and/or chromosomal rearrangements that could lead to unbalanced translocations in the offspring. MAIN RESULTS AND THE ROLE OF CHANCE Embryos were selected based on the absence or presence of the disease allele, a trisomy or other chromosomal abnormality leading to known developmental disorders. In addition, morphologically normal Day 5 embryos were prioritized for transfer based on the presence of other chromosomal imbalances and/or carrier information. LIMITATIONS, REASONS FOR CAUTION Some of the choices made and principles put forward are specific for cleavage-stage-based genetic testing. The proposed guidelines are subject to continuous update based on the accumulating knowledge from the implementation of genome-wide methods for PGD in many different centers world-wide as well as the results of ongoing scientific research. WIDER IMPLICATIONS OF THE FINDINGS Our embryo selection principles have a profound impact on the organization of PGD operations and on the information that is transferred among the genetic unit, the fertility clinic and the patients. These principles are also important for the organization of pre- and post-counseling and influence the interpretation and reporting of preimplantation genotyping results. As novel genome-wide approaches for embryo selection are revolutionizing the field of reproductive genetics, national and international discussions to set general guidelines are warranted. STUDY FUNDING/COMPETING INTEREST(S) The European Union's Research and Innovation funding programs FP7-PEOPLE-2012-IAPP SARM: 324509 and Horizon 2020 WIDENLIFE: 692065 to J.R.V., T.V., E.D. and M.Z.E. J.R.V., T.V. and M.Z.E. have patents ZL910050-PCT/EP2011/060211-WO/2011/157846 (‘Methods for haplotyping single cells’) with royalties paid and ZL913096-PCT/EP2014/068315-WO/2015/028576 (‘Haplotyping and copy-number typing using polymorphic variant allelic frequencies’) with royalties paid, licensed to Cartagenia (Agilent technologies). J.R.V. also has a patent ZL91 2076-PCT/EP20 one 3/070858 (‘High throughout genotyping by sequencing’) with royalties paid

Antigenic Conservation of an Immunodominant Invariable Region of the VlsE Lipoprotein among European Pathogenic Genospecies of Borrelia burgdorferi SL

Lyme disease is caused by genetically divergent spirochetes, including 3 pathogenic genospecies: Borrelia burgdorferi sensu stricto, B. garinii, and B. afzelii. Serodiagnosisis complicated by this genetic diversity. A synthetic peptide (C6), based on the 26-mer invariable region (IR6) of the variable surface antigen of B. burgdorferi (VlsE), was used as ELISA antigen, to test serum samples collected from mice experimentally infected with the 3 genospecies and from European patients with Lyme disease. Regardless of the infecting strains, mice produced a strong antibody response to C6, which indicates that IR6 is antigenically conserved among the pathogenic genospecies. Twenty of 23 patients with culture-confirmed erythema migrans had a detectable antibody response to C6. A sensitivity of 95.2% was achieved, with serum samples collected from patients with well-defined acrodermatitis chronica atrophicans. Fourteen of 20 patients with symptoms of late Lyme disease also had a positive anti-IR6 ELISA. Thus, it is possible that C6 may be used to serodiagnose Lyme disease universall