Determination of High-Resolution 3D Chromatin Organization Using Circular Chromosome Conformation Capture (4C-seq)

JMJD1C forms condensates to facilitate a RUNX1-dependent gene expression program shared by multiple types of AML cells

JMJD1C, a member of the lysine demethylase 3 (KDM3) family, is universally required for the survival of several types of acute myeloid leukemia (AML) cells with different genetic mutations, representing a therapeutic opportunity with broad application. Yet how JMJD1C regulates the leukemic programs of various AML cells is largely unexplored. Here we show that JMJD1C interacts with the master hematopoietic transcription factor RUNX1, which thereby recruits JMJD1C to the genome to facilitate a RUNX1-driven transcriptional program that supports leukemic cell survival. The underlying mechanism hinges on the long N-terminal disordered region of JMJD1C, which harbors two inseparable abilities: condensate formation and direct interaction with RUNX1. This dual capability of JMJD1C may influence enhancer-promoter contacts crucial for the expression of key leukemic genes regulated by RUNX1. Our findings demonstrate a previously unappreciated role for the non-catalytic function of JMJD1C in transcriptional regulation, underlying a mechanism shared by different types of leukemias.

CircR-loop: a novel RNA:DNA interaction on genome instability

CircR-loop, a recently unearthed regulatory mechanism situated at the crossroads of circular RNA and DNA interactions, constitute a subset of R-loop. This circR-loop have emerged as a crucial player in pivotal regulatory functions within both animal and plant systems. The journey into the realm of circR-loop commenced with their discovery within the human mitochondrial genome, where they serve as critical directors of mitochondrial DNA replication. In the plant kingdom, circR-loop wield influence over processes such as alternative splicing and centromere organization, impacting the intricacies of floral development and genome stability, respectively. Their significance extends to the animal domain, where circR-loop has captured attention for their roles in cancer-related phenomena, exerting control over transcription, chromatin architecture, and orchestrating responses to DNA damage. Moreover, their involvement in nuclear export anomalies further underscores their prominence in cellular regulation. This article summarizes the important regulatory mechanisms and physiological roles of circR-loop in plants and animals, and offers a comprehensive exploration of the methodologies employed for the identification, characterization, and functional analysis of circR-loop, underscoring the pressing need for innovative approaches that can effectively distinguish them from their linear RNA counterparts while elucidating their precise functions. Lastly, the article sheds light on the challenges and opportunities that lie ahead in the field of circR-loop research, emphasizing the vital importance of continued investigations to uncover their regulatory roles and potential applications in the realm of biology. In summary, circR-loop represents a captivating and novel regulatory mechanism with broad-reaching implications spanning the realms of genetics, epigenetics, and disease biology. Their exploration opens new avenues for comprehending gene regulation and holds significant promise for future therapeutic interventions.

ARID1A regulates DNA repair through chromatin organization and its deficiency triggers DNA damage-mediated anti-tumor immune response

AT-rich interaction domain protein 1A (ARID1A), a SWI/SNF chromatin remodeling complex subunit, is frequently mutated across various cancer entities. Loss of ARID1A leads to DNA repair defects. Here, we show that ARID1A plays epigenetic roles to promote both DNA double-strand breaks (DSBs) repair pathways, non-homologous end-joining (NHEJ) and homologous recombination (HR). ARID1A is accumulated at DSBs after DNA damage and regulates chromatin loops formation by recruiting RAD21 and CTCF to DSBs. Simultaneously, ARID1A facilitates transcription silencing at DSBs in transcriptionally active chromatin by recruiting HDAC1 and RSF1 to control the distribution of activating histone marks, chromatin accessibility, and eviction of RNAPII. ARID1A depletion resulted in enhanced accumulation of micronuclei, activation of cGAS-STING pathway, and an increased expression of immunomodulatory cytokines upon ionizing radiation. Furthermore, low ARID1A expression in cancer patients receiving radiotherapy was associated with higher infiltration of several immune cells. The high mutation rate of ARID1A in various cancer types highlights its clinical relevance as a promising biomarker that correlates with the level of immune regulatory cytokines and estimates the levels of tumor-infiltrating immune cells, which can predict the response to the combination of radio- and immunotherapy.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci

BACKGROUND

Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space.

RESULTS

To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4.

CONCLUSIONS

Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

