Chromatin architecture reorganization during stem cell differentiation

The Temporal Regulation of S Phase Proteins During G1

Successful DNA replication requires intimate coordination with cell-cycle progression. Prior to DNA replication initiation in S phase, a series of essential preparatory events in G₁ phase ensures timely, complete, and precise genome duplication. Among the essential molecular processes are regulated transcriptional upregulation of genes that encode replication proteins, appropriate post-transcriptional control of replication factor abundance and activity, and assembly of DNA-loaded protein complexes to license replication origins. In this chapter we describe these critical G₁ events necessary for DNA replication and their regulation in the context of both cell-cycle entry and cell-cycle progression.

A promoter interaction map for cardiovascular disease genetics

Over 500 genetic loci have been associated with risk of cardiovascular diseases (CVDs); however, most loci are located in gene-distal non-coding regions and their target genes are not known. Here, we generated high-resolution promoter capture Hi-C (PCHi-C) maps in human induced pluripotent stem cells (iPSCs) and iPSC-derived cardiomyocytes (CMs) to provide a resource for identifying and prioritizing the functional targets of CVD associations. We validate these maps by demonstrating that promoters preferentially contact distal sequences enriched for tissue-specific transcription factor motifs and are enriched for chromatin marks that correlate with dynamic changes in gene expression. Using the CM PCHi-C map, we linked 1999 CVD-associated SNPs to 347 target genes. Remarkably, more than 90% of SNP-target gene interactions did not involve the nearest gene, while 40% of SNPs interacted with at least two genes, demonstrating the importance of considering long-range chromatin interactions when interpreting functional targets of disease loci.

Novel long-range regulatory mechanisms controlling PKD2 gene expression

Background

Cis-regulatory elements control gene expression over large distances through the formation of chromatin loops, which allow contact between enhancers and gene promoters. Alterations in cis-acting regulatory systems could be linked to human genetic diseases. Here, we analyse the spatial organization of a large region spanning the polycystic kidney disease 2 (PKD2) gene, one of the genes responsible of autosomal dominant polycystic kidney disease (ADPKD).

Results

By using chromosome conformation capture carbon copy (5C) technology in primary human renal cyst epithelial cells, we identify novel contacts of the PKD2 promoter with chromatin regions, which display characteristics of regulatory elements. In parallel, by using functional analysis with a reporter assay, we demonstrate that three DNAse I hypersensitive sites regions are involved in the regulation of PKD2 gene expression.

Conclusions

Finally, through alignment of CCCTC-binding factor (CTCF) sites, we suggest that these novel enhancer elements are brought to the PKD2 promoter by chromatin looping via the recruitment of CTCF.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4892-6) contains supplementary material, which is available to authorized users.

Epigenetic maintenance of topological domains in the highly rearranged gibbon genome

The relationship between evolutionary genome remodeling and the three-dimensional structure of the genome remain largely unexplored. Here, we use the heavily rearranged gibbon genome to examine how evolutionary chromosomal rearrangements impact genome-wide chromatin interactions, topologically associating domains (TADs), and their epigenetic landscape. We use high-resolution maps of gibbon-human breaks of synteny (BOS), apply Hi-C in gibbon, measure an array of epigenetic features, and perform cross-species comparisons. We find that gibbon rearrangements occur at TAD boundaries, independent of the parameters used to identify TADs. This overlap is supported by a remarkable genetic and epigenetic similarity between BOS and TAD boundaries, namely presence of CpG islands and SINE elements, and enrichment in CTCF and H3K4me3 binding. Cross-species comparisons reveal that regions orthologous to BOS also correspond with boundaries of large (400-600 kb) TADs in human and other mammalian species. The colocalization of rearrangement breakpoints and TAD boundaries may be due to higher chromatin fragility at these locations and/or increased selective pressure against rearrangements that disrupt TAD integrity. We also examine the small portion of BOS that did not overlap with TAD boundaries and gave rise to novel TADs in the gibbon genome. We postulate that these new TADs generally lack deleterious consequences. Last, we show that limited epigenetic homogenization occurs across breakpoints, irrespective of their time of occurrence in the gibbon lineage. Overall, our findings demonstrate remarkable conservation of chromatin interactions and epigenetic landscape in gibbons, in spite of extensive genomic shuffling.

Unsupervised embedding of single-cell Hi-C data

Motivation

Single-cell Hi-C (scHi-C) data promises to enable scientists to interrogate the 3D architecture of DNA in the nucleus of the cell, studying how this structure varies stochastically or along developmental or cell-cycle axes. However, Hi-C data analysis requires methods that take into account the unique characteristics of this type of data. In this work, we explore whether methods that have been developed previously for the analysis of bulk Hi-C data can be applied to scHi-C data. We apply methods designed for analysis of bulk Hi-C data to scHi-C data in conjunction with unsupervised embedding.

Results

We find that one of these methods, HiCRep, when used in conjunction with multidimensional scaling (MDS), strongly outperforms three other methods, including a technique that has been used previously for scHi-C analysis. We also provide evidence that the HiCRep/MDS method is robust to extremely low per-cell sequencing depth, that this robustness is improved even further when high-coverage and low-coverage cells are projected together, and that the method can be used to jointly embed cells from multiple published datasets.

Three-dimensional genome organization in normal and malignant haematopoiesis

PURPOSE OF REVIEW

The three-dimensional organization of the genome inside the nucleus impacts on key aspects of genome function, including transcription, DNA replication and repair. The chromosome maintenance complex cohesin and the DNA binding protein CTCF cooperate to drive the formation of self-interacting topological domains. This facilitates transcriptional regulation via enhancer-promoter interactions, controls the distribution and release of torsional strain, and affects the frequency with which particular translocations arise, based on the spatial proximity of translocation partners. Here we discuss recent insights into the mechanisms of three-dimensional genome organization, their relationship to haematopoietic differentiation and malignant transformation.

RECENT FINDINGS

Cohesin mutations are frequently found in myeloid malignancies. Significantly, cohesin mutations can drive increased self-renewal of haematopoietic stem and progenitor cells, which may facilitate the accumulation of genetic lesions and leukaemic transformation. It is therefore important to elucidate the mechanisms that link cohesin to pathways that regulate the balance between self-renewal and differentiation. Chromosomal translocations are key to lymphoid malignancies, and recent findings link three-dimensional genome organization to the frequency and the genomic position of DNA double strand breaks.

SUMMARY

Three-dimensional genome organization can help explain genome function in normal and malignant haematopoiesis.

Identification of mutations in patients with acquired pure red cell aplasia

