indelPost: harmonizing ambiguities in simple and complex indel alignments.

SUMMARY: Small insertions and deletions (indels) in nucleotide sequence may be represented differently between mapping algorithms and variant callers, or in the flanking sequence context. Representational ambiguity is especially profound for complex indels, complicating comparisons between multiple mappings and call sets. Complex indels may additionally suffer from incomplete allele representation, potentially leading to critical misannotation of variant effect. We present indelPost, a Python library that harmonizes these ambiguities for simple and complex indels via realignment and read-based phasing. We demonstrate that indelPost enables accurate analysis of ambiguous data and can derive the correct complex indel alleles from the simple indel predictions provided by standard small variant detectors, with improved performance over a specialized tool for complex indel analysis. AVAILABILITY AND IMPLEMENTATION: indelPost is freely available at: https://github.com/stjude/indelPost. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

RNAIndel: discovering somatic coding indels from tumor RNA-Seq data.

MOTIVATION: Reliable identification of expressed somatic insertions/deletions (indels) is an unmet need due to artifacts generated in PCR-based RNA-Seq library preparation and the lack of normal RNA-Seq data, presenting analytical challenges for discovery of somatic indels in tumor transcriptome. RESULTS: We present RNAIndel, a tool for predicting somatic, germline and artifact indels from tumor RNA-Seq data. RNAIndel leverages features derived from indel sequence context and biological effect in a machine-learning framework. Except for tumor samples with microsatellite instability, RNAIndel robustly predicts 88-100% of somatic indels in five diverse test datasets of pediatric and adult cancers, even recovering subclonal (VAF range 0.01-0.15) driver indels missed by targeted deep-sequencing, outperforming the current best-practice for RNA-Seq variant calling which had 57% sensitivity but with 14 times more false positives. AVAILABILITY AND IMPLEMENTATION: RNAIndel is freely available at https://github.com/stjude/RNAIndel. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Enantioselective Denitrogenative Annulation of 1H-Tetrazoles with Styrenes Catalyzed by Rhodium.

Sulfonylation of 1H-tetrazoles with triflic anhydride in the presence of chiral rhodium(II) carboxylate dimers causes denitrogenation to generate α-azo rhodium(II) carbenoid species as new types of donor/acceptor carbenoids, which then readily react with styrenes to afford 3,5-diaryl-2-pyrazolines with a high degree of enantioselectivity.

18F-FBPA as a tumor-specific probe of L-type amino acid transporter 1 (LAT1): a comparison study with 18F-FDG and 11C-Methionine PET.

PURPOSE: The purpose of this study was to evaluate the usefulness of L-4-borono-2-18F-fluoro-phenylalanine (18F-FBPA) as a tumor-specific probe, in comparison to 18F-FDG and 11C-methionine (Met), focusing on its transport selectivity by L-type amino acid transporter 1 (LAT1), which is highly upregulated in cancers. METHODS: Cellular analyses of FBPA were performed to evaluate the transportablity and Km value. PET studies were performed in rat xenograft models of C6 glioma (n = 12) and in rat models of turpentine oil-induced subcutaneous inflammation (n = 9). The kinetic parameters and uptake values on static PET images were compared using the one-tissue compartment model (K1, k2) and maximum standardized uptake value (SUVmax). RESULTS: The cellular analyses showed that FBPA had a lower affinity to a normal cell-type transporter LAT2 and induced less efflux through LAT2 among FBPA, Met, and BPA, while the efflux through LAT1 induced by FBPA was similar among the three compounds. The Km value of 18F-FBPA for LAT1 (196.8 ± 11.4 µM) was dramatically lower than that for LAT2 (2813.8 ± 574.5 µM), suggesting the higher selectivity of 18F-FBPA for LAT1. K1 and k2 values were significantly smaller in 18F-FBPA PET (K1 = 0.04 ± 0.01 ml/ccm/min and k2 = 0.07 ± 0.01 /min) as compared to 11C-Met PET (0.22 ± 0.09 and 0.52 ± 0.10, respectively) in inflammatory lesions. Static PET analysis based on the SUVmax showed significantly higher accumulation of 18F-FDG in the tumor and inflammatory lesions (7.2 ± 2.1 and 4.6 ± 0.63, respectively) as compared to both 18F-FBPA (3.2 ± 0.40 and 1.9 ± 0.19) and 11C-Met (3.4 ± 0.43 and 1.6 ± 0.11). No significant difference was observed between 18F-FBPA and 11C-Met in the static PET images. CONCLUSION: This study shows the utility of 18F-FBPA as a tumor-specific probe of LAT1 with low accumulation in the inflammatory lesions.

