The study looked for an optimal set of genes differentiating between papillary thyroid cancer (PTC) and normal thyroid tissue and assessed the sources of variability in gene expression profiles. The ...analysis was done by oligonucleotide microarrays (GeneChip HG-U133A) in 50 tissue samples taken intraoperatively from 33 patients (23 PTC patients and 10 patients with other thyroid disease). In the initial group of 16 PTC and 16 normal samples, we assessed the sources of variability in the gene expression profile by singular value decomposition which specified three major patterns of variability. The first and the most distinct mode grouped transcripts differentiating between tumor and normal tissues. Two consecutive modes contained a large proportion of immunity-related genes. To generate a multigene classifier for tumor-normal difference, we used support vector machines-based technique (recursive feature replacement). It included the following 19 genes: DPP4, GJB3, ST14, SERPINA1, LRP4, MET, EVA1, SPUVE, LGALS3, HBB, MKRN2, MRC2, IGSF1, KIAA0830, RXRG, P4HA2, CDH3, IL13RA1, and MTMR4, and correctly discriminated 17 of 18 additional PTC/normal thyroid samples and all 16 samples published in a previous microarray study. Selected novel genes (LRP4, EVA1, TMPRSS4, QPCT, and SLC34A2) were confirmed by Q-PCR. Our results prove that the gene expression signal of PTC is easily detectable even when cancer cells do not prevail over tumor stroma. We indicate and separate the confounding variability related to the immune response. Finally, we propose a potent molecular classifier able to discriminate between PTC and nonmalignant thyroid in more than 90% of investigated samples.
We developed a computational platform including machine learning and a mechanistic mathematical model to find the optimal protocol for administration of platinum-doublet chemotherapy in a palliative ...setting. The platform has been applied to advanced metastatic non-small cell lung cancer (NSCLC). The 42 NSCLC patients treated with palliative intent at Maria Sklodowska-Curie National Research Institute of Oncology, Gliwice Branch, were collected from a retrospective cohort of patients diagnosed in 2004-2014. Patients were followed-up, for three years. Clinical data collected include complete information about the clinical course of the patients including treatment schedule, response according to RECIST classification, and survival. The core of the platform is the mathematical model, in the form of a system of ordinary differential equations, describing dynamics of platinum-sensitive and platinum-resistant cancer cells and interactions reflecting competition for space and resources. The model is simulated stochastically by sampling the parameter values from a joint probability distribution function. The machine learning model is applied to calibrate the mathematical model and to fit it to the overall survival curve. The model simulations faithfully reproduce the clinical cohort at three levels long-term response (OS), the initial response (according to RECIST criteria), and the relationship between the number of chemotherapy cycles and time between two consecutive chemotherapy cycles. In addition, we investigated the relationship between initial and long-term response. We showed that those two variables do not correlate which means that we cannot predict patient survival solely based on the initial response. We also tested several chemotherapy schedules to find the best one for patients treated with palliative intent. We found that the optimal treatment schedule depends, among others, on the strength of competition among various subclones in a tumor. The computational platform developed allows optimizing chemotherapy protocols, within admissible limits of toxicity, for palliative treatment of metastatic NSCLC. The simplicity of the method allows its application to chemotherapy optimization in different cancers.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Drug resistance is a bottleneck in cancer treatment. Commonly, a molecular treatment for cancer leads to the emergence of drug resistance in the long term. Thus, some drugs, despite their initial ...excellent response, are withdrawn from the market. Lung cancer is one of the most mutated cancers, leading to dozens of targeted therapeutics available against it. Here, we have developed a mechanistic mathematical model describing sensitization to nine groups of targeted therapeutics and the emergence of intrinsic drug resistance. As we focus only on intrinsic drug resistance, we perform the computer simulations of the model only until clinical diagnosis. We have utilized, for model calibration, the whole-exome sequencing data combined with clinical information from over 1000 non-small-cell lung cancer patients. Next, the model has been applied to find an answer to the following questions: When does intrinsic drug resistance emerge? And how long does it take for early-stage lung cancer to grow to an advanced stage? The results show that drug resistance is inevitable at diagnosis but not always detectable and that the time interval between early and advanced-stage tumors depends on the selection advantage of cancer cells.
Context: There are an increasing number of studies analyzing gene expression profiles in various benign and malignant thyroid tumors. This creates the opportunity to validate results obtained from ...one microarray study with those from other data sets. This process requires rigorous methods for accurate comparison.
Objective: The ability to compare data sets derived from different Affymetrix GeneChip generations and the influence of intra- and interindividual comparisons of gene expression data were evaluated to build multigene classifiers of benign thyroid nodules to verify a previously proposed papillary thyroid carcinoma (PTC) classifier and to look for molecular pathways essential for PTC oncogenesis.
Methods: Gene expression profile data sets from autonomously functioning and cold thyroid nodules and from PTC were analyzed by support vector machines. GenMAPP analysis was used for PTC data analysis to examine the expression patterns of biologically relevant gene sets.
Results: Only intraindividual reference samples allowed the identification of subtle changes in the expression patterns of relevant signaling cascades, such as the MAPK pathway in PTC. Using an artificial intelligence approach, the autonomously functioning and cold thyroid nodule multigene classifiers were derived and evaluated by cross-comparisons.
Conclusion: We recommend defining classifiers within one generation of gene chips and subsequently checking them across different array generations. Using this approach, we have demonstrated the specificity of a previously reported PTC classifier on an independent collection of benign tumors. Moreover, we propose multigene classifiers for different types of benign thyroid nodules.
Living cells, like whole living organisms during evolution, communicate with their neighbors, interact with the environment, divide, change their phenotypes, and eventually die. The development of ...specific ways of communication (through signaling molecules and receptors) allows some cellular subpopulations to survive better, to coordinate their physiological status, and during embryonal development to create tissues and organs or in some conditions to become tumors. Populations of cells cultured in vitro interact similarly, also competing for space and nutrients and stimulating each other to better survive or to die. The results of these intercellular interactions of different types seem to be good examples of biological evolutionary games, and have been the subjects of simulations by the methods of evolutionary game theory where individual cells are treated as players. Here we present examples of intercellular contacts in a population of living human cancer HeLa cells cultured in vitro and propose an evolutionary game theory approach to model the development of such populations. We propose a new technique termed Mixed Spatial Evolutionary Games (MSEG) which are played on multiple lattices corresponding to the possible cellular phenotypes which gives the possibility of simulating and investigating the effects of heterogeneity at the cellular level in addition to the population level. Analyses performed with MSEG suggested different ways in which cellular populations develop in the case of cells communicating directly and through factors released to the environment.
•Cells in a population communicate by signaling factors and/or by direct contacts.•Communication determines evolution of multi-phenotypic cell populations.•Evolution can be simulated by Spatial Evolutionary Game Theory (SEGT) model.•Cell signaling types are reflected by probabilistic and deterministic SEGT approaches.•Novel Mixed Evolutionary Game Theory simulates well multi-phenotypic evolution.
Purpose
Following the nuclear accidents in Chernobyl and later in Fukushima, the nuclear community has been faced with important issues concerning how to search for and diagnose biological ...consequences of low-dose internal radiation contamination. Although after the Chernobyl accident an increase in childhood papillary thyroid cancer (PTC) was observed, it is still not clear whether the molecular biology of PTCs associated with low-dose radiation exposure differs from that of sporadic PTC.
Methods
We investigated tissue samples from 65 children/young adults with PTC using DNA microarray (Affymetrix, Human Genome U133 2.0 Plus) with the aim of identifying molecular differences between radiation-induced (exposed to Chernobyl radiation, ECR) and sporadic PTC. All participants were resident in the same region so that confounding factors related to genetics or environment were minimized.
Results
There were small but significant differences in the gene expression profiles between ECR and non-ECR PTC (global test,
p
< 0.01), with 300 differently expressed probe sets (
p
< 0.001) corresponding to 239 genes. Multifactorial analysis of variance showed that besides radiation exposure history, the
BRAF
mutation exhibited independent effects on the PTC expression profile; the histological subset and patient age at diagnosis had negligible effects. Ten genes (
PPME1
,
HDAC11
,
SOCS7
,
CIC
,
THRA
,
ERBB2
,
PPP1R9A
,
HDGF
,
RAD51AP1
, and
CDK1
) from the 19 investigated with quantitative RT-PCR were confirmed as being associated with radiation exposure in an independent, validation set of samples.
Conclusion
Significant, but subtle, differences in gene expression in the post-Chernobyl PTC are associated with previous low-dose radiation exposure.
In this paper we consider a problem of controllability of discrete time linear system with randomly jumping parameters which can be described by finite state Markov chain. Necessary and sufficient ...conditions for existence of control which governs the system from any initial conditions to zero, from zero initial condition to any final value and from any initial condition to any final value in given time and with probability one are given. The relations between such kinds of controllability and stochastic stabilizability are discussed.
Evolutionary game theory (EGT) has been widely used to simulate tumour processes. In almost all studies on EGT models analysis is limited to two or three phenotypes. Our model contains four main ...phenotypes. Moreover, in a standard approach only heterogeneity of populations is studied, while cancer cells remain homogeneous. A multilayer approach proposed in this paper enables to study heterogeneity of single cells.
In the extended model presented in this paper we consider four strategies (phenotypes) that can arise by mutations. We propose multilayer spatial evolutionary games (MSEG) played on multiple 2D lattices corresponding to the possible phenotypes. It enables simulation and investigation of heterogeneity on the player-level in addition to the population-level. Moreover, it allows to model interactions between arbitrary many phenotypes resulting from the mixture of basic traits.
Different equilibrium points and scenarios (monomorphic and polymorphic populations) have been achieved depending on model parameters and the type of played game. However, there is a possibility of stable quadromorphic population in MSEG games for the same set of parameters like for the mean-field game.
The model assumes an existence of four possible phenotypes (strategies) in the population of cells that make up tumour. Various parameters and relations between cells lead to complex analysis of this model and give diverse results. One of them is a possibility of stable coexistence of different tumour cells within the population, representing almost arbitrary mixture of the basic phenotypes.
This article was reviewed by Tomasz Lipniacki, Urszula Ledzewicz and Jacek Banasiak.
Our goal was to determine the influence of sex, age and the head/brain size on the compartmental brain volumes in the radiologically verified healthy population (96 subjects; 54 women and 42 men) ...from the Upper Silesia region in Poland. The MRI examinations were done using 3T Philips Achieva with the same T1-weighted and T2-weighted protocols. The image segmentation procedures were performed with SPM (Statistical Parameter Mapping) and FSL-FIRST software. The volumes of 14 subcortical structures for the left and right hemispheres and 4 overall volumes were calculated. The General Linear Models (GLM) analysis was used with and without the Total Brain Volume (TBV) and Intracranial Volume (ICV) parameters as the covariates to study the regional vs. global brain atrophy. After the ICV/TBV adjustments, the majority of sex differences in the specific volumes of interest (VOIs) revealed to be linked to the difference in the head/brain size parameters. The analysis also confirmed the significant effect of the aging process on the brain loss. After the TBV adjustment, the age- and sex-related volumetric trends for the gray and white matter volumes were observed: the negative age dependence of the gray matter volume is more pronounced in the males, while in case of the white matter the positive age-related trend in the female group is weaker. The local losses of the left caudate nucleus and the right thalamus are more advanced than the global brain atrophy. Different head-size correction strategies are not interchangeable and may yield various volumetric results, but when used together, facilitate studies on the regional dependencies inherent to a healthy, but aging, brain.
Examination of physiological processes and the influences of the drugs on them can be efficiently supported by mathematical modeling. One of the biggest problems is related to the exact fitting of ...the parameters of a model. Conditions inside the organism change dynamically, so the rates of processes are very difficult to estimate. Perturbations in the model parameters influence the steady state so a desired therapeutic goal may not be reached. Here we investigate the effect of parameter deviation on the steady state in three simple models of the influence of a therapeutic drug on its target protein. Two types of changes in the model parameters are taken into account: small perturbations in the system parameter values, and changes in the switching time of a specific parameter. Additionally, we examine the systems response in case of a drug concentration decreasing with time.
The models which we analyze are simplified, because we want to avoid influences of complex dynamics on the results. A system with a negative feedback loop is the most robust and the most rapid, so it requires the largest drug dose but the effects are observed very quickly. On the other hand a system with positive feedback is very sensitive to changes, so small drug doses are sufficient to reach a therapeutic target. In systems without feedback or with positive feedback, perturbations in the model parameters have a bigger influence on the reachability of the therapeutic target than in systems with negative feedback. Drug degradation or inactivation in biological systems enforces multiple drug applications to maintain the level of a drug's target under the desired threshold. The frequency of drug application should be fitted to the system dynamics, because the response velocity is tightly related to the therapeutic effectiveness and the time for achieving the goal.
Systems with different types of regulation vary in their dynamics and characteristic features. Depending on the feedback loop, different types of therapy may be the most appropriate, and deviations in the model parameters have different influences on the reachability of the therapeutic target.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK