Stem cell factor and erythropoietin-independent production of cultured reticulocytes.

Cultured reticulocytes can supplement transfusion needs and offer promise for drug delivery and immune tolerization. They can be produced from induced pluripotent stem cells (iPSCs), but the 45-day culture time and cytokine costs make large-scale production prohibitive. To overcome these limitations, we have generated IPSCs that express constitutive SCF receptor and jak2 adaptor alleles. We show that iPSC lines carrying these alleles can differentiate into self-renewing erythroblast (SRE) that can proliferate for up to 70 cell-doubling in a cost-effective, chemically-defined, albumin- and cytokine-free medium. These kitjak2 SREs retain the ability to enucleate at a high rate up to senescence. Kitjak2 derived cultured reticulocytes should be safe for transfusion because they can be irradiated to eliminate residual nucleated cells. The kitjak2 cells express blood group 0 and test negative for RhD and other clinically significant RBCs antigens and have sufficient proliferation capacity to meet global RBC needs.

Translesion polymerase eta both facilitates DNA replication and promotes increased human genetic variation at common fragile sites.

Common fragile sites (CFSs) are difficult-to-replicate genomic regions that form gaps and breaks on metaphase chromosomes under replication stress. They are hotspots for chromosomal instability in cancer. Repetitive sequences located at CFS loci are inefficiently copied by replicative DNA polymerase (Pol) delta. However, translesion synthesis Pol eta has been shown to efficiently polymerize CFS-associated repetitive sequences in vitro and facilitate CFS stability by a mechanism that is not fully understood. Here, by locus-specific, single-molecule replication analysis, we identified a crucial role for Pol eta (encoded by the gene POLH) in the in vivo replication of CFSs, even without exogenous stress. We find that Pol eta deficiency induces replication pausing, increases initiation events, and alters the direction of replication-fork progression at CFS-FRA16D in both lymphoblasts and fibroblasts. Furthermore, certain replication pause sites at CFS-FRA16D were associated with the presence of non-B DNA-forming motifs, implying that non-B DNA structures could increase replication hindrance in the absence of Pol eta. Further, in Pol eta-deficient fibroblasts, there was an increase in fork pausing at fibroblast-specific CFSs. Importantly, while not all pause sites were associated with non-B DNA structures, they were embedded within regions of increased genetic variation in the healthy human population, with mutational spectra consistent with Pol eta activity. From these findings, we propose that Pol eta replicating through CFSs may result in genetic variations found in the human population at these sites.

Allele-specific genome-wide profiling in human primary erythroblasts reveal replication program organization.

We have developed a new approach to characterize allele-specific timing of DNA replication genome-wide in human primary basophilic erythroblasts. We show that the two chromosome homologs replicate at the same time in about 88% of the genome and that large structural variants are preferentially associated with asynchronous replication. We identified about 600 megabase-sized asynchronously replicated domains in two tested individuals. The longest asynchronously replicated domains are enriched in imprinted genes suggesting that structural variants and parental imprinting are two causes of replication asynchrony in the human genome. Biased chromosome X inactivation in one of the two individuals tested was another source of detectable replication asynchrony. Analysis of high-resolution TimEX profiles revealed small variations termed timing ripples, which were undetected in previous, lower resolution analyses. Timing ripples reflect highly reproducible, variations of the timing of replication in the 100 kb-range that exist within the well-characterized megabase-sized replication timing domains. These ripples correspond to clusters of origins of replication that we detected using novel nascent strands DNA profiling methods. Analysis of the distribution of replication origins revealed dramatic differences in initiation of replication frequencies during S phase and a strong association, in both synchronous and asynchronous regions, between origins of replication and three genomic features: G-quadruplexes, CpG Islands and transcription start sites. The frequency of initiation in asynchronous regions was similar in the two homologs. Asynchronous regions were richer in origins of replication than synchronous regions.

GenPlay Multi-Genome, a tool to compare and analyze multiple human genomes in a graphical interface.

SUMMARY: Parallel visualization of multiple individual human genomes is a complex endeavor that is rapidly gaining importance with the increasing number of personal, phased and cancer genomes that are being generated. It requires the display of variants such as SNPs, indels and structural variants that are unique to specific genomes and the introduction of multiple overlapping gaps in the reference sequence. Here, we describe GenPlay Multi-Genome, an application specifically written to visualize and analyze multiple human genomes in parallel. GenPlay Multi-Genome is ideally suited for the comparison of allele-specific expression and functional genomic data obtained from multiple phased genomes in a graphical interface with access to multiple-track operation. It also allows the analysis of data that have been aligned to custom genomes rather than to a standard reference and can be used as a variant calling format file browser and as a tool to compare different genome assembly, such as hg19 and hg38. AVAILABILITY AND IMPLEMENTATION: GenPlay is available under the GNU public license (GPL-3) from http://genplay.einstein.yu.edu. The source code is available at https://github.com/JulienLajugie/GenPlay.

