let-7 miRNAs repress HIC2 to regulate BCL11A transcription and hemoglobin switching

ABSTRACT

The switch from fetal hemoglobin (γ-globin, HBG) to adult hemoglobin (β-globin, HBB) gene transcription in erythroid cells serves as a paradigm for a complex and clinically relevant developmental gene regulatory program. We previously identified HIC2 as a regulator of the switch by inhibiting the transcription of BCL11A, a key repressor of HBG production. HIC2 is highly expressed in fetal cells, but the mechanism of its regulation is unclear. Here we report that HIC2 developmental expression is controlled by microRNAs (miRNAs), as loss of global miRNA biogenesis through DICER1 depletion leads to upregulation of HIC2 and HBG messenger RNA. We identified the adult-expressed let-7 miRNA family as a direct posttranscriptional regulator of HIC2. Ectopic expression of let-7 in fetal cells lowered HIC2 levels, whereas inhibition of let-7 in adult erythroblasts increased HIC2 production, culminating in decommissioning of a BCL11A erythroid enhancer and reduced BCL11A transcription. HIC2 depletion in let-7-inhibited cells restored BCL11A-mediated repression of HBG. Together, these data establish that fetal hemoglobin silencing in adult erythroid cells is under the control of a miRNA-mediated inhibitory pathway (let-7 ⊣ HIC2 ⊣ BCL11A ⊣ HBG).

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Genome folding principles revealed in condensin-depleted mitotic chromosomes

During mitosis, condensin activity interferes with interphase chromatin structures. Here, we generated condensin-free mitotic chromosomes to investigate genome folding principles. Co- depletion of condensin I and II, but neither alone, triggered mitotic chromosome compartmentalization in ways that differ from interphase. Two distinct euchromatic compartments, indistinguishable in interphase, rapidly emerged upon condensin loss with different interaction preferences and dependence on H3K27ac. Constitutive heterochromatin gradually self-aggregated and co-compartmentalized with the facultative heterochromatin, contrasting with their separation during interphase. While topologically associating domains (TADs) and CTCF/cohesin mediated structural loops remained undetectable, cis-regulatory element contacts became apparent, providing an explanation for their quick re-establishment during mitotic exit. HP1 proteins, which are thought to partition constitutive heterochromatin, were absent from mitotic chromosomes, suggesting, surprisingly, that constitutive heterochromatin can self-aggregate without HP1. Indeed, in cells traversing from M- to G1-phase in the combined absence of HP1α, HP1Π and HP1γ, re-established constitutive heterochromatin compartments normally. In sum, "clean-slate" condensin-deficient mitotic chromosomes illuminate mechanisms of genome compartmentalization not revealed in interphase cells.

Molecular basis of polycomb group protein-mediated fetal hemoglobin repression

The switch from fetal hemoglobin (HbF) to adult hemoglobin (HbA) is a paradigm for developmental gene expression control with relevance to sickle cell disease and β-thalassemia. Polycomb repressive complex (PRC) proteins regulate this switch, and an inhibitor of PRC2 has entered a clinical trial for HbF activation. Yet, how PRC complexes function in this process, their target genes, and relevant subunit composition are unknown. Here, we identified the PRC1 subunit BMI1 as a novel HbF repressor. We uncovered the RNA binding proteins LIN28B, IGF2BP1, and IGF2BP3 genes as direct BMI1 targets, and demonstrate that they account for the entirety of BMI1's effect on HbF regulation. BMI1 functions as part of the canonical PRC1 (cPRC1) subcomplex as revealed by the physical and functional dissection of BMI1 protein partners. Lastly, we demonstrate that BMI1/cPRC1 acts in concert with PRC2 to repress HbF through the same target genes. Our study illuminates how PRC silences HbF, highlighting an epigenetic mechanism involved in hemoglobin switching.

Snapshot: a package for clustering and visualizing epigenetic history during cell differentiation

BACKGROUND

Epigenetic modification of chromatin plays a pivotal role in regulating gene expression during cell differentiation. The scale and complexity of epigenetic data pose significant challenges for biologists to identify the regulatory events controlling cell differentiation.

RESULTS

To reduce the complexity, we developed a package, called Snapshot, for clustering and visualizing candidate cis-regulatory elements (cCREs) based on their epigenetic signals during cell differentiation. This package first introduces a binarized indexing strategy for clustering the cCREs. It then provides a series of easily interpretable figures for visualizing the signal and epigenetic state patterns of the cCREs clusters during the cell differentiation. It can also use different hierarchies of cell types to highlight the epigenetic history specific to any particular cell lineage. We demonstrate the utility of Snapshot using data from a consortium project for ValIdated Systematic IntegratiON (VISION) of epigenomic data in hematopoiesis.

CONCLUSION

The package Snapshot can identify all distinct clusters of genomic locations with unique epigenetic signal patterns during cell differentiation. It outperforms other methods in terms of interpreting and reproducing the identified cCREs clusters. The package of Snapshot is available at GitHub: https://github.com/guanjue/Snapshot .

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

CTCF blocks antisense transcription initiation at divergent promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

CTCF blocks antisense transcription initiation at divergent promoters

Transcription at most promoters is divergent, initiating at closely spaced oppositely oriented core promoters to produce sense transcripts along with often unstable upstream antisense transcripts (uasTrx). How antisense transcription is regulated and to what extent it is coordinated with sense transcription is not well understood. Here, by combining acute degradation of the multi-functional transcription factor CTCF and nascent transcription measurements, we find that CTCF specifically suppresses antisense but not sense transcription at hundreds of divergent promoters. Primary transcript RNA-FISH shows that CTCF lowers burst fraction but not burst intensity of uasTrx and that co-bursting of sense and antisense transcripts is disfavored. Genome editing, chromatin conformation studies and high-resolution transcript mapping revealed that precisely positioned CTCF directly suppresses the initiation of uasTrx, in a manner independent of its architectural function. In sum, CTCF shapes the transcriptional landscape in part by suppressing upstream antisense transcription.

