Machine and deep learning methods for predicting 3D genome organization.

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

Inducible deletion of Ezh2 in CD4+ T cells inhibits kidney T cell infiltration and prevents interstitial nephritis in MRL/lpr lupus-prone mice.

Systemic lupus erythematosus is a remitting relapsing autoimmune disease characterized by autoantibody production and multi-organ involvement. T cell epigenetic dysregulation plays an important role in the pathogenesis of lupus. We have previously demonstrated upregulation of the key epigenetic regulator EZH2 in CD4+ T cells isolated from lupus patients. To further investigate the role of EZH2 in the pathogenesis of lupus, we generated a tamoxifen-inducible CD4+ T cell Ezh2 conditional knockout mouse on the MRL/lpr lupus-prone background. We demonstrate that Ezh2 deletion abrogates lupus-like disease and prevents T cell differentiation. Single-cell analysis suggests impaired T cell function and activation of programed cell death pathways in EZH2-deficient mice. Ezh2 deletion in CD4+ T cells restricts TCR clonal repertoire and prevents kidney-infiltrating effector CD4+ T cell expansion and tubulointerstitial nephritis, which has been linked to end-stage renal disease in patients with lupus nephritis.

A novel multi-omics mendelian randomization method for gene set enrichment and its application to psychiatric disorders.

Genome-wide association studies (GWAS) of psychiatric disorders (PD) yield numerous loci with significant signals, but often do not implicate specific genes. Because GWAS risk loci are enriched in expression/protein/methylation quantitative loci (e/p/mQTL, hereafter xQTL), transcriptome/proteome/methylome-wide association studies (T/P/MWAS, hereafter XWAS) that integrate xQTL and GWAS information, can link GWAS signals to effects on specific genes. To further increase detection power, gene signals are aggregated within relevant gene sets (GS) by performing gene set enrichment (GSE) analyses. Often GSE methods test for enrichment of "signal" genes in curated GS while overlooking their linkage disequilibrium (LD) structure, allowing for the possibility of increased false positive rates. Moreover, no GSE tool uses xQTL information to perform mendelian randomization (MR) analysis. To make causal inference on association between PD and GS, we develop a novel MR GSE (MR-GSE) procedure. First, we generate a "synthetic" GWAS for each MSigDB GS by aggregating summary statistics for x-level (mRNA, protein or DNA methylation (DNAm) levels) from the largest xQTL studies available) of genes in a GS. Second, we use synthetic GS GWAS as exposure in a generalized summary-data-based-MR analysis of complex trait outcomes. We applied MR-GSE to GWAS of nine important PD. When applied to the underpowered opioid use disorder GWAS, none of the four analyses yielded any signals, which suggests a good control of false positive rates. For other PD, MR-GSE greatly increased the detection of GO terms signals (2,594) when compared to the commonly used (non-MR) GSE method (286). Some of the findings might be easier to adapt for treatment, e.g., our analyses suggest modest positive effects for supplementation with certain vitamins and/or omega-3 for schizophrenia, bipolar and major depression disorder patients. Similar to other MR methods, when applying MR-GSE researchers should be mindful of the confounding effects of horizontal pleiotropy on statistical inference.

Genome-wide analysis of hepatic DNA methylation reveals impact of epigenetic aging on xenobiotic metabolism and transport genes in an aged mouse model.

Hepatic xenobiotic metabolism and transport decline with age, while intact xenobiotic metabolism is associated with longevity. However, few studies have examined the genome-wide impact of epigenetic aging on these processes. We used reduced representation bisulfite sequencing (RRBS) to map DNA methylation changes in liver DNA from mice ages 4 and 24 months. We identified several thousand age-associated differentially methylated sites (a-DMS), many of which overlapped genes encoding Phase I and Phase II drug metabolizing enzymes, in addition to ABC and SLC classes of transporters. Notable genes harboring a-DMS were Cyp1a2, Cyp2d9, and Abcc2 that encode orthologs of the human drug metabolizing enzymes CYP1A2 and CYP2D6, and the multidrug resistance protein 2 (MRP2) transporter. Cyp2d9 hypermethylation with age was significantly associated with reduced gene expression, while Abcc2 expression was unchanged with age. Cyp1a2 lost methylation with age while, counterintuitively, its expression also reduced with age. We hypothesized that age-related dysregulation of the hepatic transcriptional machinery caused down-regulation of genes despite age-related hypomethylation. Bioinformatic analysis of hypomethylated a-DMS in our sample found them to be highly enriched for hepatic nuclear factor 4 alpha (HNF4α) binding sites. HNF4α promotes Cyp1a2 expression and is downregulated with age, which could explain the reduction in Cyp1a2 expression. Overall, our study supports the broad impact of epigenetic aging on xenobiotic metabolism and transport. Future work should evaluate the interplay between hepatic nuclear receptor function and epigenetic aging. These results may have implications for studies of longevity and healthy aging.

Targeting of SUMOylation leads to cBAF complex stabilization and disruption of the SS18::SSX transcriptome in Synovial Sarcoma.

Synovial Sarcoma (SS) is driven by the SS18::SSX fusion oncoprotein and is ultimately refractory to therapeutic approaches. SS18::SSX alters ATP-dependent chromatin remodeling BAF (mammalian SWI/SNF) complexes, leading to the degradation of canonical (cBAF) complex and amplified presence of an SS18::SSX-containing non-canonical BAF (ncBAF or GBAF) that drives an SS-specific transcription program and tumorigenesis. We demonstrate that SS18::SSX activates the SUMOylation program and SSs are sensitive to the small molecule SAE1/2 inhibitor, TAK-981. Mechanistically, TAK-981 de-SUMOylates the cBAF subunit SMARCE1, stabilizing and restoring cBAF on chromatin, shifting away from SS18::SSX-ncBAF-driven transcription, associated with DNA damage and cell death and resulting in tumor inhibition across both human and mouse SS tumor models. TAK-981 synergized with cytotoxic chemotherapy through increased DNA damage, leading to tumor regression. Targeting the SUMOylation pathway in SS restores cBAF complexes and blocks the SS18::SSX-ncBAF transcriptome, identifying a therapeutic vulnerability in SS, positioning the in-clinic TAK-981 to treat SS.

Targeting of SUMOylation leads to cBAF complex stabilization and disruption of the SS18::SSX transcriptome in Synovial Sarcoma.

Synovial Sarcoma (SS) is driven by the SS18::SSX fusion oncoprotein. and is ultimately refractory to therapeutic approaches. SS18::SSX alters ATP-dependent chromatin remodeling BAF (mammalian SWI/SNF) complexes, leading to the degradation of canonical (cBAF) complex and amplified presence of an SS18::SSX-containing non-canonical BAF (ncBAF or GBAF) that drives an SS-specific transcription program and tumorigenesis. We demonstrate that SS18::SSX activates the SUMOylation program and SSs are sensitive to the small molecule SAE1/2 inhibitor, TAK-981. Mechanistically, TAK-981 de-SUMOylates the cBAF subunit SMARCE1, stabilizing and restoring cBAF on chromatin, shifting away from SS18::SSX-ncBAF-driven transcription, associated with DNA damage and cell death and resulting in tumor inhibition across both human and mouse SS tumor models. TAK-981 synergized with cytotoxic chemotherapy through increased DNA damage, leading to tumor regression. Targeting the SUMOylation pathway in SS restores cBAF complexes and blocks the SS18::SSX-ncBAF transcriptome, identifying a therapeutic vulnerability in SS, positioning the in-clinic TAK-981 to treat SS.