Kras oncogene ablation prevents resistance in advanced lung adenocarcinomas

KRASG12C inhibitors have revolutionized the clinical management of patients with KRASG12C-mutant lung adenocarcinoma. However, patient exposure to these inhibitors leads to the rapid onset of resistance. In this study, we have used genetically engineered mice to compare the therapeutic efficacy and the emergence of tumor resistance between genetic ablation of mutant Kras expression and pharmacological inhibition of oncogenic KRAS activity. Whereas Kras ablation induces massive tumor regression and prevents the appearance of resistant cells in vivo, treatment of KrasG12C/Trp53-driven lung adenocarcinomas with sotorasib, a selective KRASG12C inhibitor, caused a limited antitumor response similar to that observed in the clinic, including the rapid onset of resistance. Unlike in human tumors, we did not observe mutations in components of the RAS-signaling pathways. Instead, sotorasib-resistant tumors displayed amplification of the mutant Kras allele and activation of xenobiotic metabolism pathways, suggesting that reduction of the on-target activity of KRASG12C inhibitors is the main mechanism responsible for the onset of resistance. In sum, our results suggest that resistance to KRAS inhibitors could be prevented by achieving a more robust inhibition of KRAS signaling mimicking the results obtained upon Kras ablation.

[1]  M. Barbacid,et al.  KRAS4A induces metastatic lung adenocarcinomas in vivo in the absence of the KRAS4B isoform , 2021, Proceedings of the National Academy of Sciences.

[2]  Giovanni Parmigiani,et al.  ComBat-seq: batch effect adjustment for RNA-seq count data , 2020, bioRxiv.

[3]  J. Flowers,et al.  Origins and geographic diversification of African rice (Oryza glaberrima) , 2018, bioRxiv.

[4]  M. Barbacid,et al.  c-RAF Ablation Induces Regression of Advanced Kras/Trp53 Mutant Lung Adenocarcinomas by a Mechanism Independent of MAPK Signaling. , 2018, Cancer cell.

[5]  D. Pisano,et al.  Nextpresso: Next Generation Sequencing Expression Analysis Pipeline , 2017, Current Bioinformatics.

[6]  Eric Talevich,et al.  CNVkit: Genome-Wide Copy Number Detection and Visualization from Targeted DNA Sequencing , 2016, PLoS Comput. Biol..

[7]  F. Cunningham,et al.  The Ensembl Variant Effect Predictor , 2016, bioRxiv.

[8]  J. Cigudosa,et al.  Truncated RUNX1 protein generated by a novel t(1;21)(p32;q22) chromosomal translocation impairs the proliferation and differentiation of human hematopoietic progenitors , 2016, Oncogene.

[9]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[10]  Heng Li.  Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM , 2013, 1303.3997.

[11]  Justin Guinney,et al.  GSVA: gene set variation analysis for microarray and RNA-Seq data , 2013, BMC Bioinformatics.

[12]  Jeremy Wazny,et al.  Xenome—a tool for classifying reads from xenograft samples , 2012, Bioinform..

[13]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[14]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[15]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[16]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[17]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[18]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[19]  A. Olshen,et al.  A Faster Circular Binary Segmentation Algorithm for the Analysis of Array Cgh Data , 2022 .