Quantitative theory of telomere length regulation and cellular senescence. Proc Natl Acad Sci U S A. 2010;107:5387-5392. [PMID: 20207949 DOI: 10.1073/pnas.0914502107]

The multifaceted functions of DNA-PKcs: implications for the therapy of human diseases

The DNA-dependent protein kinase (DNA-PK), catalytic subunit, also known as DNA-PKcs, is complexed with the heterodimer Ku70/Ku80 to form DNA-PK holoenzyme, which is well recognized as initiator in the nonhomologous end joining (NHEJ) repair after double strand break (DSB). During NHEJ, DNA-PKcs is essential for both DNA end processing and end joining. Besides its classical function in DSB repair, DNA-PKcs also shows multifaceted functions in various biological activities such as class switch recombination (CSR) and variable (V) diversity (D) joining (J) recombination in B/T lymphocytes development, innate immunity through cGAS-STING pathway, transcription, alternative splicing, and so on, which are dependent on its function in NHEJ or not. Moreover, DNA-PKcs deficiency has been proven to be related with human diseases such as neurological pathogenesis, cancer, immunological disorder, and so on through different mechanisms. Therefore, it is imperative to summarize the latest findings about DNA-PKcs and diseases for better targeting DNA-PKcs, which have shown efficacy in cancer treatment in preclinical models. Here, we discuss the multifaceted roles of DNA-PKcs in human diseases, meanwhile, we discuss the progresses of DNA-PKcs inhibitors and their potential in clinical trials. The most updated review about DNA-PKcs will hopefully provide insights and ideas to understand DNA-PKcs associated diseases.

Review: Are moles senescent?

Melanocytic nevi (skin moles) have been regarded as a valuable example of cell senescence occurring in vivo. However, a study of induced nevi in a mouse model reported that the nevi were arrested by cell interactions rather than a cell-autonomous process like senescence, and that size distributions of cell nests within nevi could not be accounted for by a stochastic model of oncogene-induced senescence. Moreover, others reported that some molecular markers used to identify cell senescence in human nevi are also found in melanoma cells-not senescent. It has thus been questioned whether nevi really are senescent, with potential implications for melanoma diagnosis and therapy. Here I review these areas, along with the genetic, biological, and molecular evidence supporting senescence in nevi. In conclusion, there is strong evidence that cells of acquired human benign (banal) nevi are very largely senescent, though some must contain a minor non-senescent cell subpopulation. There is also persuasive evidence that this senescence is primarily induced by dysfunctional telomeres rather than directly oncogene-induced.

Predicting the number of lifetime divisions for hematopoietic stem cells from telomere length measurements

How many times does a typical hematopoietic stem cell (HSC) divide to maintain a daily production of over 1011 blood cells over a human lifetime? It has been predicted that relatively few, slowly dividing HSCs occupy the top of the hematopoietic hierarchy. However, tracking HSCs directly is extremely challenging due to their rarity. Here, we utilize previously published data documenting the loss of telomeric DNA repeats in granulocytes, to draw inferences about HSC division rates, the timing of major changes in those rates, as well as lifetime division totals. Our method uses segmented regression to identify the best candidate representations of the telomere length data. Our method predicts that, on average, an HSC divides 56 times over an 85-year lifespan (with lower and upper bounds of 36 and 120, respectively), with half of these divisions during the first 24 years of life.

Telomeres, aging, and cancer: the big picture

The role of telomeres in human health and disease is yet to be fully understood. The limitations of mouse models for the study of human telomere biology and difficulties in accurately measuring the length of telomere repeats in chromosomes and cells have diverted attention from many important and relevant observations. The goal of this perspective is to summarize some of these observations and to discuss the antagonistic role of telomere loss in aging and cancer in the context of developmental biology, cell turnover, and evolution. It is proposed that both damage to DNA and replicative loss of telomeric DNA contribute to aging in humans, with the differences in leukocyte telomere length between humans being linked to the risk of developing specific diseases. These ideas are captured in the Telomere Erosion in Disposable Soma theory of aging proposed herein.

Frustrated folding of guanine quadruplexes in telomeric DNA

Human chromosomes terminate in long, single-stranded, DNA overhangs of the repetitive sequence (TTAGGG)n. Sets of four adjacent TTAGGG repeats can fold into guanine quadruplexes (GQ), four-stranded structures that are implicated in telomere maintenance and cell immortalization and are targets in cancer therapy. Isolated GQs have been studied in detail, however much less is known about folding in long repeat sequences. Such chains adopt an enormous number of configurations containing various arrangements of GQs and unfolded gaps, leading to a highly frustrated energy landscape. To better understand this phenomenon, we used mutagenesis, thermal melting, and global analysis to determine stability, kinetic, and cooperativity parameters for GQ folding within chains containing 8–12 TTAGGG repeats. We then used these parameters to simulate the folding of 32-repeat chains, more representative of intact telomeres. We found that a combination of folding frustration and negative cooperativity between adjacent GQs increases TTAGGG unfolding by up to 40-fold, providing an abundance of unfolded gaps that are potential binding sites for telomeric proteins. This effect was most pronounced at the chain termini, which could promote telomere extension by telomerase. We conclude that folding frustration is an important and largely overlooked factor controlling the structure of telomeric DNA.

