Chromatin architecture reorganization during stem cell differentiation

Systematically characterizing dysfunctional long intergenic non-coding RNAs in multiple brain regions of major psychosis

Schizophrenia (SZ) and bipolar disorder (BD) are severe neuropsychiatric disorders with serious impact on patients, together termed “major psychosis”. Recently, long intergenic non-coding RNAs (lincRNAs) were reported to play important roles in mental diseases. However, little was known about their molecular mechanism in pathogenesis of SZ and BD. Here, we performed RNA sequencing on 82 post-mortem brain tissues from three brain regions (orbitofrontal cortex (BA11), anterior cingulate cortex (BA24) and dorsolateral prefrontal cortex (BA9)) of patients with SZ and BD and control subjects, generating over one billion reads. We characterized lincRNA transcriptome in the three brain regions and identified 20 differentially expressed lincRNAs (DELincRNAs) in BA11 for BD, 34 and 1 in BA24 and BA9 for SZ, respectively. Our results showed that these DELincRNAs exhibited brain region-specific patterns. Applying weighted gene co-expression network analysis, we revealed that DELincRNAs together with other genes can function as modules to perform different functions in different brain regions, such as immune system development in BA24 and oligodendrocyte differentiation in BA9. Additionally, we found that DNA methylation alteration could partly explain the dysregulation of lincRNAs, some of which could function as enhancers in the pathogenesis of major psychosis. Together, we performed systematical characterization of dysfunctional lincRNAs in multiple brain regions of major psychosis, which provided a valuable resource to understand their roles in SZ and BD pathology and helped to discover novel biomarkers.

Combined Ectopic Expression of Homologous Recombination Factors Promotes Embryonic Stem Cell Differentiation

Homologous recombination (HR), which ensures accurate DNA replication and strand-break repair, is necessary to preserve embryonic stem cell (ESC) self-renewal. However, little is known about how HR factors modulate ESC differentiation and replication stress-associated DNA breaks caused by unique cell-cycle progression. Here, we report that ESCs utilize Rad51-dependent HR to enhance viability and induce rapid proliferation through a replication-coupled pathway. In addition, ESC differentiation was shown to be enhanced by ectopic expression of a subset of recombinases. Abundant expression of HR proteins throughout the ESC cycle, but not during differentiation, facilitated immediate HR-mediated repair of single-stranded DNA (ssDNA) gaps incurred during S-phase, via a mechanism that does not perturb cellular progression. Intriguingly, combined ectopic expression of two recombinases, Rad51 and Rad52, resulted in efficient ESC differentiation and diminished cell death, indicating that HR factors promote cellular differentiation by repairing global DNA breaks induced by chromatin remodeling signals. Collectively, these findings provide insight into the role of key HR factors in rapid DNA break repair following chromosome duplication during self-renewal and differentiation of ESCs.

Alterations of specific chromatin conformation affect ATRA-induced leukemia cell differentiation

Chromatin conformation plays a key role in regulating gene expression and controlling cell differentiation. However, the whole-genome chromatin conformation changes that occur during leukemia cell differentiation are poorly understood. Here, we characterized the changes in chromatin conformation, histone states, chromatin accessibility, and gene expression using an all-trans retinoic acid (ATRA)-induced HL-60 cell differentiation model. The results showed that the boundaries of topological associated domains (TADs) were stable during differentiation; however, the chromatin conformations within several specific TADs were obviously changed. By combining H3K4me3, H3K27ac, and Hi-C signals, we annotated the differential gene-regulatory chromatin interactions upon ATRA induction. The gains and losses of the gene-regulatory chromatin interactions are significantly correlated with gene expression and chromatin accessibility. Finally, we found that the loss of GATA2 expression and DNA binding are crucial for the differentiation process, and changes in the chromatin structure around the GATA2 regulate its expression upon ATRA induction. This study provided both statistical insights and experimental details regarding the relationship between chromatin conformation changes and transcription regulation during leukemia cell differentiation, and the results suggested that the chromatin conformation is a new type of potential drug target for cancer therapy.

Stratification of TAD boundaries reveals preferential insulation of super-enhancers by strong boundaries

The metazoan genome is compartmentalized in areas of highly interacting chromatin known as topologically associating domains (TADs). TADs are demarcated by boundaries mostly conserved across cell types and even across species. However, a genome-wide characterization of TAD boundary strength in mammals is still lacking. In this study, we first use fused two-dimensional lasso as a machine learning method to improve Hi-C contact matrix reproducibility, and, subsequently, we categorize TAD boundaries based on their insulation score. We demonstrate that higher TAD boundary insulation scores are associated with elevated CTCF levels and that they may differ across cell types. Intriguingly, we observe that super-enhancers are preferentially insulated by strong boundaries. Furthermore, we demonstrate that strong TAD boundaries and super-enhancer elements are frequently co-duplicated in cancer patients. Taken together, our findings suggest that super-enhancers insulated by strong TAD boundaries may be exploited, as a functional unit, by cancer cells to promote oncogenesis.

Topologically associating domains (TADs) detected by Hi-C technologies are megabase-scale areas of highly interacting chromatin. Here Gong, Lazaris et al. develop a computational approach to improve the reproducibility of Hi-C contact matrices and stratify TAD boundaries based on their insulating strength.

Sequence-based multiscale modeling for high-throughput chromosome conformation capture (Hi-C) data analysis

In this paper, we introduce sequence-based multiscale modeling for biomolecular data analysis. We employ spectral clustering method in our modeling and reveal the difference between sequence-based global scale clustering and local scale clustering. Essentially, two types of distances, i.e., Euclidean (or spatial) distance and genomic (or sequential) distance, can be used in data clustering. Clusters from sequence-based global scale models optimize spatial distances, meaning spatially adjacent loci are more likely to be assigned into the same cluster. Sequence-based local scale models, on the other hand, result in clusters that optimize genomic distances. That is to say, in these models, sequentially adjoining loci tend to be cluster together. We propose two sequence-based multiscale models (SeqMMs) for the study of chromosome hierarchical structures, including genomic compartments and topological associated domains (TADs). We find that genomic compartments are determined only by global scale information in the Hi-C data. The removal of all the local interactions within a band region as large as 10 Mb in genomic distance has almost no significant influence on the final compartment results. Further, in TAD analysis, we find that when the sequential scale is small, a tiny variation of diagonal band region in a contact map will result in a great change in the predicted TAD boundaries. When the scale value is larger than a threshold value, the TAD boundaries become very consistent. This threshold value is highly related to TAD sizes. By the comparison of our results with those previously obtained using a spectral clustering model, we find that our method is more robust and reliable. Finally, we demonstrate that almost all TAD boundaries from both clustering methods are local minimum of a TAD summation function.

