Acquired resistance to anti-PD1 therapy in patients with NSCLC associates with immunosuppressive T cell phenotype.

Immune checkpoint inhibitor treatment has the potential to prolong survival in non-small cell lung cancer (NSCLC), however, some of the patients develop resistance following initial response. Here, we analyze the immune phenotype of matching tumor samples from a cohort of NSCLC patients showing good initial response to immune checkpoint inhibitors, followed by acquired resistance at later time points. By using imaging mass cytometry and whole exome and RNA sequencing, we detect two patterns of resistance¨: One group of patients is characterized by reduced numbers of tumor-infiltrating CD8+ T cells and reduced expression of PD-L1 after development of resistance, whereas the other group shows high CD8+ T cell infiltration and high expression of PD-L1 in addition to markedly elevated expression of other immune-inhibitory molecules. In two cases, we detect downregulation of type I and II IFN pathways following progression to resistance, which could lead to an impaired anti-tumor immune response. This study thus captures the development of immune checkpoint inhibitor resistance as it progresses and deepens our mechanistic understanding of immunotherapy response in NSCLC.

Microbial peptides activate tumour-infiltrating lymphocytes in glioblastoma.

Microbial organisms have key roles in numerous physiological processes in the human body and have recently been shown to modify the response to immune checkpoint inhibitors1,2. Here we aim to address the role of microbial organisms and their potential role in immune reactivity against glioblastoma. We demonstrate that HLA molecules of both glioblastoma tissues and tumour cell lines present bacteria-specific peptides. This finding prompted us to examine whether tumour-infiltrating lymphocytes (TILs) recognize tumour-derived bacterial peptides. Bacterial peptides eluted from HLA class II molecules are recognized by TILs, albeit very weakly. Using an unbiased antigen discovery approach to probe the specificity of a TIL CD4+ T cell clone, we show that it recognizes a broad spectrum of peptides from pathogenic bacteria, commensal gut microbiota and also glioblastoma-related tumour antigens. These peptides were also strongly stimulatory for bulk TILs and peripheral blood memory cells, which then respond to tumour-derived target peptides. Our data hint at how bacterial pathogens and bacterial gut microbiota can be involved in specific immune recognition of tumour antigens. The unbiased identification of microbial target antigens for TILs holds promise for future personalized tumour vaccination approaches.

Vaccination with Designed Neopeptides Induces Intratumoral, Cross-reactive CD4+ T-cell Responses in Glioblastoma.

PURPOSE: The low mutational load of some cancers is considered one reason for the difficulty to develop effective tumor vaccines. To overcome this problem, we developed a strategy to design neopeptides through single amino acid mutations to enhance their immunogenicity. EXPERIMENTAL DESIGN: Exome and RNA sequencing as well as in silico HLA-binding predictions to autologous HLA molecules were used to identify candidate neopeptides. Subsequently, in silico HLA-anchor placements were used to deduce putative T-cell receptor (TCR) contacts of peptides. Single amino acids of TCR contacting residues were then mutated by amino acid replacements. Overall, 175 peptides were synthesized and sets of 25 each containing both peptides designed to bind to HLA class I and II molecules applied in the vaccination. Upon development of a tumor recurrence, the tumor-infiltrating lymphocytes (TIL) were characterized in detail both at the bulk and clonal level. RESULTS: The immune response of peripheral blood T cells to vaccine peptides, including natural peptides and designed neopeptides, gradually increased with repetitive vaccination, but remained low. In contrast, at the time of tumor recurrence, CD8+ TILs and CD4+ TILs responded to 45% and 100%, respectively, of the vaccine peptides. Furthermore, TIL-derived CD4+ T-cell clones showed strong responses and tumor cell lysis not only against the designed neopeptide but also against the unmutated natural peptides of the tumor. CONCLUSIONS: Turning tumor self-peptides into foreign antigens by introduction of designed mutations is a promising strategy to induce strong intratumoral CD4+ T-cell responses in a cold tumor like glioblastoma.

