Hunt for the tipping point during endocrine resistance process in breast cancer by dynamic network biomarkers

Abstract Acquired drug resistance is the major reason why patients fail to respond to cancer therapies. It is a challenging task to determine the tipping point of endocrine resistance and detect the associated molecules. Derived from new systems biology theory, the dynamic network biomarker (DNB) method is designed to quantitatively identify the tipping point of a drastic system transition and can theoretically identify DNB genes that play key roles in acquiring drug resistance. We analyzed time-course mRNA sequence data generated from the tamoxifen-treated estrogen receptor (ER)-positive MCF-7 cell line, and identified the tipping point of endocrine resistance with its leading molecules. The results show that there is interplay between gene mutations and DNB genes, in which the accumulated mutations eventually affect the DNB genes that subsequently cause the change of transcriptional landscape, enabling full-blown drug resistance. Survival analyses based on clinical datasets validated that the DNB genes were associated with the poor survival of breast cancer patients. The results provided the detection for the pre-resistance state or early signs of endocrine resistance. Our predictive method may greatly benefit the scheduling of treatments for complex diseases in which patients are exposed to considerably different drugs and may become drug resistant.

[1]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[2]  K. Aihara,et al.  Early Diagnosis of Complex Diseases by Molecular Biomarkers, Network Biomarkers, and Dynamical Network Biomarkers , 2014, Medicinal research reviews.

[3]  Susumu Goto,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 2000, Nucleic Acids Res..

[4]  Mauricio O. Carneiro,et al.  From FastQ Data to High‐Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline , 2013, Current protocols in bioinformatics.

[5]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Chen Li,et al.  Dysfunction of PLA2G6 and CYP2C44-associated network signals imminent carcinogenesis from chronic inflammation to hepatocellular carcinoma , 2017, Journal of molecular cell biology.

[7]  Kazuyuki Aihara,et al.  Identifying critical transitions and their leading biomolecular networks in complex diseases , 2012, Scientific Reports.

[8]  G. Vachtsevanos,et al.  Epileptic Seizures May Begin Hours in Advance of Clinical Onset A Report of Five Patients , 2001, Neuron.

[9]  Wanwei Zhang,et al.  Discovering a critical transition state from nonalcoholic hepatosteatosis to nonalcoholic steatohepatitis by lipidomics and dynamical network biomarkers. , 2016, Journal of molecular cell biology.

[10]  Tilo Winkler,et al.  Self-organized patchiness in asthma as a prelude to catastrophic shifts , 2005, Nature.

[11]  Rui Liu,et al.  Detecting tissue-specific early warning signals for complex diseases based on dynamical network biomarkers: study of type 2 diabetes by cross-tissue analysis , 2014, Briefings Bioinform..

[12]  Kazuyuki Aihara,et al.  Development of a mathematical model that predicts the outcome of hormone therapy for prostate cancer. , 2010, Journal of theoretical biology.

[13]  Jiang-xia Zhao,et al.  Circular RNA is enriched and stable in exosomes: a promising biomarker for cancer diagnosis , 2015, Cell Research.

[14]  Luonan Chen,et al.  Coexpression network analysis in chronic hepatitis B and C hepatic lesions reveals distinct patterns of disease progression to hepatocellular carcinoma. , 2012, Journal of molecular cell biology.

[15]  Jiang Shou,et al.  Development of resistance to targeted therapies transforms the clinically associated molecular profile subtype of breast tumor xenografts. , 2008, Cancer research.

[16]  Yongjun Li,et al.  Detecting critical state before phase transition of complex biological systems by hidden Markov model , 2016, Bioinform..

[17]  M. Zaccolo,et al.  Protein kinase A gating of a pseudopodial-located RhoA/ROCK/p38/NHE1 signal module regulates invasion in breast cancer cell lines. , 2005, Molecular biology of the cell.

[18]  Jyh‐cherng Yu,et al.  Asymptomatic pulmonary nodules in a patient with early-stage breast cancer: Cryptococcus infection. , 2010, International journal of infectious diseases : IJID : official publication of the International Society for Infectious Diseases.

[19]  U. Gunnarsson,et al.  Role of the multidisciplinary team in breast cancer management: results from a large international survey involving 39 countries. , 2012, Annals of oncology : official journal of the European Society for Medical Oncology.

[20]  M. Millward,et al.  Dynamic versus static biomarkers in cancer immune checkpoint blockade: unravelling complexity , 2017, Nature Reviews Drug Discovery.

[21]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[22]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[23]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[24]  M. DePristo,et al.  The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. , 2010, Genome research.

[25]  R. Schiff,et al.  Tamoxifen resistance in breast tumors is driven by growth factor receptor signaling with repression of classic estrogen receptor genomic function. , 2008, Cancer research.

[26]  Simak Ali,et al.  Phosphatidylinositol 3-Kinase/AKT-mediated Activation of Estrogen Receptor α , 2001, The Journal of Biological Chemistry.

[27]  I. Ellis,et al.  Phosphorylation of ERK1/2 mitogen‐activated protein kinase is associated with poor response to anti‐hormonal therapy and decreased patient survival in clinical breast cancer , 2001, International journal of cancer.

[28]  R. Riggins,et al.  Pathways to tamoxifen resistance. , 2007, Cancer letters.

[29]  Minyi Guo,et al.  Cerebrospinal fluid protein dynamic driver network: At the crossroads of brain tumorigenesis. , 2015, Methods.

[30]  Robert L. Sutherland,et al.  Biological determinants of endocrine resistance in breast cancer , 2009, Nature Reviews Cancer.

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

[32]  Patrick E. McSharry,et al.  Prediction of epileptic seizures: are nonlinear methods relevant? , 2003, Nature Medicine.

[33]  Hannah H. Chang,et al.  Cell Fate Decision as High-Dimensional Critical State Transition , 2016, bioRxiv.

[34]  R. Nicholson,et al.  Elevated levels of epidermal growth factor receptor/c-erbB2 heterodimers mediate an autocrine growth regulatory pathway in tamoxifen-resistant MCF-7 cells. , 2003, Endocrinology.

[35]  Paul Theodor Pyl,et al.  HTSeq—a Python framework to work with high-throughput sequencing data , 2014, bioRxiv.

[36]  Davide Heller,et al.  STRING v10: protein–protein interaction networks, integrated over the tree of life , 2014, Nucleic Acids Res..

[37]  Kazuyuki Aihara,et al.  Identifying critical transitions of complex diseases based on a single sample , 2014, Bioinform..

[38]  Meiyi Li,et al.  Dynamic network biomarker indicates pulmonary metastasis at the tipping point of hepatocellular carcinoma , 2018, Nature Communications.

[39]  Wanwei Zhang,et al.  Deciphering early development of complex diseases by progressive module network. , 2014, Methods.

[40]  Rudiyanto Gunawan,et al.  Single-Cell-Based Analysis Highlights a Surge in Cell-to-Cell Molecular Variability Preceding Irreversible Commitment in a Differentiation Process , 2016, PLoS biology.