Chromatin Interactions Research Articles

Mutations of PDS5 genes enhance TAD-like domain formation Arabidopsis thaliana

In eukaryotes, topologically associating domains (TADs) organize the genome into functional compartments. While TAD-like structures are common in mammals and many plants, they are challenging to detect in Arabidopsis thaliana. Here, we demonstrate that Arabidopsis PDS5 proteins play a negative role in TAD-like domain formation. Through Hi-C analysis, we show that mutations in PDS5 genes lead to the widespread emergence of enhanced TAD-like domains throughout the Arabidopsis genome, excluding pericentromeric regions. These domains exhibit increased chromatin insulation and enhanced chromatin interactions, without significant changes in gene expression or histone modifications. Our results suggest that PDS5 proteins are key regulators of genome architecture, influencing 3D chromatin organization independently of transcriptional activity. This study provides insights into the unique chromatin structure of Arabidopsis and the broader mechanisms governing plant genome folding.

Loopy virus or controlled contortionist? 3D regulation of HCMV gene expression by CTCF-driven chromatin interactions.

Three-dimensional chromatin control of eukaryotic transcription is pivotal for regulating gene expression. This additional layer of epigenetic regulation is also utilized by DNA viruses, including herpesviruses. Dynamic, spatial genomic organization often involves looping of chromatin anchored by host-encoded CCCTC-binding factor (CTCF) and other factors, which control crosstalk between promoters and enhancers. Herein, we review the contribution of CTCF-mediated looping in regulating transcription during herpesvirus infection, with a specific focus on the betaherpesvirus, human cytomegalovirus (HCMV).

Pharmacologic Induction of ERα SUMOylation Disrupts Its Chromatin Binding.

Estrogen receptor α (ERα)-positive breast cancer patients are typically treated with ERα inhibitors, including selective estrogen receptor modulators (SERMs) and selective estrogen receptor degraders (SERDs). However, the distinct pharmacological properties of various ERα inhibitors remain incompletely understood. In this study, we employed formaldehyde cross-linking followed by ERα immunoprecipitation and mass spectrometry to reveal that fulvestrant, the first FDA-approved SERD, induces the interaction between ERα and SUMO E3 ligases PIAS1 and PIAS2. Biochemical and genomic assays confirmed that fulvestrant induces SUMOylation of ERα, which inhibits ERα's binding to chromatin DNA. In addition, raloxifene (a SERM) and elacestrant (the first FDA-approved oral SERD) were identified as compounds that similarly induce ERα SUMOylation and inhibit its chromatin interaction. Our findings reveal a mechanism by which select ERα inhibitors disrupt ERα function through SUMOylation, offering insights for the development of next-generation ERα-targeted therapies.

Chromatin interaction maps identify oncogenic targets of enhancer duplications in cancer.

As a major type of structural variants, tandem duplication plays a critical role in tumorigenesis by increasing oncogene dosage. Recent work has revealed that noncoding enhancers are also affected by duplications leading to the activation of oncogenes that are inside or outside of the duplicated regions. However, the prevalence of enhancer duplication and the identity of their target genes remains largely unknown in the cancer genome. Here, by analyzing whole-genome sequencing data in a non-gene-centric manner, we identify 881 duplication hotspots in 13 major cancer types, most of which do not contain protein-coding genes. We show that the hotspots are enriched with distal enhancer elements and are highly lineage-specific. We develop a HiChIP-based methodology that navigates enhancer-promoter contact maps to prioritize the target genes for the duplication hotspots harboring enhancer elements. The methodology identifies many novel enhancer duplication events activating oncogenes such as ESR1, FOXA1, GATA3, GATA6, TP63, and VEGFA, as well as potentially novel oncogenes such as GRHL2, IRF2BP2, and CREB3L1 In particular, we identify a duplication hotspot on Chromosome 10p15 harboring a cluster of enhancers, which skips over two genes, through a long-range chromatin interaction, to activate an oncogenic isoform of the NET1 gene to promote migration of gastric cancer cells. Focusing on tandem duplications, our study substantially extends the catalog of noncoding driver alterations in multiple cancer types, revealing attractive targets for functional characterization and therapeutic intervention.