Multi-feature clustering of CTCF binding creates robustness for loop extrusion blocking and Topologically Associating Domain boundaries

Topologically Associating Domains (TADs) separate vertebrate genomes into insulated regulatory neighborhoods that focus genome-associated processes. TADs are formed by Cohesin-mediated loop extrusion, with many TAD boundaries consisting of clustered binding sites of the CTCF insulator protein. Here we determine how this clustering of CTCF binding contributes to the blocking of loop extrusion and the insulation between TADs. We identify enrichment of three features of CTCF binding at strong TAD boundaries, consisting of strongly bound and closely spaced CTCF binding peaks, with a further enrichment of DNA-binding motifs within these peaks. Using multi-contact Nano-C analysis in cells with normal and perturbed CTCF binding, we establish that individual CTCF binding sites contribute to the blocking of loop extrusion, but in an incomplete manner. When clustered, individual CTCF binding sites thus create a stepwise insulation between neighboring TADs. Based on these results, we propose a model whereby multiple instances of temporal loop extrusion blocking create strong insulation between TADs.

The Drosophila Fab-7 boundary modulates Abd-B gene activity by guiding an inversion of collinear chromatin organization and alternate promoter use

Hox genes encode transcription factors that specify segmental identities along the anteroposterior body axis. These genes are organized in clusters, where their order corresponds to their activity along the body axis, a feature known as collinearity. In Drosophila, the BX-C cluster contains the three most posterior Hox genes, where their collinear activation incorporates progressive changes in histone modifications, chromatin architecture, and use of boundary elements and cis-regulatory regions. To dissect functional hierarchies, we compare chromatin organization in cell lines and larvae, with a focus on the Abd-B gene. Our work establishes the importance of the Fab-7 boundary for insulation between 3D domains carrying different histone modifications. Interestingly, we detect a non-canonical inversion of collinear chromatin dynamics at Abd-B, with the domain of active histone modifications progressively decreasing in size. This dynamic chromatin organization differentially activates the alternative promoters of the Abd-B gene, thereby expanding the possibilities for fine-tuning of transcriptional output.

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

Promoter and enhancer RNAs regulate chromatin reorganization and activation of miR-10b/HOXD locus, and neoplastic transformation in glioma

miR-10b is silenced in normal neuroglial cells of the brain but commonly activated in glioma, where it assumes an essential tumor-promoting role. We demonstrate that the entire miR-10b-hosting HOXD locus is activated in glioma via the cis-acting mechanism involving 3D chromatin reorganization and CTCF-cohesin-mediated looping. This mechanism requires two interacting lncRNAs, HOXD-AS2 and LINC01116, one associated with HOXD3/HOXD4/miR-10b promoter and another with the remote enhancer. Knockdown of either lncRNA in glioma cells alters CTCF and cohesin binding, abolishes chromatin looping, inhibits the expression of all genes within HOXD locus, and leads to glioma cell death. Conversely, in cortical astrocytes, enhancer activation is sufficient for HOXD/miR-10b locus reorganization, gene derepression, and neoplastic cell transformation. LINC01116 RNA is essential for this process. Our results demonstrate the interplay of two lncRNAs in the chromatin folding and concordant regulation of miR-10b and multiple HOXD genes normally silenced in astrocytes and triggering the neoplastic glial transformation.

Non-coding RNA LEVER sequestration of PRC2 can mediate long range gene regulation

Polycomb Repressive Complex 2 (PRC2) is an epigenetic regulator required for gene silencing during development. Although PRC2 is a well-established RNA-binding complex, the biological function of PRC2-RNA interaction has been controversial. Here, we study the gene-regulatory role of the inhibitory PRC2-RNA interactions. We report a nuclear long non-coding RNA, LEVER, which mapped 236 kb upstream of the β-globin cluster as confirmed by Nanopore sequencing. LEVER RNA interacts with PRC2 in its nascent form, and this prevents the accumulation of the H3K27 repressive histone marks within LEVER locus. Interestingly, the accessible LEVER chromatin, in turn, suppresses the chromatin interactions between the ε-globin locus and β-globin locus control region (LCR), resulting in a repressive effect on ε-globin gene expression. Our findings validate that the nascent RNA-PRC2 interaction inhibits local PRC2 function in situ. More importantly, we demonstrate that such a local process can in turn regulate the expression of neighboring genes.

