Rapid identification of somatic mutations in colorectal and breast cancer tissues using mismatch repair detection (MRD)

Chromosomal copy number and mutational status are required to authenticate ovarian cancer cell lines as appropriate cell models

BACKGROUND

The mutational status of ovarian cancer cell line IGROV-1 is inconsistent across the literature, suggestive of multiple clonal populations of the cell line. IGROV-1 has previously been categorised as an inappropriate model for high-grade serous ovarian cancer.

METHODS

IGROV-1 cells were obtained from the Netherlands Cancer Institute (IGROV-1-NKI) and the MD Anderson Cancer Centre (IGROV-1-MDA). Cell lines were STR fingerprinted and had their chromosomal copy number analysed and BRCA1/2 genes sequenced. Mutation status of ovarian cancer-related genes were extracted from the literature.

RESULTS

The IGROV-1-NKI cell line has a tetraploid chromosomal profile. In contrast, the IGROV-1-MDA cell line has pseudo-normal chromosomes. The IGROV-1-NKI and IGROV-MDA are both STR matches (80.7% and 84.6%) to the original IGROV-1 cells isolated in 1985. However, IGROV-1-NKI and IGROV-1-MDA are not an STR match to each other (78.1%) indicating genetic drift. The BRCA1 and BRCA2 gene sequences are 100% identical between IGROV-1-MDA and IGROV-1-NKI, including a BRCA1 heterozygous deleterious mutation. The IGROV-1-MDA cells are more resistant to cisplatin and olaparib than IGROV-1-NKI. IGROV-1 has a mutational profile consistent with both Type I (PTEN, PIK3CA and ARID1A) and Type II ovarian cancer (BRCA1, TP53) and is likely to be a Type II high-grade serous carcinoma of the SET (Solid, pseudo-Endometroid and Transitional cell carcinoma-like morphology) subtype.

CONCLUSIONS

Routine testing of chromosomal copy number as well as the mutational status of ovarian cancer related genes should become the new standard alongside STR fingerprinting to ensure that ovarian cancer cell lines are appropriate models.

Stochastic phenotypes in RAS-dependent developmental diseases

Germline mutations upregulating RAS signaling are associated with multiple developmental disorders. A hallmark of these conditions is that the same mutation may present vastly different phenotypes in different individuals, even in monozygotic twins. Here, we demonstrate how the origins of such largely unexplained phenotypic variations may be dissected using highly controlled studies in Drosophila that have been gene edited to carry activating variants of MEK, a core enzyme in the RAS pathway. This allowed us to measure the small but consistent increase in signaling output of such alleles in vivo. The fraction of mutation carriers reaching adulthood was strongly reduced, but most surviving animals had normal RAS-dependent structures. We rationalize these results using a stochastic signaling model and support it by quantifying cell fate specification errors in bilaterally symmetric larval trachea, a RAS-dependent structure that allows us to isolate the effects of mutations from potential contributions of genetic modifiers and environmental differences. We propose that the small increase in signaling output shifts the distribution of phenotypes into a regime, where stochastic variation causes defects in some individuals, but not in others. Our findings shed light on phenotypic heterogeneity of developmental disorders caused by deregulated RAS signaling and offer a framework for investigating causal effects of other pathogenic alleles and mild mutations in general.

Computational screening and structural analysis of Gly201Arg and Gly201Asp missense mutations in human cyclin-dependent kinase 4 protein

The regulatory proteins, cyclins, and cyclin-dependent kinases (CDKs) control the cell cycle progression. CDK4 gene mutations are associated with certain cancers such as melanoma, breast cancer, and rhabdomyosarcoma. Therefore, understanding the mechanisms of cell cycle control and cell proliferation is essential in developing cancer treatment regimens. In this study, we obtained cancer-causing CDK4 mutations from the COSMIC database and subjected them to a series of in silico analyses to identify the most significant mutations. An overall of 238 mutations (119 missense mutations) retrieved from the COSMIC database were investigated for the pathogenic and destabilizing properties using the PredictSNP and iStable algorithms. Further, the amino acid position of the most pathogenic and destabilizing mutations were analyzed to understand the nature of amino acid conservation across the species during the evolution. We observed that the missense mutations G201R and G201D were more significant and the Glycine at position 201 was found to highly conserved. These significant mutations were subjected to molecular dynamics simulation analysis to understand the protein's structural changes. The results from molecular dynamics simulations revealed that both G201R and G201D of CDK4 are capable of altering the protein's native form. On comparison among the most significant mutations, G201R disrupted the protein structure higher than the protein with G201D.