HiCDiffusion - diffusion-enhanced, transformer-based prediction of chromatin interactions from DNA sequences

Prediction of chromatin interactions from DNA sequence has been a significant research challenge in the last couple of years. Several solutions have been proposed, most of which are based on encoder-decoder architecture, where 1D sequence is convoluted, encoded into the latent representation, and then decoded using 2D convolutions into the Hi-C pairwise chromatin spatial proximity matrix. Those methods, while obtaining high correlation scores and improved metrics, produce Hi-C matrices that are artificial - they are blurred due to the deep learning model architecture. In our study, we propose the HiCDiffusion, sequence-only model that addresses this problem. We first train the encoder-decoder neural network and then use it as a component of the diffusion model - where we guide the diffusion using a latent representation of the sequence, as well as the final output from the encoder-decoder. That way, we obtain the high-resolution Hi-C matrices that not only better resemble the experimental results - improving the Fréchet inception distance by an average of 11 times, with the highest improvement of 56 times - but also obtain similar classic metrics to current state-of-the-art encoder-decoder architectures used for the task.

Chromatin interaction maps of human arterioles reveal new mechanisms for the genetic regulation of blood pressure.

Arterioles are small blood vessels located just upstream of capillaries in nearly all tissues. The constriction and dilation of arterioles regulate tissue perfusion and are primary determinants of systemic blood pressure (BP). Abnormalities in arterioles are central to the development of major diseases such as hypertension, stroke, and microvascular complications of diabetes. Despite the broad and essential role of arterioles in physiology and disease, current knowledge of the functional genomics of arterioles is largely absent, partly because it is challenging to obtain and analyze human arteriole samples. Here, we report extensive maps of chromatin interactions, single-cell expression, and other molecular features in human arterioles and uncover new mechanisms linking human genetic variants to gene expression in vascular cells and the development of hypertension. Compared to large arteries, arterioles exhibited a higher proportion of pericytes which were strongly associated with BP traits. BP-associated single nucleotide polymorphisms (SNPs) were enriched in chromatin interaction regions in arterioles, particularly through enhancer SNP-promoter interactions, which were further linked to gene expression specificity across tissue components and cell types. Using genomic editing in animal models and human induced pluripotent stem cells, we discovered novel mechanisms linking BP-associated noncoding SNP rs1882961 to gene expression through long-range chromatin contacts and revealed remarkable effects of a 4-bp noncoding genomic segment on hypertension in vivo. We anticipate that our rich data and findings will advance the study of the numerous diseases involving arterioles. Moreover, our approach of integrating chromatin interaction mapping in trait-relevant tissues with SNP analysis and in vivo and in vitro genome editing can be applied broadly to bridge the critical gap between genetic discoveries and physiological understanding.

Drought-responsive dynamics of H3K9ac-marked 3D chromatin interactions are integrated by OsbZIP23-associated super-enhancer-like promoter regions in rice

BackgroundIn response to drought stress (DS), plants undergo complex processes that entail significant transcriptome reprogramming. However, the intricate relationship between the dynamic alterations in the three-dimensional (3D) genome and the modulation of gene co-expression in drought responses remains a relatively unexplored area.ResultsIn this study, we reconstruct high-resolution 3D genome maps based on genomic regions marked by H3K9ac, an active histone modification that dynamically responds to soil water variations in rice. We discover a genome-wide disconnection of 3D genome contact upon DS with over 10,000 chromatin loops lost, which are partially recovered in the subsequent re-watering. Loops integrating promoter–promoter interactions (PPI) contribute to gene expression in addition to basal H3K9ac modifications. Moreover, H3K9ac-marked promoter regions with high affinities in mediating PPIs, termed as super-promoter regions (SPRs), integrate spatially clustered PPIs in a super-enhancer-like manner. Interestingly, the knockout mutation of OsbZIP23, a well-defined DS-responsive transcription factor, leads to the disassociation of over 80% DS-specific PPIs and decreased expression of the corresponding genes under DS. As a case study, we show how OsbZIP23 integrates the PPI cluster formation and the co-expression of four dehydrin genes, RAB16A–D, through targeting the RAB16C SPR in a stress signaling-dependent manner.ConclusionsOur high-resolution 3D genome maps unveil the principles and details of dynamic genome folding in response to water supply variations and illustrate OsbZIP23 as an indispensable integrator of the yet unique 3D genome organization that is essential for gene co-expression under DS in rice.

