As genes with similar expression pattern are very likely having the same biological function, cluster analysis becomes an important tool to understand and predict gene functions from gene expression profi les. In many situations, each gene expression profi le only contains a few data points. Directly applying traditional clustering algorithms to such short gene expression profi les does not yield satisfactory results. Developing clustering algorithms for short gene expression profi les is necessary. In this thesis, two novel methods are developed for clustering short gene expression pro files. The fi rst method, called the network-based clustering method, deals with the defect of short gene expression profi les by generating a gene co-expression network using conditional mutual information (CMI), which measures the non-linear relationship between two genes, as well as considering indirect gene relationships in the presence of other genes. The network-based clustering method consists of two steps. A gene co-expression network is firstly constructed from short gene expression profi les using a path consistency algorithm (PCA) based on the CMI between genes. Then, a gene functional module is identi ed in terms of cluster cohesiveness. The network-based clustering method is evaluated on 10 large scale Arabidopsis thaliana short time-course gene expression profi le datasets in terms of gene ontology (GO) enrichment analysis, and compared with an existing method called Clustering with Over-lapping Neighbourhood Expansion (ClusterONE). Gene functional modules identi ed by the network-based clustering method for 10 datasets returns target GO p-values as low as 10-24, whereas the original ClusterONE yields insigni cant results. In order to more speci cally cluster gene expression profi les, a second clustering method, namely the protein-protein interaction (PPI) integrated clustering method, is developed. It is designed for clustering short gene expression profi les by integrating gene expression profi le patterns and curated PPI data. The method consists of the three following steps: (1) generate a number of prede ned profi le patterns according to the number of data points in the profi les and assign each gene to the prede fined profi le to which its expression profi le is the most similar; (2) integrate curated PPI data to refi ne the initial clustering result from (1); (3) combine the similar clusters from (2) to gradually reduce cluster numbers by a hierarchical clustering method. The PPI-integrated clustering method is evaluated on 10 large scale A. thaliana datasets using GO enrichment analysis, and by comparison with an existing method called Short Time-series Expression Miner (STEM). Target gene functional clusters identi ed by the PPI-integrated clustering method for 10 datasets returns GO p-values as low as 10-62, whereas STEM returns GO p-values as low as 10-38. In addition to the method development, obtained clusters by two proposed methods are further analyzed to identify cross-talk genes under fi ve stress conditions in root and shoot tissues. A list of potential abiotic stress tolerant genes are found.

No description available.

Indiana University-Purdue University Indianapolis (IUPUI) / New neurons are generated throughout adulthood in two regions of the brain, the dentate gyrus of the hippocampus, which is important for memory formation and cognitive functions, and the sub-ventricular zone of the olfactory bulb, which is important for the sense of smell, and are incorporated into hippocampal network circuitry. Disruption of this process has been postulated to contribute to neurodegenerative disorders including Alzheimer’s disease [1]. AD is the most common form of adult-onset dementia and the number of patients with AD escalates dramatically each year. The generation of new neurons in the dentate gyrus declines with age and in AD. Many of the molecular players in AD are also modulators of adult neurogenesis, but the genetic mechanisms influencing adult neurogenesis in AD are unclear. The overall goal of this project is to identify candidate genes and pathways that play a role in neurogenesis in the adult brain and to test the hypotheses that 1) hippocampal neurogenesis-related genes and pathways are significantly perturbed in AD and 2) neurogenesis-related pathways are significantly associated with hippocampal volume and other AD-related biomarker endophenotypes including brain deposition of amyloid-β and tau pathology. First, potential modulators of adult neurogenesis and their roles in neurodegenerative diseases were evaluated. Candidate genes that control the turnover process of neural stem cells/precursors to new functional neurons during adult neurogenesis were manually curated using a pathway-based systems biology approach. Second, a targeted neurogenesis pathway-based gene analysis was performed resulting in the identification of ADORA2A as associated with hippocampal volume and memory performance in mild cognitive impairment and AD. Third, a genome-wide gene-set enrichment analysis was conducted to discover associations between hippocampal volume and AD related endophenotypes and neurogenesis-related pathways. Within the discovered neurogenesis enriched pathways, a gene-based association analysis identified TESC and ACVR1 as significantly associated with hippocampal volume and APOE and PVLR2 as significantly associated with tau and amyloid beta levels in cerebrospinal fluid. This project identifies new genetic contributions to hippocampal neurogenesis with translational implications for novel therapeutic targets related to learning and memory and neuroprotection in AD.

