Detection and analysis of complex structural variation in human genomes across populations and in brains of donors with psychiatric disorders.

Complex structural variations (cxSVs) are often overlooked in genome analyses due to detection challenges. We developed ARC-SV, a probabilistic and machine-learning-based method that enables accurate detection and reconstruction of cxSVs from standard datasets. By applying ARC-SV across 4,262 genomes representing all continental populations, we identified cxSVs as a significant source of natural human genetic variation. Rare cxSVs have a propensity to occur in neural genes and loci that underwent rapid human-specific evolution, including those regulating corticogenesis. By performing single-nucleus multiomics in postmortem brains, we discovered cxSVs associated with differential gene expression and chromatin accessibility across various brain regions and cell types. Additionally, cxSVs detected in brains of psychiatric cases are enriched for linkage with psychiatric GWAS risk alleles detected in the same brains. Furthermore, our analysis revealed significantly decreased brain-region- and cell-type-specific expression of cxSV genes, specifically for psychiatric cases, implicating cxSVs in the molecular etiology of major neuropsychiatric disorders.

Chromatin Accessibility Plays an Important Epigenetic Role on Antibody Expression From CMV Promoter and DNA Elements Flanking the CHO TI Host Landing-Pad.

Targeted integration (TI) Chinese hamster ovary (CHO) platforms are commonly used for protein expression. However, the impact of epigenetic modifications on protein expression in TI cell lines remains elusive since almost all the epigenetic studies focus on random integration (RI) of the gene of interest and only within the promoter region. To address the impact of epigenetic modifications on TI CHO cells, we utilized a standard mAb-1 to identify and characterize TI clones with the same transgene copy numbers but different levels of transgene transcription and titer. Surprisingly, while CMV promoters were not methylated and histone acetylation/methylation was present, these epigenetic markers did not trend with mRNA transcription and protein expression in our TI model. Instead, ATAC-seq data analysis revealed that differences in chromatin accessibility within the TI site could be a major factor impacting these observed differences. However, neither chromatin accessibility nor histone acetylation/methylation profiles in early cultures were predictive of high-expressing clones early during the CLD process. Finally, modulation of the histone profiles (H3K27ac and H3K4me3) at the CMV promoters within the TI integration site using dCas9 fusion proteins was not effective in further increasing mAb titers which could have been likely due to interference of the dCas9 fusion proteins with transcription from the CMV promoters. Overall, our data suggests increasing chromatin accessibility at the TI site is the most effective way to increase mRNA transcription and hence, productivity in TI cell lines.

Overloading And unpacKing (OAK) - droplet-based combinatorial indexing for ultra-high throughput single-cell multiomic profiling.

Multiomic profiling of single cells by sequencing is a powerful technique for investigating cellular diversity. Existing droplet-based microfluidic methods produce many cell-free droplets, underutilizing bead barcodes and reagents. Combinatorial indexing on microplates is more efficient for barcoding but labor-intensive. Here we present Overloading And unpacKing (OAK), which uses a droplet-based barcoding system for initial compartmentalization followed by a second aliquoting round to achieve combinatorial indexing. We demonstrate OAK's versatility with single-cell RNA sequencing as well as paired single-nucleus RNA sequencing and accessible chromatin profiling. We further showcase OAK's performance on complex samples, including differentiated bronchial epithelial cells and primary retinal tissue. Finally, we examine transcriptomic responses of over 400,000 melanoma cells to a RAF inhibitor, belvarafenib, discovering a rare resistant cell population (0.12%). OAK's ultra-high throughput, broad compatibility, high sensitivity, and simplified procedures make it a powerful tool for large-scale molecular analysis, even for rare cells.

Single-cell long-read targeted sequencing reveals transcriptional variation in ovarian cancer.

Single-cell RNA sequencing predominantly employs short-read sequencing to characterize cell types, states and dynamics; however, it is inadequate for comprehensive characterization of RNA isoforms. Long-read sequencing technologies enable single-cell RNA isoform detection but are hampered by lower throughput and unintended sequencing of artifacts. Here we develop Single-cell Targeted Isoform Long-Read Sequencing (scTaILoR-seq), a hybridization capture method which targets over a thousand genes of interest, improving the median number of on-target transcripts per cell by 29-fold. We use scTaILoR-seq to identify and quantify RNA isoforms from ovarian cancer cell lines and primary tumors, yielding 10,796 single-cell transcriptomes. Using long-read variant calling we reveal associations of expressed single nucleotide variants (SNVs) with alternative transcript structures. Phasing of SNVs across transcripts enables the measurement of allelic imbalance within distinct cell populations. Overall, scTaILoR-seq is a long-read targeted RNA sequencing method and analytical framework for exploring transcriptional variation at single-cell resolution.