Temporally distinct 3D multi-omic dynamics in the developing human brain.

The human hippocampus and prefrontal cortex play critical roles in learning and cognition1,2, yet the dynamic molecular characteristics of their development remain enigmatic. Here we investigated the epigenomic and three-dimensional chromatin conformational reorganization during the development of the hippocampus and prefrontal cortex, using more than 53,000 joint single-nucleus profiles of chromatin conformation and DNA methylation generated bysingle-nucleus methyl-3C sequencing (snm3C-seq3)3. The remodelling of DNA methylation is temporally separated from chromatin conformation dynamics. Using single-cell profiling and multimodal single-molecule imaging approaches, we have found that short-range chromatin interactions are enriched in neurons, whereas long-range interactions are enriched in glial cells and non-brain tissues. We reconstructed the regulatory programs of cell-type development and differentiation, finding putatively causal common variants for schizophrenia strongly overlapping with chromatin loop-connected, cell-type-specific regulatory regions. Our data provide multimodal resources for studying gene regulatory dynamics in brain development and demonstrate that single-cell three-dimensional multi-omics is a powerful approach for dissecting neuropsychiatric risk loci.

8574 Single cell multiomics suggests tumor cell energy metabolism heterogeneity in non-functioning gonadotroph pituitary adenomas

Abstract Disclosure: F.M. Ruf-Zamojski: None. Z. Zhang: None. W. Cheng: None. G.R. Smith: None. M. Zamojski: None. N. Mendelev: None. Y. Ge: None. D. Marrero-Rodríguez: None. K. Taniguchi-Ponciano: None. C.L. Andoniadou: None. E. Zaslavsky: None. X. Chen: None. O.G. Troyanskaya: None. M. Mercado: None. S.C. Sealfon: None. Pituitary adenomas cause hormonal dysregulation and severe morbidity. The absence of clear biomarkers for prognosis/treatment and the high risk of recurrence make management challenging. Non-functioning gonadotroph pituitary adenomas (NFPAs) represent proliferation of gonadotroph lineage cells without an increase in secretion of gonadotropins. Detection is often delayed until mass effects cause visual defects. Higher resolution molecular study may improve classification, diagnosis, treatment, and provide insight into adenoma biology. We analyzed four NFPAs using unbiased genome-wide same-cell single nucleus 10X Genomics multiome assay (transcriptome and chromatin accessibility) to define individual cell identity/states in comparison with normal gonadotrope cell references from twelve age/sex-matched human post-mortem pituitaries. In total, we generated high quality multiomics single-cell data for 33,832 adenoma cells and 5,451 normal gonadotropes. We identified tumor, immune, vascular, and proliferative cells confirming cellular heterogeneity. Multidimensional scaling analysis using Manhattan distances between cell types showed that NFPA gonadotropes were closer in gene expression to healthy gonadotropes among all pituitary cell types. In addition, we observed differences in the macrophage and endothelial cell populations between tumors and healthy pituitaries. Next, we studied gene regulatory circuitry (regulated gene, modulating transcription factor and its chromatin interaction site) using the new Control of Regulation Extracted from Multiomics Assays (CREMA, [1]) multiome analysis method which showed high variation between normal and tumor gonadotropes. Within each adenoma, the NFPA tumor cells encompassed two main subgroups of cells, one group (GT+ marker tumor cells) expressing established gonadotrope markers and another group not (GT- marker tumor cells). GT+ and GT- tumor cells formed a cluster distinct from healthy gonadotropes. GT+ and GT- tumor cells differed by their level of expression of mitochondrial genes, suggesting decreased aerobic capacity associated with the GT- tumor subgroup. Interestingly, differential analyses between the healthy and the NFPA gonadotropes from same-cell multiome data highlighted dissimilar gene expression and chromatin accessibility patterns.Our newly collected single cell data characterize the difference between NFPA tumor cells and normal gonadotrope cells. Our results suggest a Warburg-effect like difference in aerobic respiratory capacity among tumor cells within the same pituitary adenoma. Overall, this study brings new insights into the molecular characteristics of NFPAs, including tumor cell heterogeneity, and represents a new resource for the field. [1] PMID: 38084892. Presentation: 6/1/2024

