Search

Structural variants (SVs) that alter DNA sequence emerge as a driving force involved in the reorganisation of DNA spatial folding, thus affecting gene transcription. In this work, we describe an improved version of our integrated web service for structural modeling of three-dimensional genome (3D-GNOME), which now incorporates all types of SVs to model changes to the reference 3D conformation of chromatin. In 3D-GNOME 2.0, the default reference 3D genome structure is generated using ChIA-PET data from the GM12878 cell line and SVs data are sourced from the population-scale catalogue of SVs identified by the 1000 Genomes Consortium. However, users may also submit their own structural data to set a customized reference genome structure, and/or a custom input list of SVs. 3D-GNOME 2.0 provides novel tools to inspect, visualize and compare 3D models for regions that differ in terms of their linear genomic sequence. Contact diagrams are displayed to compare the reference 3D structure with the one altered by SVs. In our opinion, 3D-GNOME 2.0 is a unique online tool for modeling and analyzing conformational changes to the human genome induced by SVs across populations. It can be freely accessed at https://3dgnome.cent.uw.edu.pl/. Polish National Science Centre [2014/15/B/ST6/05082]; Foundation for Polish Science co-financed by the European Union under the European Regional Development Fund (TEAM to D.P.); 'Nucleome Positioning System for Spatiotemporal Genome Organization and Regulation' within the 4DNucleome National Institute of Health program, and by the European Commission as European Cooperation in Science and Technology COST actions: CA18127 'International Nucleome Consortium' (INC) [1U54DK107967-01]; 'Impact of Nuclear Domains On Gene Expression and Plant Traits' [CA16212]; Horizon 2020 Marie Sklodowska-Curie ITN Enhpathy grant 'Molecular Basis of Human enhanceropathies'; D.P. and M.W. were supported by the RENOIR Project from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie [691152]; Ministry of Science and Higher Education (Poland) [W34/H2020/2016, 329025/PnH/2016]. Funding for open access charge: Foundation for Polish Science co-financed by the European Union under the European Regional Development Fund (TEAM to DP). ACKNOWLEDGEMENTS We thank Veronika Mancikova for language editing and proofreading of the manuscript and Agnieszka Bucka and Agnes Alcantara Paculdar for critical reading of the text. Plewczynski, D (corresponding author), Univ Warsaw, Ctr New Technol, PL-02097 Warsaw, Poland ; Warsaw Univ Technol, Fac Math & Informat Sci, PL-00662 Warsaw, Poland ; Jackson Lab Genom Med, Farmington, CT 06032 USA. d.plewczynski@cent.uw.edu.pl

Blocking PD-1/PD-L1 signaling transforms cancer therapy and is assumed to unleash exhausted tumor-reactive CD8

Mitochondria are essential organelles in eukaryotic cells that provide critical support for energetic and metabolic homeostasis. Mutations that accumulate in mitochondrial DNA (mtDNA) in somatic cells have been implicated in cancer, degenerative diseases, and the aging process. However, the mechanisms used by somatic cells to maintain proper functions despite their mtDNA mutation load are poorly understood. Here, we analyzed somatic mtDNA mutations in more than 30,000 human single peripheral and bone marrow mononuclear cells and observed a significant overrepresentation of homoplastic mtDNA mutations in B, T and NK lymphocytes despite their lower mutational burden than other hematopoietic cells. The characteristic mutational landscape of mtDNA in lymphocytes were validated with data from multiple platforms and individuals. Single-cell RNA-seq and computational modeling demonstrated a stringent mitochondrial bottleneck during lymphocyte development likely caused by lagging mtDNA replication relative to cell proliferation. These results illuminate a potential mechanism used by highly metabolically active immune cells for quality control of their mitochondrial genomes.

Thyroid carcinoma (TC) is the most common endocrine malignancy, and papillary TC (PTC) is the most frequent subtype of TC, accounting for 85–90% of all the cases. Aberrant histone acetylation contributes to carcinogenesis by inducing the dysregulation of certain cancer-related genes. However, the histone acetylation landscape in PTC remains elusive. Here, we interrogated the epigenomes of PTC and benign thyroid nodule (BTN) tissues by applying H3K27ac chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) along with RNA-sequencing. By comparing the epigenomic features between PTC and BTN, we detected changes in H3K27ac levels at active regulatory regions, identified PTC-specific super-enhancer-associated genes involving immune-response and cancer-related pathways, and uncovered several genes that associated with disease-free survival of PTC. In summary, our data provided a genome-wide landscape of histone modification in PTC and demonstrated the role of enhancers in transcriptional regulations associated with prognosis of PTC.

