Cancer progression modeling using static sample data.
Ontology highlight
ABSTRACT:
SUBMITTER: Sun Y
PROVIDER: S-EPMC4196119 | biostudies-literature | 2014 Aug
REPOSITORIES: biostudies-literature
Publications
Cancer progression modeling using static sample data.
Sun Yijun Y Yao Jin J Nowak Norma J NJ Goodison Steve S
Genome biology 20140826 8
As molecular profiling data continues to accumulate, the design of integrative computational analyses that can provide insights into the dynamic aspects of cancer progression becomes feasible. Here, we present a novel computational method for the construction of cancer progression models based on the analysis of static tumor samples. We demonstrate the reliability of the method with simulated data, and describe the application to breast cancer data. Our findings support a linear, branching model ...[more]
Similar Datasets
Project description:As with any biological process, cancer development is inherently dynamic. While major efforts continue to catalog the genomic events associated with human cancer, it remains difficult to interpret and extrapolate the accumulating data to provide insights into the dynamic aspects of the disease. Here, we present a computational strategy that enables the construction of a cancer progression model using static tumor sample data. The developed approach overcame many technical limitations of existing methods. Application of the approach to breast cancer data revealed a linear, branching model with two distinct trajectories for malignant progression. The validity of the constructed model was demonstrated in 27 independent breast cancer data sets, and through visualization of the data in the context of disease progression we were able to identify a number of potentially key molecular events in the advance of breast cancer to malignancy.
| S-EPMC5436003 | biostudies-literature
Project description:Gastric carcinogenesis is a multifactorial process described as a stepwise progression from non-active gastritis (NAG), chronic active gastritis (CAG), precursor lesions of gastric cancer (PLGC) and gastric adenocarcinoma. Gastric cancer (GC) 5-year survival rate is highly dependent upon stage of disease at diagnosis, which is based on endoscopy, biopsy and pathological examinations. Non-invasive GC biomarkers would facilitate its diagnosis at early stages leading to improved GC prognosis. We analyzed plasma samples collected from 80 patients diagnosed with NAG without H. pylori infection (NAG-), CAG with H. pylori infection (CAG+), PLGC and GC. A panel of 208 metabolites including acylcarnitines, amino acids and biogenic amines, sphingolipids, glycerophospholipids, hexoses, and tryptophan and phenylalanine metabolites were quantified using two complementary quantitative approaches: Biocrates AbsoluteIDQ®p180 kit and a LC-MS method designed for the analysis of 29 tryptophan pathway and phenylalanine metabolites. Significantly altered metabolic profiles were found in GC patients that allowing discrimination from NAG-, CAG+ and PLGC patients. Pathway analysis showed significantly altered tryptophan and nitrogen metabolic pathways (FDR P?<?0.01). Three metabolites (histidine, tryprophan and phenylacetylglutamine) discriminated between non-GC and GC groups. These metabolic signatures open new possibilities to improve surveillance of PLGC patients using a minimally invasive blood analysis.
| S-EPMC5736578 | biostudies-literature
Project description:Kidney diseases can cause severe morbidity, mortality, and health burden. Determining the risk factors associated with kidney damage and deterioration has become a priority for the prevention and treatment of kidney disease. This study followed 1042 chronic kidney disease (CKD) patients with Stage 3-5 kidney disease who were treated at a public veteran's hospital through the national prevention program. A total of 12.5 years of records of clinical measurements were collected and analyzed using dynamic and static Cox hazard models to predict the progression to dialysis treatment. The results showed that the statistical significance of several variables in patients with Stage 3-5 CKD was attenuated while the dynamic model was being used. The estimated glomerular filtration rate (eGFR) and urine protein to creatinine ratio (PCR) had the powerful ability to predict the progression of CKD patients with Stage 3a and Stage 3b-5 kidney disease, whereas serum calcium was also predictive for the progression of Stages 3b-5 CKD. Because these two sub-stages of Stage 3 CKD are often associated with differences in routine measurements and the risk analysis of renal dialysis, future research can use this predictive model as a reference while similar prevention programs are implemented.
| S-EPMC8306754 | biostudies-literature
Project description:Cancer is a genomic disease involving various intertwined pathways with complex cross-communication links. Conceptually, this complex interconnected system forms a network, which allows one to model the dynamic behavior of the elements that characterize it to describe the entire system's development in its various evolutionary stages of carcinogenesis. Knowing the activation or inhibition status of the genes that make up the network during its temporal evolution is necessary for the rational intervention on the critical factors for controlling the system's dynamic evolution. In this report, we proposed a methodology for building data-driven boolean networks that model breast cancer tumors. We defined the network components and topology based on gene expression data from RNA-seq of breast cancer cell lines. We used a Boolean logic formalism to describe the network dynamics. The combination of single-cell RNA-seq and interactome data enabled us to study the dynamics of malignant subnetworks of up-regulated genes. First, we used the same Boolean function construction scheme for each network node, based on canalyzing functions. Using single-cell breast cancer datasets from The Cancer Genome Atlas, we applied a binarization algorithm. The binarized version of scRNA-seq data allowed identifying attractors specific to patients and critical genes related to each breast cancer subtype. The model proposed in this report may serve as a basis for a methodology to detect critical genes involved in malignant attractor stability, whose inhibition could have potential applications in cancer theranostics.
| S-EPMC8564392 | biostudies-literature
Project description:Using laser capture microdissection to isolate over 100 specific cell populations, we report the profiling of prostate cancer progression from benign epithelium to metastatic disease. By analyzing these expression signatures in the context of over 15,000 â??molecular conceptsâ??, or sets of biologically connected genes, we generated an integrative model of prostate cancer progression. Keywords: disease state analysis Experiment Design: â?¢ Type of experiment: This study used microarray experiments to profile prostate cancer progression through the isolation of pure cell populations using laser capture microdissection (LCM) and OmniPlex Whole Transcriptome (WTA) Amplification. â?¢ Experimental factors: This study profiled various pure cell populations isolated by LCM from prostate tissue. â?¢ The number of hybridizations performed in the experiment: A total of 104 hybridizations were performed. â?¢ The type of reference used for the hybridizations: A single commercially available pool of benign prostate tissue (CPP, Clontech) was used as the reference in all hybridizations. For the reference sample, 10 ng of CPP total RNA was converted into a WTA cDNA library in a volume of 25 µl and five 5 µl aliquots were amplified by WTA PCR and purified. Five ng aliquots of the pooled products were subjected to a second WTA PCR amplification in the presence of amino-allyl dUTP and all products were pooled before labeling to ensure equal representation across all of the hybridizations. â?¢ Hybridization design: For all hybridizations, the experimental prostate sample isolated by LCM was labeled with Cy5 and the common reference of benign pooled prostate as described above was labeled with Cy3. â?¢ Quality control steps taken: In order to assess biological and technical reproducibility, we captured identical regions from two specimens (PCA_11 and MET_HR_3) and subjected one sample, PCA_30, to two hybridizations. For PCA_11, adjacent regions of the same section were captured and treated as separate samples. For MET_HR_3, two separate foci of the liver metastasis (from separate frozen blocks) were captured and treated as separate samples. For PCA_30, two aliquots of the primary WTA PCR products were separately subjected to a second round of WTA PCR amplification. Replicates from PCA_30 and MET_HR_3 were also hybridized to arrays from each of the print runs utilized in this study. Samples used, extract preparation and labeling: â?¢ The origin of the biological sample (for instance, name of the organism, the provider of the sample) and its characteristics: All samples used were of human origin. The table in the MIAME Checklist (Supplmentary Methods) summarizes the samples used. Hybridizations are named according to the class of sample profiled and the background color is the same as shown in Figure S3. Two print runs were used during the course of the study and the indicated run for each hybridization is given. For prostate cancer samples, the Gleason patterns present on the sample section used are indicated, as well as the Gleason patterns of the actual cells isolated by LCM. A unique case number for all profiled samples is given. Descriptions of samples are given where applicable, including the location of metastatic foci. â?¢ Manipulation of biological samples and protocols used: Laser Capture Microdissection (LCM) was performed from frozen tissue sections with the SL Microtest device using µCUT software (MMI). Approximately 10,000 cells were captured for each sample. Serial sections were used if cells could not be obtained from a single section. Total RNA was isolated from captured cells with the RNAqueous Micro kit (Ambion) and treated with DNAse I according to the manufacturer's instructions. RNA quantification was perfomed using Ribogreen (Molecular Probes). â?¢ Protocol for preparing the hybridization extract: For each sample, total RNA (~10ng) in a volume of 25 µl was converted into an OmniPlex WTA cDNA library and amplified by WTA PCR using reagents and protocols according to the beta-commercial TransPlex WTA kit (Rubicon Genomics, Ann Arbor, MI). For each sample, a single 10 µl aliquot of the WTA cDNA library was amplified by WTA PCR and products were purified. Four or five 5 ng aliquots of product were subjected to a second WTA PCR amplification in the presence of amino-allyl dUTP for post amplification labeling and products were pooled before proceeding with the hybridization. Yields after all WTA PCR amplifications were between 2 to 5 μg per reaction, and reaction progress was monitored in real time using SYBR green I (Molecular Probes) on the iCycler IQ (BioRad) or ABI 7300 Real Time PCR system (Applied Biosystems). Real-time PCR amplifications were terminated at or before plateau phase as measured by fluorescence incorporation to preserve maximum representation. â?¢ Labeling protocol(s): For all WTA amplified cDNA samples with amino-allyl dUTP incorporated (10-20 ï?g), nuclease free water was added to 500 ï?l, the sample was transferred to a Microcon YM-30 filter, and centrifuged at 13,000 rcf at RT for 12 minutes. The spin column was inverted in a new collection tube and centrifuged at 13,000 rcf for 2 min at RT. The amount of sample was measured and the volume was adjusted to 5 ï?l with nuclease free water. For labeling, 5 ï?l of 1 M sodium bicarbonate (pH 9.0) was added and allowed to incubate at RT for 15 min. Nine ï?l of anhydrous DMSO were added to Mono Reactive Cy3 and Cy5 dye packs (Amersham, Buckinghamshire, England), mixed thoroughly, and 1 ï?l of the proper dye was added to each sample. The labeling mixture was incubated in the dark at RT for 1 hour. The reaction was stopped by the addition of 9 ï?l of 4 M Hydroxylamine for 15 min. For each labeling mixture, RNase free water was added to 100 ï?l, 500 ï?l of PB buffer was added and Cy3 and Cy5 mixtures were added to separate Qiaquick (Qiagen, Valencia, CA) spin columns and centrifuged at 13,000 rcf. Columns were washed twice with 700 ï?l of buffer PE, and centrifuged dry at 13,000 rcf for 2 min. To elute, 60 ï?l of EB buffer was placed on the column, incubated for 5 min and centrifuged at 13,000 rcf for 1 min. The Cy3 labeled and Cy5 labeled cDNA for each hybridization were mixed and 1.5 ï?l were used for analysis on a ND-1000 spectrophotometer. â?¢ External controls (spikes): Not utilized Hybridization procedures and parameters: â?¢ The protocol and conditions used during hybridization, blocking and washing: For hybridization, 40 ï?g of human Cot-1 DNA (Invitrogen, Carlsbad, CA) was added to 120 ï?l of labeled cDNA, the probe was transferred to a Microcon YM-30 filter and centrifuged at 13,000 rcf for 6 min at RT. The spin column was inverted in a new collection tube and centrifuged at 13,000 rcf for 2 min at RT. The eluted probe volume was adjusted to 18.6 ï?l with nuclease free water and the following were added: 4 ï?l of yeast tRNA (10 ï?g/ï?l, Invitrogen, Carlsbad, CA), 4.9 ï?l of 20X SSC and 0.84 ï?l of 10% SDS. The probe was denatured at 100 Co for 3 min and centrifuged at 13,000 rcf for 45 sec. The probe was added directly to the microarray and a cover slip was added. Slides were hybridized overnight in DIE-TECH 1 or 5 slide hybridization chambers in a 65 Co water bath. After hybridization, cover slips were removed by incubation in 2x SSC / 0.05% SDS. Slides were then washed in 2x SSC / 0.05% SDS for 5 min and 0.2x SSC / 0.05% SDS for five minutes. Slides were then washed in 0.2x SSC for 10 sec and centrifuged dry at 500 rpm for 5 min. Measurement data and specifications: â?¢ Type of scanning hardware and software used: Microarrays were scanned using a GenePix 4000B scanner (Axon Instruments, Union City, CA). â?¢ Type of image analysis software used: Images were gridded and spots were quantified using GenePix Pro 4.0 software (Axon Instruments, Union City, CA). â?¢ A description of the measurements produced by the image-analysis software and a description of which measurements were used in the analysis: The measurements made by GenePix Pro 4.0 are described in the sample files. For all hybridizations, the log2 normalized Median of Ratios (as described below) was used. â?¢ Data selection and transformation procedures: Arrays were auto-gridded by GenePix 4.0. Features flagged by GenePix as not found during grid
2010-06-05 | E-GEOD-6099 | biostudies-arrayexpress
Project description:BackgroundIn the process of constructing a dynamic model of a metabolic pathway, a large number of parameters such as kinetic constants and initial metabolite concentrations are required. However, in many cases, experimental determination of these parameters is time-consuming. Therefore, for large-scale modelling, it is essential to develop a method that requires few experimental parameters. The hybrid dynamic/static (HDS) method is a combination of the conventional kinetic representation and metabolic flux analysis (MFA). Since no kinetic information is required in the static module, which consists of MFA, the HDS method may dramatically reduce the number of required parameters. However, no adequate method for developing a hybrid model from experimental data has been proposed.ResultsIn this study, we develop a method for constructing hybrid models based on metabolome data. The method discriminates enzymes into static modules and dynamic modules using metabolite concentration time series data. Enzyme reaction rate time series were estimated from the metabolite concentration time series data and used to distinguish enzymes optimally for the dynamic and static modules. The method was applied to build hybrid models of two microbial central-carbon metabolism systems using simulation results from their dynamic models.ConclusionA protocol to build a hybrid model using metabolome data and a minimal number of kinetic parameters has been developed. The proposed method was successfully applied to the strictly regulated central-carbon metabolism system, demonstrating the practical use of the HDS method, which is designed for computer modelling of metabolic systems.
| S-EPMC1892778 | biostudies-literature
Project description:PurposeTo compare the detection sensitivities of the progression of retinitis pigmentosa (RP) by automated perimetry to obtain the mean deviation (MD) and total point score and by optical coherence tomography (OCT) to determine the residual ellipsoid zone (EZ) length and thickness of retinal layers.MethodsTwenty-two eyes of 22 patients with RP who underwent annual automated perimetry (Humphrey Field Analyzer 10-2) and OCT examinations during the same period more than four times were included. Disease progression was evaluated using linear regression analysis with the least-squares method. The disease progression speed and interinspection fluctuations for the different examinations were compared using standardized values. The progression detection ability factor, defined as the average of the least squares divided by the square of annual change, was used to compare the sensitivities of the examinations for detecting the progression of RP.ResultsEZ length showed a high correlation with MD (R = 0.87; P = 1.12E-07) at baseline. Disease progression was detected more frequently using EZ length (12/22 eyes) than using MD (3/22 eyes; P = 0.004) or central retinal thickness (1/11 eyes; P = 0.012). Linear regression using standardized values showed that the EZ length had the fastest annual change, with the smallest least absolute values. EZ length was more sensitive for detecting RP progression than MD, total point score, visual acuity, or central retinal thickness.ConclusionsEZ measurement was sensitive for detecting RP progression, and the results of this study indicate that EZ length is appropriate for end points in clinical trials.Translational relevanceThe study provides a basis for conducting future clinical trials.
| S-EPMC8322706 | biostudies-literature
Project description:Apoptotic cell death can be initiated through the extrinsic and intrinsic signaling pathways. While cell cycle progression promotes the responsiveness to intrinsic apoptosis induced by genotoxic stress or spindle poisons, this has not yet been studied conclusively for extrinsic apoptosis. Here, we combined fluorescence-based time-lapse monitoring of cell cycle progression and cell death execution by long-term time-lapse microscopy with sampling-based mathematical modeling to study cell cycle dependency of TRAIL-induced extrinsic apoptosis in NCI-H460/geminin cells. In particular, we investigated the interaction of cell death timing and progression of cell cycle states. We not only found that TRAIL prolongs cycle progression, but in reverse also that cell cycle progression affects the kinetics of TRAIL-induced apoptosis: Cells exposed to TRAIL in G1 died significantly faster than cells stimulated in S/G2/M. The connection between cell cycle state and apoptosis progression was captured by developing a mathematical model, for which parameter estimation revealed that apoptosis progression decelerates in the second half of the cell cycle. Similar results were also obtained when studying HCT-116 cells. Our results therefore reject the null hypothesis of independence between cell cycle progression and extrinsic apoptosis and, supported by simulations and experiments of synchronized cell populations, suggest that unwanted escape from TRAIL-induced apoptosis can be reduced by enriching the fraction of cells in G1 phase. Besides novel insight into the interrelation of cell cycle progression and extrinsic apoptosis signaling kinetics, our findings are therefore also relevant for optimizing future TRAIL-based treatment strategies.
| S-EPMC7271993 | biostudies-literature
Project description:Cancer is a heterogeneous disease often requiring a complexity of alterations to drive a normal cell to a malignancy and ultimately to a metastatic state. Certain genetic perturbations have been implicated for initiation and progression. However, to a great extent, underlying mechanisms often remain elusive. These genetic perturbations are most likely reflected by the altered expression of sets of genes or pathways, rather than individual genes, thus creating a need for models of deregulation of pathways to help provide an understanding of the mechanisms of tumorigenesis. We introduce an integrative hierarchical analysis of tumor progression that discovers which a priori defined pathways are relevant either throughout or in particular steps of progression. Pathway interaction networks are inferred for these relevant pathways over the steps in progression. This is followed by the refinement of the relevant pathways to those genes most differentially expressed in particular disease stages. The final analysis infers a gene interaction network for these refined pathways. We apply this approach to model progression in prostate cancer and melanoma, resulting in a deeper understanding of the mechanisms of tumorigenesis. Our analysis supports previous findings for the deregulation of several pathways involved in cell cycle control and proliferation in both cancer types. A novel finding of our analysis is a connection between ErbB4 and primary prostate cancer.
| S-EPMC2242820 | biostudies-literature
Project description:Citizen-science databases have been used to develop species distribution models (SDMs), although many taxa may be only georeferenced to county. It is tacitly assumed that SDMs built from county-scale data should be less precise than those built with more accurate localities, but the extent of the bias is currently unknown. Our aims in this study were to illustrate the effects of using county-scale data on the spatial extent and accuracy of SDMs relative to true locality data and to compare potential compensatory methods (including increased sample size and using overall county environmental averages rather than point locality environmental data). To do so, we developed SDMs in maxent with PRISM-derived BIOCLIM parameters for 283 and 230 species of odonates (dragonflies and damselflies) and butterflies, respectively, for five subsets from the OdonataCentral and Butterflies and Moths of North America citizen-science databases: (1) a true locality dataset, (2) a corresponding sister dataset of county-centroid coordinates, (3) a dataset where the average environmental conditions within each county were assigned to each record, (4) a 50/50% mix of true localities and county-centroid coordinates, and (5) a 50/50% mix of true localities and records assigned the average environmental conditions within each county. These mixtures allowed us to quantify the degree of bias from county-scale data. Models developed with county centroids overpredicted the extent of suitable habitat by 15% on average compared to true locality models, although larger sample sizes (>100 locality records) reduced this disparity. Assigning county-averaged environmental conditions did not offer consistent improvement, however. Because county-level data are of limited value for developing SDMs except for species that are widespread and well collected or that inhabit regions where small, climatically uniform counties predominate, three means of encouraging more accurate georeferencing in citizen-science databases are provided.
| S-EPMC5551104 | biostudies-literature