An integrated model for predicting KRAS dependency.

The clinical approvals of KRAS G12C inhibitors have been a revolutionary advance in precision oncology, but response rates are often modest. To improve patient selection, we developed an integrated model to predict KRAS dependency. By integrating molecular profiles of a large panel of cell lines from the DEMETER2 dataset, we built a binary classifier to predict a tumor's KRAS dependency. Monte Carlo cross validation via ElasticNet within the training set was used to compare model performance and to tune parameters α and λ. The final model was then applied to the validation set. We validated the model with genetic depletion assays and an external dataset of lung cancer cells treated with a G12C inhibitor. We then applied the model to several Cancer Genome Atlas (TCGA) datasets. The final "K20" model contains 20 features, including expression of 19 genes and KRAS mutation status. In the validation cohort, K20 had an AUC of 0.94 and accurately predicted KRAS dependency in both mutant and KRAS wild-type cell lines following genetic depletion. It was also highly predictive across an external dataset of lung cancer lines treated with KRAS G12C inhibition. When applied to TCGA datasets, specific subpopulations such as the invasive subtype in colorectal cancer and copy number high pancreatic adenocarcinoma were predicted to have higher KRAS dependency. The K20 model has simple yet robust predictive capabilities that may provide a useful tool to select patients with KRAS mutant tumors that are most likely to respond to direct KRAS inhibitors.

Translating transcriptomic findings from cancer model systems to humans through joint dimension reduction.

Model systems are an essential resource in cancer research. They simulate effects that we can infer into humans, but come at a risk of inaccurately representing human biology. This inaccuracy can lead to inconclusive experiments or misleading results, urging the need for an improved process for translating model system findings into human-relevant data. We present a process for applying joint dimension reduction (jDR) to horizontally integrate gene expression data across model systems and human tumor cohorts. We then use this approach to combine human TCGA gene expression data with data from human cancer cell lines and mouse model tumors. By identifying the aspects of genomic variation joint-acting across cohorts, we demonstrate how predictive modeling and clinical biomarkers from model systems can be improved.

Rad18 mediates specific mutational signatures and shapes the genomic landscape of carcinogen-induced tumors in vivo.

The E3 ubiquitin ligase Rad18 promotes a damage-tolerant and error-prone mode of DNA replication termed trans-lesion synthesis that is pathologically activated in cancer. However, the impact of vertebrate Rad18 on cancer genomes is not known. To determine how Rad18 affects mutagenesis in vivo, we have developed and implemented a novel computational pipeline to analyze genomes of carcinogen (7, 12-Dimethylbenz[a]anthracene, DMBA)-induced skin tumors from Rad18+/+ and Rad18- / - mice. We show that Rad18 mediates specific mutational signatures characterized by high levels of A(T)>T(A) single nucleotide variations (SNVs). In Rad18- /- tumors, an alternative mutation pattern arises, which is characterized by increased numbers of deletions >4 bp. Comparison with annotated human mutational signatures shows that COSMIC signature 22 predominates in Rad18+/+ tumors whereas Rad18- / - tumors are characterized by increased contribution of COSMIC signature 3 (a hallmark of BRCA-mutant tumors). Analysis of The Cancer Genome Atlas shows that RAD18 expression is strongly associated with high SNV burdens, suggesting RAD18 also promotes mutagenesis in human cancers. Taken together, our results show Rad18 promotes mutagenesis in vivo, modulates DNA repair pathway choice in neoplastic cells, and mediates specific mutational signatures that are present in human tumors.

Genetic determinants of cellular addiction to DNA polymerase theta.

Polymerase theta (Pol θ, gene name Polq) is a widely conserved DNA polymerase that mediates a microhomology-mediated, error-prone, double strand break (DSB) repair pathway, referred to as Theta Mediated End Joining (TMEJ). Cells with homologous recombination deficiency are reliant on TMEJ for DSB repair. It is unknown whether deficiencies in other components of the DNA damage response (DDR) also result in Pol θ addiction. Here we use a CRISPR genetic screen to uncover 140 Polq synthetic lethal (PolqSL) genes, the majority of which were previously unknown. Functional analyses indicate that Pol θ/TMEJ addiction is associated with increased levels of replication-associated DSBs, regardless of the initial source of damage. We further demonstrate that approximately 30% of TCGA breast cancers have genetic alterations in PolqSL genes and exhibit genomic scars of Pol θ/TMEJ hyperactivity, thereby substantially expanding the subset of human cancers for which Pol θ inhibition represents a promising therapeutic strategy.