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.

Topoisomerase II contributes to DNA secondary structure-mediated double-stranded breaks.

DNA double-stranded breaks (DSBs) trigger human genome instability, therefore identifying what factors contribute to DSB induction is critical for our understanding of human disease etiology. Using an unbiased, genome-wide approach, we found that genomic regions with the ability to form highly stable DNA secondary structures are enriched for endogenous DSBs in human cells. Human genomic regions predicted to form non-B-form DNA induced gross chromosomal rearrangements in yeast and displayed high indel frequency in human genomes. The extent of instability in both analyses is in concordance with the structure forming ability of these regions. We also observed an enrichment of DNA secondary structure-prone sites overlapping transcription start sites (TSSs) and CCCTC-binding factor (CTCF) binding sites, and uncovered an increase in DSBs at highly stable DNA secondary structure regions, in response to etoposide, an inhibitor of topoisomerase II (TOP2) re-ligation activity. Importantly, we found that TOP2 deficiency in both yeast and human leads to a significant reduction in DSBs at structure-prone loci, and that sites of TOP2 cleavage have a greater ability to form highly stable DNA secondary structures. This study reveals a direct role for TOP2 in generating secondary structure-mediated DNA fragility, advancing our understanding of mechanisms underlying human genome instability.

Pausing sites of RNA polymerase II on actively transcribed genes are enriched in DNA double-stranded breaks.

DNA double-stranded breaks (DSBs) are strongly associated with active transcription, and promoter-proximal pausing of RNA polymerase II (Pol II) is a critical step in transcriptional regulation. Mapping the distribution of DSBs along actively expressed genes and identifying the location of DSBs relative to pausing sites can provide mechanistic insights into transcriptional regulation. Using genome-wide DNA break mapping/sequencing techniques at single-nucleotide resolution in human cells, we found that DSBs are preferentially located around transcription start sites of highly transcribed and paused genes and that Pol II promoter-proximal pausing sites are enriched in DSBs. We observed that DSB frequency at pausing sites increases as the strength of pausing increases, regardless of whether the pausing sites are near or far from annotated transcription start sites. Inhibition of topoisomerase I and II by camptothecin and etoposide treatment, respectively, increased DSBs at the pausing sites as the concentrations of drugs increased, demonstrating the involvement of topoisomerases in DSB generation at the pausing sites. DNA breaks generated by topoisomerases are short-lived because of the religation activity of these enzymes, which these drugs inhibit; therefore, the observation of increased DSBs with increasing drug doses at pausing sites indicated active recruitment of topoisomerases to these sites. Furthermore, the enrichment and locations of DSBs at pausing sites were shared among different cell types, suggesting that Pol II promoter-proximal pausing is a common regulatory mechanism. Our findings support a model in which topoisomerases participate in Pol II promoter-proximal pausing and indicated that DSBs at pausing sites contribute to transcriptional activation.

Complex DNA structures trigger copy number variation across the Plasmodium falciparum genome.

Antimalarial resistance is a major obstacle in the eradication of the human malaria parasite, Plasmodium falciparum. Genome amplifications, a type of DNA copy number variation (CNV), facilitate overexpression of drug targets and contribute to parasite survival. Long monomeric A/T tracks are found at the breakpoints of many Plasmodium resistance-conferring CNVs. We hypothesize that other proximal sequence features, such as DNA hairpins, act with A/T tracks to trigger CNV formation. By adapting a sequence analysis pipeline to investigate previously reported CNVs, we identified breakpoints in 35 parasite clones with near single base-pair resolution. Using parental genome sequence, we predicted the formation of stable hairpins within close proximity to all future breakpoint locations. Especially stable hairpins were predicted to form near five shared breakpoints, establishing that the initiating event could have occurred at these sites. Further in-depth analyses defined characteristics of these 'trigger sites' across the genome and detected signatures of error-prone repair pathways at the breakpoints. We propose that these two genomic signals form the initial lesion (hairpins) and facilitate microhomology-mediated repair (A/T tracks) that lead to CNV formation across this highly repetitive genome. Targeting these repair pathways in P. falciparum may be used to block adaptation to antimalarial drugs.

