Dynamic Transformations of Genome-wide Epigenetic Marking and Transcriptional Control Establish T Cell Identity

T cell development comprises a stepwise process of commitment from a multipotent precursor. To define molecular mechanisms controlling this progression, we probed five stages spanning the commitment process using RNA-seq and ChIP-seq to track genome-wide shifts in transcription, cohorts of active transcription factor genes, histone modifications at diverse classes of cis-regulatory elements, and binding repertoire of GATA-3 and PU.1, transcription factors with complementary roles in T cell development. The results highlight potential promoter-distal cis-regulatory elements in play and reveal both activation sites and diverse mechanisms of repression that silence genes used in alternative lineages. Histone marking is dynamic and reversible, and though permissive marks anticipate, repressive marks often lag behind changes in transcription. In vivo binding of PU.1 and GATA-3 relative to epigenetic marking reveals distinctive factor-specific rules for recruitment of these crucial transcription factors to different subsets of their potential sites, dependent on dose and developmental context

Hematopoiesis and T-cell specification as a model developmental system

The pathway to generate T cells from hematopoietic stem cells guides progenitors through a succession of fate choices while balancing differentiation progression against proliferation, stage to stage. Many elements of the regulatory system that controls this process are known, but the requirement for multiple, functionally distinct transcription factors needs clarification in terms of gene network architecture. Here, we compare the features of the T-cell specification system with the rule sets underlying two other influential types of gene network models: first, the combinatorial, hierarchical regulatory systems that generate the orderly, synchronized increases in complexity in most invertebrate embryos; second, the dueling ‘master regulator’ systems that are commonly used to explain bistability in microbial systems and in many fate choices in terminal differentiation. The T-cell specification process shares certain features with each of these prevalent models but differs from both of them in central respects. The T-cell system is highly combinatorial but also highly dose-sensitive in its use of crucial regulatory factors. The roles of these factors are not always T-lineage-specific, but they balance and modulate each other's activities long before any mutually exclusive silencing occurs. T-cell specification may provide a new hybrid model for gene networks in vertebrate developmental systems

Loss of T Cell Progenitor Checkpoint Control Underlies Leukemia Initiation in Rag1-Deficient Nonobese Diabetic Mice

NOD mice exhibit major defects in the earliest stages of T cell development in the thymus. Genome-wide genetic and transcriptome analyses were used to investigate the origins and consequences of an early T cell developmental checkpoint breakthrough in Rag1-deficient NOD mice. Quantitative trait locus analysis mapped the presence of checkpoint breakthrough cells to several known NOD diabetes susceptibility regions, particularly insulin-dependent diabetes susceptibility genes (Idd)9/11 on chromosome 4, suggesting common genetic origins for T cell defects affecting this trait and autoimmunity. Genome-wide RNA deep-sequencing of NOD and B6 Rag1-deficient thymocytes revealed the effects of genetic background prior to breakthrough, as well as the cellular consequences of the breakthrough. Transcriptome comparison between the two strains showed enrichment in differentially expressed signal transduction genes, prominently tyrosine kinase and actin-binding genes, in accord with their divergent sensitivities to activating signals. Emerging NOD breakthrough cells aberrantly expressed both stem cell–associated proto-oncogenes, such as Lmo2, Hhex, Lyl1, and Kit, which are normally repressed at the commitment checkpoint, and post–β-selection checkpoint genes, including Cd2 and Cd5. Coexpression of genes characteristic of multipotent progenitors and more mature T cells persists in the expanding population of thymocytes and in the thymic leukemias that emerge with age in these mice. These results show that Rag1-deficient NOD thymocytes have T cell defects that can collapse regulatory boundaries at two early T cell checkpoints, which may predispose them to both leukemia and autoimmunity

Bcl11b and combinatorial resolution of cell fate in the T-cell gene regulatory network