Idiopathic acquired pure red cell aplasia (PRCA) is a rare, autoimmune-related disease. This study aimed to describe the previously unidentified DNA alterations associated with PRCA. Here, next generation sequencing using a panel containing 295 critical genes was applied to detect potentially pathogenic mutations in four patients with PRCA. A total of 529 mutations were identified and further classified into three categories, namely, uncertain (n = 25), likely benign (n = 20) and benign (n = 484) mutations, based on the American College of Medical Genetics and Genomics (ACMG) 2015 guidelines and ClinVar database. The spatial proximity between two loci of the uncertain or benign mutations was evaluated using Hi-C datasets of KBM7 and K562 cell lines, respectively. Significant spatial proximity was observed in uncertain mutation pairs compared with benign mutation pairs. In addition, 17 variants were eventually identified after excluding those with mutant frequencies >0.001, including 7 newly identified variants. FANCF and LRP1B mutations existed twice in patients. FANCF and LRP1B mutations were likely to affect protein stability based on prediction analysis. Taken together, our data may provide valuable information about PRCA. FANCF and LRP1B mutations may be associated with acquired PRCA.

An atlas of chromatin accessibility in the adult human brain

Most common genetic risk variants associated with neuropsychiatric disease are noncoding and are thought to exert their effects by disrupting the function of cis regulatory elements (CREs), including promoters and enhancers. Within each cell, chromatin is arranged in specific patterns to expose the repertoire of CREs required for optimal spatiotemporal regulation of gene expression. To further understand the complex mechanisms that modulate transcription in the brain, we used frozen postmortem samples to generate the largest human brain and cell-type–specific open chromatin data set to date. Using the Assay for Transposase Accessible Chromatin followed by sequencing (ATAC-seq), we created maps of chromatin accessibility in two cell types (neurons and non-neurons) across 14 distinct brain regions of five individuals. Chromatin structure varies markedly by cell type, with neuronal chromatin displaying higher regional variability than that of non-neurons. Among our findings is an open chromatin region (OCR) specific to neurons of the striatum. When placed in the mouse, a human sequence derived from this OCR recapitulates the cell type and regional expression pattern predicted by our ATAC-seq experiments. Furthermore, differentially accessible chromatin overlaps with the genetic architecture of neuropsychiatric traits and identifies differences in molecular pathways and biological functions. By leveraging transcription factor binding analysis, we identify protein-coding and long noncoding RNAs (lncRNAs) with cell-type and brain region specificity. Our data provide a valuable resource to the research community and we provide this human brain chromatin accessibility atlas as an online database “Brain Open Chromatin Atlas (BOCA)” to facilitate interpretation.

Dendritic Cells Actively Limit Interleukin-10 Production Under Inflammatory Conditions via DC-SCRIPT and Dual-Specificity Phosphatase 4

Dendritic cell (DC)-based immunotherapy makes use of the DC’s ability to direct the adaptive immune response toward activation or inhibition. DCs perform this immune orchestration in part by secretion of selected cytokines. The most potent anti-inflammatory cytokine interleukin-10 (IL-10) is under tight regulation, as it needs to be predominantly expressed during the resolution phase of the immune response. Currently it is not clear whether there is active suppression of IL-10 by DCs at the initial pro-inflammatory stage of the immune response. Previously, knockdown of the DC-specific transcription factor DC-SCRIPT has been demonstrated to mediate an extensive increase in IL-10 production upon encounter with pro-inflammatory immune stimuli. Here, we explored how DC-SCRIPT contributes to IL-10 suppression under pro-inflammatory conditions by applying chromatin immunoprecipitation sequencing analysis of DC-SCRIPT and the epigenetic marks H3K4me3 and H3K27ac in human DCs. The data showed binding of DC-SCRIPT to a GA-rich motif at H3K27ac-marked genomic enhancers that associated with genes encoding MAPK dual-specificity phosphatases (DUSPs). Functional studies revealed that upon knockdown of DC-SCRIPT, human DCs express much less DUSP4 and exhibit increased phosphorylation of the three major MAPKs (ERK, JNK, and p38). Enhanced ERK signaling in DC-SCRIPT-knockdown-DCs led to higher production of IL-10, which was reverted by rescuing DUSP4 expression. Finally, DC-SCRIPT-knockdown-DCs induced less IFN-γ and increased IL-10 production in naïve T cells, indicative for a more anti-inflammatory phenotype. In conclusion, we have delineated a new mechanism by which DC-SCRIPT allows DCs to limit IL-10 production under inflammatory conditions and potentiate pro-inflammatory Th1 responses. These insights may be exploited to improve DC-based immunotherapies.

Understanding direct neuronal reprogramming-from pioneer factors to 3D chromatin

Cell replacement therapies aim at reestablishment of neuronal circuits after brain injury, stroke or neurodegeneration. Recently, direct reprogramming of resident glial cells into the affected neuronal subtypes has become a feasible and promising option for central nervous system regeneration. Direct reprogramming relies on the implementation of a new transcriptional program defining the desired neuronal identity in fully differentiated glial cells implying the more or less complete down-regulation of the program for the former identity of the glial cell. Despite the enormous progress achieved in this regard with highly efficient in vivo reprogramming after injury, a number of hurdles still need to be resolved. One way to further improve direct neuronal reprogramming is to understand the molecular hurdles which we discuss with the focus on chromatin states of the starting versus the reprogrammed cells.

Interpreting Chromosomal Rearrangements in the Context of 3-Dimentional Genome Organization: A Practical Guide for Medical Genetics

In this exciting era of "next-gen cytogenetics", the use of novel molecular methods such as comparative genome hybridization and whole genome and whole exome sequencing becomes more and more common in clinics. This results in generation of large amounts of high-resolution patient-specific data and challenges the development of new approaches for interpretation of obtained information. Usually, interpretation of chromosomal rearrangements is focused on alterations of linear genome sequence, underestimating the role of spatial chromatin organization. In this article, we describe the main features of 3-dimentional genome organization, emphasizing their role in normal and pathological development. We highlight some tips to help physicians estimating the impact of chromosomal rearrangements on the patient phenotype. A separate section describes available tools that can be used to visualize and analyze human genome architecture.

Low-Grade Astrocytoma Mutations in IDH1, P53, and ATRX Cooperate to Block Differentiation of Human Neural Stem Cells via Repression of SOX2

Low-grade astrocytomas (LGAs) carry neomorphic mutations in isocitrate dehydrogenase (IDH) concurrently with P53 and ATRX loss. To model LGA formation, we introduced R132H IDH1, P53 shRNA, and ATRX shRNA into human neural stem cells (NSCs). These oncogenic hits blocked NSC differentiation, increased invasiveness in vivo, and led to a DNA methylation and transcriptional profile resembling IDH1 mutant human LGAs. The differentiation block was caused by transcriptional silencing of the transcription factor SOX2 secondary to disassociation of its promoter from a putative enhancer. This occurred because of reduced binding of the chromatin organizer CTCF to its DNA motifs and disrupted chromatin looping. Our human model of IDH mutant LGA formation implicates impaired NSC differentiation because of repression of SOX2 as an early driver of gliomagenesis.

