How best to identify chromosomal interactions: a comparison of approaches

Using Functional Genomic Data in Monocytes/Macrophages and Genotyping to Nominate Disease-Driving Single Nucleotide Polymorphisms and Target Genes in Juvenile Idiopathic Arthritis

Introduction

GWAS have identified multiple regions that confer risk for juvenile idiopathic arthritis (JIA). However, identifying the single nucleotide polymorphisms (SNPs) that drive disease risk is impeded by the SNPs' that identify risk loci being in linkage disequilibrium (LD) with hundreds of other SNPs. Since the causal SNPs remain unknown, it is difficult to identify target genes and use genetic information to inform patient care. We used genotyping and functional data in primary human monocytes/macrophages to nominate disease-driving SNPs on JIA risk haplotypes and identify their likely target genes.

Methods

We identified JIA risk haplotypes using Immunochip data from Hinks et al (Nature Gen 2013) and the meta-analysis from McIntosh et al (Arthritis Rheum 2017). We used genotyping data from 3,939 children with JIA and 14,412 healthy controls to identify SNPs that: (1) were situated within open chromatin in multiple immune cell types and (2) were more common in children with JIA than the controls (p< 0.05). We intersected the chosen SNPs (n=846) with regions of bi-directional transcription initiation characteristic of non-coding regulatory regions detected using dREG to analyze GRO-seq data. Finally, we used MicroC data to identify gene promoters interacting with the regulatory regions harboring the candidate causal SNPs.

Results

We identified 190 SNPs that overlap with dREG peaks in monocytes and126 SNPs that overlap with dREG peaks in macrophages. Of these SNPs, 101 were situated within dREG peaks in both monocytes and macrophages, suggesting that these SNPs exert their effects independent of the cellular activation state. MicroC data in monocytes identified 20 genes/transcripts whose promoters interact with the enhancers harboring the SNPs of interest.

Conclusion

SNPs in JIA risk regions that are candidate causal variants can be further screened using functional data such as GRO-seq. This process identifies a finite number of candidate causal SNPs, the majority of which are likely to exert their biological effects independent of cellular activation state in monocytes. Three-dimensional chromatin data generated with MicroC identifies genes likely to be influenced by these SNPs. These studies demonstrate the importance of investigations into the role of innate immunity in JIA.

BRCAness, DNA gaps, and gain and loss of PARP inhibitor-induced synthetic lethality

Mutations in the tumor-suppressor genes BRCA1 and BRCA2 resulting in BRCA1/2 deficiency are frequently identified in breast, ovarian, prostate, pancreatic, and other cancers. Poly(ADP-ribose) polymerase (PARP) inhibitors (PARPis) selectively kill BRCA1/2-deficient cancer cells by inducing synthetic lethality, providing an effective biomarker-guided strategy for targeted cancer therapy. However, a substantial fraction of cancer patients carrying BRCA1/2 mutations do not respond to PARPis, and most patients develop resistance to PARPis over time, highlighting a major obstacle to PARPi therapy in the clinic. Recent studies have revealed that changes of specific functional defects of BRCA1/2-deficient cells, particularly their defects in suppressing and protecting single-stranded DNA gaps, contribute to the gain or loss of PARPi-induced synthetic lethality. These findings not only shed light on the mechanism of action of PARPis, but also lead to revised models that explain how PARPis selectively kill BRCA-deficient cancer cells. Furthermore, new mechanistic principles of PARPi sensitivity and resistance have emerged from these studies, generating potentially useful guidelines for predicting the PARPi response and design therapies for overcoming PARPi resistance. In this Review, we will discuss these recent studies and put them in context with the classic views of PARPi-induced synthetic lethality, aiming to stimulate the development of new therapeutic strategies to overcome PARPi resistance and improve PARPi therapy.

Regulating epithelial-mesenchymal plasticity from 3D genome organization

Epithelial-mesenchymal transition (EMT) is a dynamic process enabling polarized epithelial cells to acquire mesenchymal features implicated in development and carcinoma progression. As our understanding evolves, it is clear the reversible execution of EMT arises from complex epigenomic regulation involving histone modifications and 3-dimensional (3D) genome structural changes, leading to a cascade of transcriptional events. This review summarizes current knowledge on chromatin organization in EMT, with a focus on hierarchical structures of the 3D genome and chromatin accessibility changes.

Removing unwanted variation between samples in Hi-C experiments

Hi-C data are commonly normalized using single sample processing methods, with focus on comparisons between regions within a given contact map. Here, we aim to compare contact maps across different samples. We demonstrate that unwanted variation, of likely technical origin, is present in Hi-C data with replicates from different individuals, and that properties of this unwanted variation change across the contact map. We present band-wise normalization and batch correction, a method for normalization and batch correction of Hi-C data and show that it substantially improves comparisons across samples, including in a quantitative trait loci analysis as well as differential enrichment across cell types.

MMCT-Loop: a mix model-based pipeline for calling targeted 3D chromatin loops

Protein-specific Chromatin Conformation Capture (3C)-based technologies have become essential for identifying distal genomic interactions with critical roles in gene regulation. The standard techniques include Chromatin Interaction Analysis by Paired-End Tag (ChIA-PET), in situ Hi-C followed by chromatin immunoprecipitation (HiChIP) also known as PLAC-seq. To identify chromatin interactions from these data, a variety of computational methods have emerged. Although these state-of-art methods address many issues with loop calling, only few methods can fit different data types simultaneously, and the accuracy as well as the efficiency these approaches remains limited. Here we have generated a pipeline, MMCT-Loop, which ensures the accurate identification of strong loops as well as dynamic or weak loops through a mixed model. MMCT-Loop outperforms existing methods in accuracy, and the detected loops show higher activation functionality. To highlight the utility of MMCT-Loop, we applied it to conformational data derived from neural stem cell (NSCs) and uncovered several previously unidentified regulatory regions for key master regulators of stem cell identity. MMCT-Loop is an accurate and efficient loop caller for targeted conformation capture data, which supports raw data or pre-processed valid pairs as input, the output interactions are formatted and easily uploaded to a genome browser for visualization.

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. In this study, we examine the positioning and functional significance of human genes, grouped by their evolutionary age, within the 3D organization of the genome. We reveal that genes of different evolutionary origin have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, such associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young genes to ancient genes. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