The exosome complex establishes a barricade to erythroid maturation.

Complex genetic networks control hematopoietic stem cell differentiation into progenitors that give rise to billions of erythrocytes daily. Previously, we described a role for the master regulator of erythropoiesis, GATA-1, in inducing genes encoding components of the autophagy machinery. In this context, the Forkhead transcription factor, Foxo3, amplified GATA-1-mediated transcriptional activation. To determine the scope of the GATA-1/Foxo3 cooperativity, and to develop functional insights, we analyzed the GATA-1/Foxo3-dependent transcriptome in erythroid cells. GATA-1/Foxo3 repressed expression of Exosc8, a pivotal component of the exosome complex, which mediates RNA surveillance and epigenetic regulation. Strikingly, downregulating Exosc8, or additional exosome complex components, in primary erythroid precursor cells induced erythroid cell maturation. Our results demonstrate a new mode of controlling erythropoiesis in which multiple components of the exosome complex are endogenous suppressors of the erythroid developmental program.

High-resolution mapping of h1 linker histone variants in embryonic stem cells.

H1 linker histones facilitate higher-order chromatin folding and are essential for mammalian development. To achieve high-resolution mapping of H1 variants H1d and H1c in embryonic stem cells (ESCs), we have established a knock-in system and shown that the N-terminally tagged H1 proteins are functionally interchangeable to their endogenous counterparts in vivo. H1d and H1c are depleted from GC- and gene-rich regions and active promoters, inversely correlated with H3K4me3, but positively correlated with H3K9me3 and associated with characteristic sequence features. Surprisingly, both H1d and H1c are significantly enriched at major satellites, which display increased nucleosome spacing compared with bulk chromatin. While also depleted at active promoters and enriched at major satellites, overexpressed H1(0) displays differential binding patterns in specific repetitive sequences compared with H1d and H1c. Depletion of H1c, H1d, and H1e causes pericentric chromocenter clustering and de-repression of major satellites. These results integrate the localization of an understudied type of chromatin proteins, namely the H1 variants, into the epigenome map of mouse ESCs, and we identify significant changes at pericentric heterochromatin upon depletion of this epigenetic mark.

Erythropoiesis from human embryonic stem cells through erythropoietin-independent AKT signaling.

Unlimited self renewal capacity and differentiation potential make human pluripotent stem cells (PSC) a promising source for the ex vivo manufacture of red blood cells (RBCs) for safe transfusion. Current methods to induce erythropoiesis from PSC suffer from low yields of RBCs, most of which are immature and contain embryonic and fetal rather than adult hemoglobins. We have previously shown that homodimerization of the intracellular component of MPL (ic-MPL) induces erythropoiesis from human cord blood progenitors. The goal of this study was to investigate the potential of ic-MPL dimerization to induce erythropoiesis from human embryonic stem cells (hESCs) and to identify the signaling pathways activated by this strategy. We present here the evidence that ic-MPL dimerization induces erythropoietin (EPO)-independent erythroid differentiation from hESC by inducing the generation of erythroid progenitors and by promoting more efficient erythroid maturation with increased RBC enucleation as well as increased gamma:epsilon globin ratio and production of beta-globin protein. ic-MPL dimerization is significantly more potent than EPO in inducing erythropoiesis, and its effect is additive to EPO. Signaling studies show that dimerization of ic-MPL, unlike stimulation of the wild type MPL receptor, activates AKT in the absence of JAK2/STAT5 signaling. AKT activation upregulates GATA-1 and FOXO3 transcriptional pathways with resulting inhibition of apoptosis, modulation of cell cycle, and enhanced maturation of erythroid cells. These findings open up potential new targets for the generation of therapeutically relevant RBC products from hPSC.

Zinc-finger nuclease-mediated correction of α-thalassemia in iPS cells.

