Search

Pseudogenes may regulate expression of their parental genes as well as other protein-coding genes through various mechanisms. One such mechanism is the ability to act as competitive endogenous RNA (ceRNA) and participate in microRNA-mediated cross-regulation. Here, we outline how to predict the targets of pseudogene ceRNAs bioinformatically and how to validate them experimentally.

The identification of direct targets of transcription factors is a key problem in the study of gene regulatory networks. However, the use of high throughput experimental methods, such as ChIP-chip and ChIP-sequencing, is limited by their high cost and strong dependence on cellular type and context. We developed a computational method for the genome-wide identification of functional transcription factor binding sites based on positional weight matrices, comparative genomics, and gene expression profiling. The method was applied to Stat3, a transcription factor playing crucial roles in inflammation, immunity and oncogenesis, and able to induce distinct subsets of target genes in different cell types or conditions. A newly generated positional weight matrix enabled us to assign affinity scores of high specificity, as measured by EMSA competition assays. Phylogenetic conservation with 7 vertebrate species was used to select the binding sites most likely to be functional. Validation was carried out on predicted sites within genes identified as differentially expressed in the presence or absence of Stat3 by microarray analysis. Twelve of the fourteen sites tested were bound by Stat3 in vivo, as assessed by Chromatin Immunoprecipitation, allowing us to identify 9 Stat3 transcriptional targets. Given its high validation rate, and the availability of large transcription factor-dependent gene expression datasets obtained under diverse experimental conditions, our approach appears to be a valid alternative to high-throughput experimental assays for the discovery of novel direct targets of transcription factors.

The BioMart Community Portal (www.biomart.org) is a community-driven effort to provide a unified interface to biomedical databases that are distributed worldwide. The portal provides access to numerous database projects supported by 30 scientific organizations. It includes over 800 different biological datasets spanning genomics, proteomics, model organisms, cancer data, ontology information and more. All resources available through the portal are independently administered and funded by their host organizations. The BioMart data federation technology provides a unified interface to all the available data. The latest version of the portal comes with many new databases that have been created by our ever-growing community. It also comes with better support and extensibility for data analysis and visualization tools. A new addition to our toolbox, the enrichment analysis tool is now accessible through graphical and web service interface. The BioMart community portal averages over one million requests per day. Building on this level of service and the wealth of information that has become available, the BioMart Community Portal has introduced a new, more scalable and cheaper alternative to the large data stores maintained by specialized organizations.

Simple Summary: The metabolic alterations characteristic of cancer cells play a significant role in tumors' natural history and response to therapy. Recent technological advances have allowed the production of unprecedented amounts of data on many types of cancers. We exploited the most comprehensive collection of such data, The Cancer Genome Atlas (TCGA), to systematically investigate the associations between metabolic alterations and other tumor features. We used sets of genes known to be associated with specific metabolic pathways to classify patients into "metabolic subtypes". Then, we systematically looked for associations between the metabolic subtypes and other tumor features, including histological classification, patient survival, and genome alterations. Our results, while correlative in nature, can provide a guide to the formulation of specific mechanistic hypotheses to be tested experimentally so as to improve our understanding of the biology of cancer and our ability to tailor therapeutic interventions to the specific features of each patient. The alterations of metabolic pathways in cancer have been investigated for many years, beginning long before the discovery of the role of oncogenes and tumor suppressors, and the last few years have witnessed renewed interest in this topic. Large-scale molecular and clinical data on tens of thousands of samples allow us to tackle the problem from a general point of view. Here, we show that transcriptomic profiles of tumors can be exploited to define metabolic cancer subtypes, which can be systematically investigated for associations with other molecular and clinical data. We find thousands of significant associations between metabolic subtypes and molecular features such as somatic mutations, structural variants, epigenetic modifications, protein abundance and activation, and with clinical/phenotypic data, including survival probability, tumor grade, and histological types, which we make available to the community in a dedicated web resource. Our work provides a methodological framework and a rich database of statistical associations, which will contribute to the understanding of the role of metabolic alterations in cancer and to the development of precision therapeutic strategies. [ABSTRACT FROM AUTHOR]

Formalin fixed paraffin-embedded (FFPE) tumor specimens are the conventionally archived material in clinical practice, representing an invaluable tissue source for biomarkers development, validation and routine implementation. For many prospective clinical trials, this material has been collected allowing for a prospective-retrospective study design which represents a successful strategy to define clinical utility for candidate markers. Gene expression data can be obtained even from FFPE specimens with the broadly used Affymetrix HG-U133 Plus 2.0 microarray platform. Nevertheless, important major discrepancies remain in expression data obtained from FFPE compared to fresh-frozen samples, prompting the need for appropriate data processing which could help to obtain more consistent results in downstream analyses. In a publicly available dataset of matched frozen and FFPE expression data, the performances of different normalization methods and specifically designed Chip Description Files (CDFs) were compared. The use of an alternative CDFs together with fRMA normalization significantly improved frozen-FFPE sample correlations, frozen-FFPE probeset correlations and agreement of differential analysis between different tumor subtypes. The relevance of our optimized data processing was assessed and validated using two independent datasets. In this study we demonstrated that an appropriate data processing can significantly improve the reliability of gene expression data derived from FFPE tissues using the standard Affymetrix platform. Tools for the implementation of our data processing algorithm are made publicly available at http://www.biocut.unito.it/cdf-ffpe/. [ABSTRACT FROM AUTHOR]