T-cell development from hematopoietic progenitors depends on multiple transcription factors, mobilized and modulated by intrathymic Notch signaling. Key aspects of T-cell specification network architecture have been illuminated through recent reports defining roles of transcription factors PU.1, GATA-3, and E2A, their interactions with Notch signaling, and roles of Runx1, TCF-1, and Hes1, providing bases for a comprehensively updated model of the T-cell specification gene regulatory network presented herein. However, the role of lineage commitment factor Bcl11b has been unclear. We use self-organizing maps on 63 RNA-seq datasets from normal and perturbed T-cell development to identify functional targets of Bcl11b during commitment and relate them to other regulomes. We show that both activation and repression target genes can be bound by Bcl11b in vivo, and that Bcl11b effects overlap with E2A-dependent effects. The newly clarified role of Bcl11b distinguishes discrete components of commitment, resolving how innate lymphoid, myeloid, and dendritic, and B-cell fate alternatives are excluded by different mechanisms

GATA-3 Dose-Dependent Checkpoints in Early T Cell Commitment

GATA-3 expression is crucial for T cell development and peaks during commitment to the T cell lineage, midway through the CD4CD8 (double-negative [DN]) stages 1–3. We used RNA interference and conditional deletion to reduce GATA-3 protein acutely at specific points during T cell differentiation in vitro. Even moderate GATA-3 reduction killed DN1 cells, delayed progression to the DN2 stage, skewed DN2 gene regulation, and blocked appearance of the DN3 phenotype. Although a Bcl-2 transgene rescued DN1 survival and improved DN2 cell generation, it did not restore DN3 differentiation. Gene expression analyses (quantitative PCR, RNA sequencing) showed that GATA-3–deficient DN2 cells quickly upregulated genes, including Spi1 (PU.1) and Bcl11a, and downregulated genes, including Cpa3, Ets1, Zfpm1, Bcl11b, Il9r, and Il17rb with gene-specific kinetics and dose dependencies. These targets could mediate two distinct roles played by GATA-3 in lineage commitment, as revealed by removing wild-type or GATA-3–deficient early T lineage cells from environmental Notch signals. GATA-3 worked as a potent repressor of B cell potential even at low expression levels, so that only full deletion of GATA-3 enabled pro–T cells to reveal B cell potential. The ability of GATA-3 to block B cell development did not require T lineage commitment factor Bcl11b. In prethymic multipotent precursors, however, titration of GATA-3 activity using tamoxifen-inducible GATA-3 showed that GATA-3 inhibits B and myeloid developmental alternatives at different threshold doses. Furthermore, differential impacts of a GATA-3 obligate repressor construct imply that B and myeloid development are inhibited through distinct transcriptional mechanisms. Thus, the pattern of GATA-3 expression sequentially produces B lineage exclusion, T lineage progression, and myeloid-lineage exclusion for commitment

A far downstream enhancer for murine Bcl11b controls its T-cell specific expression

Bcl11b is a T-cell specific gene in hematopoiesis that begins expression during T-lineage commitment and is required for this process. Aberrant expression of BCL11B or proto-oncogene translocation to the vicinity of BCL11B can be a contributing factor in human T-ALL. To identify the mechanism that controls its distinctive T-lineage expression, we corrected the identified Bcl11b transcription start site and mapped a cell-type–specific differentially methylated region bracketing the Bcl11b promoter. We identified a 1.9-kb region 850 kb downstream of Bcl11b, “Major Peak,” distinguished by its dynamic histone marking pattern in development that mirrors the pattern at the Bcl11b promoter. Looping interactions between promoter-proximal elements including the differentially methylated region and downstream elements in the Major Peak are required to recapitulate the T-cell specific expression of Bcl11b in stable reporter assays. Functional dissection of the Major Peak sequence showed distinct subregions, in which TCF-1 sites and a conserved element were required for T-lineage–specific activation and silencing in non-T cells. A bacterial artificial chromosome encompassing the full Bcl11b gene still required the addition of the Major Peak to exhibit T-cell specific expression. Thus, promoter-proximal and Major Peak sequences are cis-regulatory elements that interact over 850 kb to control expression of Bcl11b in hematopoietic cells

Asynchronous combinatorial action of four regulatory factors activates Bcl11b for T cell commitment