HIV‐1 latency is a major obstacle to achieving a functional cure for AIDS. Reactivation of HIV‐1‐infected cells followed by their elimination via immune surveillance is one proposed strategy for eradicating the viral reservoir. However, current latency‐reversing agents (LRAs) show high toxicity and low efficiency, and new targets are needed to develop more promising LRAs. Here, we found that the histone chaperone CAF‐1 (chromatin assembly factor 1) is enriched on the HIV‐1 long terminal repeat (LTR) and forms nuclear bodies with liquid–liquid phase separation (LLPS) properties. CAF‐1 recruits epigenetic modifiers and histone chaperones to the nuclear bodies to establish and maintain HIV‐1 latency in different latency models and primary CD4(+) T cells. Three disordered regions of the CHAF1A subunit are important for phase‐separated CAF‐1 nuclear body formation and play a key role in maintaining HIV‐1 latency. Disruption of phase‐separated CAF‐1 bodies could be a potential strategy to reactivate latent HIV‐1.

New genome pipeline enables seamless execution of processing, analyzing, and visualizing of ChIA-PET data., ChIA-PET (chromatin interaction analysis with paired-end tags) enables genome-wide discovery of chromatin interactions involving specific protein factors, with base pair resolution. Interpretation of ChIA-PET data requires a robust analytic pipeline. Here, we introduce ChIA-PIPE, a fully automated pipeline for ChIA-PET data processing, quality assessment, visualization, and analysis. ChIA-PIPE performs linker filtering, read mapping, peak calling, and loop calling and automates quality control assessment for each dataset. To enable visualization, ChIA-PIPE generates input files for two-dimensional contact map viewing with Juicebox and HiGlass and provides a new dockerized visualization tool for high-resolution, browser-based exploration of peaks and loops. To enable structural interpretation, ChIA-PIPE calls chromatin contact domains, resolves allele-specific peaks and loops, and annotates enhancer-promoter loops. ChIA-PIPE also supports the analysis of other related chromatin-mapping data types.

Additional file 2: Supplementary figures. Figure S1. Comparison of CTCF and RNAPII binding and looping in PR9 and PR9+Zn cells. Figure S2. PML-RARα-mediated chromatin interaction domain and impacts on the myeloid genome. Figure S3. Transcription repression of myeloid specific genes by PML-RARα. Figure S4. Further details of motif analysis around PML-RARα binding sites. Figure S5. Super-enhancers of myeloid specific genes affected by PML-RARα. Figure S6. Additional examples of PML-RARα in NB4 cell.

Additional file 3. Supplementary methods. 1. ChIA-PET data processing and clustering of interaction PETs. 2. Comparison of ChIA-PET data for calling of binding peaks, loops, and chromatin domains. 3. Identification and classification of RARα and PML-RARα binding sites. 4. Identification and characterization of PML-RARα-mediated chromatin interactions. 5. Visualization of 2D contact map and BASIC browser. 6. 3D chromatin structure modeling, measurement of 3D structure condensation and visualization. 7. Identification of binding motifs at PML-RARα peak regions. 8. Comparison of RNAPII-mediated chromatin interactions between control and treatment condition. 9. RNA-Seq and differential gene expression analysis. 10. ChIP-Seq analysis. 11. Identification and characterization of super-enhancers. 12. Quantification and statistical analysis methods.