The ENCODE mouse postnatal developmental time course identifies regulatory programs of cell types and cell states.

Postnatal genomic regulation significantly influences tissue and organ maturation but is under-studied relative to existing genomic catalogs of adult tissues or prenatal development in mouse. The ENCODE4 consortium generated the first comprehensive single-nucleus resource of postnatal regulatory events across a diverse set of mouse tissues. The collection spans seven postnatal time points, mirroring human development from childhood to adulthood, and encompasses five core tissues. We identified 30 cell types, further subdivided into 69 subtypes and cell states across adrenal gland, left cerebral cortex, hippocampus, heart, and gastrocnemius muscle. Our annotations cover both known and novel cell differentiation dynamics ranging from early hippocampal neurogenesis to a new sex-specific adrenal gland population during puberty. We used an ensemble Latent Dirichlet Allocation strategy with a curated vocabulary of 2,701 regulatory genes to identify regulatory "topics," each of which is a gene vector, linked to cell type differentiation, subtype specialization, and transitions between cell states. We find recurrent regulatory topics in tissue-resident macrophages, neural cell types, endothelial cells across multiple tissues, and cycling cells of the adrenal gland and heart. Cell-type-specific topics are enriched in transcription factors and microRNA host genes, while chromatin regulators dominate mitosis topics. Corresponding chromatin accessibility data reveal dynamic and sex-specific regulatory elements, with enriched motifs matching transcription factors in regulatory topics. Together, these analyses identify both tissue-specific and common regulatory programs in postnatal development across multiple tissues through the lens of the factors regulating transcription.

The ENCODE Imputation Challenge: a critical assessment of methods for cross-cell type imputation of epigenomic profiles.

A promising alternative to comprehensively performing genomics experiments is to, instead, perform a subset of experiments and use computational methods to impute the remainder. However, identifying the best imputation methods and what measures meaningfully evaluate performance are open questions. We address these questions by comprehensively analyzing 23 methods from the ENCODE Imputation Challenge. We find that imputation evaluations are challenging and confounded by distributional shifts from differences in data collection and processing over time, the amount of available data, and redundancy among performance measures. Our analyses suggest simple steps for overcoming these issues and promising directions for more robust research.

Proinflammatory polarization of monocytes by particulate air pollutants is mediated by induction of trained immunity in pediatric asthma.

BACKGROUND: The impact of exposure to air pollutants, such as fine particulate matter (PM), on the immune system and its consequences on pediatric asthma, are not well understood. We investigated whether ambient levels of fine PM with aerodynamic diameter ≤2.5 microns (PM2.5 ) are associated with alterations in circulating monocytes in children with or without asthma. METHODS: Monocyte phenotyping was performed by cytometry time-of-flight (CyTOF). Cytokines were measured using cytometric bead array and Luminex assay. ChIP-Seq was utilized to address histone modifications in monocytes. RESULTS: Increased exposure to ambient PM2.5 was linked to specific monocyte subtypes, particularly in children with asthma. Mechanistically, we hypothesized that innate trained immunity is evoked by a primary exposure to fine PM and accounts for an enhanced inflammatory response after secondary stimulation in vitro. We determined that the trained immunity was induced in circulating monocytes by fine particulate pollutants, and it was characterized by the upregulation of proinflammatory mediators, such as TNF, IL-6, and IL-8, upon stimulation with house dust mite or lipopolysaccharide. This phenotype was epigenetically controlled by enhanced H3K27ac marks in circulating monocytes. CONCLUSION: The specific alterations of monocytes after ambient pollution exposure suggest a possible prognostic immune signature for pediatric asthma, and pollution-induced trained immunity may provide a potential therapeutic target for asthmatic children living in areas with increased air pollution.

KLF4 recruits SWI/SNF to increase chromatin accessibility and reprogram the endothelial enhancer landscape under laminar shear stress.

Physiologic laminar shear stress (LSS) induces an endothelial gene expression profile that is vasculo-protective. In this report, we delineate how LSS mediates changes in the epigenetic landscape to promote this beneficial response. We show that under LSS, KLF4 interacts with the SWI/SNF nucleosome remodeling complex to increase accessibility at enhancer sites that promote the expression of homeostatic endothelial genes. By combining molecular and computational approaches we discover enhancers that loop to promoters of KLF4- and LSS-responsive genes that stabilize endothelial cells and suppress inflammation, such as BMPR2, SMAD5, and DUSP5. By linking enhancers to genes that they regulate under physiologic LSS, our work establishes a foundation for interpreting how non-coding DNA variants in these regions might disrupt protective gene expression to influence vascular disease.