RNAs as Proximity-Labeling Media for Identifying Nuclear Speckle Positions Relative to the Genome

It remains challenging to identify all parts of the nuclear genome that are in proximity to nuclear speckles, due to physical separation between the nuclear speckle cores and chromatin. We hypothesized that noncoding RNAs including small nuclear RNA (snRNAs) and Malat1, which accumulate at the periphery of nuclear speckles (nsaRNA [nuclear speckle-associated RNA]), may extend to sufficient proximity to the genome. Leveraging a transcriptome-genome interaction assay (mapping of RNA-genome interactions [MARGI]), we identified clusters of nsaRNA-interacting genomic sequences (nsaPeaks). Posttranscriptional pre-mRNAs, which also accumulate to nuclear speckles, exhibited proximity to nsaPeaks but rarely to other genomic regions. Our combined DNA fluorescence in situ hybridization and immunofluorescence analysis in 182 single cells revealed a 3-fold increase in odds for nuclear speckles to localize near an nsaPeak than its neighboring genomic sequence. These data suggest a model that nsaRNAs are located in sufficient proximity to the nuclear genome and leave identifiable genomic footprints, thus revealing the parts of genome proximal to nuclear speckles.

•

MARGI captures interactions of nuclear speckle-associated RNAs (nsaRNA) and DNA

•

nsaRNA-interacting genomic sequences were clustered (nsaPeaks) in the genome

•

Posttranscriptional pre-mRNAs and CDK9 proteins exhibited proximity to nsaPeaks

•

Single-cell images confirmed proximity of nuclear speckles to an nsaPeak

SMCHD1 Merges Chromosome Compartments and Assists Formation of Super-Structures on the Inactive X

Mammalian chromosomes are partitioned into A/B compartments and topologically associated domains (TADs). The inactive X (Xi) chromosome, however, adopts a distinct conformation without evident compartments or TADs. Here, through exploration of an architectural protein, structural-maintenance-of-chromosomes hinge domain containing 1 (SMCHD1), we probe how the Xi is reconfigured during X chromosome inactivation. A/B compartments are first fused into "S1" and "S2" compartments, coinciding with Xist spreading into gene-rich domains. SMCHD1 then binds S1/S2 compartments and merges them to create a compartment-less architecture. Contrary to current views, TADs remain on the Xi but in an attenuated state. Ablating SMCHD1 results in a persistent S1/S2 organization and strengthening of TADs. Furthermore, loss of SMCHD1 causes regional defects in Xist spreading and erosion of heterochromatic silencing. We present a stepwise model for Xi folding, where SMCHD1 attenuates a hidden layer of Xi architecture to facilitate Xist spreading.

Forces driving the three-dimensional folding of eukaryotic genomes

The last decade has radically renewed our understanding of higher order chromatin folding in the eukaryotic nucleus. As a result, most current models are in support of a mostly hierarchical and relatively stable folding of chromosomes dividing chromosomal territories into A‐ (active) and B‐ (inactive) compartments, which are then further partitioned into topologically associating domains (TADs), each of which is made up from multiple loops stabilized mainly by the CTCF and cohesin chromatin‐binding complexes. Nonetheless, the structure‐to‐function relationship of eukaryotic genomes is still not well understood. Here, we focus on recent work highlighting the biophysical and regulatory forces that contribute to the spatial organization of genomes, and we propose that the various conformations that chromatin assumes are not so much the result of a linear hierarchy, but rather of both converging and conflicting dynamic forces that act on it.

Distinct roles of cohesin-SA1 and cohesin-SA2 in 3D chromosome organization

Two variant cohesin complexes containing SMC1, SMC3, RAD21 and either SA1 (also known as STAG1) or SA2 (also known as STAG2) are present in all cell types. We report here their genomic distribution and specific contributions to genome organization in human cells. Although both variants are found at CCCTC-binding factor (CTCF) sites, a distinct population of the SA2-containing cohesin complexes (hereafter referred to as cohesin-SA2) localize to enhancers lacking CTCF, are linked to tissue-specific transcription and cannot be replaced by the SA1-containing cohesin complex (cohesin-SA1) when SA2 is absent, a condition that has been observed in several tumors. Downregulation of each of these variants has different consequences for gene expression and genome architecture. Our results suggest that cohesin-SA1 preferentially contributes to the stabilization of topologically associating domain boundaries together with CTCF, whereas cohesin-SA2 promotes cell-type-specific contacts between enhancers and promoters independently of CTCF. Loss of cohesin-SA2 rewires local chromatin contacts and alters gene expression. These findings provide insights into how cohesin mediates chromosome folding and establish a novel framework to address the consequences of mutations in cohesin genes in cancer.

A Decade of Exploring the Mammalian Sperm Epigenome: Paternal Epigenetic and Transgenerational Inheritance

The past decade has seen a tremendous increase in interest and progress in the field of sperm epigenetics. Studies have shown that chromatin regulation during male germline development is multiple and complex, and that the spermatozoon possesses a unique epigenome. Its DNA methylation profile, DNA-associated proteins, nucleo-protamine distribution pattern and non-coding RNA set up a unique epigenetic landscape which is delivered, along with its haploid genome, to the oocyte upon fertilization, and therefore can contribute to embryogenesis and to the offspring health. An emerging body of compelling data demonstrates that environmental exposures and paternal lifestyle can change the sperm epigenome and, consequently, may affect both the embryonic developmental program and the health of future generations. This short review will attempt to provide an overview of what is currently known about sperm epigenome and the existence of transgenerational epigenetic inheritance of paternally acquired traits that may contribute to the offspring phenotype.

Computational prediction of CTCF/cohesin-based intra-TAD loops that insulate chromatin contacts and gene expression in mouse liver

CTCF and cohesin are key drivers of 3D-nuclear organization, anchoring the megabase-scale Topologically Associating Domains (TADs) that segment the genome. Here, we present and validate a computational method to predict cohesin-and-CTCF binding sites that form intra-TAD DNA loops. The intra-TAD loop anchors identified are structurally indistinguishable from TAD anchors regarding binding partners, sequence conservation, and resistance to cohesin knockdown; further, the intra-TAD loops retain key functional features of TADs, including chromatin contact insulation, blockage of repressive histone mark spread, and ubiquity across tissues. We propose that intra-TAD loops form by the same loop extrusion mechanism as the larger TAD loops, and that their shorter length enables finer regulatory control in restricting enhancer-promoter interactions, which enables selective, high-level expression of gene targets of super-enhancers and genes located within repressive nuclear compartments. These findings elucidate the role of intra-TAD cohesin-and-CTCF binding in nuclear organization associated with widespread insulation of distal enhancer activity.