Identification of a long non-coding RNA LEVER, that inhibits the Polycomb Repressive Complex 2 (PRC2) and controls nearby embryonic form of beta-globin gene, provides additional evidence for PRC2-RNA functional interaction.

A fast Myosin super enhancer dictates muscle fiber phenotype through competitive interactions with Myosin genes

The contractile properties of adult myofibers are shaped by their Myosin heavy chain isoform content. Here, we identify by snATAC-seq a 42 kb super-enhancer at the locus regrouping the fast Myosin genes. By 4C-seq we show that active fast Myosin promoters interact with this super-enhancer by DNA looping, leading to the activation of a single promoter per nucleus. A rainbow mouse transgenic model of the locus including the super-enhancer recapitulates the endogenous spatio-temporal expression of adult fast Myosin genes. In situ deletion of the super-enhancer by CRISPR/Cas9 editing demonstrates its major role in the control of associated fast Myosin genes, and deletion of two fast Myosin genes at the locus reveals an active competition of the promoters for the shared super-enhancer. Last, by disrupting the organization of fast Myosin, we uncover positional heterogeneity within limb skeletal muscles that may underlie selective muscle susceptibility to damage in certain myopathies.

The contractile properties of adult myofibers are shaped by their Myosin heavy chain isoform content. Here the authors show that a super enhancer controls the spatiotemporal expression of the genes at the fast myosin heavy chain locus by DNA looping and that this expression profile is recapitulated in a rainbow transgenic mouse model of the locus.

Screening of ETO2-GLIS2-induced Super Enhancers identifies targetable cooperative dependencies in acute megakaryoblastic leukemia

Super Enhancers (SEs) are clusters of regulatory elements associated with cell identity and disease. However, whether these elements are induced by oncogenes and can regulate gene modules cooperating for cancer cell transformation or maintenance remains elusive. To address this question, we conducted a genome-wide CRISPRi-based screening of SEs in ETO2-GLIS2⁺ acute megakaryoblastic leukemia. This approach revealed SEs essential for leukemic cell growth and survival that are induced by ETO2-GLIS2 expression. In particular, we identified a de novo SE specific of this leukemia subtype and regulating expression of tyrosine kinase-associated receptors KIT and PDGFRA. Combined expression of these two receptors was required for leukemic cell growth, and CRISPRi-mediated inhibition of this SE or treatment with tyrosine kinase inhibitors impaired progression of leukemia in vivo in patient-derived xenografts experiments. Our results show that fusion oncogenes, such as ETO2-GLIS2, can induce activation of SEs regulating essential gene modules synergizing for leukemia progression.

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.

Chromatin Accessibility of Human Mitral Valves and Functional Assessment of MVP Risk Loci

RATIONALE

Mitral valve prolapse (MVP) is a common valvopathy that leads to mitral insufficiency, heart failure, and sudden death. Functional genomic studies in mitral valves are needed to better characterize MVP-associated variants and target genes.

OBJECTIVE

To establish the chromatin accessibility profiles and assess functionality of variants and narrow down target genes at MVP loci.

METHODS AND RESULTS

We mapped the open chromatin regions in nuclei from 11 human pathogenic and 7 nonpathogenic mitral valves by an assay for transposase-accessible chromatin with high-throughput sequencing. Open chromatin peaks were globally similar between pathogenic and nonpathogenic valves. Compared with the heart tissue and cardiac fibroblasts, we found that MV-specific assay for transposase-accessible chromatin with high-throughput sequencing peaks are enriched near genes involved in extracellular matrix organization, chondrocyte differentiation, and connective tissue development. One of the most enriched motifs in MV-specific open chromatin peaks was for the nuclear factor of activated T cells family of TFs (transcription factors) involved in valve endocardial and interstitial cell formation. We also found that MVP-associated variants were significantly enriched (P<0.05) in mitral valve open chromatin peaks. Integration of the assay for transposase-accessible chromatin with high-throughput sequencing data with risk loci, extensive functional annotation, and gene reporter assay suggest plausible causal variants for rs2641440 at the SMG6/SRR locus and rs6723013 at the IGFBP2/IGFBP5/TNS1 locus. CRISPR-Cas9 deletion of the sequence including rs6723013 in human fibroblasts correlated with increased expression only for TNS1. Circular chromatin conformation capture followed by high-throughput sequencing experiments provided evidence for several target genes, including SRR, HIC1, and DPH1 at the SMG6/SRR locus and further supported TNS1 as the most likely target gene on chromosome 2.