Multiomic analysis reveals cell-type-specific molecular determinants of COVID-19 severity.

The determinants of severe COVID-19 in healthy adults are poorly understood, which limits the opportunity for early intervention. We present a multiomic analysis using machine learning to characterize the genomic basis of COVID-19 severity. We use single-cell multiome profiling of human lungs to link genetic signals to cell-type-specific functions. We discover >1,000 risk genes across 19 cell types, which account for 77% of the SNP-based heritability for severe disease. Genetic risk is particularly focused within natural killer (NK) cells and T cells, placing the dysfunction of these cells upstream of severe disease. Mendelian randomization and single-cell profiling of human NK cells support the role of NK cells and further localize genetic risk to CD56bright NK cells, which are key cytokine producers during the innate immune response. Rare variant analysis confirms the enrichment of severe-disease-associated genetic variation within NK-cell risk genes. Our study provides insights into the pathogenesis of severe COVID-19 with potential therapeutic targets.

Genome-wide identification of the genetic basis of amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a complex disease that leads to motor neuron death. Despite heritability estimates of 52%, genome-wide association studies (GWASs) have discovered relatively few loci. We developed a machine learning approach called RefMap, which integrates functional genomics with GWAS summary statistics for gene discovery. With transcriptomic and epigenetic profiling of motor neurons derived from induced pluripotent stem cells (iPSCs), RefMap identified 690 ALS-associated genes that represent a 5-fold increase in recovered heritability. Extensive conservation, transcriptome, network, and rare variant analyses demonstrated the functional significance of candidate genes in healthy and diseased motor neurons and brain tissues. Genetic convergence between common and rare variation highlighted KANK1 as a new ALS gene. Reproducing KANK1 patient mutations in human neurons led to neurotoxicity and demonstrated that TDP-43 mislocalization, a hallmark pathology of ALS, is downstream of axonal dysfunction. RefMap can be readily applied to other complex diseases.

The dynamic, combinatorial cis-regulatory lexicon of epidermal differentiation.

Transcription factors bind DNA sequence motif vocabularies in cis-regulatory elements (CREs) to modulate chromatin state and gene expression during cell state transitions. A quantitative understanding of how motif lexicons influence dynamic regulatory activity has been elusive due to the combinatorial nature of the cis-regulatory code. To address this, we undertook multiomic data profiling of chromatin and expression dynamics across epidermal differentiation to identify 40,103 dynamic CREs associated with 3,609 dynamically expressed genes, then applied an interpretable deep-learning framework to model the cis-regulatory logic of chromatin accessibility. This analysis framework identified cooperative DNA sequence rules in dynamic CREs regulating synchronous gene modules with diverse roles in skin differentiation. Massively parallel reporter assay analysis validated temporal dynamics and cooperative cis-regulatory logic. Variants linked to human polygenic skin disease were enriched in these time-dependent combinatorial motif rules. This integrative approach shows the combinatorial cis-regulatory lexicon of epidermal differentiation and represents a general framework for deciphering the organizational principles of the cis-regulatory code of dynamic gene regulation.

ALDH1A3 Coordinates Metabolism With Gene Regulation in Pulmonary Arterial Hypertension.

