Convergent evolution and multi-wave clonal invasion in H3 K27-altered diffuse midline gliomas treated with a PDGFR inhibitor.

The majority of diffuse midline gliomas, H3 K27-altered (DMG-H3 K27-a), are infiltrating pediatric brain tumors that arise in the pons with no effective treatment. To understand how clonal evolution contributes to the tumor's invasive spread, we performed exome sequencing and SNP array profiling on 49 multi-region autopsy samples from 11 patients with pontine DMG-H3 K27-a enrolled in a phase I clinical trial of PDGFR inhibitor crenolanib. For each patient, a phylogenetic tree was constructed by testing multiple possible clonal evolution models to select the one consistent with somatic mutations and copy number variations across all tumor regions. The tree was then used to deconvolute subclonal composition and prevalence at each tumor region to study convergent evolution and invasion patterns. Somatic variants in the PI3K pathway, a late event, are enriched in our cohort, affecting 70% of patients. Convergent evolution of PI3K at distinct phylogenetic branches was detected in 40% of the patients. 24 (~ 50%) of tumor regions were occupied by subclones of mixed lineages with varying molecular ages, indicating multiple waves of invasion across the pons and extrapontine. Subclones harboring a PDGFRA amplicon, including one that amplified a PDGRFAY849C mutant allele, were detected in four patients; their presence in extrapontine tumor and normal brain samples imply their involvement in extrapontine invasion. Our study expands the current knowledge on tumor invasion patterns in DMG-H3 K27-a, which may inform the design of future clinical trials.

Discovery of regulatory noncoding variants in individual cancer genomes by using cis-X.

We developed cis-X, a computational method for discovering regulatory noncoding variants in cancer by integrating whole-genome and transcriptome sequencing data from a single cancer sample. cis-X first finds aberrantly cis-activated genes that exhibit allele-specific expression accompanied by an elevated outlier expression. It then searches for causal noncoding variants that may introduce aberrant transcription factor binding motifs or enhancer hijacking by structural variations. Analysis of 13 T-lineage acute lymphoblastic leukemias identified a recurrent intronic variant predicted to cis-activate the TAL1 oncogene, a finding validated in vivo by chromatin immunoprecipitation sequencing of a patient-derived xenograft. Candidate oncogenes include the prolactin receptor PRLR activated by a focal deletion that removes a CTCF-insulated neighborhood boundary. cis-X may be applied to pediatric and adult solid tumors that are aneuploid and heterogeneous. In contrast to existing approaches, which require large sample cohorts, cis-X enables the discovery of regulatory noncoding variants in individual cancer genomes.

CICERO: a versatile method for detecting complex and diverse driver fusions using cancer RNA sequencing data.

To discover driver fusions beyond canonical exon-to-exon chimeric transcripts, we develop CICERO, a local assembly-based algorithm that integrates RNA-seq read support with extensive annotation for candidate ranking. CICERO outperforms commonly used methods, achieving a 95% detection rate for 184 independently validated driver fusions including internal tandem duplications and other non-canonical events in 170 pediatric cancer transcriptomes. Re-analysis of TCGA glioblastoma RNA-seq unveils previously unreported kinase fusions (KLHL7-BRAF) and a 13% prevalence of EGFR C-terminal truncation. Accessible via standard or cloud-based implementation, CICERO enhances driver fusion detection for research and precision oncology. The CICERO source code is available at https://github.com/stjude/Cicero.

Therapy-induced mutations drive the genomic landscape of relapsed acute lymphoblastic leukemia.

To study the mechanisms of relapse in acute lymphoblastic leukemia (ALL), we performed whole-genome sequencing of 103 diagnosis-relapse-germline trios and ultra-deep sequencing of 208 serial samples in 16 patients. Relapse-specific somatic alterations were enriched in 12 genes (NR3C1, NR3C2, TP53, NT5C2, FPGS, CREBBP, MSH2, MSH6, PMS2, WHSC1, PRPS1, and PRPS2) involved in drug response. Their prevalence was 17% in very early relapse (<9 months from diagnosis), 65% in early relapse (9-36 months), and 32% in late relapse (>36 months) groups. Convergent evolution, in which multiple subclones harbor mutations in the same drug resistance gene, was observed in 6 relapses and confirmed by single-cell sequencing in 1 case. Mathematical modeling and mutational signature analysis indicated that early relapse resistance acquisition was frequently a 2-step process in which a persistent clone survived initial therapy and later acquired bona fide resistance mutations during therapy. In contrast, very early relapses arose from preexisting resistant clone(s). Two novel relapse-specific mutational signatures, one of which was caused by thiopurine treatment based on in vitro drug exposure experiments, were identified in early and late relapses but were absent from 2540 pan-cancer diagnosis samples and 129 non-ALL relapses. The novel signatures were detected in 27% of relapsed ALLs and were responsible for 46% of acquired resistance mutations in NT5C2, PRPS1, NR3C1, and TP53. These results suggest that chemotherapy-induced drug resistance mutations facilitate a subset of pediatric ALL relapses.

CRISPR knockout screening identifies combinatorial drug targets in pancreatic cancer and models cellular drug response.