A guided-inquiry investigation of genetic variants using Oxford nanopore sequencing for an undergraduate molecular biology laboratory course

Next Generation Sequencing (NGS) has become an important tool in the biological sciences and has a growing number of applications across medical fields. Currently, few undergraduate programs provide training in the design and implementation of NGS applications. Here, we describe an inquiry-based laboratory exercise for a college-level molecular biology laboratory course that uses real-time MinION deep sequencing and bioinformatics to investigate characteristic genetic variants found in cancer cell-lines. The overall goal for students was to identify non-small cell lung cancer (NSCLC) cell-lines based on their unique genomic profiles. The units described in this laboratory highlight core principles in multiplex PCR primer design, real-time deep sequencing, and bioinformatics analysis for genetic variants. We found that the MinION device is an appropriate, feasible tool that provides a comprehensive, hands-on NGS experience for undergraduates. Student evaluations demonstrated increased confidence in using molecular techniques and enhanced understanding of NGS concepts. Overall, this exercise provides a pedagogical tool for incorporating NGS approaches in the teaching laboratory as way of enhancing students' comprehension of genomic sequence analysis. Further, this NGS lab module can easily be added to a variety of lab-based courses to help undergraduate students learn current DNA sequencing methods with limited effort and cost.

Molecular mechanisms underlying cellular effects of human MEK1 mutations

Terminal regions of Drosophila embryos are patterned by signaling through ERK, which is genetically deregulated in multiple human diseases. Quantitative studies of terminal patterning have been recently used to investigate gain-of-function variants of human MEK1, encoding the MEK kinase that directly activates ERK by dual phosphorylation. Unexpectedly, several mutations reduced ERK activation by extracellular signals, possibly through a negative feedback triggered by signal-independent activity of the mutant variants. Here we present experimental evidence supporting this model. Using a MEK variant that combines a mutation within the negative regulatory region with alanine substitutions in the activation loop, we prove that pathogenic variants indeed acquire signal-independent kinase activity. We also demonstrate that signal-dependent activation of these variants is independent of kinase suppressor of Ras, a conserved adaptor that is indispensable for activation of normal MEK. Finally, we show that attenuation of ERK activation by extracellular signals stems from transcriptional induction of Mkp3, a dual specificity phosphatase that deactivates ERK by dephosphorylation. These findings in the Drosophila embryo highlight its power for investigating diverse effects of human disease mutations.

Inference of Multisite Phosphorylation Rate Constants and Their Modulation by Pathogenic Mutations

Multisite protein phosphorylation plays a critical role in cell regulation [1-3]. It is widely appreciated that the functional capabilities of multisite phosphorylation depend on the order and kinetics of phosphorylation steps, but kinetic aspects of multisite phosphorylation remain poorly understood [4-6]. Here, we focus on what appears to be the simplest scenario, when a protein is phosphorylated on only two sites in a strict, well-defined order. This scenario describes the activation of ERK, a highly conserved cell-signaling enzyme. We use Bayesian parameter inference in a structurally identifiable kinetic model to dissect dual phosphorylation of ERK by MEK, a kinase that is mutated in a large number of human diseases [7-12]. Our results reveal how enzyme processivity and efficiencies of individual phosphorylation steps are altered by pathogenic mutations. The presented approach, which connects specific mutations to kinetic parameters of multisite phosphorylation mechanisms, provides a systematic framework for closing the gap between studies with purified enzymes and their effects in the living organism.

Phospho-proteomic analyses of B-Raf protein complexes reveal new regulatory principles

