Faculty Bibliography

Investigating how chromatin organization determines cell-type-specific gene expression remains challenging. Experimental methods for measuring three-dimensional chromatin organization, such as Hi-C, are costly and have technical limitations, restricting their broad application particularly in high-throughput genetic perturbations. We present C.Origami, a multimodal deep neural network that performs de novo prediction of cell-type-specific chromatin organization using DNA sequence and two cell-type-specific genomic features-CTCF binding and chromatin accessibility. C.Origami enables in silico experiments to examine the impact of genetic changes on chromatin interactions. We further developed an in silico genetic screening approach to assess how individual DNA elements may contribute to chromatin organization and to identify putative cell-type-specific trans-acting regulators that collectively determine chromatin architecture. Applying this approach to leukemia cells and normal T cells, we demonstrate that cell-type-specific in silico genetic screening, enabled by C.Origami, can be used to systematically discover novel chromatin regulation circuits in both normal and disease-related biological systems.

Resistance to targeted therapies is an important clinical problem in HER2-positive (HER2+) breast cancer. "Drug-tolerant persisters" (DTPs), a sub-population of cancer cells that survive via reversible, non-genetic mechanisms, are implicated in resistance to tyrosine kinase inhibitors (TKIs) in other malignancies, but DTPs following HER2 TKI exposure have not been well characterized. We found that HER2 TKIs evoke DTPs with a luminal-like or a mesenchymal-like transcriptome. Lentiviral barcoding/single cell RNA-sequencing reveal that HER2+ breast cancer cells cycle stochastically through a "pre-DTP" state, characterized by a G0-like expression signature and enriched for diapause and/or senescence genes. Trajectory analysis/cell sorting show that pre-DTPs preferentially yield DTPs upon HER2 TKI exposure. Cells with similar transcriptomes are present in HER2+ breast tumors and are associated with poor TKI response. Finally, biochemical experiments indicate that luminal-like DTPs survive via estrogen receptor-dependent induction of SGK3, leading to rewiring of the PI3K/AKT/mTORC1 pathway to enable AKT-independent mTORC1 activation.

In this interview, Professor Jane Skok speaks with Storm Johnson, commissioning editor for Epigenomics, on her work to date in the field of chromosome architecture and regulatory elements. Jane Skok's lab uses sophisticated microscopic techniques to visualize recombination in individual cells, tracing the dynamic changes in chromosome architecture and nuclear location at different stages of this complex process. This line of research unites two lifelong passions: science and art. After completing her PhD in immunology and genetics at the Imperial Cancer Research Fund in Lincoln's Inn Fields, Dr Skok took 12 years off and pursued training in art while caring for her young children. She then returned to science, joining David Gray's lab at Imperial College London as a postdoctoral fellow to study B cell biology and acquired expertise in Mandy Fisher's lab to understand how nuclear organization of the antigen receptor genes regulate V(D)J recombination and allelic exclusion. Dr Skok continued to pursue these questions in her own lab at University College London and elucidated the roles of Pax5, locus contraction and nuclear subcompartmentalization in maintaining allelic exclusion. In 2006, Dr Skok was recruited to New York University School of Medicine, where her lab has revealed the activities of several signaling factors in guiding B cell development and they made the surprising discovery that the RAG proteins and the DNA damage response factor ATM help ensure allelic exclusion at the immunoglobulin gene loci. More recently, those at the Skok lab have turned their attention to understanding how localized and long-range chromatin contacts impact gene regulation in health and disease settings.

CCCTC-binding factor (CTCF) is critical to three-dimensional genome organization. Upon differentiation, CTCF insulates active and repressed genes within Hox gene clusters. We conducted a genome-wide CRISPR knockout (KO) screen to identify genes required for CTCF-boundary activity at the HoxA cluster, complemented by biochemical approaches. Among the candidates, we identified Myc-associated zinc-finger protein (MAZ) as a cofactor in CTCF insulation. MAZ colocalizes with CTCF at chromatin borders and, similar to CTCF, interacts with the cohesin subunit RAD21. MAZ KO disrupts gene expression and local contacts within topologically associating domains. Similar to CTCF motif deletions, MAZ motif deletions lead to derepression of posterior Hox genes immediately after CTCF boundaries upon differentiation, giving rise to homeotic transformations in mouse. Thus, MAZ is a factor contributing to appropriate insulation, gene expression and genomic architecture during development.

