Dataset Information

Epigenetic maintenance of topological domains in the rearranged gibbon genome [ChIP-seq]

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

ORGANISM(S): Nomascus leucogenys

PROVIDER: GSE115064 | GEO | 2018/05/30

REPOSITORIES: GEO

ACCESS DATA

Shared Molecules

Only show the datasets with similarity scores above: 0.5

Threshold

0.5

Similar Datasets

Project description:Three-dimensional genome structure plays an important role in gene regulation. Globally chromosomes are organized into active and inactive compartments, while at the gene level looping interactions connect promoters to regulatory elements. Topologically Associating Domains (TADs), typically several hundred kilobases in size form an intermediate level of organization. Major questions include how TADs are formed and what their relation is with looping interactions between genes and regulatory elements. Here we performed a focused 5C analysis of a 2.8 Mb region on chromosome 7 surrounding CFTR in a panel of cell types. We find that the same TAD boundaries are present in all cell types, indicating that TADs represent a universal chromosome architecture. Further, we find that these TAD boundaries are present irrespective of expression and looping of genes located between them. In contrast looping interactions between promoters and regulatory elements are cell-type specific and occur mostly within TADs. This is exemplified by the CFTR promoter that in different cell types interacts with distinct sets of distal cell type-specific regulatory elements that are all located within the same TAD. Finally, we find that long-range associations between loci located in different TADs are also detected but these display much lower interaction frequencies than looping interactions within TADs. Interestingly, interactions between TADs are also highly cell type-specific and often involve loci clustered around TAD boundaries. These data point to key roles of invariant TAD boundaries in constraining as well as mediating cell type-specific long-range interactions and gene regulation. We investigated a 2.8 Mb region on Chromosome 7 (hg18 chr7: 115797757-118405450) containing the ENCODE region ENm001 42. The 5C experiment was designed to interrogate looping interactions between HindIII fragments containing transcription start sites (TSSs) and any other HindIII restriction fragment (distal fragments) in the target region. Libraries were generated for five cell lines: Caco2, Calu3, Capan1, GM12878 and HepG2, with two biological replicates for each line. 5C probes were designed at HindIII restriction sites (AAGCTT) using 5C primer design tools previously developed and made publicly available online at our My5C website (http://my5C.umassmed.edu). Probes were designed based on the ENCODE manual region 1 (ENM001) design 25 with additional probes placed throughout the region when appropriate. We also added probes to extend the analysis to include a 700 Kb gene desert region located directly adjacent to ENM001. Probe settings were: U-BLAST, 3; S-BLAST, 100; 15-MER, 3,000; MIN_FSIZE, 250; MAX_FSIZE, 20,000; OPT_TM, 65; OPT_PSIZE, 40. We designed 74 reverse 5C probes, and 605 forward 5C probes.