Integrative analysis of mutational and transcriptional profiles reveals driver mutations of metastatic breast cancers

Despite the explosion in the numbers of cancer genomic studies, metastasis is still the major cause of cancer mortality. In breast cancer, approximately one-fifth of metastatic patients survive 5 years. Therefore, detecting the patients at a high risk of developing distant metastasis at first diagnosis is critical for effective treatment strategy. We hereby present a novel systems biology approach to identify driver mutations escalating the risk of metastasis based on both exome and RNA sequencing of our collected 78 normal-paired breast cancers. Unlike driver mutations occurring commonly in cancers as reported in the literature, the mutations detected here are relatively rare mutations occurring in less than half metastatic samples. By supposing that the driver mutations should affect the metastasis gene signatures, we develop a novel computational pipeline to identify the driver mutations that affect transcription factors regulating metastasis gene signatures. We identify driver mutations in ADPGK, NUP93, PCGF6, PKP2 and SLC22A5, which are verified to enhance cancer cell migration and prompt metastasis with in vitro experiments. The discovered somatic mutations may be helpful for identifying patients who are likely to develop distant metastasis.

[1]  Tian-Li Wang,et al.  Frequent Mutations of Chromatin Remodeling Gene ARID1A in Ovarian Clear Cell Carcinoma , 2010, Science.

[2]  Shanshan Yang,et al.  Aberrant expression of δ-like ligand 4 contributes significantly to axillary lymph node metastasis and predicts postoperative outcome in breast cancer. , 2014, Human pathology.

[3]  K. Possinger,et al.  PPARgamma ligands and ATRA inhibit the invasion of human breast cancer cells in vitro. , 2003, Breast cancer research and treatment.

[4]  Steven A. Roberts,et al.  Mutational heterogeneity in cancer and the search for new cancer-associated genes , 2013 .

[5]  L. Fallowfield,et al.  Locally recurrent or metastatic breast cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. , 2010, Annals of oncology : official journal of the European Society for Medical Oncology.

[6]  M. J. Kwon,et al.  Emerging Roles of Claudins in Human Cancer , 2013, International journal of molecular sciences.

[7]  Ming Yan,et al.  Plakophilin-2 Promotes Tumor Development by Enhancing Ligand-Dependent and -Independent Epidermal Growth Factor Receptor Dimerization and Activation , 2014, Molecular and Cellular Biology.

[8]  E. Lander,et al.  Lessons from the Cancer Genome , 2013, Cell.

[9]  D. Saur,et al.  IKKα controls canonical TGFβ-SMAD signaling to regulate genes expressing SNAIL and SLUG during EMT in Panc1 cells , 2011, Development.

[10]  David Chen,et al.  ESR1 ligand binding domain mutations in hormone-resistant breast cancer , 2013, Nature Genetics.

[11]  Tao Xie,et al.  Inferring causal genomic alterations in breast cancer using gene expression data , 2011, BMC Systems Biology.

[12]  Craig D. Shriver,et al.  A gene expression signature that defines breast cancer metastases , 2008, Clinical & Experimental Metastasis.

[13]  C. Theillet,et al.  Snail and Slug Play Distinct Roles during Breast Carcinoma Progression , 2006, Clinical Cancer Research.

[14]  Douglas S. Darling,et al.  Zeb1 links epithelial-mesenchymal transition and cellular senescence , 2008, Development.

[15]  Robert C. Bast,et al.  HER2-targeting Antibodies Modulate the Cyclin-dependent Kinase Inhibitor p27Kip1 via Multiple Signaling Pathways , 2005, Cell cycle.

[16]  Haiyuan Yu,et al.  Detecting overlapping protein complexes in protein-protein interaction networks , 2012, Nature Methods.

[17]  Iannis Aifantis,et al.  Oncogenic and tumor suppressor functions of Notch in cancer: it’s NOTCH what you think , 2011, The Journal of experimental medicine.

[18]  Wei Hu,et al.  Reduced expression of Snail decreases breast cancer cell motility by downregulating the expression and inhibiting the activity of RhoA GTPase , 2013, Oncology letters.

[19]  Samy Lamouille,et al.  TGF-β-induced epithelial to mesenchymal transition , 2009, Cell Research.

