Genome-wide chemical mutagenesis screens allow unbiased saturation of the cancer genome and identification of drug resistance mutations

Base editing screens define the genetic landscape of cancer drug resistance mechanisms

Drug resistance is a principal limitation to the long-term efficacy of cancer therapies. Cancer genome sequencing can retrospectively delineate the genetic basis of drug resistance, but this requires large numbers of post-treatment samples to nominate causal variants. Here we prospectively identify genetic mechanisms of resistance to ten oncology drugs from CRISPR base editing mutagenesis screens in four cancer cell lines using a guide RNA library predicted to install 32,476 variants in 11 cancer genes. We identify four functional classes of protein variants modulating drug sensitivity and use single-cell transcriptomics to reveal how these variants operate through distinct mechanisms, including eliciting a drug-addicted cell state. We identify variants that can be targeted with alternative inhibitors to overcome resistance and functionally validate an epidermal growth factor receptor (EGFR) variant that sensitizes lung cancer cells to EGFR inhibitors. Our variant-to-function map has implications for patient stratification, therapy combinations and drug scheduling in cancer treatment.

In vitro evolution and whole genome analysis to study chemotherapy drug resistance in haploid human cells

In vitro evolution and whole genome analysis has proven to be a powerful method for studying the mechanism of action of small molecules in many haploid microbes but has generally not been applied to human cell lines in part because their diploid state complicates the identification of variants that confer drug resistance. To determine if haploid human cells could be used in MOA studies, we evolved resistance to five different anticancer drugs (doxorubicin, gemcitabine, etoposide, topotecan, and paclitaxel) using a near-haploid cell line (HAP1) and then analyzed the genomes of the drug resistant clones, developing a bioinformatic pipeline that involved filtering for high frequency alleles predicted to change protein sequence, or alleles which appeared in the same gene for multiple independent selections with the same compound. Applying the filter to sequences from 28 drug resistant clones identified a set of 21 genes which was strongly enriched for known resistance genes or known drug targets (TOP1, TOP2A, DCK, WDR33, SLCO3A1). In addition, some lines carried structural variants that encompassed additional known resistance genes (ABCB1, WWOX and RRM1). Gene expression knockdown and knockout experiments of 10 validation targets showed a high degree of specificity and accuracy in our calls and demonstrates that the same drug resistance mechanisms found in diverse clinical samples can be evolved, discovered and studied in an isogenic background.

Comparative Analysis of Proteomic of Curcumin Reversing Multidrug Resistance in HCT-8/VCR Cells

To further explore the mechanisms of curcumin reversing multidrug resistance (MDR) in HCT8/VCR cells. Here, we employed comparative analysis of proteomic of essential proteins of human colon carcinoma HCT8/VCR cells with or without treatment of curcumin by separating and quantifying the essential protein posttranslational modification through radical-free two-dimensional polyacrylamide gel electrophoresis with strong reductant. The reverse impact of curcumin on multidrug resistance of HCT8/VCR and HCT8/VCR cells was evaluated using MTT assay. After adding curcumin 25 μM for 72 h, by 2-DE and mass spectrometry, twenty proteins were certified with changed expression levels. Three protein sites were upregulated and seventeen protein sites were downregulated in curcumin-treated HCT-8/VCR. Verification analyses were conducted using RT-PCR and Western blotting for downregulated proteins including GSTP1 and PRDX6. The proteins might have a direct or indirect contact with multidrug resistance. The finding of the research would provide novel sights for systematically comprehending the mechanisms of the reversal impacts of curcumin on MDR in HCT8/VCR cells and contribute to the recognition and application of new markers in clinical practice.

Resistance to obinutuzumab-induced antibody-dependent cellular cytotoxicity caused by abnormal Fas signaling is overcome by combination therapies

Background

Obinutuzumab, a Type II anti-CD20 antibody, is used to treat follicular lymphoma. A major mode of action of obinutuzumab is antibody-dependent cellular cytotoxicity (ADCC). Knowledge of the mechanisms of resistance to obinutuzumab is important for the development of next-line strategies to follow obinutuzumab-containing therapy, including obinutuzumab retreatment. Unfortunately, the mechanisms by which tumor cells acquire resistance to ADCC are still poorly understood. To address this, we examined the mechanisms of resistance to obinutuzumab-induced ADCC and the combination efficacy of obinutuzumab and clinically available agents in the established resistant cells.