Drug screening and genome editing in human pancreatic cancer organoids identifies drug-gene interactions and candidates for off-label treatment.

Pancreatic cancer (PDAC) is a highly aggressive malignancy for which the identification of novel therapies is urgently needed. Here, we establish a human PDAC organoid biobank from 31 genetically distinct lines, covering a representative range of tumor subtypes, and demonstrate that these reflect the molecular and phenotypic heterogeneity of primary PDAC tissue. We use CRISPR-Cas9 genome editing and drug screening to characterize drug-gene interactions with ARID1A and BRCA2. We find that missense- but not frameshift mutations in the PDAC driver gene ARID1A are associated with increased sensitivity to the kinase inhibitors dasatinib (p < 0.0001) and VE-821 (p < 0.0001). We conduct an automated drug-repurposing screen with 1,172 FDA-approved compounds, identifying 26 compounds that effectively kill PDAC organoids, including 19 chemotherapy drugs currently approved for other cancer types. We validate the activity of these compounds in vitro and in vivo. The in vivo validated hits include emetine and ouabain, compounds which are approved for non-cancer indications and which perturb the ability of PDAC organoids to respond to hypoxia. Our study provides proof-of-concept for advancing precision oncology and identifying candidates for drug repurposing via genome editing and drug screening in tumor organoid biobanks.

Diverse tumorigenic consequences of human papillomavirus integration in primary oropharyngeal cancers.

Human papillomavirus (HPV) causes 5% of all cancers and frequently integrates into host chromosomes. The HPV oncoproteins E6 and E7 are necessary but insufficient for cancer formation, indicating that additional secondary genetic events are required. Here, we investigate potential oncogenic impacts of virus integration. Analysis of 105 HPV-positive oropharyngeal cancers by whole-genome sequencing detects virus integration in 77%, revealing five statistically significant sites of recurrent integration near genes that regulate epithelial stem cell maintenance (i.e., SOX2, TP63, FGFR, MYC) and immune evasion (i.e., CD274). Genomic copy number hyperamplification is enriched 16-fold near HPV integrants, and the extent of focal host genomic instability increases with their local density. The frequency of genes expressed at extreme outlier levels is increased 86-fold within ±150 kb of integrants. Across 95% of tumors with integration, host gene transcription is disrupted via intragenic integrants, chimeric transcription, outlier expression, gene breaking, and/or de novo expression of noncoding or imprinted genes. We conclude that virus integration can contribute to carcinogenesis in a large majority of HPV-positive oropharyngeal cancers by inducing extensive disruption of host genome structure and gene expression.

The Tumor Profiler Study: integrated, multi-omic, functional tumor profiling for clinical decision support.

The application and integration of molecular profiling technologies create novel opportunities for personalized medicine. Here, we introduce the Tumor Profiler Study, an observational trial combining a prospective diagnostic approach to assess the relevance of in-depth tumor profiling to support clinical decision-making with an exploratory approach to improve the biological understanding of the disease.

Tracing Clonal Dynamics Reveals that Two- and Three-dimensional Patient-derived Cell Models Capture Tumor Heterogeneity of Clear Cell Renal Cell Carcinoma.