[20]  A. Harris,et al.  Vascular Biology , Atherosclerosis and Endothelium Biolog y Expression of Vascular Notch Ligand Delta-Like 4 and Inflammatory Markers in Breast Cancer , 2010 .

[21]  J. Gray,et al.  CUL4A induces epithelial-mesenchymal transition and promotes cancer metastasis by regulating ZEB1 expression. , 2014, Cancer research.

[22]  R. Gambari,et al.  Decoy molecules based on PNA-DNA chimeras and targeting Sp1 transcription factors inhibit the activity of urokinase-type plasminogen activator receptor (uPAR) promoter. , 2005, Oncology research.

[23]  J. Baell,et al.  LIM kinase inhibition reduces breast cancer growth and invasiveness but systemic inhibition does not reduce metastasis in mice , 2012, Clinical & Experimental Metastasis.

[24]  Meng Li,et al.  Somatic Mutations in the Chromatin Remodeling Gene ARID1A Occur in Several Tumor Types , 2011, Human mutation.

[25]  Andy J. Minn,et al.  Genes that mediate breast cancer metastasis to lung , 2005, Nature.

[26]  Shi-Hua Zhang,et al.  Efficient methods for identifying mutated driver pathways in cancer , 2012, Bioinform..

[27]  Zhengyan Kan,et al.  Exome sequencing identifies frequent mutation of ARID1A in molecular subtypes of gastric cancer , 2011, Nature Genetics.

[28]  K. Possinger,et al.  PPARγ Ligands and ATRA Inhibit the Invasion of Human Breast Cancer Cells in vitro , 2003, Breast Cancer Research and Treatment.

[29]  J. Zavadil,et al.  TGF-beta and epithelial-to-mesenchymal transitions. , 2005, Oncogene.

[30]  Gary D. Bader,et al.  HyperModules: identifying clinically and phenotypically significant network modules with disease mutations for biomarker discovery , 2014, Bioinform..

[31]  B. Joseph,et al.  The Cdk inhibitor p57Kip2 controls LIM-kinase 1 activity and regulates actin cytoskeleton dynamics , 2009, Oncogene.

[32]  E. Sahai,et al.  LIM kinases are required for invasive path generation by tumor and tumor-associated stromal cells , 2010, The Journal of cell biology.

[33]  W. Gerald,et al.  Genes that mediate breast cancer metastasis to the brain , 2009, Nature.

[34]  Xing-Ming Zhao,et al.  Automatic Modeling of Signaling Pathways by Network Flow Model , 2009, J. Bioinform. Comput. Biol..

[35]  Y. L. Ramachandra,et al.  Human Proteinpedia enables sharing of human protein data , 2008, Nature Biotechnology.

[36]  P. Lobie,et al.  PAX5α Enhances the Epithelial Behavior of Human Mammary Carcinoma Cells , 2010, Molecular Cancer Research.

[37]  T. Lu,et al.  Loss of USF transcriptional activity in breast cancer cell lines , 1999, Oncogene.

[38]  A. Børresen-Dale,et al.  The landscape of cancer genes and mutational processes in breast cancer , 2012, Nature.

[39]  K. Kinzler,et al.  Cancer Genome Landscapes , 2013, Science.

[40]  Terrence S. Furey,et al.  The UCSC Table Browser data retrieval tool , 2004, Nucleic Acids Res..

[41]  E. Castanas,et al.  Immunohistochemical study of pElk-1 expression in human breast cancer: association with breast cancer biologic profile and clinicopathologic features. , 2013, Breast.

[42]  S. Natsugoe,et al.  Clinical implications of DLL4 expression in gastric cancer , 2013, Journal of experimental & clinical cancer research : CR.

[43]  Pei-Shan Wang,et al.  Troglitazone inhibits cell migration, adhesion, and spreading by modulating cytoskeletal rearrangement in human breast cancer cells , 2008, Molecular carcinogenesis.

[44]  J. Siedlecki,et al.  HAX-1 overexpression, splicing and cellular localization in tumors , 2010, BMC Cancer.

[45]  Aviv Regev,et al.  Whole exome sequencing of circulating tumor cells provides a window into metastatic prostate cancer , 2014, Nature Biotechnology.