Methods and results

We established cells resistant to obinutuzumab-induced ADCC using the non-Hodgkin lymphoma cell line RL and examined their mechanisms of resistance and the combination efficacy of obinutuzumab and clinically available agents. Comprehensive analysis by RNA sequencing of resistance mechanisms revealed that abnormal Fas signaling decreased sensitivity to ADCC in resistant clones. Combination treatment with prednisolone, a component of CHOP and CVP, was found to enhance ADCC sensitivity of RL cells and resistant clones and to significantly suppress tumor growth in xenograft models. Treatment with prednisolone upregulated expression of CD20 and an apoptosis-inducing protein BIM, which might augment perforin/granzyme B-mediated cell death. Furthermore, pretreatment of the effector cells with bendamustine enhanced ADCC activity, and treatment with obinutuzumab plus bendamustine showed significant antitumor efficacy in xenograft models. It was speculated that bendamustine upregulates ADCC activity by potentiating granules-mediated cell killing.

Conclusions

Our study revealed a novel mechanism underlying obinutuzumab-induced ADCC resistance and indicated that ADCC resistance could be overcome by combining obinutuzumab with prednisolone or bendamustine. This study provides a scientific rationale for obinutuzumab-retreatment in combination with clinically available chemotherapeutic agents for obinutuzumab resistant follicular lymphoma.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-022-07280-w.

A Fifteen-Gene Classifier to Predict Neoadjuvant Chemotherapy Responses in Patients with Stage IB to IIB Squamous Cervical Cancer

Neoadjuvant chemotherapy (NACT) remains an attractive alternative for controlling locally advanced cervical cancer. However, approximately 15-34% of women do not respond to induction therapy. To develop a risk stratification tool, 56 patients with stage IB-IIB cervical cancer are included in 2 research centers from the discovery cohort. Patient-specific somatic mutations led to NACT non-responsiveness are identified by whole-exome sequencing. Next, CRISPR/Cas9-based library screenings are performed based on these genes to confirm their biological contribution to drug resistance. A 15-gene classifier is developed by generalized linear regression analysis combined with the logistic regression model. In an independent validation cohort of 102 patients, the classifier showed good predictive ability with an area under the curve of 0.80 (95% confidence interval (CI), 0.69-0.91). Furthermore, the 15-gene classifier is significantly associated with patient responsiveness to NACT in both univariate (odds ratio, 10.8; 95% CI, 3.55-32.86; p = 2.8 × 10-5) and multivariate analysis (odds ratio, 17.34; 95% CI, 4.04-74.40; p = 1.23 × 10-4) in the validation set. In conclusion, the 15-gene classifier can accurately predict the clinical response to NACT before treatment, representing a promising approach for guiding the selection of appropriate treatment strategies for locally advanced cervical cancer.

Obinutuzumab in Combination with Chemotherapy Enhances Direct Cell Death in CD20-Positive Obinutuzumab-resistant Non-Hodgkin Lymphoma Cells

Follicular lymphoma commonly recurs and is difficult to cure. Obinutuzumab is a humanized glycoengineered type II anti-CD20 antibody with a mode of action that includes induction of antibody-dependent cellular cytotoxicity, antibody-dependent cellular phagocytosis, and direct cell death. There is no evidence on the effectiveness of retreatment with obinutuzumab in patients with prior obinutuzumab treatment. Using obinutuzumab-induced direct-cell-death-resistant cells, we investigated the efficacy of obinutuzumab retreatment in combination with chemotherapeutic agents used in follicular lymphoma treatment. Human non-Hodgkin lymphoma SU-DHL-4 cells were sustainably exposed to obinutuzumab in vitro, and 17 resistant clones expressing CD20 and showing 100-fold higher IC₅₀ of obinutuzumab than parental cells were established. The growth inhibition effect of obinutuzumab in combination with bendamustine, 4-hydroperoxy-cyclophosphamide, doxorubicin, vincristine, or prednisolone was estimated using an interaction index based on the Bliss independence model. For each clone, there were various combinations of obinutuzumab and chemotherapeutic agents that showed supra-additive effects. Obinutuzumab combined with doxorubicin enhanced caspase-dependent apoptosis and growth inhibition effect. Obinutuzumab combined with prednisolone enhanced DNA fragmentation and G₀-G₁ arrest. These combinations also had an antitumor effect in mouse xenograft models. Our results indicate that retreatment with obinutuzumab, when it is combined with chemotherapeutic agents, is effective in the CD20-positive obinutuzumab-induced direct-cell-death-resistant cells.