Transcription factors orchestrate dynamic interplay between genome topology and gene regulation during cell reprogramming

Chromosomal architecture is known to influence gene expression, yet its role in controlling cell fate remains poorly understood. Reprogramming of somatic cells into pluripotent stem cells (PSCs) by the transcription factors (TFs) OCT4, SOX2, KLF4 and MYC offers an opportunity to address this question but is severely limited by the low proportion of responding cells. We have recently developed a highly efficient reprogramming protocol that synchronously converts somatic into pluripotent stem cells. Here, we used this system to integrate time-resolved changes in genome topology with gene expression, TF binding and chromatin-state dynamics. The results showed that TFs drive topological genome reorganization at multiple architectural levels, often before changes in gene expression. Removal of locus-specific topological barriers can explain why pluripotency genes are activated sequentially, instead of simultaneously, during reprogramming. Together, our results implicate genome topology as an instructive force for implementing transcriptional programs and cell fate in mammals.

Evolutionary dynamics of 3D genome architecture following polyploidization in cotton

The formation of polyploids significantly increases the complexity of transcriptional regulation, which is expected to be reflected in sophisticated higher-order chromatin structures. However, knowledge of three-dimensional (3D) genome structure and its dynamics during polyploidization remains poor. Here, we characterize 3D genome architectures for diploid and tetraploid cotton, and find the existence of A/B compartments and topologically associated domains (TADs). By comparing each subgenome in tetraploids with its extant diploid progenitor, we find that genome allopolyploidization has contributed to the switching of A/B compartments and the reorganization of TADs in both subgenomes. We also show that the formation of TAD boundaries during polyploidization preferentially occurs in open chromatin, coinciding with the deposition of active chromatin modification. Furthermore, analysis of inter-subgenomic chromatin interactions has revealed the spatial proximity of homoeologous genes, possibly associated with their coordinated expression. This study advances our understanding of chromatin organization in plants and sheds new light on the relationship between 3D genome evolution and transcriptional regulation.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Histone H3 lysine 4 monomethylation modulates long-range chromatin interactions at enhancers

Long-range chromatin interactions between enhancers and promoters are essential for transcription of many developmentally controlled genes in mammals and other metazoans. Currently, the exact mechanisms that connect distal enhancers to their specific target promoters remain to be fully elucidated. Here, we show that the enhancer-specific histone H3 lysine 4 monomethylation (H3K4me1) and the histone methyltransferases MLL3 and MLL4 (MLL3/4) play an active role in this process. We demonstrate that in differentiating mouse embryonic stem cells, MLL3/4-dependent deposition of H3K4me1 at enhancers correlates with increased levels of chromatin interactions, whereas loss of this histone modification leads to reduced levels of chromatin interactions and defects in gene activation during differentiation. H3K4me1 facilitates recruitment of the Cohesin complex, a known regulator of chromatin organization, to chromatin in vitro and in vivo, providing a potential mechanism for MLL3/4 to promote chromatin interactions between enhancers and promoters. Taken together, our results support a role for MLL3/4-dependent H3K4me1 in orchestrating long-range chromatin interactions at enhancers in mammalian cells.

The Role of Indoleamine-2,3-Dioxygenase in Cancer Development, Diagnostics, and Therapy

Tumors are composed of abnormally transformed cell types and tissues that differ from normal tissues in their genetic and epigenetic makeup, metabolism, and immunology. Molecular compounds that modulate the immune response against neoplasms offer promising new strategies to combat cancer. Inhibitors targeting the indoleamine-2,3-dioxygenase 1 enzyme (IDO1) represent one of the most potent therapeutic opportunities to inhibit tumor growth. Herein, we assess the biochemical role of IDO1 in tumor metabolism and immune surveillance, and review current diagnostic and therapeutic approaches that are intended to increase the effectiveness of immunotherapies against highly aggressive and difficult-to-treat IDO-expressing cancers.

Genome-Wide Association Studies in Glioma

Since the first reports in 2009, genome-wide association studies (GWAS) have been successful in identifying germline variants associated with glioma susceptibility. In this review, we describe a chronological history of glioma GWAS, culminating in the most recent study comprising 12,496 cases and 18,190 controls. We additionally summarize associations at the 27 glioma-risk SNPs that have been reported so far. Future efforts are likely to be principally focused on assessing association of germline-risk SNPs with particular molecular subgroups of glioma, as well as investigating the functional basis of the risk loci in tumor formation. These ongoing studies will be important to maximize the impact of research into glioma susceptibility, both in terms of insight into tumor etiology as well as opportunities for clinical translation. Cancer Epidemiol Biomarkers Prev; 27(4); 418-28. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Three-dimensional genome architecture in health and disease

More than a decade of massive DNA sequencing efforts have generated a large body of genomic, transcriptomic and epigenomic information that has provided a more and more detailed view of the functional elements and transactions within the human genome. Considerable efforts have also focused on linking these elements with one another by mapping their interactions and by establishing 3-dimensional (3D) genomic landscapes in various cell and tissue types. In parallel, multiple studies have associated genomic deletions, duplications and other rearrangements with human pathologies. In this review, we explore recent progresses that have allowed connecting disease-causing alterations with perturbations of the 3D genome organization.

Genetic and epigenetic factors which modulate differentiation propensity in human pluripotent stem cells

BACKGROUND

Human pluripotent stem cell (hPSC) lines are known to have a bias in their differentiation. This gives individual cell lines a propensity to preferentially differentiate towards one germ layer or cell type over others. Chromosomal aberrations, mitochondrial mutations, genetic diversity and epigenetic variance are the main drivers of this phenomenon, and can lead to a wide range of phenotypes.

OBJECTIVE AND RATIONALE

Our aim is to provide a comprehensive overview of the different factors which influence differentiation propensity. Specifically, we sought to highlight known genetic variances and their mechanisms, in addition to more general observations from larger abnormalities. Furthermore, we wanted to provide an up-to-date list of a growing number of predictive indicators which are able to identify differentiation propensity before the initiation of differentiation. As differentiation propensity can lead to difficulties in both research as well as clinical translation, our thorough overview could be a useful tool.

SEARCH METHODS

Combinations of the following key words were applied as search criteria in the PubMed database: embryonic stem cells, induced pluripotent stem cells, differentiation propensity (also: potential, efficiency, capacity, bias, variability), epigenetics, chromosomal abnormalities, genetic aberrations, X chromosome inactivation, mitochondrial function, mitochondrial metabolism, genetic diversity, reprogramming, predictive marker, residual stem cell, clinic. Only studies in English were included, ranging from 2000 to 2017, with a majority ranging from 2010 to 1017. Further manuscripts were added from cross-references.