Predicting p53-dependent cell transitions from thermodynamic models.

A cell's fate involves transitions among its various states, each defined by a distinct gene expression profile governed by the topology of gene regulatory networks, which are affected by 3D genome organization. Here, we develop thermodynamic models to determine the fate of a malignant cell as governed by the tumor suppressor p53 signaling network, taking into account long-range chromatin interactions in the mean-field approximation. The tumor suppressor p53 responds to stress by selectively triggering one of the potential transcription programs that influence many layers of cell signaling. These range from p53 phosphorylation to modulation of its DNA binding affinity, phase separation phenomena, and internal connectivity among cell fate genes. We use the minimum free energy of the system as a fundamental property of biological networks that influences the connection between the gene network topology and the state of the cell. We constructed models based on network topology and equilibrium thermodynamics. Our modeling shows that the binding of phosphorylated p53 to promoters of target genes can have properties of a first order phase transition. We apply our model to cancer cell lines ranging from breast cancer (MCF-7), colon cancer (HCT116), and leukemia (K562), with each one characterized by a specific network topology that determines the cell fate. Our results clarify the biological relevance of these mechanisms and suggest that they represent flexible network designs for switching between developmental decisions.

Rational fragment-based design of compounds targeting the PWWP domain of the HRP family

Lens epithelium-derived growth factor p75 (LEDGF/p75), member of the hepatoma-derived growth-factor-related protein (HRP) family, is a transcriptional co-activator and involved in several pathologies including HIV infection and malignancies such as MLL-rearranged leukemia. LEDGF/p75 acts by tethering proteins to the chromatin through its integrase binding domain. This chromatin interaction occurs between the PWWP domain of LEDGF/p75 and nucleosomes carrying a di- or trimethylation mark on histone H3 Lys36 (H3K36me2/3). Our aim is to rationally devise small molecule drugs capable of inhibiting such interaction. To bootstrap this development, we resorted to X-ray crystallography-based fragment screening (FBS-X). Given that the LEDGF PWWP domain crystals were not suitable for FBS-X, we employed crystals of the closely related PWWP domain of paralog HRP-2. As a result, as many as 68 diverse fragment hits were identified, providing a detailed sampling of the H3K36me2/3 pocket pharmacophore. Subsequent structure-guided fragment expansion in three directions yielded multiple compound series binding to the pocket, as verified through X-ray crystallography, nuclear magnetic resonance and differential scanning fluorimetry. Our best compounds have double-digit micromolar affinity and optimally sample the interactions available in the pocket, judging by the Kd-based ligand efficiency exceeding 0.5 kcal/mol per non-hydrogen atom. Beyond π-stacking within the aromatic cage of the pocket and hydrogen bonding, the best compounds engage in a σ-hole interaction between a halogen atom and a conserved water buried deep in the pocket. Notably, the binding pocket in LEDGF PWWP is considerably smaller compared to the related PWWP1 domains of NSD2 and NSD3 which feature an additional subpocket and for which nanomolar affinity compounds have been developed recently. The absence of this subpocket in LEDGF PWWP limits the attainable affinity. Additionally, these structural differences in the H3K36me2/3 pocket across the PWWP domain family translate into a distinct selectivity of the compounds we developed. Our top ranked compounds are interacting with both homologous LEDGF and HRP-2 PWWP domains, yet they showed no affinity for the NSD2 PWWP1 and BRPF2 PWWP domains which belong to other PWWP domain subfamilies. Nevertheless, our developed compound series provide a strong foundation for future drug discovery targeting the LEDGF PWWP domain as they can further be explored through combinatorial chemistry. Given that the affinity of H3K36me2/3 nucleosomes to LEDGF/p75 is driven by interactions within the pocket as well as with the DNA-binding residues, we suggest that future compound development should target the latter region as well. Beyond drug discovery, our compounds can be employed to devise tool compounds to investigate the mechanism of LEDGF/p75 in epigenetic regulation.

