Shortcut barcoding and early pooling for scalable multiplex single-cell reduced-representation CpG methylation sequencing at single nucleotide resolution

DNA methylation is essential for a wide variety of biological processes, yet the development of a highly efficient and robust technology remains a challenge for routine single-cell analysis. We developed a multiplex scalable single-cell reduced representation bisulfite sequencing (msRRBS) technology. It allows cell-specific barcoded DNA fragments of individual cells to be pooled before bisulfite conversion, free of enzymatic modification or physical capture of the DNA ends, and achieves read mapping rates of 62.5 ± 3.9%, covering 60.0 ± 1.4% of CpG islands and 71.6 ± 1.6% of promoters in K562 cells. Its reproducibility is shown in duplicates of bulk cells with close to perfect correlation (R = 0.97-0.99). At a low 1 Mb of clean reads, msRRBS provides highly consistent coverage of CpG islands and promoters, outperforming the conventional methods with orders of magnitude reduction in cost. Here, we use this method to characterize the distinct methylation patterns and cellular heterogeneity of six cell lines, plus leukemia and hepatocellular carcinoma models. Taking 4 h of hands-on time, msRRBS offers a unique, highly efficient approach for dissecting methylation heterogeneity in a variety of multicellular systems.

Sample-multiplexing approaches for single-cell sequencing

Single-cell sequencing is widely used in biological and medical studies. However, its application with multiple samples is hindered by inefficient sample processing, high experimental costs, ambiguous identification of true single cells, and technical batch effects. Here, we introduce sample-multiplexing approaches for single-cell sequencing in transcriptomics, epigenomics, genomics, and multiomics. In single-cell transcriptomics, sample multiplexing uses variants of native or artificial features as sample markers, enabling sample pooling and decoding. Such features include: (1) natural genetic variation, (2) nucleotide-barcode anchoring on cellular or nuclear membranes, (3) nucleotide-barcode internalization to the cytoplasm or nucleus, (4) vector-based barcode expression in cells, and (5) nucleotide-barcode incorporation during library construction. Other single-cell omics methods are based on similar concepts, particularly single-cell combinatorial indexing. These methods overcome current challenges, while enabling super-loading of single cells. Finally, selection guidelines are presented that can accelerate technological application.

Spatial transcriptome profiling by MERFISH reveals fetal liver hematopoietic stem cell niche architecture

The hematopoietic stem cell (HSC) niche has been extensively studied in bone marrow, yet a more systematic investigation into the microenvironment regulation of hematopoiesis in fetal liver is necessary. Here we investigate the spatial organization and transcriptional profile of individual cells in both wild type (WT) and Tet2^−/− fetal livers, by multiplexed error robust fluorescence in situ hybridization. We find that specific pairs of fetal liver cell types are preferentially positioned next to each other. Ligand-receptor signaling molecule pairs such as Kitl and Kit are enriched in neighboring cell types. The majority of HSCs are in direct contact with endothelial cells (ECs) in both WT and Tet2^−/− fetal livers. Loss of Tet2 increases the number of HSCs, and upregulates Wnt and Notch signaling genes in the HSC niche. Two subtypes of ECs, arterial ECs and sinusoidal ECs, and other cell types contribute distinct signaling molecules to the HSC niche. Collectively, this study provides a comprehensive picture and bioinformatic foundation for HSC spatial regulation in fetal liver.

Colorectal Cancer Stem Cell States Uncovered by Simultaneous Single-Cell Analysis of Transcriptome and Telomeres

Cancer stem cells (CSCs) presumably contribute to tumor progression and drug resistance, yet their definitive features have remained elusive. Here, simultaneous measurement of telomere length and transcriptome in the same cells enables systematic assessment of CSCs in primary colorectal cancer (CRC). The in-depth transcriptome profiled by SMART-seq2 is independently validated by high-throughput scRNA-seq using 10 × Genomics. It is found that rare CSCs exist in dormant state and display plasticity toward cancer epithelial cells (EPCs) that essentially are presumptive tumor-initiating cells (TICs), while both retaining the prominent signaling pathways including WNT, TGF-β, and HIPPO/YAP. Moreover, CSCs exhibit chromosome copy number variation (CNV) pattern resembling cancer EPCs but distinct from normal stem cells, suggesting the phylogenetic relationship between CSCs and cancer EPCs. Notably, CSCs maintain shorter telomeres and possess minimal telomerase activity consistent with their nonproliferative nature, unlike cancer EPCs. Additionally, the specific signature of CSCs particularly NOTUM, SMOC2, BAMBI, PHLDA1, and TNFRSF19 correlates with the prognosis of CRC. These findings characterize the heterogeneity of CSCs and their linkage to cancer EPCs/TICs, some of which are conventionally regarded as CSCs.

Dysregulation of BRD4 Function Underlies the Functional Abnormalities of MeCP2 Mutant Neurons

Rett syndrome (RTT), mainly caused by mutations in methyl-CpG binding protein 2 (MeCP2), is one of the most prevalent intellectual disorders without effective therapies. Here, we used 2D and 3D human brain cultures to investigate MeCP2 function. We found that MeCP2 mutations cause severe abnormalities in human interneurons (INs). Surprisingly, treatment with a BET inhibitor, JQ1, rescued the molecular and functional phenotypes of MeCP2 mutant INs. We uncovered that abnormal increases in chromatin binding of BRD4 and enhancer-promoter interactions underlie the abnormal transcription in MeCP2 mutant INs, which were recovered to normal levels by JQ1. We revealed cell-type-specific transcriptome impairment in MeCP2 mutant region-specific human brain organoids that were rescued by JQ1. Finally, JQ1 ameliorated RTT-like phenotypes in mice. These data demonstrate that BRD4 dysregulation is a critical driver for RTT etiology and suggest that targeting BRD4 could be a potential therapeutic opportunity for RTT.

A versatile flow-based assay for immunocyte-mediated cytotoxicity

Cell-mediated cytotoxicity is a critical function of the immune system in mounting defense against pathogens and cancers. Current methods that allow direct evaluation of cell-mediated cytotoxicity suffer from a wide-range of drawbacks. Here, we present a novel strategy to measure cytotoxicity that is direct, sensitive, rapid, and highly adaptable. Moreover, it allows accurate measurement of viability of both target and effector cells. Target cells are fluorescently labeled with a non-toxic, cell-permeable dye that covalently binds to cell proteins, including nuclear proteins. The labeled target cells are incubated with effector cells to begin killing. Following the killing reaction, the cell mixture is incubated with another dye that specifically stains proteins of dead cells, including nuclear proteins. In the final step, cell nuclei are released by Triton X-100, and analyzed by flow cytometry. This results in four nuclear staining patterns that separate target and effector nuclei as well as nuclei of live and dead cells. Analyzing nuclei, instead of cells, greatly reduces flow cytometry errors caused by the presence of target-effector cell aggregates. Target killing time can often be reduced to 2 h and the assay can be done in a high throughput format. We have successfully validated this assay in a variety of cytotoxicity scenarios including those mediated by NK-92 cells, Chimeric Antigen Receptor (CAR)-T cells, and Tumor Infiltrating Lymphocytes (TIL). Therefore, this technique is broadly applicable, highly sensitive and easily administered, making it a powerful tool to assess immunotherapy-based, cell-mediated cytotoxicity.