OUTCOMES

Differentiation propensity is affected by a wide variety of (epi)genetic factors. These factors clearly lead to a loss of differentiation capacity, preference towards certain cell types and oftentimes, phenotypes which begin to resemble cancer. Broad changes in (epi)genetics, such as aneuploidies or wide-ranging modifications to the epigenetic landscape tend to lead to extensive, less definite changes in differentiation capacity, whereas more specific abnormalities often have precise ramifications in which certain cell types become more preferential. Furthermore, there appears to be a greater, though often less considered, contribution to differentiation propensity by factors such as mitochondria and inherent genetic diversity. Varied differentiation capacity can also lead to potential consequences in the clinical translation of hPSC, including the occurrence of residual undifferentiated stem cells, and the transplantation of potentially transformed cells.

WIDER IMPLICATIONS

As hPSC continue to advance towards the clinic, our understanding of them progresses as well. As a result, the challenges faced become more numerous, but also more clear. If the transition to the clinic is to be achieved with a minimum number of potential setbacks, thorough evaluation of the cells will be an absolute necessity. Altered differentiation propensity represents at least one such hurdle, for which researchers and eventually clinicians will need to find solutions. Already, steps are being taken to tackle the issue, though further research will be required to evaluate any long-term risks it poses.

Identifying co-opted transposable elements using comparative epigenomics

The human genome gives rise to different epigenomic landscapes that define each cell type and can be deregulated in disease. Recent efforts by ENCODE, the NIH Roadmap and the International Human Epigenome Consortium (IHEC) have made significant advances towards assembling reference epigenomic maps of various tissues. Notably, these projects have found that approximately 80% of human DNA was biochemically active in at least one epigenomic assay while only approximately 10% of the sequence displayed signs of purifying selection. Given that transposable elements (TEs) make up at least 50% of the human genome and can be actively transcribed or act as regulatory elements either for their own purposes or be co-opted for the benefit of their host; we are interested in exploring their overall contribution to the "functional" genome. Traditional methods used to identify functional DNA have relied on comparative genomics, conservation analysis and low throughput validation assays. To discover co-opted TEs, and distinguish them from noisy genomic elements, we argue that comparative epigenomic methods will also be important.

Signaling by Steroid Hormones in the 3D Nuclear Space

Initial studies showed that ligand-activated hormone receptors act by binding to the proximal promoters of individual target genes. Genome-wide studies have now revealed that regulation of transcription by steroid hormones mainly depends on binding of the receptors to distal regulatory elements. Those distal elements, either enhancers or silencers, act on the regulation of target genes by chromatin looping to the gene promoters. In the nucleus, this level of chromatin folding is integrated within dynamic higher orders of genome structures, which are organized in a non-random fashion. Terminally differentiated cells exhibit a tissue-specific three-dimensional (3D) organization of the genome that favors or restrains the activity of transcription factors and modulates the function of steroid hormone receptors, which are transiently activated upon hormone exposure. Conversely, integration of the hormones signal may require modifications of the 3D organization to allow appropriate transcriptional outcomes. In this review, we summarize the main levels of organization of the genome, review how they can modulate the response to steroids in a cell specific manner and discuss the role of receptors in shaping and rewiring the structure in response to hormone. Taking into account the dynamics of 3D genome organization will contribute to a better understanding of the pleiotropic effects of steroid hormones in normal and cancer cells.

The Epstein-Barr Virus Episome Maneuvers between Nuclear Chromatin Compartments during Reactivation

The human genome is structurally organized in three-dimensional space to facilitate functional partitioning of transcription. We learned that the latent episome of the human Epstein-Barr virus (EBV) preferentially associates with gene-poor chromosomes and avoids gene-rich chromosomes. Kaposi's sarcoma-associated herpesvirus behaves similarly, but human papillomavirus does not. Contacts on the EBV side localize to OriP, the latent origin of replication. This genetic element and the EBNA1 protein that binds there are sufficient to reconstitute chromosome association preferences of the entire episome. Contacts on the human side localize to gene-poor and AT-rich regions of chromatin distant from transcription start sites. Upon reactivation from latency, however, the episome moves away from repressive heterochromatin and toward active euchromatin. Our work adds three-dimensional relocalization to the molecular events that occur during reactivation. Involvement of myriad interchromosomal associations also suggests a role for this type of long-range association in gene regulation.

IMPORTANCE The human genome is structurally organized in three-dimensional space, and this structure functionally affects transcriptional activity. We set out to investigate whether a double-stranded DNA virus, Epstein-Barr virus (EBV), uses mechanisms similar to those of the human genome to regulate transcription. We found that the EBV genome associates with repressive compartments of the nucleus during latency and with active compartments during reactivation. This study advances our knowledge of the EBV life cycle, adding three-dimensional relocalization as a novel component to the molecular events that occur during reactivation. Furthermore, the data add to our understanding of nuclear compartments, showing that disperse interchromosomal interactions may be important for regulating transcription.

Sub-kb Hi-C in D. melanogaster reveals conserved characteristics of TADs between insect and mammalian cells

Topologically associating domains (TADs) are fundamental elements of the eukaryotic genomic structure. However, recent studies suggest that the insulating complexes, CTCF/cohesin, present at TAD borders in mammals are absent from those in Drosophila melanogaster, raising the possibility that border elements are not conserved among metazoans. Using in situ Hi-C with sub-kb resolution, here we show that the D. melanogaster genome is almost completely partitioned into >4000 TADs, nearly sevenfold more than previously identified. The overwhelming majority of these TADs are demarcated by the insulator complexes, BEAF-32/CP190, or BEAF-32/Chromator, indicating that these proteins may play an analogous role in flies as that of CTCF/cohesin in mammals. Moreover, extended regions previously thought to be unstructured are shown to consist of small contiguous TADs, a property also observed in mammals upon re-examination. Altogether, our work demonstrates that fundamental features associated with the higher-order folding of the genome are conserved from insects to mammals.

Conserved non-coding elements: developmental gene regulation meets genome organization

Comparative genomics has revealed a class of non-protein-coding genomic sequences that display an extraordinary degree of conservation between two or more organisms, regularly exceeding that found within protein-coding exons. These elements, collectively referred to as conserved non-coding elements (CNEs), are non-randomly distributed across chromosomes and tend to cluster in the vicinity of genes with regulatory roles in multicellular development and differentiation. CNEs are organized into functional ensembles called genomic regulatory blocks–dense clusters of elements that collectively coordinate the expression of shared target genes, and whose span in many cases coincides with topologically associated domains. CNEs display sequence properties that set them apart from other sequences under constraint, and have recently been proposed as useful markers for the reconstruction of the evolutionary history of organisms. Disruption of several of these elements is known to contribute to diseases linked with development, and cancer. The emergence, evolutionary dynamics and functions of CNEs still remain poorly understood, and new approaches are required to enable comprehensive CNE identification and characterization. Here, we review current knowledge and identify challenges that need to be tackled to resolve the impasse in understanding extreme non-coding conservation.