Be-1DCNN: a neural network model for chromatin loop prediction based on bagging ensemble learning

The chromatin loops in the three-dimensional (3D) structure of chromosomes are essential for the regulation of gene expression. Despite the fact that high-throughput chromatin capture techniques can identify the 3D structure of chromosomes, chromatin loop detection utilizing biological experiments is arduous and time-consuming. Therefore, a computational method is required to detect chromatin loops. Deep neural networks can form complex representations of Hi-C data and provide the possibility of processing biological datasets. Therefore, we propose a bagging ensemble one-dimensional convolutional neural network (Be-1DCNN) to detect chromatin loops from genome-wide Hi-C maps. First, to obtain accurate and reliable chromatin loops in genome-wide contact maps, the bagging ensemble learning method is utilized to synthesize the prediction results of multiple 1DCNN models. Second, each 1DCNN model consists of three 1D convolutional layers for extracting high-dimensional features from input samples and one dense layer for producing the prediction results. Finally, the prediction results of Be-1DCNN are compared to those of the existing models. The experimental results indicate that Be-1DCNN predicts high-quality chromatin loops and outperforms the state-of-the-art methods using the same evaluation metrics. The source code of Be-1DCNN is available for free at https://github.com/HaoWuLab-Bioinformatics/Be1DCNN.

Ultrafine mapping of chromosome conformation at hundred basepair resolution reveals regulatory genome architecture

The resolution limit of chromatin conformation capture methodologies (3Cs) has restrained their application in detection of fine-level chromatin structure mediated by cis-regulatory elements (CREs). Here, we report two 3C-derived methods, Tri-4C and Tri-HiC, which utilize multirestriction enzyme digestions for ultrafine mapping of targeted and genome-wide chromatin interaction, respectively, at up to one hundred basepair resolution. Tri-4C identified CRE loop interaction networks and quantitatively revealed their alterations underlying dynamic gene control. Tri-HiC uncovered global fine-gauge regulatory interaction networks, identifying >20-fold more enhancer:promoter (E:P) loops than in situ Hi-C. In addition to vastly improved identification of subkilobase-sized E:P loops, Tri-HiC also uncovered interaction stripes and contact domain insulation from promoters and enhancers, revealing their loop extrusion behaviors resembling the topologically associating domain boundaries. Tri-4C and Tri-HiC provide robust approaches to achieve the high-resolution interactome maps required for characterizing fine-gauge regulatory chromatin interactions in analysis of development, homeostasis, and disease.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

Genome architecture and selective signals compensatorily shape plastic response to a new environment

Transcriptional plasticity interacts with natural selection in complex ways and is crucial for the survival of species under rapid climate change. How 3D genome architecture affects transcriptional plasticity and its interaction with genetic adaptation are unclear. We transplanted estuarine oysters to a new environment and found that genes located in active chromatin regions exhibited greater transcriptional plasticity, and changes in these regions were negatively correlated with selective signals. This indicates a trade-off between 3D active regions and selective signals in shaping plastic responses to a new environment. Specifically, a mutation, lincRNA, and changes in the accessibility of a distal enhancer potentially affect its interaction with the ManⅡa gene, which regulates the muscle function and survival of oysters. Our findings reveal that 3D genome architecture compensates for the role of genetic adaptation in environmental response to new environments and provide insights into synergetic genetic and epigenetic interactions critical for fitness-related trait and survival in a model marine species.

Determining chromatin architecture with Micro Capture-C

Micro Capture-C (MCC) is a chromatin conformation capture (3C) method for visualizing reproducible three-dimensional contacts of specified regions of the genome at base pair resolution. These methods are an established family of techniques that use proximity ligation to assay the topology of chromatin. MCC can generate data at substantially higher resolution than previous techniques through multiple refinements of the 3C method. Using a sequence agnostic nuclease, the maintenance of cellular integrity and full sequencing of the ligation junctions, MCC achieves subnucleosomal levels of resolution, which can be used to reveal transcription factor binding sites analogous to DNAse I footprinting. Gene dense regions, close-range enhancer-promoter contacts, individual enhancers within super-enhancers and multiple other types of loci or regulatory regions that were previously challenging to assay with conventional 3C techniques, are readily observed using MCC. MCC requires training in common molecular biology techniques and bioinformatics to perform the experiment and analyze the data. The protocol can be expected to be completed in a 3 week timeframe for experienced molecular biologists.

Endogenous neural stem cells characterization using omics approaches: Current knowledge in health and disease

Neural stem cells (NSCs), an invaluable source of neuronal and glial progeny, have been widely interrogated in the last twenty years, mainly to understand their therapeutic potential. Most of the studies were performed with cells derived from pluripotent stem cells of either rodents or humans, and have mainly focused on their potential in regenerative medicine. High-throughput omics technologies, such as transcriptomics, epigenetics, proteomics, and metabolomics, which exploded in the past decade, represent a powerful tool to investigate the molecular mechanisms characterizing the heterogeneity of endogenous NSCs. The transition from bulk studies to single cell approaches brought significant insights by revealing complex system phenotypes, from the molecular to the organism level. Here, we will discuss the current literature that has been greatly enriched in the “omics era”, successfully exploring the nature and function of endogenous NSCs and the process of neurogenesis. Overall, the information obtained from omics studies of endogenous NSCs provides a sharper picture of NSCs function during neurodevelopment in healthy and in perturbed environments.

ATM suppresses c-Myc overexpression in the mammary epithelium in response to estrogen

ATM gene mutation carriers are predisposed to estrogen-receptor-positive breast cancer (BC). ATM prevents BC oncogenesis by activating p53 in every cell; however, much remains unknown about tissue-specific oncogenesis after ATM loss. Here, we report that ATM controls the early transcriptional response to estrogens. This response depends on topoisomerase II (TOP2), which generates TOP2-DNA double-strand break (DSB) complexes and rejoins the breaks. When TOP2-mediated ligation fails, ATM facilitates DSB repair. After estrogen exposure, TOP2-dependent DSBs arise at the c-MYC enhancer in human BC cells, and their defective repair changes the activation profile of enhancers and induces the overexpression of many genes, including the c-MYC oncogene. CRISPR/Cas9 cleavage at the enhancer also causes c-MYC overexpression, indicating that this DSB causes c-MYC overexpression. Estrogen treatment induced c-Myc protein overexpression in mammary epithelial cells of ATM-deficient mice. In conclusion, ATM suppresses the c-Myc-driven proliferative effects of estrogens, possibly explaining such tissue-specific oncogenesis.