B-Raf represents a critical physiological regulator of the Ras/RAF/MEK/ERK-pathway and a pharmacological target of growing clinical relevance, in particular in oncology. To understand how B-Raf itself is regulated, we combined mass spectrometry with genetic approaches to map its interactome in MCF-10A cells as well as in B-Raf deficient murine embryonic fibroblasts (MEFs) and B-Raf/Raf-1 double deficient DT40 lymphoma cells complemented with wildtype or mutant B-Raf expression vectors. Using a multi-protease digestion approach, we identified a novel ubiquitination site and provide a detailed B-Raf phospho-map. Importantly, we identify two evolutionary conserved phosphorylation clusters around T401 and S419 in the B-Raf hinge region. SILAC labelling and genetic/biochemical follow-up revealed that these clusters are phosphorylated in the contexts of oncogenic Ras, sorafenib induced Raf dimerization and in the background of the V600E mutation. We further show that the vemurafenib sensitive phosphorylation of the T401 cluster occurs in trans within a Raf dimer. Substitution of the Ser/Thr-residues of this cluster by alanine residues enhances the transforming potential of B-Raf, indicating that these phosphorylation sites suppress its signaling output. Moreover, several B-Raf phosphorylation sites, including T401 and S419, are somatically mutated in tumors, further illustrating the importance of phosphorylation for the regulation of this kinase.

Afatinib Activity in Platinum-Refractory Metastatic Urothelial Carcinoma in Patients With ERBB Alterations

PURPOSE

Somatic mutations and copy number variation in the ERBB family are frequent in urothelial carcinoma (UC) and may represent viable therapeutic targets. We studied whether afatinib (an oral, irreversible inhibitor of the ErbB family) has activity in UC and if specific ERBB molecular alterations are associated with clinical response.

PATIENTS AND METHODS

In this phase II trial, patients with metastatic platinum-refractory UC received afatinib 40 mg/day continuously until progression or intolerance. The primary end point was 3-month progression-free survival (PFS3). Prespecified tumor analysis for alterations in EGFR, HER2, ERBB3, and ERBB4 was conducted.

RESULTS

The first-stage enrollment goal of 23 patients was met. Patient demographic data included: 78% male, median age 67 years (range, 36 to 82 years), hemoglobin < 10 g/dL in 17%, liver metastases in 30%, median time from prior chemotherapy of 3.6 months, and Eastern Cooperative Oncology Group performance status ≤ 1 in 100%. No unexpected toxicities were observed; two patients required dose reduction for grade 3 fatigue and rash. Overall, five of 23 patients (21.7%) met PFS3 (two partial response, three stable disease). Notably, among the 21 tumors analyzed, five of six patients (83.3%) with HER2 and/or ERBB3 alterations achieved PFS3 (PFS = 10.3, 7.0, 6.9, 6.3, and 5.0 months, respectively) versus none of 15 patients without alterations (P < .001). Three of four patients with HER2 amplification and three of three patients with ERBB3 somatic mutations (G284R, V104M, and R103G) met PFS3. One patient with both HER2 amplification and ERBB3 mutation never progressed on therapy, but treatment was discontinued after 10.3 months as a result of depressed ejection fraction. The median time to progression/discontinuation was 6.6 months in patients with HER2/ERBB3 alterations versus 1.4 months in patients without alterations (P < .001).

CONCLUSION

Afatinib demonstrated significant activity in patients with platinum-refractory UC with HER2 or ERBB3 alterations. The potential contribution of ERBB3 to afatinib sensitivity is novel. Afatinib deserves further investigation in molecularly selected UC.

EGFR and HER2 signals play a salvage role in MEK1-mutated gastric cancer after MEK inhibition

Since the prognosis of unresectable advanced gastric cancer remains poor, novel therapeutic strategies are needed. Somatic MEK1 gene mutations have been reported as oncogenic activating mutations in gastric cancer, and MEK inhibitors can be effective against such gastric cancers. In the present study, however, activated EGFR and HER2 signals after treatment with a MEK inhibitor (trametinib) were found in a MEK1-mutated gastric cancer cell line (OCUM-1 cell line) using a phospho-receptor tyrosine kinase array. The phosphorylation of EGFR and HER2 reactivated ERK1/2, which had been inhibited by trametinib, and EGF stimulation led to resistance to trametinib in this cell line. Lapatinib, an EGFR and an HER2 inhibitor, reversed the activation of ERK1/2 by inhibiting the phosphorylation of EGFR and HER2 and cancelled the resistance. The combination of trametinib and lapatinib synergistically inhibited the cell growth of the OCUM-1 cell line and strongly induced apoptosis by inhibiting the activated EGFR and HER2 signals. These results suggest that the EGFR and HER2 signals play a salvage role and are related to resistance to MEK inhibitors in MEK1‑mutated gastric cancer. Moreover, combination therapy with trametinib and lapatinib can exhibit a synergistic effect and may contribute to overcoming the resistance to MEK inhibitors.