Phenotypic interpretation of complex chromosomal rearrangements informed by nucleotide-level resolution and structural organization of chromatin

Molecular characterization of balanced chromosomal abnormalities constitutes a powerful tool in understanding the pathogenic mechanisms of complex genetic disorders. Here we report a male with severe global developmental delay in the presence of a complex karyotype and normal microarray and exome studies. The subject, referred to as DGAP294, has two de novo apparently balanced translocations involving chromosomes 1 and 14, and chromosomes 4 and 10, disrupting several different transcripts of adhesion G protein-coupled receptor L2 (ADGRL2) and protocadherin 15 (PCDH15). In addition, a maternally inherited inversion disrupts peptidyl arginine deiminase types 3 and 4 (PADI3 and PADI4) on chromosome 1. None of these gene disruptions explain the patient's phenotype. Using genome regulatory annotations and chromosome conformation data, we predict a position effect ~370 kb upstream of a translocation breakpoint located at 14q12. The position effect involves forkhead box G1 (FOXG1), mutations in which are associated with the congenital form of Rett syndrome and FOXG1 syndrome. We believe the FOXG1 position effect largely accounts for the clinical phenotype in DGAP294, which can be classified as FOXG1 syndrome like. Our findings emphasize the significance of not only analyzing disrupted genes by chromosomal rearrangements, but also evaluating potential long-range position effects in clinical diagnoses.

Shaping epigenetic memory via genomic bookmarking

Reconciling the stability of epigenetic patterns with the rapid turnover of histone modifications and their adaptability to external stimuli is an outstanding challenge. Here, we propose a new biophysical mechanism that can establish and maintain robust yet plastic epigenetic domains via genomic bookmarking (GBM). We model chromatin as a recolourable polymer whose segments bear non-permanent histone marks (or colours) which can be modified by 'writer' proteins. The three-dimensional chromatin organisation is mediated by protein bridges, or 'readers', such as Polycomb Repressive Complexes and Transcription Factors. The coupling between readers and writers drives spreading of biochemical marks and sustains the memory of local chromatin states across replication and mitosis. In contrast, GBM-targeted perturbations destabilise the epigenetic patterns. Strikingly, we demonstrate that GBM alone can explain the full distribution of Polycomb marks in a whole Drosophila chromosome. We finally suggest that our model provides a starting point for an understanding of the biophysics of cellular differentiation and reprogramming.

Spatiotemporal allele organization by allele-specific CRISPR live-cell imaging (SNP-CLING)

Imaging and chromatin capture techniques have shed important insights into our understanding of nuclear organization. A limitation of these techniques is the inability to resolve allele-specific spatiotemporal properties of genomic loci in living cells. Here, we describe an allele-specific CRISPR live-cell DNA imaging technique (SNP-CLING) to provide the first comprehensive insights into allelic positioning across space and time in mouse embryonic stem cells and fibroblasts. In 3D-imaging, we studied alleles on different chromosomes in relation to one another and relative to nuclear substructures such as the nucleolus. We find that alleles maintain similar positions relative to each other and the nucleolus, however loci occupy different unique positions. To monitor spatiotemporal dynamics by SNP-CLING, we performed 4D-imaging, determining that alleles are either stably positioned, or fluctuating during cell state transitions, such as apoptosis. SNP-CLING is a universally applicable technique that enables dissecting allele-specific spatiotemporal genome organization in live cells.

β-Actin-dependent global chromatin organization and gene expression programs control cellular identity

During differentiation and development, cell fate and identity are established by waves of genetic reprogramming. Although the mechanisms are largely unknown, during these events, dynamic chromatin reorganization is likely to ensure that multiple genes involved in the same cellular functions are coregulated, depending on the nuclear environment. In this study, using high-content screening of embryonic fibroblasts from a β-actin knockout (KO) mouse, we found major chromatin rearrangements and changes in histone modifications, such as methylated histone (H)3-lysine-(K)9. Genome-wide H3K9 trimethylation-(Me)3 landscape changes correlate with gene up- and down-regulation in β-actin KO cells. Mechanistically, we found loss of chromatin association by the Brahma-related gene ( Brg)/Brahma-associated factor (BAF) chromatin remodeling complex subunit Brg1 in the absence of β-actin. This actin-dependent chromatin reorganization was concomitant with the up-regulation of sets of genes involved in angiogenesis, cytoskeletal organization, and myofibroblast features in β-actin KO cells. Some of these genes and phenotypes were gained in a β-actin dose-dependent manner. Moreover, reintroducing a nuclear localization signal-containing β-actin in the knockout cells affected nuclear features and gene expression. Our results suggest that, by affecting the genome-wide organization of heterochromatin through the chromatin-binding activity of the BAF complex, β-actin plays an essential role in the determination of gene expression programs and cellular identity.-Xie, X., Almuzzaini, B., Drou, N., Kremb, S., Yousif, A., Östlund Farrants, A.-K., Gunsalus, K., Percipalle, P. β-Actin-dependent global chromatin organization and gene expression programs control cellular identity.

Somatic Superenhancer Duplications and Hotspot Mutations Lead to Oncogenic Activation of the KLF5 Transcription Factor

The Krüppel-like family of transcription factors plays critical roles in human development and is associated with cancer pathogenesis. Krüppel-like factor 5 gene (KLF5) has been shown to promote cancer cell proliferation and tumorigenesis and to be genomically amplified in cancer cells. We recently reported that the KLF5 gene is also subject to other types of somatic coding and noncoding genomic alterations in diverse cancer types. Here, we show that these alterations activate KLF5 by three distinct mechanisms: (i) Focal amplification of superenhancers activates KLF5 expression in squamous cell carcinomas; (ii) Missense mutations disrupt KLF5-FBXW7 interactions to increase KLF5 protein stability in colorectal cancer; (iii) Cancer type-specific hotspot mutations within a zinc-finger DNA binding domain of KLF5 change its DNA binding specificity and reshape cellular transcription. Utilizing data from CRISPR/Cas9 gene knockout screening, we reveal that cancer cells with KLF5 overexpression are dependent on KLF5 for their proliferation, suggesting KLF5 as a putative therapeutic target.Significance: Our observations, together with previous studies that identified oncogenic properties of KLF5, establish the importance of KLF5 activation in human cancers, delineate the varied genomic mechanisms underlying this occurrence, and nominate KLF5 as a putative target for therapeutic intervention in cancer. Cancer Discov; 8(1); 108-25. ©2017 AACR.This article is highlighted in the In This Issue feature, p. 1.