Understanding How Genetic Mutations Collaborate with Genomic Instability in Cancer

Chromosomal instability is the process of mis-segregation for ongoing chromosomes, which leads to cells with an abnormal number of chromosomes, also known as an aneuploid state. Induced aneuploidy is detrimental during development and in primary cells but aneuploidy is also a hallmark of cancer cells. It is therefore believed that premalignant cells need to overcome aneuploidy-imposed stresses to become tumorigenic. Over the past decade, some aneuploidy-tolerating pathways have been identified through small-scale screens, which suggest that aneuploidy tolerance pathways can potentially be therapeutically exploited. However, to better understand the processes that lead to aneuploidy tolerance in cancer cells, large-scale and unbiased genetic screens are needed, both in euploid and aneuploid cancer models. In this review, we describe some of the currently known aneuploidy-tolerating hits, how large-scale genome-wide screens can broaden our knowledge on aneuploidy specific cancer driver genes, and how we can exploit the outcomes of these screens to improve future cancer therapy.

From Phenotypic Hit to Chemical Probe: Chemical Biology Approaches to Elucidate Small Molecule Action in Complex Biological Systems

Biologically active small molecules have a central role in drug development, and as chemical probes and tool compounds to perturb and elucidate biological processes. Small molecules can be rationally designed for a given target, or a library of molecules can be screened against a target or phenotype of interest. Especially in the case of phenotypic screening approaches, a major challenge is to translate the compound-induced phenotype into a well-defined cellular target and mode of action of the hit compound. There is no "one size fits all" approach, and recent years have seen an increase in available target deconvolution strategies, rooted in organic chemistry, proteomics, and genetics. This review provides an overview of advances in target identification and mechanism of action studies, describes the strengths and weaknesses of the different approaches, and illustrates the need for chemical biologists to integrate and expand the existing tools to increase the probability of evolving screen hits to robust chemical probes.

MAP2K1 is a potential therapeutic target in erlotinib resistant head and neck squamous cell carcinoma

Epidermal growth factor receptor (EGFR) targeted therapies have shown limited efficacy in head and neck squamous cell carcinoma (HNSCC) patients despite its overexpression. Identifying molecular mechanisms associated with acquired resistance to EGFR-TKIs such as erlotinib remains an unmet need and a therapeutic challenge. In this study, we employed an integrated multi-omics approach to delineate mechanisms associated with acquired resistance to erlotinib by carrying out whole exome sequencing, quantitative proteomic and phosphoproteomic profiling. We observed amplification of several genes including AXL kinase and transcription factor YAP1 resulting in protein overexpression. We also observed expression of constitutively active mutant MAP2K1 (p.K57E) in erlotinib resistant SCC-R cells. An integrated analysis of genomic, proteomic and phosphoproteomic data revealed alterations in MAPK pathway and its downstream targets in SCC-R cells. We demonstrate that erlotinib-resistant cells are sensitive to MAPK pathway inhibition. This study revealed multiple genetic, proteomic and phosphoproteomic alterations associated with erlotinib resistant SCC-R cells. Our data indicates that therapeutic targeting of MAPK pathway is an effective strategy for treating erlotinib-resistant HNSCC tumors.

Coupling chemical mutagenesis to next generation sequencing for the identification of drug resistance mutations in Leishmania