Cytokines critically control immune responses, but how regulatory programs are altered to allow T cells to differentially respond to distinct cytokine stimuli remains poorly understood. Here, we have globally analyzed enhancer elements bound by IL-2-activated STAT5 and IL-21-activated STAT3 in T cells and identified Il2ra as the top-ranked gene regulated by an IL-2-activated STAT5-bound superenhancer and one of the top genes regulated by STAT3-bound superenhancers. Moreover, we found that STAT5 binding was rapidly superenriched at genes highly induced by IL-2 and that IL-2-activated STAT5 binding induces new and augmented chromatin interactions within superenhancer-containing genes. Based on chromatin interaction analysis by paired-end tag (ChIA-PET) sequencing data, we used CRISPR-Cas9 gene editing to target three of the STAT5 binding sites within the Il2ra superenhancer in mice. Each mutation decreased STAT5 binding and altered IL-2-induced Il2ra gene expression, revealing that individual elements within the superenhancer were not functionally redundant and that all were required for normal gene expression. Thus, we demonstrate cooperative utilization of superenhancer elements to optimize gene expression and show that STAT5 mediates IL-2-induced chromatin looping at superenhancers to preferentially regulate highly inducible genes, thereby providing new insights into the mechanisms underlying cytokine-dependent superenhancer function.

Chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) is a robust method for capturing genome-wide chromatin interactions. Unlike other 3C-based methods, it includes a chromatin immunoprecipitation (ChIP) step that enriches for interactions mediated by specific target proteins. This unique feature allows ChIA-PET to provide the functional specificity and higher resolution needed to detect chromatin interactions, which chromosome conformation capture (3C)/Hi-C approaches have not achieved. The original ChIA-PET protocol generates short paired-end tags (2 × 20 base pairs (bp)) to detect two genomic loci that are far apart on linear chromosomes but are in spatial proximity in the folded genome. We have improved the original approach by developing long-read ChIA-PET, in which the length of the paired-end tags is increased (up to 2 × 250 bp). The longer PET reads not only improve the tag-mapping efficiency but also increase the probability of covering phased single-nucleotide polymorphisms (SNPs), which allows haplotype-specific chromatin interactions to be identified. Here, we provide the detailed protocol for long-read ChIA-PET that includes cell fixation and lysis, chromatin fragmentation by sonication, ChIP, proximity ligation with a bridge linker, Tn5 tagmentation, PCR amplification and high-throughput sequencing. For a well-trained molecular biologist, it typically takes 6 d from cell harvesting to the completion of library construction, up to a further 36 h for DNA sequencing and

The three-dimensional architecture of chromatin in the nucleus is important for genome regulation and function. Advanced high-throughput sequencing-based methods have been developed for capturing chromatin interactions (Hi-C, genome-wide chromosome conformation capture) or enriching for those involving a specific protein (ChIA-PET, chromatin interaction analysis with paired-end tag sequencing). There is widespread interest in utilizing and interpreting ChIA-PET and Hi-C. We review methods for comparative ChIA-PET and Hi-C data analysis and visualization. The topics reviewed include: downloading ChIA-PET and Hi-C data from the ENCODE and 4DN portals; processing ChIA-PET data using ChIA-PIPE; processing Hi-C data using Juicer or distiller and cooler; viewing 2D contact maps using Juicebox or Higlass; viewing peaks, loops, and domains using BASIC Browser; annotating convergent and tandem CTCF loops.

Supplemental figures. (PDF 17323 kb)

ChIA-PET enables the genome-wide discovery of chromatin interactions involving specific protein factors, with base-pair resolution. Interpreting ChIA-PET data depends on having a robust analytic pipeline. Here, we introduce ChIA-PIPE, a fully automated pipeline for ChIA-PET data processing, quality assessment, analysis, and visualization. ChIA-PIPE performs linker filtering, read mapping, peak calling, loop calling, chromatin-contact-domain calling, and can resolve allele-specific peaks and loops. ChIA-PIPE also automates quality-control assessment for each dataset. Furthermore, ChIA-PIPE generates input files for visualizing 2D contact maps with Juicebox and HiGlass, and provides a new dockerized visualization tool for high-resolution, browser-based exploration of peaks and loops. With minimal adjusting, ChIA-PIPE can also be suited for the analysis of other related chromatin-mapping data.

This genome-wide study is focused on the impact of structural variants identified in individuals from 26 human populations onto three-dimensional structures of their genomes. We assess the tendency of structural variants to accumulate in spatially interacting genomic segments and design a high-resolution computational algorithm to model the 3D conformational changes resulted by structural variations. We show that differential gene transcription is closely linked to variation in chromatin interaction networks mediated by RNA polymerase II. We also demonstrate that CTCF-mediated interactions are well conserved across population, but enriched with disease-associated SNPs. Altogether, this study assesses the critical impact of structural variants on the higher order organization of chromatin folding and provides unique insight into the mechanisms regulating gene transcription at the population scale, among which the local arrangement of chromatin loops seems to be the leading one. It is the first insight into the variability of the human 3D genome at the population scale.

