Multiomic Profiling Identifies cis-Regulatory Networks Underlying Human Pancreatic β Cell Identity and Function.

EndoC-βH1 is emerging as a critical human β cell model to study the genetic and environmental etiologies of β cell (dys)function and diabetes. Comprehensive knowledge of its molecular landscape is lacking, yet required, for effective use of this model. Here, we report chromosomal (spectral karyotyping), genetic (genotyping), epigenomic (ChIP-seq and ATAC-seq), chromatin interaction (Hi-C and Pol2 ChIA-PET), and transcriptomic (RNA-seq and miRNA-seq) maps of EndoC-βH1. Analyses of these maps define known (e.g., PDX1 and ISL1) and putative (e.g., PCSK1 and mir-375) β cell-specific transcriptional cis-regulatory networks and identify allelic effects on cis-regulatory element use. Importantly, comparison with maps generated in primary human islets and/or β cells indicates preservation of chromatin looping but also highlights chromosomal aberrations and fetal genomic signatures in EndoC-βH1. Together, these maps, and a web application we created for their exploration, provide important tools for the design of experiments to probe and manipulate the genetic programs governing β cell identity and (dys)function in diabetes

Balanced chromosomal rearrangements offer insights into coding and noncoding genomic features associated with developmental disorders

Balanced chromosomal rearrangements (BCRs), including inversions, translocations, and insertions, reorganize large sections of the genome and contribute substantial risk for developmental disorders (DDs). However, the rarity and lack of systematic screening for BCRs in the population has precluded unbiased analyses of the genomic features and mechanisms associated with risk for DDs versus normal developmental outcomes. Here, we sequenced and analyzed 1,420 BCR breakpoints across 710 individuals, including 406 DD cases and the first large-scale collection of 304 control BCR carriers. We found that BCRs were not more likely to disrupt genes in DD cases than controls, but were seven-fold more likely to disrupt genes associated with dominant DDs (21.3% of cases vs. 3.4% of controls; P = 1.60×10$^{−12}$). Moreover, BCRs that did not disrupt a known DD gene were significantly enriched for breakpoints that altered topologically associated domains (TADs) containing dominant DD genes in cases compared to controls (odds ratio [OR] = 1.43, P = 0.036). We discovered six TADs enriched for noncoding BCRs (false discovery rate < 0.1) that contained known DD genes (MEF2C, FOXG1, SOX9, BCL11A, BCL11B, and SATB2) and represent candidate pathogenic long-range positional effect (LRPE) loci. These six TADs were collectively disrupted in 7.4% of the DD cohort. Phased Hi-C analyses of five cases with noncoding BCR breakpoints localized to one of these putative LRPEs, the 5q14.3 TAD encompassing MEF2C, confirmed extensive disruption to local 3D chromatin structures and reduced frequency of contact between the MEF2C promoter and annotated enhancers. We further identified six genomic features enriched in TADs preferentially disrupted by noncoding BCRs in DD cases versus controls and used these features to build a model to predict TADs at risk for LRPEs across the genome. These results emphasize the potential impact of noncoding structural variants to cause LRPEs in unsolved DD cases, as well as the complex interaction of features associated with predicting three-dimensional chromatin structures intolerant to disruption

Mechanistic dissection of chromatin topology disruption in the 5q14.3 MEF2C locus as an indirect driver of neurodevelopmental disorders

Structural variants have the potential to create long-range positional effects, uncouple genes from regulatory elements, and facilitate aberrant 3D chromatin folding. In an independent study, we revealed a significant enrichment of intergenic balanced chromosomal abnormality (BCA) breakpoints from congenital anomaly cases at chromosome 5q14.3. All 11 5q14.3 BCA carriers had breakpoints disrupting the topologically associated domain (TAD) housing MEF2C, a known driver of neurodevelopmental disorders. In a second study, we showed, using 4C-seq, that MEF2C forms proximal and distal interaction loops within this TAD, contacting multiple neuronal enhancers. To understand the functional impact of these BCAs, we performed a mechanistic dissection of the 3D regulatory network at the 5q14.3 locus and its constituent functional elements using Cas9-based genome editing. We generated >200 cell lines, representing deletions of MEF2C alongside five MEF2C TAD boundary and interaction loop targets in iPS-derived neural stem cells (NSCs) and cortical induced neurons (iNs). We profiled changes in expression and 3D chromatin interactions within these models using RNA-seq and 4C-seq. MEF2C was variably differentially expressed upon deletion of proximal and distal MEF2C loop regions, with effects depending on cell type and variant class. Underlying chromatin interaction patterns revealed evidence of loop maintenance in these models, possibly via CTCF buffering, highlighting compensatory mechanisms against 3D chromatin disruption. In contrast, deletion of the proximal TAD boundary facilitated increased contacts with predicted enhancers in the adjacent TAD. Our results suggest novel regulatory mechanisms driving phenotypic outcomes for the 5q14.3 region, with significant implications for interpretation of pathogenic structural variation

Chromatin alternates between A and B compartments at kilobase scale for subgenic organization

Abstract Nuclear compartments are prominent features of 3D chromatin organization, but sequencing depth limitations have impeded investigation at ultra fine-scale. CTCF loops are generally studied at a finer scale, but the impact of looping on proximal interactions remains enigmatic. Here, we critically examine nuclear compartments and CTCF loop-proximal interactions using a combination of in situ Hi-C at unparalleled depth, algorithm development, and biophysical modeling. Producing a large Hi-C map with 33 billion contacts in conjunction with an algorithm for performing principal component analysis on sparse, super massive matrices (POSSUMM), we resolve compartments to 500 bp. Our results demonstrate that essentially all active promoters and distal enhancers localize in the A compartment, even when flanking sequences do not. Furthermore, we find that the TSS and TTS of paused genes are often segregated into separate compartments. We then identify diffuse interactions that radiate from CTCF loop anchors, which correlate with strong enhancer-promoter interactions and proximal transcription. We also find that these diffuse interactions depend on CTCF’s RNA binding domains. In this work, we demonstrate features of fine-scale chromatin organization consistent with a revised model in which compartments are more precise than commonly thought while CTCF loops are more protracted