During T cell development, multipotent progenitors relinquish competence for other fates and commit to the T cell lineage by turning on Bcl11b, which encodes a transcription factor. To clarify lineage commitment mechanisms, we followed developing T cells at the single-cell level using Bcl11b knock-in fluorescent reporter mice. Notch signaling and Notch-activated transcription factors collaborate to activate Bcl11b expression irrespectively of Notch-dependent proliferation. These inputs work via three distinct, asynchronous mechanisms: an early locus 'poising' function dependent on TCF-1 and GATA-3, a stochastic-permissivity function dependent on Notch signaling, and a separate amplitude-control function dependent on Runx1, a factor already present in multipotent progenitors. Despite their necessity for Bcl11b expression, these inputs act in a stage-specific manner, providing a multitiered mechanism for developmental gene regulation

Asynchronous combinatorial action of four regulatory factors activates Bcl11b for T cell commitment

During T cell development, multipotent progenitors relinquish competence for other fates and commit to the T cell lineage by turning on Bcl11b, which encodes a transcription factor. To clarify lineage commitment mechanisms, we followed developing T cells at the single-cell level using Bcl11b knock-in fluorescent reporter mice. Notch signaling and Notch-activated transcription factors collaborate to activate Bcl11b expression irrespectively of Notch-dependent proliferation. These inputs work via three distinct, asynchronous mechanisms: an early locus 'poising' function dependent on TCF-1 and GATA-3, a stochastic-permissivity function dependent on Notch signaling, and a separate amplitude-control function dependent on Runx1, a factor already present in multipotent progenitors. Despite their necessity for Bcl11b expression, these inputs act in a stage-specific manner, providing a multitiered mechanism for developmental gene regulation

Gene Expression Profile of Neuronal Progenitor Cells Derived from hESCs: Activation of Chromosome 11p15.5 and Comparison to Human Dopaminergic Neurons

BACKGROUND: We initiated differentiation of human embryonic stem cells (hESCs) into dopamine neurons, obtained a purified population of neuronal precursor cells by cell sorting, and determined patterns of gene transcription. METHODOLOGY: Dopaminergic differentiation of hESCs was initiated by culturing hESCs with a feeder layer of PA6 cells. Differentiating cells were then sorted to obtain a pure population of PSA-NCAM-expressing neuronal precursors, which were then analyzed for gene expression using Massive Parallel Signature Sequencing (MPSS). Individual genes as well as regions of the genome which were activated were determined. PRINCIPAL FINDINGS: A number of genes known to be involved in the specification of dopaminergic neurons, including MSX1, CDKN1C, Pitx1 and Pitx2, as well as several novel genes not previously associated with dopaminergic differentiation, were expressed. Notably, we found that a specific region of the genome located on chromosome 11p15.5 was highly activated. This region contains several genes which have previously been associated with the function of dopaminergic neurons, including the gene for tyrosine hydroxylase (TH), the rate-limiting enzyme in catecholamine biosynthesis, IGF2, and CDKN1C, which cooperates with Nurr1 in directing the differentiation of dopaminergic neurons. Other genes in this region not previously recognized as being involved in the functions of dopaminergic neurons were also activated, including H19, TSSC4, and HBG2. IGF2 and CDKN1C were also found to be highly expressed in mature human TH-positive dopamine neurons isolated from human brain samples by laser capture. CONCLUSIONS: The present data suggest that the H19-IGF2 imprinting region on chromosome 11p15.5 is involved in the process through which undifferentiated cells are specified to become neuronal precursors and/or dopaminergic neurons

Finishing the euchromatic sequence of the human genome

The sequence of the human genome encodes the genetic instructions for human physiology, as well as rich information about human evolution. In 2001, the International Human Genome Sequencing Consortium reported a draft sequence of the euchromatic portion of the human genome. Since then, the international collaboration has worked to convert this draft into a genome sequence with high accuracy and nearly complete coverage. Here, we report the result of this finishing process. The current genome sequence (Build 35) contains 2.85 billion nucleotides interrupted by only 341 gaps. It covers ∼99% of the euchromatic genome and is accurate to an error rate of ∼1 event per 100,000 bases. Many of the remaining euchromatic gaps are associated with segmental duplications and will require focused work with new methods. The near-complete sequence, the first for a vertebrate, greatly improves the precision of biological analyses of the human genome including studies of gene number, birth and death. Notably, the human enome seems to encode only 20,000-25,000 protein-coding genes. The genome sequence reported here should serve as a firm foundation for biomedical research in the decades ahead