BACKGROUND: Extensive DNA sequencing has led to an unprecedented view of the diversity of individual genomes and their evolution among patients with clear cell renal cell carcinoma (ccRCC). OBJECTIVE: To understand subclonal architecture and dynamics of patient-derived two-dimensional (2D) and three-dimensional (3D) ccRCC models in vitro, in order to determine whether they mirror ccRCC inter- and intratumor heterogeneity. DESIGN, SETTING, AND PARTICIPANTS: We have established a comprehensive platform of living renal cancer cell models from ccRCC surgical specimens. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: We confirmed the concordance of 2D and 3D patient-derived cell (PDC) models with the original tumor tissue in terms of histology, biomarker expression, cancer driver mutations, and copy number alterations. We addressed inter- and intrapatient heterogeneity by analyzing clonal dynamics during serial passaging. RESULTS AND LIMITATIONS: In-depth genetic characterization verified the presence of heterogeneous cell populations, and revealed a high degree of similarity between subclonal compositions of monolayer and organoid cell cultures and the corresponding parental ccRCCs. Clonal dynamics were evident during serial passaging of cells in vitro, suggesting that PDC cultures can offer insights into evolutionary potential and treatment susceptibility of ccRCC subclones in vivo. Proof-of-concept drug profiling using selected ccRCC-targeted therapy agents highlighted patient-specific vulnerabilities in PDC models that could not be anticipated by interrogating commercially available cell lines. CONCLUSIONS: We demonstrate that PDC models mirror inter- and intratumor heterogeneity of ccRCC in vitro. Based on our findings, we envision that the use of these models will advance our understanding of the trajectories that cause genetic diversity and their consequences for treatment on an individual level. PATIENT SUMMARY: In this study, we developed two- and three-dimensional patient-derived models from clear cell renal cell carcinoma (ccRCC) as "mini-tumors in a dish." We show that these cell models retain important features of the human ccRCCs such as the profound tumor heterogeneity, thus highlighting their importance for cancer research and precision medicine.

Convergent network effects along the axis of gene expression during prostate cancer progression.

BACKGROUND: Tumor-specific genomic aberrations are routinely determined by high-throughput genomic measurements. It remains unclear how complex genome alterations affect molecular networks through changing protein levels and consequently biochemical states of tumor tissues. RESULTS: Here, we investigate the propagation of genomic effects along the axis of gene expression during prostate cancer progression. We quantify genomic, transcriptomic, and proteomic alterations based on 105 prostate samples, consisting of benign prostatic hyperplasia regions and malignant tumors, from 39 prostate cancer patients. Our analysis reveals the convergent effects of distinct copy number alterations impacting on common downstream proteins, which are important for establishing the tumor phenotype. We devise a network-based approach that integrates perturbations across different molecular layers, which identifies a sub-network consisting of nine genes whose joint activity positively correlates with increasingly aggressive tumor phenotypes and is predictive of recurrence-free survival. Further, our data reveal a wide spectrum of intra-patient network effects, ranging from similar to very distinct alterations on different molecular layers. CONCLUSIONS: This study uncovers molecular networks with considerable convergent alterations across tumor sites and patients. It also exposes a diversity of network effects: we could not identify a single sub-network that is perturbed in all high-grade tumor regions.

HLA-DR15 Molecules Jointly Shape an Autoreactive T Cell Repertoire in Multiple Sclerosis.

The HLA-DR15 haplotype is the strongest genetic risk factor for multiple sclerosis (MS), but our understanding of how it contributes to MS is limited. Because autoreactive CD4+ T cells and B cells as antigen-presenting cells are involved in MS pathogenesis, we characterized the immunopeptidomes of the two HLA-DR15 allomorphs DR2a and DR2b of human primary B cells and monocytes, thymus, and MS brain tissue. Self-peptides from HLA-DR molecules, particularly from DR2a and DR2b themselves, are abundant on B cells and thymic antigen-presenting cells. Furthermore, we identified autoreactive CD4+ T cell clones that can cross-react with HLA-DR-derived self-peptides (HLA-DR-SPs), peptides from MS-associated foreign agents (Epstein-Barr virus and Akkermansia muciniphila), and autoantigens presented by DR2a and DR2b. Thus, both HLA-DR15 allomorphs jointly shape an autoreactive T cell repertoire by serving as antigen-presenting structures and epitope sources and by presenting the same foreign peptides and autoantigens to autoreactive CD4+ T cells in MS.

ClinVAP: a reporting strategy from variants to therapeutic options.