Telomere Maintenance-Associated PML Is a Potential Specific Therapeutic Target of Human Colorectal Cancer

Telomere length maintenance is essential for cell proliferation, which is particularly prominent in cancer. We validate that the primary colorectal tumors exhibit heterogeneous telomere lengths but mostly (90%) short telomeres relative to normal tissues. Intriguingly, relatively short telomeres are associated with tumor malignancy as indicated by poorly differentiated state, and these tumors contain more cancer stem-like cells (CSLCs) identified by several commonly used markers CD44, EPHB2 or LGR5. Moreover, promyelocytic leukemia (PML) and ALT-associated PML nuclear bodies (APBs) are frequently found in tumors with short telomeres and high proliferation. In contrast, distant normal tissues rarely or only minimally express PML. Inhibition of PML and APBs by an ATR inhibitor decreases proliferation of CSLCs and organoids, suggesting a potential therapeutic target to progressive colorectal tumors. Together, telomere maintenance underling tumor progression is connected with CSLCs.

hESC-Derived Thalamic Organoids Form Reciprocal Projections When Fused with Cortical Organoids

Human brain organoid techniques have rapidly advanced to facilitate investigating human brain development and diseases. These efforts have largely focused on generating telencephalon due to its direct relevance in a variety of forebrain disorders. Despite its importance as a relay hub between cortex and peripheral tissues, the investigation of three-dimensional (3D) organoid models for the human thalamus has not been explored. Here, we describe a method to differentiate human embryonic stem cells (hESCs) to thalamic organoids (hThOs) that specifically recapitulate the development of thalamus. Single-cell RNA sequencing revealed a formation of distinct thalamic lineages, which diverge from telencephalic fate. Importantly, we developed a 3D system to create the reciprocal projections between thalamus and cortex by fusing the two distinct region-specific organoids representing the developing thalamus or cortex. Our study provides a platform for understanding human thalamic development and modeling circuit organizations and related disorders in the brain.

Correction to: MISC: missing imputation for single-cell RNA sequencing data

It was highlighted that the original article [1] contained a typesetting error in the last name of Allon Canaan. This was incorrectly captured as Allon Canaann in the original article which has since been updated.

Local and global chromatin interactions are altered by large genomic deletions associated with human brain development

Large copy number variants (CNVs) in the human genome are strongly associated with common neurodevelopmental, neuropsychiatric disorders such as schizophrenia and autism. Here we report on the epigenomic effects of the prominent large deletion CNVs on chromosome 22q11.2 and on chromosome 1q21.1. We use Hi-C analysis of long-range chromosome interactions, including haplotype-specific Hi-C analysis, ChIP-Seq analysis of regulatory histone marks, and RNA-Seq analysis of gene expression patterns. We observe changes on all the levels of analysis, within the deletion boundaries, in the deletion flanking regions, along chromosome 22q, and genome wide. We detect gene expression changes as well as pronounced and multilayered effects on chromatin states, chromosome folding and on the topological domains of the chromatin, that emanate from the large CNV locus. These findings suggest basic principles of how such large genomic deletions can alter nuclear organization and affect genomic molecular activity.

Copy number variants in the human genome (CNVs) are associated with neurodevelopmental and psychiatric disorders such as schizophrenia and autism. Here the authors investigate how  the large deletion CNV on chromosome 22q11.2 alters chromatin organization.

MISC: missing imputation for single-cell RNA sequencing data

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) technology provides an effective way to study cell heterogeneity. However, due to the low capture efficiency and stochastic gene expression, scRNA-seq data often contains a high percentage of missing values. It has been showed that the missing rate can reach approximately 30% even after noise reduction. To accurately recover missing values in scRNA-seq data, we need to know where the missing data is; how much data is missing; and what are the values of these data.

METHODS

To solve these three problems, we propose a novel model with a hybrid machine learning method, namely, missing imputation for single-cell RNA-seq (MISC). To solve the first problem, we transformed it to a binary classification problem on the RNA-seq expression matrix. Then, for the second problem, we searched for the intersection of the classification results, zero-inflated model and false negative model results. Finally, we used the regression model to recover the data in the missing elements.

RESULTS

We compared the raw data without imputation, the mean-smooth neighbor cell trajectory, MISC on chronic myeloid leukemia data (CML), the primary somatosensory cortex and the hippocampal CA1 region of mouse brain cells. On the CML data, MISC discovered a trajectory branch from the CP-CML to the BC-CML, which provides direct evidence of evolution from CP to BC stem cells. On the mouse brain data, MISC clearly divides the pyramidal CA1 into different branches, and it is direct evidence of pyramidal CA1 in the subpopulations. In the meantime, with MISC, the oligodendrocyte cells became an independent group with an apparent boundary.

CONCLUSIONS

Our results showed that the MISC model improved the cell type classification and could be instrumental to study cellular heterogeneity. Overall, MISC is a robust missing data imputation model for single-cell RNA-seq data.

Ritornello: high fidelity control-free chromatin immunoprecipitation peak calling

With the advent of next generation high-throughput DNA sequencing technologies, omics experiments have become the mainstay for studying diverse biological effects on a genome wide scale. Chromatin immunoprecipitation (ChIP-seq) is the omics technique that enables genome wide localization of transcription factor (TF) binding or epigenetic modification events. Since the inception of ChIP-seq in 2007, many methods have been developed to infer ChIP-target binding loci from the resultant reads after mapping them to a reference genome. However, interpreting these data has proven challenging, and as such these algorithms have several shortcomings, including susceptibility to false positives due to artifactual peaks, poor localization of binding sites and the requirement for a total DNA input control which increases the cost of performing these experiments. We present Ritornello, a new approach for finding TF-binding sites in ChIP-seq, with roots in digital signal processing that addresses all of these problems. We show that Ritornello generally performs equally or better than the peak callers tested and recommended by the ENCODE consortium, but in contrast, Ritornello does not require a matched total DNA input control to avoid false positives, effectively decreasing the sequencing cost to perform ChIP-seq. Ritornello is freely available at https://github.com/KlugerLab/Ritornello.

