Integrated drug profiling and CRISPR screening identify BCR::ABL1-independent vulnerabilities in chronic myeloid leukemia.

BCR::ABL1-independent pathways contribute to primary resistance to tyrosine kinase inhibitor (TKI) treatment in chronic myeloid leukemia (CML) and play a role in leukemic stem cell persistence. Here, we perform ex vivo drug screening of CML CD34+ leukemic stem/progenitor cells using 100 single drugs and TKI-drug combinations and identify sensitivities to Wee1, MDM2, and BCL2 inhibitors. These agents effectively inhibit primitive CD34+CD38- CML cells and demonstrate potent synergies when combined with TKIs. Flow-cytometry-based drug screening identifies mepacrine to induce differentiation of CD34+CD38- cells. We employ genome-wide CRISPR-Cas9 screening for six drugs, and mediator complex, apoptosis, and erythroid-lineage-related genes are identified as key resistance hits for TKIs, whereas the Wee1 inhibitor AZD1775 and mepacrine exhibit distinct resistance profiles. KCTD5, a consistent TKI-resistance-conferring gene, is found to mediate TKI-induced BCR::ABL1 ubiquitination. In summary, we delineate potential mechanisms for primary TKI resistance and non-BCR::ABL1-targeting drugs, offering insights for optimizing CML treatment.

Severe asthma trajectories in adults: findings from the NORDSTAR cohort.

BACKGROUND: There is limited evidence on the pathways leading to severe asthma and we are presently unable to effectively predict the progression of the disease. We aimed to describe the longitudinal trajectories leading to severe asthma and to describe clinical events preceding disease progression in a nationwide population of patients with severe asthma. METHODS: We conducted an observational study based on Swedish data from the NORdic Dataset for aSThmA Research (NORDSTAR) research collaboration platform. We identified adult patients with severe asthma in 2018 according to the European Respiratory Society/American Thoracic Society definition and used latent class analysis to identify trajectories of asthma severity over a 10-year retrospective period from 2018. RESULTS: Among 169 128 asthma patients, we identified 4543 severe asthma patients. We identified four trajectories of severe asthma that were labelled as: trajectory 1 "consistently severe asthma" (n=389 (8.6%)), trajectory 2 "gradual onset severe asthma" (n=942 (20.7%)), trajectory 3 "intermittent severe asthma" (n=1685 (37.1%)) and trajectory 4 "sudden onset severe asthma" (n=1527 (33.6%)). "Consistently severe asthma" had a higher daily inhaled corticosteroid dose and more prevalent osteoporosis compared with the other trajectories. Patients with "gradual onset severe asthma" and "sudden onset severe asthma" developed type 2-related comorbidities concomitantly with development of severe asthma. In the latter group, this primarily occurred within 1-3 years preceding onset of severe asthma. CONCLUSIONS: Four distinct trajectories of severe asthma were identified illustrating different patterns of progression of asthma severity. This may eventually enable the development of better preventive management strategies in severe asthma.

Prevalence and management of severe asthma in the Nordic countries: findings from the NORDSTAR cohort.

Background: Real-life evidence on prevalence and management of severe asthma is limited. Nationwide population registries across the Nordic countries provide unique opportunities to describe prevalence and management patterns of severe asthma at population level. In nationwide register data from Sweden, Norway and Finland, we examined the prevalence of severe asthma and the proportion of severe asthma patients being managed in specialist care. Methods: This is a cross-sectional study based on the Nordic Dataset for Asthma Research (NORDSTAR) research collaboration platform. We identified patients with severe asthma in adults (aged ≥18 years) and in children (aged 6-17 years) in 2018 according to the European Respiratory Society/American Thoracic Society definition. Patients managed in specialist care were those with an asthma-related specialist outpatient contact (only available in Sweden and Finland). Results: Overall, we identified 598 242 patients with current asthma in Sweden, Norway and Finland in 2018. Among those, the prevalence of severe asthma was 3.5%, 5.4% and 5.2% in adults and 0.4%, 1.0%, and 0.3% in children in Sweden, Norway and Finland, respectively. In Sweden and Finland, 37% and 40% of adult patients with severe asthma and two or more exacerbations, respectively, were managed in specialist care; in children the numbers were 56% and 41%, respectively. Conclusion: In three Nordic countries, population-based nationwide data demonstrated similar prevalence of severe asthma. In children, severe asthma was a rare condition. Notably, a large proportion of patients with severe asthma were not managed by a respiratory specialist, suggesting the need for increased recognition of severe asthma in primary care.