Stationary Distribution of Telomere Lengths in Cells with Telomere Length Maintenance and its Parametric Inference

Telomeres are nucleotide caps located at the ends of each eukaryotic chromosome. Under normal physiological conditions as well as in culture, they shorten during each DNA replication round. Short telomeres initiate a proliferative arrest of cells termed 'replicative senescence'. However, cancer cells possessing limitless replication potential can avoid senescence by the telomere maintenance mechanism, which offsets telomeric loss. Therefore, cancer cells have sufficiently long telomeres even though their lengths are significantly shorter than their normal counterparts. This implies that the attrition and elongation rates play crucial roles in deciding whether and when cells ultimately become carcinogenic. In this research, we propose a concise mathematical model that shows the shortest telomere length at each cell division and prove mathematical conditions related to the attrition and elongation rates, which are necessary and sufficient for the existence of stationary distribution of telomere lengths. Moreover, we estimate the parameters of the telomere length maintenance process based on frequentist and Bayesian approaches. This study expands our knowledge of the mathematical relationship between the telomere attrition and elongation rates in cancer cells, which is important because the telomere length dynamics is a useful biomarker of cancer diagnosis and prognosis.

Senescence of Normal Human Fibroblasts Relates to the Expression of Ionotropic Glutamate Receptor GluR6/Grik2

BACKGROUND/AIM

Glutamate receptor GRIK2, previously designated as GluR6, is best described in neuronal cells. However, its biological relevance in non-neuronal cells is not well understood. We have investigated the expression of this important protein in normal human fibroblasts as a function of cell proliferation.

MATERIALS AND METHODS

We introduced expression constructs of all five isoforms (A-E) of GRIK2 in normal human fibroblasts and investigated the cells for the presence and localization of GRIK2, as well as for cell proliferation and senescence over a period of 24 days.

RESULTS

The expression of GRIK2-A isoform led to immediate cessation of cell proliferation. However, the cell numbers increased by 1.5- to 9.0-fold in 24 days upon transfection with B, C, D and E isoforms, after which they entered a state of senescence. The decreased proliferation was reflected by incorporation of BrdU in only 2-8% of transfected cells even after culturing them for 16 days.

CONCLUSION

Our results are indicative of an association between GRIK2 and aging of fibroblasts.

Replicative cellular age distributions in compartmentalized tissues

The cellular age distribution of hierarchically organized tissues can reveal important insights into the dynamics of cell differentiation and self-renewal and associated cancer risks. Here, we examine the effect of progenitor compartments with varying differentiation and self-renewal capacities on the resulting observable distributions of replicative cellular ages. We find that strongly amplifying progenitor compartments, i.e. compartments with high self-renewal capacities, substantially broaden the age distributions which become skewed towards younger cells with a long tail of few old cells. For several of these strongly amplifying compartments, the age distribution becomes virtually independent of the influx from the stem cell compartment. By contrast, if tissues are organized into many downstream compartments with low self-renewal capacity, the shape of the replicative cell distribution in more differentiated compartments is dominated by stem cell dynamics with little added variation. In the limiting case of a strict binary differentiation tree without self-renewal, the shape of the output distribution becomes indistinguishable from that of the input distribution. Our results suggest that a comparison of cellular age distributions between healthy and cancerous tissues may inform about dynamical changes within the hierarchical tissue structure, i.e. an acquired increased self-renewal capacity in certain tumours. Furthermore, we compare our theoretical results to telomere length distributions in granulocyte populations of 10 healthy individuals across different ages, highlighting that our theoretical expectations agree with experimental observations.

Mathematical model and computer simulations of telomere loss

The molecular mechanisms that control the limited number of human cell divisions has occupied researchers ever since its first description in 1961. There is evidence that this limited growth capacity, referred to as cellular or replicative senescence, is the basis for organismal ageing. Numerous studies point to the molecular mechanisms of telomere involvement in this phenomenon. A hallmark of cell senescence is high stochasticity where individual cells enter senescence in a completely random and stochastic fashion. Therefore, mathematical modelling and computational simulations of telomere dynamics are often used to explain this stochastic nature of cell ageing. Models published thus far were based on the molecular mechanisms of telomere biology and how they dictate the dynamics of cell culture proliferation. In the present work we propose an advanced model of telomere controlled cell senescence based on abrupt telomere shortening, thus explaining some important but thus far overlooked aspects of cell senescence. We test our theory by simulating the proliferative potential and two-sister experiment originally conducted by Smith and Whitney in 1980.

Divergent Aging of Isogenic Yeast Cells Revealed through Single-Cell Phenotypic Dynamics

Although genetic mutations that alter organisms' average lifespans have been identified in aging research, our understanding of the dynamic changes during aging remains limited. Here, we integrate single-cell imaging, microfluidics, and computational modeling to investigate phenotypic divergence and cellular heterogeneity during replicative aging of single S. cerevisiae cells. Specifically, we find that isogenic cells diverge early in life toward one of two aging paths, which are characterized by distinct age-associated phenotypes. We captured the dynamics of single cells along the paths with a stochastic discrete-state model, which accurately predicts both the measured heterogeneity and the lifespan of cells on each path within a cell population. Our analysis suggests that genetic and environmental factors influence both a cell's choice of paths and the kinetics of paths themselves. Given that these factors are highly conserved throughout eukaryotes, divergent aging might represent a general scheme in cellular aging of other organisms.