BACKGROUND: Metabolic alterations provide substrates that influence chromatin structure to regulate gene expression that determines cell function in health and disease. Heightened proliferation of smooth muscle cells (SMC) leading to the formation of a neointima is a feature of pulmonary arterial hypertension (PAH) and systemic vascular disease. Increased glycolysis is linked to the proliferative phenotype of these SMC. METHODS: RNA sequencing was applied to pulmonary arterial SMC (PASMC) from PAH patients with and without a BMPR2 (bone morphogenetic receptor 2) mutation versus control PASMC to uncover genes required for their heightened proliferation and glycolytic metabolism. Assessment of differentially expressed genes established metabolism as a major pathway, and the most highly upregulated metabolic gene in PAH PASMC was aldehyde dehydrogenase family 1 member 3 (ALDH1A3), an enzyme previously linked to glycolysis and proliferation in cancer cells and systemic vascular SMC. We determined if these functions are ALDH1A3-dependent in PAH PASMC, and if ALDH1A3 is required for the development of pulmonary hypertension in a transgenic mouse. Nuclear localization of ALDH1A3 in PAH PASMC led us to determine whether and how this enzyme coordinately regulates gene expression and metabolism in PAH PASMC. RESULTS: ALDH1A3 mRNA and protein were increased in PAH versus control PASMC, and ALDH1A3 was required for their highly proliferative and glycolytic properties. Mice with Aldh1a3 deleted in SMC did not develop hypoxia-induced pulmonary arterial muscularization or pulmonary hypertension. Nuclear ALDH1A3 converted acetaldehyde to acetate to produce acetyl coenzyme A to acetylate H3K27, marking active enhancers. This allowed for chromatin modification at NFYA (nuclear transcription factor Y subunit α) binding sites via the acetyltransferase KAT2B (lysine acetyltransferase 2B) and permitted NFY-mediated transcription of cell cycle and metabolic genes that is required for ALDH1A3-dependent proliferation and glycolysis. Loss of BMPR2 in PAH SMC with or without a mutation upregulated ALDH1A3, and transcription of NFYA and ALDH1A3 in PAH PASMC was ß-catenin dependent. CONCLUSIONS: Our studies have uncovered a metabolic-transcriptional axis explaining how dividing cells use ALDH1A3 to coordinate their energy needs with the epigenetic and transcriptional regulation of genes required for SMC proliferation. They suggest that selectively disrupting the pivotal role of ALDH1A3 in PAH SMC, but not endothelial cells, is an important therapeutic consideration.

PPARγ-p53-Mediated Vasculoregenerative Program to Reverse Pulmonary Hypertension.

RATIONALE: In pulmonary arterial hypertension (PAH), endothelial dysfunction and obliterative vascular disease are associated with DNA damage and impaired signaling of BMPR2 (bone morphogenetic protein type 2 receptor) via two downstream transcription factors, PPARγ (peroxisome proliferator-activated receptor gamma), and p53. OBJECTIVE: We investigated the vasculoprotective and regenerative potential of a newly identified PPARγ-p53 transcription factor complex in the pulmonary endothelium. METHODS AND RESULTS: In this study, we identified a pharmacologically inducible vasculoprotective mechanism in pulmonary arterial and lung MV (microvascular) endothelial cells in response to DNA damage and oxidant stress regulated in part by a BMPR2 dependent transcription factor complex between PPARγ and p53. Chromatin immunoprecipitation sequencing and RNA-sequencing established an inducible PPARγ-p53 mediated regenerative program regulating 19 genes involved in lung endothelial cell survival, angiogenesis and DNA repair including, EPHA2 (ephrin type-A receptor 2), FHL2 (four and a half LIM domains protein 2), JAG1 (jagged 1), SULF2 (extracellular sulfatase Sulf-2), and TIGAR (TP53-inducible glycolysis and apoptosis regulator). Expression of these genes was partially impaired when the PPARγ-p53 complex was pharmacologically disrupted or when BMPR2 was reduced in pulmonary artery endothelial cells (PAECs) subjected to oxidative stress. In endothelial cell-specific Bmpr2-knockout mice unable to stabilize p53 in endothelial cells under oxidative stress, Nutlin-3 rescued endothelial p53 and PPARγ-p53 complex formation and induced target genes, such as APLN (apelin) and JAG1, to regenerate pulmonary microvessels and reverse pulmonary hypertension. In PAECs from BMPR2 mutant PAH patients, pharmacological induction of p53 and PPARγ-p53 genes repaired damaged DNA utilizing genes from the nucleotide excision repair pathway without provoking PAEC apoptosis. CONCLUSIONS: We identified a novel therapeutic strategy that activates a vasculoprotective gene regulation program in PAECs downstream of dysfunctional BMPR2 to rehabilitate PAH PAECs, regenerate pulmonary microvessels, and reverse disease. Our studies pave the way for p53-based vasculoregenerative therapies for PAH by extending the therapeutic focus to PAEC dysfunction and to DNA damage associated with PAH progression.

Smooth Muscle Contact Drives Endothelial Regeneration by BMPR2-Notch1-Mediated Metabolic and Epigenetic Changes.