Implementing a Functional Precision Medicine Tumor Board for Acute Myeloid Leukemia.

We generated ex vivo drug-response and multiomics profiling data for a prospective series of 252 samples from 186 patients with acute myeloid leukemia (AML). A functional precision medicine tumor board (FPMTB) integrated clinical, molecular, and functional data for application in clinical treatment decisions. Actionable drugs were found for 97% of patients with AML, and the recommendations were clinically implemented in 37 relapsed or refractory patients. We report a 59% objective response rate for the individually tailored therapies, including 13 complete responses, as well as bridging five patients with AML to allogeneic hematopoietic stem cell transplantation. Data integration across all cases enabled the identification of drug response biomarkers, such as the association of IL15 overexpression with resistance to FLT3 inhibitors. Integration of molecular profiling and large-scale drug response data across many patients will enable continuous improvement of the FPMTB recommendations, providing a paradigm for individualized implementation of functional precision cancer medicine. SIGNIFICANCE: Oncogenomics data can guide clinical treatment decisions, but often such data are neither actionable nor predictive. Functional ex vivo drug testing contributes significant additional, clinically actionable therapeutic insights for individual patients with AML. Such data can be generated in four days, enabling rapid translation through FPMTB.See related commentary by Letai, p. 290.This article is highlighted in the In This Issue feature, p. 275.

Chimeric NUP98-NSD1 transcripts from the cryptic t(5;11)(q35.2;p15.4) in adult de novo acute myeloid leukemia.

The t(5;11)(q35;p15.4) is a clinically significant marker of poor prognosis in acute myeloid leukemia (AML), which is difficult to detect due to sub-telomeric localization of the breakpoints. To facilitate the detection of this rearrangement, we studied NUP98-NSD1 transcript variants in patients with the t(5;11) using paired-end RNA sequencing and standard molecular biology techniques. We discovered three NUP98-NSD1 transcripts with two fusion junctions (NUP98 exon 11-12/NSD1 exon 6), alternative 5' donor site in NUP98 exon 7, and NSD1 exon 7 skipping. Two of the transcripts were in-frame and occurred in all t(5;11) samples (N = 5). The exonic splicing events were present in all samples (N = 23) regardless of the NUP98-NSD1 suggesting that these novel splice events are unassociated with t(5;11). In conclusion, we provide evidence of two different NUP98-NSD1 fusion transcripts in adult AML, which result in functional proteins and represent suitable molecular entities for monitoring t(5;11) AML patients.

Intrinsic resistance to PIM kinase inhibition in AML through p38α-mediated feedback activation of mTOR signaling.

Although conventional therapies for acute myeloid leukemia (AML) and diffuse large B-cell lymphoma (DLBCL) are effective in inducing remission, many patients relapse upon treatment. Hence, there is an urgent need for novel therapies. PIM kinases are often overexpressed in AML and DLBCL and are therefore an attractive therapeutic target. However, in vitro experiments have demonstrated that intrinsic resistance to PIM inhibition is common. It is therefore likely that only a minority of patients will benefit from single agent PIM inhibitor treatment. In this study, we performed an shRNA-based genetic screen to identify kinases whose suppression is synergistic with PIM inhibition. Here, we report that suppression of p38α (MAPK14) is synthetic lethal with the PIM kinase inhibitor AZD1208. PIM inhibition elevates reactive oxygen species (ROS) levels, which subsequently activates p38α and downstream AKT/mTOR signaling. We found that p38α inhibitors sensitize hematological tumor cell lines to AZD1208 treatment in vitro and in vivo. These results were validated in ex vivo patient-derived AML cells. Our findings provide mechanistic and translational evidence supporting the rationale to test a combination of p38α and PIM inhibitors in clinical trials for AML and DLBCL.