Stochastic Telomere Shortening and the Route to Limitless Replicative Potential

In human tissues, the replicative potential of stem cells is limited by the shortening of telomere, limitless replicative potential is a hallmark of cancer. Telomere length changes stochastically during cell division mainly due to the competition between the end replication problem and telomerase, short telomere can lead to replicative senescence and cell apoptosis. Here, we investigate how stochastic changes of telomere length in individual cells may affect the population dynamics of clonal growth. We established a computational model that couples telomerase-regulated stochastic telomere length changes with the replicative potential of clones. Model simulations reveal qualitative dependence of clone proliferation potential with activities of telomerase; mutations in cells to alter the activities of telomerase and its inhibitors can induce abnormal tissue growth and lead to limitless replicative potential.

Redundancy principle and the role of extreme statistics in molecular and cellular biology

The paradigm of chemical activation rates in cellular biology has been shifted from the mean arrival time of a single particle to the mean of the first among many particles to arrive at a small activation site. The activation rate is set by extremely rare events, which have drastically different time scales from the mean times between activations, and depends on different structural parameters. This shift calls for reconsideration of physical processes used in deterministic and stochastic modeling of chemical reactions that are based on the traditional forward rate, especially for fast activation processes in living cells. Consequently, the biological activation time is not necessarily exponentially distributed. We review here the physical models, the mathematical analysis and the new paradigm of setting the scale to be the shortest time for activation that clarifies the role of population redundancy in selecting and accelerating transient cellular search processes. We provide examples in cellular transduction, gene activation, cell senescence activation or spermatozoa selection during fertilization, where the rate depends on numbers. We conclude that the statistics of the minimal time to activation set kinetic laws in biology, which can be very different from the ones associated to average times.

Telomere shortening in breast cancer cells (MCF7) under treatment with low doses of the benzylisoquinoline alkaloid chelidonine

Telomeres, the specialized dynamic structures at chromosome ends, regularly shrink with every replication. Thus, they function as an internal molecular clock counting down the number of cell divisions. However, most cancer cells escape this limitation by activating telomerase, which can maintain telomere length. Previous studies showed that the benzylisoquinoline alkaloid chelidonine stimulates multiple modes of cell death and strongly down-regulates telomerase. It is still unknown if down-regulation of telomerase by chelidonine boosts substantial telomere shortening. The breast cancer cell line MCF7 was sequentially treated with very low concentrations of chelidonine over several cell passages. Telomere length and telomerase activity were measured by a monochrome multiplex quantitative PCR and a q-TRAP assay, respectively. Changes in population size and doubling time correlated well with telomerase inhibition and telomere shortening. MCF7 cell growth was arrested completely after three sequential treatments with 0.1 μM chelidonine, each ending after 48 h, while telomere length was reduced to almost 10% of the untreated control. However, treatment with 0.01 μM chelidonine did not have any apparent consequence. In addition to dose and time dependent telomerase inhibition, chelidonine changed the splicing pattern of hTERT towards non-enzyme coding isoforms of the transcript. In conclusion, telomere length and telomere stability are strongly affected by chelidonine in addition to microtubule formation.

The protective function of noncoding DNA in genome defense of eukaryotic male germ cells

Peripheral and abundant noncoding DNA has been hypothesized to protect the genome and the central protein-coding sequences against DNA damage in somatic genome. In the cytosol, invading exogenous nucleic acids may first be deactivated by small RNAs encoded by noncoding DNA via mechanisms similar to the prokaryotic CRISPR-Cas system. In the nucleus, the radicals generated by radiation in the cytosol, radiation energy and invading exogenous nucleic acids are absorbed, blocked and/or reduced by peripheral heterochromatin, and damaged DNA in heterochromatin is removed and excluded from the nucleus to the cytoplasm through nuclear pore complexes. To further strengthen the hypothesis, this review summarizes the experimental evidence supporting the protective function of noncoding DNA in the genome of male germ cells. Based on these data, this review provides evidence supporting the protective role of noncoding DNA in the genome defense of sperm genome through similar mechanisms to those of the somatic genome.

Telomerase activity and telomere on stem progeny senescence

The end of linear chromosomes is formed of a special nucleoprotein heterochromatin structure with repetitive TTAGGG sequences called telomere. Telomere length is regulated by a special enzyme called telomerase, a specific DNA polymerase that adds new telomeric sequences to the chromosome ends. Telomerase consists of two parts; the central protein part and the accessory part which is a RNA component transported by the central part. Regulation of telomere length by this enzyme is a multi-stage process. Telomere length elongation is strongly influenced by the level of telomerase and has a strong correlation with the activity of telomerase enzyme. Human Telomerase Reverse Transcriptase (hTERT) gene expression plays an important role in maintaining telomere length and high proliferative property of cells. Except a low activity of telomerase enzyme in hematopoietic and few types of stem cells, most of somatic cells didn't showed telomerase activity. Moreover, cytokines are secretory proteins that control many aspects of hematopoiesis, especially immune responses and inflammation. Also, the induction of hTERT gene expression by cytokines is organized through the PI3K/AKT and NF/kB signaling pathways. In this review we have tried to talk about effects of immune cell cytokines on telomerase expression/telomere length and the induction of telomerase expression by cytokines.