Peculiar k-mer Spectra Are Correlated with 3D Contact Frequencies and Breakpoint Regions in the Human Genome.

It is widely accepted that the 3D chromatin organization in human cell nuclei is not random and recent investigations point towards an interactive relation of epigenetic functioning and chromatin (re-)organization. Although chromatin organization seems to be the result of self-organization of the entirety of all molecules available in the cell nucleus, a general question remains open as to what extent chromatin organization might additionally be predetermined by the DNA sequence and, if so, if there are characteristic differences that distinguish typical regions involved in dysfunction-related aberrations from normal ones, since typical DNA breakpoint regions involved in disease-related chromosome aberrations are not randomly distributed along the DNA sequence. Highly conserved k-mer patterns in intronic and intergenic regions have been reported in eukaryotic genomes. In this article, we search and analyze regions deviating from average spectra (ReDFAS) of k-mer word frequencies in the human genome. This includes all assembled regions, e.g., telomeric, centromeric, genic as well as intergenic regions. A positive correlation between k-mer spectra and 3D contact frequencies, obtained exemplarily from given Hi-C datasets, has been found indicating a relation of ReDFAS to chromatin organization and interactions. We also searched and found correlations of known functional annotations, e.g., genes correlating with ReDFAS. Selected regions known to contain typical breakpoints on chromosomes 9 and 5 that are involved in cancer-related chromosomal aberrations appear to be enriched in ReDFAS. Since transposable elements like ALUs are often assigned as major players in 3D genome organization, we also studied their impact on our examples but could not find a correlation between ALU regions and breakpoints comparable to ReDFAS. Our findings might show that ReDFAS are associated with instable regions of the genome and regions with many chromatin contacts which is in line with current research indicating that chromatin loop anchor points lead to genomic instability.

Virus Infection Induces Immune Gene Activation with CTCF Anchored Enhancers and Chromatin Interactions in Pig Genome.

Chromatin organization is important for gene transcription in pig genome. However, its three-dimensional (3D) structure and dynamics are much less investigated than those in human. Here we applied the long-reads chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) method to map the whole-genome chromatin interactions mediated by CCCTC-binding factor (CTCF) and RNA polymerase Ⅱ (RNAPⅡ or POLⅡ) in porcine macrophage cells before and after polyinosinic-polycytidylic acid [Poly(I:C)] induction. Our results revealed that Poly(I:C) induction impacts the 3D genome organization in the 3D4/21 cells at the fine-scale chromatin loop level rather than at the large-scale domain level. Furthermore, our findings underscored the pivotal role of CTCF anchored chromatin interactions in reshaping chromatin architecture during immune responses. Knock-out of the CTCF locus further confirmed that the CTCF anchored enhancers are associated with the activation of immune genes via long-range interactions. Notably, ChIA-PET data also supported the spatial relationship between single nucleotide polymorphisms (SNPs) and the related gene transcription in 3D genome aspect. Our findings in this study provide new clues and potential targets to explore key elements related to diseases in swine and are also likely to shed light on elucidating chromatin organization and dynamics underlying the process of mammalian infectious diseases.

Prediction and functional interpretation of inter-chromosomal genome architecture from DNA sequence with TwinC.