CONCLUSIONS

Here, we describe unprecedented genome-wide open chromatin profiles from human pathogenic and nonpathogenic MVs and report specific gene regulation profiles, compared with the heart. We also report in vitro functional evidence for potential causal variants and target genes at MVP risk loci involving established and new biological mechanisms. Graphic Abstract: A graphic abstract is available for this article.

Loop extrusion as a mechanism for formation of DNA damage repair foci

The repair of DNA double-strand breaks (DSBs) is essential for safeguarding genome integrity. When a DSB forms, the PI3K-related ATM kinase rapidly triggers the establishment of megabase-sized, chromatin domains decorated with phosphorylated histone H2AX (γH2AX), which act as seeds for the formation of DNA-damage response foci1. It is unclear how these foci are rapidly assembled to establish a 'repair-prone' environment within the nucleus. Topologically associating domains are a key feature of 3D genome organization that compartmentalize transcription and replication, but little is known about their contribution to DNA repair processes2,3. Here we show that topologically associating domains are functional units of the DNA damage response, and are instrumental for the correct establishment of γH2AX-53BP1 chromatin domains in a manner that involves one-sided cohesin-mediated loop extrusion on both sides of the DSB. We propose a model in which H2AX-containing nucleosomes are rapidly phosphorylated as they actively pass by DSB-anchored cohesin. Our work highlights the importance of chromosome conformation in the maintenance of genome integrity and demonstrates the establishment of a chromatin modification by loop extrusion.

Suppression of the NTS-CPS1 regulatory axis by AFF1 in lung adenocarcinoma cells

Upregulation of the neuropeptide neurotensin (NTS) in a subgroup of lung cancers has been linked to poor prognosis. However, the regulatory pathway centered on NTS in lung cancer remains unclear. Here we identified the NTS-specific enhancer in lung adenocarcinoma cells. The AF4/FMR2 (AFF) family protein AFF1 occupies the NTS enhancer and inhibits NTS transcription. Clustering analysis of lung adenocarcinoma gene expression data demonstrated that NTS expression is highly positively correlated with the expression of the oncogenic factor CPS1. Detailed analyses demonstrated that the IL6 pathway antagonizes NTS in regulating CPS1. Thus, our analyses revealed a novel NTS-centered regulatory axis, consisting of AFF1 as a master transcription suppressor and IL6 as an antagonist in lung adenocarcinoma cells.

SLE non-coding genetic risk variant determines the epigenetic dysfunction of an immune cell specific enhancer that controls disease-critical microRNA expression

Since most variants that impact polygenic disease phenotypes localize to non-coding genomic regions, understanding the consequences of regulatory element variants will advance understanding of human disease mechanisms. Here, we report that the systemic lupus erythematosus (SLE) risk variant rs2431697 as likely causal for SLE through disruption of a regulatory element, modulating miR-146a expression. Using epigenomic analysis, genome-editing and 3D chromatin structure analysis, we show that rs2431697 tags a cell-type dependent distal enhancer specific for miR-146a that physically interacts with the miR-146a promoter. NF-kB binds the disease protective allele in a sequence-specific manner, increasing expression of this immunoregulatory microRNA. Finally, CRISPR activation-based modulation of this enhancer in the PBMCs of SLE patients attenuates type I interferon pathway activation by increasing miR-146a expression. Our work provides a strategy to define non-coding RNA functional regulatory elements using disease-associated variants and provides mechanistic links between autoimmune disease risk genetic variation and disease etiology.

CTCF modulates allele-specific sub-TAD organization and imprinted gene activity at the mouse Dlk1-Dio3 and Igf2-H19 domains

BACKGROUND

Genomic imprinting is essential for mammalian development and provides a unique paradigm to explore intra-cellular differences in chromatin configuration. So far, the detailed allele-specific chromatin organization of imprinted gene domains has mostly been lacking. Here, we explored the chromatin structure of the two conserved imprinted domains controlled by paternal DNA methylation imprints-the Igf2-H19 and Dlk1-Dio3 domains-and assessed the involvement of the insulator protein CTCF in mouse cells.

RESULTS

