Krzysztof Fujarewicz

Gene Expression Profile of Papillary Thyroid Cancer: Sources of Variability and Diagnostic Implications

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.

A multi-gene approach to differentiate papillary thyroid carcinoma from benign lesions: gene selection using support vector machines with bootstrapping

Selection of novel molecular markers is an important goal of cancer genomics studies. The aim of our analysis was to apply the multivariate bioinformatical tools to rank the genes – potential markers of papillary thyroid cancer (PTC) according to their diagnostic usefulness. We also assessed the accuracy of benign/malignant classification, based on gene expression profiling, for PTC. We analyzed a 180-array dataset (90 HG-U95A and 90 HG-U133A oligonucleotide arrays), which included a collection of 57 PTCs, 61 benign thyroid tumors, and 62 apparently normal tissues. Gene selection was carried out by the support vector machines method with bootstrapping, which allowed us 1) ranking the genes that were most important for classification quality and appeared most frequently in the classifiers (bootstrap-based feature ranking, BBFR); 2) ranking the samples, and thus detecting cases that were most difficult to classify (bootstrap-based outlier detection). The accuracy of PTC diagnosis was 98.5% for a 20-gene classifier, its 95% confidence interval (CI) was 95.9–100%, with the lower limit of CI exceeding 95% already for five genes. Only 5 of 180 samples (2.8%) were misclassified in more than 10% of bootstrap iterations. We specified 43 genes which are most suitable as molecular markers of PTC, among them some well-known PTC markers (MET, fibronectin 1, dipeptidylpeptidase 4, or adenosine A1 receptor) and potential new ones (UDP-galactose-4-epimerase, cadherin 16, gap junction protein 3, sushi, nidogen, and EGF-like domains 1, inhibitor of DNA binding 3, RUNX1, leiomodin 1, F-box protein 9, and tripartite motif-containing 58). The highest ranking gene, metallophosphoesterase domain-containing protein 2, achieved 96.7% of the maximum BBFR score.

Analysis options for high-throughput sequencing in miRNA expression profiling

BackgroundRecently high-throughput sequencing (HTS) using next generation sequencing techniques became useful in digital gene expression profiling.Our study introduces analysis options for HTS data based on mapping to miRBase or counting and grouping of identical sequence reads. Those approaches allow a hypothesis free detection of miRNA differential expression.MethodsWe compare our results to microarray and qPCR data from one set of RNA samples. We use Illumina platforms for microarray analysis and miRNA sequencing of 20 samples from benign follicular thyroid adenoma and malignant follicular thyroid carcinoma. Furthermore, we use three strategies for HTS data analysis to evaluate miRNA biomarkers for malignant versus benign follicular thyroid tumors.ResultsHigh correlation of qPCR and HTS data was observed for the proposed analysis methods. However, qPCR is limited in the differential detection of miRNA isoforms. Moreover, we illustrate a much broader dynamic range of HTS compared to microarrays for small RNA studies. Finally, our data confirm hsa-miR-197-3p, hsa-miR-221-3p, hsa-miR-222-3p and both hsa-miR-144-3p and hsa-miR-144-5p as potential follicular thyroid cancer biomarkers.ConclusionsCompared to microarrays HTS provides a global profile of miRNA expression with higher specificity and in more detail. Summarizing of HTS reads as isoform groups (analysis pipeline B) or according to functional criteria (seed analysis pipeline C), which better correlates to results of qPCR are promising new options for HTS analysis. Finally, data opens future miRNA research perspectives for HTS and indicates that qPCR might be limited in validating HTS data in detail.

Using SVD and SVM methods for selection, classification, clustering and modeling of DNA microarray data

Abstract DNA microarray technology is the latest and the most advanced tool for parallel measuring of the activity and interactions of thousands of genes. This modern technology promises new insight into mechanisms of living systems, for example only two high-density oligonucleotide microarrays are sufficient to inspect the whole human genome. However, it provides unprecedented amount of data that require application of advanced computational methods. The appropriate choice of data analysis technique depends both on data and on goals of an experiment. In this paper we focus on two promising methods: singular value decomposition and support vector machines. We discuss the possibility of application of these methods for different purposes; particularly for clustering, classification, feature selection and modeling of dynamics of gene expression. We use for testing presented approaches existing data sets, which are widely available via Internet, and one new tumor/normal thyroid microarray data set.

Stable feature selection and classification algorithms for multiclass microarray data

BackgroundRecent studies suggest that gene expression profiles are a promising alternative for clinical cancer classification. One major problem in applying DNA microarrays for classification is the dimension of obtained data sets. In this paper we propose a multiclass gene selection method based on Partial Least Squares (PLS) for selecting genes for classification. The new idea is to solve multiclass selection problem with the PLS method and decomposition to a set of two-class sub-problems: one versus rest (OvR) and one versus one (OvO). We use OvR and OvO two-class decomposition for other recently published gene selection method. Ranked gene lists are highly unstable in the sense that a small change of the data set often leads to big changes in the obtained ordered lists. In this paper, we take a look at the assessment of stability of the proposed methods. We use the linear support vector machines (SVM) technique in different variants: one versus one, one versus rest, multiclass SVM (MSVM) and the linear discriminant analysis (LDA) as a classifier. We use balanced bootstrap to estimate the prediction error and to test the variability of the obtained ordered lists.ResultsThis paper focuses on effective identification of informative genes. As a result, a new strategy to find a small subset of significant genes is designed. Our results on real multiclass cancer data show that our method has a very high accuracy rate for different combinations of classification methods, giving concurrently very stable feature rankings.ConclusionsThis paper shows that the proposed strategies can improve the performance of selected gene sets substantially. OvR and OvO techniques applied to existing gene selection methods improve results as well. The presented method allows to obtain a more reliable classifier with less classifier error. In the same time the method generates more stable ordered feature lists in comparison with existing methods.ReviewersThis article was reviewed by Prof Marek Kimmel, Dr Hans Binder (nominated by Dr Tomasz Lipniacki) and Dr Yuriy Gusev