Exploiting sequence-based features for predicting enhancer-promoter interactions

Motivation

A large number of distal enhancers and proximal promoters form enhancer–promoter interactions to regulate target genes in the human genome. Although recent high-throughput genome-wide mapping approaches have allowed us to more comprehensively recognize potential enhancer–promoter interactions, it is still largely unknown whether sequence-based features alone are sufficient to predict such interactions.

Results

Here, we develop a new computational method (named PEP) to predict enhancer–promoter interactions based on sequence-based features only, when the locations of putative enhancers and promoters in a particular cell type are given. The two modules in PEP (PEP-Motif and PEP-Word) use different but complementary feature extraction strategies to exploit sequence-based information. The results across six different cell types demonstrate that our method is effective in predicting enhancer–promoter interactions as compared to the state-of-the-art methods that use functional genomic signals. Our work demonstrates that sequence-based features alone can reliably predict enhancer–promoter interactions genome-wide, which could potentially facilitate the discovery of important sequence determinants for long-range gene regulation.

Availability and Implementation

The source code of PEP is available at: https://github.com/ma-compbio/PEP.

Supplementary information

Supplementary data are available at Bioinformatics online.

The Post-GWAS Era: From Association to Function

During the past 12 years, genome-wide association studies (GWASs) have uncovered thousands of genetic variants that influence risk for complex human traits and diseases. Yet functional studies aimed at delineating the causal genetic variants and biological mechanisms underlying the observed statistical associations with disease risk have lagged. In this review, we highlight key advances in the field of functional genomics that may facilitate the derivation of biological meaning post-GWAS. We highlight the evidence suggesting that causal variants underlying disease risk often function through regulatory effects on the expression of target genes and that these expression effects might be modest and cell-type specific. We moreover discuss specific studies as proof-of-principle examples for current statistical, bioinformatic, and empirical bench-based approaches to downstream elucidation of GWAS-identified disease risk loci.

Diffuse gliomas classified by 1p/19q co-deletion, TERT promoter and IDH mutation status are associated with specific genetic risk loci

Recent genome-wide association studies of glioma have led to the discovery of single nucleotide polymorphisms (SNPs) at 25 loci influencing risk. Gliomas are heterogeneous, hence to investigate the relationship between risk SNPs and glioma subtype we analysed 1659 tumours profiled for IDH mutation, TERT promoter mutation and 1p/19q co-deletion. These data allowed definition of five molecular subgroups of glioma: triple-positive (IDH mutated, 1p/19q co-deletion, TERT promoter mutated); TERT-IDH (IDH mutated, TERT promoter mutated, 1p/19q-wild-type); IDH-only (IDH mutated, 1p/19q wild-type, TERT promoter wild-type); triple-negative (IDH wild-type, 1p/19q wild-type, TERT promoter wild-type) and TERT-only (TERT promoter mutated, IDH wild-type, 1p/19q wild-type). Most glioma risk loci showed subtype specificity: (1) the 8q24.21 SNP for triple-positive glioma; (2) 5p15.33, 9p21.3, 17p13.1 and 20q13.33 SNPs for TERT-only glioma; (3) 1q44, 2q33.3, 3p14.1, 11q21, 11q23.3, 14q12, and 15q24.2 SNPs for IDH mutated glioma. To link risk SNPs to target candidate genes we analysed Hi-C and gene expression data, highlighting the potential role of IDH1 at 2q33.3, MYC at 8q24.21 and STMN3 at 20q13.33. Our observations provide further insight into the nature of susceptibility to glioma.

Hypoxia and Chromatin: A Focus on Transcriptional Repression Mechanisms

Hypoxia or reduced oxygen availability has been studied extensively for its ability to activate specific genes. Hypoxia-induced gene expression is mediated by the HIF transcription factors, but not exclusively so. Despite the extensive knowledge about how hypoxia activates genes, much less is known about how hypoxia promotes gene repression. In this review, we discuss the potential mechanisms underlying hypoxia-induced transcriptional repression responses. We highlight HIF-dependent and independent mechanisms as well as the potential roles of dioxygenases with functions at the nucleosome and DNA level. Lastly, we discuss recent evidence regarding the involvement of transcriptional repressor complexes in hypoxia.

Dosage Compensation of the X Chromosome: A Complex Epigenetic Assignment Involving Chromatin Regulators and Long Noncoding RNAs

X chromosome regulation represents a prime example of an epigenetic phenomenon where coordinated regulation of a whole chromosome is required. In flies, this is achieved by transcriptional upregulation of X chromosomal genes in males to equalize the gene dosage differences in females. Chromatin-bound proteins and long noncoding RNAs (lncRNAs) constituting a ribonucleoprotein complex known as the male-specific lethal (MSL) complex or the dosage compensation complex mediate this process. MSL complex members decorate the male X chromosome, and their absence leads to male lethality. The male X chromosome is also enriched with histone H4 lysine 16 acetylation (H4K16ac), indicating that the chromatin compaction status of the X chromosome also plays an important role in transcriptional activation. How the X chromosome is specifically targeted and how dosage compensation is mechanistically achieved are central questions for the field. Here, we review recent advances, which reveal a complex interplay among lncRNAs, the chromatin landscape, transcription, and chromosome conformation that fine-tune X chromosome gene expression.

Epigenetics and epitranscriptomics in temporal patterning of cortical neural progenitor competence

Yoon et al. review epigenetic and epitranscriptomic mechanisms that regulate the lineage specification of neural progenitor cells in the developing brain.

During embryonic brain development, neural progenitor/stem cells (NPCs) sequentially give rise to different subtypes of neurons and glia via a highly orchestrated process. To accomplish the ordered generation of distinct progenies, NPCs go through multistep transitions of their developmental competence. The molecular mechanisms driving precise temporal coordination of these transitions remains enigmatic. Epigenetic regulation, including changes in chromatin structures, DNA methylation, and histone modifications, has been extensively investigated in the context of cortical neurogenesis. Recent studies of chemical modifications on RNA, termed epitranscriptomics, have also revealed their critical roles in neural development. In this review, we discuss advances in understanding molecular regulation of the sequential lineage specification of NPCs in the embryonic mammalian brain with a focus on epigenetic and epitranscriptomic mechanisms. In particular, the discovery of lineage-specific gene transcripts undergoing rapid turnover in NPCs suggests that NPC developmental fate competence is determined much earlier, before the final cell division, and is more tightly controlled than previously appreciated. We discuss how multiple regulatory systems work in harmony to coordinate NPC behavior and summarize recent findings in the context of a model of epigenetic and transcriptional prepatterning to explain NPC developmental competence.

Chromatin Architectural Changes during Cellular Senescence and Aging