Portrait of the PI3K/AKT pathway in colorectal cancer

PI3K/AKT signaling leads to reduced apoptosis, stimulates cell growth and increases proliferation. Under normal conditions, PI3K/AKT activation is tightly controlled and dependent on both extracellular growth signals and the availability of amino acids and glucose. Genetic aberrations leading to PI3K/AKT hyper-activation are observed at considerable frequency in all major nodes in most tumors. In colorectal cancer the most commonly observed pathway changes are IGF2 overexpression, PIK3CA mutations and PTEN mutations and deletions. Combined, these alterations are found in about 40% of large bowel tumors. In addition, but not mutually exclusive to these, KRAS mutations are observed at a similar frequency. There are however additional, less frequent and more poorly understood events that may also push the PI3K/AKT pathway into overdrive and thus promote malignant growth. Here we discuss aberrations of components at the genetic, epigenetic, transcriptional, post-transcriptional, translational and post-translational level where perturbations may drive excessive PI3K/AKT signaling. Integrating multiple molecular levels will advance our understanding of this cancer critical circuit and more importantly, improve our ability to pharmacologically target the pathway in view of clonal development, tumor heterogeneity and drug resistance mechanisms. In this review, we revisit the PI3K/AKT pathway cancer susceptibility syndromes, summarize the known aberrations at the different regulatory levels and the prognostic and predictive values of these alterations in colorectal cancer.

Proliferation of colorectal cancer is promoted by two signaling transduction expression patterns: ErbB2/ErbB3/AKT and MET/ErbB3/MAPK

One of the recent breakthroughs in cancer research is the identification of activating mutations in various receptor tyrosine kinase(RTK) pathways in many cancers including colorectal cancer(CRC). We hypothesize that, alternative to mutations, overexpression of various oncogenic RTKs may also underpin CRC pathogenesis, and different RTK may couple with distinct downstream signaling pathways in different subtypes of human CRC. By immunohistochemistry, we show here that RTK members ErbB2, ErbB3 and c-Met were in deed differentially overexpressed in colorectal cancer patient samples leading to constitutive activation of RTK signaling pathways. Using ErbB2 specific inhibitor Lapatinib and c-Met specific inhibitor PHA-665752, we further demonstrated that this constitutive activation of RTK signaling is necessary for the survival of colorectal cancer cells. Furthermore, we show that RTK overexpression pattern dictates the use of downstream AKT and/or MAPK pathways. Our data are important additions to current oncogenic mutation models, and further explain the clinical variation in therapeutic responses of colorectal cancer. Our findings advocate for more personalized therapy tailored to individual patients based on their type of RTK expression in addition to their mutation status.

Mutation spectrum in human colorectal cancers and potential functional relevance

Background

Somatic variants, which occur in the genome of all cells, are well accepted to play a critical role in cancer development, as their accumulation in genes could affect cell proliferations and cell cycle.

Methods

In order to understand the role of somatic mutations in human colorectal cancers, we characterized the mutation spectrum in two colorectal tumor tissues and their matched normal tissues, by analyzing deep-sequenced transcriptome data.

Results

We found a higher mutation rate of somatic variants in tumor tissues in comparison with normal tissues, but no trend was observed for mutation properties. By applying a series of stringent filters, we identified 418 genes with tumor specific disruptive somatic variants. Of these genes, three genes in mucin protein family (MUC2, MUC4, and MU12) are of particular interests. It has been reported that the expression of mucin proteins was correlated with the progression of colorectal cancer therefore somatic variants within those genes can interrupt their normal expression and thus contribute to the tumorigenesis.

