Hi-C Analysis Research Articles

CFDP1 is a neuroblastoma susceptibility gene that regulates transcription factors of the noradrenergic cell identity.

Pleiotropic genetic factors (e.g., DNA polymorphisms) may be involved in the initiation of neuroblastoma (NB) and coronary artery disease (CAD) given their common origin from defects in neural crest development. To discover novel NB susceptibility genes, we conducted a three-stage survey including a meta-analysis of NB and CAD genome-wide association data, prioritization of NB causal variants, and validation in an independent cohort of affected individuals-control subjects. The lead SNP, rs13337397 at the 16q23.1 locus, associated with both diseases in the meta-analysis and with NB in the validation study. All the SNPs in linkage disequilibrium with rs13337397 were annotated using the H3K27ac epigenetic marker of neural crest cells (NCC) and NB cell lines. Indeed, we identified the functional SNP rs13337017, mapping within an enhancer of NCCs and NB cell lines and showing long-range interactions with CFDP1 by Hi-C analysis. Luciferase assays indicated that the risk allele of rs13337017 increased CFDP1 expression in NB cell lines. Of note, CFDP1 high expression associated with unfavorable prognostic markers in an analysis including 498 NB transcriptomes. Moreover, depletion of CFDP1 markedly decreased viability and migration and increased apoptotic rates in NB cell lines. Finally, transcriptome and qPCR analyses revealed that the depletion of CFDP1 may affect noradrenergic neuron differentiation by downregulating master regulators of sympathetic noradrenergic identity, including PHOX2B, HAND2, and GATA3. Our data strongly suggest that CFDP1 acts as oncogene in NB. In addition, we provide evidence that genetic predisposition to NB can be mediated by the alteration of noradrenergic lineage-specific gene expression.

Author Correction: Chromosome-level genomes of seeded and seedless date plum based on third-generation DNA sequencing and Hi-C analysis

Author Correction: Chromosome-level genomes of seeded and seedless date plum based on third-generation DNA sequencing and Hi-C analysis

HNF1B-driven three-dimensional chromatin structure for molecular classification in pancreatic cancers.

The molecular subtypes of pancreatic cancer (PC), either classical/progenitor-like or basal/squamous-like, are currently a major topic of research because of their direct association with clinical outcomes. Some transcription factors (TFs) have been reported to be associated with these subtypes. However, the mechanisms by which these molecular signatures of PCs are established remain unknown. Epigenetic regulatory processes, supported by dynamic changes in the chromatin structure, are essential for transcriptional profiles. Previously, we reported the importance of open chromatin profiles in the biological features and transcriptional status of PCs. Here, we aimed to analyze the relationships between three-dimensional (3D) genome structures and the molecular subtypes of human PCs using Hi-C analysis. We observed a correlation of the specific elements of 3D genome modules, including compartments, topologically associating domains, and enhancer-promoter loops, with the expression of related genes. We focused on HNF1B, a TF that is implicated in the progenitor subtype. Forced expression of HNF1B in squamous-type PC organoids induced the upregulation and downregulation of genes associated with progenitor and squamous subtypes, respectively. Long-range genomic interactions induced by HNF1B were accompanied by compartment modulation and H3K27ac redistribution. We also found that these HNF1B-induced changes in subtype-related gene expression required an intrinsically disordered region, suggesting a possible involvement of phase separation in compartment modulation. Thus, mapping of 3D structural changes induced by TFs, such as HNF1B, may become a useful resource for further understanding the molecular features of PCs.

Multiscale reorganization of the genome following DNA damage facilitates chromosome translocations via nuclear actin polymerization