Chromatin 3D structure is highly dynamic and associated with many biological processes, such as cell cycle progression, cellular differentiation, cell fate reprogramming, cancer development, cellular senescence, and aging. Recently, by using chromosome conformation capture technologies, tremendous findings have been reported about the dynamics of genome architecture, their associated proteins, and the underlying mechanisms involved in regulating chromatin spatial organization and gene expression. Cellular senescence and aging, which involve multiple cellular and molecular functional declines, also undergo significant chromatin structural changes, including alternations of heterochromatin and disruption of higher-order chromatin structure. In this review, we summarize recent findings related to genome architecture, factors regulating chromatin spatial organization, and how they change during cellular senescence and aging.

DNA Methylation and Psychiatric Disorders

DNA methylation has been an important area of research in the study of molecular mechanism to psychiatric disorders. Recent evidence has suggested that abnormalities in global methylation, methylation of genes, and pathways could play a role in the etiology of many forms of mental illness. In this article, we review the mechanisms of DNA methylation, including the genetic and environmental factors affecting methylation changes. We report and discuss major findings regarding DNA methylation in psychiatric patients, both within the context of global methylation studies and gene-specific methylation studies. Finally, we discuss issues surrounding data quality improvement, the limitations of current methylation analysis methods, and the possibility of using DNA methylation-based treatment for psychiatric disorders in the future.

Srf destabilizes cellular identity by suppressing cell-type-specific gene expression programs

Multicellular organisms consist of multiple cell types. The identity of these cells is primarily maintained by cell-type-specific gene expression programs; however, mechanisms that suppress these programs are poorly defined. Here we show that serum response factor (Srf), a transcription factor that is activated by various extracellular stimuli, can repress cell-type-specific genes and promote cellular reprogramming to pluripotency. Manipulations that decrease β-actin monomer quantity result in the nuclear accumulation of Mkl1 and the activation of Srf, which downregulate cell-type-specific genes and alter the epigenetics of regulatory regions and chromatin organization. Mice overexpressing Srf exhibit various pathologies including an ulcerative colitis-like symptom and a metaplasia-like phenotype in the pancreas. Our results demonstrate an unexpected function of Srf via a mechanism by which extracellular stimuli actively destabilize cell identity and suggest Srf involvement in a wide range of diseases.

The transcription factor Srf is a central regulator of immediate-early and actin cytoskeletal genes. Here the authors show that Srf is activated by repression of β-actin, promoting iPSC reprogramming of neural progenitor cells and hepatoblasts by repressing cell-type specific genes.

Genomes in Focus: Development and Applications of CRISPR-Cas9 Imaging Technologies

The discovery of the CRISPR-Cas9 endonuclease has enabled facile genome editing in living cells and organisms. Catalytically inactive Cas9 (dCas9) retains the ability to bind DNA in an RNA-guided fashion, and has additionally been explored as a tool for transcriptional modulation, epigenetic editing, and genome imaging. This Review highlights recent progress and challenges in the development of dCas9 for imaging genomic loci. The emergence and maturation of this technology offers the potential to answer mechanistic questions about chromosome dynamics and three-dimensional genome organization in vivo.

Producing genome structure populations with the dynamic and automated PGS software

Chromosome conformation capture technologies such as Hi-C are widely used to investigate the spatial organization of genomes. Because genome structures can vary considerably between individual cells of a population, interpreting ensemble-averaged Hi-C data can be challenging, in particular for long-range and interchromosomal interactions. We pioneered a probabilistic approach for the generation of a population of distinct diploid 3D genome structures consistent with all the chromatin-chromatin interaction probabilities from Hi-C experiments. Each structure in the population is a physical model of the genome in 3D. Analysis of these models yields new insights into the causes and the functional properties of the genome's organization in space and time. We provide a user-friendly software package, called PGS, which runs on local machines (for practice runs) and high-performance computing platforms. PGS takes a genome-wide Hi-C contact frequency matrix, along with information about genome segmentation, and produces an ensemble of 3D genome structures entirely consistent with the input. The software automatically generates an analysis report, and provides tools to extract and analyze the 3D coordinates of specific domains. Basic Linux command-line knowledge is sufficient for using this software. A typical running time of the pipeline is ∼3 d with 300 cores on a computer cluster to generate a population of 1,000 diploid genome structures at topological-associated domain (TAD)-level resolution.

Chromosome positioning in interphase nuclei of hematopoietic stem cell and myeloid precursor

Human myelopoiesis is an intriguing biological process during which multipotent stem cells limit their differentiation potential generating precursors that evolve into terminally differentiated cells. The differentiation process is correlated with differential gene expression and changes in nuclear architecture. In interphase, chromosomes are distinct entities known as chromosome territories and they show a radial localization that could result in a constrain of inter-homologous distance. This element plays a role in genome stability and gene expression. Here, we provide the first experimental evidence of 3D chromosomal arrangement considering two steps of human normal myelopoiesis. Specifically, multicolor 3D-FISH and 3D image analysis revealed that, in both normal human hematopoietic stem cells and myelod precursors CD14-, chromosomal position is correlated with gene density. However, we observed that inter-homologue distances are totally different during differentiation. This could be associated with differential gene expression that we found comparing the two cell types. Our results disclose an unprecedented framework relevant for deciphering the genomic mechanisms at the base of normal human myelopoiesis.

Tethering of CHROMATOR and dCTCF proteins results in decompaction of condensed bands in the Drosophila melanogaster polytene chromosomes but does not affect their transcription and replication timing

Instulator proteins are central to domain organization and gene regulation in the genome. We used ectopic tethering of CHROMATOR (CHRIZ/CHRO) and dCTCF to pre-defined regions of the genome to dissect the influence of these proteins on local chromatin organization, to analyze their interaction with other key chromatin proteins and to evaluate the effects on transcription and replication. Specifically, using UAS-GAL4DBD system, CHRO and dCTCF were artificially recruited into highly compacted polytene chromosome bands that share the features of silent chromatin type known as intercalary heterochromatin (IH). This led to local chromatin decondensation, formation of novel DHSes and recruitment of several "open chromatin" proteins. CHRO tethering resulted in the recruitment of CP190 and Z4 (PZG), whereas dCTCF tethering attracted CHRO, CP190, and Z4. Importantly, formation of a local stretch of open chromatin did not result in the reactivation of silent marker genes yellow and mini-white immediately adjacent to the targeting region (UAS), nor did RNA polII become recruited into this chromatin. The decompacted region retained late replicated, similarly to the wild-type untargeted region.

Genome-wide analysis of replication timing by next-generation sequencing with E/L Repli-seq