CNCC: an analysis tool to determine genome-wide DNA break end structure at single-nucleotide resolution.

BACKGROUND: DNA double-stranded breaks (DSBs) are potentially deleterious events in a cell. The end structures (blunt, 3'- and 5'-overhangs) at DSB sites contribute to the fate of their repair and provide critical information concerning the consequences of the damage. Therefore, there has been a recent eruption of DNA break mapping and sequencing methods that aim to map at single-nucleotide resolution where breaks are generated genome-wide. These methods provide high resolution data for the location of DSBs, which can encode the type of end-structure present at these breaks. However, genome-wide analysis of the resulting end structures has not been investigated following these sequencing methods. RESULTS: To address this analysis gap, we develop the use of a coverage-normalized cross correlation analysis (CNCC) to process the high-precision genome-wide break mapping data, and determine genome-wide break end structure distributions at single-nucleotide resolution. We take advantage of the single-nucleotide position and the knowledge of strandness from every mapped break to analyze the relative shifts between positive and negative strand encoded break nucleotides. By applying CNCC we can identify the most abundant end structures captured by a break mapping technique, and further can make comparisons between different samples and treatments. We validate our analysis with restriction enzyme digestions of genomic DNA and establish the sensitivity of the analysis using end structures that only exist as a minor fraction of total breaks. Finally, we demonstrate the versatility of our analysis by applying CNCC to the breaks resulting after treatment with etoposide and study the variety of resulting end structures. CONCLUSION: For the first time, on a genome-wide scale, our analysis revealed the increase in the 5' to 3' end resection following etoposide treatment, and the global progression of the resection. Furthermore, our method distinguished the change in the pattern of DSB end structure with increasing doses of the drug. The ability of this method to determine DNA break end structures without a priori knowledge of break sequences or genomic position should have broad applications in understanding genome instability.

CRISPR-STOP: gene silencing through base-editing-induced nonsense mutations.

CRISPR-Cas9-induced DNA damage may have deleterious effects at high-copy-number genomic regions. Here, we use CRISPR base editors to knock out genes by changing single nucleotides to create stop codons. We show that the CRISPR-STOP method is an efficient and less deleterious alternative to wild-type Cas9 for gene-knockout studies. Early stop codons can be introduced in ∼17,000 human genes. CRISPR-STOP-mediated targeted screening demonstrates comparable efficiency to WT Cas9, which indicates the suitability of our approach for genome-wide functional screenings.

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.

Basal ganglia echogenicity in tauopathies.

Accumulating data confirm the usefulness of transcranial sonography (TCS) in the diagnosis of Parkinson's disease. The relevance of basal ganglia abnormalities depicted by TCS in atypical parkinsonian syndromes still needs further assessment. In the present study, 20 patients with progressive supranuclear palsy (PSP) and 13 patients with corticobasal syndrome (CBS) were studied with the use of transcranial sonography. Echogenicity of the substantia nigra (SN) and lenticular nucleus (LN) were assessed. 0/20 patients with PSP and 8/12 (66.6 %) patients with CBS were characterized with SN hyperechogenicity. LN hyperechogenicity was observed in 9/20 patients diagnosed with PSP and 0/11 of CBS patients. The combination of SN isoechogenicity and LN hyperechogenicity reached 100 % sensitivity and positive predictive value for the diagnosis of PSP. The results of this study point out that CBS has to be taken into consideration when SN hyperechogenicity is depicted in a patient with parkinsonian syndrome. Normal echogenicity of the SN coexisting with LN hyperechogenicity practically excludes CBS.

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.

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.

SequencErr: measuring and suppressing sequencer errors in next-generation sequencing data.