Both imprinted domains are located within overarching topologically associating domains (TADs) that are similar on both parental chromosomes. At each domain, a single differentially methylated region is bound by CTCF on the maternal chromosome only, in addition to multiple instances of bi-allelic CTCF binding. Combinations of allelic 4C-seq and DNA-FISH revealed that bi-allelic CTCF binding alone, on the paternal chromosome, correlates with a first level of sub-TAD structure. On the maternal chromosome, additional CTCF binding at the differentially methylated region adds a further layer of sub-TAD organization, which essentially hijacks the existing paternal-specific sub-TAD organization. Perturbation of maternal-specific CTCF binding site at the Dlk1-Dio3 locus, using genome editing, results in perturbed sub-TAD organization and bi-allelic Dlk1 activation during differentiation.

CONCLUSIONS

Maternal allele-specific CTCF binding at the imprinted Igf2-H19 and the Dlk1-Dio3 domains adds an additional layer of sub-TAD organization, on top of an existing three-dimensional configuration and prior to imprinted activation of protein-coding genes. We speculate that this allele-specific sub-TAD organization provides an instructive or permissive context for imprinted gene activation during development.

Generation of transgenic zebrafish with 2 populations of RFP- and GFP-labeled thrombocytes: analysis of their lipids

Zebrafish thrombocytes are similar to mammalian platelets. Mammals have young platelets (also called reticulated platelets) and mature platelets. Likewise, zebrafish have 2 populations of thrombocytes; one is DiI-C18 (DiI)+ (DP), and the other is DiI- (DN). However, the mechanism of selective thrombocyte labeling by DiI is unknown. Furthermore, there is no transgenic zebrafish line where DP and DN thrombocytes are differentially labeled with fluorescent proteins. In this study, we found that Glo fish, in which the myosin light chain 2 promoter drives the rfp gene, have a population of thrombocytes that are red fluorescent protein (RFP) labeled. We also generated transgenic GloFli fish in which DP and DN thrombocytes are labeled with RFP and green fluorescent protein (GFP), respectively. Single-cell lipid analysis showed a twofold increase in phosphatidylethanolamine (PE) and a twofold decrease in phosphatidylcholine (PC) in RFP+ thrombocytes compared with GFP+ thrombocytes, suggesting that lipid composition may be important for DiI differential labeling. Therefore, we tested liposomes prepared with different ratios of PC and PE and observed that liposomes prepared with higher amounts of PE favor DiI labeling, whereas the PC concentration had a modest effect. In liposomes prepared using only PE or PC, increased concentrations of PE resulted in increased DiI binding. These results suggest that because RFP+ thrombocytes have higher PE concentrations, DiI may bind to them efficiently, thus explaining the selective labeling of thrombocytes by DiI. This work also provides GloFli fish that should be useful in understanding the mechanism of thrombocyte maturation.

Dynamic enhancer partitioning instructs activation of a growth-related gene during exit from naïve pluripotency

During early mammalian development, the chromatin landscape undergoes profound transitions. The Zdbf2 gene—involved in growth control—provides a valuable model to study this window: upon exit from naïve pluripotency and prior to tissue differentiation, it undergoes a switch from a distal to a proximal promoter usage, accompanied by a switch from polycomb to DNA methylation occupancy. Using a mouse embryonic stem cell (ESC) system to mimic this period, we show here that four enhancers contribute to the Zdbf2 promoter switch, concomitantly with dynamic changes in chromatin architecture. In ESCs, the locus is partitioned to facilitate enhancer contacts with the distal Zdbf2 promoter. Relieving the partition enhances proximal Zdbf2 promoter activity, as observed during differentiation or with genetic mutants. Importantly, we show that 3D regulation occurs upstream of the polycomb and DNA methylation pathways. Our study reveals the importance of multi-layered regulatory frameworks to ensure proper spatio-temporal activation of developmentally important genes.

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture (4C-seq)

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and effect sizes of regulatory elements to a target gene. Some progress has been made with the bioinformatic prediction of the existence and function of proximal epigenetic modifications associated with activated gene expression using conserved transcription factor binding sites. Chromatin conformation capture studies have revolutionized our ability to discover physical chromatin contacts between sequences and even within an entire genome. Circular chromatin conformation capture coupled with next-generation sequencing (4C-seq), in particular, is designed to discover all possible physical chromatin interactions for a given sequence of interest (viewpoint), such as a target gene or a regulatory enhancer. Current 4C-seq strategies directly sequence from within the viewpoint but require numerous and diverse viewpoints to be simultaneously sequenced to avoid the technical challenges of uniform base calling (imaging) with next generation sequencing platforms. This volume of experiments may not be practical for many laboratories. Here, we report a modified approach to the 4C-seq protocol that incorporates both an additional restriction enzyme digest and qPCR-based amplification steps that are designed to facilitate a greater capture of diverse sequence reads and mitigate the potential for PCR bias, respectively. Our modified 4C method is amenable to the standard molecular biology lab for assessing chromatin architecture.