Molecular classification of Crohn's disease reveals two clinically relevant subtypes

OBJECTIVE

The clinical presentation and course of Crohn's disease (CD) is highly variable. We sought to better understand the cellular and molecular mechanisms that guide this heterogeneity, and characterise the cellular processes associated with disease phenotypes.

DESIGN

We examined both gene expression and gene regulation (chromatin accessibility) in non-inflamed colon tissue from a cohort of adult patients with CD and control patients. To support the generality of our findings, we analysed previously published expression data from a large cohort of treatment-naïve paediatric CD and control ileum.

RESULTS

We found that adult patients with CD clearly segregated into two classes based on colon tissue gene expression-one that largely resembled the normal colon and one where certain genes showed expression patterns normally specific to the ileum. These classes were supported by changes in gene regulatory profiles observed at the level of chromatin accessibility, reflective of a fundamental shift in underlying molecular phenotypes. Furthermore, gene expression from the ilea of a treatment-naïve cohort of paediatric patients with CD could be similarly subdivided into colon-like and ileum-like classes. Finally, expression patterns within these CD subclasses highlight large-scale differences in the immune response and aspects of cellular metabolism, and were associated with multiple clinical phenotypes describing disease behaviour, including rectal disease and need for colectomy.

CONCLUSIONS

Our results strongly suggest that these molecular signatures define two clinically relevant forms of CD irrespective of tissue sampling location, patient age or treatment status.

Practical Analysis of Hi-C Data: Generating A/B Compartment Profiles

Recent advances in next-generation sequencing (NGS) and chromosome conformation capture (3C) analysis have led to the development of Hi-C, a genome-wide version of the 3C method. Hi-C has identified new levels of chromosome organization such as A/B compartments, topologically associating domains (TADs) as well as large megadomains on the inactive X chromosome, while allowing the identification of chromatin loops at the genome scale. Despite its powerfulness, Hi-C data analysis is much more involved compared to conventional NGS applications such as RNA-seq or ChIP-seq and requires many more steps. This presents a significant hurdle for those who wish to implement Hi-C technology into their laboratory. On the other hand, genomics data repository sites sometimes contain processed Hi-C data sets, allowing researchers to perform further analysis without the need for high-spec workstations and servers. In this chapter, we provide a detailed description on how to calculate A/B compartment profiles from processed Hi-C data on the autosomes and the active/inactive X chromosomes.

Dynamic epigenomic landscapes during early lineage specification in mouse embryos

In mammals, all somatic development originates from lineage segregation in early embryos. However, the dynamics of transcriptomes and epigenomes acting in concert with initial cell fate commitment remains poorly characterized. Here we report a comprehensive investigation of transcriptomes and base-resolution methylomes for early lineages in peri- and postimplantation mouse embryos. We found allele-specific and lineage-specific de novo methylation at CG and CH sites that led to differential methylation between embryonic and extraembryonic lineages at promoters of lineage regulators, gene bodies, and DNA-methylation valleys. By using Hi-C experiments to define chromatin architecture across the same developmental period, we demonstrated that both global demethylation and remethylation in early development correlate with chromatin compartments. Dynamic local methylation was evident during gastrulation, which enabled the identification of putative regulatory elements. Finally, we found that de novo methylation patterning does not strictly require implantation. These data reveal dynamic transcriptomes, DNA methylomes, and 3D chromatin landscapes during the earliest stages of mammalian lineage specification.

Subnuclear distribution of proteins: Links with genome architecture

Metazoan genomes have a hierarchal 3-dimensional (3D) organization scaling from nucleosomes, loops, topologically associating domains (TADs), compartments, to chromosome territories. The 3D organization of genome has been linked with development, differentiation and disease. However, the principles governing the 3D chromatin architecture are just beginning to get unraveled. The nucleus has very high concentration of proteins and these proteins are either diffusely distributed throughout the nucleus, or aggregated in the form of foci/bodies/clusters/speckles or in combination of both. Several evidences suggest that the distribution of proteins within the nuclear space is linked to the organization and function of genome. Here, we describe advances made in understanding the relationship between subnuclear distribution of proteins and genome architecture.

Regulatory variants of FOXG1 in the context of its topological domain organisation

FOXG1 syndrome is caused by FOXG1 intragenic point mutations, or by long-range position effects (LRPE) of intergenic structural variants. However, the size of the FOXG1 regulatory landscape is uncertain, because the associated topologically associating domain (TAD) in fibroblasts is split into two domains in embryonic stem cells (hESC). Indeed, it has been suggested that the pathogenetic mechanism of deletions that remove the stem-cell-specific TAD boundary may be enhancer adoption due to ectopic activity of enhancer(s) located in the distal hESC-TAD. Herein we map three de novo translocation breakpoints to the proximal regulatory domain of FOXG1. The classical FOXG1 syndrome in these and in other translocation patients, and in a patient with an intergenic deletion that removes the hESC-specific TAD boundary, do not support the hypothesised enhancer adoption as a main contributor to the FOXG1 syndrome. Also, virtual 4 C and HiC-interaction data suggest that the hESC-specific TAD boundary may not be critical for FOXG1 regulation in a majority of human cells and tissues, including brain tissues and a neuronal progenitor cell line. Our data support the importance of a critical regulatory region (SRO) proximal to the hESC-specific TAD boundary. We further narrow this critical region by a deletion distal to the hESC-specific boundary, associated with a milder clinical phenotype. The distance from FOXG1 to the SRO ( > 500 kb) highlight a limitation of ENCODE DNase hypersensitivity data for functional prediction of LRPE. Moreover, the SRO has little overlap with a cluster of frequently associating regions (FIREs) located in the proximal hESC-TAD.

The nuclear matrix protein HNRNPU maintains 3D genome architecture globally in mouse hepatocytes