Specific transport of 3-fluoro-l-α-methyl-tyrosine by LAT1 explains its specificity to malignant tumors in imaging.

3-(18)F-l-α-methyl-tyrosine ([18F]FAMT), a PET probe for tumor imaging, has advantages of high cancer-specificity and lower physiologic background. FAMT-PET has been proved useful in clinical studies for the prediction of prognosis, the assessment of therapy response and the differentiation of malignant tumors from inflammation and benign lesions. The tumor uptake of [18F]FAMT in PET is strongly correlated with the expression of L-type amino acid transporter 1 (LAT1), an isoform of system L upregulated in cancers. In this study, to assess the transporter-mediated mechanisms in FAMT uptake by tumors, we examined amino acid transporters for FAMT transport. We synthesized [14C]FAMT and measured its transport by human amino acid transporters expressed in Xenopus oocytes. The transport of FAMT was compared with that of l-methionine, a well-studied amino acid PET probe. The significance of LAT1 in FAMT uptake by tumor cells was confirmed by siRNA knockdown. Among amino acid transporters, [14C]FAMT was specifically transported by LAT1, whereas l-[14C]methionine was taken up by most of the transporters. Km of LAT1-mediated [14C]FAMT transport was 72.7 µM, similar to that for endogenous substrates. Knockdown of LAT1 resulted in the marked reduction of [14C]FAMT transport in HeLa S3 cells, confirming the contribution of LAT1 in FAMT uptake by tumor cells. FAMT is highly specific to cancer-type amino acid transporter LAT1, which explains the cancer-specific accumulation of [18F]FAMT in PET. This, vice versa, further supports the cancer-specific expression of LAT1. This study has established FAMT as a LAT1-specific molecular probe to monitor the expression of a potential tumor biomarker LAT1.

Transport of 3-fluoro-L-α-methyl-tyrosine (FAMT) by organic ion transporters explains renal background in [(18)F]FAMT positron emission tomography.

A PET tracer for tumor imaging, 3-(18)F-l-α-methyl-tyrosine ([(18)F]FAMT), has advantages of high cancer-specificity and low physiological background. In clinical studies, FAMT-PET has been proved useful for the detection of malignant tumors and their differentiation from inflammation and benign lesions. The tumor specific uptake of FAMT is due to its high-selectivity to cancer-type amino acid transporter LAT1 among amino acid transporters. In [(18)F]FAMT PET, kidney is the only organ that shows high physiological background. To reveal transporters involved in renal accumulation of FAMT, we have examined [(14)C]FAMT uptake on the organic ion transporters responsible for the uptake into tubular epithelial cells. We have found that OAT1, OAT10 and OCTN2 transport [(14)C]FAMT. The [(14)C]FAMT uptake was inhibited by probenecid, furosemide and ethacrynic acid, consistent with the properties of the transporters. The amino acid uptake inhibitor, 2-amino-2-norbornanecarboxylic acid (BCH), also inhibited the [(14)C]FAMT uptake, whereas OCTN2-mediated [(14)C]FAMT uptake was Na(+)-dependent. We propose that FAMT uptake by OAT1, OAT10 and OCTN2 into tubular epithelial cells could contribute to the renal accumulation of FAMT. The results from this study would provide clues to the treatments to reduce renal background and enhance tumor uptake as well as to designing PET tracers with less renal accumulation.

Boronophenylalanine, a boron delivery agent for boron neutron capture therapy, is transported by ATB0,+, LAT1 and LAT2.