Adjoint Systems for Models of Cell Signaling Pathways and their Application to Parameter Fitting

The paper concerns the problem of fitting mathematical models of cell signaling pathways. Such models frequently take the form of sets of nonlinear ordinary differential equations. While the model is continuous in time, the performance index used in the fitting procedure, involves measurements taken at discrete time moments. Adjoint sensitivity analysis is a tool, which can be used for finding the gradient of a performance index in the space of parameters of the model. In the paper a structural formulation of adjoint sensitivity analysis called the Generalized Backpropagation Through Time (GBPTT) is used. The method is especially suited for hybrid, continuous-discrete time systems. As an example we use the mathematical model of the NF-kB regulatory module, which plays a major role in the innate immune response in animals.

miRNAs with the Potential to Distinguish Follicular Thyroid Carcinomas from Benign Follicular Thyroid Tumors: Results of a Meta-analysis

The detection of somatic mutations in indeterminate or follicular proliferation fine-needle aspiration cytologies (FNACs) is able to clarify only a subgroup of those FNACs. Therefore, further markers to differentiate this problematic FNAC category by the identification of mutation negative thyroid cancers and benign nodules are urgently needed. Our objective was to evaluate previously published miRNA markers and discover novel ones from all publicly available miRNA expression profiling data sets. By literature review and data repository search we gathered 3 data sets describing human miRNA expression profiles of follicular thyroid cancer (FTC) and follicular adenoma (FA) samples. Literature review summarized 27 previously published miRNAs, which were validated in the 3 available data sets. By means of uniform statistical analysis 6 further miRNAs were identified and tested in an independent, previously published microarray data set. Meta-analysis confirmed 7 out of 27 previously published, and 4 out of 6 de novo identified miRNAs. The low confirmation rate of previously published miRNA markers was induced by low numbers of samples in the analyzed studies and high false discovery rates that were higher than 0.2. Finally, miR-637, miR-181c-3p, miR-206, and miR-7-5p were discovered as de novo potential FTC markers and validated in at least one independent, previously published data set. Two out of these new identified miRNAs (miR-7-5p and miR-206) were validated by qPCR in an independent sample set of 32 FTC and 46 FA samples. Especially miR-7-5p was able to differentiate benign and malignant thyroid tumors in several datasets.

A NOTE ON CLASSIFICATION OF GENE EXPRESSION DATA USING SUPPORT VECTOR MACHINES

Microarrays provide a new technique of measuring gene expression that attracted a lot of research interest in recent years. It has been suggested that gene expression data from microarrays (biochips) can be utilized in many biomedical areas, for example in cancer classification. Whereas several, new and existing, methods of classification has been tested, a selection of proper (optimal) set of genes, which expression serves during classification, is still an open problem. In this paper we propose a heuristic method of choosing suboptimal set of genes by using support vector machines (SVM). Obtained set of genes optimizes leave-one-out cross-validation error. The method is tested on microarray gene expression data of samples of two cancer types: acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL). The results show that quality of classification is much better than for sets obtained using other methods of feature selection. In addition, we demonstrate that maximum separation in a training data set may lead to deterioration of performance in an independent validation data set, a phenomenon akin to overfitting.

Generalized Backpropagation through Time for Continuous Time Neural Networks and Discrete Time Measurements

This paper deals with the problem of identification of continuous time dynamic neural networks when the measurements are given only at discrete time moments, not necessarily uniformly distributed. It is shown that the modified adjoint system, generating the gradient of the performance index, is a continuous-time system with jumps of state variables at moments corresponding to moments of measurements.

The different radiation response and radiation-induced bystander effects in colorectal carcinoma cells differing in p53 status.

Radiation-induced bystander effect, appearing as different biological changes in cells that are not directly exposed to ionizing radiation but are under the influence of molecular signals secreted by irradiated neighbors, have recently attracted considerable interest due to their possible implication for radiotherapy. However, various cells present diverse radiosensitivity and bystander responses that depend, inter alia, on genetic status including TP53, the gene controlling the cell cycle, DNA repair and apoptosis. Here we compared the ionizing radiation and bystander responses of human colorectal carcinoma HCT116 cells with wild type or knockout TP53 using a transwell co-culture system. The viability of exposed to X-rays (0-8 Gy) and bystander cells of both lines showed a roughly comparable decline with increasing dose. The frequency of micronuclei was also comparable at lower doses but at higher increased considerably, especially in bystander TP53-/- cells. Moreover, the TP53-/- cells showed a significantly elevated frequency of apoptosis, while TP53+/+ counterparts expressed high level of senescence. The cross-matched experiments where irradiated cells of one line were co-cultured with non-irradiated cells of opposite line show that both cell lines were also able to induce bystander effects in their counterparts, however different endpoints revealed with different strength. Potential mediators of bystander effects, IL-6 and IL-8, were also generated differently in both lines. The knockout cells secreted IL-6 at lower doses whereas wild type cells only at higher doses. Secretion of IL-8 by TP53-/- control cells was many times lower than that by TP53+/+ but increased significantly after irradiation. Transcription of the NFκBIA was induced in irradiated TP53+/+ mainly, but in bystanders a higher level was observed in TP53-/- cells, suggesting that TP53 is required for induction of NFκB pathway after irradiation but another mechanism of activation must operate in bystander cells.