A population genetics perspective on the determinants of intra-tumor heterogeneity

Cancer results from the acquisition of somatic alterations in a microevolutionary process that typically occurs over many years, much of which is occult. Understanding the evolutionary dynamics that are operative at different stages of progression in individual tumors might inform the earlier detection, diagnosis, and treatment of cancer. Although these processes cannot be directly observed, the resultant spatiotemporal patterns of genetic variation amongst tumor cells encode their evolutionary histories. Such intra-tumor heterogeneity is pervasive not only at the genomic level, but also at the transcriptomic, phenotypic, and cellular levels. Given the implications for precision medicine, the accurate quantification of heterogeneity within and between tumors has become a major focus of current research. In this review, we provide a population genetics perspective on the determinants of intra-tumor heterogeneity and approaches to quantify genetic diversity. We summarize evidence for different modes of evolution based on recent cancer genome sequencing studies and discuss emerging evolutionary strategies to therapeutically exploit tumor heterogeneity. This article is part of a Special Issue entitled: Evolutionary principles - heterogeneity in cancer?, edited by Dr. Robert A. Gatenby.

Modelling the molecular mechanisms of aging

The aging process is driven at the cellular level by random molecular damage that slowly accumulates with age. Although cells possess mechanisms to repair or remove damage, they are not 100% efficient and their efficiency declines with age. There are many molecular mechanisms involved and exogenous factors such as stress also contribute to the aging process. The complexity of the aging process has stimulated the use of computational modelling in order to increase our understanding of the system, test hypotheses and make testable predictions. As many different mechanisms are involved, a wide range of models have been developed. This paper gives an overview of the types of models that have been developed, the range of tools used, modelling standards and discusses many specific examples of models that have been grouped according to the main mechanisms that they address. We conclude by discussing the opportunities and challenges for future modelling in this field.

Semen preparation methods and sperm telomere length: density gradient centrifugation versus the swim up procedure

Previous studies have shown that both density gradient centrifugation (DGC) and swim up (SU) procedures can select spermatozoa with longer telomeres for assisted reproduction techniques (ART). However, it is unknown which approach is more effective. The aim of the present study was to compare the effects of these two methods on sperm telomere length (STL). A total of 150 normozoospermic subjects were recruited. STL, DNA fragmentation index (DFI), reactive oxygen species (ROS) content and progressive motility of semen samples were detected before and after the procedures of DGC and SU. When compared to raw semen, the average length of sperm telomeres was significantly longer after the two sperm preparation methods. However, no significant difference was found between the DGC and SU procedures. We also found that semen prepared by the two methods had lower DNA fragmentation, ROS content and sperm progressive motility. However, no significant difference was found in those parameters between the two procedures. This is the first study that compares the effects of the DGC and SU procedures on STL, and the results show that both methods can recover a sperm population with longer STL and better DNA integrity for ART.

Effects of initial telomere length distribution on senescence onset and heterogeneity

Replicative senescence, induced by telomere shortening, exhibits considerable asynchrony and heterogeneity, the origins of which remain unclear. Here, we formally study how telomere shortening mechanisms impact on senescence kinetics and define two regimes of senescence, depending on the initial telomere length variance. We provide analytical solutions to the model, highlighting a non-linear relationship between senescence onset and initial telomere length distribution. This study reveals the complexity of the collective behavior of telomeres as they shorten, leading to senescence heterogeneity.

•

New stochastic model of telomere shortening and elongation by telomerase.

•

Analysis of the distribution of telomere length in the presence of telomerase.

•

Senescence onset expressed as a function of the initial telomere length distribution.

•

The initial telomere length distribution variance defines two senescence regimes.

Inhibition of YAP signaling contributes to senescence of hepatic stellate cells induced by tetramethylpyrazine

Accumulating evidence indicates that hepatic stellate cells (HSCs) are the central mediators and major effectors in the development of hepatic fibrosis. It is well-known that regulation of cell proliferation and apoptosis are potential strategies to block the activation of HSCs. Recently, several studies have revealed that induction of HSC senescence could prevent and cure the liver fibrosis. In our previous work, we have demonstrated that the natural product tetramethylpyrazine (TMP) could inhibit the activation of HSCs and ameliorate hepatic fibrosis. The aim of this study was to identify a new role of TMP in the regulation of activated HSC senescence and to elucidate the underlying mechanisms. In this study, our data showed that TMP could promote HSC senescence in vivo and in vitro. Moreover, TMP affected the cell cycle and telomerase activity. We further demonstrated that P53 siRNA or P53 pharmacological inhibitor PFT-α abrogated the TMP-induced HSC senescence in vitro. Meanwhile, similar results were obtained in vivo. Further studies indicated that TMP promoted the expression of P53 through a YAP inhibition-dependent mechanism. Moreover, silencing YAP enhanced TMP induction of activated HSC senescence. Collectively, our results suggested that TMP inhibited the activation of HSCs by inducing senescence and had therapeutic implication for the treatment of liver fibrosis.