Conclusions

Our findings provide evidence of the utility of RNA-Seq in mutation screening in cancer studies, and suggest a list of candidate genes for future colorectal cancer diagnosis and treatment.

MEK genomics in development and disease

The mitogen-activated protein kinase kinases (the MAPK/ERK kinases; MKKs or MEKs) and their downstream substrates, the extracellular-regulated kinases have been intensively studied for their roles in development and disease. Until recently, it had been assumed any mutation affecting their function would have lethal consequences. However, the identification of MEK1 and MEK2 mutations in developmental syndromes as well as chemotherapy-resistant tumors, and the discovery of genomic variants in MEK1 and MEK2 have led to the realization the extent of genomic variation associated with MEKs is much greater than had been appreciated. In this review, we will discuss these recent advances, relating them to what is currently understood about the structure and function of MEKs, and describe how they change our understanding of the role of MEKs in development and disease.

Oncogenic MAP2K1 mutations in human epithelial tumors

The scirrhous subtype of gastric cancer is a highly infiltrative tumor with a poor outcome. To identify a transforming gene in this intractable disorder, we constructed a retroviral complementary DNA (cDNA) expression library from a cell line (OCUM-1) of scirrhous gastric cancer. A focus formation assay with the library and mouse 3T3 fibroblasts led to the discovery of a transforming cDNA, encoding for MAP2K1 with a glutamine-to-proline substitution at amino acid position 56. Interestingly, treatment with a MAP2K1-specific inhibitor clearly induced cell death of OCUM-1 but not of other two cells lines of scirrhous gastric cancer that do not carry MAP2K1 mutations, revealing the essential role of MAP2K1(Q56P) in the transformation mechanism of OCUM-1 cells. By using a next-generation sequencer, we further conducted deep sequencing of the MAP2K1 cDNA among 171 human cancer specimens or cell lines, resulting in the identification of one known (D67N) and four novel (R47Q, R49L, I204T and P306H) mutations within MAP2K1. The latter four changes were further shown to confer transforming potential to MAP2K1. In our experiments, a total of six (3.5%) activating mutations in MAP2K1 were thus identified among 172 of specimens or cell lines for human epithelial tumors. Given the addiction of cancer cells to the elevated MAP2K1 activity for proliferation, human cancers with such MAP2K1 mutations are suitable targets for the treatment with MAP2K1 inhibitors.

Transactional database transformation and its application in prioritizing human disease genes

Binary (0,1) matrices, commonly known as transactional databases, can represent many application data, including genephenotype data where “1” represents a confirmed gene-phenotype relation and “0” represents an unknown relation. It is natural to ask what information is hidden behind these “0”s and “1”s. Unfortunately, recent matrix completion methods, though very effective in many cases, are less likely to infer something interesting from these (0,1)-matrices. To answer this challenge, we propose INDEVI, a very succinct and effective algorithm to perform independent-evidence-based transactional database transformation. Each entry of a (0,1)-matrix is evaluated by “independent evidence” (maximal supporting patterns) extracted from the whole matrix for this entry. The value of an entry, regardless of its value as 0 or 1, has completely no effect for its independent evidence. The experiment on a genephenotype database shows that our method is highly promising in ranking candidate genes and predicting unknown disease genes.

Exome sequencing identifies recurrent somatic MAP2K1 and MAP2K2 mutations in melanoma

We performed exome sequencing to detect somatic mutations in protein-coding regions in seven melanoma cell lines and donor-matched germline cells. All melanoma samples had high numbers of somatic mutations, which showed the hallmark of UV-induced DNA repair. Such a hallmark was absent in tumor sample-specific mutations in two metastases derived from the same individual. Two melanomas with non-canonical BRAF mutations harbored gain-of-function MAP2K1 and MAP2K2 (MEK1 and MEK2, respectively) mutations, resulting in constitutive ERK phosphorylation and higher resistance to MEK inhibitors. Screening a larger cohort of individuals with melanoma revealed the presence of recurring somatic MAP2K1 and MAP2K2 mutations, which occurred at an overall frequency of 8%. Furthermore, missense and nonsense somatic mutations were frequently found in three candidate melanoma genes, FAT4, LRP1B and DSC1.