From structure to molecular condensates: emerging mechanisms for Mediator function

Mediator is a large modular protein assembly whose function as a coactivator of transcription is conserved in all eukaryotes. The Mediator complex can integrate and relay signals from gene-specific activators bound at enhancers to activate the general transcription machinery located at promoters. It has thus been described as a bridge between these elements during initiation of transcription. Here, we review recent studies on Mediator relating to its structure, gene specificity and general requirement, roles in chromatin architecture as well as novel concepts involving phase separation and transcriptional bursting. We revisit the mechanism of action of Mediator and ultimately put forward models for its mode of action in gene activation.

Identifying differential regulatory control of APOE ɛ4 on African versus European haplotypes as potential therapeutic targets

We previously demonstrated that in Alzheimer's disease (AD) patients, European apolipoprotein E (APOE) ε4 carriers express significantly more APOE ε4 in their brains than African AD carriers. We examined single nucleotide polymorphisms near APOE with significant frequency differences between African and European/Japanese APOE ε4 haplotypes that could contribute to this difference in expression through regulation. Two enhancer massively parallel reporter assay (MPRA) approaches were performed, supplemented with single fragment reporter assays. We used Capture C analyses to support interactions with the APOE promoter. Introns within TOMM40 showed increased enhancer activity in the European/Japanese versus African haplotypes in astrocytes and microglia. This region overlaps with APOE promoter interactions as assessed by Capture C analysis. Single variant analyses pinpoints rs2075650/rs157581, and rs59007384 as functionally different on these haplotypes. Identification of the mechanisms for differential regulatory function for APOE expression between African and European/Japanese haplotypes could lead to therapeutic targets for APOE ε4 carriers.

Engineered model of t(7;12)(q36;p13) AML recapitulates patient-specific features and gene expression profiles

Acute myeloid leukaemia carrying the translocation t(7;12)(q36;p13) is an adverse-risk leukaemia uniquely observed in infants. Despite constituting up to 30% of cases in under 2-year-olds, it remains poorly understood. Known molecular features are ectopic overexpression of the MNX1 gene and generation of a fusion transcript in 50% of patients. Lack of research models has hindered understanding of t(7;12) biology, which has historically focused on MNX1 overexpression rather than the cytogenetic entity itself. Here, we employed CRISPR/Cas9 to generate t(7;12) in the human K562 cell line, and in healthy CD34+ haematopoietic progenitors where the translocation was not sustained in long-term cultures or through serial replating. In contrast, in K562 cells, t(7;12) was propagated in self-renewing clonogenic assays, with sustained myeloid bias in colony formation and baseline depletion of erythroid signatures. Nuclear localisation analysis revealed repositioning of the translocated MNX1 locus to the interior of t(7;12)-harbouring K562 nuclei — a known phenomenon in t(7;12) patients which associates with ectopic overexpression of MNX1. Crucially, the K562-t(7;12) model successfully recapitulated the transcriptional landscape of t(7;12) patient leukaemia. In summary, we engineered a clinically-relevant model of t(7;12) acute myeloid leukaemia with the potential to unravel targetable molecular mechanisms of disease.

EPI-Mind: Identifying Enhancer-Promoter Interactions Based on Transformer Mechanism

MOTIVATION

Enhancer-Promoter Interactions (EPIs) is an essential step in the gene regulation process. However, the detection of EPIs by traditional wet experimental techniques is time-consuming and expensive. Thus, computational methods would be very useful for understanding the mechanism of EPIs. A number of approaches have been proposed to address this problem. Nevertheless, there is room for exploration and improvement for the existing research methods.

METHODS

In this study, a novel deep-learning model named EPI-Mind was proposed to predict EPIs with sequences features. First, we encoded enhancers and promoters sequences with pre-trained DNA vectors. Then, the Convolutional Neural Network (CNN) was utilized to rough extract the global and local features. Finally, the transformer mechanism was introduced to further extract the feature. We first trained a model named EPI-Mind_spe which can predict EPIs in one cell line. To achieve general prediction across different cell lines and further improve the performance of the model, a second-time training was carried on. The redivided dataset were used to train a new model called EPI-Mind_gen which can predict EPIs across different cell lines. To further improve the accuracy, a new model named EPI-Mind_best was trained which used the EPI-Mind_gen as a pre-trained model.

RESULTS

EPI-Mind_spe has the ability of predict EPIs with average AUROC above 90% and average AUPR above 70% but with cell lines specificity. EPI-Mind_gen can predict EPIs across different cell lines and its average AUROC is higher than the EPI-Mind_spe about 4.8%. EPI-Mind_best is superior to the state-of-the-art predictors on benchmarking datasets. EPI-Mind_best achieved best in 5 indicators within 12 indicators consists of AUPR and AUROC which is better than pioneers.

CONCLUSION

This research proposed a method, which was called EPI-Mind, to predict EPIs only with enhancer and promoters sequences, the framework of which was based on deep learning. This manuscript may provide a new route to solve the problem.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

The 3D genome landscape: Diverse chromosomal interactions and their functional implications

Genome organization includes contacts both within a single chromosome and between distinct chromosomes. Thus, regulatory organization in the nucleus may include interplay of these two types of chromosomal interactions with genome activity. Emerging advances in omics and single-cell imaging technologies have allowed new insights into chromosomal contacts, including those of homologs and sister chromatids, and their significance to genome function. In this review, we highlight recent studies in this field and discuss their impact on understanding the principles of chromosome organization and associated functional implications in diverse cellular processes. Specifically, we describe the contributions of intra-chromosomal, inter-homolog, and inter-sister chromatid contacts to genome organization and gene expression.

Multi-layered transcriptional control of cranial neural crest development