SUMMARY The tandem duplicator phenotype (TDP) is a genome-wide instability configuration primarily observed in breast, ovarian, and endometrial carcinomas. Here, we stratify TDP tumors by classifying their tandem duplications (TDs) into three span intervals, with modal values of 11 kb, 231 kb, and 1.7 Mb, respectively. TDPs with ~11 kb TDs feature loss of TP53 and BRCA1. TDPs with ~231 kb and ~1.7 Mb TDs associate with CCNE1 pathway activation and CDK12 disruptions, respectively. We demonstrate that p53 and BRCA1 conjoint abrogation drives TDP induction by generating short-span TDP mammary tumors in genetically modified mice lacking them. Lastly, we show how TDs in TDP tumors disrupt heterogeneous combinations of tumor suppressors and chromatin topologically associating domains while duplicating oncogenes and super-enhancers., Graphical Abstract, In Brief Menghi et al. stratify tandem duplicator phenotype tumors by classifying their tandem duplications (TDs) into three span sizes associated with different pathway alterations and show how TDs disrupt tumor suppressors and chromatin topologically associating domains while duplicating oncogenes and super-enhancers.

Superenhancers regulate the expression of genes that specify cell type-specific development, but little is known regarding their function and regulation in vivo. Here, we study the cytokines IL-2 and IL-21, which critically control the immune response. These cytokines induce the binding of transcription factors STAT5 and STAT3, respectively, at superenhancers in a cytokine- and gene-specific manner. STAT5-bound superenhancers regulate genes highly induced by IL-2, with STAT5 mediating chromatin looping within such genes, including Il2ra, which mediates responsiveness to IL-2. By deleting three STAT5 binding sites that mediate IL-2–induced chromatin looping at the Il2ra locus, we demonstrate that superenhancer elements cooperatively control gene expression. Overall, we elucidate cytokine-dependent superenhancer function in general and provide detailed analysis of the Il2ra superenhancer.

Enhancers regulate spatiotemporal gene expression and impart cell-specific transcriptional outputs that drive cell identity. Super-enhancers (SEs), also known as stretch-enhancers, are a subset of enhancers especially important for genes associated with cell identity and genetic risk of disease. CD4(+) T cells are critical for host defence and autoimmunity. Here we analysed maps of mouse T-cell SEs as a non-biased means of identifying key regulatory nodes involved in cell specification. We found that cytokines and cytokine receptors were the dominant class of genes exhibiting SE architecture in T cells. Nonetheless, the locus encoding Bach2, a key negative regulator of effector differentiation, emerged as the most prominent T-cell SE, revealing a network in which SE-associated genes critical for T-cell biology are repressed by BACH2. Disease-associated single-nucleotide polymorphisms for immune-mediated disorders, including rheumatoid arthritis, were highly enriched for T-cell SEs versus typical enhancers or SEs in other cell lineages. Intriguingly, treatment of T cells with the Janus kinase (JAK) inhibitor tofacitinib disproportionately altered the expression of rheumatoid arthritis risk genes with SE structures. Together, these results indicate that genes with SE architecture in T cells encompass a variety of cytokines and cytokine receptors but are controlled by a 'guardian' transcription factor, itself endowed with an SE. Thus, enumeration of SEs allows the unbiased determination of key regulatory nodes in T cells, which are preferentially modulated by pharmacological intervention.

Oranges are an important nutritional source for human health and have immense economic value. Here we present a comprehensive analysis of the draft genome of sweet orange (Citrus sinensis). The assembled sequence covers 87.3% of the estimated orange genome, which is relatively compact, as 20% is composed of repetitive elements. We predicted 29,445 protein-coding genes, half of which are in the heterozygous state. With additional sequencing of two more citrus species and comparative analyses of seven citrus genomes, we present evidence to suggest that sweet orange originated from a backcross hybrid between pummelo and mandarin. Focused analysis on genes involved in vitamin C metabolism showed that GalUR, encoding the rate-limiting enzyme of the galacturonate pathway, is significantly upregulated in orange fruit, and the recent expansion of this gene family may provide a genomic basis. This draft genome represents a valuable resource for understanding and improving many important citrus traits in the future.