Bisulfite-independent analysis of CpG island methylation enables genome-scale stratification of single cells

Conventional DNA bisulfite sequencing has been extended to single cell level, but the coverage consistency is insufficient for parallel comparison. Here we report a novel method for genome-wide CpG island (CGI) methylation sequencing for single cells (scCGI-seq), combining methylation-sensitive restriction enzyme digestion and multiple displacement amplification for selective detection of methylated CGIs. We applied this method to analyzing single cells from two types of hematopoietic cells, K562 and GM12878 and small populations of fibroblasts and induced pluripotent stem cells. The method detected 21 798 CGIs (76% of all CGIs) per cell, and the number of CGIs consistently detected from all 16 profiled single cells was 20 864 (72.7%), with 12 961 promoters covered. This coverage represents a substantial improvement over results obtained using single cell reduced representation bisulfite sequencing, with a 66-fold increase in the fraction of consistently profiled CGIs across individual cells. Single cells of the same type were more similar to each other than to other types, but also displayed epigenetic heterogeneity. The method was further validated by comparing the CpG methylation pattern, methylation profile of CGIs/promoters and repeat regions and 41 classes of known regulatory markers to the ENCODE data. Although not every minor methylation differences between cells are detectable, scCGI-seq provides a solid tool for unsupervised stratification of a heterogeneous cell population.

Fusion of Regionally Specified hPSC-Derived Organoids Models Human Brain Development and Interneuron Migration

Organoid techniques provide unique platforms to model brain development and neurological disorders. Whereas several methods for recapitulating corticogenesis have been described, a system modeling human medial ganglionic eminence (MGE) development, a critical ventral brain domain producing cortical interneurons and related lineages, has been lacking until recently. Here, we describe the generation of MGE and cortex-specific organoids from human pluripotent stem cells that recapitulate the development of MGE and cortex domains, respectively. Population and single-cell RNA sequencing (RNA-seq) profiling combined with bulk assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) analyses revealed transcriptional and chromatin accessibility dynamics and lineage relationships during MGE and cortical organoid development. Furthermore, MGE and cortical organoids generated physiologically functional neurons and neuronal networks. Finally, fusing region-specific organoids followed by live imaging enabled analysis of human interneuron migration and integration. Together, our study provides a platform for generating domain-specific brain organoids and modeling human interneuron migration and offers deeper insight into molecular dynamics during human brain development.

Detecting heterogeneity in single-cell RNA-Seq data by non-negative matrix factorization

Single-cell RNA-Sequencing (scRNA-Seq) is a fast-evolving technology that enables the understanding of biological processes at an unprecedentedly high resolution. However, well-suited bioinformatics tools to analyze the data generated from this new technology are still lacking. Here we investigate the performance of non-negative matrix factorization (NMF) method to analyze a wide variety of scRNA-Seq datasets, ranging from mouse hematopoietic stem cells to human glioblastoma data. In comparison to other unsupervised clustering methods including K-means and hierarchical clustering, NMF has higher accuracy in separating similar groups in various datasets. We ranked genes by their importance scores (D-scores) in separating these groups, and discovered that NMF uniquely identifies genes expressed at intermediate levels as top-ranked genes. Finally, we show that in conjugation with the modularity detection method FEM, NMF reveals meaningful protein-protein interaction modules. In summary, we propose that NMF is a desirable method to analyze heterogeneous single-cell RNA-Seq data. The NMF based subpopulation detection package is available at: https://github.com/lanagarmire/NMFEM.

Evolution and heterogeneity of non-hereditary colorectal cancer revealed by single-cell exome sequencing

Recently single-cell whole-exome sequencing (scWES) has deeply expanded and sharpened our knowledge of cancer evolution and subclonality. Herein, with scWES and matched bulk whole-exome sequencing (bulk WES) on two colorectal cancer (CRC) patients with normal or adenomatous polyps, we found that both the adenoma and cancer were of monoclonal origin, and both shared partial mutations in the same signaling pathways, but each showed a specific spectrum of heterogeneous somatic mutations. In addition, the adenoma and cancer further developed intratumor heterogeneity with the accumulation of nonrandom somatic mutations specifically in GPCR, PI3K-Akt and FGFR signaling pathways. We identified novel driver mutations that developed during adenoma and cancer evolution, particularly in OR1B1 (GPCR signaling pathway) for adenoma evolution, and LAMA1 (PI3K-Akt signaling pathway) and ADCY3 (FGFR signaling pathway) for CRC evolution. In summary, we demonstrated that both colorectal adenoma and CRC are monoclonal in origin, and the CRCs further diversified into different subclones with heterogeneous mutation profiles accumulating in GPCR, PI3K-Akt and FGFR signaling pathways. ScWES provides evidence for the importance of mutations in certain pathways that would not be as apparent from bulk sequencing of tumors, and can potentially establish whether specific mutations are mutually exclusive or occur sequentially in the same subclone of cells.

Regulation of the DNA Methylation Landscape in Human Somatic Cell Reprogramming by the miR-29 Family

Reprogramming to pluripotency after overexpression of OCT4, SOX2, KLF4, and MYC is accompanied by global genomic and epigenomic changes. Histone modification and DNA methylation states in induced pluripotent stem cells (iPSCs) have been shown to be highly similar to embryonic stem cells (ESCs). However, epigenetic differences still exist between iPSCs and ESCs. In particular, aberrant DNA methylation states found in iPSCs are a major concern when using iPSCs in a clinical setting. Thus, it is critical to find factors that regulate DNA methylation states in reprogramming. Here, we found that the miR-29 family is an important epigenetic regulator during human somatic cell reprogramming. Our global DNA methylation and hydroxymethylation analysis shows that DNA demethylation is a major event mediated by miR-29a depletion during early reprogramming, and that iPSCs derived from miR-29a depletion are epigenetically closer to ESCs. Our findings uncover an important miRNA-based approach to generate clinically robust iPSCs.

A Molecular Chipper technology for CRISPR sgRNA library generation and functional mapping of noncoding regions

Clustered regularly-interspaced palindromic repeats (CRISPR)-based genetic screens using single-guide-RNA (sgRNA) libraries have proven powerful to identify genetic regulators. Applying CRISPR screens to interrogate functional elements in noncoding regions requires generating sgRNA libraries that are densely covering, and ideally inexpensive, easy to implement and flexible for customization. Here we present a Molecular Chipper technology for generating dense sgRNA libraries for genomic regions of interest, and a proof-of-principle screen that identifies novel cis-regulatory domains for miR-142 biogenesis. The Molecular Chipper approach utilizes a combination of random fragmentation and a type III restriction enzyme to derive a densely covering sgRNA library from input DNA. Applying this approach to 17 microRNAs and their flanking regions and with a reporter for miR-142 activity, we identify both the pre-miR-142 region and two previously unrecognized cis-domains important for miR-142 biogenesis, with the latter regulating miR-142 processing. This strategy will be useful for identifying functional noncoding elements in mammalian genomes.