Activation of Nrf2 contributes to the protective effect of Exendin-4 against angiotensin II-induced vascular smooth muscle cell senescence

Oxidative stress and impaired antioxidant defense are believed to be contributors to the cardiovascular aging process. The transcription factor nuclear factor-E2-related factor 2 (Nrf2) plays a key role in orchestrating cellular antioxidant defenses and maintaining redox homeostasis. Our previous study showed that Exendin-4, a glucagon-like peptide-1 analog, alleviates angiotensin II (ANG II)-induced vascular smooth muscle cell (VSMC) senescence by inhibiting Rac1 activation via cAMP/PKA (Zhao L, Li AQ, Zhou TF, Zhang MQ, Qin XM. Am J Physiol Cell Physiol 307: C1130-C1141, 2014). The objective of this study is to investigate if Nrf2 mediates the antisenescent effect of Exendin-4 in ANG II-induced VSMCs. Here we report that Exendin-4 triggered Nrf2 nuclear translocation, a downstream target of cAMP-responsive element-binding protein (CREB) and expressions of antioxidant genes heme oxygenase-1 (HO-1) and NAD(P)H quinone oxidoreductase-1 (NQO-1) in a dose- and time-dependent manner. In addition, knock-down of Nrf2 attenuated the inhibitory effects of Exendin-4 on ANG II-induced superoxidant generation and VSMC senescence. PKA/CREB pathway participated in the upregulations of HO-1 and NQO-1 induced by Exendin-4. Notably, our study revealed that Exendin-4 dose-dependently increased the acetylation of Nrf2 and the recruitment of transcriptional coactivator CREB binding protein (CBP) to Nrf2. The Exendin-4-induced Nrf2 transactivation was diminished in the presence of CBP small interfering RNA. Microscope imaging of Nrf2, as well as immunoblotting for Nrf2, showed that the Exendin-4-evoked Nrf2 acetylation favored its nuclear retention. Importantly, CBP silencing attenuated the suppressing effects of Exendin-4 on ANG II-induced VSMC senescence and superoxidant production. In conclusion, these results provide a mechanistic insight into how Nrf2 signaling mediates the antisenescent and antioxidative effects induced by Exendin-4 in VSMCs.

Aging and computational systems biology

Aging research is undergoing a paradigm shift, which has led to new and innovative methods of exploring this complex phenomenon. The systems biology approach endeavors to understand biological systems in a holistic manner, by taking account of intrinsic interactions, while also attempting to account for the impact of external inputs, such as diet. A key technique employed in systems biology is computational modeling, which involves mathematically describing and simulating the dynamics of biological systems. Although a large number of computational models have been developed in recent years, these models have focused on various discrete components of the aging process, and to date no model has succeeded in completely representing the full scope of aging. Combining existing models or developing new models may help to address this need and in so doing could help achieve an improved understanding of the intrinsic mechanisms which underpin aging.

Quantifying replicative senescence as a tumor suppressor pathway and a target for cancer therapy

To study quantitatively replicative senescence as a tumor suppressor mechanism, we investigate the distribution of a growing clonal cell population restricted by Hayflick’s limit. We find that in the biologically relevant range of parameters, if the imbalance between cell division and death is moderate or low (high death-to-birth ratio), senescence offers significant protection against cancer by halting abnormal cell proliferation at early pre-diagnostic stages of tumor development. We also find that by the time tumors are typically detected, there is a high probability that telomerase is activated, even if the cell of origin was telomerase negative. Hence, the fact that most cancers are positive for telomerase is not necessarily an indication that cancer originated in a telomerase positive cell. Finally, we discuss how the population dynamics of cells can determine the outcomes of anti-telomerase cancer therapies, and provide guidelines on how the model could potentially be applied to develop clinically useful tools to predict the response to treatment by telomerase inhibitors in individual patients.

Preventing clonal evolutionary processes in cancer: Insights from mathematical models

Clonal evolutionary processes can drive pathogenesis in human diseases, with cancer being a prominent example. To prevent or treat cancer, mechanisms that can potentially interfere with clonal evolutionary processes need to be understood better. Mathematical modeling is an important research tool that plays an ever-increasing role in cancer research. This paper discusses how mathematical models can be useful to gain insights into mechanisms that can prevent disease initiation, help analyze treatment responses, and aid in the design of treatment strategies to combat the emergence of drug-resistant cells. The discussion will be done in the context of specific examples. Among defense mechanisms, we explore how replicative limits and cellular senescence induced by telomere shortening can influence the emergence and evolution of tumors. Among treatment approaches, we consider the targeted treatment of chronic lymphocytic leukemia (CLL) with tyrosine kinase inhibitors. We illustrate how basic evolutionary mathematical models have the potential to make patient-specific predictions about disease and treatment outcome, and argue that evolutionary models could become important clinical tools in the field of personalized medicine.

A multicompartment mathematical model of cancer stem cell-driven tumor growth dynamics