Pituitary tumor-transforming gene 1 (PTTG1) encodes a multifunctional protein that is involved in many cellular processes. However, the potential role of PTTG1 in tumor formation and its prognostic function in human pan-cancer is still unknown. The analysis of gene alteration, PTTG1 expression, prognostic function, and PTTG1-related immune analysis in 33 types of tumors was performed based on various databases such as The Cancer Genome Atlas database, the Genotype-Tissue Expression database, and the Human Protein Atlas database. Additionally, PTTG1-related gene enrichment analysis was performed to investigate the potential relationship and possible molecular mechanisms between PTTG1 and tumors. Overexpression of PTTG1 may lead to tumor formation and poor prognosis in various tumors. Consequently, PTTG1 acts as a potential oncogene in most tumors. Additionally, PTTG1 is related to immune infiltration, immune checkpoints, tumor mutational burden, and microsatellite instability. Thus, PTTG1 could be potential biomarker for both prognosis and outcomes of tumor treatment and it could also be a promising target in tumor therapy.

Regulator of Chromatin Condensation 1 (RCC1) is the only known guanine nucleotide exchange factor that acts on the Ras-like G protein Ran and plays a key role in cell cycle regulation. Although there is growing evidence to support the relationship between RCC1 and cancer, detailed pancancer analyses have not yet been performed. In this genome database study, based on The Cancer Genome Atlas, Genotype-Tissue Expression and Gene Expression Omnibus databases, the potential role of RCC1 in 33 tumors’ entities was explored. The results show that RCC1 is highly expressed in most human malignant neoplasms in contrast to healthy tissues. RCC1 expression is closely related to the prognosis of a broad variety of tumor patients. Enrichment analysis showed that some tumor-related pathways such as “cell cycle” and “RNA transport” were involved in the functional mechanism of RCC1. In particular, the conducted analysis reveals the relation of RCC1 to multiple immune checkpoint genes and suggests that the regulation of RCC1 is closely related to tumor infiltration of cancer-associated fibroblasts and CD8+ T cells. Coherent data demonstrate the association of RCC1 with the tumor mutation burden and microsatellite instability in various tumors. These findings provide new insights into the role of RCC1 in oncogenesis and tumor immunology in various tumors and indicate its potential as marker for therapy prognosis and targeted treatment strategies.

The overexpression of the enzymes involved in the degradation of procollagen lysine is correlated with various tumor entities. Procollagen-lysine, 2-oxoglutarate 5-dioxygenase 3 (PLOD3) expression was found to be correlated to the progression and migration of cancer cells in gastric, lung and prostate cancer. Here, we analyzed the gene expression, protein expression, and the clinical parameters of survival across 33 cancers based on the Clinical Proteomic Tumor Analysis Consortium (CPTAC), function annotation of the mammalian genome 5 (FANTOM5), Gene Expression Omnibus (GEO), Genotype-Tissue Expression (GTEx), Human Protein Atlas (HPA) and The Cancer Genome Atlas (TCGA) databases. Genetic alteration, immune infiltration and relevant cellular pathways were analyzed in detail. PLOD3 expression negatively correlated with survival periods and the infiltration level of CD8+ T cells, but positively correlated to the infiltration of cancer associated fibroblasts in diverse cancers. Immunohistochemistry in colon carcinomas, glioblastomas, and soft tissue sarcomas further confirm PLOD 3 expression in human cancer tissue. Moreover, amplification and mutation accounted for the largest proportion in esophageal adenocarcinoma and uterine corpus endometrial carcinoma, respectively; the copy number alteration of PLOD3 appeared in all cancers from TCGA; and molecular mechanisms further proved the effect of PLOD3 on tumorigenesis. In particular, PLOD3 expression appears to have a tumor immunological effect, and is related to multiple immune cells. Furthermore, it is also associated with tumor mutation burden and microsatellite instability in various tumors. PLOD3 acts as an inducer of various cancers, and it could be a potential biomarker for prognosis and targeted treatment.