The neural crest (NC) is an emblematic population of embryonic stem-like cells with remarkable migratory ability. These distinctive attributes have inspired the curiosity of developmental biologists for over 150 years, however only recently the regulatory mechanisms controlling the complex features of the NC have started to become elucidated at genomic scales. Regulatory control of NC development is achieved through combinatorial transcription factor binding and recruitment of associated transcriptional complexes to distal cis-regulatory elements. Together, they regulate when, where and to what extent transcriptional programmes are actively deployed, ultimately shaping ontogenetic processes. Here, we discuss how transcriptional networks control NC ontogeny, with a special emphasis on the molecular mechanisms underlying specification of the cephalic NC. We also cover emerging properties of transcriptional regulation revealed in diverse developmental systems, such as the role of three-dimensional conformation of chromatin, and how they are involved in the regulation of NC ontogeny. Finally, we highlight how advances in deciphering the NC transcriptional network have afforded new insights into the molecular basis of human diseases.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

Chromatin interaction maps identify Wnt responsive cis-regulatory elements coordinating Paupar-Pax6 expression in neuronal cells

Central nervous system-expressed long non-coding RNAs (lncRNAs) are often located in the genome close to protein coding genes involved in transcriptional control. Such lncRNA-protein coding gene pairs are frequently temporally and spatially co-expressed in the nervous system and are predicted to act together to regulate neuronal development and function. Although some of these lncRNAs also bind and modulate the activity of the encoded transcription factors, the regulatory mechanisms controlling co-expression of neighbouring lncRNA-protein coding genes remain unclear. Here, we used high resolution NG Capture-C to map the cis-regulatory interaction landscape of the key neuro-developmental Paupar-Pax6 lncRNA-mRNA locus. The results define chromatin architecture changes associated with high Paupar-Pax6 expression in neurons and identify both promoter selective as well as shared cis-regulatory-promoter interactions involved in regulating Paupar-Pax6 co-expression. We discovered that the TCF7L2 transcription factor, a regulator of chromatin architecture and major effector of the Wnt signalling pathway, binds to a subset of these candidate cis-regulatory elements to coordinate Paupar and Pax6 co-expression. We describe distinct roles for Paupar in Pax6 expression control and show that the Paupar DNA locus contains a TCF7L2 bound transcriptional silencer whilst the Paupar transcript can act as an activator of Pax6. Our work provides important insights into the chromatin interactions, signalling pathways and transcription factors controlling co-expression of adjacent lncRNAs and protein coding genes in the brain.

CTCF: a misguided jack-of-all-trades in cancer cells

The emergence and progression of cancers is accompanied by a dysregulation of transcriptional programs. The three-dimensional (3D) organization of the human genome has emerged as an important multi-level mediator of gene transcription and regulation. In cancer cells, this organization can be restructured, providing a framework for the deregulation of gene activity. The CTCF protein, initially identified as the product from a tumor suppressor gene, is a jack-of-all-trades for the formation of 3D genome organization in normal cells. Here, we summarize how CTCF is involved in the multi-level organization of the human genome and we discuss emerging insights into how perturbed CTCF function and DNA binding causes the activation of oncogenes in cancer cells, mostly through a process of enhancer hijacking. Moreover, we highlight non-canonical functions of CTCF that can be relevant for the emergence of cancers as well. Finally, we provide guidelines for the computational identification of perturbed CTCF binding and reorganized 3D genome structure in cancer cells.

Analysis of sub-kilobase chromatin topology reveals nano-scale regulatory interactions with variable dependence on cohesin and CTCF

Enhancers and promoters predominantly interact within large-scale topologically associating domains (TADs), which are formed by loop extrusion mediated by cohesin and CTCF. However, it is unclear whether complex chromatin structures exist at sub-kilobase-scale and to what extent fine-scale regulatory interactions depend on loop extrusion. To address these questions, we present an MNase-based chromosome conformation capture (3C) approach, which has enabled us to generate the most detailed local interaction data to date (20 bp resolution) and precisely investigate the effects of cohesin and CTCF depletion on chromatin architecture. Our data reveal that cis-regulatory elements have distinct internal nano-scale structures, within which local insulation is dependent on CTCF, but which are independent of cohesin. In contrast, we find that depletion of cohesin causes a subtle reduction in longer-range enhancer-promoter interactions and that CTCF depletion can cause rewiring of regulatory contacts. Together, our data show that loop extrusion is not essential for enhancer-promoter interactions, but contributes to their robustness and specificity and to precise regulation of gene expression.

Super-enhancers and novel therapeutic targets in colorectal cancer

Transcription factors, cofactors, chromatin regulators, and transcription apparatuses interact with transcriptional regulatory elements, including promoters, enhancers, and super-enhancers (SEs), to coordinately regulate the transcription of target genes and thereby control cell behaviors. Among these transcriptional regulatory components and related elements, SEs often play a central role in determining cell identity and tumor initiation and progression. Therefore, oncogenic SEs, which are generated within cancer cells in oncogenes and other genes important in tumor pathogenesis, have emerged as attractive targets for novel cancer therapeutic strategies in recent years. Herein, we review the identification, formation and activation modes, and regulatory mechanisms for downstream genes and pathways of oncogenic SEs. We also review the therapeutic strategies and compounds targeting oncogenic SEs in colorectal cancer and other malignancies.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.

Accessible Region Conformation Capture (ARC-C) gives high-resolution insights into genome architecture and regulation

Nuclear organization and chromatin interactions are important for genome function, yet determining chromatin connections at high resolution remains a major challenge. To address this, we developed Accessible Region Conformation Capture (ARC-C), which profiles interactions between regulatory elements genome-wide without a capture step. Applied to Caenorhabditis elegans, ARC-C identifies approximately 15,000 significant interactions between regulatory elements at 500-bp resolution. Of 105 TFs or chromatin regulators tested, we find that the binding sites of 60 are enriched for interacting with each other, making them candidates for mediating interactions. These include cohesin and condensin II. Applying ARC-C to a mutant of transcription factor BLMP-1 detected changes in interactions between its targets. ARC-C simultaneously profiles domain-level architecture, and we observe that C. elegans chromatin domains defined by either active or repressive modifications form topologically associating domains (TADs) that interact with A/B (active/inactive) compartment-like structure. Furthermore, we discover that inactive compartment interactions are dependent on H3K9 methylation. ARC-C is a powerful new tool to interrogate genome architecture and regulatory interactions at high resolution.

Bacon: a comprehensive computational benchmarking framework for evaluating targeted chromatin conformation capture-specific methodologies