Recent advances in high-throughput chromosome conformation capture (3C) technology, such as HiC and ChIA-PET, have demonstrated the importance of 3D genome organization in development, cell differentiation and transcriptional regulation. There is now a widespread need for computational tools to generate and analyze 3D structural models from 3C data. Here we introduce our 3D GeNOme Modeling Engine (3D-GNOME), a web service which generates 3D structures from 3C data and provides tools to visually inspect and annotate the resulting structures, in addition to a variety of statistical plots and heatmaps which characterize the selected genomic region. Users submit a bedpe (paired-end BED format) file containing the locations and strengths of long range contact points, and 3D-GNOME simulates the structure and provides a convenient user interface for further analysis. Alternatively, a user may generate structures using published ChIA-PET data for the GM12878 cell line by simply specifying a genomic region of interest. 3D-GNOME is freely available at http://3dgnome.cent.uw.edu.pl/. National Leading Research Centre in Bialystok [to P.S., D.P.]; European Union under the European Social Fund [to P.S., D.P.]; Polish National Science Centre [2014/15/B/ST6/05082, 2013/09/B/NZ2/00121 to D.P., G.M., P.W.]; European Cooperation in Science and Technology action [COST BM1405, BM1408 to D.P., G.M., P.W.]; Director Innovation Fund from Jackson Laboratory and NCI [R01 CA186714 to Y.R.]. Y.R. is also supported by the Roux family as a Florine Roux Endowed Chair and Professor in Genomics and Computational Biology. Funding for open access charge: Polish National Science Centre [2014/15/B/ST6/05082].

ChIA-PET is a high-throughput mapping technology that reveals long-range chromatin interactions and provides insights into the basic principles of spatial genome organization and gene regulation mediated by specific protein factors. Recently, we showed that a single ChIA-PET experiment provides information at all genomic scales of interest, from the high-resolution locations of binding sites and enriched chromatin interactions mediated by specific protein factors, to the low resolution of nonenriched interactions that reflect topological neighborhoods of higher-order chromosome folding. This multilevel nature of ChIA-PET data offers an opportunity to use multiscale 3D models to study structural-functional relationships at multiple length scales, but doing so requires a structural modeling platform. Here, we report the development of 3D-GNOME (3-Dimensional Genome Modeling Engine), a complete computational pipeline for 3D simulation using ChIA-PET data. 3D-GNOME consists of three integrated components: a graph-distance-based heat map normalization tool, a 3D modeling platform, and an interactive 3D visualization tool. Using ChIA-PET and Hi-C data derived from human B-lymphocytes, we demonstrate the effectiveness of 3D-GNOME in building 3D genome models at multiple levels, including the entire genome, individual chromosomes, and specific segments at megabase (Mb) and kilobase (kb) resolutions of single average and ensemble structures. Further incorporation of CTCF-motif orientation and high-resolution looping patterns in 3D simulation provided additional reliability of potential biologically plausible topological structures. D.P., P.S., M.S., and M.J.P. are supported by grants from the Polish National Science Centre (2014/15/B/ST6/05082 and 2013/09/B/ NZ2/00121), and the European Cooperation in Science and Technology action (COST BM1405 and BM1408). P.S. and D.P. are supported by funds from National Leading Research Centre in Bialystok and the European Union under the European Social Fund. Y.R. is supported by the Director Innovation Fund from the Jackson Laboratory and National Cancer Institute (NCI) grant R01 CA186714. Y.R. is also supported by the Roux family as a Florine Roux Endowed Chair and Professor in Genomics and Computational Biology. All authors were supported by National Institutes of Health grant 1U54DK107967-01, “Nucleome Positioning System for Spatiotemporal Genome Organization and Regulation,” within the 4DNucleome NIH program.