Single-Cell Sequencing for Precise Cancer Research: Progress and Prospects

Advances in genomic technology have enabled the faithful detection and measurement of mutations and the gene expression profile of cancer cells at the single-cell level. Recently, several single-cell sequencing methods have been developed that permit the comprehensive and precise analysis of the cancer-cell genome, transcriptome, and epigenome. The use of these methods to analyze cancer cells has led to a series of unanticipated discoveries, such as the high heterogeneity and stochastic changes in cancer-cell populations, the new driver mutations and the complicated clonal evolution mechanisms, and the novel identification of biomarkers of variant tumors. These methods and the knowledge gained from their utilization could potentially improve the early detection and monitoring of rare cancer cells, such as circulating tumor cells and disseminated tumor cells, and promote the development of personalized and highly precise cancer therapy. Here, we discuss the current methods for single cancer-cell sequencing, with a strong focus on those practically used or potentially valuable in cancer research, including single-cell isolation, whole genome and transcriptome amplification, epigenome profiling, multi-dimensional sequencing, and next-generation sequencing and analysis. We also examine the current applications, challenges, and prospects of single cancer-cell sequencing.

A comparative encyclopedia of DNA elements in the mouse genome

The laboratory mouse shares the majority of its protein-coding genes with humans, making it the premier model organism in biomedical research, yet the two mammals differ in significant ways. To gain greater insights into both shared and species-specific transcriptional and cellular regulatory programs in the mouse, the Mouse ENCODE Consortium has mapped transcription, DNase I hypersensitivity, transcription factor binding, chromatin modifications and replication domains throughout the mouse genome in diverse cell and tissue types. By comparing with the human genome, we not only confirm substantial conservation in the newly annotated potential functional sequences, but also find a large degree of divergence of sequences involved in transcriptional regulation, chromatin state and higher order chromatin organization. Our results illuminate the wide range of evolutionary forces acting on genes and their regulatory regions, and provide a general resource for research into mammalian biology and mechanisms of human diseases.

Co-detection and sequencing of genes and transcripts from the same single cells facilitated by a microfluidics platform

Despite the recent advance of single-cell gene expression analyses, co-measurement of both genomic and transcriptional signatures at the single-cell level has not been realized. However such analysis is necessary in order to accurately delineate how genetic information is transcribed, expressed, and regulated to give rise to an enormously diverse range of cell phenotypes. Here we report on a microfluidics-facilitated approach that allows for controlled separation of cytoplasmic and nuclear contents of a single cell followed by on-chip amplification of genomic DNA and cytoplasmic mRNA. When coupled with off-chip polymerase chain reaction, gel electrophoresis and Sanger sequencing, a panel of genes and transcripts from the same single cell can be co-detected and sequenced. This platform is potentially an enabling tool to permit multiple genomic measurements performed on the same single cells and opens new opportunities to tackle a range of fundamental biology questions including non-genetic cell-to-cell variability, epigenetic regulation, and stem cell fate control. It also helps address clinical challenges such as diagnosing intra-tumor heterogeneity and dissecting complex cellular immune responses.

Long intergenic non-coding RNA HOTAIRM1 regulates cell cycle progression during myeloid maturation in NB4 human promyelocytic leukemia cells

HOTAIRM1 is a long intergenic non-coding RNA encoded in the human HOXA gene cluster, with gene expression highly specific for maturing myeloid cells. Knockdown of HOTAIRM1 in the NB4 acute promyelocytic leukemia cell line retarded all-trans retinoid acid (ATRA)-induced granulocytic differentiation, resulting in a significantly larger population of immature and proliferating cells that maintained cell cycle progression from G1 to S phases. Correspondingly, HOTAIRM1 knockdown resulted in retained expression of many otherwise ATRA-suppressed cell cycle and DNA replication genes, and abated ATRA induction of cell surface leukocyte activation, defense response, and other maturation-related genes. Resistance to ATRA-induced cell cycle arrest at the G1/S phase transition in knockdown cells was accompanied by retained expression of ITGA4 (CD49d) and decreased induction of ITGAX (CD11c). The coupling of cell cycle progression with temporal dynamics in the expression patterns of these integrin genes suggests a regulated switch to control the transit from the proliferative phase to granulocytic maturation. Furthermore, ITGAX was among a small number of genes showing perturbation in transcript levels upon HOTAIRM1 knockdown even without ATRA treatment, suggesting a direct pathway of regulation. These results indicate that HOTAIRM1 provides a regulatory link in myeloid maturation by modulating integrin-controlled cell cycle progression at the gene expression level.

Nonstochastic reprogramming from a privileged somatic cell state

Reprogramming somatic cells to induced pluripotency by Yamanaka factors is usually slow and inefficient and is thought to be a stochastic process. We identified a privileged somatic cell state, from which acquisition of pluripotency could occur in a nonstochastic manner. Subsets of murine hematopoietic progenitors are privileged whose progeny cells predominantly adopt the pluripotent fate with activation of endogenous Oct4 locus after four to five divisions in reprogramming conditions. Privileged cells display an ultrafast cell cycle of ∼8 hr. In fibroblasts, a subpopulation cycling at a similar ultrafast speed is observed after 6 days of factor expression and is increased by p53 knockdown. This ultrafast cycling population accounts for >99% of the bulk reprogramming activity in wild-type or p53 knockdown fibroblasts. Our data demonstrate that the stochastic nature of reprogramming can be overcome in a privileged somatic cell state and suggest that cell-cycle acceleration toward a critical threshold is an important bottleneck for reprogramming. PAPERCLIP:

A highly integrated and complex PPARGC1A transcription factor binding network in HepG2 cells

PPARGC1A is a transcriptional coactivator that binds to and coactivates a variety of transcription factors (TFs) to regulate the expression of target genes. PPARGC1A plays a pivotal role in regulating energy metabolism and has been implicated in several human diseases, most notably type II diabetes. Previous studies have focused on the interplay between PPARGC1A and individual TFs, but little is known about how PPARGC1A combines with all of its partners across the genome to regulate transcriptional dynamics. In this study, we describe a core PPARGC1A transcriptional regulatory network operating in HepG2 cells treated with forskolin. We first mapped the genome-wide binding sites of PPARGC1A using chromatin-IP followed by high-throughput sequencing (ChIP-seq) and uncovered overrepresented DNA sequence motifs corresponding to known and novel PPARGC1A network partners. We then profiled six of these site-specific TF partners using ChIP-seq and examined their network connectivity and combinatorial binding patterns with PPARGC1A. Our analysis revealed extensive overlap of targets including a novel link between PPARGC1A and HSF1, a TF regulating the conserved heat shock response pathway that is misregulated in diabetes. Importantly, we found that different combinations of TFs bound to distinct functional sets of genes, thereby helping to reveal the combinatorial regulatory code for metabolic and other cellular processes. In addition, the different TFs often bound near the promoters and coding regions of each other's genes suggesting an intricate network of interdependent regulation. Overall, our study provides an important framework for understanding the systems-level control of metabolic gene expression in humans.