This protocol is an extension to: Nat. Protoc. 6, 870-895 (2014); doi:10.1038/nprot.2011.328; published online 02 June 2011Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early- and late-replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and subnuclear position. Moreover, RT is regulated during development and is altered in diseases. Here, we describe E/L Repli-seq, an extension of our Repli-chip protocol. E/L Repli-seq is a rapid, robust and relatively inexpensive protocol for analyzing RT by next-generation sequencing (NGS), allowing genome-wide assessment of how cellular processes are linked to RT. Briefly, cells are pulse-labeled with BrdU, and early and late S-phase fractions are sorted by flow cytometry. Labeled nascent DNA is immunoprecipitated from both fractions and sequenced. Data processing leads to a single bedGraph file containing the ratio of nascent DNA from early versus late S-phase fractions. The results are comparable to those of Repli-chip, with the additional benefits of genome-wide sequence information and an increased dynamic range. We also provide computational pipelines for downstream analyses, for parsing phased genomes using single-nucleotide polymorphisms (SNPs) to analyze RT allelic asynchrony, and for direct comparison to Repli-chip data. This protocol can be performed in up to 3 d before sequencing, and requires basic cellular and molecular biology skills, as well as a basic understanding of Unix and R.

An international meta-analysis confirms the association of BNC2 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a common spinal deformity with the prevalence of approximately 3%. We previously conducted a genome-wide association study (GWAS) using a Japanese cohort and identified a novel locus on chromosome 9p22.2. However, a replication study using multi-population cohorts has not been conducted. To confirm the association of 9p22.2 locus with AIS in multi-ethnic populations, we conducted international meta-analysis using eight cohorts. In total, we analyzed 8,756 cases and 27,822 controls. The analysis showed a convincing evidence of association between rs3904778 and AIS. Seven out of eight cohorts had significant P value, and remaining one cohort also had the same trend as the seven. The combined P was 3.28 × 10⁻¹⁸ (odds ratio = 1.19, 95% confidence interval = 1.14–1.24). In silico analyses suggested that BNC2 is the AIS susceptibility gene in this locus.

Computational methods for analyzing genome-wide chromosome conformation capture data

In all organisms, chromatin is packed to fulfil structural constraints and functional requirements. The hierarchical model of chromatin organization in the 3D nuclear space encompasses different topologies at diverse scale lengths, with chromosomes occupying distinct volumes, further organized in compartments, inside which the chromatin fibers fold into large domains and short-range loops. In the recent years, the combination of chromosome conformation capture (3C) techniques and high-throughput sequencing allowed probing chromatin spatial organization at the whole genome-scale. 3C-based methods produce enormous amounts of genomic data that are analyzed using ad-hoc computational procedures. Here, we review the common pipelines and methods for the analysis of genome-wide chromosome conformation capture data, highlighting recent developments in key steps for the identification of chromatin structures.

3C and 3C-based techniques: the powerful tools for spatial genome organization deciphering

It is well known that the chromosomes are organized in the nucleus and this spatial arrangement of genome play a crucial role in gene regulation and genome stability. Different techniques have been developed and applied to uncover the intrinsic mechanism of genome architecture, especially the chromosome conformation capture (3C) and 3C-derived methods. 3C and 3C-derived techniques provide us approaches to perform high-throughput chromatin architecture assays at the genome scale. However, the advantage and disadvantage of current methodologies of C-technologies have not been discussed extensively. In this review, we described and compared the methodologies of C-technologies used in genome organization studies with an emphasis on Hi-C method. We also discussed the crucial challenges facing current genome architecture studies based on 3C and 3C-derived technologies and the direction of future technologies to address currently outstanding questions in the field. These latest news contribute to our current understanding of genome structure, and provide a comprehensive reference for researchers to choose the appropriate method in future application. We consider that these constantly improving technologies will offer a finer and more accurate contact profiles of entire genome and ultimately reveal specific molecular machines govern its shape and function.

Genome-wide map of proximity linkage to renin proximal promoter in rat

A challenge to understanding enhancer-gene relationships is that enhancers are not always sequentially close to the gene they regulate. Physical proximity mapping through sequencing can provide an unbiased view of the chromatin close to the proximal promoter of the renin gene ( Ren). Our objective was to determine genomic regions that physically interact with the renin proximal promoter, using two different genetic backgrounds, the Dahl salt sensitive and normotensive SS-13BN, which have been shown to have different regulation of plasma renin in vivo. The chromatin conformation capture method with sequencing focused at the Ren proximal promoter in rat-derived cardiac endothelial cells was used. Cells were fixed, chromatin close to the Ren promoter was captured, and fragments were sequenced. The clustering of mapped reads produced a genome-wide map of chromatin in contact with the Ren promoter. The largest number of contacts was found on chromosome 13, the chromosome with Ren, and contacts were found on all other chromosomes except chromosome X. These contacts were significantly enriched with genes positively correlated with Ren expression and with mapped quantitative trait loci associated with blood pressure, cardiovascular, and renal phenotypes. The results were reproducible in an independent biological replicate. The findings reported here represent the first map between a critical cardiovascular gene and physical interacting loci throughout the genome and will provide the basis for several new directions of research.

Inter-chromosomal Contact Properties in Live-Cell Imaging and in Hi-C

Imaging (fluorescence in situ hybridization [FISH]) and genome-wide chromosome conformation capture (Hi-C) are two major approaches to the study of higher-order genome organization in the nucleus. Intra-chromosomal and inter-chromosomal interactions (referred to as non-homologous chromosomal contacts [NHCCs]) have been observed by several FISH-based studies, but locus-specific NHCCs have not been detected by Hi-C. Due to crosslinking, neither of these approaches assesses spatiotemporal properties. Toward resolving the discrepancies between imaging and Hi-C, we sought to understand the spatiotemporal properties of NHCCs in living cells by CRISPR/Cas9 live-cell imaging (CLING). In mammalian cells, we find that NHCCs are stable and occur as frequently as intra-chromosomal interactions, but NHCCs occur at farther spatial distance that could explain their lack of detection in Hi-C. By revealing the spatiotemporal properties in living cells, our study provides fundamental insights into the biology of NHCCs.

Dissecting super-enhancer hierarchy based on chromatin interactions

Recent studies have highlighted super-enhancers (SEs) as important regulatory elements for gene expression, but their intrinsic properties remain incompletely characterized. Through an integrative analysis of Hi-C and ChIP-seq data, here we find that a significant fraction of SEs are hierarchically organized, containing both hub and non-hub enhancers. Hub enhancers share similar histone marks with non-hub enhancers, but are distinctly associated with cohesin and CTCF binding sites and disease-associated genetic variants. Genetic ablation of hub enhancers results in profound defects in gene activation and local chromatin landscape. As such, hub enhancers are the major constituents responsible for SE functional and structural organization.

Regulation and Dysregulation of Chromosome Structure in Cancer

