Genome-wide functional perturbation of human microsatellite repeats using engineered zinc finger transcription factors.

Repeat elements can be dysregulated at a genome-wide scale in human diseases. For example, in Ewing sarcoma, hundreds of inert GGAA repeats can be converted into active enhancers when bound by EWS-FLI1. Here we show that fusions between EWS and GGAA-repeat-targeted engineered zinc finger arrays (ZFAs) can function at least as efficiently as EWS-FLI1 for converting hundreds of GGAA repeats into active enhancers in a Ewing sarcoma precursor cell model. Furthermore, a fusion of a KRAB domain to a ZFA can silence GGAA microsatellite enhancers genome wide in Ewing sarcoma cells, thereby reducing expression of EWS-FLI1-activated genes. Remarkably, this KRAB-ZFA fusion showed selective toxicity against Ewing sarcoma cells compared with non-Ewing cancer cells, consistent with its Ewing sarcoma-specific impact on the transcriptome. These findings demonstrate the value of ZFAs for functional annotation of repeats and illustrate how aberrant microsatellite activities might be regulated for potential therapeutic applications.

Single-cell analysis and functional characterization uncover the stem cell hierarchies and developmental origins of rhabdomyosarcoma.

Rhabdomyosarcoma (RMS) is a common childhood cancer that shares features with developing skeletal muscle. Yet, the conservation of cellular hierarchy with human muscle development and the identification of molecularly defined tumor-propagating cells has not been reported. Using single-cell RNA-sequencing, DNA-barcode cell fate mapping and functional stem cell assays, we uncovered shared tumor cell hierarchies in RMS and human muscle development. We also identified common developmental stages at which tumor cells become arrested. Fusion-negative RMS cells resemble early myogenic cells found in embryonic and fetal development, while fusion-positive RMS cells express a highly specific gene program found in muscle cells transiting from embryonic to fetal development at 7-7.75 weeks of age. Fusion-positive RMS cells also have neural pathway-enriched states, suggesting less-rigid adherence to muscle-lineage hierarchies. Finally, we identified a molecularly defined tumor-propagating subpopulation in fusion-negative RMS that shares remarkable similarity to bi-potent, muscle mesenchyme progenitors that can make both muscle and osteogenic cells.

Therapeutic targeting of LCK tyrosine kinase and mTOR signaling in T-cell acute lymphoblastic leukemia.

Relapse and refractory T-cell acute lymphoblastic leukemia (T-ALL) has a poor prognosis, and new combination therapies are sorely needed. Here, we used an ex vivo high-throughput screening platform to identify drug combinations that kill zebrafish T-ALL and then validated top drug combinations for preclinical efficacy in human disease. This work uncovered potent drug synergies between AKT/mTORC1 (mammalian target of rapamycin complex 1) inhibitors and the general tyrosine kinase inhibitor dasatinib. Importantly, these same drug combinations effectively killed a subset of relapse and dexamethasone-resistant zebrafish T-ALL. Clinical trials are currently underway using the combination of mTORC1 inhibitor temsirolimus and dasatinib in other pediatric cancer indications, leading us to prioritize this therapy for preclinical testing. This combination effectively curbed T-ALL growth in human cell lines and primary human T-ALL and was well tolerated and effective in suppressing leukemia growth in patient-derived xenografts (PDX) grown in mice. Mechanistically, dasatinib inhibited phosphorylation and activation of the lymphocyte-specific protein tyrosine kinase (LCK) to blunt the T-cell receptor (TCR) signaling pathway, and when complexed with mTORC1 inhibition, induced potent T-ALL cell killing through reducing MCL-1 protein expression. In total, our work uncovered unexpected roles for the LCK kinase and its regulation of downstream TCR signaling in suppressing apoptosis and driving continued leukemia growth. Analysis of a wide array of primary human T-ALLs and PDXs grown in mice suggest that combination of temsirolimus and dasatinib treatment will be efficacious for a large fraction of human T-ALLs.

The chromatin landscape of primary synovial sarcoma organoids is linked to specific epigenetic mechanisms and dependencies.