RATIONALE: Maintaining endothelial cells (EC) as a monolayer in the vessel wall depends on their metabolic state and gene expression profile, features influenced by contact with neighboring cells such as pericytes and smooth muscle cells (SMC). Failure to regenerate a normal EC monolayer in response to injury can result in occlusive neointima formation in diseases such as atherosclerosis and pulmonary arterial hypertension. OBJECTIVE: We investigated the nature and functional importance of contact-dependent communication between SMC and EC to maintain EC integrity. METHODS AND RESULTS: We found that in SMC and EC contact cocultures, BMPR2 (bone morphogenetic protein receptor 2) is required by both cell types to produce collagen IV to activate ILK (integrin-linked kinase). This enzyme directs p-JNK (phospho-c-Jun N-terminal kinase) to the EC membrane, where it stabilizes presenilin1 and releases N1ICD (Notch1 intracellular domain) to promote EC proliferation. This response is necessary for EC regeneration after carotid artery injury. It is deficient in EC-SMC Bmpr2 double heterozygous mice in association with reduced collagen IV production, decreased N1ICD, and attenuated EC proliferation, but can be rescued by targeting N1ICD to EC. Deletion of EC- Notch1 in transgenic mice worsens hypoxia-induced pulmonary hypertension, in association with impaired EC regenerative function associated with loss of precapillary arteries. We further determined that N1ICD maintains EC proliferative capacity by increasing mitochondrial mass and by inducing the phosphofructokinase PFKFB3 (fructose-2,6-bisphosphatase 3). Chromatin immunoprecipitation sequencing analyses showed that PFKFB3 is required for citrate-dependent H3K27 acetylation at enhancer sites of genes regulated by the acetyl transferase p300 and by N1ICD or the N1ICD target MYC and necessary for EC proliferation and homeostasis. CONCLUSIONS: Thus, SMC-EC contact is required for activation of Notch1 by BMPR2, to coordinate metabolism with chromatin remodeling of genes that enable EC regeneration, and to maintain monolayer integrity and vascular homeostasis in response to injury.

NF90/ILF3 is a transcription factor that promotes proliferation over differentiation by hierarchical regulation in K562 erythroleukemia cells.

NF90 and splice variant NF110 are DNA- and RNA-binding proteins encoded by the Interleukin enhancer-binding factor 3 (ILF3) gene that have been established to regulate RNA splicing, stabilization and export. The roles of NF90 and NF110 in regulating transcription as chromatin-interacting proteins have not been comprehensively characterized. Here, chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) identified 9,081 genomic sites specifically occupied by NF90/NF110 in K562 cells. One third of NF90/NF110 peaks occurred at promoters of annotated genes. NF90/NF110 occupancy colocalized with chromatin marks associated with active promoters and strong enhancers. Comparison with 150 ENCODE ChIP-seq experiments revealed that NF90/NF110 clustered with transcription factors exhibiting preference for promoters over enhancers (POLR2A, MYC, YY1). Differential gene expression analysis following shRNA knockdown of NF90/NF110 in K562 cells revealed that NF90/NF110 activates transcription factors that drive growth and proliferation (EGR1, MYC), while attenuating differentiation along the erythroid lineage (KLF1). NF90/NF110 associates with chromatin to hierarchically regulate transcription factors that promote proliferation and suppress differentiation.

Lineage-specific dynamic and pre-established enhancer-promoter contacts cooperate in terminal differentiation.

Chromosome conformation is an important feature of metazoan gene regulation; however, enhancer-promoter contact remodeling during cellular differentiation remains poorly understood. To address this, genome-wide promoter capture Hi-C (CHi-C) was performed during epidermal differentiation. Two classes of enhancer-promoter contacts associated with differentiation-induced genes were identified. The first class ('gained') increased in contact strength during differentiation in concert with enhancer acquisition of the H3K27ac activation mark. The second class ('stable') were pre-established in undifferentiated cells, with enhancers constitutively marked by H3K27ac. The stable class was associated with the canonical conformation regulator cohesin, whereas the gained class was not, implying distinct mechanisms of contact formation and regulation. Analysis of stable enhancers identified a new, essential role for a constitutively expressed, lineage-restricted ETS-family transcription factor, EHF, in epidermal differentiation. Furthermore, neither class of contacts was observed in pluripotent cells, suggesting that lineage-specific chromatin structure is established in tissue progenitor cells and is further remodeled in terminal differentiation.

Identification of Human Neuronal Protein Complexes Reveals Biochemical Activities and Convergent Mechanisms of Action in Autism Spectrum Disorders.