Induced pluripotent stem (iPS) cell technology holds vast promises for a cure to the hemoglobinopathies. Constructs and methods to safely insert therapeutic genes to correct the genetic defect need to be developed. Site-specific insertion is a very attractive method for gene therapy because the risks of insertional mutagenesis are eliminated provided that a "safe harbor" is identified, and because a single set of validated constructs can be used to correct a large variety of mutations simplifying eventual clinical use. We report here the correction of α-thalassemia major hydrops fetalis in transgene-free iPS cells using zinc finger-mediated insertion of a globin transgene in the AAVS1 site on human chromosome 19. Homozygous insertion of the best of the 4 constructs tested led to complete correction of globin chain imbalance in erythroid cells differentiated from the corrected iPS cells.

Predictable dynamic program of timing of DNA replication in human cells.

The organization of mammalian DNA replication is poorly understood. We have produced high-resolution dynamic maps of the timing of replication in human erythroid, mesenchymal, and embryonic stem (ES) cells using TimEX, a method that relies on gaussian convolution of massive, highly redundant determinations of DNA copy-number variations during S phase to produce replication timing profiles. We first obtained timing maps of 3% of the genome using high-density oligonucleotide tiling arrays and then extended the TimEX method genome-wide using massively parallel sequencing. We show that in untransformed human cells, timing of replication is highly regulated and highly synchronous, and that many genomic segments are replicated in temporal transition regions devoid of initiation, where replication forks progress unidirectionally from origins that can be hundreds of kilobases away. Absence of initiation in one transition region is shown at the molecular level by single molecule analysis of replicated DNA (SMARD). Comparison of ES and erythroid cells replication patterns revealed that these cells replicate about 20% of their genome in different quarters of S phase. Importantly, we detected a strong inverse relationship between timing of replication and distance to the closest expressed gene. This relationship can be used to predict tissue-specific timing of replication profiles from expression data and genomic annotations. We also provide evidence that early origins of replication are preferentially located near highly expressed genes, that mid-firing origins are located near moderately expressed genes, and that late-firing origins are located far from genes.

A transgenic mouse model expressing exclusively human hemoglobin E: indications of a mild oxidative stress.

Hemoglobin (Hb) E (ß26 Glu→Lys) is the most common abnormal hemoglobin (Hb) variant in the world. Homozygotes for HbE are mildly thalassemic as a result of the alternate splice mutation and present with a benign clinical picture (microcytic and mildly anemic) with rare clinical symptoms. Given that the human red blood cell (RBC) contains both HbE and excess α-chains along with minor hemoglobins, the consequence of HbE alone on RBC pathophysiology has not been elucidated. This becomes critical for the highly morbid ß(E)-thalassemia disease. We have generated transgenic mice exclusively expressing human HbE (HbEKO) that exhibit the known aberrant splicing of ß(E) globin mRNA, but are essentially non-thalassemic as demonstrated by RBC α/ß (human) globin chain synthesis. These mice exhibit hematological characteristics similar to presentations in human EE individuals: microcytic RBC with low MCV and MCH but normal MCHC; target RBC; mild anemia with low Hb, HCT and mildly elevated reticulocyte levels and decreased osmotic fragility, indicating altered RBC surface area to volume ratio. These alterations are correlated with a mild RBC oxidative stress indicated by enhanced membrane lipid peroxidation, elevated zinc protoporphyrin levels, and by small but significant changes in cardiac function. The C57 (background) mouse and full KO mouse models expressing HbE with the presence of HbS or HbA are used as controls. In select cases, the HbA full KO mouse model is compared but found to be limited due to its RBC thalassemic characteristics. Since the HbEKO mouse RBC lacks an abundance of excess α-chains that would approximate a mouse thalassemia (or a human thalassemia), the results indicate that the observed in vivo RBC mild oxidative stress arises, at least in part, from the molecular consequences of the HbE mutation.

GenPlay, a multipurpose genome analyzer and browser.

MOTIVATION: Rapidly decreasing sequencing cost due to the emergence and improvement of massively parallel sequencing technologies has resulted in a dramatic increase in the quantity of data that needs to be analyzed. Therefore, software tools to process, visualize, analyze and integrate data produced on multiple platforms and using multiple methods are needed. RESULTS: GenPlay is a fast, easy to use and stable tool for rapid analysis and data processing. It is written in Java and runs on all major operating systems. GenPlay recognizes a wide variety of common genomic data formats from microarray- or sequencing-based platforms and offers a library of operations (normalization, binning, smoothing) to process raw data into visualizable tracks. GenPlay displays tracks adapted to summarize gene structure, gene expression, repeat families, CPG islands, etc. as well as custom tracks to show the results of RNA-Seq, ChIP-Seq, TimEX-Seq and single nucleotide polymorphism (SNP) analysis. GenPlay can generate statistics (minimum, maximum, SD, correlation, etc.). The tools provided include Gaussian filter, peak finders, signal saturation, island finders. The software also offers graphical features such as scatter plots and bar charts to depict signal repartition. The library of operations is continuously growing based on the emerging needs. AVAILABILITY: GenPlay is an open-source project available from http://www.genplay.net. The code source of the software is available at https://genplay.einstein.yu.edu/svn/GenPlay.