An encyclopedia of mouse DNA elements (Mouse ENCODE)

To complement the human Encyclopedia of DNA Elements (ENCODE) project and to enable a broad range of mouse genomics efforts, the Mouse ENCODE Consortium is applying the same experimental pipelines developed for human ENCODE to annotate the mouse genome.

A multiprotein complex necessary for both transcription and DNA replication at the β-globin locus

DNA replication, repair, transcription and chromatin structure are intricately associated nuclear processes, but the molecular links between these events are often obscure. In this study, we have surveyed the protein complexes that bind at β-globin locus control region, and purified and characterized the function of one such multiprotein complex from human erythroleukemic K562 cells. We further validated the existence of this complex in human CD34+ cell-derived normal erythroid cells. This complex contains ILF2/ILF3 transcription factors, p300 acetyltransferase and proteins associated with DNA replication, transcription and repair. RNAi knockdown of ILF2, a DNA-binding component of this complex, abrogates the recruitment of the complex to its cognate DNA sequence and inhibits transcription, histone acetylation and usage of the origin of DNA replication at the β-globin locus. These results imply a direct link between mammalian DNA replication, transcription and histone acetylation mediated by a single multiprotein complex.

Variation in transcription factor binding among humans

Differences in gene expression may play a major role in speciation and phenotypic diversity. We examined genome-wide differences in transcription factor (TF) binding in several humans and a single chimpanzee by using chromatin immunoprecipitation followed by sequencing. The binding sites of RNA polymerase II (PolII) and a key regulator of immune responses, nuclear factor kappaB (p65), were mapped in 10 lymphoblastoid cell lines, and 25 and 7.5% of the respective binding regions were found to differ between individuals. Binding differences were frequently associated with single-nucleotide polymorphisms and genomic structural variants, and these differences were often correlated with differences in gene expression, suggesting functional consequences of binding variation. Furthermore, comparing PolII binding between humans and chimpanzee suggests extensive divergence in TF binding. Our results indicate that many differences in individuals and species occur at the level of TF binding, and they provide insight into the genetic events responsible for these differences.

Dynamics of alpha-globin locus chromatin structure and gene expression during erythroid differentiation of human CD34(+) cells in culture

OBJECTIVE

The aim of the present study has been to establish serum-free culture conditions for ex vivo expansion and differentiation of human CD34(+) cells into erythroid lineage and to study the chromatin structure, gene expression, and transcription factor recruitment at the alpha-globin locus in the developing erythron.

MATERIALS AND METHODS

A basal Iscove's modified Dulbecco's medium cell culture medium with 1% bovine serum albumin as a serum replacement and a combination of cytokines and growth factors was used for expansion and differentiation of the CD34(+) cells. Expression patterns of the alpha- and beta-like genes at various stages of erythropoiesis was studied by reverse transcriptase quantitative polymerase chain reaction analysis, profile of key erythroid transcription factors was investigated by Western blotting, and the chromatin structure and transcription factor recruitment at the alpha-globin locus was investigated by chromatin immunoprecipitation quantitative polymerase chain reaction analysis.

RESULTS

Human CD34(+) cells in the serum-free medium undergo near synchronous erythroid differentiation to yield large amount of cells at different differentiation stages. We observe distinct patterns of the histone modifications and transcription factor binding at the alpha-globin locus during erythroid differentiation of CD34(+) cells. Nuclear factor erythroid-derived 2 (NF-E2) was present at upstream activator sites even before addition of erythropoietin (EPO), while bound GATA-1 was only detectable after EPO treatment. After 7 days of EPO treatment, H3K4Me2 modification uniformly increases throughout the alpha-globin locus. Acetylation at H3K9 and binding of Pol II, NF-E2, and GATA-1 were restricted to certain hypersensitive sites of the enhancer and theta gene, and were conspicuously low at the alpha-like globin promoters. Rearrangement of the insulator binding factor CTCF took place at and around the alpha-globin locus as CD34(+) cells differentiated into erythroid pathway.

CONCLUSION

Our results indicate that remodeling of the upstream elements may be the primary event in activation of alpha-globin gene expression. Activation of alpha-globin genes upon EPO treatment involves initial binding of Pol II, downregulation of pre-existing factors like NF-E2, removal of CTCF from the locus, then rebinding of CTCF in an altered pattern, and concurrent or subsequent binding of transcription factors like GATA-1.

Genome-wide screen of promoter methylation identifies novel markers in melanoma

DNA methylation is an important component of epigenetic modifications, which influences the transcriptional machinery aberrant in many human diseases. In this study we present the first genome-wide integrative analysis of promoter methylation and gene expression for the identification of methylation markers in melanoma. Genome-wide promoter methylation and gene expression of eight early-passage human melanoma cell strains were compared with newborn and adult melanocytes. We used linear mixed effect models (LME) in combination with a series of filters based on the localization of promoter methylation relative to the transcription start site, overall promoter CpG content, and differential gene expression to discover DNA methylation markers. This approach identified 76 markers, of which 68 were hyper- and eight hypomethylated (LME, P < 0.05). Promoter methylation and differential gene expression of five markers (COL1A2, NPM2, HSPB6, DDIT4L, MT1G) were validated by sequencing of bisulfite-modified DNA and real-time reverse transcriptase PCR, respectively. Importantly, the incidence of promoter methylation of the validated markers increased moderately in early and significantly in advanced-stage melanomas, using early-passage cell strains and snap-frozen tissues (n = 18 and n = 24, respectively) compared with normal melanocytes and nevi (n = 11 and n = 9, respectively). Our approach allows robust identification of methylation markers that can be applied to other studies involving genome-wide promoter methylation. In conclusion, this study represents the first unbiased systematic effort to determine methylation markers in melanoma and revealed several novel genes regulated by promoter methylation that were not described in cancer cells before.

Chimeric receptor mRNA transfection as a tool to generate antineoplastic lymphocytes