Predicting the response and identifying additional targets that will improve the efficacy of chemotherapy is a major goal in cancer research. Through large-scale in vivo and in vitro CRISPR knockout screens in pancreatic ductal adenocarcinoma cells, we identified genes whose genetic deletion or pharmacologic inhibition synergistically increase the cytotoxicity of MEK signaling inhibitors. Furthermore, we show that CRISPR viability scores combined with basal gene expression levels could model global cellular responses to the drug treatment. We develop drug response evaluation by in vivo CRISPR screening (DREBIC) method and validated its efficacy using large-scale experimental data from independent experiments. Comparative analyses demonstrate that DREBIC predicts drug response in cancer cells from a wide range of tissues with high accuracy and identifies therapeutic vulnerabilities of cancer-causing mutations to MEK inhibitors in various cancer types.

Genomic analysis of DNA repair genes and androgen signaling in prostate cancer.

BACKGROUND: The cellular effects of androgen are transduced through the androgen receptor, which controls the expression of genes that regulate biosynthetic processes, cell growth, and metabolism. Androgen signaling also impacts DNA damage signaling through mechanisms involving gene expression and transcription-associated DNA damaging events. Defining the contributions of androgen signaling to DNA repair is important for understanding androgen receptor function, and it also has translational implications. METHODS: We generated RNA-seq data from multiple prostate cancer lines and used bioinformatic analyses to characterize androgen-regulated gene expression. We compared the results from cell lines with gene expression data from prostate cancer xenografts, and patient samples, to query how androgen signaling and prostate cancer progression influences the expression of DNA repair genes. We performed whole genome sequencing to help characterize the status of the DNA repair machinery in widely used prostate cancer lines. Finally, we tested a DNA repair enzyme inhibitor for effects on androgen-dependent transcription. RESULTS: Our data indicates that androgen signaling regulates a subset of DNA repair genes that are largely specific to the respective model system and disease state. We identified deleterious mutations in the DNA repair genes RAD50 and CHEK2. We found that inhibition of the DNA repair enzyme MRE11 with the small molecule mirin inhibits androgen-dependent transcription and growth of prostate cancer cells. CONCLUSIONS: Our data supports the view that crosstalk between androgen signaling and DNA repair occurs at multiple levels, and that DNA repair enzymes in addition to PARPs, could be actionable targets in prostate cancer.

Live cell imaging of low- and non-repetitive chromosome loci using CRISPR-Cas9.

Imaging chromatin dynamics is crucial to understand genome organization and its role in transcriptional regulation. Recently, the RNA-guidable feature of CRISPR-Cas9 has been utilized for imaging of chromatin within live cells. However, these methods are mostly applicable to highly repetitive regions, whereas imaging regions with low or no repeats remains as a challenge. To address this challenge, we design single-guide RNAs (sgRNAs) integrated with up to 16 MS2 binding motifs to enable robust fluorescent signal amplification. These engineered sgRNAs enable multicolour labelling of low-repeat-containing regions using a single sgRNA and of non-repetitive regions with as few as four unique sgRNAs. We achieve tracking of native chromatin loci throughout the cell cycle and determine differential positioning of transcriptionally active and inactive regions in the nucleus. These results demonstrate the feasibility of our approach to monitor the position and dynamics of both repetitive and non-repetitive genomic regions in live cells.

The role of iron in neurodegeneration--Mössbauer spectroscopy, electron microscopy, enzyme-linked immunosorbent assay and neuroimaging studies.

The possible role of iron in neurodegeneration was studied by various techniques: electron microscopy, enzyme-linked immunosorbent assay, Mössbauer spectroscopy, atomic absorption, ultrasonography and magnetic resonance imaging. The measurements were made on human tissues extracted from liver and from brain structures involved in diseases of the human brain: substantia nigra (Parkinson's, PD), hippocampal cortex (Alzheimer's, AD) and globus pallidus (progressive supranuclear palsy, PSP). The sizes of the iron cores of ferritin, the main iron storage compound in tissues, were found to be smaller in brain than in liver. Brain ferritin has a higher proportion of H to L chains compared to liver. A significant decrease of the concentration of L chains in PD compared to control was found. No increase in the concentration of iron in PD versus control was detected; however, there was an increase of labile iron, which constitutes only 2‰ of brain iron. In AD an increase in the concentration of ferritin was noticed, without a significant increase in iron concentration. In PSP an increase of total iron was observed. Our findings suggest that the mechanisms leading to the death of nerve cells in these three diseases may be different, although all may be related to iron mediated oxidative stress.

Brain tissue echogenicity--implications for substantia nigra studies in parkinsonian patients.

The aim of the present study was to assess the origin of the substantia nigra hyperechogenicity in Parkinson disease patients. The cause of hyperechogenicity was tested on an animal model. Fresh porcine brains were injected consecutively with ferritin, apoferritin and water. Then, glioma samples were inserted into animal model. The echogenicity of the region of interest was assessed before and after experimental procedures. We observed the same echogenicity of porcine brain before and after injections of iron-loaded ferritin, apoferritin and water. Increased echogenicity of glioma samples compared to surrounding porcine brain tissue could be clearly seen. We postulate that the relative gliosis might be, at least partially, responsible for the increased echogenicity of the substantia nigra in Parkinson disease patients. Keeping in mind all limitations and inaccuracies of animal model used, it seems that hyperechogenicity of substantia nigra is caused rather by structural changes within the brain tissue than by increased iron concentration.