Subtype-specific 3D genome alteration in acute myeloid leukaemia

Acute myeloid leukaemia (AML) represents a set of heterogeneous myeloid malignancies, and hallmarks include mutations in epigenetic modifiers, transcription factors and kinases1-5. The extent to which mutations in AML drive alterations in chromatin 3D structure and contribute to myeloid transformation is unclear. Here we use Hi-C and whole-genome sequencing to analyse 25 samples from patients with AML and 7 samples from healthy donors. Recurrent and subtype-specific alterations in A/B compartments, topologically associating domains and chromatin loops were identified. RNA sequencing, ATAC with sequencing and CUT&Tag for CTCF, H3K27ac and H3K27me3 in the same AML samples also revealed extensive and recurrent AML-specific promoter-enhancer and promoter-silencer loops. We validated the role of repressive loops on their target genes by CRISPR deletion and interference. Structural variation-induced enhancer-hijacking and silencer-hijacking events were further identified in AML samples. Hijacked enhancers play a part in AML cell growth, as demonstrated by CRISPR screening, whereas hijacked silencers have a downregulating role, as evidenced by CRISPR-interference-mediated de-repression. Finally, whole-genome bisulfite sequencing of 20 AML and normal samples revealed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Treatment of AML cells with a DNA hypomethylating agent and triple knockdown of DNMT1, DNMT3A and DNMT3B enabled the manipulation of DNA methylation to revert 3D genome organization and gene expression. Overall, this study provides a resource for leukaemia studies and highlights the role of repressive loops and hijacked cis elements in human diseases.

HIC2 controls developmental hemoglobin switching by repressing BCL11A transcription

The fetal-to-adult switch in hemoglobin production is a model of developmental gene control with relevance to the treatment of hemoglobinopathies. The expression of transcription factor BCL11A, which represses fetal β-type globin (HBG) genes in adult erythroid cells, is predominantly controlled at the transcriptional level but the underlying mechanism is unclear. We identify HIC2 as a repressor of BCL11A transcription. HIC2 and BCL11A are reciprocally expressed during development. Forced expression of HIC2 in adult erythroid cells inhibits BCL11A transcription and induces HBG expression. HIC2 binds to erythroid BCL11A enhancers to reduce chromatin accessibility and binding of transcription factor GATA1, diminishing enhancer activity and enhancer-promoter contacts. DNA-binding and crystallography studies reveal direct steric hindrance as one mechanism by which HIC2 inhibits GATA1 binding at a critical BCL11A enhancer. Conversely, loss of HIC2 in fetal erythroblasts increases enhancer accessibility, GATA1 binding and BCL11A transcription. HIC2 emerges as an evolutionarily conserved regulator of hemoglobin switching via developmental control of BCL11A.

Abstract 2953: Subtype-specific and structure variation induced 3D genome alteration in acute myeloid leukemia

Acute myeloid leukemia (AML) represents a set of heterogeneous myeloid malignancies hallmarked by mutations in epigenetic modifiers, transcription factors, and kinases that can cause epigenetic reshaping. It is unclear to what extent AML mutations drive chromatin 3D structure alteration and contribute to myeloid transformation. We first performed Hi-C and whole-genome sequencing in 25 AML patient samples and seven healthy donor samples, and identified recurrent alterations of A/B compartments, TADs, and chromatin loops that are unique to different subtypes. To investigate how altered chromatin organization contributes to transcriptional misregulation, we performed RNA-Seq, ATAC-Seq and CUT&ag for CTCF, H3K27ac, and H3K27me3 in the same AML samples. We identified extensive and recurrent AML-specific promoter-enhancer and promoter-repressor loops. We performed both CRISPR deletion and interference experiments and validated two repressor loops that downregulated cancer related genes IKZF2 and RTTN. Furthermore, by using our recently developed algorithm, we identified structural variation-induced enhancer-hijacking and repressor-hijacking events in AML samples. We further demonstrated the role of hijacked enhancers in AML cell growth by CRISPR screening, and the role of hijacked repressors by CRISPR de-repression. We performed whole-genome bisulfite sequencing in 20 AML and normal samples, and showed the delicate relationship between DNA methylation, CTCF binding and 3D genome structure. Finally, by treating the AML cells with the DNA hypomethylating agent and performing triple knockdown of DNMT1/3A/3B, we demonstrated the impact of altered DNA methylation on gene expression and 3D genome organization. Overall this study provides an invaluable resource for leukemia studies and also highlighted the role of repressor-loops and hijacked cis-elements in gene regulation and human diseases.

Citation Format: Jie Xu, Fan Song, Baozhen Zhang, Huijue Lyu, Mikoto Kobayashi, Ziyu Zhao, Ye Hou, Xiaotao Wang, Yu Luan, Bei Jia, Lena Stasiak, Qixuan Wang, Qi Jin, Qiushi Jin, Yihao Fu, Ross C. Hardison, Sinisa Dovat, Leonidas C. Platanias, Yue Yang, Tomoko Yamada, Aaron D. Viny, Ross L. Levine, David F. Claxton, James R. Broach, Hong Zheng, Feng Yue. Subtype-specific and structure variation induced 3D genome alteration in acute myeloid leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2953.

Dual function NFI factors control fetal hemoglobin silencing in adult erythroid cells

The mechanisms by which the fetal-type β-globin-like genes HBG1 and HBG2 are silenced in adult erythroid precursor cells remain a fundamental question in human biology and have therapeutic relevance to sickle cell disease (SCD) and β-thalassemia. Here, we identify via a CRISPR-Cas9 genetic screen two members of the NFI transcription factor family – NFIA and NFIX – as HBG1/2 repressors. NFIA and NFIX are expressed at elevated levels in adult erythroid cells compared to fetal cells, and function cooperatively to repress HBG1/2 in cultured cells and in human-to-mouse xenotransplants. Genomic profiling, genome editing, and DNA binding assays demonstrate that the potent concerted activity of NFIA and NFIX is explained in part by their ability to stimulate the expression of BCL11A, a known silencer of the HBG1/2 genes, and in part by directly repressing the HBG1/2 genes. Thus, NFI factors emerge as versatile regulators of the fetal-to-adult switch in β-globin production.

