WO2013062505A1 - Pathway recognition algorithm using data integration on genomic models (paradigm) - Google Patents
Pathway recognition algorithm using data integration on genomic models (paradigm) Download PDFInfo
- Publication number
- WO2013062505A1 WO2013062505A1 PCT/US2011/001844 US2011001844W WO2013062505A1 WO 2013062505 A1 WO2013062505 A1 WO 2013062505A1 US 2011001844 W US2011001844 W US 2011001844W WO 2013062505 A1 WO2013062505 A1 WO 2013062505A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- pathway
- network
- cell
- protein
- disease
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B45/00—ICT specially adapted for bioinformatics-related data visualisation, e.g. displaying of maps or networks
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/22—Indexing; Data structures therefor; Storage structures
- G06F16/2219—Large Object storage; Management thereof
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B5/00—ICT specially adapted for modelling or simulations in systems biology, e.g. gene-regulatory networks, protein interaction networks or metabolic networks
- G16B5/20—Probabilistic models
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B99/00—Subject matter not provided for in other groups of this subclass
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/10—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/30—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to physical therapies or activities, e.g. physiotherapy, acupressure or exercising
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/40—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mechanical, radiation or invasive therapies, e.g. surgery, laser therapy, dialysis or acupuncture
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/60—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to nutrition control, e.g. diets
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Definitions
- PARADIGM GENOMIC MODELS
- the present invention relates to a method for identifying components of biological pathways in an individual or subject and determining if the individual or subject is a candidate for a clinical regimen or treatment.
- the invention also relates to using the methods to diagnose whether a subject is susceptible to cancer, autoimmune diseases, cell cycle disorders, or other disorders.
- a central premise in modern cancer treatment is that patient diagnosis, prognosis, risk assessment, and treatment response prediction can be improved by stratification of cancers based on genomic, transcriptional and epigenomic characteristics of the tumor alongside relevant clinical information gathered at the time of diagnosis (for example, patient history, tumor histology and stage) as well as subsequent clinical follow-up data (for example, treatment regimens and disease recurrence events).
- Cancer is a disease of the genome that is associated with aberrant alterations that lead to disregulation of the cellular system. What is not clear is how genomic changes feed into genetic pathways that underlie cancer phenotypes.
- High-throughput functional genomics investigations have made tremendous progress in the past decade (Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A, Boldrick JC, Sabet H, Tran T, Yu X, Powell JI, Yang L, Marti GE, Moore T, Hudson J, Lu L, Lewis DB, Tibshirani R, SHERLOCK G, Chan WC, Greiner TC, Weisenburger DD, Armitage JO, Warnke R, Levy R, Wilson W, Grever MR, Byrd JC, Botstein D, Brown PO, Staudt LM.
- GSEA gene set enrichment analysis
- Gene set enrichment analysis a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. U.S.A. 2005 Oct.; 102(43): 15545- 15550.). Other approaches include using a hypergeometric test- based method to identify Gene Ontology (Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, SHERLOCK G. Gene ontology: tool for the unification of biology.
- Overrepresentation analyses are limited in their efficacy because they do not incorporate known interdependencies among genes in a pathway that can increase the detection signal for pathway relevance. In addition, they treat all gene alterations as equal, which is not expected to be valid for many biological systems.
- PGMs probabilistic graphical models
- approaches in probabilistic graphical models have been developed for learning causal networks compatible with multiple levels of observations. Efficient algorithms are available to learn pathways automatically from data (Friedman N, Goldszmidt M. (1997) Sequential Update of Bayesian Network Structure. In: Proceedings of the Thirteenth Conference on Uncertainty in Artificial Intelligence (UAI'97), Morgan Kaufmann Publishers, pp. 165-174; Murphy K, Weiss Y. Loopy belief propagation for approximate inference: An empirical study. In: Proceedings of Uncertainty in AI. 1999) and are well adapted to problems in genetic network inference (Friedman N. Inferring cellular networks using probabilistic graphical models.
- Breast cancer is clinically and genomically heterogeneous and is composed of several pathologically and molecular!y distinct subtypes.
- Patient responses to conventional and targeted therapeutics differ among subtypes motivating the development of marker guided therapeutic strategies.
- Collections of breast cancer cell lines mirror many of the molecular subtypes and pathways found in tumors, suggesting that treatment of cell lines with candidate therapeutic compounds can guide identification of associations between molecular subtypes, pathways and drug response.
- nearly all drugs show differential responses across these cell lines and approximately half show subtype-, pathway and/or genomic aberration-specific responses.
- the MicMa set has been used in parallel pilot studies of whole genome mRNA expression (1 Naume, B. et al., (2007), Presence of bone marrow micrometastasis is associated with different recurrence risk within molecular subtypes of breast cancer, 1 : 160-171), arrayCGH (Russnes HG, Vollan HKM, Lingjaerde OC, Krasnitz A, Lundin P, Naume B, S0rlie T, Borgen E, Rye 1H, Langer0d A, Chin S, Teschendorff AE, Stephens PJ, Miner S, Schlichting E, Baumbusch LO, Karesen R, Stratton MP, Wigler M, Caldas C, Zetterberg A, Hicks J, B0rresen-Dale A.
- Genomic architecture characterizes tumor progression paths and fate in breast cancer patients. Sci Transl Med 2010 Jun.;2(38):38ra47), DNA methylation (R0nneberg JA, Fleischer T, Solvang HK, Nordgard SH, Edvardsen H, Potapenko I, Nebdal D, Daviaud C, Gut I, Bukholm I, Naume B, B0rresen-Dale A, Tost J, Kristensen V.
- Methylation profiling with a panel of cancer related genes association with estrogen receptor, TP53 mutation status and expression subtypes in sporadic breast cancer. Mol Oncol 2011 Feb.;5( l):61-76), whole genome S P and S P-CGH (Van, Loo P. et al., (2010), Allele-specific copy number analysis of tumors, 107: 16910- 169154), whole genome miRNA expression analyses (Enerly, E. et al., (2011 ), miRNA-mRNA Integrated Analysis Reveals Roles for miRNAs in Primary Breast Tumors, 6: el6915-), TP53 mutation status dependent pathways and high throughput paired end sequencing (Stephens, P. J.
- a topic of great importance in cancer research is the identification of genomic aberrations that drive the development of cancer.
- Utilizing whole-genome copy number and expression profiles from the MicMa cohort we defined several filtering steps, each designed to identify the most promising candidates among the genes selected in the previous step.
- the first two steps involve identification of commonly aberrant and in-cis correlated to expression genes, i.e. genes for which copy number changes have substantial effect on expression.
- the method considers in-trans effects of the selected genes to further narrow down the potential novel candidate driver genes (Miriam Ragle Aure, Israel Steinfeld Lars Oliver Baumbusch nut Liest0l Doron Lipson Bj0rn Naume Vessela N. Kristensen Anne-Lise B0rresen-Dale Ole-Christian Lingjserde and Zohar Yakhini, (2011), A robust novel method for the integrated analysis of copy number and expression reveals new candidate driver genes in breast cancer).
- the inventors have discovered various systems and methods of pathway analysis that allow for integration of multiple attributes of multiple pathway elements (typically of one or more pathways) where at least one pathway element has an a priori known attribute, where at least another pathway element has an assumed attribute, where the pathway elements are cross-correlated and assigned specific influence levels on at least one pathway to so construct a probabilistic pathway model (PPM).
- PPM probabilistic pathway model
- Measured attributes for multiple elements of a patient sample are then used in conjunction with the PPM to so produce a patient sample specific dynamic pathway map (DPM).
- DPM patient sample specific dynamic pathway map
- the inventors contemplate a method of generating a dynamic pathway map (DPM) in which in one step access is provided to a pathway element database that stores a plurality of pathway elements, wherein each pathway element is characterized by its involvement in at least one pathway.
- access is provided to a modification engine that is coupled to the pathway element database, and the modification engine is used to associate a first pathway element with at least one a priori known attribute.
- DPM dynamic pathway map
- the modification engine is used to associate a second pathway element with at least one assumed attribute, and in yet a further step, the modification engine is used to cross-correlate and assign an influence level of the first and second pathway elements for at least one pathway using the known and assumed attributes, respectively, to thereby form a probabilistic pathway model.
- the probabilistic pathway model is used, via an analysis engine, to derive from a plurality of measured attributes for a plurality of elements of a patient sample the DPM having reference pathway activity information for a particular pathway.
- the pathway is within a regulatory pathway network, and particularly contemplated regulatory pathway networks include an ageing pathway network, an apoptosis pathway network, a homeostasis pathway network, a metabolic pathway network, a replication pathway network, and an immune response pathway network.
- the pathway may also be within a signaling pathway network and/or within a network of distinct pathway networks.
- suitable signaling pathway networks include
- calcium/calmodulin dependent signaling pathway network a cytokine mediated signaling pathway network, a chemokine mediated signaling pathway network, a growth factor signaling pathway network, a hormone signaling pathway network, a MAP kinase signaling pathway network, a phosphatase mediated signaling pathway network, a Ras superfamily mediated signaling pathway network, and a transcription factor mediated signaling pathway network.
- preferred pathway elements are proteins.
- preferred proteins include a receptor, a hormone binding protein, a kinase, a transcription factor, a methylase, a histone acetylase, and a histone deacetylase.
- preferred pathway elements are nucleic acids, such nucleic acids will typically include a protein coding sequence, a genomic regulatory sequence, a regulatory RNA, and a trans-activating sequence.
- the reference pathway activity information is specific with respect to a normal tissue, a diseased tissue, an ageing tissue, and/or a recovering tissue.
- Known and assumed attributes are typically and independently a compound attribute, a class attribute, a gene copy number, a transcription level, a translation level, or a protein activity, while the measured attributes are preferably a mutation, a differential genetic sequence object, a gene copy number, a transcription level, a translation level, a protein activity, and/or a protein interaction.
- the inventors contemplate a method of generating a dynamic pathway map (DPM) in which in one step access to a model database is provided that stores a probabilistic pathway model that comprises a plurality of pathway elements.
- DPM dynamic pathway map
- a first number of the plurality of pathway elements are cross-correlated and assigned an influence level for at least one pathway on the basis of known attributes
- a second number of the plurality of pathway elements are cross-correlated and assigned an influence level for at least one pathway on the basis of assumed attributes.
- a plurality of measured attributes for a plurality of elements of a patient sample is used, via an analysis engine, to modify the probabilistic pathway model to so obtain the DPM, wherein the DPM has reference pathway activity information for a particular pathway.
- the pathway is within a regulatory pathway network, a signaling pathway network, and/or a network of distinct pathway networks, and/or that the pathway element is a protein (for example, a receptor, a hormone binding protein, a kinase, a transcription factor, a methylase, a histone acetylase, a histone deacetylase, etc.) or a nucleic acid (for example, a genomic regulatory sequence, regulatory RNA, trans-activating sequence, etc.).
- the reference pathway activity information the known attribute, the assumed attribute, and the measured attribute, the same considerations as outlined above apply.
- a method of analyzing biologically relevant information may include a step of providing access to a model database that stores a dynamic pathway map (DPM), wherein the DPM is generated by modification of a probabilistic pathway model with a plurality of measured attributes for a plurality of elements of a first cell or patient sample.
- DPM dynamic pathway map
- a plurality of measured attributes for a plurality of elements of a second cell or patient sample is obtained, and the DPM and the plurality of measured attributes for the plurality of elements of the second cell or patient sample are used, via an analysis engine, to determine a predicted pathway activity information for the second cell or patient sample.
- the measured attributes for the plurality of elements of the first cell or patient sample are characteristic for a healthy cell or tissue, a specific age of a cell or tissue, a specific disease of a cell or tissue, a specific disease stage of a diseased cell or tissue, a specific gender, a specific ethnic group, a specific occupational group, and/or a specific species.
- the measured attributes for the plurality of elements of the second cell or patient sample will include information about a mutation, a differential genetic sequence object, a gene copy number, a transcription level, a translation level, a protein activity, and/or a protein interaction.
- the first and second samples are obtained from the same cell or patient, and it should be appreciated that treatment (for example, radiation, administration of a pharmaceutical) may be provided to the cell or patient before obtaining the plurality of measured attributes for the plurality of elements of the second cell or patient sample.
- treatment for example, radiation, administration of a pharmaceutical
- the treatment includes administration of a candidate molecule to the cell (for example, where the candidate molecule is a member of a library of candidate molecule).
- the predicted pathway activity information identifies an element as a hierarchical- dominant element in at least one pathway, and/or identifies the element as a disease-determinant element in at least one pathway with respect to a disease.
- a graphical representation of predicted pathway activity information may be provided, and/or a treatment recommendation may be generated that is at least in part based on the predicted pathway activity information.
- the predicted pathway activity information may be used to formulate a diagnosis, a prognosis for a disease, or a recommendation selected from the group consisting of a selection of a treatment option, and/or a dietary guidance, or to identify an epigenetic factor, a stress adaptation, a state of an organism, and/or a state of repair or healing.
- the invention provides a method of generating a dynamic pathway map (DPM), the method comprising: providing access to a pathway element database storing a plurality of pathway elements, each pathway element being characterized by its involvement in at least one pathway; providing access to a modification engine coupled to the pathway element database; using the modification engine to associate a first pathway element with at least one a priori known attribute; using the modification engine to associate a second pathway element with at least one assumed attribute; using the modification engine to cross-correlate and assign an influence level of the first and second pathway elements for at least one pathway using the known and assumed attributes, respectively, to form a probabilistic pathway model; and using the probabilistic pathway model, via an analysis engine, to derive from a plurality of measured attributes for a plurality of elements of a patient sample the DPM having reference pathway activity information for a particular pathway.
- DPM dynamic pathway map
- the pathway element is a protein.
- the protein is selected from the group consisting of a receptor, a hormone binding protein, a kinase, a transcription factor, a methylase, a histone acetylase, and a histone deacetylase.
- the pathway element is a nucleic acid.
- the nucleic acid is selected from the group consisting of a protein coding sequence, a genomic regulatory sequence, a regulatory RNA, and a irarcs-activating sequence.
- the reference pathway activity information is specific with respect to a normal tissue, a diseased tissue, an ageing tissue, or a recovering tissue.
- the known attribute is selected from the group consisting of a compound attribute, a class attribute, a gene copy number, a transcription level, a translation level, and a protein activity.
- the assumed attribute is selected from the group consisting of a compound attribute, a class attribute, a gene copy number, a transcription level, a translation level, and a protein activity.
- the measured attributes are selected from the group consisting of a mutation, a differential genetic sequence object, a gene copy number, a transcription level, a translation level, a protein activity, and a protein interaction.
- the pathway is within a regulatory pathway network.
- the regulatory pathway network is selected from the group consisting of an ageing pathway network, an apoptosis pathway network, a homeostasis pathway network, a metabolic pathway network, a replication pathway network, and an immune response pathway network.
- the pathway is within a signaling pathway network.
- the pathway is within a network of distinct pathway networks.
- the signaling pathway network is selected from the group consisting of a calcium/calmodulin dependent signaling pathway network, a cytokine mediated signaling pathway network, a chemokine mediated signaling pathway network, a growth factor signaling pathway network, a hormone signaling pathway network, a MAP kinase signaling pathway network, a phosphatase mediated signaling pathway network, a Ras superfamily mediated signaling pathway network, and a transcription factor mediated signaling pathway network.
- the invention also provides a method of generating a dynamic pathway map (DPM), the method comprising: providing access to a model database that stores a probabilistic pathway model that comprises a plurality of pathway elements; wherein a first number of the plurality of pathway elements are cross-correlated and assigned an influence level for at least one pathway on the basis of known attributes; wherein a second number of the plurality of pathway elements are cross-correlated and assigned an influence level for at least one pathway on the basis of assumed attributes; and using a plurality of measured attributes for a plurality of elements of a patient sample, via an analysis engine, to modify the probabilistic pathway model to obtain the DPM, wherein the DPM has reference pathway activity information for a particular pathway.
- DPM dynamic pathway map
- the pathway is within a regulatory pathway network, a signaling pathway network, or a network of distinct pathway networks.
- the pathway element is a protein selected from the group consisting of a receptor, a hormone binding protein, a kinase, a transcription factor, a methylase, a histone acetylase, and a histone deacetylase or a nucleic acid is selected from the group consisting of a genomic regulatory sequence, a regulatory RNA, and a /rani-activating sequence.
- the reference pathway activity information is specific with respect to a normal tissue, a diseased tissue, an ageing tissue, or a recovering tissue.
- the known attribute is selected from the group consisting of a compound attribute, a class attribute, a gene copy number, a transcription level, a translation level, and a protein activity.
- the assumed attribute is selected from the group consisting of a compound attribute, a class attribute, a gene copy number, a transcription level, a translation level, and a protein activity.
- the measured attributes are selected from the group consisting of a mutation, a differential genetic sequence object, a gene copy number, a transcription level, a translation level, a protein activity, and a protein interaction.
- the invention further provides a method of analyzing biologically relevant information, comprising: providing access to a model database that stores a dynamic pathway map (DPM), wherein the DPM is generated by modification of a probabilistic pathway model with a plurality of measured attributes for a plurality of elements of a first cell or patient sample; obtaining a plurality of measured attributes for a plurality of elements of a second cell or patient sample; and using the DPM and the plurality of measured attributes for the plurality of elements of the second cell or patient sample, via an analysis engine, to determine a predicted pathway activity information for the second cell or patient sample.
- DPM dynamic pathway map
- the measured attributes for the plurality of elements of the first cell or patient sample are characteristic for a healthy cell or tissue, a specific age of a cell or tissue, a specific disease of a cell or tissue, a specific disease stage of a diseased cell or tissue, a specific gender, a specific ethnic group, a specific occupational group, and a specific species.
- the measured attributes for the plurality of elements of the second cell or patient sample are selected from the group consisting of a mutation, a differential genetic sequence object, a gene copy number, a transcription level, a translation level, a protein activity, and a protein interaction.
- the first and second samples are obtained from the same cell or patient, and further comprising providing a treatment to the cell or patient before obtaining the plurality of measured attributes for the plurality of elements of the second cell or patient sample.
- the treatment is selected from the group consisting of radiation, administration of a pharmaceutical to the patient, and administration of a candidate molecule to the cell.
- the candidate molecule is a member of a library of candidate molecules.
- the predicted pathway activity information identifies an element as a hierarchical-dominant element in at least one pathway.
- the predicted pathway activity information identifies an element as a disease-determinant element in at least one pathway with respect to a disease.
- the method further comprises a step of generating a graphical representation of predicted pathway activity information.
- the method further comprises a step of generating a treatment recommendation that is at least in part based on the predicted pathway activity information.
- the method further comprises a step of using the predicted pathway activity information to formulate a diagnosis, a prognosis for a disease, or a recommendation selected from the group consisting of a selection of a treatment option, and a dietary guidance.
- the method further comprises a step of using the predicted pathway activity information to identify an epigenetic factor, a stress adaptation, a state of an organism, and a state of repair or healing.
- the invention provides a transformation method for creating a matrix of integrated pathway activities (IPAs) for predicting a clinical outcome for an individual in need, the method comprising the steps of (i) providing a set of curated pathways, wherein the pathways comprise a plurality of entities; (ii) converting each curated pathway into a distinct probabilistic graphical model (PGM), wherein the PGM is derived from factor graphs of each curated pathway, (iii) providing a biological sample from the individual wherein the biological sample comprises at least one endogenous entity comprised in one of the curated pathways; (iv) determining the levels of endogenous entity in the biological sample; (v) comparing the levels of the endogenous entity with those levels of the entity in a previously determined control sample from another individual; (vi) determining whether the levels of the endogenous entity relative to the control entity levels are activated, nominal, or inactivated; (vii) assigning the endogenous entity a numeric state, wherein the state representing activated is +1
- the method for creating a matrix of IPAs comprises predicting a clinical outcome, providing a diagnosis, providing a treatment, delivering a treatment, administering a treatment, conducting a treatment, managing a treatment, or dispensing a treatment to an individual in need.
- the set of curated pathways is from an analysis of human biology.
- the set of curated pathways is from an analysis of non-human biology.
- the determining of the levels of the endogenous entity relative to the control entity levels is performed using Student's r-test.
- the determining of the levels of the endogenous entity relative to the control entity levels is performed using ANOVA.
- the transforming method comprise the steps of wherein a plurality of matrices of integrated pathway activities from more than one individual are combined, the combined plurality of matrices resulting in a cluster, and where the distances between the individuals' matrices of the resulting cluster are determined.
- the determined distances are analysed using K-means cluster analysis.
- the determined distances are analysed using 2 -means cluster analysis.
- the transforming method comprises the step of determining the levels of endogenous entity in the biological sample comprises detecting the endogenous entity with an antibody and thereby determining the levels of endogenous entity.
- the step of determining the levels of endogenous entity in the biological sample comprises detecting the endogenous entity with a nucleic acid probe and thereby determining the levels of endogenous entity.
- the step of determining the levels of endogenous entity in the biological sample comprises detecting the endogenous entity with an organic reagent, wherein the organic reagent binds to the endogenous entity thereby resulting in a detectable signal and thereby determining the levels of endogenous entity.
- the step of determining the levels of endogenous entity in the biological sample comprises detecting the endogenous entity with an inorganic reagent, wherein the inorganic reagent binds to the endogenous entity thereby resulting in a detectable signal and thereby determining the levels of endogenous entity.
- the step of determining the levels of endogenous entity in the biological sample comprises detecting the endogenous entity with an organic reagent, wherein the organic reagent reacts with the endogenous entity thereby resulting in a detectable signal and thereby determining the levels of endogenous entity.
- the step of determining the levels of endogenous entity in the biological sample comprises detecting the endogenous entity with an inorganic reagent, wherein the inorganic reagent reacts with the endogenous entity thereby resulting in a detectable signal and thereby determining the levels of endogenous entity.
- the step of determining the levels of endogenous entity in the biological sample comprises measuring the absorbance of the endogenous entity at the optimal wavelength for the endogenous entity and thereby determining the levels of endogenous entity.
- the step of determining the levels of endogenous entity in the biological sample comprises measuring the fluorescence of the endogenous entity at the optimal wavelength for the endogenous entity and thereby determining the levels of endogenous entity.
- the step of determining the levels of endogenous entity in the biological sample comprises reacting the endogenous entity with an enzyme, wherein the enzyme selectively digests the endogenous entity to create at least one product, detecting the at least one product, and thereby determining the levels of endogenous entity.
- the step of reacting the endogenous entity with an enzyme results in creating at least two products.
- the step of reacting the endogenous entity with an enzyme resulting at least two products is followed by a step of treating the products with another enzyme, wherein the enzyme selectively digests at least one of the products to create at least a third product, detecting the at least a third product, and thereby determining the levels of endogenous entity.
- the individual is selected from the group of a healthy individual, an asymptomatic individual, and a symptomatic individual.
- the individual is selected from the group consisting of an individual diagnosed with a condition, the condition selected from the group consisting of a disease and a disorder.
- the condition is selected from the group consisting of acquired immunodeficiency syndrome (AIDS), Addison's disease, adult respiratory distress syndrome, allergies, ankylosing spondylitis, amyloidosis, anemia, asthma, atherosclerosis, autoimmune hemolytic anemia, autoimmune thyroiditis, benign prostatic hyperplasia, bronchitis, Chediak-Higashi syndrome, cholecystitis, Crohn's disease, atopic dermatitis, dermnatomyositis, diabetes mellitus, emphysema,
- AIDS acquired immunodeficiency syndrome
- Addison's disease AIDS
- adult respiratory distress syndrome e.g., allergies, ankylosing spondylitis, amyloidosis, anemia, asthma, atherosclerosis, autoimmune hemolytic anemia, autoimmune thyroiditis, benign prostatic hyperplasia, bronchitis, Chediak-Higashi syndrome, cholecystitis, Crohn
- erythroblastosis fetalis, erythema nodosum, atrophic gastritis, glomerulonephritis, Goodpasture's syndrome, gout, chronic granulomatous diseases, Graves' disease, Hashimoto's thyroiditis, hypereosinophilia, irritable bowel syndrome, multiple sclerosis, myasthenia gravis, myocardial or pericardial inflammation, osteoarthritis, osteoporosis, pancreatitis, polycystic ovary syndrome, polymyositis, psoriasis, Reiter's syndrome, rheumatoid arthritis, scleroderma, severe combined immunodeficiency disease (SCID), Sjogren's syndrome, systemic anaphylaxis, systemic lupus erythematosus, systemic sclerosis, thrombocytopenic purpura, ulcerative colitis, uveitis, Werner syndrome, complications
- the condition is selected from the group consisting of cancers such as adenocarcinoma, leukemia, lymphoma, melanoma, myeloma, sarcoma, teratocarcinoma, and, in particular, cancers of the adrenal gland, bladder, bone, bone marrow, brain, breast, cervix, gall bladder, ganglia, gastrointestinal tract, heart, kidney, liver, lung, muscle, ovary, pancreas, parathyroid, penis, prostate, salivary glands, skin, spleen, testis, thymus, thyroid, and uterus; immune disorders such as acquired immunodeficiency syndrome (AIDS), Addison's disease, adult respiratory distress syndrome, allergies, ankylosing spondylitis, amyloidosis, anemia, asthma, atherosclerosis, autoimmune hemolytic anemia, autoimmune thyroiditis, bronchitis, cholecystitis, contact dermatitis, Crohn'
- cancers such
- the condition is selected from the group consisting of endocrinological disorders such as disorders associated with hypopituitarism including hypogonadism, Sheehan syndrome, diabetes insipidus, Kallman's disease, Hand-Schuller-Christian disease, Letterer-Siwe disease, sarcoidosis, empty sella syndrome, and dwarfism; hyperpituitarism including acromegaly, giantism, and syndrome of inappropriate antidiuretic hormone (ADH) secretion (SIADH); and disorders associated with hypothyroidism including goiter, myxedema, acute thyroiditis associated with bacterial infection, subacute thyroiditis associated with viral infection, autoimmune thyroiditis (Hashimoto's disease), and cretinism; disorders associated with hyperthyroidism including thyrotoxicosis and its various forms, Grave's disease, pretibial myxedema, toxic multinodular goiter, thyroid carcinoma, and Plummer's disease; and disorders associated with hyper
- the Invention also provides the transforming method as disclosed herein wherein matrix A can then be used in place of the original constituent datasets to identify associations with clinical outcomes.
- the curated pathways are selected from the group consisting of biochemical pathways, genetic pathways, metabolic pathways, gene regulatory pathways, gene transcription pathways, gene translation pathways.
- the entities are selected from the group consisting of nucleic acids, peptides, proteins, peptide nucleic acids, carbohydrates, lipids, proteoglycans, factors, co-factors, biochemical metabolites, organic compositions, inorganic compositions, and salts.
- the biological sample is selected from the group consisting of patient samples, control samples, experimentally-treated animal samples, experimentally-treated tissue culture samples, experimentally-treated cell culture samples, and experimentally-treated in vitro biochemical composition samples.
- the biological sample is a patient sample.
- the invention also provides a probabilistic graphical model (PGM) framework having an output that infers the molecular pathways altered in a patient sample, the PGM comprising a plurality of factor graphs, wherein the factor graphs represent integrated biological datasets, and wherein the inferred molecular pathways that are altered in a patient sample comprise molecular pathways known from data and wherein said molecular pathways effect a clinical or non-clinical condition, wherein the inferred molecular pathways are known to be modulated by a clinical regimen or treatment, and wherein the output indicates a clinical regimen.
- the data is selected from experimental data, clinical data, epidemiological data, and
- the condition is selected from the group consisting of a disease and a disorder.
- the condition is selected from the group consisting of acquired immunodeficiency syndrome (AIDS), Addison's disease, adult respiratory distress syndrome, allergies, ankylosing spondylitis, amyloidosis, anemia, asthma, atherosclerosis, autoimmune hemolytic anemia, autoimmune thyroiditis, benign prostatic hyperplasia, bronchitis, Chediak-Higashi syndrome, cholecystitis, Crohn's disease, atopic dermatitis, dermnatomyositis, diabetes mellitus, emphysema, erythroblastosis fetalis, erythema nodosum, atrophic gastritis, glomerulonephritis, Goodpasture's syndrome, gout, chronic granulomatous diseases, Graves' disease, Hashimoto's thyroiditis
- AIDS acquired immunodeficiency syndrome
- the condition is selected from the group consisting of cancers such as adenocarcinoma, leukemia, lymphoma, melanoma, myeloma, sarcoma, teratocarcinoma, and, in particular, cancers of the adrenal gland, bladder, bone, bone marrow, brain, breast, cervix, gall bladder, ganglia, gastrointestinal tract, heart, kidney, liver, lung, muscle, ovary, pancreas, parathyroid, penis, prostate, salivary glands, skin, spleen, testis, thymus, thyroid, and uterus; immune disorders such as acquired immunodeficiency syndrome (AIDS), Addison's disease, adult respiratory distress syndrome, allergies, ankylosing spondylitis, amyloidosis, anemia, asthma, atherosclerosis, autoimmune hemolytic anemia, autoimmune thyroiditis, bronchitis, cholecystitis, contact dermatitis, Crohn
- cancers such as
- the condition is selected from the group consisting of endocrinological disorders such as disorders associated with hypopituitarism including hypogonadism, Sheehan syndrome, diabetes insipidus, Kallman's disease, Hand-Schuller-Christian disease, Letterer-Siwe disease, sarcoidosis, empty sella syndrome, and dwarfism; hyperpituitarism including acromegaly, giantism, and syndrome of inappropriate antidiuretic hormone (ADH) secretion (SIADH); and disorders associated with hypothyroidism including goiter, myxedema, acute thyroiditis associated with bacterial infection, subacute thyroiditis associated with viral infection, autoimmune thyroiditis (Hashimoto's disease), and cretinism; disorders associated with hyperthyroidism including thyrotoxicosis and its various forms, Grave's disease, pretibial myxedema, toxic multinodular goiter, thyroid carcinoma, and Plummer's disease; and
- FIG. 1 illustrates an overview of the PARADIGM method.
- PARADIGM uses a pathway schematic with functional genomic data to infer genetic activities that can be used for further downstream analysis.
- FIG. 2 illustrates the conversion of a genetic pathway diagram into a PARADIGM model.
- PARADIGM uses a pathway schematic with functional genomic data to infer genetic activities that can be used for further downstream analysis.
- A. Data on a single patient is integrated for a single gene using a set of four different biological entities for the gene describing the DNA copies, mRNA and protein levels, and activity of the protein.
- B. PARADIGM models various types of interactions across genes including transcription factors to targets (upper-left), subunits aggregating in a complex (upper- right), post- translational modification (lower-left), and sets of genes in a family performing redundant functions (lower-right).
- C. Toy example of a small sub-pathway involving P53, an inhibitor MDM2, and the high level process, apoptosis as represented in the model.
- FIG 3 illustrates exemplary NCI pathway interactions in The Cancer Genome Atlas (TCGA) project (http://cancergenome.nih.gov) glioblastoma multiform (GMB) data.
- TCGA Cancer Genome Atlas
- GMB glioblastoma multiform
- a histogram of correlations between randomly paired genes is shown for C2E (dashed red) and E2E (dashed blue). Arrows point to the enrichment of positive correlations found for the C2E (red) and E2E (blue) correlation.
- Figure 4 illustarates exemplary learning parameters for the anti-apoptotic serine-threonine kinase 1 (AKTI).
- IPAs Integrated Pathway Activities
- EM Expectation-Maximization
- Dots show IPAs from permuted samples and circles show IPAs fdrom real samples.
- the red line denotes the mean IPA in real samples and the green line denotes the man ⁇ of null samples.
- Figure 5 illustrates distinguishing decoy from real pathways with PARADIGM and Signaling Pathway Impact Analysis (SPIA). Decoy pathways werecreated by assigning a new gene name to each gene in a pathway. PARADIGM and SPIA were then used to compute the perturbation of every pathway. Each line shows the receiver-operator characteristic for distinguishing real from decoy pathways using the perturbation ranking.
- AUCs areas under the curve
- GBM glioblastoma multiform
- the AUCs are 0.642 and 0.604, respectively.
- Figure 6 illustrates exemplary patient sample IPAs compared with within permutations for Class I phosphatidylinositol-3-kinase (PI3K) signaling events mediated by Akt in breast cancer.
- PI3K phosphatidylinositol-3-kinase
- Biological entities were sorted by mean EPA in the patient samples (red) and compared with the mean IPA for the peruted samples.
- the colored areas around each mean denote the standard deviation (SD) of each set.
- the IPAs of the right include AKTI , CHUK, and MDM2.
- Figure 7 illustrates an exemplary CIRCLEMAP display of the ErbB2 pathway.
- estrogen receptor (ER) status For each node, estrogen receptor (ER) status, IPAs, expression data, and copy-number data are displayed as concentric circles, from innermost to outermost respectively.
- the apoptosis node and the ErbB2/ErbB3/neuregulin 2 complex node have circles only for ER status and for IPAs, as there are no direct observations of these entities.
- Each patient's data is displayed along one angle from the circle center to edge.
- Figure 8 illustarates exemplary clustering of IPAs for TCGA GBM.Each column corresponds to a single sample, and each row to a biomolecular entity. Color bars beneath the hierarchical clustering tree denote clusters used for Figure 9.
- Figure 9 illustrates Kaplan-Meier survival plots for the clusters from Figure 8.
- Figure 10 illustrates that cell lines show a broad range of responses to therapeutic compounds.
- A Luminal and ERBB2AMP cell lines preferentially respond to AKT inhibition. Each bar represents the response of a single breast cancer cell line to the Sigma AKT1-2 inhibitor. Cell lines are ordered by increasing sensitivity (-logio(GI 50 )) and colored according to subtype.
- B GI50 values for compounds with similar mechanisms are highly correlated. Heatmap shows hierarchical clustering of correlations between responses breast cancer cell lines treated with various compounds.
- C. Compounds with similar modes of action show similar patterns of response across the panel of cell lines. Each column represents one cell line, each row represents a compound tested. GI50 values are hierarchically clustered. Only compounds with a significant subtype effect are included.
- CNAs are associated sensitivity. Boxplots show distribution of response sensitivity for cell lines with aberrant (A) and normal (N) copy number at the noted genomic locus. FDR p values for the association between drug response and CNA are noted, a. 9p21 (CDKN2A) deletion is associated with response to ixabepilone, vinerolbine and fascaplysin. b. 20ql3 (STK15/AURKA) amplification is associated with VX-680 and GSK1070916. c. Amplification at 1 lq 13 (CCND1) is associated with response to carboplatin and GSK1070916.
- Figure 11 shows a heatmap of non-redundant PARADIGM activities both cell line and TCGA samples.
- Cluster dendrogram represents Euclidian distance between samples and was created using Eisen Cluster and drawn using Java Treeview. Colored bars below dendrogram represent sample subtype (top) and sample cohort (bottom).
- FIG 12 illustrates that cell line subtypes have unique network features.
- each node in the graph represents a different pathway "concept" corresponding to either a protein (circles), a multimeric complex (hexagons), or a an abstract cellular process (squares).
- the size of the nodes were drawn in proportion to the differential activity score such that larger nodes correspond to pathway concepts with activities more correlated with basal versus non-basal cell lines. Color indicates whether the concept is positively correlated (red) or negatively correlated (blue) with the basal subtype.
- Links represent different interactions including protein- protein level interactions (dashed lines) and transcriptional (solid lines).
- the MYC/MAX and ERK1/2 subnet is preferentially activated in basal breast cancer cell lines.
- the CTTNB l network is activated in claudin-low cell lines.
- a FOXA 1/FOXA2 network is upregulated in the luminal subtype.
- the ERBB2AMP subtype shows down-regulation of the RPS6KB 1 pathway.
- Figure 13 Illustrates how pathway diagrams can be used to predict response to therapies.
- Figure 14 illustrates exemplary genomic and transcriptional profiles of the breast cancer cell lines.
- A. DNA copy number aberrations for 43 breast cancer cell lines are plotted with log l0 (FDR) of GISTIC analysis on the y-axis and chromosome position on the x-axis. Copy number gains are shown in red with positive log, 0 (FDR) and losses are shown in green with negative log, 0 (FDR).
- Figure 15 illustrates that GI50 calculations are highly reproducible.
- A Each bar a count of the frequency of replicated drug/cell line combinations. Most cell lines were tested only one time against a particular compound, but some drug/cell line combinations were tested multiple times.
- B Each boxplot represents the distribution of median average deviations for drug/cell line pairs with 3 or 4 replicates.
- Figure 16 shows that doubling time varies across cell line subtype.
- A. Growth rate, computed as the median doubling time in hours, of the breast cancer cell lines subtypes are shown as box-plots. The basal and claudin-low subtypes have shorter median doubling time as compared to lumi nal and ERBB2 AMP subtypes, ruskal-Wallis p value (p 0.006).
- B. The ANCOVA model shows strong effects of both subtype and growth rate on response to 5'FU. Luminal (black) and basal/claudin-low (red) breast cancer lines each show significant associations to growth rate but have distinct slopes.
- Figure 17 shows that inferred pathway activities are more strongly correlated within subtypes than within cohorts. Shown is a histogram of t-statistics derived from Pearson correlations computed between cell lines and TCGA samples of the same subtype (red) compared to t-statistics of Pearson correlations between cell lines of different subtypes (black).
- X-axis corresponds to the Pearson correlation t-statistic;
- y-axis shows the density of (cell-line, cell-line) or (cell-line, TCGA sample) pairs.
- K-S test P ⁇ l x lO -22 ) indicates cell lines and TCGA samples of the same subtype are more alike than cell lines of other subtypes.
- Supplementary Figures 18-21 illustrate an exemplary network architecture for each of the four subnetworks identified from the SuperPathway.
- Figure 18 illustrates a network diagram of basal pathway markers.
- Each node in the graph represents a different pathway "concept" corresponding to either a protein (circles), a multimeric complex (hexagons), or a an abstract cellular process (squares).
- the size of the nodes are drawn in proportion to the differential activity score such that larger nodes correspond to pathway concepts with activities more correlated with basal versus non-basal cell lines. Color indicates whether the concept is positively correlated (red) or negatively correlated (blue) with the basal subtype.
- Links represent different interactions including protein-protein level interactions (dashed lines) and transcriptional (solid lines). Interactions were included in the map only if they interconnect concepts whose absolute level of differential activity is higher than the mean absolute level.
- Figure 19 illustrates an exemplary network diagram of claudin-low pathway markers. Convention as in Figure 18.
- Figure 20 illustrates an exemplary network diagram of luminal pathway markers. Convention as in Figure 18.
- Figure 21 illustrates an exemplary network diagram of ERBB2AMP pathway markers. Convention as in Figure 18.
- Figure 22 illustrates exemplary URKB-FOXMl -CCNB l networks in luminal, claudin-low and basal cell lines.
- A Network surrounding AURKB and FOXMl in luminal cell lines. CCNB l was not significantly downregulated and therefore does not appear on the pathway map.
- B In claudin-low cell lines, AURKB and FOXMl both up-regulated; activity for CCNB l was not significant.
- C AURKB, FOXMl and CCNB l are all up- regulated in basal cell lines. Convention as in Figure 18.
- Figure 23 illustrates an exemplary distribution of unsupervised clusters and survival curves of the patients of the MicMa cohort according to CNA, mRNA expression, DNA methylation and miRNA expression. For each type of genomic level the size of each cluster are plotted on the left, and to the right, survival curves are shown. Significance of differential survival are assessed by two methods (see Examples).
- Figure 24 illustrates an exemplary distribution of indentified PARADIGM clusters and survival.
- Figure 25 illustrates an exemplary heatmaps of Paradigm IPLs for each dataset. Each row shows the IPL of a gene or complex across all three cohorts. The colored bar across the top shows the MicMa-derived Paradigm clusters, as in Figure 2. Members of pathways of interest are labeled by their pathway. Red represents an activated IPL, blue a deactivated IPL.
- Figure 26 illustrates the FOXM1 Transcription Factor Network.
- the upper network diagram summarizes data from cluster pdgm.3, whereas the lower cluster summarizes the data from other clusters.
- Nodes shapes denote the data type that was most frequently perturbed within each cluster, and node color denote the direction of perturbation.
- Edge arrows denote the sign of interactions, and color denotes the type of interaction.
- Figure 27 illustrates a toy example of a small fragment of the p53 apoptosis pathway.
- a pathway diagram from NCI was converted into a factor graph that includes both hidden and observed states.
- Figure 28 illustrates an exemplary heatmap of Inferred Pathway Activities (IPAs).
- IP As representing 1598 inferences of molecular entities (rows) inferred to be activated (red) or inactivated (blue) are plotted for each of 316 patient tumor samples (columns).
- IPAs were hierarchically clustered by pathway entity and tumor sample, and labels on the right show sections of the heatmap enriched with entities of individual pathways.
- the colorbar legend is in log base 10.
- Figure 29 summarises FOXM1 integrated pathway activities (IPAs) across all samples.
- the arithmetic mean of IPAs across tumor samples for each entity in the FOXM1 transcription factor network is shown in red, with heavier red shading indicating two standard deviations. Gray line and shading indicates the mean and two standard deviations for IPAs derived from the 1000 "null" samples.
- Figure 30 shows a comparison of IPAs of FOXM1 to those of other tested transcription factors (TFs) in NCI Pathway Interaction Database.
- A Histogram of IPAs with non-active (zero- valued) IPAs removed. FOXM1 targets are significantly more activated than other NCI TFs (P ⁇ 10 '267 ; Kolmogorov-Smirnov ( S) test).
- B Histogram of all IPAs including non-active IPAs. Using all IPAs, FOXMl 's activity relative to other TFs is interpreted with somewhat higher significance (P ⁇ 10 '301 ; KS test).
- FIG 31 illustrates that FOXM1 is not expressed in fallopian epithelium compared to serous ovarian carcinoma.
- FOXMl 's expression levels in fallopian tube was compared to its levels in serous ovarian carcinoma using the data from Tone et al (PMTD: 18593983).
- FOXMl 's expression is much lower in fallopian tube, including in samples carrying BRCA 1/2 mutations, indicating that FOXMl 's elevated expression observed in the TCGA serous ovarian cancers is not simply due to an epithelial signature.
- Figure 32 shows expression of FOXM1 transcription factor network genes in high grade versus low grade carcinoma.
- Expression levels for FOXM1 and nine selected FOXM1 targets were plotted for both low-grade (I; tan boxes; 26 samples) and high-grade (II/III; blue boxes; 296 samples) ovarian carcinomas. Seven out of the nine targets were showed to have significantly high expression of FOXM 1 in the high-grade carcinomas (Student's t-test; p-values noted under boxplots).
- Figure 33 shows that the cell lines show a broad range of responses to therapeutic compounds.
- C. Compounds with similar modes of action show similar patterns of response across the panel of cell lines. Each column represents one cell line, each row represents a compound tested. GI50 values are hierarchically clustered.
- CNAs are associated sensitivity. Boxplots show distribution of response sensitivity for cell lines with aberrant (A) and normal (N) copy number at the noted genomic locus. FDR p values for the association between drug response and CNA are noted, a. 9p21 (CDK 2A) deletion is associated with response to ixabepilone, vinerolbine and fascaplysin. b. 20ql3 (STK15/AURKA) amplification is associated with VX-680 and GSK1070916. c. Amplification at 1 lql 3 (CCND1) is associated with response to carboplatin and GSK 1070916.
- FIG. 34 A. Heatmap of non-redundant PARADIGM activities both cell line and TCGA samples.
- Cluster dendrogram represents Euclidian distance between samples and was created using Eisen Cluster and drawn using Java Treeview. Colored bars below dendrogram represent sample subtype (top) and sample cohort (bottom).
- FIG 35 shows that the cell line subtypes have unique network features.
- each node in the graph represents a different pathway "concept" corresponding to either a protein (circles), a multimeric complex (hexagons), or a an abstract cellular process (squares).
- the size of the nodes were drawn in proportion to the differential activity score such that larger nodes correspond to pathway concepts with activities more correlated with basal versus non-basal cell lines. Color indicates whether the concept is positively correlated (red) or negatively correlated (blue) with the basal subtype.
- Links represent different interactions including protein- protein level interactions (dashed lines) and transcriptional (solid lines).
- the MYC/MAX and ERK1/2 subnet is preferentially activated in basal breast cancer cell lines.
- the CTTNB l network is activated in claudin-low cell lines.
- a FOXA1/FOXA2 network is upregulated in the luminal subtype.
- the ERBB2AMP subtype shows down-regulation of the RPS6KB 1 pathway.
- Figure 36 shows that the pathway diagrams can be used to predict response to therapies.
- Figure 37 illustrates genome copy number abnormalities, (a) Copy-number profiles of 489 HGS- OvCa, compared to profiles of 197 glioblastoma multiforme (GBM) tumors46. Copy number increases (red) and decreases (blue) are plotted as a function of distance along the normal genome, (b) Significant, focally amplified (red) and deleted (blue) regions are plotted along the gnome. Annotations include the 20 most significant amplified and deleted regions, well-localized regions with 8 or fewer genes, and regions with known cancer genes or genes identified by genome-wide loss-of-function screens. The number of genes included in each region is given in brackets, (c) Significantly amplified (red) and deleted (blue) chromosome arms.
- Figure 38 illustrates gene and miRNA expression patterns of molecular subtype and outcome prediction in HGS- OvCa.
- a prognostic gene signature was defined and applied to a test dataset.
- Figure 39 illustartes altered Pathways in HGS-OvCa.
- NOTCH pathway, identified by HotNet analysis are commonly altered. Alterations are defined by somatic mutations, DNA copy-number changes, or in some cases by significant up- or down- regulation compared to expression in diploid tumors. Alteration frequencies are in percentage of all cases;
- activated genes are red, inactivated genes are blue, (c) Genes in the HR pathway are altered in up to 49% of cases. Survival analysis of BRCA status shows divergent outcome for BRCA mutated cases (exhibiting better overall survival) than BRCA wild-type, and BRCA / epigenetically silenced cases exhibiting worse survival, (d) The FOX 1 transcription factor network is activated in 87% of cases.
- Each gene is depicted as a multi-ring circle in which its copy number (outer ring) and gene expression (inner ring) are plotted such that each "spoke" in the ring represents a single patient sample, with samples sorted in increasing order of FOXMl expression.
- Excitatory (red arrows) and inhibitory interactions (blue lines) were taken from the NCI Pathway Interaction Database. Dashed lines indicate transcriptional regulation.
- Fig. 40 is a schematic of an exemplary computer system to produce a dynamic pathway map according to the inventive subject matter.
- the term "curated” means the relationships between a set of biological molecules and/or non-biological molecules that has been tested, analyzed, and identified according to scientific and/or clinical principles using methods well known in the art, such as molecular biological, biochemical, physiological, anatomical, genomic, transcriptomic, proteomic, metabolomic, ADME, and bioinformatic techniques, and the like.
- the relationships may be biochemical such as biochemical pathways, genetic pathways, metabolic pathways, gene regulatory pathways, gene transcription pathways, gene translation pathways, miRNA-regulated pathways, pseudogene-regulated pathways, and the like.
- High-throughput data is providing a comprehensive view of the molecular changes in cancer tissues.
- New technologies allow for the simultaneous genome-wide assay of the state of genome copy number variation, gene expression, DNA methylation, and epigenetics of tumor samples and cancer cell lines.
- the inventors have developed systems and methods where multiple attributes of multiple pathway elements are integrated into a probabilistic pathway model that is then modified using patient data to produce a dynamic pathway map.
- the attributes for pathway elements within a pathway need not be known a priori. Indeed, at least some of the attributes of at least some pathway elements are assumed.
- the pathway elements are then cross-correlated and assigned specific influence levels on or more pathways to so construct the probabilistic pathway model, which is preferably representative of a particular reference state (for example, healthy or diseased).
- Measured attributes for multiple elements of a patient sample are then used in conjunction with the probabilistic pathway model to so produce a patient sample specific dynamic pathway map that provides reference pathway activity information for one or more particular pathways.
- computing devices comprise a processor configured to execute software instructions stored on a tangible, non-transitory computer readable storage medium (for example, hard drive, solid state drive, RAM, flash, ROM, etc.).
- the software instructions preferably configure the computing device to provide the roles, responsibilities, or other functionality as discussed below with respect to the disclosed apparatus.
- the various servers, systems, databases, or interfaces exchange data using standardized protocols or algorithms, possibly based on HTTP, HTTPS, AES, public-private key exchanges, web service APIs, known financial transaction protocols, or other electronic information exchanging methods.
- Data exchanges preferably are conducted over a packet-switched network, the Internet, LAN, WAN, VPN, or other type of packet switched network.
- inventive subject matter is considered to include all possible combinations of the disclosed elements.
- inventive subject matter is also considered to include other remaining combinations of A, B, C, or D, even if not explicitly disclosed.
- a gene is modeled by a factor graph as a set of interconnected variables encoding the expression and known activity of a gene and its products, allowing the incorporation of many types of -omic data as evidence.
- the method predicts the degree to which a pathway's activities (for example, internal gene states, interactions, or high-level ""outputs") are altered in the patient using probabilistic inference.
- a pathway's activities for example, internal gene states, interactions, or high-level ""outputs"
- SPIA a competing pathway activity inference approach
- our method identifies altered activities in cancer- related pathways with fewer false-positives in, but not limted to, both a glioblastoma multiform (GBM) and a breast cancer dataset.
- GBM glioblastoma multiform
- Pathway Recognition Algorithm using Data integration on Genomic Models identified consistent pathway-level activities for subsets of the GBM patients that are overlooked when genes are considered in isolation. Further, grouping GBM patients based on their significant pathway perturbations using the algorithm divides them into clinically-relevant subgroups having significantly different survival outcomes.
- PGM probabilistic graphical model
- PARADIGM PAthway Recognition Algorithm using Data Integration on Genomic Models
- Figure 1 illustrates the overview of the approach. Multiple genome-scale measurements on a single patient sample are combined to infer the activities of genes, products, and abstract process inputs and outputs for a single National Cancer Institute (NCI) pathway.
- PARADIGM produces a matrix of integrated pathway activities (IPAs) A where A u represents the inferred activity of entity i in patient sample j. The matrix A can then be used in place of the original constituent datasets to identify associations with clinical outcomes.
- IPAs integrated pathway activities
- the method may be used to provide clinical information that can be used in a variety of diagnostic and therapeutic applications, such as detection of cancer tissue, staging of cancer tissue, detection of metastatic tissue, and the like; detection of neurological disorders, such as, but not limited to, Alzheimer's disease, amyotrophic lateral sclerosis (ALS), Parkinson's disease, schizophrenia, epilepsy, and their complications; developmental disorders such as DiGeorge Syndrome, autism, autoimmune disorders such as multiple sclerosis, diabetes, and the like; treatment of an infection, such as, but not limited to, viral infection, bacterial infection, fungal infection, leishmania, schistosomiasis, malaria, tape-worm, elephantiasis, infections by nematodes, nematines, and the like.
- diagnostic and therapeutic applications such as detection of cancer tissue, staging of cancer tissue, detection of metastatic tissue, and the like
- detection of neurological disorders such as, but not limited to, Alzheimer's disease, amyotrophic lateral sclerosis (ALS), Parkinson's disease, schizophrenia, epilepsy, and
- the method may be used to provide clinical information to detect and quantify altered gene expression, absence/presence versus excess, expression of mRNAs or to monitor mRNA levels during therapeutic intervention.
- Conditions, diseases or disorders associated with altered expression include acquired immunodeficiency syndrome (AIDS), Addison's disease, adult respiratory distress syndrome, allergies, ankylosing spondylitis, amyloidosis, anemia, asthma, atherosclerosis, autoimmune hemolytic anemia, autoimmune thyroiditis, benign prostatic hyperplasia, bronchitis, Chediak-Higashi syndrome, cholecystitis, Crohn's disease, atopic dermatitis, dermnatomyositis, diabetes mellitus, emphysema, erythroblastosis fetalis, erythema nodosum, atrophic gastritis, glomerulonephritis, Goodpasture's syndrome, gout, chronic granulomatous diseases,
- AIDS acquired immuno
- the method may be used to provide clinical information to detect and quantify altered gene expression; absence, presence, or excess expression of mRNAs; or to monitor mRNA levels during therapeutic intervention.
- Disorders associated with altered expression include akathesia, Alzheimer's disease, amnesia, amyotrophic lateral sclerosis (ALS), ataxias, bipolar disorder, catatonia, cerebral palsy, cerebrovascular disease Creutzfeldt-Jakob disease, dementia, depression, Down's syndrome, tardive dyskinesia, dystonias, epilepsy, Huntington's disease, multiple sclerosis, muscular dystrophy, neuralgias, neurofibromatosis, neuropathies, Parkinson's disease, Pick's disease, retinitis pigmentosa, schizophrenia, seasonal affective disorder, senile dementia, stroke, Tourette's syndrome and cancers including adenocarcinomas, melanomas, and teratocarcinomas, particularly of the brain.
- the method may be used to provide clinical information for a condition associated with altered expression or activity of the mammalian protein.
- conditions include, but are not limited to, acquired immunodeficiency syndrome (AIDS), Addison's disease, adult respiratory distress syndrome, allergies, ankylosing spondylitis, amyloidosis, anemia, asthma, atherosclerosis, autoimmune hemolytic anemia, autoimmune thyroiditis, benign prostatic hyperplasia, bronchitis, Chediak-Higashi syndrome, cholecystitis, Crohn's disease, atopic dermatitis, dermatomyositis, diabetes mellitus, emphysema, erythroblastosis fetalis, erythema nodosum, atrophic gastritis, glomerulonephritis, Goodpasture's syndrome, gout, chronic granulomatous diseases, Graves' disease, Hashimoto's thyroiditis,
- AIDS acquired immunodefic
- the methods disclosed erein may be used to detect, stage, diagnose, and/or treat a disorder associated with decreased expression or activity of the nucleic acid sequences.
- disorders include, but are not limited to, cancers such as adenocarcinoma, leukemia, lymphoma, melanoma, myeloma, sarcoma, teratocarcinoma, and, in particular, cancers of the adrenal gland, bladder, bone, bone marrow, brain, breast, cervix, gall bladder, ganglia, gastrointestinal tract, heart, kidney, liver, lung, muscle, ovary, pancreas, parathyroid, penis, prostate, salivary glands, skin, spleen, testis, thymus, thyroid, and uterus; immune disorders such as acquired immunodeficiency syndrome (AIDS), Addison's disease, adult respiratory distress syndrome, allergies, ankylosing spondylitis, amyloidosis, anemia, asthma, at
- the methods disclosed erein may be used to detect, stage, diagnose, and/or treat a disorder associated with expression of the nucleic acid sequences.
- a disorder associated with expression of the nucleic acid sequences include, but are not limited to, endocrinological disorders such as disorders associated with hypopituitarism including hypogonadism, Sheehan syndrome, diabetes insipidus, Kallman's disease, Hand-Schuller-Christian disease, Letterer-Siwe disease, sarcoidosis, empty sella syndrome, and dwarfism; hyperpituitarism including acromegaly, giantism, and syndrome of inappropriate antidiuretic hormone (ADH) secretion (S1ADH); and disorders associated with hypothyroidism including goiter, myxedema, acute thyroiditis associated with bacterial infection, subacute thyroiditis associated with viral infection, autoimmune thyroiditis (Hashimoto's disease), and cretinism; disorders associated with hyperthyroidism including thyrotoxico
- the polynucleotide sequences may be used in Southern or Northern analysis, dot blot, or other membrane-based technologies; in PCR technologies; in dipstick, pin, and ELISA assays; and in microarrays utilizing fluids or tissues from patients to detect altered nucleic acid sequence expression. Such qualitative or quantitative methods are well known in the art.
- PARADIGM Inference of patient-specific pathway activities from multi-dimensional cancer genomics data using PARADIGM.
- pathway-based approaches the genetic interactions found in pathway databases carry information for interpreting correlations between gene expression changes detected in cancer. For example, if a cancer-related pathway includes a link from a transcriptional activator A to a target gene T, we expect the expression of A to be positively correlated with the expression of T (E2E correlation). Likewise, we also expect a positive correlation between A's copy number and T's expression (C2E correlation). Further, we expect C2E correlation to be weaker than E2E correlation because amplification in A does not necessarily imply A is expressed at higher levels, which in turn is necessary to upregulate B. In this way, each link in a pathway provides an expectation about the data; pathways with many consistent links may be relevant for further consideration.
- the PARADIGM method integrates diverse high-throughput genomics information with known signaling pathways to provide patient-specific genomic inferences on the state of gene activities, complexes, and cellular processes.
- the core of the method uses a factor graph to leverage inference for combining the various data sources.
- the use of such inferences in place of, or in conjunction with, the original high-throughput datasets improves our ability to classify samples into clinically relevant subtypes.
- Clustering the GBM patients based on the PARADIGM-integrated activities revealed patient subtypes correlated with different survival profiles. In contrast, clustering the samples either using the expression data or the copy-number data did not reveal any significant clusters in the dataset.
- PARADIGM produces pathway inferences of significantly altered gene activities in tumor samples from both GBM and breast cancer. Compared to a competing pathway activity inference approach called SPIA, our method identifies altered activities in cancer-related pathways with fewer false-positives. For computational efficiency, PARADIGM currently uses the NCI pathways as is.
- the NCI's Discovery Therapeutic Program has pursued this approach on large scale, identifying associations between molecular features and responses to >100,000 compounds in a collection of -60 cancer cell lines (Weinstein, J. N. Spotlight on molecular profiling: "Integromic" analysis of the NCI-60 cancer cell lines. Mol Cancer Ther 5, 2601 -2605 (2006); Bussey, K. J. et al. Integrating data on DNA copy number with gene expression levels and drug sensitivities in the NCI-60 cell line panel. Mol Cancer Ther 5, 853-867 (2006)). Although useful for detecting compounds with diverse responses, the NCI60 panel is arguably of limited power in detecting subtype specific responses because of the relatively sparse representation of specific cancer subtypes in the collection.
- the collection carries only 6 breast cancer cell lines, which is not enough to adequately represent the known heterogeneity.
- the MicMa set has been used in parallel pilot studies of whole genome mRNA expression ( Naume, B. et al., (2007), Presence of bone marrow micrometastasis is associated with different recurrence risk within molecular subtypes of breast cancer, 1 : 160- 17), arrayCGH ( Russnes, H. G. et al., (2010), Genomic architecture characterizes tumor progression paths and fate in breast cancer patients, 2: 38ra472), DNA methylation (Ronneberg, J. A.
- KEGG pathways Over 40 KEGG pathways were identified showing differential enrichment according to TP53 mutation status at the p- value cut-off level of 10e-6 in this cohort of breast cancer patients. The differential enrichment of pathways was also observed on the cross-platform dataset consisting of 187 breast cancer samples, based on two different microarray platforms. Differentially enriched pathways included several known cancer pathways such as TP53 signaling and cell cycle, signaling pathways including immune response and cytokine activation and metabolic pathways including fatty acid metabolism (Joshi et al, 2011 supra).
- TP53/mRNA TP53/mRNA
- PARADIGM Pathway Recognition Algorithm using Data integration on Genomic Models
- the clusters and deregulated pathways identified by PARADIGM were then validated in another dataset (Chin, S. F. et al., (2007), Using array-comparative genomic hybridization to define molecular portraits of primary breast cancers, 26: 1959- 1970), and also studied in a dataset of premalignant neoplasia such as DCIS, (ductal carcinoma in situ) (Muggerud, A. A. et al., (2010), Molecular diversity in ductal carcinoma in situ (DCIS) and early invasive breast cancer, 4: 357- 368).
- DCIS ductal carcinoma in situ
- NCI-PID National Cancer Institutes' Pathway Interaction Database
- Table 3 shows the pathways altered by at least three standard deviations with respect to permuted samples found by PARADIGM.
- the FOXMl transcription factor network was altered in the largest number of samples among all pathways tested - 67% of entities with altered activities when averaged across samples.
- pathways with the next highest level of altered activities in the ovarian cohort included PLK1 signaling events (27%), Aurora B signaling (24%), and Thromboxane A2 receptor signaling (20%).
- the FOXMl network harbors significantly more altered activities than other pathways with respect to the ovarian samples.
- FOXMl is a multifunctional transcription factor with three known dominant splice forms, each regulating distinct subsets of genes with a variety of roles in cell proliferation and DNA repair.
- the FOXM lc isoform directly regulates several targets with known roles in cell proliferation including AUKB, PLK1 , CDC25, and BIRC5 (PMID: 15671063).
- the FOXMlb isoform regulates a completely different subset of genes that include the DNA repair genes BRCA2 and XRCC 1 (PMTD: 17101782).
- CHEK2 which is under indirect control of ATM, directly regulates FOXMls expression level.
- the FOXM 1 network was found to have 13 of its genes (52%) differentially regulated, which is a significant proportion based on the hypergeometric test (P ⁇ 3.8* 10 '12 ). Thus, high expression of the FOXM1 network genes does appear to be specifically associated with high-grade disease when compared to the expression of typical genes in the genome.
- FOXM 1 is a multifunctional transcription factor with three known splice forms, each regulating distinct subsets of genes with a variety of roles in cell proliferation and DNA repair.
- An excerpt of FOXM l 's interaction network relevant to this analysis is shown in Figure 27.
- the FOXM la isoform directly regulates several targets with known roles in cell proliferation including AUKB, PLKl , CDC25, and BIRC5.
- the FOXMlb isoform regulates a completely different subset of genes that include the DNA repair genes BRCA2 and XRCC1.
- CHEK2 which is under indirect control of ATM, directly regulates FOXM l 's expression level.
- FOXM1 in most of the ovarian patients, a small subset also have increased copy number amplifications detected by CBS (19% with copy number increases in the top 5% quantile of all genes in the genome measured).
- CBS copy number amplifications detected by CBS (19% with copy number increases in the top 5% quantile of all genes in the genome measured).
- FOXMl may be involved in the control switch between DNA repair and cell proliferation.
- future high- throughput sequencing of the mRNA of these samples may help determine the differential levels of the FOXM 1 isoforms.
- PARADIGM detected the highest level of altered activity centered on this transcription factor suggests that FOXM1 resides at a critical regulatory point in the cell.
- the methods herein described may be used to detect and quantify altered gene expression, absence/presence versus excess, expression of mRNAs or to monitor mRNA levels during therapeutic intervention.
- Conditions, diseases or disorders associated with altered expression include idiopathic pulmonary arterial hypertension, secondary pulmonary hypertension, a cell proliferative disorder, particularly anaplastic oligodendroglioma, astrocytoma, oligoastrocytoma, glioblastoma, meningioma, ganglioneuroma, neuronal neoplasm, multiple sclerosis, Huntington's disease, breast adenocarcinoma, prostate adenocarcinoma, stomach adenocarcinoma, metastasizing neuroendocrine carcinoma, nonproliferative fibrocystic and proliferative fibrocystic breast disease, gallbladder cholecystitis and cholelithiasis, osteoarthritis, and rheumatoid
- gynecomastia actinic keratosis, arteriosclerosis, bursitis, cirrhosis, hepatitis, mixed connective tissue disease (MCTD), myelofibrosis, paroxysmal nocturnal hemoglobinuria, polycythemia vera, primary thrombocythemia, complications of cancer, cancers including adenocarcinoma, leukemia, lymphoma, melanoma, myeloma, sarcoma, teratocarcinoma, and, in particular, cancers of the adrenal gland, bladder, bone, bone marrow, brain, breast, cervix, gall bladder, ganglia, gastrointestinal tract, heart, kidney, liver, lung, muscle, ovary, pancreas, parathyroid, penis, prostate, salivary glands, skin, spleen, testis, thymus, thyroid, and uterus.
- the nucleic acid of the invention in another aspect, the
- the methods described herein may be used to detect and quantify altered gene expression; absence, presence, or excess expression of mRNAs; or to monitor mRNA levels during therapeutic intervention.
- Disorders associated with altered expression include akathesia, Alzheimer's disease, amnesia, amyotrophic lateral sclerosis, ataxias, bipolar disorder, catatonia, cerebral palsy, cerebrovascular disease Creutzfeldt-Jakob disease, dementia, depression, Down's syndrome, tardive dyskinesia, dystonias, epilepsy, Huntington's disease, multiple sclerosis, muscular dystrophy, neuralgias, neurofibromatosis, neuropathies, Parkinson's disease, Pick's disease, retinitis pigmentosa, schizophrenia, seasonal affective disorder, senile dementia, stroke, Tourette's syndrome and cancers including adenocarcinomas, melanomas, and teratocarcinomas, particularly of the brain.
- a normal or standard expression profile is established. This may be accomplished by combining a biological sample taken from normal subjects, either animal or human, with a probe under conditions for hybridization or amplification. Standard hybridization may be quantified by comparing the values obtained using normal subjects with values from an experiment in which a known amount of a substantially purified target sequence is used. Standard values obtained in this manner may be compared with values obtained from samples from patients who are symptomatic for a particular condition, disease, or disorder. Deviation from standard values toward those associated with a particular condition is used to diagnose that condition.
- Such assays may also be used to evaluate the efficacy of a particular therapeutic treatment regimen in animal studies and in clinical trial or to monitor the treatment of an individual patient. Once the presence of a condition is established and a treatment protocol is initiated, diagnostic assays may be repeated on a regular basis to determine if the level of expression in the patient begins to approximate the level that is observed in a normal subject. The results obtained from successive assays may be used to show the efficacy of treatment over a period ranging from several days to months.
- Animal models may be used as bioassays where they exhibit a toxic response similar to that of humans and where exposure conditions are relevant to human exposures. Mammals are the most common models, and most toxicity studies are performed on rodents such as rats or mice because of low cost, availability, and abundant reference toxicology. Inbred rodent strains provide a convenient model for investigation of the physiological consequences of under- or over-expression of genes of interest and for the development of methods for diagnosis and treatment of diseases. A mammal inbred to over-express a particular gene (for example, secreted in milk) may also serve as a convenient source of the protein expressed by that gene.
- Toxicology is the study of the effects of agents on living systems. The majority of toxicity studies are performed on rats or mice to help predict the effects of these agents on human health. Observation of qualitative and quantitative changes in physiology, behavior, homeostatic processes, and lethality are used to generate a toxicity profile and to assess the consequences on human health following exposure to the agent.
- Genetic toxicology identifies and analyzes the ability of an agent to produce genetic mutations.
- Genotoxic agents usually have common chemical or physical properties that facilitate interaction with nucleic acids and are most harmful when chromosomal aberrations are passed along to progeny.
- Toxicological studies may identify agents that increase the frequency of structural or functional abnormalities in progeny if administered to either parent before conception, to the mother during pregnancy, or to the developing organism. Mice and rats are most frequently used in these tests because of their short reproductive cycle that produces the number of organisms needed to satisfy statistical requirements.
- Acute toxicity tests are based on a single administration of the agent to the subject to determine the symptomology or lethality of the agent. Three experiments are conducted: (a) an initial dose-range-finding experiment, (b) an experiment to narrow the range of effective doses, and (c) a final experiment for establishing the dose-response curve.
- Prolonged toxicity tests are based on the repeated administration of the agent. Rats and dog are commonly used in these studies to provide data from species in different families. With the exception of carcinogenesis, there is considerable evidence that daily administration of an agent at high-dose concentrations for periods of three to four months will reveal most forms of toxicity in adult animals.
- Chronic toxicity tests with a duration of a year or more, are used to demonstrate either the absence of toxicity or the carcinogenic potential of an agent.
- studies are conducted on rats, a minimum of three test groups plus one control group are used, and animals are examined and monitored at the outset and at intervals throughout the experiment.
- Transgenic rodents which over-express or under-express a gene of interest may be inbred and used to model human diseases or to test therapeutic or toxic agents.
- the introduced gene may be activated at a specific time in a specific tissue type during fetal development or postnatally. Expression of the transgene is monitored by analysis of phenotype or tissue-specific mRNA expression in transgenic animals before, during, and after challenge with experimental drug therapies.
- Embryonic stem cells isolated from rodent embryos retain the potential to form an embryo. When ES cells are placed inside a carrier embryo, they resume normal development and contribute to all tissues of the live-born animal. ES cells are the preferred cells used in the creation of experimental knockout and knockin rodent strains.
- Mouse ES cells such as the mouse 129/SvJ cell line, are derived from the early mouse embryo and are grown under culture conditions well known in the art.
- Vectors for knockout strains contain a disease gene candidate modified to include a marker gene that disrupts transcription and/or translation in vivo. The vector is introduced into ES cells by transformation methods such as electroporation, liposome delivery, microinjection, and the like which are well known in the art.
- the endogenous rodent gene is replaced by the disrupted disease gene through homologous recombination and integration during cell division.
- Transformed ES cells are identified, and preferably microinjected into mouse cell blastocysts such as those from the C57BL 6 mouse strain.
- the blastocysts are surgically transferred to pseudopregnant dams and the resulting chimeric progeny are genotyped and bred to produce heterozygous or homozygous strains.
- ES cells are also used to study the differentiation of various cell types and tissues in vitro, such as neural cells, hematopoietic lineages, and cardiomyocytes (Bain et al. (1995) Dev. Biol. 168: 342-357; Wiles and Keller (1991 ) Development 111 : 259-267; and Klug et al. (1996) J. Clin. Invest. 98: 216-224).
- ES cells derived from human blastocysts may also be manipulated in vitro to differentiate into eight separate cell lineages, including endoderm, mesoderm, and ectodermnal cell types (Thomson (1998) Science 282: 1145- 1147).
- a region of a human disease gene candidate is enzymatically modified to include a non-mammalian gene such as the neomycin phosphotransferase gene (neo; see, for example, Capecchi (1989) Science 244: 1288- 1292).
- the inserted coding sequence disrupts transcription and translation of the targeted gene and prevents biochemical synthesis of the disease candidate protein.
- the modified gene is transformed into cultured embryonic stem cells (described above), the transformed cells are injected into rodent blastulae, and the blastulae are implanted into pseudopregnant dams. Transgenic progeny are crossbred to obtain homozygous inbred lines.
- Totipotent ES cells present in the early stages of embryonic development, can be used to create knockin humanized animals (pigs) or transgenic animal models (mice or rats) of human diseases.
- knockin technology a region of a human gene is injected into animal ES cells, and the human sequence integrates into the animal cell genome by recombination.
- Totipotent ES cells that contain the integrated human gene are handled as described above. Inbred animals are studied and treated to obtain information on the analogous human condition. These methods have been used to model several human diseases. (See, for example, Lee et al. (1998) Proc. Natl. Acad. Sci. 95: 11371-11376; Baudoin et al. (1998) Genes Dev. 12: 1202- 1216; and Zhuang et al. (1998) Mol. Cell Biol. 18: 3340-3349).
- NHPs are the first choice test animal.
- NHPs and individual humans exhibit differential sensitivities to many drugs and toxins and can be classified as "extensive metabolizers" and “poor metabolizers” of these agents.
- CTLs cytotoxic T lymphocytes
- NK natural killer
- This disclosure also provides the first large scale integrative view of the aberrations in HGS-OvCa.
- the mutational spectrum was surprisingly simple. Mutations in TP53 predominated, occurring in at least 96% of HGS-OvCa while BRCA l /2 were mutated in 22% of tumors due to a combination of germline and somatic mutations. Seven other significantly mutated genes were identified, but only in 2-6% of HGS-OvCa. In contrast, HGS-OvCa demonstrates a remarkable degree of genomic disarray.
- clear-cell OvCa have few TP53 mutations but have recurrent ARID] A and PIK3CA47-49 mutations; endometrioid OvCa have frequent CTTNB1, A RID I A, and P1K3CA mutations and a lower rate of TP5348.49 while mucinous OvCa have prevalent KRAS mutations50.
- These differences between ovarian cancer subtypes likely reflect a combination of etiologic and lineage effects, and represent an opportunity to improve ovarian cancer outcomes through subtype-stratified care.
- Ecosystem 100 can include pathway element database 120 preferably storing a plurality of pathway elements 125A through 125N, collectively referred to as pathway elements 125.
- Each of pathway elements 125 can be characterized by its involvement with one or more pathways.
- Elements 125 can be considered separately manageable data objects comprising one or more properties or values describing the characteristics of the element.
- elements 125 can be considered an n-tuple of properties or values, where each property member of an element 125 tuple can be compared, analyzed, contrasted, or otherwise evaluated against other property members in other element tuples.
- Modification engine 110 communicatively couples with pathway element database 120, possibly over a network link (for example, LAN, WAN, Internet, VPN, etc.).
- pathway element database 120 could be local to modification engine 110, while in other embodiments, pathway element database 120 could be remote from modification engine 110.
- pathway element database 120 could be accessed via the National Lambda Rail (see URL www.nlr.net) or the Internet.
- modification engine 110, or ecosystem 100 for that matter can be accessed by users over the network, possibly in exchange for a fee.
- Modification engine 110 obtains one or more of elements 125 from pathway element database 120 for analysis.
- modification engine 110 associates at least one of elements 125 (for example, elements 125A) with at least one a priori known attribute 133.
- modification 110 also associates another element, element 125N for example, with assumed attribute 137.
- modification engine 110 can make the associations automatically based on inference rules, programmatic instructions, or other techniques. For example, known attributes 137 could be obtained from known research while assumed attributes 137 could be mapped out according to an attribute parameterized space where modification engine 110 serially, or in parallel, walks through the assumed attribute space.
- a user can manually associate attributes 133 or 137 as desired through one or more user interfaces (not shown), possibly operating through an HTTP server or other suitable interfacing technology.
- Modification engine 110 further cross-correlates pathway elements 125 for one or more pathways using known attributes 133 and assumed attributes 137. Further, modification engine 110 assigns one or more influence levels 145 to elements 125. Through cross-correlation and assignment of influence levels 145, modification engine 110 constructs probabilistic pathway model 140 outlining how pathways might be influenced by assumed attributes 137 or other factors.
- probabilistic pathway model 140 can be stored within pathway model database 150 for archival purposes, or for analysis as indicated. As with elements 125, probabilistic pathway model 140 can also be stored as a distinct manageable data object having properties or values describing the characteristics of the model, possibly as an n-tuple. Models 145, or even elements 125, can be stored according to any desirable schema.
- Example suitable database that can be used to construct element database 120 or model database 140 include MySQL, PostgreSQL, Oracle, or other suitable databases.
- the data objects for example, elements 125, probabilistic pathway model 145, etc.
- Ecosystem 100 preferably includes analysis engine 160 configured to further analyze probabilistic pathway model 150 with respect to actual data.
- analysis engine 160 obtains probabilistic pathway model 150, possibly under direction of a user or researcher, to derive dynamic pathway model 165.
- dynamic pathway model 160 is derived by comparing one or more measured attributes 173 from a patient sample with the attributes associated with probabilistic pathway model 140.
- analysis engine 160 seeks to modify, update, correct, or otherwise validate probabilistic pathway model 140 to form dynamic pathway model 165.
- dynamic pathway model 165 can be stored within a model database.
- analysis engine 160 can configure one or more output devices (for example, a display, a printer, a web server, etc.) to present dynamic pathway model 165.
- a system will therefore typically include a pathway element database.
- the database may be physically located on a single computer, however, distributed databases are also deemed suitable for use herein.
- the particular format of the database is not limiting to the inventive subject matter so long as such database is capable of storing and retrieval of multiple pathway elements, and so long as each pathway element can be characterized by its involvement in at least one pathway.
- suitable pathway elements will include one or more proteins (which may or may not be modified, for example, via glycosylation, myristoylation, etc.), alone or in complex with other cellular components, various nucleic acids (genomic DNA, extrachromosomal DNA, hnRNA, siRNA, mRNA, rRNA, etc) which may be a native nucleic acid or a recombinant nucleic acid, lipids, hormones, second messengers, and pharmaceutically active agents provided as therapeutic or preventive agent.
- contemplated pathway elements may have a variety of functions, and especially preferred functions include various enzymatic functions.
- suitable functions are examples of suitable functions.
- kinases/phosphatases polymerases/hydrolases, proteases, hydrolases (and especially GTPase), hydroxylases, methyl transferases/methylases, etc.
- pathway element is a protein
- suitable pathway elements include various receptors, hormone binding proteins, kinases, transcription factors, initiation factors, methylases and methyl transferases, histone acetylases, and histone deacetylases.
- the pathway element is a nucleic acid
- contemplated pathway elements will include those that encode a protein sequence, one or more genomic regulatory sequences, regulatory RNA, and a trans-activating sequences.
- pathway element Depending on the particular pathway element, it should therefore be appreciated that the nature of the pathway may vary considerably, and all known pathways are deemed suitable for use herein.
- contemplated pathways may be involved in signal transduction, in cell cycling, in cell growth and/or metabolism, in repair mechanisms (and especially in DNA repair), and in neural signaling.
- especially preferred pathways include calcium/calmodulin dependent signaling pathways and functionally associated pathway networks, a cytokine mediated signaling pathway and functionally associated pathway networks, a chemokine mediated signaling pathway and functionally associated pathway networks, a growth factor signaling pathway and functionally associated pathway networks, a hormone signaling pathway and functionally associated pathway networks, a MAP kinase signaling pathway and functionally associated pathway networks, a phosphatase mediated signaling pathway and functionally associated pathway networks, a Ras superfamily mediated signaling pathway and functionally associated pathway networks, and a transcription factor mediated signaling pathway and functionally associated pathway networks. Therefore, it should be appreciated that the pathways may be individual pathways as well as pathways within a pathway network, and even within a network of distinct pathway networks.
- contemplated herein may be within a regulatory pathway network.
- contemplated pathway networks include an ageing pathway network, an apoptosis pathway network, a homeostasis pathway network, a metabolic pathway network, a replication pathway network, and an immune response pathway network.
- the type and numerical value of the attribute of the pathway element may vary considerably, and that the particular pathway element will in large part determine the type and numerical value of the attribute.
- the attribute may be a copy number, a particular haplotype or mutation, strength of a regulatory element (for example, promoter, repressor, etc.), transcription level or translation level.
- contemplated attributes will also include class attributes (for example, gene is activatable by particular transcription factor, or sensitive to particular hormone response element, etc.) or may be a compound attribute (for example, representative of at least two different attributes).
- the pathway element is a protein
- the attribute may be quantity of translation, protein activity, requirement for a cofactor, requirement for formation of a multi-protein complex to provide activity, etc.).
- the attributes for at least some of the pathway elements will be known from prior study and publication and can therefore be used in contemplated systems and methods as a priori known attributes for the specific pathway element.
- numerous attributes are not known a priori, however, that a large variety of such unknown attributes can be assumed with a reasonably good expectation of accuracy.
- the pathway element is a genomic sequence for a receptor, and where that sequence is preceded by a trans-activator binding sequence element, it can be reasonably assumed that one attribute of the pathway element is a requirement for binding of a trans activator.
- the strength of trans-activation is known for similarly controlled sequences, the transcription level of the pathway element can be reasonably inferred.
- the assumed attributes are not arbitrarily assumed values, but that the assumption is based on at least partially known information.
- the kind and value of the assumed attribute is also a function of a reference pathway.
- the reference pathway is a pathway of a healthy cell.
- the numerical range and kind of attribute will typically be reflective of that of a normal cell.
- non-normal cells may also be used to establish a reference pathways.
- Such pathway patterns are typically produced using a modification engine that is coupled to the pathway element database, wherein the modification engine is used (1 ) to associate a first pathway element with at least one a priori known attribute, (2) to associate a second pathway element with at least one assumed attribute, and (3) to cross-correlate and assign an influence level of the first and second pathway elements for at least one pathway using the known and assumed attributes, respectively, to ultimately form a probabilistic pathway model.
- association of the first pathway element with at least one a priori known attribute can be done in numerous manners. However, it is particularly preferred that the attribute is expressed as one of an n-tuple of attributes that is directly associated with the pathway element.
- the known attribute is derived from a peer-reviewed publication.
- secondary information sources for example, compiled and publicly available information from various databases such as SWISSPROT, EMBL, OMIM, NCI-PID, Reactome, Biocarta, KEGG, etc.
- assumed attributes can be manually associated with the pathway element, and more preferably in an at least semi-automated manner.
- Cross-correlation can be achieved through numerous techniques.
- pathway elements can be cross-correlated manually.
- elements can be cross-correlated through one or more automated techniques.
- numerous elements can be analyzed with respect to their properties via a modification engine that seeks to find possible correlation.
- the modification engine can be configured to seek such correlations via multi-variate analysis, genetic algorithms, inference reasoning, or other techniques.
- inference reason could include application of various forms of logic including deductive logic, abductive logic, inductive logic, or other forms of logic.
- contemplated engines are capable of discovering possible correlations that a researcher might otherwise overlook.
- Another example of inference reasoning can include applications using inference on probabilistic models such as belief propagation, loopy belief propagation, junction trees, variable elimination or other inference methods.
- Influence levels represent a quantitative value that an assumed attribute has on a pathway comprising elements with known attributes.
- Influence levels can comprise single values or multiple values.
- Example of a single value could include a weighting factor, possibly as an absolute value or a normalized value relative to other known influences within the pathway system under evaluation.
- Example multi-valued influence levels can include a range of values with a possible distribution width.
- initial values of an influence level can be established through various techniques including being manually set. In more preferred embodiments, the initial value can be established through a manual estimation formulated by the modification engine.
- the relative "distance" according to one or more element or pathway properties can be used to weight an influence level. The distance can be the exact distance or can be the square of the distance.
- the influence levels can be determined by maximizing the likelihood of the influence levels between all of the other values within the pathway system.
- an analysis engine will then employ the probabilistic pathway model to derive a dynamic pathway map from a plurality of measured attributes for a plurality of elements of a patient sample.
- a patient sample may be derived from a biological fluid, a biopsy, or surgical specimen, and will typically analyzed using methods well known in the art. Therefore, and among other suitable attributes, measured attributes will include mutations, differential genetic sequence object, gene copy number for one or more particular gene, transcription level for one or more particular gene, translation level for one or more particular protein, protein activity, protein interaction, presence and/or quantity of an analyte (for example, metabolite) or marker of a disease, etc.
- an analyte for example, metabolite
- the measured attributes are fed into the probabilistic pathway model to so arrive at a dynamic pathway map that can indicate deviations from the probabilistic pathway model.
- the dynamic pathway map will provide a user with reference pathway activity information for a particular pathway (which can be specific with respect to a normal tissue, a diseased tissue, an ageing tissue, or a recovering tissue, etc.). Consequently, and viewed form a different perspective, the dynamic pathway map will allow a user to readily identify information related to one or more pathways in a patient sample based on a relatively limited number of measured attributes.
- the inventors also contemplate a method of generating a dynamic pathway map in which a user is provided access to a model database that stores a probabilistic pathway model that comprises a plurality of pathway elements.
- access may be controlled in a variety of manners as the particular access protocol will be at least in part determined on the particular use.
- the access is a pay-per-use access or a pre-authorized access.
- the model database may also be accessible via a publicly accessible network.
- the plurality of pathway elements are cross-correlated and assigned an influence level for at least one pathway on the basis of known attributes, and that another number of the plurality of pathway elements are cross-correlated and assigned an influence level for at least one pathway on the basis of assumed attributes, and that an analysis engine modifies the probabilistic pathway model with a plurality of measured attributes for a plurality of elements of a patient sample to obtain the dynamic pathway map, wherein the dynamic pathway map includes most preferably reference pathway activity information for a particular pathway.
- contemplated systems and methods are not only suitable for analysis of a first sample relative to a standard pathway model (for example, representing healthy donor), but that such systems and methods also allow intra-patient analysis of diseased tissue vis-a-vis healthy tissue to so predict a pathway activity information for a tissue. Therefore, using two samples from the same patient (i.e., from diseased tissue and non-diseased tissue), susceptibility of a diseased tissue to certain pharmaceuticals can be predicted.
- the inventors also contemplate a method of analyzing biologically relevant information in which access to a model database is provided that stores a dynamic pathway map, wherein the DPM is generated by modification of a probabilistic pathway model with a plurality of measured attributes for a plurality of elements of a first cell or patient sample. Subsequently, a plurality of measured attributes are obtained for a plurality of elements of a second cell or patient sample, and the dynamic pathway map and the plurality of measured attributes for the plurality of elements of the second cell or patient sample are then used by an analysis engine to determine a predicted pathway activity information for the second cell or patient sample.
- the measured attributes for the plurality of elements of the first cell or patient sample may be characteristic for a healthy cell or tissue, a specific age of a cell or tissue, a specific disease of a cell or tissue, a specific disease stage of a diseased cell or tissue, a specific gender, a specific ethnic group, a specific occupational group, and even a specific species. So computed information may provide valuable information about actual or likely pathway differences with respect to occupation, pharmaceutical treatment, predisposition to a disease, etc. Thus, first and second samples may obtained from the same cell or patient concurrently, or at different times (most typically after treatment has commenced).
- a patient or patient sample may be subjected to a treatment (typically surgery, radiation, and administration of a pharmaceutical), and then receive a second pharmaceutical of potential therapeutic value.
- a treatment typically surgery, radiation, and administration of a pharmaceutical
- the predicted pathway activity information may be able to identify a pathway element as a hierarchical-dominant element in at least one pathway, and/or as a disease-determinant element in at least one pathway with respect to a disease. Consequently, pharmaceutical intervention can be used in a targeted fashion with high chances of achieving the desired outcome.
- the predicted pathway activity information is provided to a physician, it is generally preferred to generate a graphical representation of predicted pathway activity information to render the information more relevant to the needs of a practitioner.
- the predicted pathway activity information may be used by the system and/or user to formulate a diagnosis, a prognosis for a disease, or a recommendation (for example, selection of a treatment option or dietary guidance).
- the predicted pathway activity information can also be used to identify an epigenetic factor, a stress adaptation, a state of an organism, and/or a state of repair or healing.
- polynucleotide nucleic acids may be used in any molecular biology techniques that have yet to be developed, provided the new techniques rely on properties of nucleic acid molecules that are currently known, including, but not limited to, such properties as the triplet genetic code and specific base pair interactions.
- Probe annotations were converted to BED15 format for display in the UCSC Cancer Genomics Browser (Zhu:2009, supra) and subsequent analysis. Array data were mapped to probe annotations via probe ID. Matched expression data from Naderi (2007, supra) was obtained from MIAMIExpress at EBI using accession number E-UCon- 1. Platform annotation information for Human lA (V2) was obtained from the Agilent website.Expression data was probe-level median-normalized and mapped via probe ID to HUGO gene names.
- the glioblastoma data from TCGA was obtained from the TCGA Data Portal providing gene expression for 230 patient samples and 10 adjacent normal tissues on the Affymetrix U133A platform.
- the probes for the patient samples were normalized to the normal tissue by subtracting the median normal value of each probe.
- CBS segmented (Olshen:2004 supra pl618) copy number data for the same set of patients were obtained. Both datasets were non-parametrically normalized using the same procedure as the breast cancer data.
- each NCI pathway into a distinct probabilistic model.
- a toy example of a small fragment of the p53 apoptosis pathway is shown in Figure 2.
- a pathway diagram from NCI was converted into a factor graph that includes both hidden and observed states.
- the factor graph integrates observations on gene- and biological process-related state information with a structure describing known interactions among the entities.
- the factor graph encodes the state of a cell using a random variable for each entity X - ⁇ und*,,...., *, preparation ⁇ and a set of m non-negative functions, or factors, that constrain the entities to take on biologically meaningful values as functions of one another.
- the -th factor ⁇ defines a probability distribution over a subset of entities X, c .
- Each entity can take on one of three states corresponding to activated, nominal, or deactivated relative to a control level (for example, as measured in normal tissue) and encoded as 1 , 0, or -1 respectively.
- the states may be interpreted differently depending on the type of entity (for example, gene, protein, etc).
- an activated mRNA entity represents overexpression, while an activated genomic copy entity represents more than two copies are present in the genome.
- Figure 2 shows the conceptual model of the factor graph for a single protein-coding gene.
- entities are introduced to represent the copy number of the genome (C DNA ), mRNA expression (G mRNA ), protein level (G protein ), and protein activity (G protein ) (ovals labeled "DNA”, “mRNA”, “protein”, and “active” in Figure 2).
- C DNA copy number of the genome
- G mRNA mRNA expression
- G protein protein level
- protein activity G protein
- a single variable with molecular type For every compound, protein complex, gene family, and abstract process in the pathway, we include a single variable with molecular type " active.”
- the expected value was set to the majority vote of the parent variables. If a parent is connected by a positive edge it contributes a vote of +1 times its own state to the value of the factor. Conversely, if the parent is connected by a negative edge, then the variable votes -1 times its own state.
- by an edge labeled "minimum” get a single vote, and that vote's value is the minimum value of these variables, creating an AND-like connection.
- by an edge labeled "maximum” get a single vote, and that vote's value is the maximum value of these variables, creating an OR-like connection. Votes of zero are treated as abstained votes. If there are no votes the expected state is zero.
- Array CGH/SNP estimates of copy number alteration are associated with the 'genome' type.
- Gene expression data is associated with the 'mRNA' type.
- future expansion will include DNA methylation data with the 'mRNA' type, and proteomics and gene-resequencing data with the 'protein' and 'active' types.
- Each observation variable is also ternary valued. The factors associated with each observed type of data are shared across all entities and learned from the data, as described next.
- EM Expectation-Maximization
- the IPA score reflects a signed analog of the log-likelihood ratio, L.
- the IPA is set to L.
- the ⁇ is set to the negative of the log likelihood ratio. If the gene is most likely unchanged, the IPA is set to zero.
- Each pathway is analyzed independently of other pathways. Therefore, a gene can be associated with multiple inferences, one for each pathway in which it appears. Differing inferences for the same gene can be viewed as alternative interpretations of the data as a function of the gene's pathway context.
- IPA scores we assess the significance of IPA scores by two different permutations of the data.
- a permuted data sample is created by choosing a new tuple of data (i.e. matched gene expression and gene copy number) first by choosing a random real sample, and then choosing a random gene from within the same pathway, until tuples have been chosen for each gene in the pathway.
- the procedure is the same, but the random gene selection step could choose a gene from anywhere in the genome.
- 1 ,000 permuted samples are created, and the perturbation scores for each permuted sample is calculated. The distribution of perturbation scores from permuted samples is used as a null distribution to estimate the significance of true samples.
- SPIA Signaling Pathway Impact Analysis
- a set of decoy pathways was created for each cancer dataset.
- Each NCI pathway was used to create a decoy pathway which consisted of the same structure but where every gene in the pathway was substituted for a random gene in RefGene. All complexes and abstract processes were kept the same and the significance analysis for both PARADIGM and SPIA was run on the set of pathways containing both real and decoy pathways. The pathways were ranked within each method and the fraction of real versus total pathways was computed and visualized.
- the anti-apoptotic AKTl serine-threonine kinase is known to be involved in breast cancer and interacts with the ERBB2 pathway (Ju X, Katiyar S, Wang C, Liu M, Jiao X, Li S, Zhou J, Turner J, Lisanti MP, Russell RG, Mueller SC, Ojeifo J, Chen WS, Hay N, Pestell RG. Aktl governs breast cancer progression in vivo. Proc. Natl. Acad. Sci. U.S.A. 2007 May; 104(18):7438-7443).
- CircleMap visualization to display multiple datasets centered around each gene in a pathway ( Figure 7).
- each gene is associated with all of its data across the cohort by plotting concentric rings around the gene, where each ring corresponds to a single type of measurement or computational inference.
- Each tick in the ring corresponds to a single patient sample while the color corresponds to activated (red), deactivated (blue), or unchanged (white) levels of activity.
- CircleMaps for a subset of the ErbB2 pathway and included ER status, IP As, expression, and copy number data from the breast cancer cohort.
- the fourth cluster exhibits clear downregulation of HIF- 1 -alpha transcription factor network as well as overexpression of the E2F transcription factor network.
- HIF- 1 -alpha is a master transcription factor involved in regulation of the response to hypoxic conditions.
- two of the first three clusters have elevated EGFR signatures and an inactive MAP kinase cascade involving the GATA interleukin transcriptional cascade.
- mutations and amplifications in EGFR have been associated with high grade gliomas as well as glioblastomas (Kuan CT, Wikstrand CJ, Bigner DD.
- EGF mutant receptor vIII as a molecular target in cancer therapy. Endocr. Relat. Cancer 2001 Jun.;8(2):83-96). Amplifications and certain mutations can create a constitutively active EGFR either through self stimulation of the dimer or through ligand-independent activation. The constitutive activation of EGFR may promote oncogenesis and progression of solid tumors. Gefitinib, a molecule known to target EGFR, is currently being investigated for its efficacy in other EGFR-driven cancers. Thus, qualitatively, the clusters appeared to be honing in on biologically meaningful themes that can stratify patients.
- Cluster 4 The survival analysis revealed that the patients in cluster 4 have a significantly better survival profile.
- Cluster 4 was found to have an up-regulation of E2F, which acts with the retinoblastoma tumor suppressor. Up- regulation of E2F is therefore consistent with an active suppression of cell cycle progression in the tumor samples from the patients in cluster 4.
- cluster 4 was associated with an inactivity of the HIF- 1 -alpha transcription factor. The inactivity in the fourth cluster may be a marker that the tumors are more oxygenated, suggesting that they may be smaller or newer tumors.
- PARADIGM IP As provide a meaningful set of profiles for delineating subtypes with markedly different survival outcomes.
- Example XIII Integrated Genomic Analyses of Ovarian Carcinoma: Samples and clinical data.
- This report covers analysis of 489 clinically annotated stage II-IV HGS-OvCa and corresponding normal DNA. Patients reflected the age at diagnosis, stage, tumor grade, and surgical outcome of individuals diagnosed with HGS-OvCa. Clinical data were current as of August 25, 2010. HGS-OvCa specimens were surgically resected before systemic treatment but all patients received a platinum agent and 94% received a taxane. The median progression-free and overall survival of the cohort is similar to previously published trials 11 , 12. Twenty five percent of the patients remained free of disease and 45% were alive at the time of last follow-up, while 31% progressed within 6 months after completing platinum-based therapy. Median follow up was 30 months (range 0 to 179). Samples for TCGA analysis were selected to have > 70% tumor cell nuclei and ⁇ 20% necrosis.
- Exome capture and sequencing was performed on DNA isolated from 316 HGS-OvCa samples and matched normal samples for each individual. Capture reagents targeted -180,000 exons from -18,500 genes totaling -33 megabases of non-redundant sequence. Massively parallel sequencing on the Illumina GAIIx platform (236 sample pairs) or ABI SOLiD 3 platform (80 sample pairs) yielded -14 gigabases per sample (-9x 109 bases total). On average, 76% of coding bases were covered in sufficient depth in both the tumor and matched normal samples to allow confident mutation detection. 19,356 somatic mutations (-61 per tumor) were annotated and classified in Table 4.
- Mutations that may be important in HGS-OvCa pathophysiology were identified by (a) searching for non-synonymous or splice site mutations present at significantly increased frequencies relative to background, (b) comparing mutations in this study to those in COSMIC and OMIM and (c) predicting impact on protein function.
- Example XV Copy number analysis.
- SCNAs Somatic copy number alterations present in the 489 HGS-OvCa genomes were identified and compared with glioblastome multiforme data in Figure 37A. SCNAs were divided into regional aberrations that affected extended chromosome regions and smaller focal aberrations. A statistical analysis of regional aberrations identified 8 recurrent gains and 22 losses, all of which have been reported previously22 ( Figure 37B). Five of the gains and 18 of the losses occurred in more than 50% of tumors.
- GISTIC was used to identify recurrent focal SCNAs. This yielded 63 regions of focal amplification ( Figure 37C) including 26 that encoded 8 or fewer genes.
- the most common focal amplifications encoded CCNE1, MYC, and MECOM ( Figure 37C) each highly amplified in greater than 20% of tumors.
- New tightly- localized amplification peaks in HGS-OvCa encoded the receptor for activated C-kinase, ZMYND8; the p53 target gene, IRF2BP2; the DNA-binding protein inhibitor, ID4; the embryonic development gene, PAX8; and the telomerase catalytic subunit, TERT.
- GISTIC also identified 50 focal deletions.
- the known tumor suppressor genes PTEN, RBI, and NF1 were in regions of homozygous deletions in at least 2% of tumors.
- RBI and NF1 also were among the significantly mutated genes.
- One deletion contained only three genes, including the essential cell cycle control gene, CREBBP, which has 5 non-synonymous and 2 frameshift mutations.
- Example XVI mR A and miRNA expression and DNA methylation analysis.
- T-cell chemokine ligands CXCLI I and CXCL10, and the receptor, CXCR3, characterized the Immunoreactive subtype.
- High expression of transcription factors such as HMGA2 and SOXI 1, low expression of ovarian tumor markers (MUCl, MUCl 6) and high expression of proliferation markers such as MCM2 and PCNA defined the Proliferative subtype.
- MUCl, MUCl 6 ovarian tumor markers
- MCM2 and PCNA defined the Proliferative subtype.
- the Differentiated subtype was associated with high expression of MUCl 6 and MUCl and with expression of the secretory fallopian tube maker SLP1, suggesting a more mature stage of development.
- High expression of HOX genes and markers suggestive of increased stromal components such as for myofibroblasts (FAP) and microvascular pericytes (ANGPTL2, ANGPTL1) characterized the Mesenchymal subtype.
- FAP myofibroblasts
- Elevated DNA methylation and reduced tumor expression implicated 168 genes as epigenetically silenced in HGS-OvCa compared to fallopian tube controls26. DNA methylation was correlated with reduced gene expression across all samples. AMT, CCL21 and SPARCL1 were noteworthy because they showed promoter hypermethylation in the vast majority of the tumors. Curiously, RAB25, previously reported to be amplified and over-expressed in ovarian cancer, also appeared to be epigenetically silenced in a subset of tumors. The BRCA 1 promoter was hypermethylated and silenced in 56 of 489 (11.5%) tumors as previously reported.
- Consensus clustering of variable DNA methylation across tumors identified four subtypes that were significantly associated with differences in age, BRCA inactivation events, and survival. However, the clusters demonstrated only modest stability. [00234] Survival duration did not differ significantly for transcriptional subtypes in the TCGA dataset. The Proliferative group showed a decrease in the rate of MYC amplification and RB I deletion, whereas the
- Immunoreactive subtype showed an increased frequency of 3q26.2 ( ECOM) amplification.
- ECOM ECOM
- a 193 gene transcriptional signature predictive of overall survival was defined using the integrated expression data set from 215 samples. After univariate Cox regression analysis, 108 genes were correlated with poor survival, and 85 were correlated with good survival (p-value cutoff of 0.01). The predictive power was validated on an independent set of 255 TCGA samples as well as three independent expression data sets25,29,30. Each of the validation samples was assigned a prognostic gene score, reflecting the similarity between its expression profile and the prognostic gene signature31 (Figure 38c). Kaplan- Meier survival analysis of this signature showed statistically significant association with survival in all validation data sets ( Figure 38d).
- NMF consensus clustering of miRNA expression data identified three subtypes. Interestingly, miRNA subtype 1 overlapped the mRNA Proliferative subtype and miRNA subtype 2 overlaped the mKNAMesenchymal subtype (Figure 38d). Survival duration differed significantly between iRNA subtypes with patients in miRNA subtype 1 tumors surviving significantly longer ( Figure 38e).
- Example XVII Pathways influencing disease.
- a probabilistic graphical model searched for altered pathways in the NCI Pathway Interaction Database identifying the FOXM1 transcription factor network (Figure 39D) as significantly altered in 87% of cases.
- FOXM1 and its proliferation-related target genes A URB, CCNBl, BIRC5, CDC25, and PLK1, were consistently over-expressed but not altered by DNA copy number changes, indicative of transcriptional regulation.
- TP53 represses FOXM1 following DNA damage42, suggesting that the high rate of TP53 mutation in HGS-OvCa contributes to FOXM1 overexpression.
- the FOXM1 pathway is significantly activated in tumors relative to adjacent epithelial tissue and is associated with HGS-OvCa.
- Example XVIII Frequently altered pathways in ovarian serous carcinomas
- Table 3 shows the pathways altered by at least three standard deviations with respect to permuted samples found by PARADIGM.
- the FOXM1 transcription factor network was altered in the largest number of samples among all pathways tested - 67% of entities with altered activities when averaged across samples.
- pathways with the next highest level of altered activities in the ovarian cohort included PLK1 signaling events (27%), Aurora B signaling (24%), and Thromboxane A2 receptor signaling (20%).
- the FOXM1 network harbors significantly more altered activities than other pathways with respect to the ovarian samples.
- FOXM1 transcription factor network was found to be differentially altered in the tumor samples compared to the normal controls in the highest proportion of the patient samples (Figure 29).
- FOXM1 is a multifunctional transcription factor with three known dominant splice forms, each regulating distinct subsets of genes with a variety of roles in cell proliferation and DNA repair.
- the FOXMlc isoform directly regulates several targets with known roles in cell proliferation including AUKB, PLKl , CDC25, and BIRC5.
- the FOXMlb isoform regulates a completely different subset of genes that include the DNA repair genes BRCA2 and XRCC1.
- CHEK2 which is under indirect control of ATM, directly regulates FOXMls expression level.
- FOXMl 's role in many different cancers including breast and lung has been well documented but its role in ovarian cancer has not been investigated.
- FOXMl is a multifunctional transcription factor with three known splice variants, each regulating distinct subsets of genes with a variety of roles in cell proliferation and DNA repair.
- An excerpt of FOXMl 's interaction network relevant to this analysis is shown as Figure 27.
- the FOXMla isoform directly regulates several targets with known roles in cell proliferation including AUKB, PLKl , CDC25, and BIRC5.
- the FOXM lb isoform regulates a completely different subset of genes that include the DNA repair genes BRCA2 and XRCC 1.
- CHEK2 which is under indirect control of ATM, directly regulates FOXMl 's expression level.
- FOXMl In addition to increased expression of FOXMl in most of the ovarian patients, a small subset also have increased copy number amplifications detected by CBS ( 19% with copy number increases in the top 5% quantile of all genes in the genome measured).
- CBS copy number amplifications detected by CBS ( 19% with copy number increases in the top 5% quantile of all genes in the genome measured).
- the alternative splicing regulation of FOXMl may be involved in the control switch between DNA repair and cell proliferation.
- future high- throughput sequencing of the mRNA of these samples may help determine the differential levels of the FOXMl isoforms.
- PARADIGM detected the highest level of altered activity centered on this transcription factor suggests that FOXMl resides at a critical regulatory point in the cell.
- a collection of pathways was obtained from NCI-PID containing 131 pathways, 11 ,563 interactions, and 7,204 entities.
- An entity is molecule, complex, small molecule, or abstract concept represented as "nodes” in PARADIGM' S graphical model.
- the abstract concepts correspond to general cellular processes (such as "apoptosis” or "absorption of light,") and families of genes that share functional activity such as the RAS family of signal transducers.
- We collected interactions including protein-protein interactions, transcriptional regulatory interactions, protein modifications such as phosphorylation and ubiquitinylation interactions.
- Example XX Inference of integrated molecular activities in pathway context.
- PARADIGM which assigns an integrated pathway activity (IPA) reflecting the copy number, gene expression, and pathway context of each entity.
- IPAs The significance of IPAs was assessed using permutations of gene- and patient-specific cross-sections of data. Data for 1000 "null" patients was created by randomly selecting a gene-expression and copy number pair of values for each gene in the genome. To assess the significance of the PARADIGM IPAs, we constructed a null distribution by assigning random genes to pathways while preserving the pathway structure.
- Example XXII Clustering [00253] The use of inferred activities, which represent a change in probability of activity and not activity directly, it enables entities of various types to be clustered together into one heatmap. To globally visualize the results of PARADIGM inference, Eisen Cluster 3.0 was used to perform feature filtering and clustering. A standard deviation filtering of 0.1 resulted in 1598 out of 7204 pathway entities remaining, and average linkage, uncentered correlation hierarchical cluster was performed on both the entities and samples.
- Example XXIII Cell lines model many important tumor subtypes and features.
- HCC hierarchical consensus clustering
- FIG. 14B A refined high-resolution SNP copy number analysis (Figure 14B) confirms that the cell line panel models regions of recurrent amplification at 8q24 (MYC), l lql 3 (CCND1), 17ql2 (ERBB2), 20q l3 ⁇ STK15/AURKA), and homozygous deletion at 9p21 (CDKN2A) found in primary tumors.
- MYC l lql 3
- ERBB2 l lql 3
- ERBB2 17ql2
- CDKN2A homozygous deletion at 9p21
- FIG. 10A A representative waterfall plot illustrating the variation in response to the Sigma AKT1 -2 inhibitor along with associated transcriptional subtypes is shown in Figure 10A. Sensitivity to this compound is highest in luminal and ERBB2 AMP and lower in basal and claudin-low breast cancer cell lines. Waterfall plots showing the distribution of GI 50 values among the cell lines for all compounds are plotted. We established the reproducibility of the overall data set by computing the median absolute deviation of GI 50 values for 229 compound/cell line combinations with 3 or 4 replicates. The median average deviation was 0.15 across these replicates ( Figure 15). We assessed concordance of response to 8 compounds by computing the pairwise Pearson's correlation between sets of GI50 values ( Figure 15B. Sensitivities for pairs of drugs with similar mechanisms of action were highly correlated, suggesting similar modes of action.
- Example XV Many compounds were preferentially effective in subsets of the cell lines.
- a central premise of this study is that associations between responses and molecular subtypes observed in preclinical cell line analyses will be recapitulated in the clinic in instances where the predictive molecular features in the cell lines are mirrored in human tumors.
- Figure IOC shows a hierarchical clustering of the 34 agents with significant associations with one or more of the luminal, basal, claudin-low and ERBB2 AMP subtypes.
- the 11 agents most strongly associated with subtype were inhibitors of receptor tyrosine kinase signaling and histone deacetylase and had the highest efficacy in luminal and/or ERBB2 AMP cell lines.
- the three next most subtype-specific agents - etoposide, cisplatin, and docetaxel - show preferential activity in basal and/or claudin-low cell lines as observed clinically.
- Agents targeting the mitotic apparatus including ixabepilone, GSK461364 (polo kinase inhibitor) and
- GSK1070916 also were more active against basal and claudin-low cell lines. AG1478, BIBW2992 and gefitinib, all of which target EGFR and/or ERBB2 were positively associated with ERBB2 amplification.
- Geldanamycin, an inhibitor of HSP90 also was positively associated with ERBB2 amplification.
- VX-680 aurora kinase inhibitor
- CGC-11 144 polyamine analogue
- FIG. 10D shows that (a) Homozygous deletion at 9p21 (CDKN2A and CDKN2B) was associated with response to vinorelbine, ixabepilone and fascalypsin. Fascalypsin inhibited CDK4 and this specificity is consistent with the role of the p l6 ,NK4A product of CDKN2A in inhibiting CDK4 20 .
- Example XVI Subtype specificity dominates growth rate effects.
- Example XVII Integration of copy number and transcription measurements identifies pathways of subtype specific responses.
- PARADIGM 24 was used to identify differences in pathway activity among the subtypes in the cell line panel.
- the analysis is complicated by the fact that the curated pathways are partially overlapping. For example EGFR, PI3 kinase and MEK are often curated as separate pathways when in fact they are components of a single larger pathway.
- PARADIGM merges approximately 1400 curated signal transduction, transcriptional and metabolic pathways into a single superimposed pathway (SuperPathway) to eliminate such redundancies.
- SuperPathway superimposed pathway
- PARADIGM uses the pathway interactions to infer integrated pathway levels (IPLs) for every gene, complex, and cellular process.
- Example XVIII Identification of subtype-specific pathway markers.
- the FOXM1 and DNA damage subnetworks also were markedly upregulated in the basal cell lines. Comparison of the claudin-low subtype with all others showed upregulation of many of the same subnetworks as in basal cell lines with some exceptions, including upregulation of the beta-catenin (CTNNB 1 ) network in claudin low cell lines as compared to the basal cells ( Figure 35B). Beta-catennin has been implicated in tumorigenesis, and is associated with poor prognosis.
- ERBB2 AMP cell lines were resistant to the aurora kinase inhibitor VX-680 ( Figure 36C, upper), and further that sensitivity to this compound was not associated with amplification at 20ql 3 (AURKA). This raises the possibility that this resistance may be mediated through CCNB l , which is co-regulated with AUR B by FOXM 1.
- ERBB2 AMP is the only one that shows substantial down-regulation of CCNB 1 ( Figure 36C and Figure 22. This proposed mechanism is supported by the observation that in primary tumors, CCNB l gene expression is significantly correlated with AURKB gene expression.
- Example XVIX Cell growth inhibition assay and growth rate
- GI 50 compound concentration required to inhibit growth by 50%
- TGI Total Growth Inihibition
- LC 50 concentration necessary to reduce the population by 50%
- the drug response data was filtered to meet the following criteria: 1) median standard deviation across the 9 triplicate datapoints ⁇ 0.20; 2) DT +/- 2SD of the median DT for a particular cell line; 3) slope of the fitted curve > 0.25; 4) growth inhibition at the maximum concentration ⁇ 50% for datasets with no clear response. Approximately 80% of the drug plates pass all filtering requirements. We used the median absolute deviation (MAD), a robust version of standard deviation, to assess the reliability of our replicate measures of GI50. Curve fitting and filtering were performed with custom-written R packages.
- MAD median absolute deviation
- Each drug included in the statistical analysis satisfied the following screening criteria for data quality: 1) Missing values: No more than 40% of GI 50 values can be missing across the entire set of cell lines; 2) Variability: For at least 3 cell lines, either GI 50 > 1.5. mGI 50 or GI 50 ⁇ 0.5. mGI 50 , where mGI 50 is the median GI 50 for a given drug. Compounds failing these criteria were excluded from analysis.
- Affymetrix Genome- Wide Human SNP Array 6.0 was used to measure DNA copy number data. The array quality and data processing was performed using the R statistical framework (http://www.r-project.org) based aroma. affymetrix.
- the breast cancer cell line SNP arrays were normalized using 20 normal sample arrays as described (Bengtsson, H., Irizarry, R., Carvalho, B. & Speed, T. P. Estimation and assessment of raw copy numbers at the single locus level. Bioinformatics (Oxford, England) 24, 759-767 (2008)). Data were segmented using circular binary segmentation (CBS) from the bioconductor package DNAcopy (Olshen, A. B.,
- Gene expression data for the cell lines were derived from Affymetrix GeneChip Human Gene 1.0 ST exon arrays. Gene-level summaries of expression were computed using the aroma. affymetrix R package, with quantile normalization and a log-additive probe-level model (PLM) based on the "HuEx- l_0-st-v2,core” chip type. Transcript identifiers were converted to HGNC gene symbols by querying the Ensembl database using the BioMart R package. The resulting expression profiles were subsequently filtered to capture only those genes expressing a standard deviation greater than 1.0 on the log 2 -scale across all cell lines. The raw data are available in ArrayExpress (E-MTAB-181).
- PLM log-additive probe-level model
- Example XXIV Associations of clinically relevant subtypes and response to therapeutic agents
- Example XXV Association of genomic changes and response to therapeutic agents
- Example XXVIII TCGA and cell line clustering
- Pathways were obtained in BioPax Level 2 format from http://pid.nci.nih.gov/ and included NCI-PID, Reactome, and BioCarta databases. Interactions were combined into a merged Superimposed Pathway (SuperPathway). Genes, complexes, and abstract processes (for example, "cell cycle") were retained as pathway concepts. Before merging gene concepts, all gene identifiers were translated into HUGO nomenclature. All interactions were included and no attempt was made to resolve conflicting influences. A breadth-first undirected traversal starting from P53 (the most connected component) was performed to build one single component. The resulting merged pathway structure contained a total of 8768 concepts representing 3491 proteins, 4757 complexes, and 520 processes.
- Expectation-Maximization parameters for PARADIGM were trained on the cell line data and then applied to the TCGA samples. Data from the cell lines and tumor samples were then combined into a single data matrix. Any entry without at least 1 value above 0.5 IPL in either the data from cell lines or tumor samples was removed from further analysis.
- Example XXXI TCGA and cell line clustering
- C s be the set of cell lines of subtype s.
- T s be the set of TCGA tumor samples of subtype s.
- C bllsa i and T huml are the set of all basal cell lines and basal tumor samples respectively.
- the first distribution was made up of t-statistics derived from the Pearson correlations between every possible pair containing a cell line and tumor sample of the same subtype; i.e.
- every pairwise correlation t-statistics was computed between a pair (a, b) such that a e C s and b e T s .
- the second distribution was made of correlation t-statistics between cell lines of different subtypes; that is, computed over pairs (a, b) such that a e C s and b e C 5 and s ⁇ s '.
- Pathways were obtained in BioPax Level 2 format on October 13, 2010 from http://pid.nci.nih.gov/ and included NCI-PID, Reactome, and BioCarta databases. Interactions were combined into a merged Superimposed Pathway (SuperPathway). Genes, complexes, and abstract processes (for example, "cell cycle") were retained as pathway concepts. Before merging gene concepts, all gene identifiers were translated into HUGO nomenclature. All interactions were included and no attempt was made to resolve conflicting influences. A breadth-first undirected traversal starting from P53 (the most connected component) was performed to build one single component. The resulting merged pathway structure contained a total of 8768 concepts representing 3491 proteins, 4757 complexes, and 520 processes.
- Expectation-Maximization parameters for PARADIGM were trained on the cell line data and then applied to the TCGA samples. Data from the cell lines and tumor samples were then combined into a single data matrix. Any entry without at least 1 value above 0.5 IPL in either the data from cell lines or tumor samples was removed from further analysis.
- Example XXXIII TCGA and cell line clustering
- the matrix of inferred pathway activities for both the 47 cell lines and 183 TCGA tumor samples was clustered using complete linkage hierarchical agglomerative clustering implemented in the Eisen Cluster software package version 3.0 45 Uncentered Pearson correlation was used as the metric for the pathway concepts and Euclidean distance was used for sample metric.
- C s be the set of cell lines of subtype s.
- T s be the set of TCGA tumor samples of subtype s.
- C basal and T ba ⁇ l are the set of all basal cell lines and basal tumor samples respectively.
- the first distribution was made up of t-statistics derived from the Pearson correlations between every possible pair containing a cell line and tumor sample of the same subtype; i.e.
- every pairwise correlation t-statistics was computed between a pair (a, b) such that a e C s and b e T s .
- the second distribution was made of correlation t-statistics between cell lines of different subtypes; i.e. computed over pairs (a, b) such that a e C s and b e C 3 - and s ⁇ s '.
- Example XXXIV Molecular subtypes of tumors at various genetic molecular levels.
- Example XXXV PARADIGM for classification of invasive cancers with prognostic significance
- the method uses techniques from probabilistic graphical models (PGM) to integrated functional genomics data onto a known pathway structure. It has previously been applied to analysis of copy number and mRNA expression data from the TCGA glioblastoma and ovarian datasets. PARADIGM analysis can also be used to connect genomic alterations at multiple levels such as DNA methylation or copy number, mRNA and miRNA expression and can thus integrate any number of omics layers of data in each individual sample.
- PGM probabilistic graphical models
- pdgm.4 high ERBB4, low Angiopoietin signaling
- IL4 signaling is strongly down-regulated in conjunction with STAT6, which has been shown in human breast cancer cells to prevent growth inhibition (16 Gooch, J. L., Christy, B., and Yee, D., (2002), STAT6 mediates interleukin-4 growth inhibition in human breast cancer cells, 4: 324-331).
- Down-regulation of IL4 signaling has also promoted mast cell activation which can support greater tumor growth (17 de Visser, . E., Eichten, A., and Coussens, L. M., (2006), Paradoxical roles of the immune system during cancer development, 6: 24-37).
- IL4 is a Th-2 derived cytokine that stimulates B cells differentiation and chronic inflammation in cancer cells. Further Th-2 cells secrete IL10 that mediates immunosuppression in these cancers. This immunosuppression was shown to occur predominantly in basal and ERBB2 cancers. In support to this, it has been shown recently that "antitumor acquired immune programs can be usurped in pro-tumor microenvironments and instead promote malignancy by engaging cellular components of the innate immune system functionally involved in regulating epithelial cell behavior" ( DeNardo, D. G. et al., (2009), CD4(+) T cells regulate pulmonary metastasis of mammary carcinomas by enhancing protumor properties of macrophages, 16: 91-102).
- Example XXXVI Identified pathways whose perturbation specifically influences the PARADIGM clustering.
- FOXM 1 is a key regulator of cell cycle progression and its endogenous FOXM1 expression oscillates according to the phases of the cell cycle.
- FOXM1 confirmed as a human proto-oncogene is found upregulated in the majority of solid human cancers including liver, breast, lung, prostate, cervix of uterus, colon, pancreas, brain as well as basal cell carcinoma, the most common human cancer.
- FOXM1 is thought to promote oncogenesis through its multiple roles in cell cycle and chromosomal/genomic maintenance (Wonsey, D. R. and Follettie, M.
- Figure 26 illustrates the opposite activation modus of this pathway (red as activated vs blue inactivated) for cluster pdgm 3 (best survival) as opposed to the rest of the clusters with worse survival and the molecular levels that contribute to it (mRNA, CNA, miRNA or DNA methylation according to the shape of the figures).
- Angiopoietin receptor tie2-mediated signaling is angiopoietin receptor tie2-mediated signaling.
- Ang2 plays an important role in angiogenesis during the development and growth of human cancers.
- Ang2's role in angiogenesis generally is considered as an antagonist for Angl , inhibiting Angl -promoted Tie2 signaling, which is critical for blood vessel maturation and stabilization(23).
- Ang2 modulates angiogenesis in a cooperative manner with another important angiogenic factor, vascular endothelial growth factor A (VEGFA) (Hashizume, H. et al., (2010), Complementary actions of inhibitors of angiopoietin-2 and VEGF on tumor angiogenesis and growth, 70: 2213-2223).
- VAGFA vascular endothelial growth factor A
- Ang2 angiopoietin family members can activate Tie l , for example, Angl induces Tiel phosphorylation in endothelial cells (2 Yuan, H. T. et al., (2007), Activation of the orphan endothelial receptor Tiel modifies Tie2-mediated intracellular signaling and cell survival, 21 : 3171-3183).
- Tiel phosphorylation is, however, Tie2 dependent because Angl fails to induce Tiel phosphorylation when Tie2 is down-regulated in endothelial cells and Tiel phosphorylation is induced in the absence of Angl by either a constitutively active form of Tie2 or a Tie2 agonistic antibody (25 Yuan et al. (2007) supra).
- Angl-mediated AKT and 42/44MAPK phosphorylation is predominantly Tie2 mediated, and Tie l down- regulates this pathway.
- the main role for Tiel is to modulate blood vessel morphogenesis due to its ability to down-regulate Tie2-driven signaling and endothelial survival. Both Tie2 mediated signaling as well as VEGFR 1 and 2 mediated signaling and specific signals were observed in this dataset.
- ERBB4 ERBB4
- ERBB4 contributes to proliferation and cell movements in mammary morphogenesis and the directional cell movements of Erbb4-expressing mammary primordial epithelia while promoting mammary cell fate.
- Candidate effectors of Nrg3/Erbb4 signaling have been identified and shown here to interacts with other signalling pathways relevant to early mammary gland development and cancer.
- One of the primary functions of ErbB4 in vivo is in the maturation of mammary glands during pregnancy and lactation induction. Pregnancy and extended lactation durations have been correlated with reduced risk of breast cancer, and the role of ErbB4 in tumor suppression may therefore be linked with its role in lactation.
- Example XXXVII PARADIGM for classification in ductal carcinoma in situ (DCIS)
- DCIS Ductal carcinoma in situ
- TCR signaling in naive CD8+ T cells was on top of the list alongside with a large number of chemokines that are known to recruit CD8+ T cells.
- One is IL- 12, produced by the antigen presenting cells that was shown to stimulate IFN-gamma production from NK and T cells.
- IFN-gamma pathway was one of the deregulated pathways, higher up on the list in DCIS. IFNgamma is produced from the Th l cells and the NK cells and was shown to initiate an antitumor immune response.
- trastuzumab (herceptin) is potentiated by the co-administration of EL- 12 to patients with HER2-overexpressing tumors, and this effect is mediated by the stimulation of IFNgamma production in the NK cells (29).
- other most strong contributor (Table 8) was 84_NOX4.
- NOX4 an oxygen-sensing NAPHD oxidase, and a phagocyte-type A oxidase, is similar to that responsible for the production of large amounts of reactive oxygen species (ROS) in neutrophil granulocytes, primary immune response.
- ROS reactive oxygen species
- FNl fibronectin
- PDGFRB the platelet-derived growth factor receptor
- FNl Fibronectin- 1 belongs to a family of high molecular weight glycoproteins that are present on cell surfaces, in extracellular fluids, connective tissues, and basement membranes. Fibronectins interact with other extracellular matrix proteins and cellular ligands, such as collagen, fibrin, and integrins. Fibronectins are involved in adhesive and migratory processes of cells.
- PDGFR the platelet- derived growth factor receptor
- EGF Epidermal growth factor
- PDGF receptors which are important receptor tyrosine kinases (RTKs).
- RTKs receptor tyrosine kinases
- CNA copy number alterations
- miRNA profiling from total RNA was performed using Agilent Technologies "Human miRNA Microarray Kit (V2)" according to manufacturer's protocol. Scanning on Agilent Scanner G2565A and Feature Extraction (FE) v9.5 was used to extract signals. Experiments were performed using duplicate hybridizations (99 samples) on different arrays and time points. Two samples were profiled only once. miRNA signal intensities for replicate probes were averaged across the platform, log2 transformed and normalized to the 75 percentile. miRNA expression status was scored as present or absent for each gene in each sample by default settings in FE v9.5.
- DNA methylation One microgram of DNA was bisulphite treated using the EpiTect 96 Bisulfite Kit (Qiagen GmbH, Germany). 500 ng of bisulphite treated DNA was analyzed using the GoldenGate Methylation Cancer Panel I ( lumina Inc, CA, USA) that simultaneously analyses 1505 CpG sites in 807 cancer related genes. At least 2 CpG sites were analyzed per gene were one CpG site is in the promoter region and one CpG site is in the 1 st exon Bead studio software was used for the initial processing of the methylation data according to the manufacturer's protocol. The detection p- value for each CpG site was used to validate sample performance and the dataset was filtered based on the detection p- value were CpG sites with a detection p-value> 0.05 was omitted from further analysis.
- Copy number was segmented using CBS, then mapped to gene-level measurements by taking the median of all segments that span a RefSeq gene's coordinates in hgl8.
- measurements were first probe-normalized by subtracting the median expression value for each probe. The manufacturer's genomic location for each probe was converted from hgl 7 to hgl 8 using UCSCs liftOver tool. Per-gene measurements were then obtained by taking the median value of all probes overlapping a RefSeq gene. Methylation probes were matched to genes using manufacturers description.
- Pathway files were from the PID (36) as previously parsed.
- Figure 26 shows summaries of discretized input data, and not IPL values, by counting the fraction of observations in either an up or down bin in each datatype, and then labeling each node with the bin with the highest fraction of observations in any datatype.
- Clusters were derived using the HOPACH R implementation version 2.10 (37) running on R version 2.12.
- the correlation distance metric was used with all data types, except for Paradigm IPLs, which used cosangle due to the non-normal distribution and prevalence of zero values. For any cluster of samples that contained fewer than 5 samples, each sample was mapped to the same cluster as the most similar sample in a larger cluster.
- Paradigm clusters in the MicMa dataset were mapped to other datatypes by determining each cluster's mediod (using the median function) in the MicMa dataset, then assigning each sample in another dataset to whichever cluster mediod was closest by cosangle distance.
- Kaplain-Meier, Cluster enrichments Kaplan-Meier statistics, plots, and cluster enrichments were determined using R version 2.12. Cox p-values were determined using the Wald test from the coxph() proportional hazards model, and log-rank p-values from a chi-square test from the survdiff() function. Overall enrichment of a gene's or pathway member's values for a clustering were determined by ANOVA, and enrichment of a gene for a particular cluster label were determined by a T-test of a gene's values in a particular cluster vs. the gene's values in all other clusters. FDR was determined using the Benjamini &Hochberg method of p. adjust.
- a collection of pathways was obtained from NCI-PID containing 131 pathways, 11 ,563 interactions, and 7,204 entities.
- An entity is molecule, complex, small molecule, or abstract concept represented as "nodes” in PARADIGM' S graphical model.
- the abstract concepts correspond to general cellular processes (such as "apoptosis” or "absorption of light,") and families of genes that share functional activity such as the RAS family of signal transducers.
- We collected interactions including protein-protein interactions, transcriptional regulatory interactions, protein modifications such as phosphorylation and ubiquitinylation interactions.
- Example XL Inference of integrated molecular activities in pathway context.
- PARADIGM which assigns an integrated pathway activity (IPA) reflecting the copy number, gene expression, and pathway context of each entity.
- IPAs The significance of IPAs was assessed using permutations of gene- and patient-specific cross-sections of data. Data for 1000 "null" patients was created by randomly selecting a gene-expression and copy number pair of values for each gene in the genome. To assess the significance of the PARADIGM IPAs, we constructed a null distribution by assigning random genes to pathways while preserving the pathway structure.
- Example XLIV Identification of SNPs [00316] A region of a gene from a patient's DNA sample is amplified by PCR using the primers specifically designed for the region. The PCR products are sequenced using methods well known to those of skill in the art, as disclosed above. SNPs identified in the sequence traces are verified using Phred/Phrap/Consed software and compared with known SNPs deposited in the NCBI SNP databank.
- Linguistics Georgetown University, Washington DC is used to assess differences between genotype frequencies in normal subjects and patients with a disorder.
- One-way ANOVA with post-hoc analysis is performed as indicated to compare hemodynamics between different patient groups.
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Evolutionary Biology (AREA)
- Biotechnology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Genetics & Genomics (AREA)
- Physiology (AREA)
- Data Mining & Analysis (AREA)
- Public Health (AREA)
- Databases & Information Systems (AREA)
- Probability & Statistics with Applications (AREA)
- Pathology (AREA)
- Primary Health Care (AREA)
- Biomedical Technology (AREA)
- Epidemiology (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
Description
Claims
Priority Applications (9)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CA2853702A CA2853702A1 (en) | 2011-10-26 | 2011-10-31 | Pathway recognition algorithm using data integration on genomic models (paradigm) |
KR1020157010056A KR20150048909A (en) | 2011-10-26 | 2011-10-31 | Pathway recognition algorithm using data integration on genomic models (paradigm) |
CN201180075918.7A CN104011726A (en) | 2011-10-26 | 2011-10-31 | Pathway recognition algorithm using data integration on genomic models (PARADIGM) |
JP2014538758A JP6061941B2 (en) | 2011-10-26 | 2011-10-31 | Pathway recognition algorithm (PARADIGM) using data integration for genome models |
AU2011380013A AU2011380013B2 (en) | 2011-10-26 | 2011-10-31 | Pathway recognition algorithm using data integration on genomic models (paradigm) |
EP11874850.8A EP2771830A4 (en) | 2011-10-26 | 2011-10-31 | Pathway recognition algorithm using data integration on genomic models (paradigm) |
KR1020147013980A KR101554530B1 (en) | 2011-10-26 | 2011-10-31 | Pathway recognition algorithm using data integration on genomic models (paradigm) |
IL232243A IL232243A (en) | 2011-10-26 | 2014-04-24 | Pathway recognition algorithm using data integration on genomic models (paradigm) |
IL248380A IL248380B (en) | 2011-10-26 | 2016-10-18 | Pathway recognition algorithm using data integration on genomic models (paradigm) |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/317,769 | 2011-10-26 | ||
US13/317,769 US10192641B2 (en) | 2010-04-29 | 2011-10-26 | Method of generating a dynamic pathway map |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013062505A1 true WO2013062505A1 (en) | 2013-05-02 |
Family
ID=48168183
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2011/001844 WO2013062505A1 (en) | 2011-10-26 | 2011-10-31 | Pathway recognition algorithm using data integration on genomic models (paradigm) |
Country Status (9)
Country | Link |
---|---|
US (2) | US10192641B2 (en) |
EP (1) | EP2771830A4 (en) |
JP (4) | JP6061941B2 (en) |
KR (2) | KR20150048909A (en) |
CN (2) | CN108664762A (en) |
AU (1) | AU2011380013B2 (en) |
CA (2) | CA2853702A1 (en) |
IL (2) | IL232243A (en) |
WO (1) | WO2013062505A1 (en) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014059036A1 (en) | 2012-10-09 | 2014-04-17 | Five3 Genomics, Llc | Systems and methods for learning and identification of regulatory interactions in biological pathways |
WO2014193982A1 (en) | 2013-05-28 | 2014-12-04 | Five3 Genomics, Llc | Paradigm drug response networks |
WO2014210611A1 (en) | 2013-06-28 | 2014-12-31 | Nantomics, Llc | Pathway analysis for identification of diagnostic tests |
WO2015060994A2 (en) | 2013-07-26 | 2015-04-30 | Nant Holdings Ip, Llc | Discovery routing systems and engines |
EP3149641A4 (en) * | 2014-05-30 | 2018-01-24 | Nantomics, LLC | Systems and methods for comprehensive analysis of molecular profiles across multiple tumor and germline exomes |
US20180046751A1 (en) * | 2010-04-29 | 2018-02-15 | The Regents Of The University Of California | Pathway recognition algorithm using data integration on genomic models (paradigm) |
US9953137B2 (en) | 2012-07-06 | 2018-04-24 | Nant Holdings Ip, Llc | Healthcare analysis stream management |
WO2018093658A1 (en) | 2016-11-17 | 2018-05-24 | Nantbio, Inc. | Validation of inferred anticancer pathways |
WO2018148381A1 (en) | 2017-02-07 | 2018-08-16 | Nantcell, Inc. | Maximizing t-cell memory and compositions and methods therefor |
US10192027B2 (en) | 2014-03-25 | 2019-01-29 | Five3 Genomics, Llc | Systems and methods for RNA analysis in functional confirmation of cancer mutations |
CN109416925A (en) * | 2016-05-05 | 2019-03-01 | 南托米克斯有限责任公司 | It checks point failure and makes the method for checking point failure |
US10532089B2 (en) | 2015-10-12 | 2020-01-14 | Nantomics, Llc | Iterative discovery of neoepitopes and adaptive immunotherapy and methods therefor |
US20200292515A1 (en) * | 2018-03-01 | 2020-09-17 | Obschestvo s ogranichennoi otvetstvennostyu "Oncobox" | Genetic Information Analysis Platform Oncobox |
US10984890B2 (en) | 2016-06-30 | 2021-04-20 | Nantomics, Llc | Synthetic WGS bioinformatics validation |
US11071774B2 (en) | 2016-06-30 | 2021-07-27 | Nantcell, Inc. | Nant cancer vaccine |
US11101038B2 (en) | 2015-01-20 | 2021-08-24 | Nantomics, Llc | Systems and methods for response prediction to chemotherapy in high grade bladder cancer |
US11154597B2 (en) | 2016-03-24 | 2021-10-26 | Nantcell, Inc. | Sequence arrangements and sequences for neoepitope presentation |
US11361844B2 (en) | 2013-09-26 | 2022-06-14 | Five3 Genomics, Llc | Systems, methods, and compositions for viral-associated tumors |
US11564944B2 (en) | 2016-11-21 | 2023-01-31 | Nant Holdings Ip, Llc | Fractal combination therapy |
US11564980B2 (en) | 2018-04-23 | 2023-01-31 | Nantcell, Inc. | Tumor treatment method with an individualized peptide vaccine |
US11585805B2 (en) | 2016-02-19 | 2023-02-21 | Nantcell, Inc. | Methods of immunogenic modulation |
US11590217B2 (en) | 2018-04-23 | 2023-02-28 | Nantcell, Inc. | Neoepitope vaccine and immune stimulant combinations and methods |
US11773187B2 (en) | 2018-01-17 | 2023-10-03 | Nantbio, Inc. | Enhanced immunogenicity for GPI-anchored antigens |
US11779637B2 (en) | 2017-04-24 | 2023-10-10 | Nantcell, Inc. | Targeted neoepitope vectors and methods therefor |
US11823773B2 (en) | 2018-04-13 | 2023-11-21 | Nant Holdings Ip, Llc | Nant cancer vaccine strategies |
US11875877B2 (en) | 2016-08-25 | 2024-01-16 | Nantomics, Llc | Immunotherapy markers and uses therefor |
Families Citing this family (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011079034A1 (en) | 2011-07-12 | 2013-01-17 | Siemens Aktiengesellschaft | Control of a technical system |
US11361867B2 (en) * | 2012-10-05 | 2022-06-14 | H. Lee Moffitt Cancer Center And Research Institute, Inc. | Pathways for treating patients |
AU2013368945B2 (en) | 2012-12-26 | 2020-01-23 | Innosign B.V. | Assessment of cellular signaling pathway activity using linear combination(s) of target gene expressions |
CA2989815A1 (en) * | 2015-06-15 | 2016-12-22 | Nantomics, Llc | Systems and methods for patient-specific prediction of drug responses from cell line genomics |
KR102375973B1 (en) * | 2015-11-05 | 2022-03-17 | 김훈 | Security server using case based reasoning engine and storage medium for installing security function |
CN105447337B (en) * | 2015-11-13 | 2018-01-26 | 大连理工大学 | A kind of time series data processing method based on dynamic network map analysis |
WO2017191076A1 (en) | 2016-05-01 | 2017-11-09 | Genome Research Limited | Method of characterising a dna sample |
US20170329914A1 (en) * | 2016-05-11 | 2017-11-16 | International Business Machines Corporation | Predicting Personalized Cancer Metastasis Routes, Biological Mediators of Metastasis and Metastasis Blocking Therapies |
WO2018027076A1 (en) * | 2016-08-03 | 2018-02-08 | Nantomics, Llc | Dasatinib response prediction models and methods therefor |
EP3535422A2 (en) * | 2016-11-07 | 2019-09-11 | Grail, Inc. | Methods of identifying somatic mutational signatures for early cancer detection |
WO2018094204A1 (en) * | 2016-11-17 | 2018-05-24 | Arivale, Inc. | Determining relationships between risks for biological conditions and dynamic analytes |
US10629292B2 (en) | 2016-12-19 | 2020-04-21 | Seven Bridges Genomics, Inc. | Generation and use of simulated genomic data |
US20180251849A1 (en) * | 2017-03-03 | 2018-09-06 | General Electric Company | Method for identifying expression distinguishers in biological samples |
AU2018244307A1 (en) * | 2017-03-28 | 2019-10-24 | Nant Holdings Ip, Llc | Modeling miRNA induced silencing in breast cancer with paradigm |
US10686669B2 (en) * | 2017-06-16 | 2020-06-16 | Cisco Technology, Inc. | Collecting network models and node information from a network |
WO2019178546A1 (en) | 2018-03-16 | 2019-09-19 | Scipher Medicine Corporation | Methods and systems for predicting response to anti-tnf therapies |
CN108795934B (en) * | 2018-05-23 | 2022-06-21 | 中国农业科学院上海兽医研究所(中国动物卫生与流行病学中心上海分中心) | SiRNA of schistosoma japonicum SjELAV-like 2 gene and application thereof |
US11189361B2 (en) * | 2018-06-28 | 2021-11-30 | International Business Machines Corporation | Functional analysis of time-series phylogenetic tumor evolution tree |
US11211148B2 (en) | 2018-06-28 | 2021-12-28 | International Business Machines Corporation | Time-series phylogenetic tumor evolution trees |
JP7404648B2 (en) | 2019-04-25 | 2023-12-26 | 富士通株式会社 | Therapeutic drug presentation method, therapeutic drug presentation device, and therapeutic drug presentation program |
JP2022541125A (en) | 2019-06-27 | 2022-09-22 | サイファー メディシン コーポレーション | Developing classifiers to stratify patients |
CN112816275B (en) * | 2019-11-18 | 2022-04-05 | 中国科学院深圳先进技术研究院 | Brain repair device |
CN111161796B (en) * | 2019-12-30 | 2024-04-16 | 中南大学 | Method and system for predicting PD potential gene and miRNA |
RU2756883C1 (en) * | 2020-06-09 | 2021-10-06 | ФГАУ "Военный инновационный технополис "ЭРА" | Apparatus for probabilistic modelling of the process of functioning of a telecommunication network |
CN111883203B (en) * | 2020-07-03 | 2023-12-29 | 上海厦维医学检验实验室有限公司 | Construction method of model for predicting PD-1 curative effect |
US11481235B2 (en) * | 2021-01-11 | 2022-10-25 | Evicore Healthcare MSI, LLC | Database framework model transformation for pathway identification |
US11439159B2 (en) | 2021-03-22 | 2022-09-13 | Shiru, Inc. | System for identifying and developing individual naturally-occurring proteins as food ingredients by machine learning and database mining combined with empirical testing for a target food function |
CN113434693B (en) * | 2021-06-23 | 2023-02-21 | 重庆邮电大学工业互联网研究院 | Data integration method based on intelligent data platform |
CN113436729A (en) * | 2021-07-08 | 2021-09-24 | 湖南大学 | Synthetic lethal interaction prediction method based on heterogeneous graph convolution neural network |
CN115203172B (en) * | 2022-06-30 | 2023-11-07 | 北京亚控科技发展有限公司 | Model construction and model data subscription method and device, electronic equipment and medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001050950A2 (en) * | 2000-01-11 | 2001-07-19 | Zycare, Inc. | Apparatus and methods for monitoring and modifying anticoagulation therapy of remotely located patients |
US20030113761A1 (en) * | 2001-08-16 | 2003-06-19 | Patrick Tan | Modelling of biochemical pathways |
US20040167763A1 (en) * | 2002-05-31 | 2004-08-26 | Liebman Michael N | Information processing method for evaluating biochemical pathway models using clinical data |
US20100100331A1 (en) * | 2003-08-01 | 2010-04-22 | Dna Twopointo Inc. | Systems and methods for biopolymer engineering |
US20110093244A1 (en) * | 2005-04-13 | 2011-04-21 | Venture Gain LLC | Analysis of Transcriptomic Data Using Similarity Based Modeling |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5968770A (en) | 1995-02-10 | 1999-10-19 | Millennium Pharmaceuticals, Inc. | Compositions and methods for the treatment and diagnosis of cardiovascular disease using rchd523 as a target |
US7709454B2 (en) | 1997-06-20 | 2010-05-04 | New York University | Methods and compositions for inhibiting tumorigenesis |
AU2001229744A1 (en) | 2000-01-25 | 2001-08-07 | Cellomics, Inc. | Method and system for automated inference of physico-chemical interaction knowl edge |
EP1739584A4 (en) | 2004-03-30 | 2008-07-23 | Shigeo Ihara | Document information processing system |
JP2006185412A (en) | 2004-12-03 | 2006-07-13 | Kazusa Dna Kenkyusho | Information processor, information processing method and program thereof |
US7848890B2 (en) | 2004-12-08 | 2010-12-07 | Electronics And Telecommunications Research Institute | Method and system for predicting gene pathway using gene expression pattern data and protein interaction data |
EP1910564A1 (en) | 2005-05-13 | 2008-04-16 | Duke University | Gene expression signatures for oncogenic pathway deregulation |
US8642512B2 (en) * | 2005-05-17 | 2014-02-04 | Cellact Pharma Gmbh | Method of detecting early immune activation |
JP2007011996A (en) | 2005-07-04 | 2007-01-18 | Fujitsu Ltd | Analysis method for expression information and system therefor |
JP2007052766A (en) * | 2005-07-22 | 2007-03-01 | Mathematical Systems Inc | Pathway display method, information processing device, and pathway display program |
WO2008115497A2 (en) * | 2007-03-16 | 2008-09-25 | Gene Security Network | System and method for cleaning noisy genetic data and determining chromsome copy number |
US10262103B2 (en) | 2008-11-18 | 2019-04-16 | Raphael LEHRER | Individualized cancer treatment |
JP5429567B2 (en) | 2009-06-29 | 2014-02-26 | 株式会社ジェイテクト | Resin cage for angular contact ball bearings |
GB2474694B (en) * | 2009-10-23 | 2011-11-02 | Innovia Films Ltd | Biodegradable composites |
JP5670911B2 (en) * | 2009-11-09 | 2015-02-18 | 大日本住友製薬株式会社 | Support device, support method, and computer program |
CA2796272C (en) | 2010-04-29 | 2019-10-01 | The Regents Of The University Of California | Pathway recognition algorithm using data integration on genomic models (paradigm) |
-
2011
- 2011-10-26 US US13/317,769 patent/US10192641B2/en active Active
- 2011-10-31 KR KR1020157010056A patent/KR20150048909A/en not_active Application Discontinuation
- 2011-10-31 WO PCT/US2011/001844 patent/WO2013062505A1/en active Application Filing
- 2011-10-31 CN CN201810729146.8A patent/CN108664762A/en active Pending
- 2011-10-31 JP JP2014538758A patent/JP6061941B2/en active Active
- 2011-10-31 EP EP11874850.8A patent/EP2771830A4/en not_active Ceased
- 2011-10-31 CN CN201180075918.7A patent/CN104011726A/en active Pending
- 2011-10-31 AU AU2011380013A patent/AU2011380013B2/en active Active
- 2011-10-31 CA CA2853702A patent/CA2853702A1/en not_active Abandoned
- 2011-10-31 CA CA3021833A patent/CA3021833A1/en not_active Abandoned
- 2011-10-31 KR KR1020147013980A patent/KR101554530B1/en active IP Right Grant
-
2014
- 2014-04-24 IL IL232243A patent/IL232243A/en active IP Right Grant
- 2014-12-19 US US14/577,522 patent/US20150142465A1/en not_active Abandoned
-
2016
- 2016-01-04 JP JP2016000304A patent/JP6129998B2/en active Active
- 2016-10-18 IL IL248380A patent/IL248380B/en active IP Right Grant
-
2017
- 2017-04-12 JP JP2017078768A patent/JP6396532B2/en active Active
-
2018
- 2018-08-29 JP JP2018160366A patent/JP6854792B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001050950A2 (en) * | 2000-01-11 | 2001-07-19 | Zycare, Inc. | Apparatus and methods for monitoring and modifying anticoagulation therapy of remotely located patients |
US20030113761A1 (en) * | 2001-08-16 | 2003-06-19 | Patrick Tan | Modelling of biochemical pathways |
US20040167763A1 (en) * | 2002-05-31 | 2004-08-26 | Liebman Michael N | Information processing method for evaluating biochemical pathway models using clinical data |
US20100100331A1 (en) * | 2003-08-01 | 2010-04-22 | Dna Twopointo Inc. | Systems and methods for biopolymer engineering |
US20110093244A1 (en) * | 2005-04-13 | 2011-04-21 | Venture Gain LLC | Analysis of Transcriptomic Data Using Similarity Based Modeling |
Non-Patent Citations (1)
Title |
---|
See also references of EP2771830A4 * |
Cited By (58)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180046751A1 (en) * | 2010-04-29 | 2018-02-15 | The Regents Of The University Of California | Pathway recognition algorithm using data integration on genomic models (paradigm) |
US11355218B2 (en) | 2010-04-29 | 2022-06-07 | Thf Regents Of The University Of California | Patient-specific cellular pathway activity inference computer system |
US10916329B2 (en) | 2010-04-29 | 2021-02-09 | The Regents Of The University Of California | Method for treating an individual patient |
US11043282B2 (en) | 2010-04-29 | 2021-06-22 | The Regents Of The Unviersity Of California | Pathway analysis computing system and method |
US10991448B2 (en) | 2010-04-29 | 2021-04-27 | The Regents Of The University Of California | Pathway recognition algorithm using data integration on genomic models (paradigm) |
US20190341123A1 (en) * | 2010-04-29 | 2019-11-07 | The Regents Of The University Of California | Pathway recognition algorithm using data integration on genomic models (paradigm) |
US10770169B2 (en) | 2010-04-29 | 2020-09-08 | The Regents Of The University Of California | Method of providing a treatment to a subject based on a dynamic pathway map |
US10957429B2 (en) | 2012-07-06 | 2021-03-23 | Nant Holdings Ip, Llc | Healthcare analysis stream management |
US9953137B2 (en) | 2012-07-06 | 2018-04-24 | Nant Holdings Ip, Llc | Healthcare analysis stream management |
US10580523B2 (en) | 2012-07-06 | 2020-03-03 | Nant Holdings Ip, Llc | Healthcare analysis stream management |
US10055546B2 (en) | 2012-07-06 | 2018-08-21 | Nant Holdings Ip, Llc | Healthcare analysis stream management |
US10095835B2 (en) | 2012-07-06 | 2018-10-09 | Nant Holdings Ip, Llc | Healthcare analysis stream management |
WO2014059036A1 (en) | 2012-10-09 | 2014-04-17 | Five3 Genomics, Llc | Systems and methods for learning and identification of regulatory interactions in biological pathways |
EP3471103A1 (en) | 2012-10-09 | 2019-04-17 | Five3 Genomics, LLC | Systems and methods for learning and identification of regulatory interactions in biological pathways |
JP2018022511A (en) * | 2013-05-28 | 2018-02-08 | ファイヴ3 ゲノミクス,エルエルシー | Paradigm drug response networks |
EP3564958A1 (en) | 2013-05-28 | 2019-11-06 | Five3 Genomics, LLC | Paradigm drug response networks |
AU2014274224B2 (en) * | 2013-05-28 | 2016-06-09 | Five3 Genomics, Llc | Paradigm drug response networks |
WO2014193982A1 (en) | 2013-05-28 | 2014-12-04 | Five3 Genomics, Llc | Paradigm drug response networks |
US11011273B2 (en) | 2013-06-28 | 2021-05-18 | Nantomics, Llc | Pathway analysis for identification of diagnostic tests |
WO2014210611A1 (en) | 2013-06-28 | 2014-12-31 | Nantomics, Llc | Pathway analysis for identification of diagnostic tests |
US11017884B2 (en) | 2013-07-26 | 2021-05-25 | Nant Holdings Ip, Llc | Discovery routing systems and engines |
US10114925B2 (en) | 2013-07-26 | 2018-10-30 | Nant Holdings Ip, Llc | Discovery routing systems and engines |
US10332618B2 (en) | 2013-07-26 | 2019-06-25 | Nant Holdings Ip, Llc | Discovery routing systems and engines |
US12051485B2 (en) | 2013-07-26 | 2024-07-30 | Nant Holdings Ip, Llc | Discovery routing systems and engines |
WO2015060994A2 (en) | 2013-07-26 | 2015-04-30 | Nant Holdings Ip, Llc | Discovery routing systems and engines |
US11361844B2 (en) | 2013-09-26 | 2022-06-14 | Five3 Genomics, Llc | Systems, methods, and compositions for viral-associated tumors |
US10192027B2 (en) | 2014-03-25 | 2019-01-29 | Five3 Genomics, Llc | Systems and methods for RNA analysis in functional confirmation of cancer mutations |
EP3149641A4 (en) * | 2014-05-30 | 2018-01-24 | Nantomics, LLC | Systems and methods for comprehensive analysis of molecular profiles across multiple tumor and germline exomes |
US10614910B2 (en) | 2014-05-30 | 2020-04-07 | Nantomics, Llc | Systems and methods for comprehensive analysis of molecular profiles across multiple tumor and germline exomes |
US11101038B2 (en) | 2015-01-20 | 2021-08-24 | Nantomics, Llc | Systems and methods for response prediction to chemotherapy in high grade bladder cancer |
US10532089B2 (en) | 2015-10-12 | 2020-01-14 | Nantomics, Llc | Iterative discovery of neoepitopes and adaptive immunotherapy and methods therefor |
US11717564B2 (en) | 2015-10-12 | 2023-08-08 | Nantomics, Llc | Iterative discovery of neoepitopes and adaptive immunotherapy and methods therefor |
EP3855444A1 (en) | 2015-10-12 | 2021-07-28 | Nantomics, LLC | Iterative discovery of neoepitopes and adaptive immunotherapy and methods therefor |
US11585805B2 (en) | 2016-02-19 | 2023-02-21 | Nantcell, Inc. | Methods of immunogenic modulation |
US11154597B2 (en) | 2016-03-24 | 2021-10-26 | Nantcell, Inc. | Sequence arrangements and sequences for neoepitope presentation |
CN109416925A (en) * | 2016-05-05 | 2019-03-01 | 南托米克斯有限责任公司 | It checks point failure and makes the method for checking point failure |
EP3452936A4 (en) * | 2016-05-05 | 2020-01-15 | Nantomics, LLC | Checkpoint failure and methods therefor |
US11071774B2 (en) | 2016-06-30 | 2021-07-27 | Nantcell, Inc. | Nant cancer vaccine |
US10984890B2 (en) | 2016-06-30 | 2021-04-20 | Nantomics, Llc | Synthetic WGS bioinformatics validation |
US11439697B2 (en) | 2016-06-30 | 2022-09-13 | Nantcell, Inc. | Nant cancer vaccine |
US11207392B2 (en) | 2016-06-30 | 2021-12-28 | Nantcell Inc. | Coordinated multi-prong cancer therapy |
EP4101447A1 (en) | 2016-06-30 | 2022-12-14 | Nant Holdings IP, LLC | Nant cancer vaccine |
US11875877B2 (en) | 2016-08-25 | 2024-01-16 | Nantomics, Llc | Immunotherapy markers and uses therefor |
WO2018093658A1 (en) | 2016-11-17 | 2018-05-24 | Nantbio, Inc. | Validation of inferred anticancer pathways |
AU2017360887B2 (en) * | 2016-11-17 | 2023-01-19 | Nant Holdings Ip, Llc | Validation of inferred anticancer pathways |
JP2020506450A (en) * | 2016-11-17 | 2020-02-27 | ナントバイオ,インコーポレイテッド | Validation of inferred anticancer pathway |
KR102608922B1 (en) | 2016-11-17 | 2023-12-01 | 난트바이오 인코포레이티드 | Validation of the inferred anticancer pathway |
KR20190090805A (en) * | 2016-11-17 | 2019-08-02 | 난트바이오 인코포레이티드 | Validation of Inferred Anticancer Pathways |
US11564944B2 (en) | 2016-11-21 | 2023-01-31 | Nant Holdings Ip, Llc | Fractal combination therapy |
US11229668B2 (en) | 2017-02-07 | 2022-01-25 | Nantcell, Inc. | Maximizing T-cell memory and compositions and methods therefor |
WO2018148381A1 (en) | 2017-02-07 | 2018-08-16 | Nantcell, Inc. | Maximizing t-cell memory and compositions and methods therefor |
US11779637B2 (en) | 2017-04-24 | 2023-10-10 | Nantcell, Inc. | Targeted neoepitope vectors and methods therefor |
US11773187B2 (en) | 2018-01-17 | 2023-10-03 | Nantbio, Inc. | Enhanced immunogenicity for GPI-anchored antigens |
US11614434B2 (en) * | 2018-03-01 | 2023-03-28 | Micsway Corp. | Genetic information analysis platform oncobox |
US20200292515A1 (en) * | 2018-03-01 | 2020-09-17 | Obschestvo s ogranichennoi otvetstvennostyu "Oncobox" | Genetic Information Analysis Platform Oncobox |
US11823773B2 (en) | 2018-04-13 | 2023-11-21 | Nant Holdings Ip, Llc | Nant cancer vaccine strategies |
US11590217B2 (en) | 2018-04-23 | 2023-02-28 | Nantcell, Inc. | Neoepitope vaccine and immune stimulant combinations and methods |
US11564980B2 (en) | 2018-04-23 | 2023-01-31 | Nantcell, Inc. | Tumor treatment method with an individualized peptide vaccine |
Also Published As
Publication number | Publication date |
---|---|
IL232243A0 (en) | 2014-06-30 |
AU2011380013B2 (en) | 2014-06-05 |
JP6396532B2 (en) | 2018-09-26 |
KR20140078766A (en) | 2014-06-25 |
JP2019008812A (en) | 2019-01-17 |
CA2853702A1 (en) | 2013-05-02 |
US20120158391A1 (en) | 2012-06-21 |
US10192641B2 (en) | 2019-01-29 |
CA3021833A1 (en) | 2013-05-02 |
CN104011726A (en) | 2014-08-27 |
KR101554530B1 (en) | 2015-09-21 |
AU2011380013A1 (en) | 2014-05-08 |
JP2016122451A (en) | 2016-07-07 |
CN108664762A (en) | 2018-10-16 |
JP6129998B2 (en) | 2017-05-17 |
JP6061941B2 (en) | 2017-01-18 |
JP6854792B2 (en) | 2021-04-07 |
JP2014532859A (en) | 2014-12-08 |
IL248380A0 (en) | 2016-11-30 |
IL232243A (en) | 2016-11-30 |
IL248380B (en) | 2019-07-31 |
EP2771830A4 (en) | 2015-03-25 |
KR20150048909A (en) | 2015-05-07 |
US20150142465A1 (en) | 2015-05-21 |
JP2017182813A (en) | 2017-10-05 |
EP2771830A1 (en) | 2014-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10916329B2 (en) | Method for treating an individual patient | |
JP6854792B2 (en) | Pathway Recognition Algorithm Using Data Integration for Genome Models (PARADIGM) | |
AU2014221248B2 (en) | Pathway recognition algorithm using data integration on genomic models (paradigm) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11874850 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 232243 Country of ref document: IL |
|
ENP | Entry into the national phase |
Ref document number: 2853702 Country of ref document: CA |
|
ENP | Entry into the national phase |
Ref document number: 2014538758 Country of ref document: JP Kind code of ref document: A |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2011380013 Country of ref document: AU Date of ref document: 20111031 Kind code of ref document: A |
|
ENP | Entry into the national phase |
Ref document number: 20147013980 Country of ref document: KR Kind code of ref document: A |
|
REEP | Request for entry into the european phase |
Ref document number: 2011874850 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011874850 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 248380 Country of ref document: IL |