Nuclear actin-based movements have been shown to orchestrate clustering of DNA double-strand breaks (DSBs) into homology-directed repair domains. Here we describe multiscale three-dimensional genome reorganization following DNA damage and analyze the contribution of the nuclear WASP-ARP2/3-actin pathway toward chromatin topology alterations and pathologic repair. Hi-C analysis reveals genome-wide, DNA damage-induced chromatin compartment flips facilitated by ARP2/3 that enrich for open, A compartments. Damage promotes interactions between DSBs, which in turn facilitate aberrant, actin-dependent intra- and inter-chromosomal rearrangements. Our work establishes that clustering of resected DSBs into repair domains by nuclear actin assembly is coordinated with multiscale alterations in genome architecture that enable homology-directed repair while also increasing nonhomologous end-joining-dependent translocation frequency.

Chromosome arm length, and a species-specific determinant, define chromosome arm width.

Mitotic chromosomes in different organisms adopt various dimensions. What defines these dimensions is scarcely understood. Here, we compare mitotic chromosomes in budding and fission yeasts harboring similarly sized genomes distributed among 16 or 3 chromosomes, respectively. Hi-C analyses and superresolution microscopy reveal that budding yeast chromosomes are characterized by shorter-ranging mitotic chromatin contacts and are thinner compared with the thicker fission yeast chromosomes that contain longer-ranging mitotic contacts. These distinctions persist even after budding yeast chromosomes are fused to form three fission-yeast-length entities, revealing a species-specific organizing principle. Species-specific widths correlate with the known binding site intervals of the chromosomal condensin complex. Unexpectedly, within each species, we find that longer chromosome arms are always thicker and harbor longer-ranging contacts, a trend that we also observe with human chromosomes. Arm length as a chromosome width determinant informs mitotic chromosome formation models.

Deep generative modeling and clustering of single cell Hi-C data.

Deciphering 3D genome conformation is important for understanding gene regulation and cellular function at a spatial level. The recent advances of single cell Hi-C technologies have enabled the profiling of the 3D architecture of DNA within individual cell, which allows us to study the cell-to-cell variability of 3D chromatin organization. Computational approaches are in urgent need to comprehensively analyze the sparse and heterogeneous single cell Hi-C data. Here, we proposed scDEC-Hi-C, a new framework for single cell Hi-C analysis with deep generative neural networks. scDEC-Hi-C outperforms existing methods in terms of single cell Hi-C data clustering and imputation. Moreover, the generative power of scDEC-Hi-C could help unveil the differences of chromatin architecture across cell types. We expect that scDEC-Hi-C could shed light on deepening our understanding of the complex mechanism underlying the formation of chromatin contacts.

Turnover of mammal sex chromosomes in the Sry-deficient Amami spiny rat is due to male-specific upregulation of Sox9

Mammalian sex chromosomes are highly conserved, and sex is determined by SRY on the Y chromosome. Two exceptional rodent groups in which some species lack a Y chromosome and Sry offer insights into how novel sex genes can arise and replace Sry, leading to sex chromosome turnover. However, intensive study over three decades has failed to reveal the identity of novel sex genes in either of these lineages. We here report our discovery of a male-specific duplication of an enhancer of Sox9 in the Amami spiny rat Tokudaia osimensis, in which males and females have only a single X chromosome (XO/XO) and the Y chromosome and Sry are completely lost. We performed a comprehensive survey to detect sex-specific genomic regions in the spiny rat. Sex-related genomic differences were limited to a male-specific duplication of a 17-kb unit located 430 kb upstream of Sox9 on an autosome. Hi-C analysis using male spiny rat cells showed the duplicated region has potential chromatin interaction with Sox9. The duplicated unit harbored a 1,262-bp element homologous to mouse enhancer 14 (Enh14), a candidate Sox9 enhancer that is functionally redundant in mice. Transgenic reporter mice showed that the spiny rat Enh14 can function as an embryonic testis enhancer in mice. Embryonic gonads of XX mice in which Enh14 was replaced by the duplicated spiny rat Enh14 showed increased Sox9 expression and decreased Foxl2 expression. We propose that male-specific duplication of this Sox9 enhancer substituted for Sry function, defining a novel Y chromosome in the spiny rat.