MicroRNAs are single-stranded non-coding RNAs that simultaneously down-modulate the expression of multiple genes post-transcriptionally by binding to the 3′UTRs of target mRNAs. Here we used computational methods to predict microRNAs relevant in breast cancer progression. Specifically, we applied different microRNA target prediction algorithms to various groups of differentially expressed protein-coding genes obtained from four breast cancer datasets. Six potential candidates were identified, among them miR-223, previously described to be highly expressed in the tumor microenvironment and known to be actively transferred into breast cancer cells. To investigate the function of miR-223 in tumorigenesis and to define its molecular mechanism, we overexpressed miR-223 in breast cancer cells in a transient or stable manner. Alternatively we overexpressed miR-223 in mouse embryonic fibroblasts or HEK293 cells and used their conditioned medium to treat tumor cells. With both approaches, we obtained elevated levels of miR-223 in tumor cells and observed decreased migration, increased cell death in anoikis conditions and augmented sensitivity to chemotherapy but no effect on adhesion and proliferation. The analysis of miR-223 predicted targets revealed enrichment in cell death and survival-related genes and in pathways frequently altered in breast cancer. Among these genes, we showed that protein levels for STAT5A, ITGA3 and NRAS were modulated by miR-223. In addition, we proved that STAT5A is a direct miR-223 target and highlighted a possible correlation between miR-223 and STAT5A in migration and chemotherapy response. Our investigation revealed that a computational analysis of cancer gene expression datasets can be a relevant tool to identify microRNAs involved in cancer progression and that miR-223 has a prominent role in breast malignancy that could potentially be exploited therapeutically. [ABSTRACT FROM AUTHOR]

Among thousands of long non-coding RNAs (lncRNAs) only a small subset is functionally characterized and the functional annotation of lncRNAs on the genomic scale remains inadequate. In this study we computationally characterized two functionally different parts of human lncRNAs transcriptome based on their ability to bind the polycomb repressive complex, PRC2. This classification is enabled by the fact that while all lncRNAs constitute a diverse set of sequences, the classes of PRC2-binding and PRC2 non-binding lncRNAs possess characteristic combinations of sequence-structure patterns and, therefore, can be separated within the feature space. Based on the specific combination of features, we built several machine-learning classifiers and identified the SVM-based classifier as the best performing. We further showed that the SVM-based classifier is able to generalize on the independent data sets. We observed that this classifier, trained on the human lncRNAs, can predict up to 59.4% of PRC2-binding lncRNAs in mice. This suggests that, despite the low degree of sequence conservation, many lncRNAs play functionally conserved biological roles. [ABSTRACT FROM AUTHOR]

Data normalization is a crucial preliminary step in analyzing genomic datasets. The goal of normalization is to remove global variation to make readings across different experiments comparable. In addition, most genomic loci have non-uniform sensitivity to any given assay because of variation in local sequence properties. In microarray experiments, this non-uniform sensitivity is due to different DNA hybridization and cross-hybridization efficiencies, known as the probe effect. In this paper we introduce a new scheme, called Group Normalization (GN), to remove both global and local biases in one integrated step, whereby we determine the normalized probe signal by finding a set of reference probes with similar responses. Compared to conventional normalization methods such as Quantile normalization and physically motivated probe effect models, our proposed method is general in the sense that it does not require the assumption that the underlying signal distribution be identical for the treatment and control, and is flexible enough to correct for nonlinear and higher order probe effects. The Group Normalization algorithm is computationally efficient and easy to implement. We also describe a variant of the Group Normalization algorithm, called Cross Normalization, which efficiently amplifies biologically relevant differences between any two genomic datasets. [ABSTRACT FROM AUTHOR]

The genome of the fruitfly Drosophila melanogaster contains a single p53-like protein, phylogenetically related to the ancestor of the mammalian p53 family of tumor suppressors. We reasoned that a comprehensive map of the protein interaction profile of Drosophila p53 (Dmp53) might help identify conserved interactions of the entire p53 family in man. Using a genome-scale in vitro expression cloning approach, we identified 91 previously unreported Dmp53 interactors. considerably expanding the current Drosophila p53 interactome. Looking for evolutionary conservation of these interactions, we tested 41 mammalian orthologs and found that 37 bound to one or more p53-family members when overexpressed in human cells. An RNAi-based functional assay for modulation of the p53 pathway returned five positive hits, validating the biological relevance of these interactions. One p53 interactor is GTPBP4, a nucleolar protein involved in 605 ribosome biogenesis. We demonstrate that GTPBP4 knockdown induces p53 accumulation and activation in the absence of nucleolardisruption. In breasttumors with wild-type p53. increased expression of GTPBP4 correlates with reduced patient survival. emphasizing a potential relevance of this regulatory axis in cancer. [ABSTRACT FROM AUTHOR]

The identification of direct targets of transcription factors is a key problem in the study of gene regulatory networks. However, the use of high throughput experimental methods, such as ChIP-chip and ChIP-sequencing, is limited by their high cost and strong dependence on cellular type and context. We developed a computational method for the genome-wide identification of functional transcription factor binding sites based on positional weight matrices, comparative genomics, and gene expression profiling. The method was applied to Stat3, a transcription factor playing crucial roles in inflammation, immunity and oncogenesis, and able to induce distinct subsets of target genes in different cell types or conditions. A newly generated positional weight matrix enabled us to assign affinity scores of high specificity, as measured by EMSA competition assays. Phylogenetic conservation with 7 vertebrate species was used to select the binding sites most likely to be functional. Validation was carried out on predicted sites within genes identified as differentially expressed in the presence or absence of Stat3 by microarray analysis. Twelve of the fourteen sites tested were bound by Stat3 in vivo, as assessed by Chromatin Immunoprecipitation, allowing us to identify 9 Stat3 transcriptional targets. Given its high validation rate, and the availability of large transcription factor-dependent gene expression datasets obtained under diverse experimental conditions, our approach appears to be a valid alternative to high-throughput experimental assays for the discovery of novel direct targets of transcription factors. [ABSTRACT FROM AUTHOR]