The male sterility of thermosensitive genic male sterile (TGMS) lines of wheat (Triticum aestivum) is strictly controlled by temperature. The early phase of anther development is especially susceptible to cold stress. MicroRNAs (miRNAs) play an important role in plant development and in responses to environmental stress. In this study, deep sequencing of small RNA (smRNA) libraries obtained from spike tissues of the TGMS line under cold and control conditions identified a total of 78 unique miRNA sequences from 30 families and trans-acting small interfering RNAs (tasiRNAs) derived from two TAS3 genes. To identify smRNA targets in the wheat TGMS line, we applied the degradome sequencing method, which globally and directly identifies the remnants of smRNA-directed target cleavage. We identified 26 targets of 16 miRNA families and three targets of tasiRNAs. Comparing smRNA sequencing data sets and TaqMan quantitative polymerase chain reaction results, we identified six miRNAs and one tasiRNA (tasiRNA-ARF [for Auxin-Responsive Factor]) as cold stress-responsive smRNAs in spike tissues of the TGMS line. We also determined the expression profiles of target genes that encode transcription factors in response to cold stress. Interestingly, the expression of cold stress-responsive smRNAs integrated in the auxin-signaling pathway and their target genes was largely noncorrelated. We investigated the tissue-specific expression of smRNAs using a tissue microarray approach. Our data indicated that miR167 and tasiRNA-ARF play roles in regulating the auxin-signaling pathway and possibly in the developmental response to cold stress. These data provide evidence that smRNA regulatory pathways are linked with male sterility in the TGMS line during cold stress.

Gossypium hirsutum contributes the most production of cotton fibre, but G. barbadense is valued for its better comprehensive resistance and superior fibre properties. However, the allotetraploid genome of G. barbadense has not been comprehensively analysed. Here we present a high-quality assembly of the 2.57 gigabase genome of G. barbadense, including 80,876 protein-coding genes. The double-sized genome of the A (or At) (1.50 Gb) against D (or Dt) (853 Mb) primarily resulted from the expansion of Gypsy elements, including Peabody and Retrosat2 subclades in the Del clade and the Athila subclade in the Athila/Tat clade. Substantial gene expansion and contraction were observed and rich homoeologous gene pairs with biased expression patterns were identified, suggesting abundant gene sub-functionalization occurred by allopolyploidization. More specifically, the CesA gene family has adapted differentially temporal expression patterns, suggesting an integrated regulatory mechanism of CesA genes from At and Dt subgenomes for the primary and secondary cellulose biosynthesis of cotton fibre in a “relay race”-like fashion. We anticipate that the G. barbadense genome sequence will advance our understanding the mechanism of genome polyploidization and underpin genome-wide comparison research in this genus.

Thermo-photoperiod-sensitive genic male sterile (TPGMS) wheat is important in utilization of heterosis. To facilitate the use of such wheat line in agriculture, more knowledge about molecular mechanisms of TPGMS genes is required. In this study, we set up a normalized complementary DNA (cDNA) library based on the strategy of saturation hybridization with genomic DNA using TPGMS wheat line. This normalized cDNA library consists of cDNA from six directionally cloned cDNA libraries constructed with spike and anther tissues from spike developmental stages. From the normalized cDNA library, 3,264 single-pass expressed sequence tag (EST) were obtained. Exclusion of sequences shorter than 100 bp resulted in 3,223 vector-trimmed ESTs with a mean length of 926 bp. Clustering and assembly analysis resulted in 2,175 unique ESTs from 423 contigs and 1,752 singletons. Taking advantage of various tools and database, gene function classification showed that 60% of the ESTs were predicted to have putative gene function. Of the 2,175 unique ESTs, 264 (12%) displayed significant homology (BlastX E values

International audience; Spatial genome organization and its effect on transcription remains a fundamental question. We applied an advanced chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) strategy to comprehensively map higher-order chromosome folding and specific chromatin interactions mediated by CCCTC-binding factor (CTCF) and RNA polymerase II (RNAPII) with haplotype specificity and nucleotide resolution in different human cell lineages. We find that CTCF/cohesin-mediated interaction anchors serve as structural foci for spatial organization of constitutive genes concordant with CTCF-motif orientation, whereas RNAPII interacts within these structures by selectively drawing cell-type-specific genes toward CTCF foci for coordinated transcription. Furthermore, we show that haplotype variants and allelic interactions have differential effects on chromosome configuration, influencing gene expression, and may provide mechanistic insights into functions associated with disease susceptibility. 3D genome simulation suggests a model of chromatin folding around chromosomal axes, where CTCF is involved in defining the interface between condensed and open compartments for structural regulation. Our 3D genome strategy thus provides unique insights in the topological mechanism of human variations and diseases.

