The mechanisms that control the location and timing of firing of replication origins are poorly understood. Using a novel functional genomic approach based on the analysis of SNPs and indels in phased human genomes, we observe that replication asynchrony is associated with small cumulative variations in the initiation efficiency of multiple origins between the chromosome homologues, rather than with the activation of dormant origins. Allele-specific measurements demonstrate that the presence of G-quadruplex-forming sequences does not correlate with the efficiency of initiation. Sequence analysis reveals that the origins are highly enriched in sequences with profoundly asymmetric G/C and A/T nucleotide distributions and are almost completely depleted of antiparallel triplex-forming sequences. We therefore propose that although G4-forming sequences are abundant in replication origins, an asymmetry in nucleotide distribution, which increases the propensity of origins to unwind and adopt non-B DNA structure, rather than the ability to form G4, is directly associated with origin activity.

The molecular mechanisms governing γ-globin expression in a subset of fetal hemoglobin (α2γ2: HbF) expressing red blood cells (F-cells) and the mechanisms underlying the variability of response to hydroxyurea induced γ-globin expression in the treatment of sickle cell disease are not completely understood. Here we analyzed intra-person clonal populations of basophilic erythroblasts (baso-Es) derived from bone marrow common myeloid progenitors in serum free cultures and report the level of fetal hemoglobin production in F-cells negatively correlates with expression of BCL11A, KLF1 and TAL1. We then examined the effects of hydroxyurea on these three transcription factors and conclude that a successful induction of γ-globin includes a reduction in BCL11A, KLF1 and TAL1 expression. These data suggests that expression changes in this transcription factor network modulate γ-globin expression in F-cells during steady state erythropoiesis and after induction with hydroxyurea.

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.

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.

As cell culture methods and stem cell biology have progressed, the in vitro production of cultured RBCs (cRBCs) has emerged as a viable option to produce cells for transfusion or to carry therapeutic cargoes. RBCs produced in culture can be quality-tested either by xeno-transfusion of human cells into immuno-deficient animals, or by transfusion of autologous cells in immuno-competent models. Although murine xeno-transfusion methods have improved, they must be complemented by studies in immuno-competent models. Non-human primates (NHPs) are important pre-clinical, large animal models due to their high biological and developmental similarities with humans, including their comparable hematopoietic and immune systems. Among NHPs, baboons are particularly attractive to validate cRBCs because of the wealth of data available on the characteristics of RBCs in this species that have been generated by past blood transfusion studies. We report here that we have developed a method to produce enucleated cRBCs by differentiation of baboon induced pluripotent stem cells (iPSCs). This method will enable the use of baboons to evaluate therapeutic cRBCs and generate essential pre-clinical data in an immuno-competent, large animal model. Production of the enucleated baboon cRBCs was achieved by adapting the PSC-RED protocol that we previously developed for human cells. Baboon-PSC-RED is an efficient chemically-defined method to differentiate iPSCs into cRBCs that are about 40% to 50% enucleated. PSC-RED is relatively low cost because it requires no albumin and only small amounts of recombinant transferrin.

Eukaryotic genome duplication starts at discrete sequences (replication origins) that coordinate cell cycle progression, ensure genomic stability and modulate gene expression. Origins share some sequence features, but their activity also responds to changes in transcription and cellular differentiation status.

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.

To study transcriptional regulation by distant enhancers, we devised a system of easily modified reporter plasmids for integration into single-copy targeting cassettes in clones of HuH7, a human hepatocellular carcinoma. The plasmid constructs tested transcriptional function of a 35-kb region that contained the rat albumin gene and its upstream flanking region. Expression of integrants was analyzed in two orientations, and compared to transient expression of non-integrated plasmids. Enhancers were studied in their natural positions relative to the promoter and localized by deletion. All constructs were also analyzed by transient transfection assays. In addition to the known albumin gene enhancer (E1 at -10 kb), we demonstrated two new enhancers, E2 at -13, and E4 at +1.2 kb. All three enhancers functioned in both transient assays and integrated constructs. However, chromosomal integration demonstrated several differences from transient expression. For example, analysis of E2 showed that enhancer function within the chromosome required a larger gene region than in transient assays. Another conserved region, E3 at -0.7 kb, functioned as an enhancer in transient assays but inhibited the function of E1 and E2 when chromosomally integrated. The enhancers did not show additive or synergistic behavior,an effect consistent with competition for the promoter or inhibitory interactions among enhancers. Growth arrest by serum starvation strongly stimulated the function of some integrated enhancers, consistent with the expected disruption of enhancer-promoter looping during the cell cycle.