Cancer arises from genetic alterations that produce dysregulated gene expression programs. Normal gene regulation occurs in the context of chromosome loop structures called insulated neighborhoods, and recent studies have shown that these structures are altered and can contribute to oncogene dysregulation in various cancer cells. We review the types of genetic and epigenetic alterations that influence neighborhood structures and contribute to gene dysregulation in cancer, present models for insulated neighborhoods associated with the most prominent human oncogenes, and discuss how such models may lead to further advances in cancer diagnosis and therapy.

Precision medicine based on surgical oncology in the era of genome-scale analysis and genome editing technology

Accumulated evidence suggests that multiple molecular and cellular interactions promote cancer evolution in vivo. Surgical oncology is of growing significance to a comprehensive understanding of the malignant diseases for therapeutic application. We have analyzed more than 1000 clinical samples from surgically resected tissue to identify molecular biomarkers and therapeutic targets for advanced malignancies. Cancer stemness and mitotic instability were then determined as the essential predictors of aggressive phenotype with poor prognosis. Recently, whole genome/exome sequencing showed a mutational landscape underlying phenotype heterogeneity in caners. In addition, integrated genomic, epigenomic, transcriptomic, metabolic, proteomic and phenomic analyses elucidated several molecular subtypes that cluster in liver, pancreatic, biliary, esophageal and gastroenterological cancers. Identification of each molecular subtype is expected to realize the precise medicine targeting subtype-specific molecules; however, there are obstacle limitations to determine matching druggable targets or synthetic lethal interactions. Current breakthroughs in genome editing technology can provide us with unprecedented opportunity to recapitulate subtype-specific pathophysiology in vitro and in vivo. Given a great potential, on-demand editing system can design actionable strategy and revolutionize precision cancer medicine based on surgical oncology.

Chromatin organisation and cancer prognosis: a pan-cancer study

BACKGROUND

Chromatin organisation affects gene expression and regional mutation frequencies and contributes to carcinogenesis. Aberrant organisation of DNA has been correlated with cancer prognosis in analyses of the chromatin component of tumour cell nuclei using image texture analysis. As yet, the methodology has not been sufficiently validated to permit its clinical application. We aimed to define and validate a novel prognostic biomarker for the automatic detection of heterogeneous chromatin organisation.

METHODS

Machine learning algorithms analysed the chromatin organisation in 461 000 images of tumour cell nuclei stained for DNA from 390 patients (discovery cohort) treated for stage I or II colorectal cancer at the Aker University Hospital (Oslo, Norway). The resulting marker of chromatin heterogeneity, termed Nucleotyping, was subsequently independently validated in six patient cohorts: 442 patients with stage I or II colorectal cancer in the Gloucester Colorectal Cancer Study (UK); 391 patients with stage II colorectal cancer in the QUASAR 2 trial; 246 patients with stage I ovarian carcinoma; 354 patients with uterine sarcoma; 307 patients with prostate carcinoma; and 791 patients with endometrial carcinoma. The primary outcome was cancer-specific survival.

FINDINGS

In all patient cohorts, patients with chromatin heterogeneous tumours had worse cancer-specific survival than patients with chromatin homogeneous tumours (univariable analysis hazard ratio [HR] 1·7, 95% CI 1·2-2·5, in the discovery cohort; 1·8, 1·0-3·0, in the Gloucester validation cohort; 2·2, 1·1-4·5, in the QUASAR 2 validation cohort; 3·1, 1·9-5·0, in the ovarian carcinoma cohort; 2·5, 1·8-3·4, in the uterine sarcoma cohort; 2·3, 1·2-4·6, in the prostate carcinoma cohort; and 4·3, 2·8-6·8, in the endometrial carcinoma cohort). After adjusting for established prognostic patient characteristics in multivariable analyses, Nucleotyping was prognostic in all cohorts except for the prostate carcinoma cohort (HR 1·7, 95% CI 1·1-2·5, in the discovery cohort; 1·9, 1·1-3·2, in the Gloucester validation cohort; 2·6, 1·2-5·6, in the QUASAR 2 cohort; 1·8, 1·1-3·0, for ovarian carcinoma; 1·6, 1·0-2·4, for uterine sarcoma; 1·43, 0·68-2·99, for prostate carcinoma; and 1·9, 1·1-3·1, for endometrial carcinoma). Chromatin heterogeneity was a significant predictor of cancer-specific survival in microsatellite unstable (HR 2·9, 95% CI 1·0-8·4) and microsatellite stable (1·8, 1·2-2·7) stage II colorectal cancer, but microsatellite instability was not a significant predictor of outcome in chromatin homogeneous (1·3, 0·7-2·4) or chromatin heterogeneous (0·8, 0·3-2·0) stage II colorectal cancer.

INTERPRETATION

The consistent prognostic prediction of Nucleotyping in different biological and technical circumstances suggests that the marker of chromatin heterogeneity can be reliably assessed in routine clinical practice and could be used to objectively assist decision making in a range of clinical settings. An immediate application would be to identify high-risk patients with stage II colorectal cancer who might have greater absolute benefit from adjuvant chemotherapy. Clinical trials are warranted to evaluate the survival benefit and cost-effectiveness of using Nucleotyping to guide treatment decisions in multiple clinical settings.

FUNDING

The Research Council of Norway, the South-Eastern Norway Regional Health Authority, the National Institute for Health Research, and the Wellcome Trust.

FIND: difFerential chromatin INteractions Detection using a spatial Poisson process

Polymer-based simulations and experimental studies indicate the existence of a spatial dependency between the adjacent DNA fibers involved in the formation of chromatin loops. However, the existing strategies for detecting differential chromatin interactions assume that the interacting segments are spatially independent from the other segments nearby. To resolve this issue, we developed a new computational method, FIND, which considers the local spatial dependency between interacting loci. FIND uses a spatial Poisson process to detect differential chromatin interactions that show a significant difference in their interaction frequency and the interaction frequency of their neighbors. Simulation and biological data analysis show that FIND outperforms the widely used count-based methods and has a better signal-to-noise ratio.

Impact of broad regulatory regions on Gdf5 expression and function in knee development and susceptibility to osteoarthritis

OBJECTIVES

Given the role of growth and differentiation factor 5 (GDF5) in knee development and osteoarthritis risk, we sought to characterise knee defects resulting from Gdf5 loss of function and how its regulatory regions control knee formation and morphology.

METHODS

The brachypodism (bp) mouse line, which harbours an inactivating mutation in Gdf5, was used to survey how Gdf5 loss of function impacts knee morphology, while two transgenic Gdf5 reporter bacterial artificial chromosome mouse lines were used to assess the spatiotemporal activity and function of Gdf5 regulatory sequences in the context of clinically relevant knee anatomical features.