Three-dimensional nuclear DNA architecture comprises well-studied intra-chromosomal (cis) folding and less characterized inter-chromosomal (trans) interfaces. Current predictive models of 3D genome folding can effectively infer pairwise cis-chromatin interactions from the primary DNA sequence but generally ignore trans contacts. There is an unmet need for robust models of trans-genome organization that provide insights into their underlying principles and functional relevance. We present TwinC, an interpretable convolutional neural network model that reliably predicts trans contacts measurable through genome-wide chromatin conformation capture (Hi-C). TwinC uses a paired sequence design from replicate Hi-C experiments to learn single base pair relevance in trans interactions across two stretches of DNA. The method achieves high predictive accuracy (AUROC=0.80) on a cross-chromosomal test set from Hi-C experiments in heart tissue. Mechanistically, the neural network learns the importance of compartments, chromatin accessibility, clustered transcription factor binding and G-quadruplexes in forming trans contacts. In summary, TwinC models and interprets trans genome architecture, shedding light on this poorly understood aspect of gene regulation.

Super-silencer perturbation by EZH2 and REST inhibition leads to large loss of chromatin interactions and reduction in cancer growth

Super-silencer perturbation by EZH2 and REST inhibition leads to large loss of chromatin interactions and reduction in cancer growth

TFEB controls expression of human syncytins during cell-cell fusion.

During human development, a temporary organ is formed, the placenta, which invades the uterine wall to support nutrient, oxygen, and waste exchange between the mother and fetus until birth. Most of the human placenta is formed by a syncytial villous structure lined by syncytialized trophoblasts, a specialized cell type that forms via cell-cell fusion of underlying progenitor cells. Genetic and functional studies have characterized the membrane protein fusogens Syncytin-1 and Syncytin-2, both of which are necessary and sufficient for human trophoblast cell-cell fusion. However, identification and characterization of upstream transcriptional regulators regulating their expression have been limited. Here, using CRISPR knockout in an in vitro cellular model of syncytiotrophoblast development (BeWo cells), we found that the transcription factor TFEB, mainly known as a regulator of autophagy and lysosomal biogenesis, is required for cell-cell fusion of syncytiotrophoblasts. TFEB translocates to the nucleus, exhibits increased chromatin interactions, and directly binds the Syncytin-1 and Syncytin-2 promoters to control their expression during differentiation. Although TFEB appears to play a critical role in syncytiotrophoblast differentiation, ablation of TFEB largely does not affect lysosomal gene expression or lysosomal biogenesis in differentiating BeWo cells, suggesting a previously uncharacterized role for TFEB in controlling the expression of human syncytins.

Methods for Genome-Wide Chromatin Interaction Analysis.

Recent analyses revealed the essential function of chromatin structure in maintaining and regulating genomic information. Advancements in microscopy, nuclear structure observation techniques, and the development of methods utilizing next-generation sequencers (NGSs) have significantly progressed these discoveries. Methods utilizing NGS enable genome-wide analysis, which is challenging with microscopy, and have elucidated concepts of important chromatin structures such as a loop structure, a domain structure called topologically associating domains (TADs), and compartments. In this chapter, I introduce chromatin interaction techniques using NGS and outline the principles and features of each method.

Abstract C058: A novel approach to understanding pancreas cancer risk through enhancer- promoter interactions