Chromatin conformation capture (3C)-based technologies have enabled the accurate detection of topological genomic interactions, and the adoption of ChIP techniques to 3C-based protocols makes it possible to identify long-range interactions. To analyze these large and complex datasets, computational methods are undergoing rapid and expansive evolution. Thus, a thorough evaluation of these analytical pipelines is necessary to identify which commonly used algorithms and processing pipelines need to be improved. Here we present a comprehensive benchmark framework, Bacon, to evaluate the performance of several computational methods. Finally, we provide practical recommendations for users working with HiChIP and/or ChIA-PET analyses.

Detection of Allele-Specific 3D Chromatin Interactions Using High-Resolution In-Nucleus 4C-seq

Chromosome conformation capture techniques are a set of methods used to determine 3D genome organization through the capture and identification of physical contacts between pairs of genomic loci. Among them, 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing) allows for the identification and quantification of the sequences interacting with a preselected locus of interest. 4C-seq has been widely used in the literature, mainly to study chromatin loops between enhancers and promoters or between CTCF binding sites and to identify chromatin domain boundaries. As 3D-contacts may be established in an allele-specific manner, we describe an up-to-date allele-specific 4C-seq protocol, starting from the selection of allele-specific viewpoints to Illumina sequencing. This protocol has mainly been optimized for cultured mammalian cells, but can be adapted for other cell types with relatively minor changes in initial steps.

Non-coding regulatory elements: Potential roles in disease and the case of epilepsy

Non-coding DNA (ncDNA) refers to the portion of the genome that does not code for proteins and accounts for the greatest physical proportion of the human genome. ncDNA includes sequences that are transcribed into RNA molecules, such as ribosomal RNAs (rRNAs), microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and un-transcribed sequences that have regulatory functions, including gene promoters and enhancers. Variation in non-coding regions of the genome have an established role in human disease, with growing evidence from many areas, including several cancers, Parkinson's disease and autism. Here, we review the features and functions of the regulatory elements that are present in the non-coding genome and the role that these regions have in human disease. We then review the existing research in epilepsy and emphasise the potential value of further exploring non-coding regulatory elements in epilepsy. In addition, we outline the most widely used techniques for recognising regulatory elements throughout the genome, current methodologies for investigating variation and the main challenges associated with research in the field of non-coding DNA.

Major Reorganization of Chromosome Conformation During Muscle Development in Pig

The spatial organization of the genome in the nucleus plays a crucial role in eukaryotic cell functions, yet little is known about chromatin structure variations during late fetal development in mammals. We performed in situ high-throughput chromosome conformation capture (Hi-C) sequencing of DNA from muscle samples of pig fetuses at two late stages of gestation. Comparative analysis of the resulting Hi-C interaction matrices between both groups showed widespread differences of different types. First, we discovered a complex landscape of stable and group-specific Topologically Associating Domains (TADs). Investigating the nuclear partition of the chromatin into transcriptionally active and inactive compartments, we observed a genome-wide fragmentation of these compartments between 90 and 110 days of gestation. Also, we identified and characterized the distribution of differential cis- and trans-pairwise interactions. In particular, trans-interactions at chromosome extremities revealed a mechanism of telomere clustering further confirmed by 3D Fluorescence in situ Hybridization (FISH). Altogether, we report major variations of the three-dimensional genome conformation during muscle development in pig, involving several levels of chromatin remodeling and structural regulation.

From GWAS to Gene: Transcriptome-Wide Association Studies and Other Methods to Functionally Understand GWAS Discoveries

Since their inception, genome-wide association studies (GWAS) have identified more than a hundred thousand single nucleotide polymorphism (SNP) loci that are associated with various complex human diseases or traits. The majority of GWAS discoveries are located in non-coding regions of the human genome and have unknown functions. The valley between non-coding GWAS discoveries and downstream affected genes hinders the investigation of complex disease mechanism and the utilization of human genetics for the improvement of clinical care. Meanwhile, advances in high-throughput sequencing technologies reveal important genomic regulatory roles that non-coding regions play in the transcriptional activities of genes. In this review, we focus on data integrative bioinformatics methods that combine GWAS with functional genomics knowledge to identify genetically regulated genes. We categorize and describe two types of data integrative methods. First, we describe fine-mapping methods. Fine-mapping is an exploratory approach that calibrates likely causal variants underneath GWAS signals. Fine-mapping methods connect GWAS signals to potentially causal genes through statistical methods and/or functional annotations. Second, we discuss gene-prioritization methods. These are hypothesis generating approaches that evaluate whether genetic variants regulate genes via certain genetic regulatory mechanisms to influence complex traits, including colocalization, mendelian randomization, and the transcriptome-wide association study (TWAS). TWAS is a gene-based association approach that investigates associations between genetically regulated gene expression and complex diseases or traits. TWAS has gained popularity over the years due to its ability to reduce multiple testing burden in comparison to other variant-based analytic approaches. Multiple types of TWAS methods have been developed with varied methodological designs and biological hypotheses over the past 5 years. We dive into discussions of how TWAS methods differ in many aspects and the challenges that different TWAS methods face. Overall, TWAS is a powerful tool for identifying complex trait-associated genes. With the advent of single-cell sequencing, chromosome conformation capture, gene editing technologies, and multiplexing reporter assays, we are expecting a more comprehensive understanding of genomic regulation and genetically regulated genes underlying complex human diseases and traits in the future.

New Data on Organization and Spatial Localization of B-Chromosomes in Cell Nuclei of the Yellow-Necked Mouse Apodemus flavicollis

The gene composition, function and evolution of B-chromosomes (Bs) have been actively discussed in recent years. However, the additional genomic elements are still enigmatic. One of Bs mysteries is their spatial organization in the interphase nucleus. It is known that heterochromatic compartments are not randomly localized in a nucleus. The purpose of this work was to study the organization and three-dimensional spatial arrangement of Bs in the interphase nucleus. Using microdissection of Bs and autosome centromeric heterochromatic regions of the yellow-necked mouse (Apodemus flavicollis) we obtained DNA probes for further two-dimensional (2D)- and three-dimensional (3D)- fluorescence in situ hybridization (FISH) studies. Simultaneous in situ hybridization of obtained here B-specific DNA probes and autosomal C-positive pericentromeric region-specific probes further corroborated the previously stated hypothesis about the pseudoautosomal origin of the additional chromosomes of this species. Analysis of the spatial organization of the Bs demonstrated the peripheral location of B-specific chromatin within the interphase nucleus and feasible contact with the nuclear envelope (similarly to pericentromeric regions of autosomes and sex chromosomes). It is assumed that such interaction is essential for the regulation of nuclear architecture. It also points out that Bs may follow the same mechanism as sex chromosomes to avoid a meiotic checkpoint.