Current genome-wide screens allow system-wide study of drug resistance but detecting small nucleotide variants (SNVs) is challenging. Here, we use chemical mutagenesis, drug selection and next generation sequencing to characterize miltefosine and paromomycin resistant clones of the parasite Leishmania. We highlight several genes involved in drug resistance by sequencing the genomes of 41 resistant clones and by concentrating on recurrent SNVs. We associate genes linked to lipid metabolism or to ribosome/translation functions with miltefosine or paromomycin resistance, respectively. We prove by allelic replacement and CRISPR-Cas9 gene-editing that the essential protein kinase CDPK1 is crucial for paromomycin resistance. We have linked CDPK1 in translation by functional interactome analysis, and provide evidence that CDPK1 contributes to antimonial resistance in the parasite. This screen is powerful in exploring networks of drug resistance in an organism with diploid to mosaic aneuploid genome, hence widening the scope of its applicability.

Here, Bhattacharya et al. chemically mutagenize Leishmania and identify genes associated with resistance to miltefosine and paromomycin by next generation sequencing. The study shows that a protein kinase (CDPK1) can mediate resistance to paromomycin by affecting translation.

Chemogenomic Screen for Imipenem Resistance in Gram-Negative Bacteria

Gram-negative carbapenem-resistant bacteria are a major threat to global health. The use of genome-wide screening approaches to probe for genes or mutations enabling resistance can lead to identification of molecular markers for diagnostics applications. We describe an approach called Mut-Seq that couples chemical mutagenesis and next-generation sequencing for studying resistance to imipenem in the Gram-negative bacteria Escherichia coli, Klebsiella pneumoniae, and Pseudomonas aeruginosa. The use of this approach highlighted shared and species-specific responses, and the role in resistance of a number of genes involved in membrane biogenesis, transcription, and signal transduction was functionally validated. Interestingly, some of the genes identified were previously considered promising therapeutic targets. Our genome-wide screen has the potential to be extended outside drug resistance studies and expanded to other organisms.

Carbapenem-resistant Gram-negative bacteria are considered a major threat to global health. Imipenem (IMP) is used as a last line of treatment against these pathogens, but its efficacy is diminished by the emergence of resistance. We applied a whole-genome screen in Escherichia coli, Klebsiella pneumoniae, and Pseudomonas aeruginosa isolates that were submitted to chemical mutagenesis, selected for IMP resistance, and characterized by next-generation sequencing. A comparative analysis of IMP-resistant clones showed that most of the highly mutated genes shared by the three species encoded proteins involved in transcription or signal transduction. Of these, the rpoD gene was one of the most prevalent and an E. coli strain disrupted for rpoD displayed a 4-fold increase in resistance to IMP. E. coli and K. pneumoniae also specifically shared several mutated genes, most involved in membrane/cell envelope biogenesis, and the contribution in IMP susceptibility was experimentally proven for amidases, transferases, and transglycosidases. P. aeruginosa differed from the two Enterobacteriaceae isolates with two different resistance mechanisms, with one involving mutations in the oprD porin or, alternatively, in two-component systems. Our chemogenomic screen performed with the three species has highlighted shared and species-specific responses to IMP.

IMPORTANCE Gram-negative carbapenem-resistant bacteria are a major threat to global health. The use of genome-wide screening approaches to probe for genes or mutations enabling resistance can lead to identification of molecular markers for diagnostics applications. We describe an approach called Mut-Seq that couples chemical mutagenesis and next-generation sequencing for studying resistance to imipenem in the Gram-negative bacteria Escherichia coli, Klebsiella pneumoniae, and Pseudomonas aeruginosa. The use of this approach highlighted shared and species-specific responses, and the role in resistance of a number of genes involved in membrane biogenesis, transcription, and signal transduction was functionally validated. Interestingly, some of the genes identified were previously considered promising therapeutic targets. Our genome-wide screen has the potential to be extended outside drug resistance studies and expanded to other organisms.

Sensitivity and Resistance of MET Exon 14 Mutations in Lung Cancer to Eight MET Tyrosine Kinase Inhibitors In Vitro

BACKGROUND

MNNG HOS transforming gene (MET) exon 14 mutations in lung cancer, including exon 14 skipping and point mutations, have been attracting the attention of thoracic oncologists as new therapeutic targets. Tumors with these mutations almost always acquire resistance, which also occurs in other oncogene-addicted lung cancers. However, the resistance mechanisms and treatment strategies are not fully understood.

METHODS