Individualized systems medicine strategy to tailor treatments for patients with chemorefractory acute myeloid leukemia.

UNLABELLED: We present an individualized systems medicine (ISM) approach to optimize cancer drug therapies one patient at a time. ISM is based on (i) molecular profiling and ex vivo drug sensitivity and resistance testing (DSRT) of patients' cancer cells to 187 oncology drugs, (ii) clinical implementation of therapies predicted to be effective, and (iii) studying consecutive samples from the treated patients to understand the basis of resistance. Here, application of ISM to 28 samples from patients with acute myeloid leukemia (AML) uncovered five major taxonomic drug-response subtypes based on DSRT profiles, some with distinct genomic features (e.g., MLL gene fusions in subgroup IV and FLT3-ITD mutations in subgroup V). Therapy based on DSRT resulted in several clinical responses. After progression under DSRT-guided therapies, AML cells displayed significant clonal evolution and novel genomic changes potentially explaining resistance, whereas ex vivo DSRT data showed resistance to the clinically applied drugs and new vulnerabilities to previously ineffective drugs. SIGNIFICANCE: Here, we demonstrate an ISM strategy to optimize safe and effective personalized cancer therapies for individual patients as well as to understand and predict disease evolution and the next line of therapy. This approach could facilitate systematic drug repositioning of approved targeted drugs as well as help to prioritize and de-risk emerging drugs for clinical testing.

The HER2 amplicon includes several genes required for the growth and survival of HER2 positive breast cancer cells.

About 20% of breast cancers are characterized by amplification and overexpression of the HER2 oncogene. Although significant progress has been achieved for treating such patients with HER2 inhibitor trastuzumab, more than half of the patients respond poorly or become resistant to the treatment. Since the HER2 amplicon at 17q12 contains multiple genes, we have systematically explored the role of the HER2 co-amplified genes in breast cancer cell growth and their relation to trastuzumab resistance. We integrated aCGH data of the HER2 amplicon from 71 HER2 positive breast tumors and 10 cell lines with systematic functional RNA interference analysis of 23 core amplicon genes with several phenotypic endpoints in a panel of trastuzumab responding and non-responding HER2 positive breast cancer cells. Silencing of HER2 caused a greater growth arrest and apoptosis in the responding compared to the non-responding cell lines, indicating that the resistant cells are inherently less dependent on the HER2 pathway. Several other genes in the amplicon also showed a more pronounced effect when silenced; indicating that expression of HER2 co-amplified genes may be needed to sustain the growth of breast cancer cells. Importantly, co-silencing of STARD3, GRB7, PSMD3 and PERLD1 together with HER2 led to an additive inhibition of cell viability as well as induced apoptosis. These studies indicate that breast cancer cells may become addicted to the amplification of several genes that reside in the HER2 amplicon. The simultaneous targeting of these genes may increase the efficacy of the anti-HER2 therapies and possibly also counteract trastuzumab resistance. The observed additive effects seem to culminate to both apoptosis and cell proliferation pathways indicating that these pathways may be interesting targets for combinatorial treatment of HER2+ breast cancers.

Integrative genomic, transcriptomic, and RNAi analysis indicates a potential oncogenic role for FAM110B in castration-resistant prostate cancer.