Human Accelerated Regions and Other Human-Specific Sequence Variations in the Context of Evolution and Their Relevance for Brain Development

The review discusses, in a format of a timeline, the studies of different types of genetic variants, present in Homo sapiens, but absent in all other primate, mammalian, or vertebrate species, tested so far. The main characteristic of these variants is that they are found in regions of high evolutionary conservation. These sequence variations include single nucleotide substitutions (called human accelerated regions), deletions, and segmental duplications. The rationale for finding such variations in the human genome is that they could be responsible for traits, specific to our species, of which the human brain is the most remarkable. As became obvious, the vast majority of human-specific single nucleotide substitutions are found in noncoding, likely regulatory regions. A number of genes, associated with these human-specific alleles, often through novel enhancer activity, were in fact shown to be implicated in human-specific development of certain brain areas, including the prefrontal cortex. Human-specific deletions may remove regulatory sequences, such as enhancers. Segmental duplications, because of their large size, create new coding sequences, like new functional paralogs. Further functional study of these variants will shed light on evolution of our species, as well as on the etiology of neurodevelopmental disorders.

Krox20 hindbrain regulation incorporates multiple modes of cooperation between cis-acting elements

Developmental genes can harbour multiple transcriptional enhancers that act simultaneously or in succession to achieve robust and precise spatiotemporal expression. However, the mechanisms underlying cooperation between cis-acting elements are poorly documented, notably in vertebrates. The mouse gene Krox20 encodes a transcription factor required for the specification of two segments (rhombomeres) of the developing hindbrain. In rhombomere 3, Krox20 is subject to direct positive feedback governed by an autoregulatory enhancer, element A. In contrast, a second enhancer, element C, distant by 70 kb, is active from the initiation of transcription independent of the presence of the KROX20 protein. Here, using both enhancer knock-outs and investigations of chromatin organisation, we show that element C possesses a dual activity: besides its classical enhancer function, it is also permanently required in cis to potentiate the autoregulatory activity of element A, by increasing its chromatin accessibility. This work uncovers a novel, asymmetrical, long-range mode of cooperation between cis-acting elements that might be essential to avoid promiscuous activation of positive autoregulatory elements.

The formation of multicellular organisms from the egg to the adult stage is largely under genetic control. The activation of specific genes is governed by regulatory DNA sequences present nearby on the chromosome. Most of these sequences promote activation and are called enhancers. In this paper, we study two enhancers governing the expression of a gene involved in the formation of the posterior brain in vertebrates. One of these enhancers is involved in a positive feedback loop: it is itself activated by the protein product of the gene that it regulates. The other enhancer was thought to be only involved in the initial accumulation of the protein, necessary for the subsequent activation of the feedback loop. Here we show that the second enhancer directly cooperates with the autoregulatory enhancer to increase its accessibility and its activity. Our work uncovers a novel, long-range mode of cooperation between enhancers that restricts the domain of action of autoregulatory enhancers within embryos and might be essential to avoid their inappropriate activation.

Interplay between genetic and epigenetic mechanisms in rheumatoid arthritis

Genetic and environmental factors contribute to the risk for rheumatoid arthritis (RA), with epigenetics serving as a possible interface through which risk factors contribute to RA. High-throughput technologies for interrogating genome and epigenome, and the availability of genetic and epigenetic datasets across a diversity of cell types, enable the identification of candidate causal genetic variants for RA to study their function in core RA processes. To date, RA risk variants were studied in the immune cells but not joint resident cells, for example, synovial fibroblasts. Synovial fibroblasts from different joints are distinct, anatomically specialized cells, defined by joint-specific transcriptomes, epigenomes and phenotypes. Cell type-specific analysis of epigenetic changes, together with genetic fine mapping and interrogation of chromatin 3D interactions may identify new disease relevant pathways, potential therapeutic targets and biomarkers for RA progression or therapy response.