Eukaryotic chromosomes are folded into higher-order conformations to coordinate genome functions. In addition to long-range chromatin loops, recent chromosome conformation capture (3C)-based studies have indicated higher levels of chromatin structures including compartments and topologically associating domains (TADs), which may serve as units of genome organization and functions. However, the molecular machinery underlying these hierarchically three-dimensional (3D) chromatin architectures remains poorly understood. Via high-throughput assays, including in situ Hi-C, DamID, ChIP-seq, and RNA-seq, we investigated roles of the Heterogeneous Nuclear Ribonucleoprotein U (HNRNPU), a nuclear matrix (NM)-associated protein, in 3D genome organization. Upon the depletion of HNRNPU in mouse hepatocytes, the coverage of lamina-associated domains (LADs) in the genome increases from 53.1% to 68.6%, and a global condensation of chromatin was observed. Furthermore, disruption of HNRNPU leads to compartment switching on 7.5% of the genome, decreases TAD boundary strengths at borders between A (active) and B (inactive) compartments, and reduces chromatin loop intensities. Long-range chromatin interactions between and within compartments or TADs are also significantly remodeled upon HNRNPU depletion. Intriguingly, HNRNPU mainly associates with active chromatin, and 80% of HNRNPU peaks coincide with the binding of CTCF or RAD21. Collectively, we demonstrated that HNRNPU functions as a major factor maintaining 3D chromatin architecture, suggesting important roles of NM-associated proteins in genome organization.

Chromatin state changes during neural development revealed by in vivo cell-type specific profiling

A key question in developmental biology is how cellular differentiation is controlled during development. While transitions between trithorax-group (TrxG) and polycomb-group (PcG) chromatin states are vital for the differentiation of ES cells to multipotent stem cells, little is known regarding the role of chromatin states during development of the brain. Here we show that large-scale chromatin remodelling occurs during Drosophila neural development. We demonstrate that the majority of genes activated during neuronal differentiation are silent in neural stem cells (NSCs) and occupy black chromatin and a TrxG-repressive state. In neurons, almost all key NSC genes are switched off via HP1-mediated repression. PcG-mediated repression does not play a significant role in regulating these genes, but instead regulates lineage-specific transcription factors that control spatial and temporal patterning in the brain. Combined, our data suggest that forms of chromatin other than canonical PcG/TrxG transitions take over key roles during neural development.

While transitions between active and repressive chromatin states are essential for differentiation, little is known regarding their role during development of the brain in Drosophila. Here, the authors investigate the large scale chromatin remodelling taking place during fly neural development.

The emerging landscape of in vitro and in vivo epigenetic allelic effects

Epigenetic mechanisms that cause maternally and paternally inherited alleles to be expressed differently in offspring have the potential to radically change our understanding of the mechanisms that shape disease susceptibility, phenotypic variation, cell fate, and gene expression. However, the nature and prevalence of these effects in vivo have been unclear and are debated. Here, I consider major new studies of epigenetic allelic effects in cell lines and primary cells and in vivo. The emerging picture is that these effects take on diverse forms, and this review attempts to clarify the nature of the different forms that have been uncovered for genomic imprinting and random monoallelic expression (RME). I also discuss apparent discrepancies between in vitro and in vivo studies. Importantly, multiple studies suggest that allelic effects are prevalent and can be developmental stage- and cell type-specific. I propose some possible functions and consider roles for allelic effects within the broader context of gene regulatory networks, cellular diversity, and plasticity. Overall, the field is ripe for discovery and is in need of mechanistic and functional studies.

3D genome of multiple myeloma reveals spatial genome disorganization associated with copy number variations

The Hi-C method is widely used to study the functional roles of the three-dimensional (3D) architecture of genomes. Here, we integrate Hi-C, whole-genome sequencing (WGS) and RNA-seq to study the 3D genome architecture of multiple myeloma (MM) and how it associates with genomic variation and gene expression. Our results show that Hi-C interaction matrices are biased by copy number variations (CNVs) and can be used to detect CNVs. Also, combining Hi-C and WGS data can improve the detection of translocations. We find that CNV breakpoints significantly overlap with topologically associating domain (TAD) boundaries. Compared to normal B cells, the numbers of TADs increases by 25% in MM, the average size of TADs is smaller, and about 20% of genomic regions switch their chromatin A/B compartment types. In summary, we report a 3D genome interaction map of aneuploid MM cells and reveal the relationship among CNVs, translocations, 3D genome reorganization, and gene expression regulation.

Chromosome conformation capture techniques enable the study of genome organization in cancer cells. Here, the authors use Hi-C, WGS, and RNA-seq to study the 3D genome of multiple myeloma and find that genome disorganization is associated with copy number variations and changes in gene expression.

3D Chromatin Architecture of Large Plant Genomes Determined by Local A/B Compartments

The spatial organization of the genome plays an important role in the regulation of gene expression. However, the core structural features of animal genomes, such as topologically associated domains (TADs) and chromatin loops, are not prominent in the extremely compact Arabidopsis genome. In this study, we examine the chromatin architecture, as well as their DNA methylation, histone modifications, accessible chromatin, and gene expression, of maize, tomato, sorghum, foxtail millet, and rice with genome sizes ranging from 0.4 to 2.4 Gb. We found that these plant genomes can be divided into mammalian-like A/B compartments. At higher resolution, the chromosomes of these plants can be further partitioned to local A/B compartments that reflect their euchromatin, heterochromatin, and polycomb status. Chromatins in all these plants are organized into domains that are not conserved across species. They show similarity to the Drosophila compartment domains, and are clustered into active, polycomb, repressive, and intermediate types based on their transcriptional activities and epigenetic signatures, with domain border overlaps with the local A/B compartment junctions. In the large maize and tomato genomes, we observed extensive chromatin loops. However, unlike the mammalian chromatin loops that are enriched at the TAD border, plant chromatin loops are often formed between gene islands outside the repressive domains and are closely associated with active compartments. Our study indicates that plants have complex and unique 3D chromatin architectures, which require further study to elucidate their biological functions.

Nanog Expression in Embryonic Stem Cells - An Ideal Model System to Dissect Enhancer Function

Embryonic stem cells (ESCs) are derived from the preimplantation embryo and can differentiate into virtually any other cell type (termed pluripotency), which is governed by lineage specific transcriptions factors (TFs) binding to cis regulatory elements (CREs) to mediate changes in gene expression. The reliance on transcriptional regulation to maintain pluripotency makes ESCs a valuable model to study the role of distal CREs such as enhancers in modulating gene expression to affect cell fate decisions. This review will highlight recent advance on transcriptional enhancers, focusing on studies performed in ESCs. In addition, we argue that the Nanog locus, which encodes for an ESC-critical TF, is particularly informative because it contains multiple co-regulated genes and enhancers in close proximity to one another. The unique landscape at Nanog permits the study of ongoing questions including whether multiple enhancers function additively versus synergistically, determinants of gene specificity, and cell-to-cell variability in gene expression.

HiTAD: detecting the structural and functional hierarchies of topologically associating domains from chromatin interactions