Background: Relapse from acute lymphoblastic leukemia (ALL) is one of the most common causes of pediatric cancer-related death. Early relapse of ALL is associated with recurrent mutations of histone methyltransferase NSD2 (nuclear receptor binding SET domain protein 2), which is specific for H3K36me2, suggesting a link to therapy resistance or other mechanisms underlying relapse such as central neural system (CNS) infiltration. NSD2 p.E1099K affects gene expression through disturbing the balance of H3K36me2/H3K27me3. Using CRISPR/Cas9-edited isogenic ALL cell lines and PDX cells, we found that NSD2 p.E1099K drives oncogenic programming, CNS infiltration and glucocorticoid (GC) resistance. However, the molecular mechanisms underlying the relapse of this subtype of ALL are still under investigation.
Aim(s): To elucidate the epigenetic landscape underlying the mechanism of the relapse of pediatric ALL with NSD2 mutation.
Method(s): We performed in vivo experiments to observe tumor burden, leukemia cell infiltration and survival of the NOD/SCID mice injected with a NSD2 p.E1099K mutation knock-out SEM cell line and knock-in CEM cell line. We determined transcriptome (RNA-Seq), chromatin accessibility (ATAC-Seq) in isogenic RCH-ACV, SEM, RPMI-8402 and CEM cell lines, transcription factor binding and histone modification (ChIP-Seq) and 3D organization (Hi-C) in RCH-ACV cells. Finally, we integrated analysis of RNA-Seq, ATAC-Seq, ChIP-Seq and Hi-C to comprehensively disclose the epigenetic landscape in ALL with NSD2 mutation.
Result(s): NOD/SCID mice xenografted with NSD2 mutant cells developed high tumor burden and infiltration to spleen, liver and brain while the mice injected with WT cells accumulated significant less tumor cells and had extended survival. RNA-Seq analysis showed that reversion of NSD2 mutation to WT caused more genes downregulated while insertion of NSD2 mutation to WT cells led to more genes upregulated (Mutant vs WT: RCH-ACV 838 vs 494, SEM 1567 vs 1158, RPMI-8402 1922 vs 1745, CEM 1809 vs 1031). 50 upregulated genes and 3 downregulated genes were in common in B-ALL and T-ALL with NSD2 mutation. Most of the upregulated genes correlated with neural development and adhesion which might contribute to CNS infiltration (e.g., NCAM1 and NEO1). A few genes were associated with GC resistance such as decreased NR3C1 and increased NR3C2. Accordingly, ATAC-Seq showed that NSD2 mutant cells had more open chromatin peaks than those of WT (RCH-ACV 4853 vs 3212, SEM 10052 vs 7595, RPMI-8402 20392 vs 12133, CEM 10155 vs 6437). ChIP-Seq revealed general large gains of H3K36me2 in intergenic regions in NSD2 mutant cells. Importantly, genes upregulated with NSD2 mutation (e.g., NCAM1 and NEO1) lost H3K27me3 at promoters but gained H3K36me2 at promoters and whole gene bodies, accompanied with increased H3K27ac at enhancers. Conversely, a small subset of genes gained H3K27me3 and lost H3K36me2 in their promoters. Concentrated H3K36me2 in gene bodies diffused and broadened was less prominent and H3K27me3 accumulation became dominant. This for example was associated with repression of NR3C1 to drive GC resistance of NSD2 mutant cells. Genes upregulated in NSD2 mutant cells were enriched for binding sites for lymphoid transcriptional activators such as EBF1 and IRF2. The promoters of the downregulated genes had motifs for transcription factors poorly expressed in lymphoid cells and were enriched for binding sites for the BCL6 transcriptional repressor. Hi-C analysis revealed 430 topologically associated domains (TADs) with increased loop interactions while 136 TADs with decreased interactions were in NSD2 mutant cells compared to WT cells. Overall, 491 regions switched from compartment A to B and 444 regions switched from B to A in NSD2 mutant cells compared to WT cells. Compartment switching from inactive B to active A correlated with upregulated gene expression levels while the reverse was true for switching from A to B. Increased intra-TAD interactions were linked to upregulated genes while decreased intra-TAD interactions were linked to downregulated genes.
Conclusion(s): The NSD2 mutation led to increased tumor burden, CNS infiltration and glucocorticoid resistance due to dysregulation of epigenetic patterns and 3D chromatin architecture, indicating mechanisms underlying the relapse of pediatric ALL and potential therapeutic targets in ALL with NSD2 mutation. Disclosures: Licht: Epizyme: Research Funding.
Copyright