mRNA transfection is a useful approach for temporal cell reprogramming with minimal risk of transgene-mediated mutagenesis. We applied this to redirect lymphocyte cytotoxicity toward malignant cells. Using the chimeric immune receptor (CIR) constructs anti-CD19 CIR and 8H9 CIR, we achieved uniform expression of CIRs on virtually the entire population of lymphocytes. We reprogrammed CD3+ CD8+, CD3+ CD4+, and natural killer (NK ) cells toward autologous and allogeneic targets such as B cells, Daudi lymphoma, primary melanoma, breast ductal carcinoma, breast adenocarcinoma, and rhabdomyosarcoma. The reprogramming procedure is fast. Although most of the experiments were performed on lymphocytes obtained after 7-day activation, only 1-day activation of T cells with anti-CD3, anti-CD28 antibodies, and interleukin-2 is sufficient to develop both lymphocyte cytotoxicity and competence for mRNA transfer. The entire procedure, which includes lymphocyte activation and reprogramming, can be completed in 2 days. The efficiency of mRNA-modified human T cells was tested in a murine xenograft model. Human CD3+CD8+ lymphocytes expressing anti-CD19 CIR mRNA inhibited Daudi lymphoma growth in NOD=SCID mice. These results demonstrate that a mixed population of cytotoxic lymphocytes, including T cells together with NK cells, can be quickly and simultaneously reprogrammed by mRNA against autologous malignancies. With relatively minor modifications the described method of lymphocyte reprogramming can be scaled up for cancer therapy.

Genome-wide occupancy of SREBP1 and its partners NFY and SP1 reveals novel functional roles and combinatorial regulation of distinct classes of genes

The sterol regulatory element-binding protein (SREBP) family member SREBP1 is a critical transcriptional regulator of cholesterol and fatty acid metabolism and has been implicated in insulin resistance, diabetes, and other diet-related diseases. We globally identified the promoters occupied by SREBP1 and its binding partners NFY and SP1 in a human hepatocyte cell line using chromatin immunoprecipitation combined with genome tiling arrays (ChIP-chip). We find that SREBP1 occupies the promoters of 1,141 target genes involved in diverse biological pathways, including novel targets with roles in lipid metabolism and insulin signaling. We also identify a conserved SREBP1 DNA-binding motif in SREBP1 target promoters, and we demonstrate that many SREBP1 target genes are transcriptionally activated by treatment with insulin and glucose using gene expression microarrays. Finally, we show that SREBP1 cooperates extensively with NFY and SP1 throughout the genome and that unique combinations of these factors target distinct functional pathways. Our results provide insight into the regulatory circuitry in which SREBP1 and its network partners coordinate a complex transcriptional response in the liver with cues from the diet.

Transcription factors (TFs) are DNA-binding proteins that regulate the transcription of their target genes. TFs typically bind in proximity to the transcription start sites of their target genes in a region called the promoter. SREBP1 is a TF that increases the transcription of numerous genes involved in cholesterol and fat metabolism and has been linked to diet-related diseases such as insulin resistance and type 2 diabetes. Using microarray technology, we identified all of the promoters in the human genome that are bound by SREBP1 and two associated TFs called NFY and SP1 in a human liver cell line. Our findings greatly expand the number of genes and biological pathways that may be regulated by SREBP1 and reveal that different combinations of SREBP1 and its partners preferentially target genes involved in different pathways. Thus, in contrast to traditional studies that focus on individual genes, we have used a genomics approach to provide a novel global view of the regulatory circuitry in which SREBP1 and its partners coordinate a transcriptional response in the liver with cues from the diet.

A genomic analysis of RNA polymerase II modification and chromatin architecture related to 3' end RNA polyadenylation

Genomic analyses have been applied extensively to analyze the process of transcription initiation in mammalian cells, but less to transcript 3' end formation and transcription termination. We used a novel approach to prepare 3' end fragments from polyadenylated RNA, and mapped the position of the poly(A) addition site using oligonucleotide arrays tiling 1% of the human genome. This approach revealed more 3' ends than had been annotated. The distribution of these ends relative to RNA polymerase II (PolII) and di- and trimethylated lysine 4 and lysine 36 of histone H3 was compared. A substantial fraction of unannotated 3' ends of RNA are intronic and antisense to the embedding gene. Poly(A) ends of annotated messages lie on average 2 kb upstream of the end of PolII binding (termination). Near the termination sites, and in some internal sites, unphosphorylated and C-terminal domain (CTD) serine 2 phosphorylated PolII (POLR2A) accumulate, suggesting pausing of the polymerase and perhaps dephosphorylation prior to release. Lysine 36 trimethylation occurs across transcribed genes, sometimes alternating with stretches of DNA in which lysine 36 dimethylation is more prominent. Lysine 36 methylation decreases at or near the site of polyadenylation, sometimes disappearing before disappearance of phosphorylated RNA PolII or release of PolII from DNA. Our results suggest that transcription termination loss of histone 3 lysine 36 methylation and later release of RNA polymerase. The latter is often associated with polymerase pausing. Overall, our study reveals extensive sites of poly(A) addition and provides insights into the events that occur during 3' end formation.

Control of beta globin genes

The developmental changes in expression of the beta like genes from embryonic to adult stages of human life are controlled at least partially at the level of the promoter sequences of these genes and their binding factors, and competition for promoter specific interactions with the locus control region (LCR). In recent years, the control of beta globin genes has also been investigated at the level of chromatin structure involving the chemical modification of histones and their remodelling by DNA dependent ATPases (SMARCA) containing protein complexes. The role of intergenic RNA is also being investigated with renewed interest. Although a wealth of information on the structure/function relationship of the LCR and globin promoters has been gathered over more than two decades, the fundamental nature of the control of these genes at the molecular level is still not completely understood. In the following pages, we intend to briefly describe the progress made in the field and discuss future directions.

Paired-end mapping reveals extensive structural variation in the human genome

Structural variation of the genome involves kilobase- to megabase-sized deletions, duplications, insertions, inversions, and complex combinations of rearrangements. We introduce high-throughput and massive paired-end mapping (PEM), a large-scale genome-sequencing method to identify structural variants (SVs) approximately 3 kilobases (kb) or larger that combines the rescue and capture of paired ends of 3-kb fragments, massive 454 sequencing, and a computational approach to map DNA reads onto a reference genome. PEM was used to map SVs in an African and in a putatively European individual and identified shared and divergent SVs relative to the reference genome. Overall, we fine-mapped more than 1000 SVs and documented that the number of SVs among humans is much larger than initially hypothesized; many of the SVs potentially affect gene function. The breakpoint junction sequences of more than 200 SVs were determined with a novel pooling strategy and computational analysis. Our analysis provided insights into the mechanisms of SV formation in humans.

Systematic prediction and validation of breakpoints associated with copy-number variants in the human genome

