Marlous Hoogstraat

Chromothripsis is a common mechanism driving genomic rearrangements in primary and metastatic colorectal cancer

BackgroundStructural rearrangements form a major class of somatic variation in cancer genomes. Local chromosome shattering, termed chromothripsis, is a mechanism proposed to be the cause of clustered chromosomal rearrangements and was recently described to occur in a small percentage of tumors. The significance of these clusters for tumor development or metastatic spread is largely unclear.ResultsWe used genome-wide long mate-pair sequencing and SNP array profiling to reveal that chromothripsis is a widespread phenomenon in primary colorectal cancer and metastases. We find large and small chromothripsis events in nearly every colorectal tumor sample and show that several breakpoints of chromothripsis clusters and isolated rearrangements affect cancer genes, including NOTCH2, EXO1 and MLL3. We complemented the structural variation studies by sequencing the coding regions of a cancer exome in all colorectal tumor samples and found somatic mutations in 24 genes, including APC, KRAS, SMAD4 and PIK3CA. A pairwise comparison of somatic variations in primary and metastatic samples indicated that many chromothripsis clusters, isolated rearrangements and point mutations are exclusively present in either the primary tumor or the metastasis and may affect cancer genes in a lesion-specific manner.ConclusionsWe conclude that chromothripsis is a prevalent mechanism driving structural rearrangements in colorectal cancer and show that a complex interplay between point mutations, simple copy number changes and chromothripsis events drive colorectal tumor development and metastasis.

Ovarian Cancer Cell Line Panel (OCCP): Clinical Importance of In Vitro Morphological Subtypes

Epithelial ovarian cancer is a highly heterogeneous disease and remains the most lethal gynaecological malignancy in the Western world. Therapeutic approaches need to account for inter-patient and intra-tumoural heterogeneity and detailed characterization of in vitro models representing the different histological and molecular ovarian cancer subtypes is critical to enable reliable preclinical testing. There are approximately 100 publicly available ovarian cancer cell lines but their cellular and molecular characteristics are largely undescribed. We have characterized 39 ovarian cancer cell lines under uniform conditions for growth characteristics, mRNA/microRNA expression, exon sequencing, drug response for clinically-relevant therapeutics and collated all available information on the original clinical features and site of origin. We tested for statistical associations between the cellular and molecular features of the lines and clinical features. Of the 39 ovarian cancer cell lines, 14 were assigned as high-grade serous, four serous-type, one low-grade serous and 20 non-serous type. Three morphological subtypes: Epithelial (n = 21), Round (n = 7) and Spindle (n = 12) were identified that showed distinct biological and molecular characteristics, including overexpression of cell movement and migration-associated genes in the Spindle subtype. Comparison with the original clinical data showed association of the spindle-like tumours with metastasis, advanced stage, suboptimal debulking and poor prognosis. In addition, the expression profiles of Spindle, Round and Epithelial morphologies clustered with the previously described C1-stromal, C5-mesenchymal and C4 ovarian subtype expression profiles respectively. Comprehensive profiling of 39 ovarian cancer cell lines under controlled, uniform conditions demonstrates clinically relevant cellular and genomic characteristics. This data provides a rational basis for selecting models to develop specific treatment approaches for histological and molecular subtypes of ovarian cancer.

Targeted next-generation sequencing: A novel diagnostic tool for primary immunodeficiencies

BACKGROUND Primary immunodeficiency (PID) disorders are a heterogeneous group of inherited disorders caused by a variety of monogenetic immune defects. Thus far, mutations in more than 170 different genes causing PIDs have been described. A clear genotype-phenotype correlation is often not available, which makes a genetic diagnosis in patients with PIDs complex and laborious. OBJECTIVE We sought to develop a robust, time-effective, and cost-effective diagnostic method to facilitate a genetic diagnosis in any of 170 known PID-related genes by using next-generation sequencing (NGS). METHODS We used both targeted array-based and in-solution enrichment combined with a SOLiD sequencing platform and a bioinformatic pipeline developed in house to analyze genetic changes in the DNA of 41 patients with PIDs with known mutations and 26 patients with undiagnosed PIDs. RESULTS This novel NGS-based method accurately detected point mutations (sensitivity and specificity >99% in covered regions) and exonic deletions (100% sensitivity and specificity). For the 170 genes of interest, the DNA coverage was greater than 20× in 90% to 95%. Nine PID-related genes proved not eligible for evaluation by using this NGS-based method because of inadequate coverage. The NGS method allowed us to make a genetic diagnosis in 4 of 26 patients who lacked a genetic diagnosis despite routine functional and genetic testing. Three of these patients proved to have an atypical presentation of previously described PIDs. CONCLUSION This novel NGS tool facilitates accurate simultaneous detection of mutations in 161 of 170 known PID-related genes. In addition, these analyses will generate more insight into genotype-phenotype correlations for the different PID disorders.