Defining genome architecture at base-pair resolution

In higher eukaryotes, many genes are regulated by enhancers that are 104-106 base pairs (bp) away from the promoter. Enhancers contain transcription-factor-binding sites (which are typically around 7-22 bp), and physical contact between the promoters and enhancers is thought to be required to modulate gene expression. Although chromatin architecture has been mapped extensively at resolutions of 1 kilobase and above; it has not been possible to define physical contacts at the scale of the proteins that determine gene expression. Here we define these interactions in detail using a chromosome conformation capture method (Micro-Capture-C) that enables the physical contacts between different classes of regulatory elements to be determined at base-pair resolution. We find that highly punctate contacts occur between enhancers, promoters and CCCTC-binding factor (CTCF) sites and we show that transcription factors have an important role in the maintenance of the contacts between enhancers and promoters. Our data show that interactions between CTCF sites are increased when active promoters and enhancers are located within the intervening chromatin. This supports a model in which chromatin loop extrusion1 is dependent on cohesin loading at active promoters and enhancers, which explains the formation of tissue-specific chromatin domains without changes in CTCF binding.

CSynth: an interactive modelling and visualization tool for 3D chromatin structure

MOTIVATION

The 3D structure of chromatin in the nucleus is important for gene expression and regulation. Chromosome conformation capture techniques, such as Hi-C, generate large amounts of data showing interaction points on the genome but these are hard to interpret using standard tools.

RESULTS

We have developed CSynth, an interactive 3D genome browser and real-time chromatin restraint-based modeller to visualize models of any chromosome conformation capture (3C) data. Unlike other modelling systems, CSynth allows dynamic interaction with the modelling parameters to allow experimentation and effects on the model. It also allows comparison of models generated from data in different tissues/cell states and the results of third-party 3D modelling outputs. In addition, we include an option to view and manipulate these complicated structures using Virtual Reality (VR) so scientists can immerse themselves in the models for further understanding. This VR component has also proven to be a valuable teaching and a public engagement tool.

AVAILABILITYAND IMPLEMENTATION

CSynth is web based and available to use at csynth.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A unified framework for inferring the multi-scale organization of chromatin domains from Hi-C

Chromosomes are giant chain molecules organized into an ensemble of three-dimensional structures characterized with its genomic state and the corresponding biological functions. Despite the strong cell-to-cell heterogeneity, the cell-type specific pattern demonstrated in high-throughput chromosome conformation capture (Hi-C) data hints at a valuable link between structure and function, which makes inference of chromatin domains (CDs) from the pattern of Hi-C a central problem in genome research. Here we present a unified method for analyzing Hi-C data to determine spatial organization of CDs over multiple genomic scales. By applying statistical physics-based clustering analysis to a polymer physics model of the chromosome, our method identifies the CDs that best represent the global pattern of correlation manifested in Hi-C. The multi-scale intra-chromosomal structures compared across different cell types uncover the principles underlying the multi-scale organization of chromatin chain: (i) Sub-TADs, TADs, and meta-TADs constitute a robust hierarchical structure. (ii) The assemblies of compartments and TAD-based domains are governed by different organizational principles. (iii) Sub-TADs are the common building blocks of chromosome architecture. Our physically principled interpretation and analysis of Hi-C not only offer an accurate and quantitative view of multi-scale chromatin organization but also help decipher its connections with genome function.

An array of square blocks and checkerboard patterns characteristic to Hi-C data reflects the multi-scale organization of the chromatin chain. Deciphering the structures of chromatin domains from Hi-C and associating them with genome function are open problems of great importance in genome research. However, most existing methods are specialized in finding domains at different scales, making it difficult to integrate the solutions. Here we develop a unified framework for modeling and inferring domain structures over multiple scales, based on a physical model of the chromosome that reflects its nature as a three-dimensional object. Our method efficiently explores the space of domain solutions at different genomic scales, and systematically infers the chromatin domains over multiple scales from Hi-C data by employing a single tuning parameter. Our principled interpretation of Hi-C not only offers a quantitative view of multi-scale chromatin organization but also helps understand its connections with genome function.

High-resolution targeted 3C interrogation of cis-regulatory element organization at genome-wide scale

Chromosome conformation capture (3C) provides an adaptable tool for studying diverse biological questions. Current 3C methods generally provide either low-resolution interaction profiles across the entire genome, or high-resolution interaction profiles at limited numbers of loci. Due to technical limitations, generation of reproducible high-resolution interaction profiles has not been achieved at genome-wide scale. Here, to overcome this barrier, we systematically test each step of 3C and report two improvements over current methods. We show that up to 30% of reporter events generated using the popular in situ 3C method arise from ligations between two individual nuclei, but this noise can be almost entirely eliminated by isolating intact nuclei after ligation. Using Nuclear-Titrated Capture-C, we generate reproducible high-resolution genome-wide 3C interaction profiles by targeting 8055 gene promoters in erythroid cells. By pairing high-resolution 3C interaction calls with nascent gene expression we interrogate the role of promoter hubs and super-enhancers in gene regulation.

BET inhibition disrupts transcription but retains enhancer-promoter contact

Enhancers are DNA sequences that enable complex temporal and tissue-specific regulation of genes in higher eukaryotes. Although it is not entirely clear how enhancer-promoter interactions can increase gene expression, this proximity has been observed in multiple systems at multiple loci and is thought to be essential for the maintenance of gene expression. Bromodomain and Extra-Terminal domain (BET) and Mediator proteins have been shown capable of forming phase condensates and are thought to be essential for super-enhancer function. Here, we show that targeting of cells with inhibitors of BET proteins or pharmacological degradation of BET protein Bromodomain-containing protein 4 (BRD4) has a strong impact on transcription but very little impact on enhancer-promoter interactions. Dissolving phase condensates reduces BRD4 and Mediator binding at enhancers and can also strongly affect gene transcription, without disrupting enhancer-promoter interactions. These results suggest that activation of transcription and maintenance of enhancer-promoter interactions are separable events. Our findings further indicate that enhancer-promoter interactions are not dependent on high levels of BRD4 and Mediator, and are likely maintained by a complex set of factors including additional activator complexes and, at some sites, CTCF and cohesin.