Epigenetic control of chromosome-associated lncRNA genes essential for replication and stability.

ASARs are long noncoding RNA genes that control replication timing of entire human chromosomes in cis. The three known ASAR genes are located on human chromosomes 6 and 15, and are essential for chromosome integrity. To identify ASARs on all human chromosomes we utilize a set of distinctive ASAR characteristics that allow for the identification of hundreds of autosomal loci with epigenetically controlled, allele-restricted behavior in expression and replication timing of coding and noncoding genes, and is distinct from genomic imprinting. Disruption of noncoding RNA genes at five of five tested loci result in chromosome-wide delayed replication and chromosomal instability, validating their ASAR activity. In addition to the three known essential cis-acting chromosomal loci, origins, centromeres, and telomeres, we propose that all mammalian chromosomes also contain "Inactivation/Stability Centers" that display allele-restricted epigenetic regulation of protein coding and noncoding ASAR genes that are essential for replication and stability of each chromosome.

Efficient generation of lens progenitor cells and lentoid bodies from human embryonic stem cells in chemically defined conditions.

The eye lens is an encapsulated avascular organ whose function is to focus light on the retina. Lens comprises a single progenitor cell lineage in multiple states of differentiation. Disruption of lens function leading to protein aggregation and opacity results in age-onset cataract. Cataract is a complex disease involving genetic and environmental factors. Here, we report the development of a new 3-stage system that differentiates human embryonic stem cells (hESCs) into large quantities of lens progenitor-like cells and differentiated 3-dimensional lentoid bodies. Inhibition of BMP signaling by noggin triggered differentiation of hESCs toward neuroectoderm. Subsequent reactivation of BMP and activation of FGF signaling stimulated formation of lens progenitor cells marked by the expression of PAX6 and alpha-crystallins. The formation of lentoid bodies was most efficient in the presence of FGF2 and Wnt-3a, yielding approximately 1000 lentoid bodies/30-mm well. Lentoid bodies expressed and accumulated lens-specific markers including alphaA-, alphaB-, beta-, and gamma-crystallins, filensin, CP49, and MIP/aquaporin 0. Collectively, these studies identify a novel procedure to generate lens cells from hESCs that can be applied for studies of lens differentiation and cataractogenesis using induced pluripotent stem (iPS) cells derived from various cataract patients.

High-resolution genome-wide cytosine methylation profiling with simultaneous copy number analysis and optimization for limited cell numbers.

Many genome-wide assays involve the generation of a subset (or representation) of the genome following restriction enzyme digestion. The use of enzymes sensitive to cytosine methylation allows high-throughput analysis of this epigenetic regulatory process. We show that the use of a dual-adapter approach allows us to generate genomic representations that includes fragments of <200 bp in size, previously not possible when using the standard approach of using a single adapter. By expanding the representation to smaller fragments using HpaII or MspI, we increase the representation by these isoschizomers to more than 1.32 million loci in the human genome, representing 98.5% of CpG islands and 91.1% of refSeq promoters. This advance allows the development of a new, high-resolution version of our HpaII-tiny fragment Enrichment by Ligation-mediated PCR (HELP) assay to study cytosine methylation. We also show that the MspI representation generates information about copy-number variation, that the assay can be used on as little as 10 ng of DNA and that massively parallel sequencing can be used as an alternative to microarrays to read the output of the assay, making this a powerful discovery platform for studies of genomic and epigenomic abnormalities.

Clonal origin in normal adults of all blood lineages and circulating hematopoietic stem cells.