Genomic landscape of TCRαβ and TCRγδ T-large granular lymphocyte leukemia

Large granular lymphocyte (LGL) leukemia comprises a group of rare lymphoproliferative disorders whose molecular landscape is incompletely defined. We leveraged paired whole-exome and transcriptome sequencing in the largest LGL leukemia cohort to date, which included 105 patients (93 T-cell receptor αβ [TCRαβ] T-LGL and 12 TCRγδ T-LGL). Seventy-six mutations were observed in 3 or more patients in the cohort, and out of those, STAT3, KMT2D, PIK3R1, TTN, EYS, and SULF1 mutations were shared between both subtypes. We identified ARHGAP25, ABCC9, PCDHA11, SULF1, SLC6A15, DDX59, DNMT3A, FAS, KDM6A, KMT2D, PIK3R1, STAT3, STAT5B, TET2, and TNFAIP3 as recurrently mutated putative drivers using an unbiased driver analysis approach leveraging our whole-exome cohort. Hotspot mutations in STAT3, PIK3R1, and FAS were detected, whereas truncating mutations in epigenetic modifying enzymes such as KMT2D and TET2 were observed. Moreover, STAT3 mutations co-occurred with mutations in chromatin and epigenetic modifying genes, especially KMT2D and SETD1B (P < .01 and P < .05, respectively). STAT3 was mutated in 50.5% of the patients. Most common Y640F STAT3 mutation was associated with lower absolute neutrophil count values, and N647I mutation was associated with lower hemoglobin values. Somatic activating mutations (Q160P, D170Y, L287F) in the STAT3 coiled-coil domain were characterized. STAT3-mutant patients exhibited increased mutational burden and enrichment of a mutational signature associated with increased spontaneous deamination of 5-methylcytosine. Finally, gene expression analysis revealed enrichment of interferon-γ signaling and decreased phosphatidylinositol 3-kinase-Akt signaling for STAT3-mutant patients. These findings highlight the clinical and molecular heterogeneity of this rare disorder.

Domain adaptive neural networks improve cross-species prediction of transcription factor binding

The intrinsic DNA sequence preferences and cell type–specific cooperative partners of transcription factors (TFs) are typically highly conserved. Hence, despite the rapid evolutionary turnover of individual TF binding sites, predictive sequence models of cell type–specific genomic occupancy of a TF in one species should generalize to closely matched cell types in a related species. To assess the viability of cross-species TF binding prediction, we train neural networks to discriminate ChIP-seq peak locations from genomic background and evaluate their performance within and across species. Cross-species predictive performance is consistently worse than within-species performance, which we show is caused in part by species-specific repeats. To account for this domain shift, we use an augmented network architecture to automatically discourage learning of training species–specific sequence features. This domain adaptation approach corrects for prediction errors on species-specific repeats and improves overall cross-species model performance. Our results show that cross-species TF binding prediction is feasible when models account for domain shifts driven by species-specific repeats.

Strategies to improve pharmacogenomic-guided treatment options for patients with β-hemoglobinopathies

Drug efficacy and toxicity are closely related to the unique genetic profile of individuals, or pharmacogenomics. Despite the fact that cardiology, psychiatry and oncology are among the clinical specialties in which pharmacogenomics has become a clinical reality, the utility of pharmacogenomics has yet to be demonstrated for several other medical specialties. Over the last 15 years, genomic variants in a number of loci have been shown to be significantly associated with the fetal hemoglobin (HbF) response to hydroxyurea, the only approved drug for HbF induction for sickle cell disease. Here, we provide an update and discuss future challenges to the application of pharmacogenomics to improve therapies for β-hemoglobinopathies in relation to the current pharmacological treatment modalities for those disorders.

Frequent somatic TET2 mutations in chronic NK-LGL leukemia with distinct patterns of cytopenias

Chronic natural killer large granular lymphocyte (NK-LGL) leukemia, also referred to as chronic lymphoproliferative disorder of NK cells, is a rare disorder defined by prolonged expansion of clonal NK cells. Similar prevalence of STAT3 mutations in chronic T-LGL and NK-LGL leukemia is suggestive of common pathogenesis. We undertook whole-genome sequencing to identify mutations unique to NK-LGL leukemia. The results were analyzed to develop a resequencing panel that was applied to 58 patients. Phosphatidylinositol 3-kinase pathway gene mutations (PIK3CD/PIK3AP1) and TNFAIP3 mutations were seen in 5% and 10% of patients, respectively. TET2 was exceptional in that mutations were present in 16 (28%) of 58 patient samples, with evidence that TET2 mutations can be dominant and exclusive to the NK compartment. Reduced-representation bisulfite sequencing revealed that methylation patterns were significantly altered in TET2 mutant samples. The promoter of TET2 and that of PTPRD, a negative regulator of STAT3, were found to be methylated in additional cohort samples, largely confined to the TET2 mutant group. Mutations in STAT3 were observed in 19 (33%) of 58 patient samples, 7 of which had concurrent TET2 mutations. Thrombocytopenia and resistance to immunosuppressive agents were uniquely observed in those patients with only TET2 mutation (Games-Howell post hoc test, P = .0074; Fisher's exact test, P = .00466). Patients with STAT3 mutation, inclusive of those with TET2 comutation, had lower hematocrit, hemoglobin, and absolute neutrophil count compared with STAT3 wild-type patients (Welch's t test, P ≤ .015). We present the discovery of TET2 mutations in chronic NK-LGL leukemia and evidence that it identifies a unique molecular subtype.

