We sought to identify drugs that could counteract cytarabine resistance in acute myeloid leukemia (AML) by generating eight resistant variants from MOLM-13 and SHI-1 AML cell lines by long-term drug treatment. These cells were compared with 66 ex vivo chemorefractory samples from cytarabine-treated AML patients. The models and patient cells were subjected to genomic and transcriptomic profiling and high-throughput testing with 250 emerging and clinical oncology compounds. Genomic profiling uncovered deletion of the deoxycytidine kinase (DCK) gene in both MOLM-13- and SHI-1-derived cytarabine-resistant variants and in an AML patient sample. Cytarabine-resistant SHI-1 variants and a subset of chemorefractory AML patient samples showed increased sensitivity to glucocorticoids that are often used in treatment of lymphoid leukemia but not AML. Paired samples taken from AML patients before treatment and at relapse also showed acquisition of glucocorticoid sensitivity. Enhanced glucocorticoid sensitivity was only seen in AML patient samples that were negative for the FLT3 mutation (P=0.0006). Our study shows that development of cytarabine resistance is associated with increased sensitivity to glucocorticoids in a subset of AML, suggesting a new therapeutic strategy that should be explored in a clinical trial of chemorefractory AML patients carrying wild-type FLT3.

Renal cell cancer (RCC) has become a prototype example of the extensive intratumor heterogeneity and clonal evolution of human cancers. However, there is little direct evidence on how the genetic heterogeneity impacts on drug response profiles of the cancer cells. Our goal was to determine how genomic clonal evolution impacts drug responses. Finding from our study could help to define the challenge that clonal evolution poses on cancer therapy. We established multiple patient-derived cells (PDCs) from different tumor regions of four RCC patients, verified their clonal relationship to each other and to the uncultured tumor tissue by genome sequencing. Furthermore, comprehensive drug-sensitivity testing with 460 oncological drugs was performed on all PDC clones. The PDCs retained many cancer-specific copy number alterations and mutations in driver genes such as VHL, PBRM1, PIK3C2A, KMD5C and TSC2 genes. The drug testing highlighted vulnerability in the PDCs toward approved RCC drugs, such as the mTOR-inhibitor temsirolimus, but also novel sensitivities were uncovered. The individual PDC clones from different tumor regions in a patient showed distinct drug-response profiles, suggesting that genomic heterogeneity contributes to the variability in drug responses. Studies of multiple PDCs from a patient with cancer are informative for elucidating cancer heterogeneity and for the determination on how the genomic evolution is manifested in cancer drug responsiveness. This approach could facilitate tailoring of drugs and drug combinations to individual patients.

In our individualized systems medicine program, personalized treatment options are identified and administered to chemorefractory acute myeloid leukemia (AML) patients based on exome sequencing and ex vivo drug sensitivity and resistance testing data. Here, we analyzed how clonal heterogeneity affects the responses of 13 AML patients to chemotherapy or targeted treatments using ultra-deep (average 68 000 × coverage) amplicon resequencing. Using amplicon resequencing, we identified 16 variants from 4 patients (frequency 0.54-2%) that were not detected previously by exome sequencing. A correlation-based method was developed to detect mutation-specific responses in serial samples across multiple time points. Significant subclone-specific responses were observed for both chemotherapy and targeted therapy. We detected subclonal responses in patients where clinical European LeukemiaNet (ELN) criteria showed no response. Subclonal responses also helped to identify putative mechanisms underlying drug sensitivities, such as sensitivity to azacitidine in DNMT3A mutated cell clones and resistance to cytarabine in a subclone with loss of NF1 gene. In summary, ultra-deep amplicon resequencing method enables sensitive quantification of subclonal variants and their responses to therapies. This approach provides new opportunities for designing combinatorial therapies blocking multiple subclones as well as for real-time assessment of such treatments.