Characterization of human cells that sustain blood cell production lifelong has historically been inferred from phenotypically defined subsets of cells assayed in vitro, in transplanted immunodeficient mice, or in patients transplanted with genetically marked cells. These approaches have led to the concept of a persistent complex hierarchical process of differentiation divisions originating from a rare population of CD34+CD38-CD45RA-CD90+CD49f+ cells with an average self-renewal potential of >0.5 and an ability to produce some or all blood cell types for >1 year. However, the role of these "49f" cells in the unperturbed adult has remained poorly understood. To address this gap, somatic single-nucleotide polymorphisms (SNVs) have recently been exploited as lineage tracing markers to enumerate and characterize active hematopoietic clones in normal adults using a capture and recapture approach. We show here that the use of somatic transversions to identify somatically acquired variant alleles enabled their detection in bulk populations at frequencies of approximately 1 in 80,000 cells. We then applied this method to blood cells isolated from two normal adults (aged 31 and 53 years) over a 1- to 3-year period. The results revealed in both donors a continued clonal output of both T- and B-lymphoid cells as well as myeloid cells identified by the same unique transversions found to distinguish single 49f cells isolated from the same donors' initial blood samples. These findings provide the first evidence of a continuing hematopoietic stem cell-derived source of all mature blood cell types in normal (unperturbed) adult humans.

HbF Levels in Sickle Cell Disease Are Associated with Proportion of Circulating Hematopoietic Stem and Progenitor Cells and CC-Chemokines.

The concentration of circulating hematopoietic stem and progenitor cells has not been studied longitudinally. Here, we report that the proportions of Lin-CD34+38- hematopoietic multipotent cells (HMCs) and of Lin-CD34+CD38+ hematopoietic progenitors cells (HPCs) are highly variable between individuals but stable over long periods of time, in both healthy individuals and sickle cell disease (SCD) patients. This suggests that these proportions are regulated by genetic polymorphisms or by epigenetic mechanisms. We also report that in SCD patients treated with hydroxyurea, the proportions of circulating HMCs and HPCs show a strong positive and negative correlation with fetal hemoglobin (HbF) levels, respectively. Titration of 65 cytokines revealed that the plasma concentration of chemokines CCL2, CCL11, CCL17, CCL24, CCL27, and PDGF-BB were highly correlated with the proportion of HMCs and HPCs and that a subset of these cytokines were also correlated with HbF levels. A linear model based on four of these chemokines could explain 80% of the variability in the proportion of circulating HMCs between individuals. The proportion of circulating HMCs and HPCs and the concentration of these chemokines might therefore become useful biomarkers for HbF response to HU in SCD patients. Such markers might become increasingly clinically relevant, as alternative treatment modalities for SCD are becoming available.

Characterization of Hematopoiesis in Sickle Cell Disease by Prospective Isolation of Stem and Progenitor Cells.

The consequences of sickle cell disease (SCD) include ongoing hematopoietic stress, hemolysis, vascular damage, and effect of chronic therapies, such as blood transfusions and hydroxyurea, on hematopoietic stem and progenitor cell (HSPC) have been poorly characterized. We have quantified the frequencies of nine HSPC populations by flow cytometry in the peripheral blood of pediatric and adult patients, stratified by treatment and control cohorts. We observed broad differences between SCD patients and healthy controls. SCD is associated with 10 to 20-fold increase in CD34dim cells, a two to five-fold increase in CD34bright cells, a depletion in Megakaryocyte-Erythroid Progenitors, and an increase in hematopoietic stem cells, when compared to controls. SCD is also associated with abnormal expression of CD235a as well as high levels CD49f antigen expression. These findings were present to varying degrees in all patients with SCD, including those on chronic therapy and those who were therapy naive. HU treatment appeared to normalize many of these parameters. Chronic stress erythropoiesis and inflammation incited by SCD and HU therapy have long been suspected of causing premature aging of the hematopoietic system, and potentially increasing the risk of hematological malignancies. An important finding of this study was that the observed concentration of CD34bright cells and of all the HSPCs decreased logarithmically with time of treatment with HU. This correlation was independent of age and specific to HU treatment. Although the number of circulating HSPCs is influenced by many parameters, our findings suggest that HU treatment may decrease premature aging and hematologic malignancy risk compared to the other therapeutic modalities in SCD.

Preventing gene silencing with human replicators.

Transcriptional silencing, one of the major impediments to gene therapy in humans, is often accompanied by replication during late S-phase. We report that transcriptional silencing and late replication were prevented by DNA sequences that can initiate DNA replication (replicators). When replicators were included in silencing-prone transgenes, they did not undergo transcriptional silencing, replicated early and maintained histone acetylation patterns characteristic of euchromatin. A mutant replicator, which could not initiate replication, could not prevent gene silencing and replicated late when included in identical transgenes and inserted at identical locations. These observations suggest that replicators introduce epigenetic chromatin changes that facilitate initiation of DNA replication and affect gene silencing. Inclusion of functional replicators in gene therapy vectors may provide a tool for stabilizing gene expression patterns.