Colorectal cancer is the third most common cancer and the leading cause of cancer deaths in Western industrialised countries. Despite recent advances in the screening, diagnosis, and treatment of colorectal cancer, an estimated 608,000 people die every year due to colon cancer. Our current knowledge of colorectal carcinogenesis indicates a multifactorial and multi-step process that involves various genetic alterations and several biological pathways. The identification of molecular markers with early diagnostic and precise clinical outcome in colon cancer is a challenging task because of tumour heterogeneity. This Ph.D.-thesis presents the molecular and cellular mechanisms leading to colorectal cancer. A systematical review of the literature is conducted on Microarray Gene expression profiling, gene ontology enrichment analysis, microRNA and system Biology and various bioinformatics tools. We aimed this study to stratify a colon tumour into molecular distinct subtypes, identification of novel diagnostic targets and prediction of reliable prognostic signatures for clinical practice using microarray expression datasets. We performed an integrated analysis of gene expression data based on genetic, epigenetic and extensive clinical information using unsupervised learning, correlation and functional network analysis. As results, we identified 267-gene and 124-gene signatures that can distinguish normal, primary and metastatic tissues, and also involved in important regulatory functions such as immune-response, lipid metabolism and peroxisome proliferator-activated receptors (PPARs) signalling pathways. For the first time, we also identify miRNAs that can differentiate between primary colon from metastatic and a prognostic signature of grade and stage levels, which can be a major contributor to complex transcriptional phenotypes in a colon tumour.

No description available.

In my thesis, I focused on integrative analysis of high-throughput oncogenomic data. This was done in two parts: In the first part, I describe IntOGen, an integrative data mining tool for the study of cancer. This system collates, annotates, pre-processes and analyzes large-scale data for transcriptomic, copy number aberration and mutational profiling of a large number of tumors in multiple cancer types. All oncogenomic data is annotated with ICD-O terms. We perform analysis at different levels of complexity: at the level of genes, at the level of modules, at the level of studies and finally combination of studies. The results are publicly available in a web service. I also present the Biomart interface of IntOGen for bulk download of data. In the final part, I propose a methodology based on sample-level enrichment analysis to identify patient subgroups from high-throughput profiling of tumors. I also apply this approach to a specific biological problem and characterize properties of worse prognosis tumor in multiple cancer types. This methodology can be used in the translational version of IntOGen.