Single-nucleotide-level mapping of DNA regulatory elements that control fetal hemoglobin expression

Pinpointing functional noncoding DNA sequences and defining their contributions to health-related traits is a major challenge for modern genetics. We developed a high-throughput framework to map noncoding DNA functions with single-nucleotide resolution in four loci that control erythroid fetal hemoglobin (HbF) expression, a genetically determined trait that modifies sickle cell disease (SCD) phenotypes. Specifically, we used the adenine base editor ABEmax to introduce 10,156 separate A•T to G•C conversions in 307 predicted regulatory elements and quantified the effects on erythroid HbF expression. We identified numerous regulatory elements, defined their epigenomic structures and linked them to low-frequency variants associated with HbF expression in an SCD cohort. Targeting a newly discovered γ-globin gene repressor element in SCD donor CD34+ hematopoietic progenitors raised HbF levels in the erythroid progeny, inhibiting hypoxia-induced sickling. Our findings reveal previously unappreciated genetic complexities of HbF regulation and provide potentially therapeutic insights into SCD.

Effects of sheared chromatin length on ChIP-seq quality and sensitivity

Chromatin immunoprecipitation followed by massively parallel, high throughput sequencing (ChIP-seq) is the method of choice for genome-wide identification of DNA segments bound by specific transcription factors or in chromatin with particular histone modifications. However, the quality of ChIP-seq datasets varies widely, with a substantial fraction being of intermediate to poor quality. Thus, it is important to discern and control the factors that contribute to variation in ChIP-seq. In this study, we focused on sonication, a user-controlled variable, to produce sheared chromatin. We systematically varied the amount of shearing of fixed chromatin from a mouse erythroid cell line, carefully measuring the distribution of resultant fragment lengths prior to ChIP-seq. This systematic study was complemented with a retrospective analysis of additional experiments. We found that the level of sonication had a pronounced impact on the quality of ChIP-seq signals. Over-sonication consistently reduced quality, while the impact of under-sonication differed among transcription factors, with no impact on sites bound by CTCF but frequently leading to the loss of sites occupied by TAL1 or bound by POL2. The bound sites not observed in low quality datasets were inferred to be a mix of both direct and indirect binding. We leveraged these findings to produce a set of CTCF ChIP-seq datasets in rare, primary hematopoietic progenitor cells. Our observation that the amount of chromatin sonication is a key variable in success of ChIP-seq experiments indicates that monitoring the level of sonication can improve ChIP-seq quality and reproducibility and facilitate ChIP-seq in rare cell types.

S3V2-IDEAS: a package for normalizing, denoising and integrating epigenomic datasets across different cell types

SUMMARY

Epigenetic modifications reflect key aspects of transcriptional regulation, and many epigenomic datasets have been generated under different biological contexts to provide insights into regulatory processes. However, the technical noise in epigenomic datasets and the many dimensions (features) examined make it challenging to effectively extract biologically meaningful inferences from these datasets. We developed a package that reduces noise while normalizing the epigenomic data by a novel normalization method, followed by integrative dimensional reduction by learning and assigning epigenetic states. This package, called S3V2-IDEAS, can be used to identify epigenetic states for multiple features, or identify discretized signal intensity levels and a master peak list across different cell types for a single feature. We illustrate the outputs and performance of S3V2-IDEAS using 137 epigenomics datasets from the VISION project that provides ValIdated Systematic IntegratiON of epigenomic data in hematopoiesis.

AVAILABILITY AND IMPLEMENTATION

S3V2-IDEAS pipeline is freely available as open source software released under an MIT license at: https://github.com/guanjue/S3V2_IDEAS_ESMP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Distinct properties and functions of CTCF revealed by a rapidly inducible degron system

CCCTC-binding factor (CTCF) is a conserved zinc finger transcription factor implicated in a wide range of functions, including genome organization, transcription activation, and elongation. To explore the basis for CTCF functional diversity, we coupled an auxin-induced degron system with precision nuclear run-on. Unexpectedly, oriented CTCF motifs in gene bodies are associated with transcriptional stalling in a manner independent of bound CTCF. Moreover, CTCF at different binding sites (CBSs) displays highly variable resistance to degradation. Motif sequence does not significantly predict degradation behavior, but location at chromatin boundaries and chromatin loop anchors, as well as co-occupancy with cohesin, are associated with delayed degradation. Single-molecule tracking experiments link chromatin residence time to CTCF degradation kinetics, which has ramifications regarding architectural CTCF functions. Our study highlights the heterogeneity of CBSs, uncovers properties specific to architecturally important CBSs, and provides insights into the basic processes of genome organization and transcription regulation.

Clinically relevant updates of the HbVar database of human hemoglobin variants and thalassemia mutations