Synovial sarcoma (SyS) is an aggressive mesenchymal malignancy invariably associated with the chromosomal translocation t(X:18; p11:q11), which results in the in-frame fusion of the BAF complex gene SS18 to one of three SSX genes. Fusion of SS18 to SSX generates an aberrant transcriptional regulator, which, in permissive cells, drives tumor development by initiating major chromatin remodeling events that disrupt the balance between BAF-mediated gene activation and polycomb-dependent repression. Here, we developed SyS organoids and performed genome-wide epigenomic profiling of these models and mesenchymal precursors to define SyS-specific chromatin remodeling mechanisms and dependencies. We show that SS18-SSX induces broad BAF domains at its binding sites, which oppose polycomb repressor complex (PRC) 2 activity, while facilitating recruitment of a non-canonical (nc)PRC1 variant. Along with the uncoupling of polycomb complexes, we observed H3K27me3 eviction, H2AK119ub deposition and the establishment of de novo active regulatory elements that drive SyS identity. These alterations are completely reversible upon SS18-SSX depletion and are associated with vulnerability to USP7 loss, a core member of ncPRC1.1. Using the power of primary tumor organoids, our work helps define the mechanisms of epigenetic dysregulation on which SyS cells are dependent.

hPCL3S promotes proliferation and migration of androgen-independent prostate cancer cells.

Polycomb repressive complex 2 (PRC2) allows the deposition of H3K27me3. PRC2 facultative subunits modulate its activity and recruitment such as hPCL3/PHF19, a human ortholog of Drosophila Polycomb-like protein (PCL). These proteins contain a TUDOR domain binding H3K36me3, two PHD domains and a "Winged-helix" domain involved in GC-rich DNA binding. The human PCL3 locus encodes the full-length hPCL3L protein and a shorter isoform, hPCL3S containing the TUDOR and PHD1 domains only. In this study, we demonstrated by RT-qPCR analyses of 25 prostate tumors that hPCL3S is frequently up-regulated. In addition, hPCL3S is overexpressed in the androgen-independent DU145 and PC3 cells, but not in the androgen-dependent LNCaP cells. hPCL3S knockdown decreased the proliferation and migration of DU145 and PC3 whereas its forced expression into LNCaP increased these properties. A mutant hPCL3S unable to bind H3K36me3 (TUDOR-W50A) increased proliferation and migration of LNCaP similarly to wt hPCL3S whereas inactivation of its PHD1 domain decreased proliferation. These effects partially relied on the up-regulation of genes known to be important for the proliferation and/or migration of prostate cancer cells such as S100A16, PlexinA2, and Spondin1. Collectively, our results suggest hPCL3S as a new potential therapeutic target in castration resistant prostate cancers.

Two faces of bivalent domain regulate VEGFA responsiveness and angiogenesis.

The bivalent domain (BD) at promoter region is an unique epigenetic feature poised for activation or repression during cell differentiation in embryonic stem cell. However, the function of BDs in already differentiated cells remains exclusive. By profiling the epigenetic landscape of endothelial cells during VEGFA (vascular endothelial growth factor A) stimulation, we discovered that BDs are widespread in endothelial cells and preferentially marked genes responsive to VEGFA. The BDs responsive to VEGFA have more permissive chromatin environment comparing to other BDs. The initial activation of bivalent genes depends on RNAPII pausing release induced by EZH1 rather than removal of H3K27me3. The later suppression of bivalent gene expression depended on KDM5A recruitment by its interaction with PRC2. Importantly, EZH1 promoted both in vitro and in vivo angiogenesis by upregulating EGR3, whereas KDM5A dampened angiogenesis. Collectively, this study demonstrates a novel dual function of BDs in endothelial cells to control VEGF responsiveness and angiogenesis.

High levels of AAV vector integration into CRISPR-induced DNA breaks.