Folding Keratin Gene Clusters during Skin Regional Specification

Regional specification is critical for skin development, regeneration, and evolution. The contribution of epigenetics in this process remains unknown. Here, using avian epidermis, we find two major strategies regulate β-keratin gene clusters. (1) Over the body, macro-regional specificities (scales, feathers, claws, etc.) established by typical enhancers control five subclusters located within the epidermal differentiation complex on chromosome 25; (2) within a feather, micro-regional specificities are orchestrated by temporospatial chromatin looping of the feather β-keratin gene cluster on chromosome 27. Analyses suggest a three-factor model for regional specification: competence factors (e.g., AP1) make chromatin accessible, regional specifiers (e.g., Zic1) target specific genome regions, and chromatin regulators (e.g., CTCF and SATBs) establish looping configurations. Gene perturbations disrupt morphogenesis and histo-differentiation. This chicken skin paradigm advances our understanding of how regulation of big gene clusters can set up a two-dimensional body surface map.

MZF1 regulates α-globin gene transcription via long-range interactions in erythroid differentiation

Precise spatiotemporal gene expression regulation is crucial for human erythropoiesis. However, dramatic changes in the chromatin structure and transcriptome involved in α-globin gene expression during erythropoiesis still not fully understand. To identify candidate regulators for α-globin gene regulation, we carried out an integrated approach by integrating publicly available transcriptomic and epigenomic data. We computed active enhancers by overlapping enriched regions marked with H3K4me1 and H3K27ac and correlated their activity with mRNA expression. Next, we cataloged potential transcription factors via de novo motif analysis. We highlighted the discovery of potential novel transcription factor MZF1 of the α-globin gene in erythroid differentiation. To validate the role of MZF1, we quantified the expression level of MZF1 and α-globin gene in HSPCs, early erythroid progenitors and late erythroid precursors cells. Both the mRNA and protein expression patterns of MZF1 were consistent with the α-globin gene. Also, the qPCR result showed that the expression of the α-globin gene was significantly increased by the MZF1 overexpression. To further investigate the role of MZF1 regulating α-globin gene transcriptional activity during erythroid differentiation, we performed ChIP-qPCR at the α-globin locus. Our results showed that MZF1 recruitment both at 4 upstream HS sites and α-globin gene promoter in erythroid precursor cells. To determine the importance of the MZF1 to enhancer-promoter interaction at the α-globin locus, we compared interaction frequency before and after knockdown of MZF1 by chromosome conformation capture (3C) assay. Upon MZF1 depletion, both the expression of the α-globin gene and all 3C signals were significantly decreased. Taken together, MZF1 plays an important role in regulating α-globin gene expression by binding to long-region enhancers and α-globin gene promoter and facilitates the organization of specific 3D chromatin architecture in erythroid differentiation.

Reflections on the organization and the physical state of chromatin in eukaryotic cells

In recent years, our perception of chromatin structure and organization in the cell nucleus has changed in fundamental ways. The 30 nm chromatin fiber has lost its status as an essential in vivo structure. Hi-C and related biochemical methods, advanced electron and super-resolved fluorescence microscopy, together with concepts from soft matter physics, have revolutionized the field. A comprehensive understanding of the structural and functional interactions that regulate cell cycle and cell type specific nuclear functions appears within reach, but it requires the integration of top-down and bottom-up approachs. In this review, I present an update on nuclear architecture studies with an emphasis on organization and the controversy regarding the physical state of chromatin in cells.

The emerging role of super enhancer-derived noncoding RNAs in human cancer

Super enhancers (SEs) are large clusters of adjacent enhancers that drive the expression of genes which regulate cellular identity; SE regions can be enriched with a high density of transcription factors, co-factors, and enhancer-associated epigenetic modifications. Through enhanced activation of their target genes, SEs play an important role in various diseases and conditions, including cancer. Recent studies have shown that SEs not only activate the transcriptional expression of coding genes to directly regulate biological functions, but also drive the transcriptional expression of non-coding RNAs (ncRNAs) to indirectly regulate biological functions. SE-derived ncRNAs play critical roles in tumorigenesis, including malignant proliferation, metastasis, drug resistance, and inflammatory response. Moreover, the abnormal expression of SE-derived ncRNAs is closely related to the clinical and pathological characterization of tumors. In this review, we summarize the functions and roles of SE-derived ncRNAs in tumorigenesis and discuss their prospective applications in tumor therapy. A deeper understanding of the potential mechanism underlying the action of SE-derived ncRNAs in tumorigenesis may provide new strategies for the early diagnosis of tumors and targeted therapy.

A modified protocol of Capture-C allows affordable and flexible high-resolution promoter interactome analysis

Large-scale epigenomic projects have mapped hundreds of thousands of potential regulatory sites in the human genome, but only a small proportion of these elements are proximal to transcription start sites. It is believed that the majority of these sequences are remote promoter-activating genomic sites scattered within several hundreds of kilobases from their cognate promoters and referred to as enhancers. It is still unclear what principles, aside from relative closeness in the linear genome, determine which promoter(s) is controlled by a given enhancer; however, this understanding is of great fundamental and clinical relevance. In recent years, C-methods (chromosome conformation capture-based methods) have become a powerful tool for the identification of enhancer-promoter spatial contacts that, in most cases, reflect their functional link. Here, we describe a new hybridisation-based promoter Capture-C protocol that makes use of biotinylated dsDNA probes generated by PCR from a custom pool of long oligonucleotides. The described protocol allows high-resolution promoter interactome description, providing a flexible and cost-effective alternative to the existing promoter Capture-C modifications. Based on the obtained data, we propose several tips on probe design that could potentially improve the results of future experiments.

Advances in optical mapping for genomic research