A current question in the high-order organization of chromatin is whether topologically associating domains (TADs) are distinct from other hierarchical chromatin domains. However, due to the unclear TAD definition in tradition, the structural and functional uniqueness of TAD is not well studied. In this work, we refined TAD definition by further constraining TADs to the optimal separation on global intra-chromosomal interactions. Inspired by this constraint, we developed a novel method, called HiTAD, to detect hierarchical TADs from Hi-C chromatin interactions. HiTAD performs well in domain sensitivity, replicate reproducibility and inter cell-type conservation. With a novel domain-based alignment proposed by us, we defined several types of hierarchical TAD changes which were not systematically studied previously, and subsequently used them to reveal that TADs and sub-TADs differed statistically in correlating chromosomal compartment, replication timing and gene transcription. Finally, our work also has the implication that the refinement of TAD definition could be achieved by only utilizing chromatin interactions, at least in part. HiTAD is freely available online.

Developmentally regulated higher-order chromatin interactions orchestrate B cell fate commitment

Genome organization in 3D nuclear-space is important for regulation of gene expression. However, the alterations of chromatin architecture that impinge on the B cell-fate choice of multi-potent progenitors are still unclear. By integrating in situ Hi-C analyses with epigenetic landscapes and genome-wide expression profiles, we tracked the changes in genome architecture as the cells transit from a progenitor to a committed state. We identified the genomic loci that undergo developmental switch between A and B compartments during B-cell fate determination. Furthermore, although, topologically associating domains (TADs) are stable, a significant number of TADs display structural alterations that are associated with changes in cis-regulatory interaction landscape. Finally, we demonstrate the potential roles for Ebf1 and its downstream factor, Pax5, in chromatin reorganization and transcription regulation. Collectively, our studies provide a general paradigm of the dynamic relationship between chromatin reorganization and lineage-specific gene expression pattern that dictates cell-fate determination.

PTBP1 and PTBP2 Serve Both Specific and Redundant Functions in Neuronal Pre-mRNA Splicing

Families of alternative splicing regulators often contain multiple paralogs presumed to fulfill different functions. Polypyrimidine tract binding proteins PTBP1 and PTBP2 reprogram developmental pre-mRNA splicing in neurons, but how their regulatory networks differ is not understood. To compare their targeting, we generated a knockin allele that conditionally expresses PTBP1. Bred to a Ptbp2 knockout, the transgene allowed us to compare the developmental and molecular phenotypes of mice expressing only PTBP1, only PTBP2, or neither protein in the brain. This knockin Ptbp1 rescued a forebrain-specific, but not a pan-neuronal, Ptbp2 knockout, demonstrating both redundant and distinct roles for the proteins. Many developmentally regulated exons exhibited different sensitivities to PTBP1 and PTBP2. Nevertheless, the two paralogs displayed similar RNA binding across the transcriptome, indicating that their differential targeting does not derive from their RNA interactions, but from possible different cofactor interactions.

Allelic reprogramming of 3D chromatin architecture during early mammalian development

In mammals, chromatin organization undergoes drastic reprogramming after fertilization. However, the three-dimensional structure of chromatin and its reprogramming in preimplantation development remain poorly understood. Here, by developing a low-input Hi-C (genome-wide chromosome conformation capture) approach, we examined the reprogramming of chromatin organization during early development in mice. We found that oocytes in metaphase II show homogeneous chromatin folding that lacks detectable topologically associating domains (TADs) and chromatin compartments. Strikingly, chromatin shows greatly diminished higher-order structure after fertilization. Unexpectedly, the subsequent establishment of chromatin organization is a prolonged process that extends through preimplantation development, as characterized by slow consolidation of TADs and segregation of chromatin compartments. The two sets of parental chromosomes are spatially separated from each other and display distinct compartmentalization in zygotes. Such allele separation and allelic compartmentalization can be found as late as the 8-cell stage. Finally, we show that chromatin compaction in preimplantation embryos can partially proceed in the absence of zygotic transcription and is a multi-level hierarchical process. Taken together, our data suggest that chromatin may exist in a markedly relaxed state after fertilization, followed by progressive maturation of higher-order chromatin architecture during early development.

The HoxD cluster is a dynamic and resilient TAD boundary controlling the segregation of antagonistic regulatory landscapes

The mammalian HoxD cluster lies between two topologically associating domains (TADs) matching distinct enhancer-rich regulatory landscapes. During limb development, the telomeric TAD controls the early transcription of Hoxd genes in forearm cells, whereas the centromeric TAD subsequently regulates more posterior Hoxd genes in digit cells. Therefore, the TAD boundary prevents the terminal Hoxd13 gene from responding to forearm enhancers, thereby allowing proper limb patterning. To assess the nature and function of this CTCF-rich DNA region in embryos, we compared chromatin interaction profiles between proximal and distal limb bud cells isolated from mutant stocks where various parts of this boundary region were removed. The resulting progressive release in boundary effect triggered inter-TAD contacts, favored by the activity of the newly accessed enhancers. However, the boundary was highly resilient, and only a 400-kb deletion, including the whole-gene cluster, was eventually able to merge the neighboring TADs into a single structure. In this unified TAD, both proximal and distal limb enhancers nevertheless continued to work independently over a targeted transgenic reporter construct. We propose that the whole HoxD cluster is a dynamic TAD border and that the exact boundary position varies depending on both the transcriptional status and the developmental context.

The three-dimensional genome: regulating gene expression during pluripotency and development

The precise expression of genes in time and space during embryogenesis is largely influenced by communication between enhancers and promoters, which is propagated and governed by the physical proximity of these elements in the nucleus. Here, we review how chromatin domains organize the genome by guiding enhancers to their target genes thereby preventing non-specific interactions with other neighboring regions. We also discuss the dynamics of chromatin interactions between enhancers and promoters, as well as the consequent changes in gene expression, that occur in pluripotent cells and during development. Finally, we evaluate how genomic changes such as deletions, inversions and duplications affect 3D chromatin configuration overall and lead to ectopic enhancer-promoter contacts, and thus gene misexpression, which can contribute to abnormal development and disease.

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.

Multiscale 3D Genome Rewiring during Mouse Neural Development

Chromosome conformation capture technologies have revealed important insights into genome folding. Yet, how spatial genome architecture is related to gene expression and cell fate remains unclear. We comprehensively mapped 3D chromatin organization during mouse neural differentiation in vitro and in vivo, generating the highest-resolution Hi-C maps available to date. We found that transcription is correlated with chromatin insulation and long-range interactions, but dCas9-mediated activation is insufficient for creating TAD boundaries de novo. Additionally, we discovered long-range contacts between gene bodies of exon-rich, active genes in all cell types. During neural differentiation, contacts between active TADs become less pronounced while inactive TADs interact more strongly. An extensive Polycomb network in stem cells is disrupted, while dynamic interactions between neural transcription factors appear in vivo. Finally, cell type-specific enhancer-promoter contacts are established concomitant to gene expression. This work shows that multiple factors influence the dynamics of chromatin interactions in development.