BRAF p.Val600Glu (V600E) Testing for Assessment of Treatment Options in Metastatic Colorectal Cancer

Colon and rectal cancer (CRC) are the third most common cancer in the United States and cause approximately 50,000 deaths per year. The anti-epidermal growth factor receptor (EGFR) monoclonal antibodies cetuximab (Erbitux®) and panitumumab (Vectibix®) have been recently introduced to treat CRC. However, the response rate with these agents is low and they are associated with serious adverse effects. Accordingly biomarkers that can predict those patients that will respond to treatment may have clinical utility. The p.Val600Glu sequence variant (often called V600E) in the BRAF gene has been investigated as a biomarker to predict patients that will not respond to treatment with the anti-EGFR monoclonal antibodies.

Clinical relevance of KRAS in human cancers

The KRAS gene (Ki-ras2 Kirsten rat sarcoma viral oncogene homolog) is an oncogene that encodes a small GTPase transductor protein called KRAS. KRAS is involved in the regulation of cell division as a result of its ability to relay external signals to the cell nucleus. Activating mutations in the KRAS gene impair the ability of the KRAS protein to switch between active and inactive states, leading to cell transformation and increased resistance to chemotherapy and biological therapies targeting epidermal growth factor receptors. This review highlights some of the features of the KRAS gene and the KRAS protein and summarizes current knowledge of the mechanism of KRAS gene regulation. It also underlines the importance of activating mutations in the KRAS gene in relation to carcinogenesis and their importance as diagnostic biomarkers, providing clues regarding human cancer patients' prognosis and indicating potential therapeutic approaches.

High-throughput, high-accuracy array-based resequencing

Although genomewide association studies have successfully identified associations of many common single-nucleotide polymorphisms (SNPs) with common diseases, the SNPs implicated so far account for only a small proportion of the genetic variability of tested diseases. It has been suggested that common diseases may often be caused by rare alleles missed by genomewide association studies. To identify these rare alleles we need high-throughput, high-accuracy resequencing technologies. Although array-based genotyping has allowed genomewide association studies of common SNPs in tens of thousands of samples, array-based resequencing has been limited for 2 main reasons: the lack of a fully multiplexed pipeline for high-throughput sample processing, and failure to achieve sufficient performance. We have recently solved both of these problems and created a fully multiplexed high-throughput pipeline that results in high-quality data. The pipeline consists of target amplification from genomic DNA, followed by allele enrichment to generate pools of purified variant (or nonvariant) DNA and ends with interrogation of purified DNA on resequencing arrays. We have used this pipeline to resequence approximately 5 Mb of DNA (on 3 arrays) corresponding to the exons of 1,500 genes in >473 samples; in total >2,350 Mb were sequenced. In the context of this large-scale study we obtained a false positive rate of approximately 1 in 500,000 bp and a false negative rate of approximately 10%.

A comprehensive assay for targeted multiplex amplification of human DNA sequences

We developed a robust and reproducible methodology to amplify human sequences in parallel for use in downstream multiplexed sequence analyses. We call the methodology SMART (Spacer Multiplex Amplification Reaction), and it is based, in part, on padlock probe technology. As a proof of principle, we used SMART technology to simultaneously amplify 485 human exons ranging from 100 to 500 bp from human genomic DNA. In multiple repetitions, >90% of the targets were successfully amplified with a high degree of uniformity, with 70% of targets falling within a 10-fold range and all products falling within a 100-fold range of each other in abundance. We used long padlock probes (LPPs) >300 bases in length for the assay, and the increased length of these probes allowed for the capture of human sequences up to 500 bp in length, which is optimal for capturing most human exons. To engineer the LPPs, we developed a method that generates ssDNA molecules with precise ends, using an appropriately designed dsDNA template. The template has appropriate restriction sites engineered into it that can be digested to generate nucleotide overhangs that are suitable for lambda exonuclease digestion, producing a single-stranded probe from dsDNA. The SMART technology is flexible and can be easily adapted to multiplex tens of thousands of target sequences in a single reaction.