In the last decade, molecular biology has expanded from a reductionist view to a systems-wide view that tries to unravel the complex interactions of cellular components. Owing to the emergence of high-throughput technology it is now possible to interrogate entire genomes at an unprecedented resolution. The dimension and unstructured nature of these data made it evident that new methodologies and tools are needed to turn data into biological knowledge. To contribute to this challenge we exploited the wealth of publicly available high-throughput genomics data and developed bioinformatics methodologies focused on extracting information at the pathway rather than the single gene level. First, we developed Gene Set Variation Analysis (GSVA), a method that facilitates the organization and condensation of gene expression profiles into gene sets. GSVA enables pathway-centric downstream analyses of microarray and RNA-seq gene expression data. The method estimates sample-wise pathway variation over a population and allows for the integration of heterogeneous biological data sources with pathway-level expression measurements. To illustrate the features of GSVA, we applied it to several use-cases employing different data types and addressing biological questions. GSVA is made available as an R package within the Bioconductor project. Secondly, we developed a pathway-centric genome-based strategy to reposition drugs in type 2 diabetes (T2D). This strategy consists of two steps, first a regulatory network is constructed that is used to identify disease driving modules and then these modules are searched for compounds that might target them. Our strategy is motivated by the observation that disease genes tend to group together in the same neighborhood forming disease modules and that multiple genes might have to be targeted simultaneously to attain an effect on the pathophenotype. To find potential compounds, we used compound exposed genomics data deposited in public databases. We collected about 20,000 samples that have been exposed to about 1,800 compounds. Gene expression can be seen as an intermediate phenotype reflecting underlying dysregulatory pathways in a disease. Hence, genes contained in the disease modules that elicit similar transcriptional responses upon compound exposure are assumed to have a potential therapeutic effect. We applied the strategy to gene expression data of human islets from diabetic and healthy individuals and identified four potential compounds, methimazole, pantoprazole, bitter orange extract and torcetrapib that might have a positive effect on insulin secretion. This is the first time a regulatory network of human islets has been used to reposition compounds for T2D. In conclusion, this thesis contributes with two pathway-centric approaches to important bioinformatic problems, such as the assessment of biological function and in silico drug repositioning. These contributions demonstrate the central role of pathway-based analyses in interpreting high-throughput genomics data. / En l'última dècada, la biologia molecular ha evolucionat des d'una perspectiva reduccionista cap a una perspectiva a nivell de sistemes que intenta desxifrar les complexes interaccions entre els components cel•lulars. Amb l'aparició de les tecnologies d'alt rendiment actualment és possible interrogar genomes sencers amb una resolució sense precedents. La dimensió i la naturalesa desestructurada d'aquestes dades ha posat de manifest la necessitat de desenvolupar noves eines i metodologies per a convertir aquestes dades en coneixement biològic. Per contribuir a aquest repte hem explotat l'abundància de dades genòmiques procedents d'instruments d'alt rendiment i disponibles públicament, i hem desenvolupat mètodes bioinformàtics focalitzats en l'extracció d'informació a nivell de via molecular en comptes de fer-ho al nivell individual de cada gen. En primer lloc, hem desenvolupat GSVA (Gene Set Variation Analysis), un mètode que facilita l'organització i la condensació de perfils d'expressió dels gens en conjunts. GSVA possibilita anàlisis posteriors en termes de vies moleculars amb dades d'expressió gènica provinents de microarrays i RNA-seq. Aquest mètode estima la variació de les vies moleculars a través d'una població de mostres i permet la integració de fonts heterogènies de dades biològiques amb mesures d'expressió a nivell de via molecular. Per il•lustrar les característiques de GSVA, l'hem aplicat a diversos casos usant diferents tipus de dades i adreçant qüestions biològiques. GSVA està disponible com a paquet de programari lliure per R dins el projecte Bioconductor. En segon lloc, hem desenvolupat una estratègia centrada en vies moleculars basada en el genoma per reposicionar fàrmacs per la diabetis tipus 2 (T2D). Aquesta estratègia consisteix en dues fases: primer es construeix una xarxa reguladora que s'utilitza per identificar mòduls de regulació gènica que condueixen a la malaltia; després, a partir d'aquests mòduls es busquen compostos que els podrien afectar. La nostra estratègia ve motivada per l'observació que els gens que provoquen una malaltia tendeixen a agrupar-se, formant mòduls patogènics, i pel fet que podria caldre una actuació simultània sobre múltiples gens per assolir un efecte en el fenotipus de la malaltia. Per trobar compostos potencials, hem usat dades genòmiques exposades a compostos dipositades en bases de dades públiques. Hem recollit unes 20.000 mostres que han estat exposades a uns 1.800 compostos. L'expressió gènica es pot interpretar com un fenotip intermedi que reflecteix les vies moleculars desregulades subjacents a una malaltia. Per tant, considerem que els gens d'un mòdul patològic que responen, a nivell transcripcional, d'una manera similar a l'exposició del medicament tenen potencialment un efecte terapèutic. Hem aplicat aquesta estratègia a dades d'expressió gènica en illots pancreàtics humans corresponents a individus sans i diabètics, i hem identificat quatre compostos potencials (methimazole, pantoprazole, extracte de taronja amarga i torcetrapib) que podrien tenir un efecte positiu sobre la secreció de la insulina. Aquest és el primer cop que una xarxa reguladora d'illots pancreàtics humans s'ha utilitzat per reposicionar compostos per a T2D. En conclusió, aquesta tesi aporta dos enfocaments diferents en termes de vies moleculars a problemes bioinformàtics importants, com ho son el contrast de la funció biològica i el reposicionament de fàrmacs "in silico". Aquestes contribucions demostren el paper central de les anàlisis basades en vies moleculars a l'hora d'interpretar dades genòmiques procedents d'instruments d'alt rendiment.