The efficacy of boron neutron capture therapy relies on the selective delivery of boron carriers to malignant cells. p-Boronophenylalanine (BPA), a boron delivery agent, has been proposed to be localized to cells through transporter-mediated mechanisms. In this study, we screened aromatic amino acid transporters to identify BPA transporters. Human aromatic amino acid transporters were functionally expressed in Xenopus oocytes and examined for BPA uptake and kinetic parameters. The roles of the transporters in BPA uptake were characterized in cancer cell lines. For the quantitative assessment of BPA uptake, HPLC was used throughout the study. Among aromatic amino acid transporters, ATB(0,+), LAT1 and LAT2 were found to transport BPA with Km values of 137.4 ± 11.7, 20.3 ± 0.8 and 88.3 ± 5.6 µM, respectively. Uptake experiments in cancer cell lines revealed that the LAT1 protein amount was the major determinant of BPA uptake at 100 µM, whereas the contribution of ATB(0,+) became significant at 1000 µM, accounting for 20-25% of the total BPA uptake in MCF-7 breast cancer cells. ATB(0,+), LAT1 and LAT2 transport BPA at affinities comparable with their endogenous substrates, suggesting that they could mediate effective BPA uptake in vivo. The high and low affinities of LAT1 and ATB(0,+), respectively, differentiate their roles in BPA uptake. ATB(0,+), as well as LAT1, could contribute significantly to the tumor accumulation of BPA at clinical dose.

Clinical Response to a PARP Inhibitor and Chemotherapy in a Child with BARD1-Mutated Refractory Neuroblastoma: A Case Report.

Despite advances in the treatment of high-risk neuroblastoma, approximately half of these patients die from the disease. Targeted therapy based on synthetic lethality associated with homologous recombination deficiency (HRD) caused by germline mutations in homologous recombination repair genes has shown great efficacy in several adult solid tumors. Here we report the first successful treatment of a pediatric patient with refractory neuroblastoma and a germline pathogenic mutation in BARD1 using a PARP inhibitor, talazoparib, in combination with cytotoxic chemotherapy and radiation therapy. Allele-specific expression in RNA-seq indicates bi-allelic loss of BARD1 in tumor; however, the HRD score was below the threshold currently used for HRD classification in adult cancers. Our study demonstrates that the use of PARP inhibition in combination with DNA-damaging agents should be considered in children with BARD1-mutated neuroblastoma and cautions against the use of HRD score alone as a biomarker for this pediatric population.

Dynamics of Age- versus Therapy-Related Clonal Hematopoiesis in Long-term Survivors of Pediatric Cancer.

We present the first comprehensive investigation of clonal hematopoiesis (CH) in 2,860 long-term survivors of pediatric cancer with a median follow-up time of 23.5 years. Deep sequencing over 39 CH-related genes reveals mutations in 15% of the survivors, significantly higher than the 8.5% in 324 community controls. CH in survivors is associated with exposures to alkylating agents, radiation, and bleomycin. Therapy-related CH shows significant enrichment in STAT3, characterized as a CH gene specific to survivors of Hodgkin lymphoma, and TP53. Single-cell profiling of peripheral blood samples revealed STAT3 mutations predominantly present in T cells and contributed by SBS25, a mutational signature associated with procarbazine exposure. Serial sample tracking reveals that larger clone size is a predictor for future expansion of age-related CH clones, whereas therapy-related CH remains stable decades after treatment. These data depict the distinct dynamics of these CH subtypes and support the need for longitudinal monitoring to determine the potential contribution to late effects. SIGNIFICANCE: This first comprehensive CH analysis in long-term survivors of pediatric cancer presents the elevated prevalence and therapy exposures/diagnostic spectrum associated with CH. Due to the contrasting dynamics of clonal expansion for age-related versus therapy-related CH, longitudinal monitoring is recommended to ascertain the long-term effects of therapy-induced CH in pediatric cancer survivors. See related commentary by Collord and Behjati, p. 811. This article is highlighted in the In This Issue feature, p. 799.

Malignant Progression of an Ancestral Bone Marrow Clone Harboring a CIC-NUTM2A Fusion in Isolated Myeloid Sarcoma.

