Stein Aerts

Gene prioritization through genomic data fusion.

The identification of genes involved in health and disease remains a challenge. We describe a bioinformatics approach, together with a freely accessible, interactive and flexible software termed Endeavour, to prioritize candidate genes underlying biological processes or diseases, based on their similarity to known genes involved in these phenomena. Unlike previous approaches, ours generates distinct prioritizations for multiple heterogeneous data sources, which are then integrated, or fused, into a global ranking using order statistics. In addition, it offers the flexibility of including additional data sources. Validation of our approach revealed it was able to efficiently prioritize 627 genes in disease data sets and 76 genes in biological pathway sets, identify candidates of 16 mono- or polygenic diseases, and discover regulatory genes of myeloid differentiation. Furthermore, the approach identified a novel gene involved in craniofacial development from a 2-Mb chromosomal region, deleted in some patients with DiGeorge-like birth defects. The approach described here offers an alternative integrative method for gene discovery.

Discovery of functional elements in 12 Drosophila genomes using evolutionary signatures

Sequencing of multiple related species followed by comparative genomics analysis constitutes a powerful approach for the systematic understanding of any genome. Here, we use the genomes of 12 Drosophila species for the de novo discovery of functional elements in the fly. Each type of functional element shows characteristic patterns of change, or ‘evolutionary signatures’, dictated by its precise selective constraints. Such signatures enable recognition of new protein-coding genes and exons, spurious and incorrect gene annotations, and numerous unusual gene structures, including abundant stop-codon readthrough. Similarly, we predict non-protein-coding RNA genes and structures, and new microRNA (miRNA) genes. We provide evidence of miRNA processing and functionality from both hairpin arms and both DNA strands. We identify several classes of pre- and post-transcriptional regulatory motifs, and predict individual motif instances with high confidence. We also study how discovery power scales with the divergence and number of species compared, and we provide general guidelines for comparative studies.

Exome sequencing identifies mutation in CNOT3 and ribosomal genes RPL5 and RPL10 in T-cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is caused by the cooperation of multiple oncogenic lesions. We used exome sequencing on 67 T-ALLs to gain insight into the mutational spectrum in these leukemias. We detected protein-altering mutations in 508 genes, with an average of 8.2 mutations in pediatric and 21.0 mutations in adult T-ALL. Using stringent filtering, we predict seven new oncogenic driver genes in T-ALL. We identify CNOT3 as a tumor suppressor mutated in 7 of 89 (7.9%) adult T-ALLs, and its knockdown causes tumors in a sensitized Drosophila melanogaster model. In addition, we identify mutations affecting the ribosomal proteins RPL5 and RPL10 in 12 of 122 (9.8%) pediatric T-ALLs, with recurrent alterations of Arg98 in RPL10. Yeast and lymphoid cells expressing the RPL10 Arg98Ser mutant showed a ribosome biogenesis defect. Our data provide insights into the mutational landscape of pediatric versus adult T-ALL and identify the ribosome as a potential oncogenic factor.

Endeavour update: a web resource for gene prioritization in multiple species