We generated Ba/F3 cells expressing MET exon 14 mutations by retroviral gene transfer. The sensitivities of these cells to eight MET-tyrosine kinase inhibitors (TKIs) were determined using a colorimetric assay. In addition, using N-ethyl-N-nitrosourea mutagenesis, we generated resistant clones, searched for secondary MET mutations, and then examined the sensitivities of these resistant cells to different TKIs.

RESULTS

Ba/F3 cells transfected with MET mutations grew in the absence of interleukin-3, indicating their oncogenic activity. These cells were sensitive to all MET-TKIs except tivantinib. We identified a variety of secondary mutations. D1228 and Y1230 were common sites for resistance mutations for type I TKIs, which bind the active form of MET, whereas L1195 and F1200 were common sites for type II TKIs, which bind the inactive form. In general, resistance mutations against type I were sensitive to type II, and vice versa.

CONCLUSIONS

MET-TKIs inhibited the growth of cells with MET exon 14 mutations. We also identified mutation sites specific for TKI types as resistance mechanisms and complementary activities between type I and type II inhibitors against those mutations. These finding should provide relevant clinical implication for treating patients with lung cancer harboring MET exon 14 mutations.

Diagnostic Value of MiR-125b as a Potential Biomarker for Stage I Lung Adenocarcinoma

Background:

We aimed at exploring biological functions of differentially expressed miRNAs during carcinogenesis, to identify miRNAs dysegulations involved in DNA repair mechanisms, and to evaluate potential of miRNAs as prognostic and diagnostic biomarkers for early lung adenocarcinomas (LAC).

Methods:

We obtained 21 LAC and paired adjacent normal formalin-fixed, paraffinembedded lung tissues from patients who underwent curative resection for stage I LAC. We compared expression levels of eight miRNAs involved in the DNA repair mechanism between LAC and adjacent tissues.

Results:

Expressions of Hsa-miR-9-5p, hsa-miR-24-3p, hsa-miR-125a-3p, hsa-miR- 125b-5p, hsa-miR-155-5p, and hsa-let-7a-5p were significantly up-regulated in stage I LAC tissues compared with those in the adjacent tissues. In addition, expressions of hsa-mir-9-5p, hsa-mir-24-3p, hsa-mir-125a-3p, hsa-mir-125b-5p, and hsa-mir-155-5p were significantly up-regulated in stage Ia LAC tissues, whereas expressions of hsa-mir- 125a-3p and hsa-mir-125b-5p were significantly up-regulated in stage Ib LAC tissues. Receiver operating characteristic (ROC) analysis revealed that AUROC of hsa-mir-125b- 5p was 0.875 (P < 0.001).

Conclusion:

Expression of hsa-mir-125b-5p could be used to distinguish LAC from adjacent tissues. Our result suggests that hsa-mir125b-5p can be a prognostic and diagnostic biomarker for LAC.

Gain- and Loss-of-Function Screens Coupled to Next-Generation Sequencing for Antibiotic Mode of Action and Resistance Studies in Streptococcus pneumoniae

Two whole-genome screening approaches are described for studying the mode of action and the mechanisms of resistance to trimethoprim (TMP) in the Gram-positive Streptococcus pneumoniae The gain-of-function approach (Int-Seq) relies on a genomic library of DNA fragments integrated into a fucose-inducible cassette. The second approach, leading to both gain- and loss-of-function mutation, is based on chemical mutagenesis coupled to next-generation sequencing (Mut-Seq). Both approaches pointed at the drug target dihydrofolate reductase (DHFR) as a major resistance mechanism to TMP. Resistance was achieved by dhfr overexpression either through the addition of fucose (Int-Seq) or by mutations upstream of the gene (Mut-Seq). Three types of mutations increased expression by disrupting a predicted Rho-independent terminator upstream of dhfr Known and novel DHFR mutations were also detected by Mut-Seq, and these were functionally validated for TMP resistance. The two approaches also suggested that an increase in the metabolic flux from purine synthesis to GTP and then to folate can modulate the susceptibility to TMP. Finally, we provide evidence for a novel role of the ABC transporter PatAB in TMP susceptibility. Our genomic screens highlighted novel aspects on the mode of action and mechanisms of resistance to antibiotics.

Evolution of Relapse-Proficient Subclones Constrained by Collateral Sensitivity to Oncogene Overdose in Wnt-Driven Mammary Cancer