Myeloid sarcoma is a rare condition consisting of extramedullary myeloid blasts found in association with acute myeloid leukemia or, in the absence of bone marrow involvement. We identified an infant with isolated myeloid sarcoma whose bone marrow was negative for involvement by flow cytometry. Sequencing revealed the fusion oncogene CIC-NUTM2A and identified the sarcoma to be clonally evolved from the bone marrow, which carried the fusion despite the absence of pathology. Murine modeling confirmed the ability of the fusion to transform hematopoietic cells and identified receptor tyrosine kinase (RTK) signaling activation consistent with disruption of the CIC transcriptional repressor. These findings extend the definition of CIC-rearranged malignancies to include hematologic disease, provide insight into the mechanism of oncogenesis, and demonstrate the importance of molecular analysis and tracking of bone marrow involvement over the course of treatment in myeloid sarcoma, including patients that lack flow cytometric evidence of leukemia at diagnosis. IMPLICATIONS: This study illustrates molecular involvement of phenotypically normal bone marrow in myeloid sarcoma, which has significant implications in clinical care. Further, it extends the definition of CIC-rearrangements to include hematologic malignancies and shows evidence of RTK activation that may be exploited therapeutically in cancer(s) driven by these fusions.

Integrated Genomic Analysis Identifies UBTF Tandem Duplications as a Recurrent Lesion in Pediatric Acute Myeloid Leukemia.

The genetics of relapsed pediatric acute myeloid leukemia (AML) has yet to be comprehensively defined. Here, we present the spectrum of genomic alterations in 136 relapsed pediatric AMLs. We identified recurrent exon 13 tandem duplications (TD) in upstream binding transcription factor (UBTF) in 9% of relapsed AML cases. UBTF-TD AMLs commonly have normal karyotype or trisomy 8 with cooccurring WT1 mutations or FLT3-ITD but not other known oncogenic fusions. These UBTF-TD events are stable during disease progression and are present in the founding clone. In addition, we observed that UBTF-TD AMLs account for approximately 4% of all de novo pediatric AMLs, are less common in adults, and are associated with poor outcomes and MRD positivity. Expression of UBTF-TD in primary hematopoietic cells is sufficient to enhance serial clonogenic activity and to drive a similar transcriptional program to UBTF-TD AMLs. Collectively, these clinical, genomic, and functional data establish UBTF-TD as a new recurrent mutation in AML. SIGNIFICANCE: We defined the spectrum of mutations in relapsed pediatric AML and identified UBTF-TDs as a new recurrent genetic alteration. These duplications are more common in children and define a group of AMLs with intermediate-risk cytogenetic abnormalities, FLT3-ITD and WT1 alterations, and are associated with poor outcomes. See related commentary by Hasserjian and Nardi, p. 173. This article is highlighted in the In This Issue feature, p. 171.

The genomic landscape of pediatric acute lymphoblastic leukemia.

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Here, using whole-genome, exome and transcriptome sequencing of 2,754 childhood patients with ALL, we find that, despite a generally low mutation burden, ALL cases harbor a median of four putative somatic driver alterations per sample, with 376 putative driver genes identified varying in prevalence across ALL subtypes. Most samples harbor at least one rare gene alteration, including 70 putative cancer driver genes associated with ubiquitination, SUMOylation, noncoding transcripts and other functions. In hyperdiploid B-ALL, chromosomal gains are acquired early and synchronously before ultraviolet-induced mutation. By contrast, ultraviolet-induced mutations precede chromosomal gains in B-ALL cases with intrachromosomal amplification of chromosome 21. We also demonstrate the prognostic significance of genetic alterations within subtypes. Intriguingly, DUX4- and KMT2A-rearranged subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of the ALL genomic landscape and associated outcomes.

Profiling chromatin accessibility in pediatric acute lymphoblastic leukemia identifies subtype-specific chromatin landscapes and gene regulatory networks.