Dynamic changes in chromatin looping structure are identified to reflect interactions of enhancers to target promoters. Such changes modulate hematopoietic lineage differentiation and are known to be mediated by some critical transcription factor, such as GATA1 and MYC. How chromatin topology changes in cancer such as leukemia is poorly understood on a global genomic view. Moreover, how oncogenic transcription factors, such as PML-RARA, AML1-ETO and others, manipulate gene regulation by globally altering the topology of chromatin, which contribute to leukemogenesis in acute myeloid leukemia (AML) remains to be described. Here we identified PML-RARA involved chromosomal interactions spanning hundreds of kilobases between promoters and distal regulatory elements. The RNAPII mediated chromatin connectivity between transcriptionally active genes and its distal regulatory elements were lost in the reorganized chromatin structure caused by PML-RARA, while the chromatin topological associated domain mediated by CTCF remained stable. PML-RARA mediated chromatin loops eliminated the occupancy of myeloid specific transcription factors (such as PU.1, IRF1 and CEBPB et al.) as well as coordinately general transcriptional factors (such as P300 and RNAPII) at the distal regulatory elements and respective promoters. Finally we show that interacting loci designated as super enhancers become depleted by PML-RARA looping to suppress myeloid differentiation regulomes. These distinct results based on 3D genome architecture uncover novel modes of how PML-RARA disrupts the regulome of myeloid differentiation to contribute to leukemogenesis. This study provides new insights and accessible tools to delineate the aberrant regulome originated from the chaotic gene regulatory interactions leading to cancer genesis. Disclosures No relevant conflicts of interest to declare.

SummaryA key finding of the ENCODE project is that the enhancer landscape of mammalian cells undergoes marked alterations during ontogeny. However, the nature and extent of these changes are unclear. As part of the NIH Mouse Regulome Project, we here combined DNaseI hypersensitivity, ChIP-seq, and ChIA-PET technologies to map the promoter-enhancer interactomes of pluripotent ES cells and differentiated B lymphocytes. We confirm that enhancer usage varies widely across tissues. Unexpectedly, we find that this feature extends to broadly transcribed genes, including Myc and Pim1 cell-cycle regulators, which associate with an entirely different set of enhancers in ES and B cells. By means of high-resolution CpG methylomes, genome editing, and digital footprinting, we show that these enhancers recruit lineage-determining factors. Furthermore, we demonstrate that the turning on and off of enhancers during development correlates with promoter activity. We propose that organisms rely on a dynamic enhancer landscape to control basic cellular functions in a tissue-specific manner.

Summary The antibody gene mutator activation-induced cytidine deaminase (AID) promiscuously damages oncogenes, leading to chromosomal translocations and tumorigenesis. Why nonimmunoglobulin loci are susceptible to AID activity is unknown. Here, we study AID-mediated lesions in the context of nuclear architecture and the B cell regulome. We show that AID targets are not randomly distributed across the genome but are predominantly grouped within super-enhancers and regulatory clusters. Unexpectedly, in these domains, AID deaminates active promoters and eRNA + enhancers interconnected in some instances over megabases of linear chromatin. Using genome editing, we demonstrate that 3D-linked targets cooperate to recruit AID-mediated breaks. Furthermore, a comparison of hypermutation in mouse B cells, AID-induced kataegis in human lymphomas, and translocations in MEFs reveals that AID damages different genes in different cell types. Yet, in all cases, the targets are predominantly associated with topological complex, highly transcribed super-enhancers, demonstrating that these compartments are key mediators of AID recruitment.

Additional file 5. Review history.

Additional file 1: Supplementary tables. Table S1. Dataset used in this study. Table S2. Reagents used in this study. Table S3. Software and algorithms used in this study.

Additional file 5. Review history.

Additional file 1: Supplementary tables. Table S1. Dataset used in this study. Table S2. Reagents used in this study. Table S3. Software and algorithms used in this study.