Tumors are appreciated to be an intrinsically heterogeneous population of cells with varying proliferation capacities and tumorigenic potentials. As a central tenet of the so-called cancer stem cell hypothesis, most cancer cells have only a limited lifespan, and thus cannot initiate or reinitiate tumors. Longevity and clonogenicity are properties unique to the subpopulation of cancer stem cells. To understand the implications of the population structure suggested by this hypothesis--a hierarchy consisting of cancer stem cells and progeny non-stem cancer cells which experience a reduction in their remaining proliferation capacity per division--we set out to develop a mathematical model for the development of the aggregate population. We show that overall tumor progression rate during the exponential growth phase is identical to the growth rate of the cancer stem cell compartment. Tumors with identical stem cell proportions, however, can have different growth rates, dependent on the proliferation kinetics of all participating cell populations. Analysis of the model revealed that the proliferation potential of non-stem cancer cells is likely to be small to reproduce biologic observations. Furthermore, a single compartment of non-stem cancer cell population may adequately represent population growth dynamics only when the compartment proliferation rate is scaled with the generational hierarchy depth.

Mathematical model of alternative mechanism of telomere length maintenance

Biopolymer length regulation is a complex process that involves a large number of biological, chemical, and physical subprocesses acting simultaneously across multiple spatial and temporal scales. An illustrative example important for genomic stability is the length regulation of telomeres-nucleoprotein structures at the ends of linear chromosomes consisting of tandemly repeated DNA sequences and a specialized set of proteins. Maintenance of telomeres is often facilitated by the enzyme telomerase but, particularly in telomerase-free systems, the maintenance of chromosomal termini depends on alternative lengthening of telomeres (ALT) mechanisms mediated by recombination. Various linear and circular DNA structures were identified to participate in ALT, however, dynamics of the whole process is still poorly understood. We propose a chemical kinetics model of ALT with kinetic rates systematically derived from the biophysics of DNA diffusion and looping. The reaction system is reduced to a coagulation-fragmentation system by quasi-steady-state approximation. The detailed treatment of kinetic rates yields explicit formulas for expected size distributions of telomeres that demonstrate the key role played by the J factor, a quantitative measure of bending of polymers. The results are in agreement with experimental data and point out interesting phenomena: an appearance of very long telomeric circles if the total telomere density exceeds a critical value (excess mass) and a nonlinear response of the telomere size distributions to the amount of telomeric DNA in the system. The results can be of general importance for understanding dynamics of telomeres in telomerase-independent systems as this mode of telomere maintenance is similar to the situation in tumor cells lacking telomerase activity. Furthermore, due to its universality, the model may also serve as a prototype of an interaction between linear and circular DNA structures in various settings.

The amyloid hypothesis, time to move on: Amyloid is the downstream result, not cause, of Alzheimer's disease

The "amyloid hypothesis" has dominated Alzheimer research for more than 20 years, and proposes that amyloid is the toxic cause of neural/synaptic damage and dementia. If correct, decreasing the formation or removing amyloid should be therapeutic. Despite discrepancies in the proposed mechanism, and failed clinical trials, amyloid continues to be considered the cause of a degenerative cascade. Alternative hypotheses must explain three features: (i) why amyloid toxicity is not the etiology of Alzheimer's disease (AD), (ii) what alternative mechanisms cause the degeneration and dementia of AD, and (iii) why increased amyloid accumulates in the brain in AD. We propose that AD, which occurs in elderly, already vulnerable brains, with multiple age-related changes, is precipitated by impaired microvascular function, resulting primarily from decreased Notch-related angiogenesis. With impaired microvasculature, a lack of vascular endothelial-derived trophic factors and decreased cerebral blood flow cause the atrophy of neural structures. Therapeutic strategies should focus on supporting normal angiogenesis.

Cancer-associated mutations in healthy individuals: assessing the risk of carcinogenesis

Mutations associated with hematopoietic malignancies have been repeatedly identified in healthy individuals. For certain cases, such as the t(14;18) translocation and monoclonal B-cell lymphocytosis, no clear link between the presence of aberrant cells and the later development of cancer has been established. Intriguingly, longitudinal studies suggest that these abnormalities persist for long periods of time in some individuals, but in others are transient in which they disappear completely. Here, we present a mathematical model, based on cellular replication limits, that provides a possible explanation for these seemingly contradictory findings. It proposes that the transient and persistent nature of the phenotypes depends on the stage in the differentiation pathway of a given lineage in which the mutation originates. Our work suggests that cellular replication limits may not only prevent cancer by aborting clonal expansion of cells, but also by influencing the fate of altered but nonneoplastic cells in healthy tissue.

Diagnostic value of telomerase activity in peripheral blood in hepatocellular carcinoma: A Meta-analysis

AIM: To explore the diagnostic value of telomerase activity in peripheral blood in hepatocellular carcinoma (HCC).

METHODS: We searched the literature for studies on the relationship between the detection of telomerase activity in peripheral blood and HCC published in recent ten years. Meta-Disc1.4 statistical software was applied to perform comprehensive quantitative analysis. We discussed the overall diagnostic features of this assay by plotting the summary ROC curve (SROC) and calculating the area under the curve (AUC).