Phased genome maps are important to understand genetic and epigenetic regulation and disease mechanisms, particularly parental imprinting defects. Phasing is also critical to assess the functional consequences of genetic variants, and to allow precise definition of haplotype blocks which is useful to understand gene-flow and genotype-phenotype association at the population level. Transmission phasing by analysis of a family quartet allows the phasing of 95% of all variants as the uniformly heterozygous positions cannot be phased. Here, we report a phasing method based on a combination of transmission analysis, physical phasing by pair-end sequencing of libraries of staggered sizes and population-based analysis. Sequencing of a healthy Caucasians quartet at 120x coverage and combination of physical and transmission phasing yielded the phased genotypes of about 99.8% of the SNPs, indels and structural variants present in the quartet, a phasing rate significantly higher than what can be achieved using any single phasing method. A false positive SNP error rate below 10*E-7 per genome and per base was obtained using a combination of filters. We provide a complete list of SNPs, indels and structural variants, an analysis of haplotype block sizes, and an analysis of the false positive and negative variant calling error rates. Improved genome phasing and family sequencing will increase the power of genome-wide sequencing as a clinical diagnosis tool and has myriad basic science applications.

In the mammalian genome, early- and late-replicating domains are often separated by temporal transition regions (TTRs) with novel properties and unknown functions. We identified a TTR in the mouse immunoglobulin heavy chain (Igh) locus, which contains replication origins that are silent in embryonic stem cells but activated during B cell development. To investigate which factors contribute to origin activation during B cell development, we systematically modified the genetic and epigenetic status of the endogenous Igh TTR and used a single-molecule approach to analyze DNA replication. Introduction of a transcription unit into the Igh TTR, activation of gene transcription, and enhancement of local histone modifications characteristic of active chromatin did not lead to origin activation. Moreover, very few replication initiation events were observed when two ectopic replication origin sequences were inserted into the TTR. These findings indicate that the Igh TTR represents a repressive compartment that inhibits replication initiation, thus maintaining the boundaries between early and late replication domains.

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.

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.

We have previously shown that human embryonic stem cells can be differentiated into embryonic and fetal type of red blood cells that sequentially express three types of hemoglobins recapitulating early human erythropoiesis. We report here that we have produced iPS from three somatic cell types: adult skin fibroblasts as well as embryonic and fetal mesenchymal stem cells. We show that regardless of the age of the donor cells, the iPS produced are fully reprogrammed into a pluripotent state that is undistinguishable from that of hESCs by low and high-throughput expression and detailed analysis of globin expression patterns by HPLC. This suggests that reprogramming with the four original Yamanaka pluripotency factors leads to complete erasure of all functionally important epigenetic marks associated with erythroid differentiation regardless of the age or the tissue type of the donor cells, at least as detected in these assays. The ability to produce large number of erythroid cells with embryonic and fetal-like characteristics is likely to have many translational applications.

Efficiency of reprogramming of human cells into induced pluripotent stem cells (iPSCs) has remained low. We report that individual adult human CD49f+ long-term hematopoietic stem cells (LT-HSCs) can be reprogrammed into iPSCs at close to 50% efficiency using Sendai virus transduction. This exquisite sensitivity to reprogramming is specific to LT-HSCs, since it progressively decreases in committed progenitors. LT-HSC reprogramming can follow multiple paths and is most efficient when transduction is performed after the cells have exited G0. Sequencing of 75 paired skin fibroblasts/LT-HSC samples collected from nine individuals revealed that LT-HSCs contain a lower load of somatic single-nucleotide variants (SNVs) and indels than skin fibroblasts and accumulate about 12 SNVs/year. Mutation analysis revealed that LT-HSCs and fibroblasts have very different somatic mutation signatures and that somatic mutations in iPSCs generally exist prior to reprogramming. LT-HSCs may become the preferred cell source for the production of clinical-grade iPSCs.

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.

Cytosines at cytosine-guanine (CG) dinucleotides are the near-exclusive target of DNA methyltransferases in mammalian genomes. Spontaneous deamination of methylcytosine to thymine makes methylated cytosines unusually susceptible to mutation and consequent depletion. The loci where CG dinucleotides remain relatively enriched, presumably due to their unmethylated status during the germ cell cycle, have been referred to as CpG islands. Currently, CpG islands are solely defined by base compositional criteria, allowing annotation of any sequenced genome. Using a novel bioinformatic approach, we show that CG clusters can be identified as an inherent property of genomic sequence without imposing a base compositional a priori assumption. We also show that the CG clusters co-localize in the human genome with hypomethylated loci and annotated transcription start sites to a greater extent than annotations produced by prior CpG island definitions. Moreover, this new approach allows CG clusters to be identified in a species-specific manner, revealing a degree of orthologous conservation that is not revealed by current base compositional approaches. Finally, our approach is able to identify methylating genomes (such as Takifugu rubripes) that lack CG clustering entirely, in which it is inappropriate to annotate CpG islands or CG clusters.