Faculty Bibliography

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

Epilepsy is a heterogenous group of disorders defined by recurrent seizure activity due to abnormal synchronized activity of neurons. A growing number of epilepsy cases are believed to be caused by genetic factors and copy number variants (CNV) contribute to up to 5% of epilepsy cases. However, CNVs in epilepsy are usually large deletions or duplications involving multiple neurodevelopmental genes. In patients who underwent seizure focus resection for treatment-resistant epilepsy, whole genome DNA methylation profiling identified 3 main clusters of which one showed strong association with receptor tyrosine kinase (RTK) genes. We identified focal copy number gains involving epidermal growth factor receptor (EGFR) and PDGFRA loci. The dysplastic neurons of cases with amplifications showed marked overexpression of EGFR and PDGFRA, while glial and endothelial cells were negative. Targeted sequencing of regulatory regions and DNA methylation analysis revealed that only enhancer regions of EGFR and gene promoter of PDGFRA were amplified, while coding regions did not show copy number abnormalities or somatic mutations. Somatic focal copy number gains of noncoding regulatory represent a previously unrecognized genetic driver in epilepsy and a mechanism of abnormal activation of RTK genes. Upregulated RTKs provide a potential avenue for therapy in seizure disorders.

DNA methylation is thought to regulate accessibility of chromatin and binding of regulatory elements; however, it is difficult to determine if chromatin accessibility or transcription factor (TF) binding overlap with methylated or unmethylated DNA if the assays are performed separately. In order to examine accessibility or TF binding simultaneously with methylation on the same DNA molecule, we developed EpiMethylTag which combines ATAC-Seq or ChIP-Seq (M-ATAC or M-ChIP) with bisulfite conversion. Our approach provides a fast, low-input, low sequencing depth method to determine whether DNAme and accessibility/TF binding are mutually exclusive or can coexist in certain locations.

Inhibiting membrane association of RAS has long been considered a rational approach to anticancer therapy, which led to the development of farnesyltransferase inhibitors (FTIs). However, FTIs proved ineffective against KRAS-driven tumors. To reveal alternative therapeutic strategies, we carried out a genome-wide CRISPR-Cas9 screen designed to identify genes required for KRAS4B membrane association. We identified five enzymes in the prenylation pathway and SAFB, a nuclear protein with both DNA and RNA binding domains. Silencing SAFB led to marked mislocalization of all RAS isoforms as well as RAP1A but not RAB7A, a pattern that phenocopied silencing FNTA, the prenyltransferase α subunit shared by farnesyltransferase and geranylgeranyltransferase type I. We found that SAFB promoted RAS membrane association by controlling FNTA expression. SAFB knockdown decreased GTP loading of RAS, abrogated alternative prenylation, and sensitized RAS-mutant cells to growth inhibition by FTI. Our work establishes the prenylation pathway as paramount in KRAS membrane association, reveals a regulator of prenyltransferase expression, and suggests that reduction in FNTA expression may enhance the efficacy of FTIs.