MOTIVATION: Next-generation sequencing has become routine in oncology and opens up new avenues of therapies, particularly in personalized oncology setting. An increasing number of cases also implies a need for a more robust, automated and reproducible processing of long lists of variants for cancer diagnosis and therapy. While solutions for the large-scale analysis of somatic variants have been implemented, existing solutions often have issues with reproducibility, scalability and interoperability. RESULTS: Clinical Variant Annotation Pipeline (ClinVAP) is an automated pipeline which annotates, filters and prioritizes somatic single nucleotide variants provided in variant call format. It augments the variant information with documented or predicted clinical effect. These annotated variants are prioritized based on driver gene status and druggability. ClinVAP is available as a fully containerized, self-contained pipeline maximizing reproducibility and scalability allowing the analysis of larger scale data. The resulting JSON-based report is suited for automated downstream processing, but ClinVAP can also automatically render the information into a user-defined template to yield a human-readable report. AVAILABILITY AND IMPLEMENTATION: ClinVAP is available at https://github.com/PersonalizedOncology/ClinVAP. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comparative analysis of mRNA and protein degradation in prostate tissues indicates high stability of proteins.

Deterioration of biomolecules in clinical tissues is an inevitable pre-analytical process, which affects molecular measurements and thus potentially confounds conclusions from cohort analyses. Here, we investigate the degradation of mRNA and protein in 68 pairs of adjacent prostate tissue samples using RNA-Seq and SWATH mass spectrometry, respectively. To objectively quantify the extent of protein degradation, we develop a numerical score, the Proteome Integrity Number (PIN), that faithfully measures the degree of protein degradation. Our results indicate that protein degradation only affects 5.9% of the samples tested and shows negligible correlation with mRNA degradation in the adjacent samples. These findings are confirmed by independent analyses on additional clinical sample cohorts and across different mass spectrometric methods. Overall, the data show that the majority of samples tested are not compromised by protein degradation, and establish the PIN score as a generic and accurate indicator of sample quality for proteomic analyses.

Bioinformatics for precision oncology.

Molecular profiling of tumor biopsies plays an increasingly important role not only in cancer research, but also in the clinical management of cancer patients. Multi-omics approaches hold the promise of improving diagnostics, prognostics and personalized treatment. To deliver on this promise of precision oncology, appropriate bioinformatics methods for managing, integrating and analyzing large and complex data are necessary. Here, we discuss the specific requirements of bioinformatics methods and software that arise in the setting of clinical oncology, owing to a stricter regulatory environment and the need for rapid, highly reproducible and robust procedures. We describe the workflow of a molecular tumor board and the specific bioinformatics support that it requires, from the primary analysis of raw molecular profiling data to the automatic generation of a clinical report and its delivery to decision-making clinical oncologists. Such workflows have to various degrees been implemented in many clinical trials, as well as in molecular tumor boards at specialized cancer centers and university hospitals worldwide. We review these and more recent efforts to include other high-dimensional multi-omics patient profiles into the tumor board, as well as the state of clinical decision support software to translate molecular findings into treatment recommendations.

Increasing Dietary Fiber Intake Is Associated with a Distinct Esophageal Microbiome.