Targeted cancer therapeutics select for drug-resistant rescue subclones (RSCs), which typically carry rescue mutations that restore oncogenic signaling. Whereas mutations underlying antibiotic resistance frequently burden drug-naive microbes with a fitness cost, it remains unknown whether and how rescue mutations underlying cancer relapse encounter negative selection prior to targeted therapy. Here, using mouse models of reversible, Wnt-driven mam-mary cancer, we uncovered stringent counter-selection against Wnt signaling overdose during the clonal evolution of RSCs. Analyzing recurrent tumors emerging during simulated targeted therapy (Wnt withdrawal) by multi-region DNA sequencing revealed polyclonal relapses comprised of multiple RSCs, which bear distinct but functionally equivalent rescue mutations that converge on sub-maximal Wnt pathway activation. When superimposed on native (i.e., undrugged) signaling, these rescue mutations faced negative selection, indicating that they burden RSCs with a fitness cost before Wnt withdrawal unmasks their selective advantage. Exploiting collateral sensitivity to oncogene overdose may help eliminate RSCs and prevent cancer relapse.

Keller and Gunther show that Wnt-driven mammary cancers challenged with simulated targeted therapy (Wnt withdrawal) undergo clonal evolution, which stringently selects for mutations that restore a “just right” level of oncogenic signaling. Therefore, cancer relapses emerge from rare subclones that are encumbered by an untapped vulnerability to oncogene overdose.

Mouse models in the era of large human tumour sequencing studies

Cancer is a complex disease in which cells progressively accumulate mutations disrupting their cellular processes. A fraction of these mutations drive tumourigenesis by affecting oncogenes or tumour suppressor genes, but many mutations are passengers with no clear contribution to tumour development. The advancement of DNA and RNA sequencing technologies has enabled in-depth analysis of thousands of human tumours from various tissues to perform systematic characterization of their (epi)genomes and transcriptomes in order to identify (epi)genetic changes associated with cancer. Combined with considerable progress in algorithmic development, this expansion in scale has resulted in the identification of many cancer-associated mutations, genes and pathways that are considered to be potential drivers of tumour development. However, it remains challenging to systematically identify drivers affected by complex genomic rearrangements and drivers residing in non-coding regions of the genome or in complex amplicons or deletions of copy-number driven tumours. Furthermore, functional characterization is challenging in the human context due to the lack of genetically tractable experimental model systems in which the effects of mutations can be studied in the context of their tumour microenvironment. In this respect, mouse models of human cancer provide unique opportunities for pinpointing novel driver genes and their detailed characterization. In this review, we provide an overview of approaches for complementing human studies with data from mouse models. We also discuss state-of-the-art technological developments for cancer gene discovery and validation in mice.

Targeted inhibition of histone H3K27 demethylation is effective in high-risk neuroblastoma

High-risk neuroblastoma is often distinguished by amplification of MYCN and loss of differentiation potential. We performed high-throughput drug screening of epigenetic-targeted therapies across a large and diverse tumor cell line panel and uncovered the hypersensitivity of neuroblastoma cells to GSK-J4, a small-molecule dual inhibitor of lysine 27 of histone 3 (H3K27) demethylases ubiquitously transcribed tetratricopeptide repeat, X chromosome (UTX), and histone demethylase Jumonji D3 (JMJD3). Mechanistically, GSK-J4 induced neuroblastoma differentiation and endoplasmic reticulum (ER) stress, with accompanying up-regulation of p53 up-regulated modulator of apoptosis (PUMA) and induction of cell death. Retinoic acid (RA)-resistant neuroblastoma cells were sensitive to GSK-J4. In addition, GSK-J4 was effective at blocking the growth of chemorefractory and patient-derived xenograft models of high-risk neuroblastoma in vivo. Furthermore, GSK-J4 and RA combination increased differentiation and ER stress over GSK-J4 effects and limited the growth of neuroblastomas resistant to either drug alone. In MYCN-amplified neuroblastoma, PUMA induction by GSK-J4 sensitized tumors to the B cell lymphoma 2 (BCL-2) inhibitor venetoclax, demonstrating that epigenetic-targeted therapies and BCL-2 homology domain 3 mimetics can be rationally combined to treat this high-risk subset of neuroblastoma. Therefore, H3K27 demethylation inhibition is a promising therapeutic target to treat high-risk neuroblastoma, and H3K27 demethylation can be part of rational combination therapies to induce robust antineuroblastoma activity.