HbVar (http://globin.bx.psu.edu/hbvar) is a widely-used locus-specific database (LSDB) launched 20 years ago by a multi-center academic effort to provide timely information on the numerous genomic variants leading to hemoglobin variants and all types of thalassemia and hemoglobinopathies. Here, we report several advances for the database. We made clinically relevant updates of HbVar, implemented as additional querying options in the HbVar query page, allowing the user to explore the clinical phenotype of compound heterozygous patients. We also made significant improvements to the HbVar front page, making comparative data querying, analysis and output more user-friendly. We continued to expand and enrich the regular data content, involving 1820 variants, 230 of which are new entries. We also increased the querying potential and expanded the usefulness of HbVar database in the clinical setting. These several additions, expansions and updates should improve the utility of HbVar both for the globin research community and in a clinical setting.

ZNF410 Uniquely Activates the NuRD Component CHD4 to Silence Fetal Hemoglobin Expression

Metazoan transcription factors typically regulate large numbers of genes. Here we identify via a CRISPR-Cas9 genetic screen ZNF410, a pentadactyl DNA-binding protein that in human erythroid cells directly activates only a single gene, the NuRD component CHD4. Specificity is conveyed by two highly evolutionarily conserved clusters of ZNF410 binding sites near the CHD4 gene with no counterparts elsewhere in the genome. Loss of ZNF410 in adult-type human erythroid cell culture systems and xenotransplantation settings diminishes CHD4 levels and derepresses the fetal hemoglobin genes. While previously known to be silenced by CHD4, the fetal globin genes are exposed here as among the most sensitive to reduced CHD4 levels.. In vitro DNA binding assays and crystallographic studies reveal the ZNF410-DNA binding mode. ZNF410 is a remarkably selective transcriptional activator in erythroid cells, and its perturbation might offer new opportunities for treatment of hemoglobinopathies.

A map of cis-regulatory elements and 3D genome structures in zebrafish

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species¹. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Alteration of genome folding via contact domain boundary insertion

Animal chromosomes are partitioned into contact domains. Pathogenic domain disruptions can result from chromosomal rearrangements or perturbation of architectural factors. However, such broad-scale alterations are insufficient to define the minimal requirements for domain formation. Moreover, to what extent domains can be engineered is just beginning to be explored. In an attempt to create contact domains, we inserted a 2-kb DNA sequence underlying a tissue-invariant domain boundary-containing a CTCF-binding site (CBS) and a transcription start site (TSS)-into 16 ectopic loci across 11 chromosomes, and characterized its architectural impact. Depending on local constraints, this fragment variably formed new domains, partitioned existing ones, altered compartmentalization and initiated contacts reflecting chromatin loop extrusion. Deletions of the CBS or the TSS individually or in combination within inserts revealed its distinct contributions to genome folding. Altogether, short DNA insertions can suffice to shape the spatial genome in a manner influenced by chromatin context.

Expanded encyclopaedias of DNA elements in the human and mouse genomes

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE1 and Roadmap Epigenomics2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system1. Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis-element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.

Perspectives on ENCODE

The Encylopedia of DNA Elements (ENCODE) Project launched in 2003 with the long-term goal of developing a comprehensive map of functional elements in the human genome. These included genes, biochemical regions associated with gene regulation (for example, transcription factor binding sites, open chromatin, and histone marks) and transcript isoforms. The marks serve as sites for candidate cis-regulatory elements (cCREs) that may serve functional roles in regulating gene expression1. The project has been extended to model organisms, particularly the mouse. In the third phase of ENCODE, nearly a million and more than 300,000 cCRE annotations have been generated for human and mouse, respectively, and these have provided a valuable resource for the scientific community.

S3norm: simultaneous normalization of sequencing depth and signal-to-noise ratio in epigenomic data

Quantitative comparison of epigenomic data across multiple cell types or experimental conditions is a promising way to understand the biological functions of epigenetic modifications. However, differences in sequencing depth and signal-to-noise ratios in the data from different experiments can hinder our ability to identify real biological variation from raw epigenomic data. Proper normalization is required prior to data analysis to gain meaningful insights. Most existing methods for data normalization standardize signals by rescaling either background regions or peak regions, assuming that the same scale factor is applicable to both background and peak regions. While such methods adjust for differences in sequencing depths, they do not address differences in the signal-to-noise ratios across different experiments. We developed a new data normalization method, called S3norm, that normalizes the sequencing depths and signal-to-noise ratios across different data sets simultaneously by a monotonic nonlinear transformation. We show empirically that the epigenomic data normalized by our method, compared to existing methods, can better capture real biological variation, such as impact on gene expression regulation.

An integrative view of the regulatory and transcriptional landscapes in mouse hematopoiesis

Thousands of epigenomic data sets have been generated in the past decade, but it is difficult for researchers to effectively use all the data relevant to their projects. Systematic integrative analysis can help meet this need, and the VISION project was established for validated systematic integration of epigenomic data in hematopoiesis. Here, we systematically integrated extensive data recording epigenetic features and transcriptomes from many sources, including individual laboratories and consortia, to produce a comprehensive view of the regulatory landscape of differentiating hematopoietic cell types in mouse. By using IDEAS as our integrative and discriminative epigenome annotation system, we identified and assigned epigenetic states simultaneously along chromosomes and across cell types, precisely and comprehensively. Combining nuclease accessibility and epigenetic states produced a set of more than 200,000 candidate cis-regulatory elements (cCREs) that efficiently capture enhancers and promoters. The transitions in epigenetic states of these cCREs across cell types provided insights into mechanisms of regulation, including decreases in numbers of active cCREs during differentiation of most lineages, transitions from poised to active or inactive states, and shifts in nuclease accessibility of CTCF-bound elements. Regression modeling of epigenetic states at cCREs and gene expression produced a versatile resource to improve selection of cCREs potentially regulating target genes. These resources are available from our VISION website to aid research in genomics and hematopoiesis.

Systematic integration of GATA transcription factors and epigenomes via IDEAS paints the regulatory landscape of hematopoietic cells

Members of the GATA family of transcription factors play key roles in the differentiation of specific cell lineages by regulating the expression of target genes. Three GATA factors play distinct roles in hematopoietic differentiation. In order to better understand how these GATA factors function to regulate genes throughout the genome, we are studying the epigenomic and transcriptional landscapes of hematopoietic cells in a model-driven, integrative fashion. We have formed the collaborative multi-lab VISION project to conduct ValIdated Systematic IntegratiON of epigenomic data in mouse and human hematopoiesis. The epigenomic data included nuclease accessibility in chromatin, CTCF occupancy, and histone H3 modifications for 20 cell types covering hematopoietic stem cells, multilineage progenitor cells, and mature cells across the blood cell lineages of mouse. The analysis used the Integrative and Discriminative Epigenome Annotation System (IDEAS), which learns all common combinations of features (epigenetic states) simultaneously in two dimensions-along chromosomes and across cell types. The result is a segmentation that effectively paints the regulatory landscape in readily interpretable views, revealing constitutively active or silent loci as well as the loci specifically induced or repressed in each stage and lineage. Nuclease accessible DNA segments in active chromatin states were designated candidate cis-regulatory elements in each cell type, providing one of the most comprehensive registries of candidate hematopoietic regulatory elements to date. Applications of VISION resources are illustrated for the regulation of genes encoding GATA1, GATA2, GATA3, and Ikaros. VISION resources are freely available from our website http://usevision.org.

OnTAD: hierarchical domain structure reveals the divergence of activity among TADs and boundaries

The spatial organization of chromatin in the nucleus has been implicated in regulating gene expression. Maps of high-frequency interactions between different segments of chromatin have revealed topologically associating domains (TADs), within which most of the regulatory interactions are thought to occur. TADs are not homogeneous structural units but appear to be organized into a hierarchy. We present OnTAD, an optimized nested TAD caller from Hi-C data, to identify hierarchical TADs. OnTAD reveals new biological insights into the role of different TAD levels, boundary usage in gene regulation, the loop extrusion model, and compartmental domains. OnTAD is available at https://github.com/anlin00007/OnTAD.

Chromatin structure dynamics during the mitosis-to-G1 phase transition

Features of higher-order chromatin organization-such as A/B compartments, topologically associating domains and chromatin loops-are temporarily disrupted during mitosis1,2. Because these structures are thought to influence gene regulation, it is important to understand how they are re-established after mitosis. Here we examine the dynamics of chromosome reorganization by Hi-C after mitosis in highly purified, synchronous mouse erythroid cell populations. We observed rapid establishment of A/B compartments, followed by their gradual intensification and expansion. Contact domains form from the 'bottom up'-smaller subTADs are formed initially, followed by convergence into multi-domain TAD structures. CTCF is partially retained on mitotic chromosomes and immediately resumes full binding in ana/telophase. By contrast, cohesin is completely evicted from mitotic chromosomes and regains focal binding at a slower rate. The formation of CTCF/cohesin co-anchored structural loops follows the kinetics of cohesin positioning. Stripe-shaped contact patterns-anchored by CTCF-grow in length, which is consistent with a loop-extrusion process after mitosis. Interactions between cis-regulatory elements can form rapidly, with rates exceeding those of CTCF/cohesin-anchored contacts. Notably, we identified a group of rapidly emerging transient contacts between cis-regulatory elements in ana/telophase that are dissolved upon G1 entry, co-incident with the establishment of inner boundaries or nearby interfering chromatin loops. We also describe the relationship between transcription reactivation and architectural features. Our findings indicate that distinct but mutually influential forces drive post-mitotic chromatin reconfiguration.

Transcriptional Burst Initiation and Polymerase Pause Release Are Key Control Points of Transcriptional Regulation

Transcriptional regulation occurs via changes to rates of different biochemical steps of transcription, but it remains unclear which rates are subject to change upon biological perturbation. Biochemical studies have suggested that stimuli predominantly affect the rates of RNA polymerase II (Pol II) recruitment and polymerase release from promoter-proximal pausing. Single-cell studies revealed that transcription occurs in discontinuous bursts, suggesting that features of such bursts like frequency and intensity could also be regulated. We combined Pol II chromatin immunoprecipitation sequencing (ChIP-seq) and single-cell transcriptional measurements to show that an independently regulated burst initiation step is required before polymerase recruitment can occur. Using a number of global and targeted transcriptional regulatory perturbations, we showed that biological perturbations regulated both burst initiation and polymerase pause release rates but seemed not to regulate polymerase recruitment rate. Our results suggest that transcriptional regulation primarily acts by changing the rates of burst initiation and polymerase pause release.

Integrative detection and analysis of structural variation in cancer genomes

Structural variants (SVs) can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of SVs in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole-genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify SVs in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of noncoding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel three-dimensional chromatin structural domains. Our results indicate that noncoding SVs may be underappreciated mutational drivers in cancer genomes.

Domain-focused CRISPR screen identifies HRI as a fetal hemoglobin regulator in human erythroid cells

Increasing fetal hemoglobin (HbF) levels in adult red blood cells provides clinical benefit to patients with sickle cell disease and some forms of β-thalassemia. To identify potentially druggable HbF regulators in adult human erythroid cells, we employed a protein kinase domain-focused CRISPR-Cas9-based genetic screen with a newly optimized single-guide RNA scaffold. The screen uncovered the heme-regulated inhibitor HRI (also known as EIF2AK1), an erythroid-specific kinase that controls protein translation, as an HbF repressor. HRI depletion markedly increased HbF production in a specific manner and reduced sickling in cultured erythroid cells. Diminished expression of the HbF repressor BCL11A accounted in large part for the effects of HRI depletion. Taken together, these results suggest HRI as a potential therapeutic target for hemoglobinopathies.

Selenoproteins regulate stress erythroid progenitors and spleen microenvironment during stress erythropoiesis

Micronutrient selenium (Se) plays a key role in redox regulation through its incorporation into selenoproteins as the 21st amino acid selenocysteine (Sec). Because Se deficiency appears to be a cofactor in the anemia associated with chronic inflammatory diseases, we reasoned that selenoproteins may contribute to erythropoietic recovery from anemia, referred to as stress erythropoiesis. Here, we report that loss of selenoproteins through Se deficiency or by mutation of the Sec tRNA (tRNA[Sec]) gene (Trsp) severely impairs stress erythropoiesis at 2 stages. Early stress erythroid progenitors failed to expand and properly differentiate into burst-forming unit-erythroid cells , whereas late-stage erythroid progenitors exhibited a maturation defect that affected the transition of proerythroblasts to basophilic erythroblasts. These defects were, in part, a result of the loss of selenoprotein W (SelenoW), whose expression was reduced at both transcript and protein levels in Se-deficient erythroblasts. Mutation of SelenoW in the bone marrow cells significantly decreased the expansion of stress burst-forming unit-erythroid cell colonies, which recapitulated the phenotypes induced by Se deficiency or mutation of Trsp Similarly, mutation of SelenoW in murine erythroblast (G1E) cell line led to defects in terminal differentiation. In addition to the erythroid defects, the spleens of Se-deficient mice contained fewer red pulp macrophages and exhibited impaired development of erythroblastic island macrophages, which make up the niche supporting erythroblast development. Taken together, these data reveal a critical role of selenoproteins in the expansion and development of stress erythroid progenitors, as well as the erythroid niche during acute anemia recovery.

Establishment of regulatory elements during erythro-megakaryopoiesis identifies hematopoietic lineage-commitment points

BACKGROUND

Enhancers and promoters are cis-acting regulatory elements associated with lineage-specific gene expression. Previous studies showed that different categories of active regulatory elements are in regions of open chromatin, and each category is associated with a specific subset of post-translationally marked histones. These regulatory elements are systematically activated and repressed to promote commitment of hematopoietic stem cells along separate differentiation paths, including the closely related erythrocyte (ERY) and megakaryocyte (MK) lineages. However, the order in which these decisions are made remains unclear.

RESULTS

To characterize the order of cell fate decisions during hematopoiesis, we collected primary cells from mouse bone marrow and isolated 10 hematopoietic populations to generate transcriptomes and genome-wide maps of chromatin accessibility and histone H3 acetylated at lysine 27 binding (H3K27ac). Principle component analysis of transcriptional and open chromatin profiles demonstrated that cells of the megakaryocyte lineage group closely with multipotent progenitor populations, whereas erythroid cells form a separate group distinct from other populations. Using H3K27ac and open chromatin profiles, we showed that 89% of immature MK (iMK)-specific active regulatory regions are present in the most primitive hematopoietic cells, 46% of which contain active enhancer marks. These candidate active enhancers are enriched for transcription factor binding site motifs for megakaryopoiesis-essential proteins, including ERG and ETS1. In comparison, only 64% of ERY-specific active regulatory regions are present in the most primitive hematopoietic cells, 20% of which containing active enhancer marks. These regions were not enriched for any transcription factor consensus sequences. Incorporation of genome-wide DNA methylation identified significant levels of de novo methylation in iMK, but not ERY.

CONCLUSIONS

Our results demonstrate that megakaryopoietic profiles are established early in hematopoiesis and are present in the majority of the hematopoietic progenitor population. However, megakaryopoiesis does not constitute a "default" differentiation pathway, as extensive de novo DNA methylation accompanies megakaryopoietic commitment. In contrast, erythropoietic profiles are not established until a later stage of hematopoiesis, and require more dramatic changes to the transcriptional and epigenetic programs. These data provide important insights into lineage commitment and can contribute to ongoing studies related to diseases associated with differentiation defects.

Evolution of hemoglobin loci and their regulatory elements

Across the expanse of vertebrate evolution, each species produces multiple forms of hemoglobin in erythroid cells at appropriate times and in the proper amounts. The multiple hemoglobins are encoded in two globin gene clusters in almost all species. One globin gene cluster, linked to the gene NPRL3, is preserved in all vertebrates, including a gene cluster encoding the highly divergent globins from jawless vertebrates. This preservation of synteny may reflect the presence of a powerful enhancer of globin gene expression in the NPRL3 gene. Despite substantial divergence in noncoding DNA sequences among mammals, several epigenetic features of the globin gene regulatory regions are preserved across vertebrates. The preserved features include multiple DNase hypersensitive sites, at least one of which is an enhancer, and binding by key lineage-restricted transcription factors such as GATA1 and TAL1, which in turn recruit coactivators such as P300 that catalyze acetylation of histones. The maps of epigenetic features are strongly correlated with activity in gene regulation, and resources for accessing and visualizing such maps are readily available to the community of researchers and students.

Accurate and reproducible functional maps in 127 human cell types via 2D genome segmentation

The Roadmap Epigenomics Consortium has published whole-genome functional annotation maps in 127 human cell types by integrating data from studies of multiple epigenetic marks. These maps have been widely used for studying gene regulation in cell type-specific contexts and predicting the functional impact of DNA mutations on disease. Here, we present a new map of functional elements produced by applying a method called IDEAS on the same data. The method has several unique advantages and outperforms existing methods, including that used by the Roadmap Epigenomics Consortium. Using five categories of independent experimental datasets, we compared the IDEAS and Roadmap Epigenomics maps. While the overall concordance between the two maps is high, the maps differ substantially in the prediction details and in their consistency of annotation of a given genomic position across cell types. The annotation from IDEAS is uniformly more accurate than the Roadmap Epigenomics annotation and the improvement is substantial based on several criteria. We further introduce a pipeline that improves the reproducibility of functional annotation maps. Thus, we provide a high-quality map of candidate functional regions across 127 human cell types and compare the quality of different annotation methods in order to facilitate biomedical research in epigenomics.

Between form and function: the complexity of genome folding

It has been known for over a century that chromatin is not randomly distributed within the nucleus. However, the question of how DNA is folded and the influence of such folding on nuclear processes remain topics of intensive current research. A longstanding, unanswered question is whether nuclear organization is simply a reflection of nuclear processes such as transcription and replication, or whether chromatin is folded by independent mechanisms and this per se encodes function? Evidence is emerging that both may be true. Here, using the α-globin gene cluster as an illustrative model, we provide an overview of the most recent insights into the layers of genome organization across different scales and how this relates to gene activity.

Comparative analysis of three-dimensional chromosomal architecture identifies a novel fetal hemoglobin regulatory element

In this study, Huang et al. compared the chromosomal architectures of fetal and adult human erythroblasts and found that, globally, chromatin structures and compartments A/B are highly similar at both developmental stages. Their results uncover a new critical regulatory region as a potential target for therapeutic genome editing for hemoglobinopathies and highlight the power of chromosome conformation analysis in discovering new cis control elements.

Chromatin structure is tightly intertwined with transcription regulation. Here we compared the chromosomal architectures of fetal and adult human erythroblasts and found that, globally, chromatin structures and compartments A/B are highly similar at both developmental stages. At a finer scale, we detected distinct folding patterns at the developmentally controlled β-globin locus. Specifically, new fetal stage-specific contacts were uncovered between a region separating the fetal (γ) and adult (δ and β) globin genes (encompassing the HBBP1 and BGLT3 noncoding genes) and two distal chromosomal sites (HS5 and 3′HS1) that flank the locus. In contrast, in adult cells, the HBBP1–BGLT3 region contacts the embryonic ε-globin gene, physically separating the fetal globin genes from the enhancer (locus control region [LCR]). Deletion of the HBBP1 region in adult cells alters contact landscapes in ways more closely resembling those of fetal cells, including increased LCR–γ-globin contacts. These changes are accompanied by strong increases in γ-globin transcription. Notably, the effects of HBBP1 removal on chromatin architecture and gene expression closely mimic those of deleting the fetal globin repressor BCL11A, implicating BCL11A in the function of the HBBP1 region. Our results uncover a new critical regulatory region as a potential target for therapeutic genome editing for hemoglobinopathies and highlight the power of chromosome conformation analysis in discovering new cis control elements.

HiCRep: assessing the reproducibility of Hi-C data using a stratum-adjusted correlation coefficient

Hi-C is a powerful technology for studying genome-wide chromatin interactions. However, current methods for assessing Hi-C data reproducibility can produce misleading results because they ignore spatial features in Hi-C data, such as domain structure and distance dependence. We present HiCRep, a framework for assessing the reproducibility of Hi-C data that systematically accounts for these features. In particular, we introduce a novel similarity measure, the stratum adjusted correlation coefficient (SCC), for quantifying the similarity between Hi-C interaction matrices. Not only does it provide a statistically sound and reliable evaluation of reproducibility, SCC can also be used to quantify differences between Hi-C contact matrices and to determine the optimal sequencing depth for a desired resolution. The measure consistently shows higher accuracy than existing approaches in distinguishing subtle differences in reproducibility and depicting interrelationships of cell lineages. The proposed measure is straightforward to interpret and easy to compute, making it well-suited for providing standardized, interpretable, automatable, and scalable quality control. The freely available R package HiCRep implements our approach.

A hyperactive transcriptional state marks genome reactivation at the mitosis-G1 transition

Hsiung et al. tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. During the earliest rounds of transcription at the mitosis–G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression.

During mitosis, RNA polymerase II (Pol II) and many transcription factors dissociate from chromatin, and transcription ceases globally. Transcription is known to restart in bulk by telophase, but whether de novo transcription at the mitosis–G1 transition is in any way distinct from later in interphase remains unknown. We tracked Pol II occupancy genome-wide in mammalian cells progressing from mitosis through late G1. Unexpectedly, during the earliest rounds of transcription at the mitosis–G1 transition, ∼50% of active genes and distal enhancers exhibit a spike in transcription, exceeding levels observed later in G1 phase. Enhancer–promoter chromatin contacts are depleted during mitosis and restored rapidly upon G1 entry but do not spike. Of the chromatin-associated features examined, histone H3 Lys27 acetylation levels at individual loci in mitosis best predict the mitosis–G1 transcriptional spike. Single-molecule RNA imaging supports that the mitosis–G1 transcriptional spike can constitute the maximum transcriptional activity per DNA copy throughout the cell division cycle. The transcriptional spike occurs heterogeneously and propagates to cell-to-cell differences in mature mRNA expression. Our results raise the possibility that passage through the mitosis–G1 transition might predispose cells to diverge in gene expression states.

A guide to translation of research results from model organisms to human

A new study helps resolve a controversy about determinants of gene expression variability and might facilitate the effective translation of research results across species.

Jointly characterizing epigenetic dynamics across multiple human cell types

Advanced sequencing technologies have generated a plethora of data for many chromatin marks in multiple tissues and cell types, yet there is lack of a generalized tool for optimal utility of those data. A major challenge is to quantitatively model the epigenetic dynamics across both the genome and many cell types for understanding their impacts on differential gene regulation and disease. We introduce IDEAS, an integrative and discriminative epigenome annotation system, for jointly characterizing epigenetic landscapes in many cell types and detecting differential regulatory regions. A key distinction between our method and existing state-of-the-art algorithms is that IDEAS integrates epigenomes of many cell types simultaneously in a way that preserves the position-dependent and cell type-specific information at fine scales, thereby greatly improving segmentation accuracy and producing comparable annotations across cell types.

Unlinking an lncRNA from Its Associated cis Element

Long non-coding (lnc) RNAs can regulate gene expression and protein functions. However, the proportion of lncRNAs with biological activities among the thousands expressed in mammalian cells is controversial. We studied Lockd (lncRNA downstream of Cdkn1b), a 434-nt polyadenylated lncRNA originating 4 kb 3' to the Cdkn1b gene. Deletion of the 25-kb Lockd locus reduced Cdkn1b transcription by approximately 70% in an erythroid cell line. In contrast, homozygous insertion of a polyadenylation cassette 80 bp downstream of the Lockd transcription start site reduced the entire lncRNA transcript level by >90% with no effect on Cdkn1b transcription. The Lockd promoter contains a DNase-hypersensitive site, binds numerous transcription factors, and physically associates with the Cdkn1b promoter in chromosomal conformation capture studies. Therefore, the Lockd gene positively regulates Cdkn1b transcription through an enhancer-like cis element, whereas the lncRNA itself is dispensable, which may be the case for other lncRNAs.