INTRODUCTION: There is increasing evidence that the microbiome contributes to esophageal disease. Diet, especially fiber and fat intake, is a known potent modifier of the colonic microbiome, but its impact on the esophageal microbiome is not well described. We hypothesized that dietary fiber and fat intake would be associated with a distinct esophageal microbiome. METHODS: We collected esophageal samples from 47 ambulatory patients scheduled to undergo endoscopy who completed a validated food frequency questionnaire quantifying dietary fiber and fat intake. Using 16S high-throughput sequencing, we determined composition of the esophageal microbiome and predicted functional capacity of microbiota based on fiber and fat intake. RESULTS: Among all samples, the most abundant phyla were Firmicutes (54.0%), Proteobacteria (19.0%), Bacteroidetes (17.0%), Actinobacteria (5.2%), and Fusobacteria (4.3%). Increasing fiber intake was significantly associated with increasing relative abundance of Firmicutes (p = 0.04) and decreasing relative abundance of Gram-negative bacteria overall (p = 0.03). Low fiber intake was associated with increased relative abundance of several Gram-negative bacteria, including Prevotella, Neisseria, and Eikenella. Several predicted metabolic pathways differed between highest and lowest quartile of fiber intake. Fat intake was associated with altered relative abundance of few taxa, with no alterations at the phylum level and no changes in microbiome functional composition. CONCLUSIONS: Dietary fiber, but not fat, intake was associated with a distinct esophageal microbiome. Diet should be considered an important modifier of the esophageal microbiome in future studies. Studies are also needed to elucidate how the effects of dietary fiber on the esophageal microbiome may contribute to esophageal disease.

SwissMTB: establishing comprehensive molecular cancer diagnostics in Swiss clinics.

BACKGROUND: Molecular precision oncology is an emerging practice to improve cancer therapy by decreasing the risk of choosing treatments that lack efficacy or cause adverse events. However, the challenges of integrating molecular profiling into routine clinical care are manifold. From a computational perspective these include the importance of a short analysis turnaround time, the interpretation of complex drug-gene and gene-gene interactions, and the necessity of standardized high-quality workflows. In addition, difficulties faced when integrating molecular diagnostics into clinical practice are ethical concerns, legal requirements, and limited availability of treatment options beyond standard of care as well as the overall lack of awareness of their existence. METHODS: To the best of our knowledge, we are the first group in Switzerland that established a workflow for personalized diagnostics based on comprehensive high-throughput sequencing of tumors at the clinic. Our workflow, named SwissMTB (Swiss Molecular Tumor Board), links genetic tumor alterations and gene expression to therapeutic options and clinical trial opportunities. The resulting treatment recommendations are summarized in a clinical report and discussed in a molecular tumor board at the clinic to support therapy decisions. RESULTS: Here we present results from an observational pilot study including 22 late-stage cancer patients. In this study we were able to identify actionable variants and corresponding therapies for 19 patients. Half of the patients were analyzed retrospectively. In two patients we identified resistance-associated variants explaining lack of therapy response. For five out of eleven patients analyzed before treatment the SwissMTB diagnostic influenced treatment decision. CONCLUSIONS: SwissMTB enables the analysis and clinical interpretation of large numbers of potentially actionable molecular targets. Thus, our workflow paves the way towards a more frequent use of comprehensive molecular diagnostics in Swiss hospitals.

Comprehensive Analysis of Alternative Splicing Across Tumors from 8,705 Patients.

Our comprehensive analysis of alternative splicing across 32 The Cancer Genome Atlas cancer types from 8,705 patients detects alternative splicing events and tumor variants by reanalyzing RNA and whole-exome sequencing data. Tumors have up to 30% more alternative splicing events than normal samples. Association analysis of somatic variants with alternative splicing events confirmed known trans associations with variants in SF3B1 and U2AF1 and identified additional trans-acting variants (e.g., TADA1, PPP2R1A). Many tumors have thousands of alternative splicing events not detectable in normal samples; on average, we identified ≈930 exon-exon junctions ("neojunctions") in tumors not typically found in GTEx normals. From Clinical Proteomic Tumor Analysis Consortium data available for breast and ovarian tumor samples, we confirmed ≈1.7 neojunction- and ≈0.6 single nucleotide variant-derived peptides per tumor sample that are also predicted major histocompatibility complex-I binders ("putative neoantigens").

Genomic and Geographic Context for the Evolution of High-Risk Carbapenem-Resistant Enterobacter cloacae Complex Clones ST171 and ST78.