BACKGROUND: Castration-resistant prostate cancer (CRPC) represents a therapeutic challenge for current medications. METHODS: In order to explore the molecular mechanisms involved in CRPC progression and to identify new therapeutic targets, we analyzed a unique sample set of 11 CRPCs and 7 advanced tumors by array-CGH and gene expression microarrays. The genome-wide DNA and RNA data were integrated to identify genes whose overexpression was driven by their amplification. To assess the functional role of these genes, their expression was analyzed in a transcriptional data set of 329 clinical prostate cancers and the corresponding gene products were silenced using RNA interference in prostate cancer cells. RESULTS: Six recurrent genetic targets were identified in the CRPCs; ATP1B1, AR, FAM110B, LAS1L, MYC, and YIPF6. In addition to AR and MYC, FAM110B emerged as a potential key gene involved in CRPC progression in a subset of the tumors. FAM110B was able to regulate AR signaling in prostate cancer cells and FAM110B itself was regulated by androgens. FAM110B siRNA inhibited the growth of prostate cancer cells in vitro, and this effect was substantially enhanced in androgen deficient conditions. Ectopic FAM110B expression in non-cancerous epithelial prostate cells induced aneuploidy and impaired antigen presentation. CONCLUSIONS: The DNA/RNA gene outlier detection combined with siRNA cell proliferation assay identified FAM110B as a potential growth promoting key gene for CRPC. FAM110B appears to have a key role in the androgen signaling and progression of CRPC impacting multiple cancer hallmarks and therefore highlighting a potential therapeutic target.

The miR-15a-miR-16-1 locus is homozygously deleted in a subset of prostate cancers.

MicroRNAs (miRNAs) are small, non-coding RNAs that negatively regulate the expression of protein coding genes. In this study, we screened highly informative prostate cancer cell lines and xenografts (n = 42) for miRNA gene copy number and expression changes. The expression profiling showed distinction between cell lines and xenografts as well as between androgen sensitive and independent models. Only a few copy number alterations that were associated with expression changes were identified. Most importantly, the miR-15a-miR-16-1 locus was found to be homozygously deleted in two samples leading to the abolishment of miR-15a, but not miR-16, expression. miR-16 is also expressed from another genomic locus. Mutation screening of the miR-15a-miR-16-1 gene in the model systems as well as clinical samples (n = 50) revealed no additional mutations. In conclusion, our data indicate that putative tumor suppressors, miR-15a and miR-16-1, are homozygously deleted in a subset of prostate cancers, further suggesting that these miRNAs could be important in the development of prostate cancer.

Identification of fusion genes in breast cancer by paired-end RNA-sequencing.