Endeavour (http://www.esat.kuleuven.be/endeavourweb; this web site is free and open to all users and there is no login requirement) is a web resource for the prioritization of candidate genes. Using a training set of genes known to be involved in a biological process of interest, our approach consists of (i) inferring several models (based on various genomic data sources), (ii) applying each model to the candidate genes to rank those candidates against the profile of the known genes and (iii) merging the several rankings into a global ranking of the candidate genes. In the present article, we describe the latest developments of Endeavour. First, we provide a web-based user interface, besides our Java client, to make Endeavour more universally accessible. Second, we support multiple species: in addition to Homo sapiens, we now provide gene prioritization for three major model organisms: Mus musculus, Rattus norvegicus and Caenorhabditis elegans. Third, Endeavour makes use of additional data sources and is now including numerous databases: ontologies and annotations, protein–protein interactions, cis-regulatory information, gene expression data sets, sequence information and text-mining data. We tested the novel version of Endeavour on 32 recent disease gene associations from the literature. Additionally, we describe a number of recent independent studies that made use of Endeavour to prioritize candidate genes for obesity and Type II diabetes, cleft lip and cleft palate, and pulmonary fibrosis.

ORegAnno: an open-access community-driven resource for regulatory annotation

ORegAnno is an open-source, open-access database and literature curation system for community-based annotation of experimentally identified DNA regulatory regions, transcription factor binding sites and regulatory variants. The current release comprises 30 145 records curated from 922 publications and describing regulatory sequences for over 3853 genes and 465 transcription factors from 19 species. A new feature called the ‘publication queue’ allows users to input relevant papers from scientific literature as targets for annotation. The queue contains 4438 gene regulation papers entered by experts and another 54 351 identified by text-mining methods. Users can enter or ‘check out’ papers from the queue for manual curation using a series of user-friendly annotation pages. A typical record entry consists of species, sequence type, sequence, target gene, binding factor, experimental outcome and one or more lines of experimental evidence. An evidence ontology was developed to describe and categorize these experiments. Records are cross-referenced to Ensembl or Entrez gene identifiers, PubMed and dbSNP and can be visualized in the Ensembl or UCSC genome browsers. All data are freely available through search pages, XML data dumps or web services at: http://www.oreganno.org.

Alteration of the microRNA network during the progression of Alzheimer's disease

An overview of miRNAs altered in Alzheimers disease (AD) was established by profiling the hippocampus of a cohort of 41 late‐onset AD (LOAD) patients and 23 controls, showing deregulation of 35 miRNAs. Profiling of miRNAs in the prefrontal cortex of a second independent cohort of 49 patients grouped by Braak stages revealed 41 deregulated miRNAs. We focused on miR‐132‐3p which is strongly altered in both brain areas. Downregulation of this miRNA occurs already at Braak stages III and IV, before loss of neuron‐specific miRNAs. Next‐generation sequencing confirmed a strong decrease of miR‐132‐3p and of three family‐related miRNAs encoded by the same miRNA cluster on chromosome 17. Deregulation of miR‐132‐3p in AD brain appears to occur mainly in neurons displaying Tau hyper‐phosphorylation. We provide evidence that miR‐132‐3p may contribute to disease progression through aberrant regulation of mRNA targets in the Tau network. The transcription factor (TF) FOXO1a appears to be a key target of miR‐132‐3p in this pathway.

iRegulon: from a gene list to a gene regulatory network using large motif and track collections.

Identifying master regulators of biological processes and mapping their downstream gene networks are key challenges in systems biology. We developed a computational method, called iRegulon, to reverse-engineer the transcriptional regulatory network underlying a co-expressed gene set using cis-regulatory sequence analysis. iRegulon implements a genome-wide ranking-and-recovery approach to detect enriched transcription factor motifs and their optimal sets of direct targets. We increase the accuracy of network inference by using very large motif collections of up to ten thousand position weight matrices collected from various species, and linking these to candidate human TFs via a motif2TF procedure. We validate iRegulon on gene sets derived from ENCODE ChIP-seq data with increasing levels of noise, and we compare iRegulon with existing motif discovery methods. Next, we use iRegulon on more challenging types of gene lists, including microRNA target sets, protein-protein interaction networks, and genetic perturbation data. In particular, we over-activate p53 in breast cancer cells, followed by RNA-seq and ChIP-seq, and could identify an extensive up-regulated network controlled directly by p53. Similarly we map a repressive network with no indication of direct p53 regulation but rather an indirect effect via E2F and NFY. Finally, we generalize our computational framework to include regulatory tracks such as ChIP-seq data and show how motif and track discovery can be combined to map functional regulatory interactions among co-expressed genes. iRegulon is available as a Cytoscape plugin from http://iregulon.aertslab.org.

TOUCAN 2: the all-inclusive open source workbench for regulatory sequence analysis

We present the second and improved release of the TOUCAN workbench for cis-regulatory sequence analysis. TOUCAN implements and integrates fast state-of-the-art methods and strategies in gene regulation bioinformatics, including algorithms for comparative genomics and for the detection of cis-regulatory modules. This second release of TOUCAN has become open source and thereby carries the potential to evolve rapidly. The main goal of TOUCAN is to allow a user to come to testable hypotheses regarding the regulation of a gene or of a set of co-regulated genes. TOUCAN can be launched from this location: .

Melanoma addiction to the long non-coding RNA SAMMSON

Focal amplifications of chromosome 3p13–3p14 occur in about 10% of melanomas and are associated with a poor prognosis. The melanoma-specific oncogene MITF resides at the epicentre of this amplicon. However, whether other loci present in this amplicon also contribute to melanomagenesis is unknown. Here we show that the recently annotated long non-coding RNA (lncRNA) gene SAMMSON is consistently co-gained with MITF. In addition, SAMMSON is a target of the lineage-specific transcription factor SOX10 and its expression is detectable in more than 90% of human melanomas. Whereas exogenous SAMMSON increases the clonogenic potential in trans, SAMMSON knockdown drastically decreases the viability of melanoma cells irrespective of their transcriptional cell state and BRAF, NRAS or TP53 mutational status. Moreover, SAMMSON targeting sensitizes melanoma to MAPK-targeting therapeutics both in vitro and in patient-derived xenograft models. Mechanistically, SAMMSON interacts with p32, a master regulator of mitochondrial homeostasis and metabolism, to increase its mitochondrial targeting and pro-oncogenic function. Our results indicate that silencing of the lineage addiction oncogene SAMMSON disrupts vital mitochondrial functions in a cancer-cell-specific manner; this silencing is therefore expected to deliver highly effective and tissue-restricted anti-melanoma therapeutic responses.

p53 induces formation of NEAT1 lncRNA-containing paraspeckles that modulate replication stress response and chemosensitivity

In a search for mediators of the p53 tumor suppressor pathway, which induces pleiotropic and often antagonistic cellular responses, we identified the long noncoding RNA (lncRNA) NEAT1. NEAT1 is an essential architectural component of paraspeckle nuclear bodies, whose pathophysiological relevance remains unclear. Activation of p53, pharmacologically or by oncogene-induced replication stress, stimulated the formation of paraspeckles in mouse and human cells. Silencing Neat1 expression in mice, which prevents paraspeckle formation, sensitized preneoplastic cells to DNA-damage-induced cell death and impaired skin tumorigenesis. We provide mechanistic evidence that NEAT1 promotes ATR signaling in response to replication stress and is thereby engaged in a negative feedback loop that attenuates oncogene-dependent activation of p53. NEAT1 targeting in established human cancer cell lines induced synthetic lethality with genotoxic chemotherapeutics, including PARP inhibitors, and nongenotoxic activation of p53. This study establishes a key genetic link between NEAT1 paraspeckles, p53 biology and tumorigenesis and identifies NEAT1 as a promising target to enhance sensitivity of cancer cells to both chemotherapy and p53 reactivation therapy.