Adeno-associated virus (AAV) vectors have shown promising results in preclinical models, but the genomic consequences of transduction with AAV vectors encoding CRISPR-Cas nucleases is still being examined. In this study, we observe high levels of AAV integration (up to 47%) into Cas9-induced double-strand breaks (DSBs) in therapeutically relevant genes in cultured murine neurons, mouse brain, muscle and cochlea. Genome-wide AAV mapping in mouse brain shows no overall increase of AAV integration except at the CRISPR/Cas9 target site. To allow detailed characterization of integration events we engineer a miniature AAV encoding a 465 bp lambda bacteriophage DNA (AAV-λ465), enabling sequencing of the entire integrated vector genome. The integration profile of AAV-465λ in cultured cells display both full-length and fragmented AAV genomes at Cas9 on-target sites. Our data indicate that AAV integration should be recognized as a common outcome for applications that utilize AAV for genome editing.

Allele-specific gene editing prevents deafness in a model of dominant progressive hearing loss.

Since most dominant human mutations are single nucleotide substitutions1,2, we explored gene editing strategies to disrupt dominant mutations efficiently and selectively without affecting wild-type alleles. However, single nucleotide discrimination can be difficult to achieve3 because commonly used endonucleases, such as Streptococcus pyogenes Cas9 (SpCas9), can tolerate up to seven mismatches between guide RNA (gRNA) and target DNA. Furthermore, the protospacer-adjacent motif (PAM) in some Cas9 enzymes can tolerate mismatches with the target DNA3,4. To circumvent these limitations, we screened 14 Cas9/gRNA combinations for specific and efficient disruption of a nucleotide substitution that causes the dominant progressive hearing loss, DFNA36. As a model for DFNA36, we used Beethoven mice5, which harbor a point mutation in Tmc1, a gene required for hearing that encodes a pore-forming subunit of mechanosensory transduction channels in inner-ear hair cells6. We identified a PAM variant of Staphylococcus aureus Cas9 (SaCas9-KKH) that selectively and efficiently disrupted the mutant allele, but not the wild-type Tmc1/TMC1 allele, in Beethoven mice and in a DFNA36 human cell line. Adeno-associated virus (AAV)-mediated SaCas9-KKH delivery prevented deafness in Beethoven mice up to one year post injection. Analysis of current ClinVar entries revealed that ~21% of dominant human mutations could be targeted using a similar approach.

Transcriptome-wide off-target RNA editing induced by CRISPR-guided DNA base editors.

CRISPR-Cas base-editor technology enables targeted nucleotide alterations, and is being increasingly used for research and potential therapeutic applications1,2. The most widely used cytosine base editors (CBEs) induce deamination of DNA cytosines using the rat APOBEC1 enzyme, which is targeted by a linked Cas protein-guide RNA complex3,4. Previous studies of the specificity of CBEs have identified off-target DNA edits in mammalian cells5,6. Here we show that a CBE with rat APOBEC1 can cause extensive transcriptome-wide deamination of RNA cytosines in human cells, inducing tens of thousands of C-to-U edits with frequencies ranging from 0.07% to 100% in 38-58% of expressed genes. CBE-induced RNA edits occur in both protein-coding and non-protein-coding sequences and generate missense, nonsense, splice site, and 5' and 3' untranslated region mutations. We engineered two CBE variants bearing mutations in rat APOBEC1 that substantially decreased the number of RNA edits (by more than 390-fold and more than 3,800-fold) in human cells. These variants also showed more precise on-target DNA editing than the wild-type CBE and, for most guide RNAs tested, no substantial reduction in editing efficiency. Finally, we show that an adenine base editor7 can also induce transcriptome-wide RNA edits. These results have implications for the use of base editors in both research and clinical settings, illustrate the feasibility of engineering improved variants with reduced RNA editing activities, and suggest the need to more fully define and characterize the RNA off-target effects of deaminase enzymes in base editor platforms.

A dynamic and integrated epigenetic program at distal regions orchestrates transcriptional responses to VEGFA.