RESULTS: According to selection and exclusion criteria, 12 independent studies were finally selected in this meta-analysis. The random effects model was used for the meta-analysis, which showed that the comprehensive sensitivity was 0.85 (95%CI: 0.82-0.87), specificity was 0.84 (95%CI: 0.81-0.86), comprehensive diagnostic odds ratio was 28.25 (95%CI: 17.90-44.56), average value of the maximum intersection point of the sensitivity and the specificity was 0.8419, AUC was 0.9098, combined positive likelihood ratio was 5.17 (95%CI: 3.96-6.74), and negative likelihood ratio was 0.20 (95%CI: 0.14-0.28).

CONCLUSION: Telomerase activity in peripheral blood can be used as an important index to help diagnose HCC, but more high-quality, large-sized, multi-center studies are still required to evaluate its clinical value.

Computing the length of the shortest telomere in the nucleus

The telomere length can either be shortened or elongated by an enzyme called telomerase after each cell division. Interestingly, the shortest telomere is involved in controlling the ability of a cell to divide. Yet, its dynamics remains elusive. We present here a stochastic approach where we model this dynamics using a Markov jump process. We solve the forward Fokker-Planck equation to obtain the steady state distribution and the statistical moments of telomere lengths. We focus specifically on the shortest one and we estimate its length difference with the second shortest telomere. After extracting key parameters such as elongation and shortening dynamics from experimental data, we compute the length of telomeres in yeast and obtain as a possible prediction the minimum concentration of telomerase required to ensure a proper cell division.

Minimizing the risk of cancer: tissue architecture and cellular replication limits

Normal somatic cells are capable of only a limited number of divisions, which prevents unlimited cell proliferation and the onset of tumours. Cancer cells find ways to circumvent this obstacle, typically by expressing the enzyme telomerase and less often by alternative recombination strategies. Given this fundamental link between cellular replication limits and cancer, it is important to understand how a tissue's architecture affects the replicative capacity of a cell population. We define this as the average number of remaining divisions at equilibrium. The lower the replication capacity, the lower the chances to escape the replication limit during abnormal growth when a tumour develops. In this paper, we examine how the replication capacity is influenced by defining characteristics of cell lineages, such as the number of intermediate cell compartments, self-renewal capability of cells and division rates. We describe an optimal tissue architecture that minimizes the replication capacity of dividing cells and thus the risk of cancer. Interestingly, some of the features that define an optimal tissue architecture have been documented in a variety of tissues, suggesting that they may have evolved as a cancer-protecting strategy in multicellular organisms.

Derivation of putative porcine embryonic germ cells and analysis of their multi-lineage differentiation potential

Embryonic germ (EG) cells are cultured pluripotent stem cells derived from the primordial germ cells (PGCs) that migrate from the dorsal mesentery of the hindgut to the developing genital ridge. In this study, the morphology of the porcine genital ridge was assessed in embryos harvested on days 22-30 of pregnancy. PGCs from embryos at these stages were cultured to obtain porcine EG cell lines, and EG-like cells were derived from PGCs from embryos harvested on days 24-28 of pregnancy. The EG-like cells expressed Oct4, Sox2, Nanog, SSEA-3, SSEA-4 and alkaline phosphatase (AP). These cells were able to form embryoid bodies (EBs) in suspension culture and differentiate into cells representative of the three germ layers as verified by a-fetoprotein (AFP), α-smooth muscle actin (α-SMA), and Nestin expression. Spontaneous differentiation from the porcine EG-like cells of delayed passage in vitro showed that they could differentiate into epithelial-like cells, mesenchymal-like cells and neuron-like cells. In vitro directed differentiation generated osteocytes, adipocytes and a variety of neural lineage cells, as demonstrated by alizarin red staining, oil red O staining, and immunofluorescence for neuronal class Ⅲ β-tubulin (Tuj1), glial fibrillary protein (GFAP) and galactosylceramidase (GALC), respectively. These results indicate that porcine EG-like cells have the potential for multi-lineage differentiation and are useful for basic porcine stem cell research.

The length of the shortest telomere as the major determinant of the onset of replicative senescence

The absence of telomerase in many eukaryotes leads to the gradual shortening of telomeres, causing replicative senescence. In humans, this proliferation barrier constitutes a tumor suppressor mechanism and may be involved in cellular aging. Yet the heterogeneity of the senescence phenotype has hindered the understanding of its onset. Here we investigated the regulation of telomere length and its control of senescence heterogeneity. Because the length of the shortest telomeres can potentially regulate cell fate, we focus on their dynamics in Saccharomyces cerevisiae. We developed a stochastic model of telomere dynamics built on the protein-counting model, where an increasing number of protein-bound telomeric repeats shift telomeres into a nonextendable state by telomerase. Using numerical simulations, we found that the length of the shortest telomere is well separated from the length of the others, suggesting a prominent role in triggering senescence. We evaluated this possibility using classical genetic analyses of tetrads, combined with a quantitative and sensitive assay for senescence. In contrast to mitosis of telomerase-negative cells, which produces two cells with identical senescence onset, meiosis is able to segregate a determinant of senescence onset among the telomerase-negative spores. The frequency of such segregation is in accordance with this determinant being the length of the shortest telomere. Taken together, our results substantiate the length of the shortest telomere as being the key genetic marker determining senescence onset in S. cerevisiae.

Modelling the regulation of telomere length: the effects of telomerase and G-quadruplex stabilising drugs