Abstract Genetic or epigenetic variations in regulatory enhancer elements increase susceptibility to a range of pathologies, including pancreas cancer. Pinpointing genes affecting pancreas cancer risk holds promise for early detection, prevention, and effective therapies. Despite recent advances, linking enhancer elements to target genes and predicting the transcriptional outcome of enhancer dysfunction remain significant challenges. Using 3D chromatin assays, we generated an extensive enhancer interaction dataset for the human pancreas, spanning more than 20 donors and five major cell types, including both the exocrine and endocrine compartments. We employed a network approach to parse chromatin interactions into enhancer-promoter tree models, facilitating quantitative, genome-wide analysis of enhancer connectivity. Using the tree models, we developed a machine learning algorithm capable of estimating the impact of enhancer perturbations on cell-type specific gene expression in the human pancreas. Complementing this computational approach, we streamlined an experimental platform using CRISPR-interference and RNA-FISH coupled with high-throughput imaging to quantitatively measure the effect size of enhancer perturbations in primary human pancreas cells. Our enhancer tree models enabled functional annotation of genetic variants associated with pancreas diseases, revealing novel target genes in specific cell types. For pancreas cancer (PDAC), our data suggests a stronger association of disease susceptibility variants with acinar cells, even though ductal cells historically have been the focus of PDAC. This discovery linked three potentially functional variants associated with PDAC risk to acinar cell enhancers, revealing novel target genes to further investigation. By integrating genomics with mechanistic findings into enhancer function, our work offers a framework for understanding genetic basis of pancreas cancer disease risk. Citation Format: Li Wang, Songjoon Baek, Gauri Prasad, Thucnhi Truongvo, Jason Hoskins, Laufey Amundadottir, Efsun Arda. A novel approach to understanding pancreas cancer risk through enhancer- promoter interactions [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Advances in Pancreatic Cancer Research; 2024 Sep 15-18; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2024;84(17 Suppl_2):Abstract nr C058.

Chromatin conformation capture in the clinic: 4C-seq/HiC distinguishes pathogenic from neutral duplications at the GPR101 locus

BackgroundX-linked acrogigantism (X-LAG; MIM: 300942) is a severe form of pituitary gigantism caused by chromosome Xq26.3 duplications involving GPR101. X-LAG-associated duplications disrupt the integrity of the topologically associating domain (TAD) containing GPR101 and lead to the formation of a neo-TAD that drives pituitary GPR101 misexpression and gigantism. As X-LAG is fully penetrant and heritable, duplications involving GPR101 identified on prenatal screening studies, like amniocentesis, can pose an interpretation challenge for medical geneticists and raise important concerns for patients and families. Therefore, providing robust information on the functional genomic impact of such duplications has important research and clinical value with respect to gene regulation and triplosensitivity traits.MethodsWe employed 4C/HiC-seq as a clinical tool to determine the functional impact of incidentally discovered GPR101 duplications on TAD integrity in three families. After defining duplications and breakpoints around GPR101 by clinical-grade and high-density aCGH, we constructed 4C/HiC chromatin contact maps for our study population and compared them with normal and active (X-LAG) controls.ResultsWe showed that duplications involving GPR101 that preserved the centromeric invariant TAD boundary did not generate a pathogenic neo-TAD and that ectopic enhancers were not adopted. This allowed us to discount presumptive/suspected X-LAG diagnoses and GPR101 misexpression, obviating the need for intensive clinical follow-up.ConclusionsThis study highlights the importance of TAD boundaries and chromatin interactions in determining the functional impact of copy number variants and provides proof-of-concept for using 4C/HiC-seq as a clinical tool to acquire crucial information for genetic counseling and to support clinical decision-making in cases of suspected TADopathies.

Predictive Prioritization of Enhancers Associated with Pancreas Disease Risk.

Genetic and epigenetic variations in regulatory enhancer elements increase susceptibility to a range of pathologies. Despite recent advances, linking enhancer elements to target genes and predicting transcriptional outcomes of enhancer dysfunction remain significant challenges. Using 3D chromatin conformation assays, we generated an extensive enhancer interaction dataset for the human pancreas, encompassing more than 20 donors and five major cell types, including both exocrine and endocrine compartments. We employed a network approach to parse chromatin interactions into enhancer-promoter tree models, facilitating a quantitative, genome-wide analysis of enhancer connectivity. With these tree models, we developed a machine learning algorithm to estimate the impact of enhancer perturbations on cell type-specific gene expression in the human pancreas. Orthogonal to our computational approach, we perturbed enhancer function in primary human pancreas cells using CRISPR interference and quantified the effects at the single-cell level through RNA FISH coupled with high-throughput imaging. Our enhancer tree models enabled the annotation of common germline risk variants associated with pancreas diseases, linking them to putative target genes in specific cell types. For pancreatic ductal adenocarcinoma, we found a stronger enrichment of disease susceptibility variants within acinar cell regulatory elements, despite ductal cells historically being assumed as the primary cell-of-origin. Our integrative approach-combining cell type-specific enhancer-promoter interaction mapping, computational models, and single-cell enhancer perturbation assays-produced a robust resource for studying the genetic basis of pancreas disorders.