Acute lymphoblastic leukemia (ALL) is a hematopoietic malignancy comprised of molecular subtypes largely characterized by aneuploidy or recurring chromosomal rearrangements. Despite extensive information on the ALL transcriptome and methylome, there is limited understanding of the ALL chromatin landscape. We therefore mapped accessible chromatin in 24 primary ALL cell biospecimens comprising three common molecular subtypes (DUX4/ERG, ETV6-RUNX1 and hyperdiploid) from patients treated at St. Jude Children's Research Hospital. Our findings highlight extensive chromatin reprogramming in ALL, including the identification ALL subtype-specific chromatin landscapes that are additionally modulated by genetic variation. Chromatin accessibility differences between ALL and normal B-cells implicate the activation of B-cell repressed chromatin domains and detail the disruption of normal B-cell development in ALL. Among ALL subtypes, we uncovered roles for basic helix-loop-helix, homeodomain and activator protein 1 transcription factors in promoting subtype-specific chromatin accessibility and distinct gene regulatory networks. In addition to chromatin subtype-specificity, we further identified over 3500 DNA sequence variants that alter the ALL chromatin landscape and contribute to inter-individual variability in chromatin accessibility. Collectively, our data suggest that subtype-specific chromatin landscapes and gene regulatory networks impact ALL biology and contribute to transcriptomic differences among ALL subtypes.

Network-based systems pharmacology reveals heterogeneity in LCK and BCL2 signaling and therapeutic sensitivity of T-cell acute lymphoblastic leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive hematological malignancy, and novel therapeutics are much needed. Profiling patient leukemia' drug sensitivities ex vivo, we discovered that 44.4% of childhood and 16.7% of adult T-ALL cases exquisitely respond to dasatinib. Applying network-based systems pharmacology analyses to examine signal circuitry, we identified preTCR-LCK activation as the driver of dasatinib sensitivity, and T-ALL-specific LCK dependency was confirmed in genome-wide CRISPR-Cas9 screens. Dasatinib-sensitive T-ALLs exhibited high BCL-XL and low BCL2 activity and venetoclax resistance. Discordant sensitivity of T-ALL to dasatinib and venetoclax is strongly correlated with T-cell differentiation, particularly with the dynamic shift in LCK vs. BCL2 activation. Finally, single-cell analysis identified leukemia heterogeneity in LCK and BCL2 signaling and T-cell maturation stage, consistent with dasatinib response. In conclusion, our results indicate that developmental arrest in T-ALL drives differential activation of preTCR-LCK and BCL2 signaling in this leukemia, providing unique opportunities for targeted therapy.

Pan-neuroblastoma analysis reveals age- and signature-associated driver alterations.

Neuroblastoma is a pediatric malignancy with heterogeneous clinical outcomes. To better understand neuroblastoma pathogenesis, here we analyze whole-genome, whole-exome and/or transcriptome data from 702 neuroblastoma samples. Forty percent of samples harbor at least one recurrent driver gene alteration and most aberrations, including MYCN, ATRX, and TERT alterations, differ in frequency by age. MYCN alterations occur at median 2.3 years of age, TERT at 3.8 years, and ATRX at 5.6 years. COSMIC mutational signature 18, previously associated with reactive oxygen species, is the most common cause of driver point mutations in neuroblastoma, including most ALK and Ras-activating variants. Signature 18 appears early and is continuous throughout disease evolution. Signature 18 is enriched in neuroblastomas with MYCN amplification, 17q gain, and increased expression of mitochondrial ribosome and electron transport-associated genes. Recurrent FGFR1 variants in six patients, and ALK N-terminal structural alterations in five samples, identify additional patients potentially amenable to precision therapy.

Integrative genomic analyses reveal mechanisms of glucocorticoid resistance in acute lymphoblastic leukemia.

Identification of genomic and epigenomic determinants of drug resistance provides important insights for improving cancer treatment. Using agnostic genome-wide interrogation of mRNA and miRNA expression, DNA methylation, SNPs, CNAs and SNVs/Indels in primary human acute lymphoblastic leukemia cells, we identified 463 genomic features associated with glucocorticoid resistance. Gene-level aggregation identified 118 overlapping genes, 15 of which were confirmed by genome-wide CRISPR screen. Collectively, this identified 30 of 38 (79%) known glucocorticoid-resistance genes/miRNAs and all 38 known resistance pathways, while revealing 14 genes not previously associated with glucocorticoid-resistance. Single cell RNAseq and network-based transcriptomic modelling corroborated the top previously undiscovered gene, CELSR2. Manipulation of CELSR2 recapitulated glucocorticoid resistance in human leukemia cell lines and revealed a synergistic drug combination (prednisolone and venetoclax) that mitigated resistance in mouse xenograft models. These findings illustrate the power of an integrative genomic strategy for elucidating genes and pathways conferring drug resistance in cancer cells.