BACKGROUND: There is currently no method to precisely measure the errors that occur in the sequencing instrument/sequencer, which is critical for next-generation sequencing applications aimed at discovering the genetic makeup of heterogeneous cellular populations. RESULTS: We propose a novel computational method, SequencErr, to address this challenge by measuring the base correspondence between overlapping regions in forward and reverse reads. An analysis of 3777 public datasets from 75 research institutions in 18 countries revealed the sequencer error rate to be ~ 10 per million (pm) and 1.4% of sequencers and 2.7% of flow cells have error rates > 100 pm. At the flow cell level, error rates are elevated in the bottom surfaces and > 90% of HiSeq and NovaSeq flow cells have at least one outlier error-prone tile. By sequencing a common DNA library on different sequencers, we demonstrate that sequencers with high error rates have reduced overall sequencing accuracy, and removal of outlier error-prone tiles improves sequencing accuracy. We demonstrate that SequencErr can reveal novel insights relative to the popular quality control method FastQC and achieve a 10-fold lower error rate than popular error correction methods including Lighter and Musket. CONCLUSIONS: Our study reveals novel insights into the nature of DNA sequencing errors incurred on DNA sequencers. Our method can be used to assess, calibrate, and monitor sequencer accuracy, and to computationally suppress sequencer errors in existing datasets.

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.

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.

Alternative DNA secondary structure formation affects RNA polymerase II promoter-proximal pausing in human.

BACKGROUND: Alternative DNA secondary structures can arise from single-stranded DNA when duplex DNA is unwound during DNA processes such as transcription, resulting in the regulation or perturbation of these processes. We identify sites of high propensity to form stable DNA secondary structure across the human genome using Mfold and ViennaRNA programs with parameters for analyzing DNA. RESULTS: The promoter-proximal regions of genes with paused transcription are significantly and energetically more favorable to form DNA secondary structure than non-paused genes or genes without RNA polymerase II (Pol II) binding. Using Pol II ChIP-seq, GRO-seq, NET-seq, and mNET-seq data, we arrive at a robust set of criteria for Pol II pausing, independent of annotation, and find that a highly stable secondary structure is likely to form about 10-50 nucleotides upstream of a Pol II pausing site. Structure probing data confirm the existence of DNA secondary structures enriched at the promoter-proximal regions of paused genes in human cells. Using an in vitro transcription assay, we demonstrate that Pol II pausing at HSPA1B, a human heat shock gene, is affected by manipulating DNA secondary structure upstream of the pausing site. CONCLUSIONS: Our results indicate alternative DNA secondary structure formation as a mechanism for how GC-rich sequences regulate RNA Pol II promoter-proximal pausing genome-wide.

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.

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.

Crystal structure of equine serum albumin in complex with cetirizine reveals a novel drug binding site.

Serum albumin (SA) is the main transporter of drugs in mammalian blood plasma. Here, we report the first crystal structure of equine serum albumin (ESA) in complex with antihistamine drug cetirizine at a resolution of 2.1Å. Cetirizine is bound in two sites--a novel drug binding site (CBS1) and the fatty acid binding site 6 (CBS2). Both sites differ from those that have been proposed in multiple reports based on equilibrium dialysis and fluorescence studies for mammalian albumins as cetirizine binding sites. We show that the residues forming the binding pockets in ESA are highly conserved in human serum albumin (HSA), and suggest that binding of cetirizine to HSA will be similar. In support of that hypothesis, we show that the dissociation constants for cetirizine binding to CBS2 in ESA and HSA are identical using tryptophan fluorescence quenching. Presence of lysine and arginine residues that have been previously reported to undergo nonenzymatic glycosylation in CBS1 and CBS2 suggests that cetirizine transport in patients with diabetes could be altered. A review of all available SA structures from the PDB shows that in addition to the novel drug binding site we present here (CBS1), there are two pockets on SA capable of binding drugs that do not overlap with fatty acid binding sites and have not been discussed in published reviews.

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.