Genomic and transcriptomic plasticity in treatment-naïve ovarian cancer

Intra-tumor heterogeneity is a hallmark of many cancers and may lead to therapy resistance or interfere with personalized treatment strategies. Here, we combined topographic mapping of somatic breakpoints and transcriptional profiling to probe intra-tumor heterogeneity of treatment-naïve stage IIIC/IV epithelial ovarian cancer. We observed that most substantial differences in genomic rearrangement landscapes occurred between metastases in the omentum and peritoneum versus tumor sites in the ovaries. Several cancer genes such as NF1, CDKN2A, and FANCD2 were affected by lesion-specific breakpoints. Furthermore, the intra-tumor variability involved different mutational hallmarks including lesion-specific kataegis (local mutation shower coinciding with genomic breakpoints), rearrangement classes, and coding mutations. In one extreme case, we identified two independent TP53 mutations in ovary tumors and omentum/peritoneum metastases, respectively. Examination of gene expression dynamics revealed up-regulation of key cancer pathways including WNT, integrin, chemokine, and Hedgehog signaling in only subsets of tumor samples from the same patient. Finally, we took advantage of the multilevel tumor analysis to understand the effects of genomic breakpoints on qualitative and quantitative gene expression changes. We show that intra-tumor gene expression differences are caused by site-specific genomic alterations, including formation of in-frame fusion genes. These data highlight the plasticity of ovarian cancer genomes, which may contribute to their strong capacity to adapt to changing environmental conditions and give rise to the high rate of recurrent disease following standard treatment regimes.

Simultaneous Detection of Clinically Relevant Mutations and Amplifications for Routine Cancer Pathology

In routine cancer molecular pathology, various independent experiments are required to determine mutation and amplification status of clinically relevant genes. Most of these tests are designed to identify a limited number of genetic aberrations, most likely in a given tumor type. We present a modified version of a multiplexed PCR and IonTorrent-based sequencing approach that can replace a large number of existing assays. The test allows for the simultaneous detection of point mutations and gene amplifications in 40 genes, including known hotspot regions in oncogenes (KRAS, BRAF), inactivating mutations in tumor suppressors (TP53, PTEN), and oncogene amplifications (ERBB2, EGFR). All point mutations were confirmed using certified diagnostic assays, and a sensitivity and specificity of 100% (95% CI, 0.875-1.0) and 99% (95% CI, 0.960-0.999), respectively, were determined for amplifications in FFPE material. Implementation of a single assay to effectively detect mutations and amplifications in clinically relevant genes not only improves the efficiency of the workflow within diagnostic laboratories but also increases the chance of detecting (rare) actionable variants for a given tumor type that are typically missed in routine pathology. The ability to obtain comprehensive and rapid mutational overviews is key for improving the efficiency of cancer patient care through tailoring treatments based on the genetic characteristics of individual tumors.

Detailed imaging and genetic analysis reveal a secondary BRAF(L505H) resistance mutation and extensive intrapatient heterogeneity in metastatic BRAF mutant melanoma patients treated with vemurafenib

Resistance to treatment is the main problem of targeted treatment for cancer. We followed ten patients during treatment with vemurafenib, by three‐dimensional imaging. In all patients, only a subset of lesions progressed. Next‐generation DNA sequencing was performed on sequential biopsies in four patients to uncover mechanisms of resistance. In two patients, we identified mutations that explained resistance to vemurafenib; one of these patients had a secondary BRAF L505H mutation. This is the first observation of a secondary BRAF mutation in a vemurafenib‐resistant patient‐derived melanoma sample, which confirms the potential importance of the BRAF L505H mutation in the development of therapy resistance. Moreover, this study hints toward an important role for tumor heterogeneity in determining the outcome of targeted treatments.

Predicting clinical benefit from everolimus in patients with advanced solid tumors, the CPCT-03 study