[46]  Jos Jonkers,et al.  Somatic inactivation of E-cadherin and p53 in mice leads to metastatic lobular mammary carcinoma through induction of anoikis resistance and angiogenesis. , 2006, Cancer cell.

[47]  E. Lander,et al.  Loss of E-cadherin promotes metastasis via multiple downstream transcriptional pathways. , 2008, Cancer research.

[48]  S. Fox,et al.  Expression of delta-like ligand 4 (Dll4) and markers of hypoxia in colon cancer , 2009, British Journal of Cancer.

[49]  H. Brauch,et al.  A novel functional polymorphism in the transforming growth factor-beta2 gene promoter and tumor progression in breast cancer. , 2006, Cancer research.

[50]  C. Sander,et al.  Genome-wide analysis of non-coding regulatory mutations in cancer , 2014, Nature Genetics.

[51]  Elias Campo Guerri,et al.  International network of cancer genome projects , 2010 .

[52]  Kazuyuki Aihara,et al.  Detecting early-warning signals for sudden deterioration of complex diseases by dynamical network biomarkers , 2012, Scientific Reports.

[53]  Ahmed H. Mekkawy,et al.  Novel protein interactors of urokinase-type plasminogen activator receptor. , 2010, Biochemical and biophysical research communications.

[54]  M. Kazanietz,et al.  Transcriptional Regulation of Oncogenic Protein Kinase Cϵ (PKCϵ) by STAT1 and Sp1 Proteins* , 2014, The Journal of Biological Chemistry.

[55]  Steven A. Roberts,et al.  Mutational heterogeneity in cancer and the search for new cancer genes , 2014 .

[56]  R. Klopfleisch,et al.  The metastatic cascade is reflected in the transcriptome of metastatic canine mammary carcinomas. , 2011, Veterinary journal.

[57]  Lior Pachter,et al.  Sequence Analysis , 2020, Definitions.

[58]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[59]  I. Shih,et al.  Functional analysis of in-frame indel ARID1A mutations reveals new regulatory mechanisms of its tumor suppressor functions. , 2012, Neoplasia.

[60]  Nur Aishah Taib,et al.  Predicting Survival of De Novo Metastatic Breast Cancer in Asian Women: Systematic Review and Validation Study , 2014, PloS one.

[61]  C. Heldin,et al.  Non-Smad TGF-beta signals. , 2005, Journal of cell science.

[62]  F. Kittrell,et al.  Overexpression of Separase induces aneuploidy and mammary tumorigenesis , 2008, Proceedings of the National Academy of Sciences.

[63]  Eduard Batlle,et al.  Snail Induction of Epithelial to Mesenchymal Transition in Tumor Cells Is Accompanied by MUC1 Repression andZEB1 Expression* , 2002, The Journal of Biological Chemistry.

[64]  Jerzy Adamski,et al.  IKK&agr; controls canonical TGF&bgr;–SMAD signaling to regulate genes expressing SNAIL and SLUG during EMT in Panc1 cells , 2013, Journal of Cell Science.

[65]  E. Leygue,et al.  Claudin 1 in Breast Tumorigenesis: Revelation of a Possible Novel “Claudin High” Subset of Breast Cancers , 2010, Journal of biomedicine & biotechnology.

[66]  C. Heldin,et al.  Non-Smad TGF-β signals , 2005, Journal of Cell Science.

[67]  O. Bernard,et al.  A role for LIM kinase in cancer invasion , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[68]  P. Brown,et al.  SLC22A5/OCTN2 expression in breast cancer is induced by estrogen via a novel intronic estrogen-response element (ERE) , 2012, Breast Cancer Research and Treatment.

[69]  Eli Upfal,et al.  De Novo Discovery of Mutated Driver Pathways in Cancer , 2011, RECOMB.

[70]  A. Sahin,et al.  LIM kinase 1 increases tumor metastasis of human breast cancer cells via regulation of the urokinase‐type plasminogen activator system , 2006, International journal of cancer.

[71]  G. Watkins,et al.  The localisation and reduction of nuclear staining of PPARgamma and PGC-1 in human breast cancer. , 2004, Oncology reports.