Telomeres are guanine-rich sequences at the end of chromosomes which shorten during each replication event and trigger cell cycle arrest and/or controlled death (apoptosis) when reaching a threshold length. The enzyme telomerase replenishes the ends of telomeres and thus prolongs the life span of cells, but also causes cellular immortalisation in human cancer. G-quadruplex (G4) stabilising drugs are a potential anticancer treatment which work by changing the molecular structure of telomeres to inhibit the activity of telomerase. We investigate the dynamics of telomere length in different conformational states, namely t-loops, G-quadruplex structures and those being elongated by telomerase. By formulating deterministic differential equation models we study the effects of various levels of both telomerase and concentrations of a G4-stabilising drug on the distribution of telomere lengths, and analyse how these effects evolve over large numbers of cell generations. As well as calculating numerical solutions, we use quasicontinuum methods to approximate the behaviour of the system over time, and predict the shape of the telomere length distribution. We find those telomerase and G4-concentrations where telomere length maintenance is successfully regulated. Excessively high levels of telomerase lead to continuous telomere lengthening, whereas large concentrations of the drug lead to progressive telomere erosion. Furthermore, our models predict a positively skewed distribution of telomere lengths, that is, telomeres accumulate over lengths shorter than the mean telomere length at equilibrium. Our model results for telomere length distributions of telomerase-positive cells in drug-free assays are in good agreement with the limited amount of experimental data available.

SNP rs1533428 at 2p16.3 as a marker for late-onset primary open-angle glaucoma

PURPOSE

To investigate the associations between gene variants in cholesterol 24S-hydroxylase (CYP46A1), LIM homeobox transcription factor 1-beta (LMX1B), plexin domain containing 2 (PLXDC2), toll-like receptor 4 (TLR4), transmembrane and tetratricopeptide repeat containing 2 (TMTC2), zona pellucida glycoprotein 4 (ZP4), chromosome 2p16.3, and primary open-angle glaucoma (POAG).

METHODS

We studied 462 POAG patients and 577 controls from three cohorts (Hong Kong, Shantou, and Beijing, China). Twelve single-nucleotide polymorphisms (SNPs) were genotyped in the Hong Kong cohort using TaqMan genotyping assay. Significant associations were validated in the Shantou and Beijing cohorts.

RESULTS

Association of POAG with TLR4 rs7037117, in a recessive model, was identified in the Hong Kong and Shantou cohorts (both southern Chinese, p(rec)=0.0019) but not the Beijing cohort (northern Chinese). rs1533428 at chromosome 2p16.3 showed a consistent trend of age-specific association in all three cohorts. Genotypes TT + CT conferred a 2.16 fold of significantly increased risk to late-onset POAG (p(dom)=0.00025), but no significant risk to POAG of younger ages of onset in the combined cohort. A joint effect was found between rs7037117 and rs1533428, with carriers of both higher-risk genotypes having a 4.53 fold of increased disease risk (p=0.00028).

CONCLUSIONS

Our study reveals discrepant association patterns of 12 candidate SNPs in 7 genes/loci with POAG in Chinese, provides positive replications for POAG markers rs1533428 at 2p16.3 and TLR4 rs7037117, and suggests that rs1533428 is a putative risk variant for late-onset POAG. The identification of an age-specific association between rs1533428 and late-onset POAG highlights a new genotype-phenotype association in POAG. Further studies are warranted to confirm the age-specific association.

Isolation, proliferation, cytogenetic, and molecular characterization and in vitro differentiation potency of canine stem cells from foetal adnexa: a comparative study of amniotic fluid, amnion, and umbilical cord matrix

The possibility to isolate canine mesenchymal stem cells (MSCs) from foetal adnexa is interesting since several canine genetic disorders are reported to resemble similar dysfunctions in humans. In this study, we successfully isolated, cytogenetically and molecularly characterized, and followed the differentiation potency of canine MSCs from foetal adnexa, such as amniotic fluid (AF), amniotic membrane (AM), and umbilical cord matrix (UCM). In the three types of cell lines, the morphology of proliferating cells typically appeared fibroblast-like, and the population doubling time (DT) significantly increased with passage number. For AF- and AM-MSCs, cell viability did not change with passages. In UCM-MSCs, cell viability remained at approximately constant levels up to P6 and significantly decreased from P7 (P < 0.05). Amnion and UCM-MSCs expressed embryonic and MSC markers, such as Oct-4 CD44, CD184, and CD29, whereas AF-MSCs expressed Oct-4, CD44. Expression of the hematopoietic markers CD34 and CD45 was not found. Dog leucocyte antigens (DLA-DRA1 and DLA-79) were expressed only in AF-MSCs at P1. Isolated cells of the three cell lines at P3 showed multipotent capacity, and differentiated in vitro into neurocyte, adipocyte, osteocyte, and chondrocyte, as demonstrated by specific stains and expression of molecular markers. Cells at P4 showed normal chromosomal number, structure, and telomerase activity. These results demonstrate that, in dog, MSCs can be successfully isolated from foetal adnexa and grown in vitro. Their proven stemness and chromosomal stability indicated that MSCs could be used as a model to study stem cell biology and have an application in therapeutic programs.