Gastric cancer is a leading cause of cancer deaths, but analysis of its molecular and clinical characteristics has been complicated by histological and aetiological heterogeneity. Here we describe a comprehensive molecular evaluation of 295 primary gastric adenocarcinomas as part of The Cancer Genome Atlas (TCGA) project. We propose a molecular classification dividing gastric cancer into four subtypes: tumours positive for Epstein–Barr virus, which display recurrent PIK3CA mutations, extreme DNA hypermethylation, and amplification of JAK2, CD274 (also known as PD-L1) and PDCD1LG2 (also knownasPD-L2); microsatellite unstable tumours, which show elevated mutation rates, including mutations of genes encoding targetable oncogenic signalling proteins; genomically stable tumours, which are enriched for the diffuse histological variant and mutations of RHOA or fusions involving RHO-family GTPase-activating proteins; and tumours with chromosomal instability, which show marked aneuploidy and focal amplification of receptor tyrosine kinases. Identification of these subtypes provides a roadmap for patient stratification and trials of targeted therapies.
Aberrant AKT activation is prevalent across multiple human cancer lineages providing an important new target for therapy. Twenty-two independent phosphorylation sites have been identified on specific AKT isoforms likely contributing to differential isoform regulation. However, the mechanisms regulating phosphorylation of individual AKT isoform molecules have not been elucidated due to the lack of robust approaches able to assess phosphorylation of multiple sites on a single AKT molecule. Using a nanofluidic proteomic immunoassay (NIA), consisting of isoelectric focusing followed by sensitive chemiluminescence detection, we demonstrate that under basal and ligand-induced conditions that the pattern of phosphorylation events is markedly different between AKT1 and AKT2. Indeed, there are at least 12 AKT1 peaks and at least 5 AKT2 peaks consistent with complex combinations of phosphorylation of different sites on individual AKT molecules. Following insulin stimulation, AKT1 was phosphorylated at Thr308 in the T-loop and Ser473 in the hydrophobic domain. In contrast, AKT2 was only phosphorylated at the equivalent sites (Thr309 and Ser474) at low levels. Further, Thr308 and Ser473 phosphorylation occurred predominantly on the same AKT1 molecules, whereas Thr309 and Ser474 were phosphorylated primarily on different AKT2 molecules. While basal AKT2 phosphorylation was sensitive to inhibition of PI3K, basal AKT1 phosphorylation was essentially resistant. PI3K inhibition decreased pThr451 on AKT2 but not pThr450 on AKT1. Thus NIA technology provides an ability to characterize coordinate phosphorylation of individual AKT molecules providing important information about AKT isoform-specific phosphorylation, which is required for optimal development and implementation of drugs targeting aberrant AKT activation.
AKT phosphorylation NIA
Traditionally, scientific research has focused on studying individual events, such as single mutations, gene function or the effect of the manipulation of one protein on a biological phenotype. A range of technologies, combined with the ability to develop robust and predictive mathematical models, is beginning to provide information that will enable a holistic view of how the genomic and epigenetic aberrations in cancer cells can alter the homeostasis of signalling networks within these cells, between cancer cells and the local microenvironment, at the organ and organism level. This systems biology process needs to be integrated with an iterative approach wherein hypotheses and predictions that arise from modelling are refined and constrained by experimental evaluation. Systems biology approaches will be vital for developing and implementing effective strategies to deliver personalized cancer therapy. Specifically, these approaches will be important to select those patients most likely to benefit from targeted therapies as well as for the development and implementation of rational combinatorial therapies. Systems biology can help to increase therapy efficacy or bypass the emergence of resistance, thus converting the current (often short term) effects of targeted therapies into durable responses, ultimately to improve quality of life and provide a cure.
The serine/threonine kinase AKT is a key mediator of cancer cell survival. We demonstrate that transient glucose deprivation modestly induces AKT phosphorylation at both Thr308 and Ser473. In contrast, prolonged glucose deprivation induces selective AKTThr308 phosphorylation and phosphorylation of a distinct subset of AKT downstream targets leading to cell survival under metabolic stress. Glucose deprivation-induced AKTThr308 phosphorylation is dependent on PDK1 and PI3K but not EGFR or IGF1R. Prolonged glucose deprivation induces the formation of a complex of AKT, PDK1, and the GRP78 chaperone protein, directing phosphorylation of AKTThr308 but AKTSer473. Our results reveal a novel mechanism of AKT activation under prolonged glucose deprivation that protects cells from metabolic stress. The selective activation of AKTThr308 phosphorylation that occurs during prolonged nutrient deprivation may provide an unexpected opportunity for the development and implementation of drugs targeting cell metabolism and aberrant AKT signaling.
glucose deprivation; site specific AKT phosphorylation; substrate specific AKT activation; cell survival
Lysophosphatidic acid (LPA, 1- or 2-acyl-sn-glycerol 3-phosphate) mediates a plethora of physiological and pathological activities via interactions with a series of high affinity G protein-coupled receptors (GPCR). Both LPA receptor family members and autotaxin (ATX/LysoPLD), the primary LPA-producing enzyme, are aberrantly expressed in many human breast cancers and several other cancer lineages. Using transgenic mice expressing either an LPA receptor or ATX, we recently demonstrated that the ATX-LPA receptor axis plays a causal role in breast tumorigenesis and cancer-related inflammation, further validating the ATX-LPA receptor axis as a rich therapeutic target in cancer.
Breast cancer; ATX; LPA; G protein-coupled receptor; inflammation; cytokines; target therapy
Lysophosphatidic acid (LPA) acts through high affinity G protein-coupled receptors to mediate a plethora of physiological and pathological activities associated with tumorigenesis. LPA receptors and autotaxin (ATX/LysoPLD), the primary enzyme producing LPA, are aberrantly expressed in multiple cancer lineages. However, the role of ATX and LPA receptors in the initiation and progression of breast cancer has not been evaluated. We demonstrate that expression of ATX or each Edg-family LPA receptor in mammary epithelium of transgenic mice is sufficient to induce a high frequency of late-onset, estrogen receptor (ER) positive, invasive and metastatic mammary cancer. Thus ATX and LPA receptors can contribute to the initiation and progression of breast cancer.
LPA; ATX; Transgenic mouse model; Breast cancer; Metastasis
Mutations in BRCA1/2 increase the risk of developing breast and ovarian cancer. Germline BRCA1/2 mutations occur in 8.6-13.7% of unselected epithelial ovarian cancers, somatic mutations are also frequent. BRCA1/2 mutated or dysfunctional cells may be sensitive to PARP inhibition by synthetic lethality. The aim of this study is to comprehensively characterise the BRCA1/2 status of a large panel of ovarian cancer cell lines available to the research community to assist in biomarker studies of novel drugs and in particular of PARP inhibitors.
The BRCA1/2 genes were sequenced in 41 ovarian cell lines, mRNA expression of BRCA1/2 and gene methylation status of BRCA1 was also examined. The cytotoxicity of PARP inhibitors olaparib and veliparib was examined in 20 cell lines.
The cell line SNU-251 has a deleterious BRCA1 mutation at 5564G>A, and is the only deleterious BRCA1/2 mutant in the panel. Two cell lines (UPN-251 and PEO1) had deleterious mutations as well as additional reversion mutations that restored the protein functionality. Heterozygous mutations in BRCA1/2 were relatively common, found in 14.6% of cell lines. BRCA1 was methylated in two cell lines (OVCAR8, A1847) and there was a corresponding decrease in gene expression. The BRCA1 methylated cell lines were more sensitive to PARP inhibition than wild-type cells. The SNU-251 deleterious mutant was more sensitive to PARP inhibition, but only in a long-term exposure to correct for its slow growth rate. Cell lines derived from metastatic disease are significantly more resistant to veliparib (2.0 fold p = 0.03) compared to those derived from primary tumours. Resistance to olaparib and veliparib was correlated Pearsons-R 0.5393, p = 0.0311.
The incidence of BRCA1/2 deleterious mutations 1/41 cell lines derived from 33 different patients (3.0%) is much lower than the population incidence. The reversion mutations and high frequency of heterozygous mutations suggest that there is a selective pressure against BRCA1/2 in cell culture similar to the selective pressure seen in the clinic after treatment with chemotherapy. PARP inhibitors may be useful in patients with BRCA1 deleterious mutations or gene methylation.
BRCA1/2; ovarian; mutation; methylation; parp inhibitor; olaparib; veliparib
Although breast cancers are known to be molecularly heterogeneous, their metabolic phenotype is less well understood and may predict response to chemotherapy. This study aimed to evaluate metabolic genes as individual predictive biomarkers in breast cancer.
mRNA microarray data from breast cancer cell lines were used to identify bimodal genes – those with highest potential for robust high/low classification in clinical assays. Metabolic function was evaluated in vitro for the highest scoring metabolic gene, lactate dehydrogenase B (LDHB). Its expression was associated with neoadjuvant chemotherapy response and relapse within clinical and PAM50-derived subtypes.
LDHB was highly expressed in cell lines with glycolytic, basal-like phenotypes. Stable knockdown of LDHB in cell lines reduced glycolytic dependence, linking LDHB expression directly to metabolic function. Using patient datasets, LDHB was highly expressed in basal-like cancers and could predict basal-like subtype within clinical groups (odds ratio = 21 for hormone-receptor (HR)-positive/HER2-negative; odds ratio = 10 for triple-negative). Furthermore, high LDHB predicted pathological complete response (pCR) to neoadjuvant chemotherapy for both HR-positive/HER2-negative (odds ratio = 4.1, P < .001) and triple-negative (odds ratio = 3.0, P = .003) cancers. For triple-negative tumors without pCR, high LDHB post-treatment also identified proliferative tumors with increased risk of recurrence (hazard ratio = 2.2, P = .006).
Expression of LDHB predicted response to neoadjuvant chemotherapy within clinical subtypes independently of standard prognostic markers and PAM50-subtyping. These observations support prospective clinical evaluation of LDHB as a predictive marker of response for breast cancer patients receiving neoadjuvant chemotherapy.
breast cancer; LDHB; lactate; glycolysis; Warburg
The AMP-activated protein kinase (AMPK) functions to monitor and maintain energy homeostasis at the cellular and organismal level. AMPK was perceived historically primarily as a component of the LKB1/STK11 tumor suppressor (LKB1 mutations cause the Peutz Jegher’s cancer predisposition syndrome, PJS) cascade upstream of the TSC1/2/mTOR pathway and thus likely to be a tumor suppressor. However, AMPK has recently been demonstrated to promote cancer cell survival in the face of extrinsic and intrinsic stressors including bioenergetic, growth factor and oncogene stress(1-3), compatible with studies showing that AMPK is required for oncogenic transformation(4). Thus whether AMPK acts as a bona fide tumor suppressor or a contextual oncogene and, of particular importance, whether AMPK should be targeted for activation or inhibition during cancer therapy is controversial and requires clarification. We aim to initiate discussions of these critical questions by reviewing the role of AMPK with an emphasis on cancer cell adaptation to microenvironment stress and therapeutic intervention. Overall AMPK function is the topic of several comprehensive reviews(5-7).
The management of melanoma has evolved due to improved understanding of its molecular drivers. To augment the current understanding of the prevalence, patterns, and associations of mutations in this disease, the results of clinical testing of 699 advanced melanoma patients using a pan-cancer next generation sequencing (NGS) panel of hotspot regions in 46 genes were reviewed. Mutations were identified in 43 of the 46 genes on the panel. The most common mutations were BRAFV600 (36%), NRAS (21%), TP53 (16%), BRAFNon-V600 (6%), and KIT (4%). Approximately one-third of melanomas had >1 mutation detected, and the number of mutations per tumor was associated with melanoma subtype. Concurrent TP53 mutations were the most frequent event in tumors with BRAFV600 and NRAS mutations. Melanomas with BRAFNon-V600 mutations frequently harbored concurrent NRAS mutations (18%), which were rare in tumors with BRAFV600 mutations (1.6%). The prevalence of BRAFV600 and KIT mutations were significantly associated with melanoma subtypes, and BRAFV600 and TP53 mutations were significantly associated with cutaneous primary tumor location. Multiple potential therapeutic targets were identified in metastatic unknown primary and cutaneous melanomas that lacked BRAFV600 and NRAS mutations. These results enrich our understanding of the patterns and clinical associations of oncogenic mutations in melanoma.
melanoma; sequencing; mutation; BRAF; NRAS; TP53
We report here on experimental and theoretical efforts to determine how best to combine drugs that inhibit HER2 and AKT in HER2+ breast cancers. We accomplished this by measuring cellular and molecular responses to lapatinib and the AKT inhibitors (AKTi) GSK690693 and GSK2141795 in a panel of 22 HER2+ breast cancer cell lines carrying wild type or mutant PIK3CA. We observed that combinations of lapatinib plus AKTi were synergistic in HER2+/PIK3CAmut cell lines but not in HER2+/PIK3CAwt cell lines. We measured changes in phospho-protein levels in 15 cell lines after treatment with lapatinib, AKTi or lapatinib + AKTi to shed light on the underlying signaling dynamics. This revealed that p-S6RP levels were less well attenuated by lapatinib in HER2+/PIK3CAmut cells compared to HER2+/PIK3CAwt cells and that lapatinib + AKTi reduced p-S6RP levels to those achieved in HER2+/PIK3CAwt cells with lapatinib alone. We also found that that compensatory up-regulation of p-HER3 and p-HER2 is blunted in PIK3CAmut cells following lapatinib + AKTi treatment. Responses of HER2+ SKBR3 cells transfected with lentiviruses carrying control or PIK3CAmut sequences were similar to those observed in HER2+/PIK3CAmut cell lines but not in HER2+/PIK3CAwt cell lines. We used a nonlinear ordinary differential equation model to support the idea that PIK3CA mutations act as downstream activators of AKT that blunt lapatinib inhibition of downstream AKT signaling and that the effects of PIK3CA mutations can be countered by combining lapatinib with an AKTi. This combination does not confer substantial benefit beyond lapatinib in HER2+/PIK3CAwt cells.
PI3K/Akt/mTOR signaling is being actively pursued as a therapeutic target for breast cancer. We sought to determine if tumor heterogeneity and biospecimen variables affect the evaluation of PI3K/Akt/mTOR pathway markers.
Intraoperative image-guided core-needle biopsies (CNB), and central and peripheral surgical tumor specimens were prospectively collected in 53 patients with invasive breast cancer. Specimens were assessed with reverse phase protein arrays (RPPA) and immunohistochemistry (IHC).
There was a moderate or strong correlation between the expression of 149 (97%) of the 154 different RPPA markers in the center and periphery. Correlation was higher for smaller tumors, in patients who did not undergo neoadjuvant therapy, and with shorter cold ischemia time. Of 154 markers, 132 (86%) were not statistically different between the center and periphery, and 97 (63%) were not different between the CNB and the surgical specimen (average of the central and peripheral specimen). pAkt S473 and PTEN had a significant correlation between central and peripheral specimens, and between CNB and surgical specimen. However, pAkt S473, pS6 S235/236 and pS6 240/244 levels were significantly higher in CNB than the central specimens both by RPPA and by IHC.
Most individual proteomic biomarkers studied do not have significant intratumoral heterogeneity. However, protein and phosphoprotein levels are affected by biospecimen type and other preanalytic variables. PI3K pathway activation is greater in CNB compared to post-excision surgical samples suggesting a potential loss of phosphorylation during surgical manipulation, or with cold ischemia of surgical specimens.
biopsy; proteomics; tumor heterogeneity; Akt
High grade serous ovarian carcinoma (HGSOC) and basal-like breast cancer (BLBC) share many features including TP53 mutations, genomic instability and poor prognosis. We recently reported that Elafin is overexpressed by HGSOC and is associated with poor overall survival. Here, we confirmed that Elafin overexpression is associated with shorter survival in 1000 HGSOC patients. Elafin confers a proliferative advantage to tumor cells through activation of the MAP kinase pathway. This mitogenic effect can be neutralized by RNA interference, specific antibodies, and a MEK inhibitor. Elafin expression in patient-derived samples was also associated with chemoresistance and strongly correlates with bcl-xL expression. We extended these findings into examination of 1100 primary breast tumors and six breast cancer cell lines. We observed that Elafin is overexpressed and secreted specifically by BLBC tumors and cell lines, leading to a similar mitogenic effect through activation of the MAP kinase pathway. Here too, Elafin overexpression is associated with poor overall survival, suggesting that it may serve as a biomarker and therapeutic target in this setting.
Elafin; ovarian cancer; basal-like breast cancer; mitogen; MAP kinase
Autophagy, a programmed process in which cell contents are delivered to lysosomes for degradation, appears to have both tumor-suppressive and tumor-promoting functions; both stimulation and inhibition of autophagy have been reported to induce cancer cell death, and particular genes and proteins have been associated both positively and negatively with autophagy. To provide a basis for incisive analysis of those complexities and ambiguities and to guide development of new autophagy-targeted treatments for cancer, we have compiled a comprehensive, curated inventory of autophagy modulators by integrating information from published siRNA screens, multiple pathway analysis algorithms, and extensive, manually curated text-mining of the literature. The resulting inventory includes 739 proteins and 385 chemicals (including drugs, small molecules, and metabolites). Because autophagy is still at an early stage of investigation, we provide extensive analysis of our sources of information and their complex relationships with each other. We conclude with a discussion of novel strategies that could potentially be used to target autophagy for cancer therapy.
autophagy; cancer; high-throughput screening; L-asparaginase; natural language processing; pathway analysis; RNAi; siRNA; text-mining
Multiple mechanisms have been described that confer BRAF inhibitor resistance to melanomas, yet the basis of this resistance remains undefined in a sizable portion of patient samples. Here, we characterized samples from a set of patients with melanoma that included individuals at baseline diagnosis, on BRAF inhibitor treatment, and with resistant tumors at both the protein and RNA levels. Using RNA and DNA sequencing, we identified known resistance-conferring mutations in 50% (6 of 12) of the resistant samples. In parallel, targeted proteomic analysis by protein array categorized the resistant samples into 3 stable groups, 2 of which were characterized by reactivation of MAPK signaling to different levels and 1 that was MAPK independent. The molecular relevance of these classifications identified in patients was supported by both mutation data and the similarity of resistance patterns that emerged during a co-clinical trial in a genetically engineered mouse (GEM) model of melanoma that recapitulates the development of BRAF inhibitor resistance. Additionally, we defined candidate biomarkers in pre- and early-treatment patient samples that have potential for predicting clinical responses. On the basis of these observations, we suggest that BRAF inhibitor–resistant melanomas can be actionably classified using protein expression patterns, even without identification of the underlying genetic alteration.
Many mutant p53 proteins (mutp53s) exert oncogenic gain-of-function (GOF) properties, but the mechanisms mediating these functions remain poorly defined. We show here that GOF mutp53s inhibit AMP-activated protein kinase (AMPK) signaling in head and neck cancer cells. Conversely, downregulation of GOF mutp53s enhances AMPK activation under energy stress, decreasing the activity of the anabolic factors acetyl-CoA carboxylase and ribosomal protein S6 and inhibiting aerobic glycolytic potential and invasive cell growth. Under conditions of energy stress, GOF mutp53s, but not wild-type p53, preferentially bind to the AMPKα subunit and inhibit AMPK activation. Given the importance of AMPK as an energy sensor and tumor suppressor that inhibits anabolic metabolism, our findings reveal that direct inhibition of AMPK activation is an important mechanism through which mutp53s can gain oncogenic function.
Aberrant von Hippel Lindau (VHL) protein function is the underlying driver of VHL-related diseases, including both sporadic and inherited clear cell renal cell carcinoma (ccRCC). About one third of VHL mutations are missense point mutations, with R167Q being the most common VHL point mutation in hereditary VHL disease. Although it has been studied extensively, the ability of VHL-R167Q to downregulate hypoxia inducible factor 2α (HIF2α) is still controversial. In addition, the manner in which the mutation contributes to tumorigenesis is not fully understood. No therapeutic approach is available to target VHL-R167Q and similar missense point mutations. We analyzed VHL-R167Q proteostasis and function at normoxia, at hypoxia with different oxygen pressure, and in a xenograft mouse model. We showed that the protein levels of VHL-R167Q dictate its ability to downregulate HIF2α and suppress tumor growth. Strikingly, the proteasome inhibitors bortezomib and carfilzomib, which are currently in clinical use, stabilize VHL-R167Q and increase its ability to downregulate HIF2α. VHL-R167Q binds elongin C and elongin B with considerably less avidity than wild-type VHL does but retains residual capacity to generate a VHL-elongin C-elongin B complex, downregulate HIF2α, and suppress tumorigenesis, which could be rescued by increase VHL-R167Q levels. Finally, we used in silico approaches and identified other missense VHL mutants in addition to VHL-R167Q that might be rescued by similar strategies. Thus, our studies revealed detailed information describing how VHL-R167Q contributes to tumorigenesis and identified a potential targeted therapy for ccRCC and other VHL-related disease in patients carrying VHL-R167Q or similar missense mutations.
VHL; ccRCC; Proteasome Inhibition; Point Mutations; Functional Restoration
During the type-setting of the final version of the article  some of the additional files were swapped. The correct files are republished in this Erratum.
The online version of the original article can be found under doi:10.1186/s13059-015-0658-5.
HPV-associated (HPV+) oropharyngeal squamous cell carcinomas (OPSCC) have different molecular and biological characteristics and clinical behavior compared to HPV-negative (HPV−) OPSCC. PIK3CA mutations are more common in HPV(+)OPSCC. To define molecular differences and tumor subsets, protein expression and phosphorylation were compared between HPV(+) and HPV(−) OPSCC and between tumors with and without PIK3CA mutations.
Expression of 137 total and phosphorylated proteins was evaluated by reverse phase protein array(RPPA) in 29 HPV(+) and 13 HPV(−)prospectively collected OPSCCs. 47 OPSCCs were tested for hotspot activating mutations in PIK3CA and AKT. Activation of PIK3CA downstream targets and sensitivity to pathway inhibitors were determined in HPV(+) head and neck cancer cells overexpressing wild-type or mutant PIK3CA.
Analyses revealed forty-one differentially expressed proteins between HPV(+) and HPV(−) OPSCC categorized into functional groups: DNA repair, cell cycle, apoptosis, PI3K/AKT/mTOR, and receptor kinase pathways. All queried DNA repair proteins were significantly upregulated in HPV(+) samples. 8 of 33 HPV(+) and 0 of 14 HPV(−) tumors contained activating PIK3CA mutations. Despite all activating PIK3CA mutations occurring in HPV(+) samples, HPV(+) tumors had lower mean levels of activated AKT and downstream AKT target phosphorylation. Ectopic expression of mutant PIK3CA in HPV(+) cells increase dmTOR, but not AKT activity. HPV E6/E7 overexpression inhibited AKT phosphorylation in HPV-negative cells. Mutant PIK3CA overexpressing cells were more sensitive to a dual PI3K/mTOR inhibitor compared to an AKT inhibitor.
Protein expression analyses suggest that HPV(+) and HPV(−) OPSCC differentially activate DNA repair, cell cycle, apoptosis, PI3K/AKT/mTOR, and receptor kinase pathways. PIK3CA mutations are more common in HPV(+) OPSCC and are associated with activation of mTOR, but not AKT. These data suggest that inhibitors for mTOR may have activity against HPV(+) PIK3CA mutant oropharyngeal cancers.
head and neck cancer; PIK3CA; AKT; HPV; RPPA; prognosis
There is growing interest in delivering genomically-informed cancer therapy. Our aim was to determine the concordance of genomic alterations between primary and recurrent breast cancer. Targeted next generation sequencing was performed on formalin-fixed paraffin embedded (FFPE) samples, profiling 3320 exons of 182 cancer-related genes plus 37 introns from 14 genes often rearranged in cancer. Point mutations, indels, copy number alterations and select rearrangements were assessed in 74 tumors from 43 patients (36 primary, 38 recurrence/metastases). Alterations potentially targetable with established or investigational therapeutics were considered “actionable”. Alterations were detected in 55 genes (mean 3.95 alterations/sample, range 1-12), including mutations in PIK3CA, TP53, ARID1A, PTEN, AKT1, NF1, FBXW7 and FGFR3 and amplifications in MCL1, CCND1, FGFR1, MYC, IGF1R, MDM2, MDM4, AKT3, CDK4, AKT2. In 33 matched primary and recurrent tumors, 97 of 112 (86.6%) somatic mutations were concordant. Of identified copy number alterations, 136 of 159 (85.5%) were concordant: 37 (23.3%) were concordant, but below the reporting threshold in one of the matched samples, and 23 (14.5%) discordant. There was an increased frequency of CDK4/MDM2 amplifications in recurrences, as well as gains and losses of other actionable alterations. 40 of 43 (93%) patients had actionable alterations that could inform targeted treatment options. In conclusion, deep genomic profiling of cancer-related genes reveals potentially actionable alterations in most breast cancer patients. Overall there was high concordance between primary and recurrent tumors. Analysis of recurrent tumors prior to treatment may provide additional insights, as both gains and losses of targets are observed.
breast cancer; genomics; next generation sequencing; genomic evolution
Motivation: Identifying subclonal mutations and their implications requires accurate estimation of mutant allele fractions from possibly duplicated sequencing reads. Removing duplicate reads assumes that polymerase chain reaction amplification from library constructions is the primary source. The alternative—sampling coincidence from DNA fragmentation—has not been systematically investigated.
Results: With sufficiently high-sequencing depth, sampling-induced read duplication is non-negligible, and removing duplicate reads can overcorrect read counts, causing systemic biases in variant allele fraction and copy number variation estimations. Minimal overcorrection occurs when duplicate reads are identified accounting for their mate reads, inserts are of a variety of lengths and samples are sequenced in separate batches. We investigate sampling-induced read duplication in deep sequencing data with 500× to 2000× duplicates-removed sequence coverage. We provide a quantitative solution to overcorrection and guidance for effective designs of deep sequencing platforms that facilitate accurate estimation of variant allele fraction and copy number variation.
Availability and implementation: A Python implementation is freely available at https://bitbucket.org/wanding/duprecover/overview.
Contact: email@example.com, firstname.lastname@example.org
Supplementary information: Supplementary data are available at Bioinformatics online.
In this issue of Cancer Discovery, AI-Ahmadie and colleagues identify a somatic mutation in the Rad50 gene as a likely contributing factor to an unusual curative response to systemic combination therapy employing the DNA-damaging agent, irinotecan, and a checkpoint kinase 1(Chk1) inhibitor in a patient with recurrent, metastatic small-cell cancer. This study highlights the importance of in-depth analysis of exceptional responders to chemo and targeted therapy in early phase clinical trials and opens new avenues for developing cancer genome-based combination therapy to improve the efficacy of traditional chemotherapy through synthetically lethal interactions.
Estrogen receptor–positive (ER+) breast cancers adapt to hormone deprivation and become resistant to antiestrogen therapy. Here, we performed deep sequencing on ER+ tumors that remained highly proliferative after treatment with the aromatase inhibitor letrozole and identified a D189Y mutation in the inhibitory SH2 domain of the SRC family kinase (SFK) LYN. Evaluation of 463 breast tumors in The Cancer Genome Atlas revealed four LYN mutations, two of which affected the SH2 domain. In addition, LYN was upregulated in multiple ER+ breast cancer lines resistant to long-term estrogen deprivation (LTED). An RNAi-based kinome screen revealed that LYN is required for growth of ER+ LTED breast cancer cells. Kinase assays and immunoblot analyses of SRC substrates in transfected cells indicated that LYND189Y has higher catalytic activity than WT protein. Further, LYND189Y exhibited reduced phosphorylation at the inhibitory Y507 site compared with LYNWT. Other SH2 domain LYN mutants, E159K and K209N, also exhibited higher catalytic activity and reduced inhibitory site phosphorylation. LYND189Y overexpression abrogated growth inhibition by fulvestrant and/or the PI3K inhibitor BKM120 in 3 ER+ breast cancer cell lines. The SFK inhibitor dasatinib enhanced the antitumor effect of BKM120 and fulvestrant against estrogen-deprived ER+ xenografts but not LYND189Y-expressing xenografts. These results suggest that LYN mutations mediate escape from antiestrogens in a subset of ER+ breast cancers.
Neoangiogenesis is an important feature in tumor growth and progression, and combining chemotherapy and antiangiogenic drugs have demonstrated clinical efficacy. However, as treatment induced resistance often develops our goal was to identify pathways indicating response and/or evolving resistance to treatment, and inhibit these pathways to optimize the treatment strategies.
To identify markers of response and/or resistance Reverse Phase Protein Array (RPPA) was utilized to characterize treatment-induced changes in a bevacizumab responsive and a nonresponsive human breast cancer xenograft. Results were combined with bioinformatic modeling to predict druggable targets for optimization of the treatment.
RPPA analysis showed that both tumor models responded to bevacizumab with an early (day 3) upregulation of growth factor receptors and downstream signaling pathways, with persistent mTOR signaling until the end of the in vivo experiment. Adding doxorubicin to bevacizumab showed significant and superior growth inhibition of basal-like tumors, whereas no additive effect was seen in the luminal-like model. The combination treatment corresponded to a continuous late attenuation of mTOR signaling in the basal-like model, while the inhibition was temporary in the luminal-like model. Integrating the bevacizumab-induced dynamic changes in protein levels with bioinformatic modeling predicted inhibition of PI3K-pathway to increase the efficacy of bevacizumab monotherapy. In vivo experiments combining bevacizumab and the PI3K/mTOR inhibitor BEZ235 confirmed their significant and additive growth inhibitory effect in the basal-like model.
Treatment with bevacizumab caused compensatory upregulation of several signaling pathways. Targeting such pathways increased the efficacy of antiangiogenic therapy.