BACKGROUND: Until recently, chromosomal translocations and fusion genes have been an underappreciated class of mutations in solid tumors. Next-generation sequencing technologies provide an opportunity for systematic characterization of cancer cell transcriptomes, including the discovery of expressed fusion genes resulting from underlying genomic rearrangements. RESULTS: We applied paired-end RNA-seq to identify 24 novel and 3 previously known fusion genes in breast cancer cells. Supported by an improved bioinformatic approach, we had a 95% success rate of validating gene fusions initially detected by RNA-seq. Fusion partner genes were found to contribute promoters (5' UTR), coding sequences and 3' UTRs. Most fusion genes were associated with copy number transitions and were particularly common in high-level DNA amplifications. This suggests that fusion events may contribute to the selective advantage provided by DNA amplifications and deletions. Some of the fusion partner genes, such as GSDMB in the TATDN1-GSDMB fusion and IKZF3 in the VAPB-IKZF3 fusion, were only detected as a fusion transcript, indicating activation of a dormant gene by the fusion event. A number of fusion gene partners have either been previously observed in oncogenic gene fusions, mostly in leukemias, or otherwise reported to be oncogenic. RNA interference-mediated knock-down of the VAPB-IKZF3 fusion gene indicated that it may be necessary for cancer cell growth and survival. CONCLUSIONS: In summary, using RNA-sequencing and improved bioinformatic stratification, we have discovered a number of novel fusion genes in breast cancer, and identified VAPB-IKZF3 as a potential fusion gene with importance for the growth and survival of breast cancer cells.

Enhanced serine production by bone metastatic breast cancer cells stimulates osteoclastogenesis.

Since bone metastatic breast cancer is an incurable disease, causing significant morbidity and mortality, an understanding of the underlying molecular mechanisms would be highly valuable. Here, we describe in vitro and in vivo evidences for the importance of serine biosynthesis in the metastasis of breast cancer to bone. We first characterized the bone metastatic propensity of the MDA-MB-231(SA) cell line variant as compared to the parental MDA-MB-231 cells by radiographic and histological observations in the inoculated mice. Genome-wide gene expression profiling of this isogenic cell line pair revealed that all the three genes involved in the L: -serine biosynthesis pathway, phosphoglycerate dehydrogenase (PHGDH), phosphoserine aminotransferase 1 (PSAT1), and phosphoserine phosphatase (PSPH) were upregulated in the highly metastatic variant. This pathway is the primary endogenous source for L: -serine in mammalian tissues. Consistently, we observed that the proliferation of MDA-MB-231(SA) cells in serine-free conditions was dependent on PSAT1 expression. In addition, we observed that L: -serine is essential for the formation of bone resorbing human osteoclasts and may thus contribute to the vicious cycle of osteolytic bone metastasis. High expression of PHGDH and PSAT1 in primary breast cancer was significantly associated with decreased relapse-free and overall survival of patients and malignant phenotypic features of breast cancer. In conclusion, high expression of serine biosynthesis genes in metastatic breast cancer cells and the stimulating effect of L: -serine on osteoclastogenesis and cancer cell proliferation indicate a functionally critical role for serine biosynthesis in bone metastatic breast cancer and thereby an opportunity for targeted therapeutic interventions.

Integrative functional genomics analysis of sustained polyploidy phenotypes in breast cancer cells identifies an oncogenic profile for GINS2.

Aneuploidy is among the most obvious differences between normal and cancer cells. However, mechanisms contributing to development and maintenance of aneuploid cell growth are diverse and incompletely understood. Functional genomics analyses have shown that aneuploidy in cancer cells is correlated with diffuse gene expression signatures and aneuploidy can arise by a variety of mechanisms, including cytokinesis failures, DNA endoreplication, and possibly through polyploid intermediate states. To identify molecular processes contributing to development of aneuploidy, we used a cell spot microarray technique to identify genes inducing polyploidy and/or allowing maintenance of polyploid cell growth in breast cancer cells. Of 5760 human genes screened, 177 were found to induce severe DNA content alterations on prolonged transient silencing. Association with response to DNA damage stimulus and DNA repair was found to be the most enriched cellular processes among the candidate genes. Functional validation analysis of these genes highlighted GINS2 as the highest ranking candidate inducing polyploidy, accumulation of endogenous DNA damage, and impairing cell proliferation on inhibition. The cell growth inhibition and induction of polyploidy by suppression of GINS2 was verified in a panel of breast cancer cell lines. Bioinformatic analysis of published gene expression and DNA copy number studies of clinical breast tumors suggested GINS2 to be associated with the aggressive characteristics of a subgroup of breast cancers in vivo. In addition, nuclear GINS2 protein levels distinguished actively proliferating cancer cells suggesting potential use of GINS2 staining as a biomarker of cell proliferation as well as a potential therapeutic target.

Array-based gene expression, CGH and tissue data defines a 12q24 gain in neuroblastic tumors with prognostic implication.

BACKGROUND: Neuroblastoma has successfully served as a model system for the identification of neuroectoderm-derived oncogenes. However, in spite of various efforts, only a few clinically useful prognostic markers have been found. Here, we present a framework, which integrates DNA, RNA and tissue data to identify and prioritize genetic events that represent clinically relevant new therapeutic targets and prognostic biomarkers for neuroblastoma. METHODS: A single-gene resolution aCGH profiling was integrated with microarray-based gene expression profiling data to distinguish genetic copy number alterations that were strongly associated with transcriptional changes in two neuroblastoma cell lines. FISH analysis using a hotspot tumor tissue microarray of 37 paraffin-embedded neuroblastoma samples and in silico data mining for gene expression information obtained from previously published studies including up to 445 healthy nervous system samples and 123 neuroblastoma samples were used to evaluate the clinical significance and transcriptional consequences of the detected alterations and to identify subsequently activated gene(s). RESULTS: In addition to the anticipated high-level amplification and subsequent overexpression of MYCN, MEIS1, CDK4 and MDM2 oncogenes, the aCGH analysis revealed numerous other genetic alterations, including microamplifications at 2p and 12q24.11. Most interestingly, we identified and investigated the clinical relevance of a previously poorly characterized amplicon at 12q24.31. FISH analysis showed low-level gain of 12q24.31 in 14 of 33 (42%) neuroblastomas. Patients with the low-level gain had an intermediate prognosis in comparison to patients with MYCN amplification (poor prognosis) and to those with no MYCN amplification or 12q24.31 gain (good prognosis) (P = 0.001). Using the in silico data mining approach, we identified elevated expression of five genes located at the 12q24.31 amplicon in neuroblastoma (DIABLO, ZCCHC8, RSRC2, KNTC1 and MPHOSPH9). Among these, DIABLO showed the strongest activation suggesting a putative role in neuroblastoma progression. CONCLUSIONS: The presented systematic and rapid framework, which integrates aCGH, gene expression and tissue data to obtain novel targets and biomarkers for cancer, identified a low-level gain of the 12q24.31 as a potential new biomarker for neuroblastoma progression. Furthermore, results of in silico data mining suggest a new neuroblastoma target gene, DIABLO, within this region, whose functional and therapeutic role remains to be elucidated in follow-up studies.

Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line.

BACKGROUND: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. METHODS: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. RESULTS: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). CONCLUSION: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.