BACKGROUND In this study, our aim was to identify molecular aberrations predictive for response to everolimus, an mTOR inhibitor, regardless of tumor type. METHODS To generate hypotheses about potential markers for sensitivity to mTOR inhibition, drug sensitivity and genomic profiles of 835 cell lines were analyzed. Subsequently, a multicenter study was conducted. Patients with advanced solid tumors lacking standard of care treatment options were included and underwent a pre-treatment tumor biopsy to enable DNA sequencing of 1,977 genes, derive copy number profiles and determine activation status of pS6 and pERK. Treatment benefit was determined according to TTP ratio and RECIST. We tested for associations between treatment benefit and single molecular aberrations, clusters of aberrations and pathway perturbation. RESULTS Cell line screens indicated several genes, such as PTEN (P = 0.016; Wald test), to be associated with sensitivity to mTOR inhibition. Subsequently 73 patients were included, of which 59 started treatment with everolimus. Response and molecular data were available from 43 patients. PTEN aberrations, i.e. copy number loss or mutation, were associated with treatment benefit (P = 0.046; Fishers exact test). CONCLUSION Loss-of-function aberrations in PTEN potentially represent a tumor type agnostic biomarker for benefit from everolimus and warrants further confirmation in subsequent studies.Background In this study, our aim was to identify molecular aberrations predictive for response to everolimus, an mTOR inhibitor, regardless of tumor type. Methods To generate hypotheses about potential markers for sensitivity to mTOR inhibition, drug sensitivity and genomic profiles of 835 cell lines were analyzed. Subsequently, a multicenter study was conducted. Patients with advanced solid tumors lacking standard of care treatment options were included and underwent a pre-treatment tumor biopsy to enable DNA sequencing of 1,977 genes, derive copy number profiles and determine activation status of pS6 and pERK. Treatment benefit was determined according to TTP ratio and RECIST. We tested for associations between treatment benefit and single molecular aberrations, clusters of aberrations and pathway perturbation. Results Cell line screens indicated several genes, such as PTEN (P = 0.016; Wald test), to be associated with sensitivity to mTOR inhibition. Subsequently 73 patients were included, of which 59 started treatment with everolimus. Response and molecular data were available from 43 patients. PTEN aberrations, i.e. copy number loss or mutation, were associated with treatment benefit (P = 0.046; Fishers exact test). Conclusion Loss-of-function aberrations in PTEN potentially represent a tumor type agnostic biomarker for benefit from everolimus and warrants further confirmation in subsequent studies.

TP53 mutated glioblastoma stem-like cell cultures are sensitive to dual mTORC1/2 inhibition while resistance in TP53 wild type cultures can be overcome by combined inhibition of mTORC1/2 and Bcl-2

Background Glioblastoma is the most malignant tumor of the central nervous system and still lacks effective treatment. This study explores mutational biomarkers of 11 drugs targeting either the RTK/Ras/PI3K, the p53 or the Rb pathway using 25 patient-derived glioblastoma stem-like cell cultures (GSCs). Results We found that TP53 mutated GSCs were approximately 3.5 fold more sensitive to dual inhibition of mammalian target of rapamycin complex 1 and 2 (mTORC1/2) compared to wild type GSCs. We identified that Bcl-2(Thr56/Ser70) phosphorylation contributed to the resistance of TP53 wild type GSCs against dual mTORC1/2 inhibition. The Bcl-2 inhibitor ABT-263 (navitoclax) increased sensitivity to the mTORC1/2 inhibitor AZD8055 in TP53 wild type GSCs, while sensitivity to AZD8055 in TP53 mutated GSCs remained unchanged. Conclusion Our data suggest that Bcl-2 confers resistance to mTORC1/2 inhibitors in TP53 wild type GSCs and that combined inhibition of both mTORC1/2 and Bcl-2 is worthwhile to explore further in TP53 wild type glioblastomas, whereas in TP53 mutated glioblastomas dual mTORC1/2 inhibitors should be explored.

RUBIC identifies driver genes by detecting recurrent DNA copy number breaks

The frequent recurrence of copy number aberrations across tumour samples is a reliable hallmark of certain cancer driver genes. However, state-of-the-art algorithms for detecting recurrent aberrations fail to detect several known drivers. In this study, we propose RUBIC, an approach that detects recurrent copy number breaks, rather than recurrently amplified or deleted regions. This change of perspective allows for a simplified approach as recursive peak splitting procedures and repeated re-estimation of the background model are avoided. Furthermore, we control the false discovery rate on the level of called regions, rather than at the probe level, as in competing algorithms. We benchmark RUBIC against GISTIC2 (a state-of-the-art approach) and RAIG (a recently proposed approach) on simulated copy number data and on three SNP6 and NGS copy number data sets from TCGA. We show that RUBIC calls more focal recurrent regions and identifies a much larger fraction of known cancer genes.

GPCR-SAS: A web application for statistical analyses on G protein-coupled receptors sequences

G protein-coupled receptors (GPCRs) are one of the largest protein families in mammals. They mediate signal transduction across cell membranes and are important targets for the pharmaceutical industry. The G Protein-Coupled Receptors—Sequence Analysis and Statistics (GPCR-SAS) web application provides a set of tools to perform comparative analysis of sequence positions between receptors, based on a curated structural-informed multiple sequence alignment. The analysis tools include: (i) percentage of occurrence of an amino acid or motif and entropy at a position or range of positions, (ii) covariance of two positions, (iii) correlation between two amino acids in two positions (or two sequence motifs in two ranges of positions), and (iv) snake-plot representation for a specific receptor or for the consensus sequence of a group of selected receptors. The analysis of conservation of residues and motifs across transmembrane (TM) segments may guide the design of more selective ligands or help to rationalize activation mechanisms, among others. As an example, here we analyze the amino acids of the “transmission switch”, that initiates receptor activation following ligand binding. The tool is freely accessible at http://lmc.uab.cat/gpcrsas/.