RESULTS

Knees from homozygous bp mice (bp/bp) exhibit underdeveloped femoral condyles and tibial plateaus, no cruciate ligaments, and poorly developed menisci. Secondary ossification is also delayed in the distal femur and proximal tibia. bp/bp mice have significantly narrower femoral condyles, femoral notches and tibial plateaus, and curvier medial femoral condyles, shallower trochlea, steeper lateral tibial slopes and smaller tibial spines. Regulatory sequences upstream from Gdf5 were weakly active in the prenatal knee, while downstream regulatory sequences were active throughout life. Importantly, downstream but not upstream Gdf5 regulatory sequences fully restored all the key morphological features disrupted in the bp/bp mice.

CONCLUSIONS

Knee morphology is profoundly affected by Gdf5 absence, and downstream regulatory sequences mediate its effects by controlling Gdf5 expression in knee tissues. This downstream region contains numerous enhancers harbouring human variants that span the osteoarthritis association interval. We posit that subtle alterations to morphology driven by changes in downstream regulatory sequence underlie this locus' role in osteoarthritis risk.

The long-range interaction map of ribosomal DNA arrays

The repeated rDNA array gives rise to the nucleolus, an organelle that is central to cellular processes as varied as stress response, cell cycle regulation, RNA modification, cell metabolism, and genome stability. The rDNA array is also responsible for the production of more than 70% of all cellular RNAs (the ribosomal RNAs). The rRNAs are produced from two sets of loci: the 5S rDNA array resides exclusively on human chromosome 1 while the 45S rDNA arrays reside on the short arm of five human acrocentric chromosomes. These critical genome elements have remained unassembled and have been excluded from all Hi-C analyses to date. Here we built the first high resolution map of 5S and 45S rDNA array contacts with the rest of the genome combining over 15 billion Hi-C reads from several experiments. The data enabled sufficiently high coverage to map rDNA-genome interactions with 1MB resolution and identify rDNA-gene contacts. The map showed that the 5S and 45S arrays display preferential contact at common sites along the genome but are not themselves sufficiently close to yield 5S-45S Hi-C contacts. Ribosomal DNA contacts are enriched in segments of closed, repressed, and late replicating chromatin, as well as CTCF binding sites. Finally, we identified functional categories whose dispersed genes coalesced in proximity to the rDNA arrays or instead avoided proximity with the rDNA arrays. The observations further our understanding of the spatial localization of rDNA arrays and their contribution to the architecture of the cell nucleus.

Endothelial cell differentiation is encompassed by changes in long range interactions between inactive chromatin regions

Endothelial cells (ECs) differentiate from mesodermal progenitors during vasculogenesis. By comparing changes in chromatin interactions between human umbilical vein ECs, embryonic stem cells and mesendoderm cells, we identified regions exhibiting EC-specific compartmentalization and changes in the degree of connectivity within topologically associated domains (TADs). These regions were characterized by EC-specific transcription, binding of lineage-determining transcription factors and cohesin. In addition, we identified 1200 EC-specific long-range interactions (LRIs) between TADs. Most of the LRIs were connected between regions enriched for H3K9me3 involving pericentromeric regions, suggesting their involvement in establishing compartmentalization of heterochromatin during differentiation. Second, we provide evidence that EC-specific LRIs correlate with changes in the hierarchy of chromatin aggregation. Despite these rearrangements, the majority of chromatin domains fall within a pre-established hierarchy conserved throughout differentiation. Finally, we investigated the effect of hypoxia on chromatin organization. Although hypoxia altered the expression of hundreds of genes, minimal effect on chromatin organization was seen. Nevertheless, 70% of hypoxia-inducible genes situated within a TAD bound by HIF1α suggesting that transcriptional responses to hypoxia largely depend on pre-existing chromatin organization. Collectively our results show that large structural rearrangements establish chromatin architecture required for functional endothelium and this architecture remains largely unchanged in response to hypoxia.

Into the Fourth Dimension: Dysregulation of Genome Architecture in Aging and Alzheimer's Disease

Alzheimer's disease (AD) is a progressive neurodegenerative disease characterized by synapse dysfunction and cognitive impairment. Understanding the development and progression of AD is challenging, as the disease is highly complex and multifactorial. Both environmental and genetic factors play a role in AD pathogenesis, highlighted by observations of complex DNA modifications at the single gene level, and by new evidence that also implicates changes in genome architecture in AD patients. The four-dimensional structure of chromatin in space and time is essential for context-dependent regulation of gene expression in post-mitotic neurons. Dysregulation of epigenetic processes have been observed in the aging brain and in patients with AD, though there is not yet agreement on the impact of these changes on transcription. New evidence shows that proteins involved in genome organization have altered expression and localization in the AD brain, suggesting that the genomic landscape may play a critical role in the development of AD. This review discusses the role of the chromatin organizers and epigenetic modifiers in post-mitotic cells, the aging brain, and in the development and progression of AD. How these new insights can be used to help determine disease risk and inform treatment strategies will also be discussed.

Enhancing Hi-C data resolution with deep convolutional neural network HiCPlus

Although Hi-C technology is one of the most popular tools for studying 3D genome organization, due to sequencing cost, the resolution of most Hi-C datasets are coarse and cannot be used to link distal regulatory elements to their target genes. Here we develop HiCPlus, a computational approach based on deep convolutional neural network, to infer high-resolution Hi-C interaction matrices from low-resolution Hi-C data. We demonstrate that HiCPlus can impute interaction matrices highly similar to the original ones, while only using 1/16 of the original sequencing reads. We show that the models learned from one cell type can be applied to make predictions in other cell or tissue types. Our work not only provides a computational framework to enhance Hi-C data resolution but also reveals features underlying the formation of 3D chromatin interactions.

Despite its popularity for measuring the spatial organization of mammalian genomes, the resolution of most Hi-C datasets is coarse due to sequencing cost. Here, Zhang et al. develop HiCPlus, a computational approach based on deep convolutional neural network, to infer high-resolution Hi-C interaction matrices from low-resolution Hi-C data.

Insights about genome function from spatial organization of the genome

Over the last 15 years, development of chromosome conformation capture (3C) and its subsequent high-throughput variants in conjunction with the fast development of sequencing technology has allowed investigators to generate large volumes of data giving insights into the spatial three-dimensional (3D) architecture of the genome. This huge data has been analyzed and validated using various statistical, mathematical, genomics, and biophysical tools in order to examine the chromosomal interaction patterns, understand the organization of the chromosome, and find out functional implications of the interactions. This review summarizes the data generated by several large-scale high-throughput chromosome conformation capture studies and the functional implications obtained from the data analyses. We also discuss emerging results on factors (both CCCTC binding factor (CTCF) related and CTCF independent) that could contribute to looping interactions.