Abstract A069: Hi-C analysis in patient-derived organoids reveals the impact of three-dimensional chromatin structures on pancreatic cancer transcriptional subtypes

Abstract Three-dimensional (3D) chromatin structures connect genes to regulatory elements, thereby contributing to cell-type-specific gene regulation. However, its role in pancreatic cancer remains unclear. Inspired by our previous finding that open chromatin profiles analyzed by ATAC-seq well characterize the biological features and transcriptional status of pancreatic neoplasms, we hereby demonstrate the previously overlooked impact of 3D epigenetic regulation in pancreatic cancers. Through integrative analysis of Hi-C, ChIP-seq, and RNA-seq using patient-derived organoids, 3D genome modules, including compartments, topologically associating domains, and enhancer-promoter loops were found to play important roles in characterizing pancreatic tumors. Interestingly, intraductal papillary mucinous neoplasm (IPMN) lineages tended to have unique 3D chromatin structure profiles. We further found that transcriptional subtypes in pancreatic cancer, either classical/progenitor-like or basal/squamous-like, were largely underpinned by 3D structures. Compartment was responsible for supporting many subtype-specific genes, by synchronously activating genes that were in the vicinity, as represented by CYP2C family genes. Meanwhile, enhancer-promoter loop fine-tuned some key subtype genes, such as GATA6. By exploring the driving forces behind these 3D structures, we found that HNF1B, a transcription factor implicated in the progenitor subtype, was sufficient to induce subtype switch, together with their underlying 3D genome structure alterations. We demonstrate that HNF1B overexpression in squamous-type pancreatic cancer organoids induced the upregulation and downregulation of genes associated with the progenitor and squamous subtype, respectively. Mechanistically, these dynamic transcriptional and 3D structural modulation were accompanied by HNF1B-induced long-range genomic interactions and H3K27ac modification at HNF1B bound sites, implying HNF1B as a driving force. We also provide evidence suggesting that this HNF1B-induced subtype switch was mediated through the phase-separating property of HNF1B, which was enabled by its intrinsically disordered region (IDR). Induction of mutant HNF1B, whose IDR was either deleted (HNF1B-dIDR) or substituted by transcriptional phase-separation forming IDR from another peptide AKAP95 (HNF1B-AKAP95), showed that subtype switching effects of HNF1B were diminished by HNF1B-dIDR or recapitulated by HNF1B-AKAP95. Thus, HNF1B contributes to the transcriptional subtype formation in pancreatic cancers through its phase-separating structural ability. Collectively, the mapping of 3D structural changes induced by transcription factors, such as HNF1B, is useful for further understanding the features of pancreatic cancers. Our study may imply the potential of 3D structures to effectively suppress or unwind the malignant transformation process and to discover subtype-specific therapeutic targets. Citation Format: Hiroyuki Kato, Dosuke Iwadate, Keisuke Yamamoto, Hiroaki Fujiwara, Hideaki Ijichi, Mitsuhiro Fujishiro, Keisuke Tateishi. Hi-C analysis in patient-derived organoids reveals the impact of three-dimensional chromatin structures on pancreatic cancer transcriptional subtypes [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr A069.

EPCO-23. DYSFUNCTION OF LARGE H3K9ME3 DOMAINS IN ATRX DEFICIENT GLIOMAS INDUCES GENETIC REARRANGEMENTS AND LATENT DEVELOPMENTAL SIGNALING NETWORKS THROUGH SUPER-ENHANCER LANDSCAPES

Abstract Somatic alterations in ATRX (Alpha Thalassemia/Mental Retardation Syndrome X-linked), a member of SWI/SNF family chromatin regulator has been found to be frequently mutated in diffuse gliomas and defines molecular subtypes with aggressive behavior. Mechanistically, ATRX regulates incorporation of histone H3.3 into chromatin sites across the genome, maintains alternative lengthening of telomeres and establishes genomic distribution of polycomb responsive genes. To understand in depth role of ATRX in gliomas, we performed ChIP-seq and/or Cut-and-tag for histone marks that define active transcription, enhancers, repressors, gene bodies and CTCF on Atrx intact and deficient mNPCs (mouse Neural Progenitor Cells). Our integrated analysis reports depletion of H3K9me3 loci’s that coincidently enriched with Atrx binding sites and Lamina-Associated Domains with genes enriched for cell cycle, motility and chromosome organization. This chromatin perturbations at heterochromatin domains was further confirmed in our Hi-C analysis with switching of A-B and B-A compartments, reorganization of TADs with occasional CTCF marks and gain of novel interacting loops that showed gene expression leakage required for gliomagenesis using Capture Hi-C. Notably, we observed aberrant levels of endogenous retroviral elements (ERVs) family members, including upregulation of Line-1 elements in mNPCs and patient derived glioma stem cells (GSCs) and our analysis shows increased copy number variations in ATRX deficient gliomas as a consequence of Line-1 activation in these subsets of tumors. Finally, our integrated omics- analysis demonstrates enrichment of super-enhancers in Atrx deficient background with several putative druggable candidates for clinical benefits. As an example, we show presence of H3K27ac super-enhancer over HOXA locus and targeting pan-HOXA genes using small inhibitor peptides diminished proliferation and migration of mNPCs and GSCs with increased in cell apoptosis with alterations in downstream developmental signaling pathways. To summarize,our data establishes tangible links between Atrx deficiency and multiscale dysregulated cellular phenotype in gliomas with identifying novel targets for therapy.

STEM-24. SINGLE-CELL ATLAS OF HUMAN FETAL CEREBELLUM DECODING MEDULLOBLASTOMA ORIGIN AND ONCOGENESIS

Abstract Medulloblastomas (MBs) are the most common malignant childhood brain tumors, yet the origin of the most aggressive subgroup-3 form remains elusive, impeding development of effective targeted treatment strategies. Previous cell-type analyses of mouse cerebella or human counterparts from frozen tissue nuclei have not fully defined the compositional heterogeneity of MBs. Here, we undertook an unprecedented single-cell profiling of freshly-isolated human fetal cerebella at different developmental stages to establish a reference map for delineating the hierarchical cellular states in MBs. We identified a unique transitional cerebellar progenitor connecting neural stem cells to neuronal lineages in human fetal cerebella. Intersectional analysis revealed that the transitional progenitors were enriched in aggressive MB subgroups, including group-3 and metastatic tumors. Integrated single-cell multi-omic profiling revealed unique regulatory networks in the transitional progenitor populations, including transcriptional determinants HNRNPH1 and SOX11, which are correlated with clinical prognosis in aggressive group-3 MBs. Genomic profiling and Hi-C analyses identified de novo long-range chromatin loops juxtaposing HNRNPH1/SOX11-targeted super-enhancers to cis-regulatory elements of MYC, an oncogenic driver for group-3 MBs. Targeting the transitional progenitor regulators inhibited MYC expression and MYC-driven group-3 MB growth. Together, our integrated single-cell atlases of human fetal cerebella and MBs reveal important cell populations predisposed to transformation and regulatory circuitries underlying tumor cell state evolution and oncogenesis, highlighting hitherto unrecognized transitional progenitor intermediates predictive of disease prognosis and potential therapeutic vulnerabilities.

DcHiC detects differential compartments across multiple Hi-C datasets

The compartmental organization of mammalian genomes and its changes play important roles in distinct biological processes. Here, we introduce dcHiC, which utilizes a multivariate distance measure to identify significant changes in compartmentalization among multiple contact maps. Evaluating dcHiC on four collections of bulk and single-cell contact maps from in vitro mouse neural differentiation (n = 3), mouse hematopoiesis (n = 10), human LCLs (n = 20) and post-natal mouse brain development (n = 3 stages), we show its effectiveness and sensitivity in detecting biologically relevant changes, including those orthogonally validated. dcHiC reported regions with dynamically regulated genes associated with cell identity, along with correlated changes in chromatin states, subcompartments, replication timing and lamin association. With its efficient implementation, dcHiC enables high-resolution compartment analysis as well as standalone browser visualization, differential interaction identification and time-series clustering. dcHiC is an essential addition to the Hi-C analysis toolbox for the ever-growing number of bulk and single-cell contact maps. Available at: https://github.com/ay-lab/dcHiC.

Depletion or cleavage of cohesin during anaphase differentially affects chromatin structure and segregation.

Chromosome segregation requires both the separation of sister chromatids and the sustained condensation of chromatids during anaphase. In yeast cells, cohesin is not only required for sister chromatid cohesion but also plays a major role determining the structure of individual chromatids in metaphase. Separase cleavage is thought to remove all cohesin complexes from chromosomes to initiate anaphase. It is thus not clear how the length and organisation of segregating chromatids is maintained during anaphase in the absence of cohesin. Here, we show that degradation of cohesin at the anaphase onset causes aberrant chromatid segregation. Hi-C analysis on segregating chromatids demonstrates that cohesin depletion causes loss of intrachromatid organisation. Surprisingly, tobacco etch virus (TEV)-mediated cleavage of cohesin does not dramatically disrupt chromatid organisation in anaphase, explaining why bulk segregation is achieved. In addition, we identified a small pool of cohesin complexes bound to telophase chromosomes in wild-type cells and show that they play a role in the organisation of centromeric regions. Our data demonstrates that in yeast cells cohesin function is not over in metaphase, but extends to the anaphase period when chromatids are segregating.

Decoding the spatial chromatin organization and dynamic epigenetic landscapes of macrophage cells during differentiation and immune activation

Immunocytes dynamically reprogram their gene expression profiles during differentiation and immunoresponse. However, the underlying mechanism remains elusive. Here, we develop a single-cell Hi-C method and systematically delineate the 3D genome and dynamic epigenetic atlas of macrophages during these processes. We propose “degree of disorder” to measure genome organizational patterns inside topologically-associated domains, which is correlated with the chromatin epigenetic states, gene expression, and chromatin structure variability in individual cells. Furthermore, we identify that NF-κB initiates systematic chromatin conformation reorganization upon Mycobacterium tuberculosis infection. The integrated Hi-C, eQTL, and GWAS analysis depicts the atlas of the long-range target genes of mycobacterial disease susceptible loci. Among these, the SNP rs1873613 is located in the anchor of a dynamic chromatin loop with LRRK2, whose inhibitor AdoCbl could be an anti-tuberculosis drug candidate. Our study provides comprehensive resources for the 3D genome structure of immunocytes and sheds insights into the order of genome organization and the coordinated gene transcription during immunoresponse.

NODULIN HOMEOBOX is required for heterochromatin homeostasis in Arabidopsis

Arabidopsis NODULIN HOMEOBOX (NDX) is a nuclear protein described as a regulator of specific euchromatic genes within transcriptionally active chromosome arms. Here we show that NDX is primarily a heterochromatin regulator that functions in pericentromeric regions to control siRNA production and non-CG methylation. Most NDX binding sites coincide with pericentromeric het-siRNA loci that mediate transposon silencing, and are antagonistic with R-loop structures that are prevalent in euchromatic chromosomal arms. Inactivation of NDX leads to differential siRNA accumulation and DNA methylation, of which CHH/CHG hypomethylation colocalizes with NDX binding sites. Hi-C analysis shows significant chromatin structural changes in the ndx mutant, with decreased intrachromosomal interactions at pericentromeres where NDX is enriched in wild-type plants, and increased interchromosomal contacts between KNOT-forming regions, similar to those observed in DNA methylation mutants. We conclude that NDX is a key regulator of heterochromatin that is functionally coupled to het-siRNA loci and non-CG DNA methylation pathways.

A chromosome-level genome assembly provides insights into the environmental adaptability and outbreaks of Chlorops oryzae

Chlorops oryzae is a pest of rice that has caused severe damage to crops in major rice-growing areas in recent years. We generated a 447.60 Mb high-quality chromosome-level genome with contig and scaffold N50 values of 1.17 Mb and 117.57 Mb, respectively. Hi-C analysis anchored 93.22% scaffolds to 4 chromosomes. The relatively high expression level of Heat Shock Proteins (HSPs) and antioxidant genes in response to thermal stress suggests these genes may play a role in the environmental adaptability of C. oryzae. The identification of multiple pathways that regulate reproductive development (juvenile hormone, 20-hydroxyecdsone, and insulin signaling pathways) provides evidence that these pathways also play an important role in vitellogenesis and thus insect population maintenance. These findings identify possible reasons for the increased frequency of outbreaks of C. oryzae in recent years. Our chromosome-level genome assembly may provide a basis for further genetic studies of C. oryzae, and promote the development of novel, sustainable strategies to control this pest.

Super-enhancer hijacking LINC01977 promotes malignancy of early-stage lung adenocarcinoma addicted to the canonical TGF-β/SMAD3 pathway

BackgroundLung adenocarcinoma (LUAD) is the leading cause of death worldwide. However, the roles of long noncoding RNAs (lncRNAs) hijacked by super-enhancers (SEs), vital regulatory elements of the epigenome, remain elusive in the progression of LUAD metastasis.MethodsSE-associated lncRNA microarrays were used to identify the dysregulated lncRNAs in LUAD. ChIP-seq, Hi-C data analysis, and luciferase reporter assays were utilized to confirm the hijacking of LINC01977 by SE. The functions and mechanisms of LINC01977 in LUAD were explored by a series of in vitro and in vivo assays.ResultsWe found that LINC01977, a cancer-testis lncRNA, was hijacked by SE, which promoted proliferation and invasion both in vitro and in vivo. LINC01977 interacted with SMAD3 to induce its nuclear transport, which facilitated the interaction between SMAD3 and CBP/P300, thereby regulating the downstream target gene ZEB1. Additionally, SMAD3 up-regulated LINC09177 transcription by simultaneously binding the promoter and SE, which was induced by the infiltration of M2-like tumor-associated macrophages (TAM2), subsequently activating the TGF-β/SMAD3 pathway. Moreover, LINC01977 expression was positively correlated with TAM2 infiltration and SMAD3 expression, especially in early-stage LUAD. Higher chromatin accessibility in the SE region of LINC01977 was observed with high expression of TGF-β. Early-stage LUAD patients with high LIN01977 expression had a shorter disease-free survival.ConclusionsTAM2 infiltration induced a rich TGF-β microenvironment, activating SMAD3 to bind the promoter and the SE of LINC01977, which up-regulated LINC01977 expression. LINC01977 also promoted malignancy via the canonical TGF-β/SMAD3 pathway. LINC01977 hijacked by SE could be a valuable therapeutic target, especially for the treatment of early-stage LUAD.

Imputation-free reconstructions of three-dimensional chromosome architectures in human diploid single-cells using allele-specified contacts

Single-cell Hi-C analysis of diploid human cells is difficult because of the lack of dense chromosome contact information and the presence of homologous chromosomes with very similar nucleotide sequences. Thus here, we propose a new algorithm to reconstruct the three-dimensional (3D) chromosomal architectures from the Hi-C dataset of single diploid human cells using allele-specific single-nucleotide variations (SNVs). We modified our recurrence plot-based algorithm, which is suitable for the estimation of the 3D chromosome structure from sparse Hi-C datasets, by newly incorporating a function of discriminating SNVs specific to each homologous chromosome. Here, we eventually regard a contact map as a recurrence plot. Importantly, the proposed method does not require any imputation for ambiguous segment information, but could efficiently reconstruct 3D chromosomal structures in single human diploid cells at a 1-Mb resolution. Datasets of segments without allele-specific SNVs, which were considered to be of little value, can also be used to validate the estimated chromosome structure. Introducing an additional mathematical measure called a refinement further improved the resolution to 40-kb or 100-kb. The reconstruction data supported the notion that human chromosomes form chromosomal territories and take fractal structures where the dimension for the underlying chromosome structure is a non-integer value.

HiC-TE: a computational pipeline for Hi-C data analysis to study the role of repeat family interactions in the genome 3D organization.

The role of repetitive DNA in the 3D organization of the interphase nucleus is a subject of intensive study. In studies of 3D nucleus organization, mutual contacts of various loci can be identified by Hi-C sequencing. Typical analyses use binning of read pairs by location to reduce noise. We use binning by repeat families instead to make similar conclusions about repeat regions. To achieve this, we combined Hi-C data, reference genome data and tools for repeat analysis into a Nextflow pipeline identifying and quantifying the contacts of specific repeat families. As an output, our pipeline produces heatmaps showing contact frequency and circular diagrams visualizing repeat contact localization. Using our pipeline with tomato data, we revealed the preferential homotypic interactions of ribosomal DNA, centromeric satellites and some LTR retrotransposon families and, as expected, little contact between organellar and nuclear DNA elements. While the pipeline can be applied to any eukaryotic genome, results in plants provide better coverage, since the built-in TE-greedy-nester software only detects tandems and LTR retrotransposons. Other repeats can be fed via GFF3 files. This pipeline represents a novel and reproducible way to analyze the role of repetitive elements in the 3D organization of genomes. https://gitlab.fi.muni.cz/lexa/hic-te/. Supplementary data are available at Bioinformatics online.

DeepLoop robustly maps chromatin interactions from sparse allele-resolved or single-cell Hi-C data at kilobase resolution.

Mapping chromatin loops from noisy Hi-C heatmaps remains a major challenge. Here we present DeepLoop, which performs rigorous bias correction followed by deep-learning-based signal enhancement for robust chromatin interaction mapping from low-depth Hi-C data. DeepLoop enables loop-resolution, single-cell Hi-C analysis. It also achieves a cross-platform convergence between different Hi-C protocols and micrococcal nuclease (micro-C). DeepLoop allowed us to map the genetic and epigenetic determinants of allele-specific chromatin interactions in the human genome. We nominate new loci with allele-specific interactions governed by imprinting or allelic DNA methylation. We also discovered that, in the inactivated X chromosome (Xi), local loops at the DXZ4 'megadomain' boundary escape X-inactivation but the FIRRE 'superloop' locus does not. Importantly, DeepLoop can pinpoint heterozygous single-nucleotide polymorphisms and large structure variants that cause allelic chromatin loops, many of which rewire enhancers with transcription consequences. Taken together, DeepLoop expands the use of Hi-C to provide loop-resolution insights into the genetics of the three-dimensional genome.

Genome-wide association and multi-omics studies identify MGMT as a novel risk gene for Alzheimer's disease among women.

Variants in the tau gene (MAPT) region are associated with breast cancer in women and Alzheimer's disease (AD) among persons lacking apolipoprotein E ε4 (ε4-). To identify novel genes associated with tau-related pathology, we conducted two genome-wide association studies (GWAS) for AD, one among 10,340 ε4- women in the Alzheimer's Disease Genetics Consortium (ADGC) and another in 31 members (22 women) of a consanguineous Hutterite kindred. We identified novel associations of AD with MGMT variants in the ADGC (rs12775171, odds ratio [OR] = 1.4, P=4.9×10-8) and Hutterite (rs12256016 and rs2803456, OR = 2.0, P=1.9×10-14) datasets. Multi-omics analyses showed that the most significant and largest number of associations among the single nucleotide polymorphisms (SNPs), DNA-methylated CpGs, MGMT expression, and AD-related neuropathological traits were observed among women. Furthermore, promoter capture Hi-C analyses revealed long-range interactions of the MGMT promoter with MGMT SNPs and CpG sites. These findings suggest that epigenetically regulated MGMT expression is involved in AD pathogenesis, especially in women.