Recent advances in optical mapping have allowed the construction of improved genome assemblies with greater contiguity. Optical mapping also enables genome comparison and identification of large-scale structural variations. Association of these large-scale genomic features with biological functions is an important goal in plant and animal breeding and in medical research. Optical mapping has also been used in microbiology and still plays an important role in strain typing and epidemiological studies. Here, we review the development of optical mapping in recent decades to illustrate its importance in genomic research. We detail its applications and algorithms to show its specific advantages. Finally, we discuss the challenges required to facilitate the optimization of optical mapping and improve its future development and application.

Cancer systems immunology

Tumor immunology is undergoing a renaissance due to the recent profound clinical successes of tumor immunotherapy. These advances have coincided with an exponential growth in the development of -omics technologies. Armed with these technologies and their associated computational and modeling toolsets, systems biologists have turned their attention to tumor immunology in an effort to understand the precise nature and consequences of interactions between tumors and the immune system. Such interactions are inherently multivariate, spanning multiple time and size scales, cell types, and organ systems, rendering systems biology approaches particularly amenable to their interrogation. While in its infancy, the field of 'Cancer Systems Immunology' has already influenced our understanding of tumor immunology and immunotherapy. As the field matures, studies will move beyond descriptive characterizations toward functional investigations of the emergent behavior that govern tumor-immune responses. Thus, Cancer Systems Immunology holds incredible promise to advance our ability to fight this disease.

Blood disease-causing and -suppressing transcriptional enhancers: general principles and GATA2 mechanisms

Intensive scrutiny of human genomes has unveiled considerable genetic variation in coding and noncoding regions. In cancers, including those of the hematopoietic system, genomic instability amplifies the complexity and functional consequences of variation. Although elucidating how variation impacts the protein-coding sequence is highly tractable, deciphering the functional consequences of variation in noncoding regions (genome reading), including potential transcriptional-regulatory sequences, remains challenging. A crux of this problem is the sheer abundance of gene-regulatory sequence motifs (cis elements) mediating protein-DNA interactions that are intermixed in the genome with thousands of look-alike sequences lacking the capacity to mediate functional interactions with proteins in vivo. Furthermore, transcriptional enhancers harbor clustered cis elements, and how altering a single cis element within a cluster impacts enhancer function is unpredictable. Strategies to discover functional enhancers have been innovated, and human genetics can provide vital clues to achieve this goal. Germline or acquired mutations in functionally critical (essential) enhancers, for example at the GATA2 locus encoding a master regulator of hematopoiesis, have been linked to human pathologies. Given the human interindividual genetic variation and complex genetic landscapes of hematologic malignancies, enhancer corruption, creation, and expropriation by new genes may not be exceedingly rare mechanisms underlying disease predisposition and etiology. Paradigms arising from dissecting essential enhancer mechanisms can guide genome-reading strategies to advance fundamental knowledge and precision medicine applications. In this review, we provide our perspective of general principles governing the function of blood disease-linked enhancers and GATA2-centric mechanisms.

Modern epigenetics methods in biological research

The definition of epigenetics refers that molecular modifications on DNA that can regulate gene activity are independent of DNA sequence and mitotically stable. Notably, epigenetics studies have grown exponentially in the past few years. Recent progresses that lead to exciting discoveries and groundbreaking nature of this area demand thorough methodologies and advanced technologies to move epigenetics to the forefront of molecular biology. The most recognized epigenetic regulations are DNA methylation, histone modifications, and non-coding RNAs (ncRNAs). This review will discuss the modern techniques that are available to detect locus-specific and genome-wide changes for all epigenetic codes. Furthermore, updated analysis of technologies, newly developed methods, recent breakthroughs and bioinformatics pipelines in epigenetic analysis will be presented. These methods, as well as many others presented in this specific issue, provide comprehensive guidelines in the area of epigenetics that facilitate further developments in this promising and rapidly developing field.

Are We Meeting the Promise of Endotypes and Precision Medicine in Asthma?

While the term asthma has long been known to describe heterogeneous groupings of patients, only recently have data evolved which enable a molecular understanding of the clinical differences. The evolution of transcriptomics (and other 'omics platforms) and improved statistical analyses in combination with large clinical cohorts opened the door for molecular characterization of pathobiologic processes associated with a range of asthma patients. When linked with data from animal models and clinical trials of targeted biologic therapies, emerging distinctions arose between patients with and without elevations in type 2 immune and inflammatory pathways, leading to the confirmation of a broad categorization of type 2-Hi asthma. Differences in the ratios, sources, and location of type 2 cytokines and their relation to additional immune pathway activation appear to distinguish several different (sub)molecular phenotypes, and perhaps endotypes of type 2-Hi asthma, which respond differently to broad and targeted anti-inflammatory therapies. Asthma in the absence of type 2 inflammation is much less well defined, without clear biomarkers, but is generally linked with poor responses to corticosteroids. Integration of "big data" from large cohorts, over time, using machine learning approaches, combined with validation and iterative learning in animal (and human) model systems is needed to identify the biomarkers and tightly defined molecular phenotypes/endotypes required to fulfill the promise of precision medicine.

Zooming in on chromosome dynamics

Until recently, our understanding of chromosome organization in higher eukaryotic cells has been based on analyses of large-scale, low-resolution changes in chromosomes structure. More recently, CRISPR-Cas9 technologies have allowed us to "zoom in" and visualize specific chromosome regions in live cells so that we can begin to examine in detail the dynamics of chromosome organization in individual cells. In this review, we discuss traditional methods of chromosome locus visualization and look at how CRISPR-Cas9 gene-targeting methodologies have helped improve their application. We also describe recent developments of the CRISPR-Cas9 technology that enable visualization of specific chromosome regions without the requirement for complex genetic manipulation.

Analysis of 4C-seq data: A comparison of methods

The circular chromosome conformation capture technique followed by sequencing (4C-seq) has been used in a number of studies to investigate chromosomal interactions between DNA fragments. Computational pipelines have been developed and published that offer various possibilities of 4C-seq data processing and statistical analysis. Here, we present an overview of four of such pipelines (fourSig, FourCSeq, 4C-ker and w4Cseq) taking into account the most important stages of computations. We provide comparisons of the methods and discuss their advantages and possible weaknesses. We illustrate the results with the use of data obtained for two different species, in a study devoted to vernalization control in Arabidopsis thaliana by the FLOWERING LOCUS C (FLC) gene and to long-range chromatin interactions in mouse embryonic stem cells.