Copy-number variants (CNVs) are an abundant form of genetic variation in humans. However, approaches for determining exact CNV breakpoint sequences (physical deletion or duplication boundaries) across individuals, crucial for associating genotype to phenotype, have been lacking so far, and the vast majority of CNVs have been reported with approximate genomic coordinates only. Here, we report an approach, called BreakPtr, for fine-mapping CNVs (available from http://breakptr.gersteinlab.org). We statistically integrate both sequence characteristics and data from high-resolution comparative genome hybridization experiments in a discrete-valued, bivariate hidden Markov model. Incorporation of nucleotide-sequence information allows us to take into account the fact that recently duplicated sequences (e.g., segmental duplications) often coincide with breakpoints. In anticipation of an upcoming increase in CNV data, we developed an iterative, "active" approach to initially scoring with a preliminary model, performing targeted validations, retraining the model, and then rescoring, and a flexible parameterization system that intuitively collapses from a full model of 2,503 parameters to a core one of only 10. Using our approach, we accurately mapped >400 breakpoints on chromosome 22 and a region of chromosome 11, refining the boundaries of many previously approximately mapped CNVs. Four predicted breakpoints flanked known disease-associated deletions. We validated an additional four predicted CNV breakpoints by sequencing. Overall, our results suggest a predictive resolution of approximately 300 bp. This level of resolution enables more precise correlations between CNVs and across individuals than previously possible, allowing the study of CNV population frequencies. Further, it enabled us to demonstrate a clear Mendelian pattern of inheritance for one of the CNVs.

Synthetic messenger RNA as a tool for gene therapy

Transfection of human cells with DNA in biomedical applications carries the risk of insertional mutagenesis. Transfection with mRNA avoids this problem; however, in vitro production of mRNA, based on preliminary DNA template cloning in special vectors, is a laborious and time-consuming procedure. We report an efficient vectorfree method of mRNA production from polymerase chain reaction-generated DNA templates. For all cell types tested mRNA was transfected more readily than DNA, and its expression was highly uniform in cell populations. Even cell types relatively resistant to transfection with DNA could express transfected mRNA well. The level of mRNA expression could be controlled over a wide range by changing the amount of input RNA. Cells could be efficiently and simultaneously loaded with several different transcripts. To test a potential clinical application of this method, we transfected human T lymphocytes with mRNA encoding a chimeric immune receptor directed against CD19, a surface antigen widely expressed in leukemia and lymphoma. The transfected mRNA conferred powerful cytotoxicity to T cells against CD19+ targets from the same donor. These results demonstrate that this method can be applied to generate autologous T lymphocytes directed toward malignant cells.

Differentiated cells are more efficient than adult stem cells for cloning by somatic cell nuclear transfer

Since the creation of Dolly via somatic cell nuclear transfer (SCNT), more than a dozen species of mammals have been cloned using this technology. One hypothesis for the limited success of cloning via SCNT (1%-5%) is that the clones are likely to be derived from adult stem cells. Support for this hypothesis comes from the findings that the reproductive cloning efficiency for embryonic stem cells is five to ten times higher than that for somatic cells as donors and that cloned pups cannot be produced directly from cloned embryos derived from differentiated B and T cells or neuronal cells. The question remains as to whether SCNT-derived animal clones can be derived from truly differentiated somatic cells. We tested this hypothesis with mouse hematopoietic cells at different differentiation stages: hematopoietic stem cells, progenitor cells and granulocytes. We found that cloning efficiency increases over the differentiation hierarchy, and terminally differentiated postmitotic granulocytes yield cloned pups with the greatest cloning efficiency.

FAT10/diubiquitin-like protein-deficient mice exhibit minimal phenotypic differences

The FAT10 gene encodes a diubiquitin-like protein containing two tandem head-to-tail ubiquitin-like domains. There is a high degree of similarity between murine and human FAT10 sequences at both the mRNA and protein levels. In various cell lines, FAT10 expression was shown to be induced by gamma interferon or by tumor necrosis factor alpha. In addition, FAT10 expression was found to be up-regulated in some Epstein-Barr virus-infected B-cell lines, in activated dendritic cells, and in several epithelial tumors. However, forced expression of FAT10 in cultured cells was also found to produce apoptotic cell death. Overall, these findings suggest that FAT10 may modulate cellular growth or cellular viability. Here we describe the steps to generate, by genetic targeting, a FAT10 gene knockout mouse model. The FAT10 knockout homozygous mice are viable and fertile. No gross lesions or obvious histological differences were found in these mutated mice. Examination of lymphocyte populations from spleen, thymus, and bone marrow did not reveal any abnormalities. However, flow cytometry analysis demonstrated that the lymphocytes of FAT10 knockout mice were, on average, more prone to spontaneous apoptotic death. Physiologically, these mice demonstrated a high level of sensitivity toward endotoxin challenge. These findings indicate that FAT10 may function as a survival factor.

Identification of coding single nucleotide polymorphisms and mutations by combination of genome tiling arrays and enrichment/depletion of mismatch cDNAs

Genome tiling array technology combined with a method for both enrichment and depletion of mismatch-containing cDNA fragments offers a useful approach for detecting coding single nucleotide polymorphisms (cSNPs) and mutations in pooled cDNA samples. Enriched mismatch and perfect match cDNA samples from human primary melanoma cells and normal melanocytes were obtained by selection using mismatch repair thymine DNA glycosylase-bound beads. These cDNA samples were then labeled and hybridized to Encyclopedia of DNA Elements genome tiling arrays. The results revealed that the hybridization intensity values of potential cDNA variation regions of the enriched mismatch samples increased, whereas the hybridization intensity values of corresponding regions of the enriched perfect match samples decreased. Six potential mutations were confirmed by polymerase chain reaction product sequencing, including two novel heterozygous mutations in melanoma cells. We suggest that this strategy should increase the efficiency of both cSNP and mutation detection throughout the entire human genome and decrease the cost and complexity of genomewide analysis of cDNA variations.

High-resolution mapping of DNA copy alterations in human chromosome 22 using high-density tiling oligonucleotide arrays

Deletions and amplifications of the human genomic sequence (copy number polymorphisms) are the cause of numerous diseases and a potential cause of phenotypic variation in the normal population. Comparative genomic hybridization (CGH) has been developed as a useful tool for detecting alterations in DNA copy number that involve blocks of DNA several kilobases or larger in size. We have developed high-resolution CGH (HR-CGH) to detect accurately and with relatively little bias the presence and extent of chromosomal aberrations in human DNA. Maskless array synthesis was used to construct arrays containing 385,000 oligonucleotides with isothermal probes of 45-85 bp in length; arrays tiling the beta-globin locus and chromosome 22q were prepared. Arrays with a 9-bp tiling path were used to map a 622-bp heterozygous deletion in the beta-globin locus. Arrays with an 85-bp tiling path were used to analyze DNA from patients with copy number changes in the pericentromeric region of chromosome 22q. Heterozygous deletions and duplications as well as partial triploidies and partial tetraploidies of portions of chromosome 22q were mapped with high resolution (typically up to 200 bp) in each patient, and the precise breakpoints of two deletions were confirmed by DNA sequencing. Additional peaks potentially corresponding to known and novel additional CNPs were also observed. Our results demonstrate that HR-CGH allows the detection of copy number changes in the human genome at an unprecedented level of resolution.

Multi-protein complexes at the beta-globin locus

Biochemical analysis of the beta-globin gene function has led to the identification of several multi-protein complexes at the locus control region (LCR), insulator and promoters. This review briefly summarizes these multi-protein complexes and discusses their contribution towards the regulation of the beta-globin gene expression.

A limited number of genes are involved in the differentiation of germinal center B cells

Mature B cells, upon activation, progressively differentiate through centroblasts into centrocytes and finally to plasmacytes that express large amounts of selected immunoglobulins. A significant part of this maturation is thought to involve induction of the unfolded protein response (UPR). We have compared gene expression in normal germinal center centroblasts, centrocytes, lymphoblastoid cells undergoing induced UPR, and the CCL155 plasmacytoma cell line. In the centroblast to centrocyte transition there is a change in the expression of a relatively small number of genes. These include a limited subset of the genes upregulated by a fully activated UPR as well as a small number of other transcription factors, some disulphide isomerases, and other genes. This is consistent with a model in which this transition is mediated by changes in the levels of expression of transcription factor B-lymphocyte-induced maturation protein 1 (Blimp1) (PRDM1), BACH2, X-box binding protein 1 (XBP1), interferon regulatory factor 4 (IRF4), and possibly vitamin D receptor (VDR) expression, together with post-transcriptional changes and a limited induction of aspects of the UPR.

Two types of precursor cells in a multipotential hematopoietic cell line

The biochemistry of early stages of hematopoietic differentiation is difficult to study because only relatively small numbers of precursor cells are available. The murine EML cell line is a multipotential cell line that can be used to model some of these steps. We found that the lineage- EML precursor cells can be separated into two populations based on cell surface markers including CD34. Both populations contain similar levels of stem cell factor (SCF) receptor (c-Kit) but only the CD34+ population shows a growth response when treated with SCF. Conversely, the CD34- population will grow in the presence of the cytokine IL-3. The human beta-globin locus control region hypersensitive site 2 plays different roles on beta-globin transcription in the CD34+ and CD34- populations. The two populations are present in about equal amounts in culture, and the CD34+ population rapidly regenerates the mixed population when grown in the presence of SCF. We suggest that this system may mimic a normal developmental transition in hematopoiesis.

Werner protein protects nonproliferating cells from oxidative DNA damage

Werner syndrome, caused by mutations of the WRN gene, mimics many changes of normal aging. Although roles for WRN protein in DNA replication, recombination, and telomere maintenance have been suggested, the pathology of rapidly dividing cells is not a feature of Werner syndrome. To identify cellular events that are specifically vulnerable to WRN deficiency, we used RNA interference (RNAi) to knockdown WRN or BLM (the RecQ helicase mutated in Bloom syndrome) expression in primary human fibroblasts. Withdrawal of WRN or BLM produced accelerated cellular senescence phenotype and DNA damage response in normal fibroblasts, as evidenced by induction of gammaH2AX and 53BP1 nuclear foci. After WRN depletion, the induction of these foci was seen most prominently in nondividing cells. Growth in physiological (3%) oxygen or in the presence of an antioxidant prevented the development of the DNA damage foci in WRN-depleted cells, whereas acute oxidative stress led to inefficient repair of the lesions. Furthermore, WRN RNAi-induced DNA damage was suppressed by overexpression of the telomere-binding protein TRF2. These conditions, however, did not prevent the DNA damage response in BLM-ablated cells, suggesting a distinct role for WRN in DNA homeostasis in vivo. Thus, manifestations of Werner syndrome may reflect an impaired ability of slowly dividing cells to limit oxidative DNA damage.

Heterogeneous nuclear ribonucleoprotein C1/C2, MeCP1, and SWI/SNF form a chromatin remodeling complex at the beta-globin locus control region

Locus control regions (LCRs) are regulatory DNA sequences that are situated many kilobases away from their cognate promoters. LCRs protect transgenes from position effect variegation and heterochromatinization and also promote copy-number dependence of the levels of transgene expression. In this work, we describe the biochemical purification of a previously undescribed LCR-associated remodeling complex (LARC) that consists of heterogeneous nuclear ribonucleoprotein C1/C2, nucleosome remodeling SWI/SNF, and nucleosome remodeling and deacetylating (NuRD)/MeCP1 as a single homogeneous complex. LARC binds to the hypersensitive 2 (HS2)-Maf recognition element (MARE) DNA in a sequence-specific manner and remodels nucleosomes. Heterogeneous nuclear ribonucleoprotein C1/C2, previously known as a general RNA binding protein, provides a sequence-specific DNA recognition element for LARC, and the LARC DNA-recognition sequence is essential for the enhancement of transcription by HS2. Independently of the initiation of transcription, LARC becomes associated with beta-like globin promoters.

Lineage specificity of gene expression patterns

The hematopoietic system offers many advantages as a model for understanding general aspects of lineage choice and specification. Using oligonucleotide microarrays, we compared gene expression patterns of multiple purified hematopoietic cell populations, including neutrophils, monocytes, macrophages, resting, centrocytic, and centroblastic B lymphocytes, dendritic cells, and hematopoietic stem cells. Some of these cells were studied under both resting and stimulated conditions. We studied the collective behavior of subsets of genes derived from the Biocarta database of functional pathways, hand-tuned groupings of genes into broad functional categories based on the Gene Ontology database, and the metabolic pathways in the Kyoto Encyclopedia of Genes and Genomes database. Principal component analysis revealed strikingly pervasive differences in relative levels of gene expression among cell lineages that involve most of the subsets examined. These results indicate that many processes in these cells behave differently in different lineages. Much of the variation among lineages was captured by the first few principal components. Principal components biplots were found to provide a convenient visual display of the contributions of the various genes within the subsets in lineage discrimination. Moreover, by applying tree-constructing methodologies borrowed from phylogenetics to the expression data from differentiated cells and stem cells, we reconstructed a tree of relationships that resembled the established hematopoietic program of lineage development. Thus, the mRNA expression data implicitly contained information about developmental relationships among cell types.