The prevalence of autism spectrum disorders (ASDs) is rapidly growing, yet its molecular basis is poorly understood. We used a systems approach in which ASD candidate genes were mapped onto the ubiquitous human protein complexes and the resulting complexes were characterized. The studies revealed the role of histone deacetylases (HDAC1/2) in regulating the expression of ASD orthologs in the embryonic mouse brain. Proteome-wide screens for the co-complexed subunits with HDAC1 and six other key ASD proteins in neuronal cells revealed a protein interaction network, which displayed preferential expression in fetal brain development, exhibited increased deleterious mutations in ASD cases, and were strongly regulated by FMRP and MECP2 causal for Fragile X and Rett syndromes, respectively. Overall, our study reveals molecular components in ASD, suggests a shared mechanism between the syndromic and idiopathic forms of ASDs, and provides a systems framework for analyzing complex human diseases.

Integrated systems analysis reveals a molecular network underlying autism spectrum disorders.

Autism is a complex disease whose etiology remains elusive. We integrated previously and newly generated data and developed a systems framework involving the interactome, gene expression and genome sequencing to identify a protein interaction module with members strongly enriched for autism candidate genes. Sequencing of 25 patients confirmed the involvement of this module in autism, which was subsequently validated using an independent cohort of over 500 patients. Expression of this module was dichotomized with a ubiquitously expressed subcomponent and another subcomponent preferentially expressed in the corpus callosum, which was significantly affected by our identified mutations in the network center. RNA-sequencing of the corpus callosum from patients with autism exhibited extensive gene mis-expression in this module, and our immunochemical analysis showed that the human corpus callosum is predominantly populated by oligodendrocyte cells. Analysis of functional genomic data further revealed a significant involvement of this module in the development of oligodendrocyte cells in mouse brain. Our analysis delineates a natural network involved in autism, helps uncover novel candidate genes for this disease and improves our understanding of its molecular pathology.

Architecture of the human regulatory network derived from ENCODE data.

Transcription factors bind in a combinatorial fashion to specify the on-and-off states of genes; the ensemble of these binding events forms a regulatory network, constituting the wiring diagram for a cell. To examine the principles of the human transcriptional regulatory network, we determined the genomic binding information of 119 transcription-related factors in over 450 distinct experiments. We found the combinatorial, co-association of transcription factors to be highly context specific: distinct combinations of factors bind at specific genomic locations. In particular, there are significant differences in the binding proximal and distal to genes. We organized all the transcription factor binding into a hierarchy and integrated it with other genomic information (for example, microRNA regulation), forming a dense meta-network. Factors at different levels have different properties; for instance, top-level transcription factors more strongly influence expression and middle-level ones co-regulate targets to mitigate information-flow bottlenecks. Moreover, these co-regulations give rise to many enriched network motifs (for example, noise-buffering feed-forward loops). Finally, more connected network components are under stronger selection and exhibit a greater degree of allele-specific activity (that is, differential binding to the two parental alleles). The regulatory information obtained in this study will be crucial for interpreting personal genome sequences and understanding basic principles of human biology and disease.

Tcf7 is an important regulator of the switch of self-renewal and differentiation in a multipotential hematopoietic cell line.

A critical problem in biology is understanding how cells choose between self-renewal and differentiation. To generate a comprehensive view of the mechanisms controlling early hematopoietic precursor self-renewal and differentiation, we used systems-based approaches and murine EML multipotential hematopoietic precursor cells as a primary model. EML cells give rise to a mixture of self-renewing Lin-SCA+CD34+ cells and partially differentiated non-renewing Lin-SCA-CD34- cells in a cell autonomous fashion. We identified and validated the HMG box protein TCF7 as a regulator in this self-renewal/differentiation switch that operates in the absence of autocrine Wnt signaling. We found that Tcf7 is the most down-regulated transcription factor when CD34+ cells switch into CD34- cells, using RNA-Seq. We subsequently identified the target genes bound by TCF7, using ChIP-Seq. We show that TCF7 and RUNX1 (AML1) bind to each other's promoter regions and that TCF7 is necessary for the production of the short isoforms, but not the long isoforms of RUNX1, suggesting that TCF7 and the short isoforms of RUNX1 function coordinately in regulation. Tcf7 knock-down experiments and Gene Set Enrichment Analyses suggest that TCF7 plays a dual role in promoting the expression of genes characteristic of self-renewing CD34+ cells while repressing genes activated in partially differentiated CD34- state. Finally a network of up-regulated transcription factors of CD34+ cells was constructed. Factors that control hematopoietic stem cell (HSC) establishment and development, cell growth, and multipotency were identified. These studies in EML cells demonstrate fundamental cell-intrinsic properties of the switch between self-renewal and differentiation, and yield valuable insights for manipulating HSCs and other differentiating systems.