Systematic bioinformatic analysis of expression levels of 17,330 human genes across 9,783 samples from 175 types of healthy and pathological tissues.

Our knowledge on tissue- and disease-specific functions of human genes is rather limited and highly context-specific. Here, we have developed a method for the comparison of mRNA expression levels of most human genes across 9,783 Affymetrix gene expression array experiments representing 43 normal human tissue types, 68 cancer types, and 64 other diseases. This database of gene expression patterns in normal human tissues and pathological conditions covers 113 million datapoints and is available from the GeneSapiens website.

Integrin trafficking regulated by Rab21 is necessary for cytokinesis.

Adherent cells undergo remarkable changes in shape during cell division. However, the functional interplay between cell adhesion turnover and the mitotic machinery is poorly understood. The endo/exocytic trafficking of integrins is regulated by the small GTPase Rab21, which associates with several integrin alpha subunits. Here, we show that targeted trafficking of integrins to and from the cleavage furrow is required for successful cytokinesis, and that this is regulated by Rab21. Rab21 activity, integrin-Rab21 association, and integrin endocytosis are all necessary for normal cytokinesis, which becomes impaired when integrin-mediated adhesion at the cleavage furrow fails. We also describe a chromosomal deletion and loss of Rab21 gene expression in human cancer, which leads to the accumulation of multinucleate cells. Importantly, reintroduction of Rab21 rescued this phenotype. In conclusion, Rab21-regulated integrin trafficking is essential for normal cell division, and its defects may contribute to multinucleation and genomic instability, which are hallmarks of cancer.

A new look towards BAC-based array CGH through a comprehensive comparison with oligo-based array CGH.

BACKGROUND: Currently, two main technologies are used for screening of DNA copy number; the BAC (Bacterial Artificial Chromosome) and the recently developed oligonucleotide-based CGH (Chromosomal Comparative Genomic Hybridization) arrays which are capable of detecting small genomic regions with amplification or deletion. The correlation as well as the discriminative power of these platforms has never been compared statistically on a significant set of human patient samples. RESULTS: In this paper, we present an exhaustive comparison between the two CGH platforms, undertaken at two independent sites using the same batch of DNA from 19 advanced prostate cancers. The comparison was performed directly on the raw data and a significant correlation was found between the two platforms. The correlation was greatly improved when the data were averaged over large chromosomic regions using a segmentation algorithm. In addition, this analysis has enabled the development of a statistical model to discriminate BAC outliers that might indicate microevents. These microevents were validated by the oligo platform results. CONCLUSION: This article presents a genome-wide statistical validation of the oligo array platform on a large set of patient samples and demonstrates statistically its superiority over the BAC platform for the Identification of chromosomic events. Taking advantage of a large set of human samples treated by the two technologies, a statistical model has been developed to show that the BAC platform could also detect microevents.