Cell behaviors are dictated by epigenetic and transcriptional programs. Little is known about how extracellular stimuli modulate these programs to reshape gene expression and control cell behavioral responses. Here, we interrogated the epigenetic and transcriptional response of endothelial cells to VEGFA treatment and found rapid chromatin changes that mediate broad transcriptomic alterations. VEGFA-responsive genes were associated with active promoters, but changes in promoter histone marks were not tightly linked to gene expression changes. VEGFA altered transcription factor occupancy and the distal epigenetic landscape, which profoundly contributed to VEGFA-dependent changes in gene expression. Integration of gene expression, dynamic enhancer, and transcription factor occupancy changes induced by VEGFA yielded a VEGFA-regulated transcriptional regulatory network, which revealed that the small MAF transcription factors are master regulators of the VEGFA transcriptional program and angiogenesis. Collectively these results revealed that extracellular stimuli rapidly reconfigure the chromatin landscape to coordinately regulate biological responses.

tp53 deficiency causes a wide tumor spectrum and increases embryonal rhabdomyosarcoma metastasis in zebrafish.

The TP53 tumor-suppressor gene is mutated in >50% of human tumors and Li-Fraumeni patients with germ line inactivation are predisposed to developing cancer. Here, we generated tp53 deleted zebrafish that spontaneously develop malignant peripheral nerve-sheath tumors, angiosarcomas, germ cell tumors, and an aggressive Natural Killer cell-like leukemia for which no animal model has been developed. Because the tp53 deletion was generated in syngeneic zebrafish, engraftment of fluorescent-labeled tumors could be dynamically visualized over time. Importantly, engrafted tumors shared gene expression signatures with predicted cells of origin in human tissue. Finally, we showed that tp53del/del enhanced invasion and metastasis in kRASG12D-induced embryonal rhabdomyosarcoma (ERMS), but did not alter the overall frequency of cancer stem cells, suggesting novel pro-metastatic roles for TP53 loss-of-function in human muscle tumors. In summary, we have developed a Li-Fraumeni zebrafish model that is amenable to large-scale transplantation and direct visualization of tumor growth in live animals.

Vangl2/RhoA Signaling Pathway Regulates Stem Cell Self-Renewal Programs and Growth in Rhabdomyosarcoma.

Tumor growth and relapse are driven by tumor propagating cells (TPCs). However, mechanisms regulating TPC fate choices, maintenance, and self-renewal are not fully understood. Here, we show that Van Gogh-like 2 (Vangl2), a core regulator of the non-canonical Wnt/planar cell polarity (Wnt/PCP) pathway, affects TPC self-renewal in rhabdomyosarcoma (RMS)-a pediatric cancer of muscle. VANGL2 is expressed in a majority of human RMS and within early mononuclear progenitor cells. VANGL2 depletion inhibited cell proliferation, reduced TPC numbers, and induced differentiation of human RMS in vitro and in mouse xenografts. Using a zebrafish model of embryonal rhabdomyosarcoma (ERMS), we determined that Vangl2 expression enriches for TPCs and promotes their self-renewal. Expression of constitutively active and dominant-negative isoforms of RHOA revealed that it acts downstream of VANGL2 to regulate proliferation and maintenance of TPCs in human RMS. Our studies offer insights into pathways that control TPCs and identify new potential therapeutic targets.

Transcriptional diversity during lineage commitment of human blood progenitors.

Blood cells derive from hematopoietic stem cells through stepwise fating events. To characterize gene expression programs driving lineage choice, we sequenced RNA from eight primary human hematopoietic progenitor populations representing the major myeloid commitment stages and the main lymphoid stage. We identified extensive cell type-specific expression changes: 6711 genes and 10,724 transcripts, enriched in non-protein-coding elements at early stages of differentiation. In addition, we found 7881 novel splice junctions and 2301 differentially used alternative splicing events, enriched in genes involved in regulatory processes. We demonstrated experimentally cell-specific isoform usage, identifying nuclear factor I/B (NFIB) as a regulator of megakaryocyte maturation-the platelet precursor. Our data highlight the complexity of fating events in closely related progenitor populations, the understanding of which is essential for the advancement of transplantation and regenerative medicine.