Computational characterization of chromatin domain boundary-associated genomic elements

Topologically associated domains (TADs) are 3D genomic structures with high internal interactions that play important roles in genome compaction and gene regulation. Their genomic locations and their association with CCCTC-binding factor (CTCF)-binding sites and transcription start sites (TSSs) were recently reported. However, the relationship between TADs and other genomic elements has not been systematically evaluated. This was addressed in the present study, with a focus on the enrichment of these genomic elements and their ability to predict the TAD boundary region. We found that consensus CTCF-binding sites were strongly associated with TAD boundaries as well as with the transcription factors (TFs) Zinc finger protein (ZNF)143 and Yin Yang (YY)1. TAD boundary-associated genomic elements include DNase I-hypersensitive sites, H3K36 trimethylation, TSSs, RNA polymerase II, and TFs such as Specificity protein 1, ZNF274 and SIX homeobox 5. Computational modeling with these genomic elements suggests that they have distinct roles in TAD boundary formation. We propose a structural model of TAD boundaries based on these findings that provides a basis for studying the mechanism of chromatin structure formation and gene regulation.

Epigenomics of Retinal Development in Mice and Humans

In this issue of Neuron, Aldiri et al. (2017) present an analysis of epigenetic changes during retinal development, and use these data to probe reprogramming of retinal iPSC cells, as well as the origin of retinoblastoma cells.

Utilizing networks for differential analysis of chromatin interactions

Chromatin conformation capture with high-throughput sequencing (Hi-C) is a powerful technique to detect genome-wide chromatin interactions. In this paper, we introduce two novel approaches to detect differentially interacting genomic regions between two Hi-C experiments using a network model. To make input data from multiple experiments comparable, we propose a normalization strategy guided by network topological properties. We then devise two measurements, using local and global connectivity information from the chromatin interaction networks, respectively, to assess the interaction differences between two experiments. When multiple replicates are present in experiments, our approaches provide the flexibility for users to either pool all replicates together to therefore increase the network coverage, or to use the replicates in parallel to increase the signal to noise ratio. We show that while the local method works better in detecting changes from simulated networks, the global method performs better on real Hi-C data. The local and global methods, regardless of pooling, are always superior to two existing methods. Furthermore, our methods work well on both unweighted and weighted networks and our normalization strategy significantly improves the performance compared with raw networks without normalization. Therefore, we believe our methods will be useful for identifying differentially interacting genomic regions.

Analysis of primary microRNA loci from nascent transcriptomes reveals regulatory domains governed by chromatin architecture

Changes in mature microRNA (miRNA) levels that occur downstream of signaling cascades play an important role during human development and disease. However, the regulation of primary microRNA (pri-miRNA) genes remains to be dissected in detail. To address this, we followed a data-driven approach and developed a transcript identification, validation and quantification pipeline for characterizing the regulatory domains of pri-miRNAs. Integration of 92 nascent transcriptomes and multilevel data from cells arising from ecto-, endo- and mesoderm lineages reveals cell type-specific expression patterns, allows fine-resolution mapping of transcription start sites (TSS) and identification of candidate regulatory regions. We show that inter- and intragenic pri-miRNA transcripts span vast genomic regions and active TSS locations differ across cell types, exemplified by the mir-29a∼29b-1, mir-100∼let-7a-2∼125b-1 and miR-221∼222 clusters. Considering the presence of multiple TSS as an important regulatory feature at miRNA loci, we developed a strategy to quantify differential TSS usage. We demonstrate that the TSS activities associate with cell type-specific super-enhancers, differential stimulus responsiveness and higher-order chromatin structure. These results pave the way for building detailed regulatory maps of miRNA loci.

The Role of Chromatin Density in Cell Population Heterogeneity during Stem Cell Differentiation

We incorporate three-dimensional (3D) conformation of chromosome (Hi-C) and single-cell RNA sequencing data together with discrete stochastic simulation, to explore the role of chromatin reorganization in determining gene expression heterogeneity during development. While previous research has emphasized the importance of chromatin architecture on activation and suppression of certain regulatory genes and gene networks, our study demonstrates how chromatin remodeling can dictate gene expression distribution by folding into distinct topological domains. We hypothesize that the local DNA density during differentiation accentuate transcriptional bursting due to the crowding effect of chromatin. This phenomenon yields a heterogeneous cell population, thereby increasing the potential of differentiation of the stem cells.

Cell Type-Specific Chromatin Signatures Underline Regulatory DNA Elements in Human Induced Pluripotent Stem Cells and Somatic Cells

RATIONALE

Regulatory DNA elements in the human genome play important roles in determining the transcriptional abundance and spatiotemporal gene expression during embryonic heart development and somatic cell reprogramming. It is not well known how chromatin marks in regulatory DNA elements are modulated to establish cell type-specific gene expression in the human heart.

OBJECTIVE

We aimed to decipher the cell type-specific epigenetic signatures in regulatory DNA elements and how they modulate heart-specific gene expression.

METHODS AND RESULTS

We profiled genome-wide transcriptional activity and a variety of epigenetic marks in the regulatory DNA elements using massive RNA-seq (n=12) and ChIP-seq (chromatin immunoprecipitation combined with high-throughput sequencing; n=84) in human endothelial cells (CD31⁺CD144⁺), cardiac progenitor cells (Sca-1⁺), fibroblasts (DDR2⁺), and their respective induced pluripotent stem cells. We uncovered 2 classes of regulatory DNA elements: class I was identified with ubiquitous enhancer (H3K4me1) and promoter (H3K4me3) marks in all cell types, whereas class II was enriched with H3K4me1 and H3K4me3 in a cell type-specific manner. Both class I and class II regulatory elements exhibited stimulatory roles in nearby gene expression in a given cell type. However, class I promoters displayed more dominant regulatory effects on transcriptional abundance regardless of distal enhancers. Transcription factor network analysis indicated that human induced pluripotent stem cells and somatic cells from the heart selected their preferential regulatory elements to maintain cell type-specific gene expression. In addition, we validated the function of these enhancer elements in transgenic mouse embryos and human cells and identified a few enhancers that could possibly regulate the cardiac-specific gene expression.

CONCLUSIONS

Given that a large number of genetic variants associated with human diseases are located in regulatory DNA elements, our study provides valuable resources for deciphering the epigenetic modulation of regulatory DNA elements that fine-tune spatiotemporal gene expression in human cardiac development and diseases.