Recent reports have established the escalating threat of carbapenem-resistant Enterobacter cloacae complex (CREC). Here, we demonstrate that CREC has evolved as a highly antibiotic-resistant rather than highly virulent nosocomial pathogen. Applying genomics and Bayesian phylogenetic analyses to a 7-year collection of CREC isolates from a northern Manhattan hospital system and to a large set of publicly available, geographically diverse genomes, we demonstrate clonal spread of a single clone, ST171. We estimate that two major clades of epidemic ST171 diverged prior to 1962, subsequently spreading in parallel from the Northeastern to the Mid-Atlantic and Midwestern United States and demonstrating links to international sites. Acquisition of carbapenem and fluoroquinolone resistance determinants by both clades preceded widespread use of these drugs in the mid-1980s, suggesting that antibiotic pressure contributed substantially to its spread. Despite a unique mobile repertoire, ST171 isolates showed decreased virulence in vitro While a second clone, ST78, substantially contributed to the emergence of CREC, it encompasses diverse carbapenemase-harboring plasmids, including a potentially hypertransmissible IncN plasmid, also present in other sequence types. Rather than heightened virulence, CREC demonstrates lineage-specific, multifactorial adaptations to nosocomial environments coupled with a unique potential to acquire and disseminate carbapenem resistance genes. These findings indicate a need for robust surveillance efforts that are attentive to the potential for local and international spread of high-risk CREC clones.IMPORTANCE Carbapenem-resistant Enterobacter cloacae complex (CREC) has emerged as a formidable nosocomial pathogen. While sporadic acquisition of plasmid-encoded carbapenemases has been implicated as a major driver of CREC, ST171 and ST78 clones demonstrate epidemic potential. However, a lack of reliable genomic references and rigorous statistical analyses has left many gaps in knowledge regarding the phylogenetic context and evolutionary pathways of successful CREC. Our reconstruction of recent ST171 and ST78 evolution represents a significant addition to current understanding of CREC and the directionality of its spread from the Eastern United States to the northern Midwestern United States with links to international collections. Our results indicate that the remarkable ability of E. cloacae to acquire and disseminate cross-class antibiotic resistance rather than virulence determinants, coupled with its ability to adapt under conditions of antibiotic pressure, likely led to the wide dissemination of CREC.

Genomic Surveillance Reveals Diversity of Multidrug-Resistant Organism Colonization and Infection: A Prospective Cohort Study in Liver Transplant Recipients.

Background: Multidrug-resistant organisms (MDROs) are an important cause of morbidity and mortality after solid organ transplantation. We aimed to characterize MDRO colonization dynamics and infection in liver transplant (LT) recipients through innovative use of active surveillance and whole-genome sequencing (WGS). Methods: We prospectively enrolled consecutive adult patients undergoing LT from March 2014 to March 2016. Fecal samples were collected at multiple timepoints from time of enrollment to 12 months posttransplant. Samples were screened for carbapenem-resistant Enterobacteriaceae (CRE), Enterobacteriaceae resistant to third-generation cephalosporins (Ceph-RE), and vancomycin-resistant enterococci. We performed WGS of CRE and selected Ceph-RE isolates. We also collected clinical data including demographics, transplant characteristics, and infection data. Results: We collected 998 stool samples and 119 rectal swabs from 128 patients. MDRO colonization was detected in 86 (67%) patients at least once and was significantly associated with subsequent MDRO infection (0 vs 19.8%, P = .002). Child-Turcotte-Pugh score at LT and duration of post-LT hospitalization were independent predictors of both MDRO colonization and infection. Temporal dynamics differed between MDROs with respect to onset of colonization, clearance, and infections. We detected an unexpected diversity of CRE colonizing isolates and previously unrecognized transmission that spanned Ceph-RE and CRE phenotypes, as well as a cluster of mcr-1-producing isolates. Conclusions: Active surveillance and WGS showed that MDRO colonization is a highly dynamic and complex process after LT. Understanding that complexity is crucial for informing decisions regarding MDRO infection control, use of therapeutic decolonization, and empiric treatment regimens.