Pathway-based dissection of the genomic heterogeneity of cancer hallmarks' acquisition with SLAPenrich

Cancer hallmarks are evolutionary traits required by a tumour to develop. While extensively characterised, the way these traits are achieved through the accumulation of somatic mutations in key biological pathways is not fully understood. To shed light on this subject, we characterised the landscape of pathway alterations associated with somatic mutations observed in 4,415 patients across ten cancer types, using 374 orthogonal pathway gene-sets mapped onto canonical cancer hallmarks. Towards this end, we developed SLAPenrich: a computational method based on population-level statistics, freely available as an open source R package. Assembling the identified pathway alterations into sets of hallmark signatures allowed us to connect somatic mutations to clinically interpretable cancer mechanisms. Further, we explored the heterogeneity of these signatures, in terms of ratio of altered pathways associated with each individual hallmark, assuming that this is reflective of the extent of selective advantage provided to the cancer type under consideration. Our analysis revealed the predominance of certain hallmarks in specific cancer types, thus suggesting different evolutionary trajectories across cancer lineages. Finally, although many pathway alteration enrichments are guided by somatic mutations in frequently altered high-confidence cancer genes, excluding these driver mutations preserves the hallmark heterogeneity signatures, thus the detected hallmarks' predominance across cancer types. As a consequence, we propose the hallmark signatures as a ground truth to characterise tails of infrequent genomic alterations and identify potential novel cancer driver genes and networks.

Identification of Mutation Accumulation as Resistance Mechanism Emerging in First-Line Osimertinib Treatment

INTRODUCTION

The survival of patients with EGFR mutation-positive lung cancer has dramatically improved since the introduction of EGFR tyrosine kinase inhibitors (EGFR-TKIs). Recently, osimertinib showed significantly prolonged progression-free survival than first-generation EGFR-TKI in first-line treatment, suggesting that a paradigm change that would move osimetinib to first-line treatment is indicated. We performed N-ethyl-N-nitrosourea (ENU) mutagenesis screening to uncover the resistant mechanism in first- and second-line osimertinib treatment.

METHODS

Ba/F3 cells harboring EGFR activating-mutation with or without secondary resistant mutation were exposed to ENU for 24 hours to introduce random mutations and selected with gefitinib, afatinib, or osimertinib. Mutations of emerging resistant cells were assessed.

RESULTS

The resistance of T790M and C797S to gefitinib and osimertinib, respectively, was prevalent in the mutagenesis screening with the Ba/F3 cells harboring activating-mutation alone. From C797S/activating-mutation expressing Ba/F3, the additional T790M was a major resistant mechanism in gefitinib and afatinib selection and the additional T854A and L792H were minor resistance mechanisms only in afatinib selection. However, the additional T854A or L792H mediated resistance to all classes of EGFR-TKI. Surprisingly, no resistant clone due to secondary mutation emerged from activating-mutation alone in the gefitinib + osimertinib selection.

CONCLUSIONS

We showed the resistance mechanism to EGFR-TKI focusing on first- and second-line osimertinib using ENU mutagenesis screening. Additional T854A and L792H on C797S/activating-mutation were found as afatinib resistance and not as gefitinib resistance. Thus, compared to afatinib, the first-generation EGFR-TKI might be preferable as second-line treatment to C797S/activating-mutation emerging after first-line osimertinib treatment.

The Pursuit of Therapeutic Biomarkers with High-Throughput Cancer Cell Drug Screens

In the last decade we have witnessed tremendous advances in our understanding of the landscape of the molecular alterations that underpin many of the most prevalent cancers, in the use of automated high-throughput platforms for high-throughput drug screens in cancer cells, in the creation of more clinically relevant cancer cell models, and lastly in the development of more useful computational approaches in the pursuit of biomarkers of drug response. Separately, each of these improvements will undoubtedly lead to improvements in the treatment of cancer patients but to fulfill the promise of truly personalized cancer medicine, we must bring these disciplines together in a truly multidisciplinary fashion.