Barrett's esophagus is associated with a distinct oral microbiome.

OBJECTIVES: The esophageal microbiome is composed of predominantly oral flora and is altered in reflux-related conditions including Barrett's esophagus (BE). Changes to the esophageal microbiome may be reflected in the oral cavity. Assessing the oral microbiome thus represents a potential non-invasive method to identify patients with BE. METHODS: Patients with and without BE undergoing upper endoscopy were prospectively enrolled. Demographics, clinical data, medications, and dietary intake were assessed. 16S rRNA gene sequencing was performed on saliva samples collected prior to endoscopy. Taxonomic differences between groups were assessed via linear discriminant analysis effect size (LEfSe). Logit models were used to develop microbiome signatures to distinguish BE from non-BE, assessed by area under the receiver operating curve (AUROC). RESULTS: A total of 49 patients were enrolled (control = 17, BE = 32). There was no significant difference in alpha diversity comparing all BE patients vs. CONTROLS: At the phylum level, the oral microbiome in BE patients had significantly increased relative abundance of Firmicutes (p = 0.005) and decreased Proteobacteria (p = 0.02). There were numerous taxonomic differences in the oral microbiome between BE and controls. A model including relative abundance of Lautropia, Streptococcus, and a genus in the order Bacteroidales distinguished BE from controls with an AUROC 0.94 (95% CI: 0.85-1.00). The optimal cutoff identified BE patients with 96.9% sensitivity and 88.2% specificity. CONCLUSIONS: The oral microbiome in BE patients was markedly altered and distinguished BE with relatively high accuracy. The oral microbiome represents a potential screening marker for BE, and validation studies in larger and distinct populations are warranted.

NGS-pipe: a flexible, easily extendable and highly configurable framework for NGS analysis.

Motivation: Next-generation sequencing is now an established method in genomics, and massive amounts of sequencing data are being generated on a regular basis. Analysis of the sequencing data is typically performed by lab-specific in-house solutions, but the agreement of results from different facilities is often small. General standards for quality control, reproducibility and documentation are missing. Results: We developed NGS-pipe, a flexible, transparent and easy-to-use framework for the design of pipelines to analyze whole-exome, whole-genome and transcriptome sequencing data. NGS-pipe facilitates the harmonization of genomic data analysis by supporting quality control, documentation, reproducibility, parallelization and easy adaptation to other NGS experiments. Availability and implementation: https://github.com/cbg-ethz/NGS-pipe. Contact: niko.beerenwinkel@bsse.ethz.ch.

Combined mutation in Vhl, Trp53 and Rb1 causes clear cell renal cell carcinoma in mice.

Clear cell renal cell carcinomas (ccRCCs) frequently exhibit inactivation of the von Hippel-Lindau tumor-suppressor gene, VHL, and often harbor multiple copy-number alterations in genes that regulate cell cycle progression. We show here that modeling these genetic alterations by combined deletion of Vhl, Trp53 and Rb1 specifically in renal epithelial cells in mice caused ccRCC. These tumors arose from proximal tubule epithelial cells and shared molecular markers and mRNA expression profiles with human ccRCC. Exome sequencing revealed that mouse and human ccRCCs exhibit recurrent mutations in genes associated with the primary cilium, uncovering a mutational convergence on this organelle and implicating a subset of ccRCCs as genetic ciliopathies. Different mouse tumors responded differently to standard therapies for advanced human ccRCC, mimicking the range of clinical behaviors in the human disease. Inhibition of hypoxia-inducible factor (HIF)-α transcription factors with acriflavine as third-line therapy had therapeutic effects in some tumors, providing preclinical evidence for further investigation of HIF-α inhibition as a ccRCC treatment. This autochthonous mouse ccRCC model represents a tool to investigate the biology of ccRCC and to identify new treatment strategies.