Integrative analysis of RUNX1 downstream pathways and target genes.

BACKGROUND: The RUNX1 transcription factor gene is frequently mutated in sporadic myeloid and lymphoid leukemia through translocation, point mutation or amplification. It is also responsible for a familial platelet disorder with predisposition to acute myeloid leukemia (FPD-AML). The disruption of the largely unknown biological pathways controlled by RUNX1 is likely to be responsible for the development of leukemia. We have used multiple microarray platforms and bioinformatic techniques to help identify these biological pathways to aid in the understanding of why RUNX1 mutations lead to leukemia. RESULTS: Here we report genes regulated either directly or indirectly by RUNX1 based on the study of gene expression profiles generated from 3 different human and mouse platforms. The platforms used were global gene expression profiling of: 1) cell lines with RUNX1 mutations from FPD-AML patients, 2) over-expression of RUNX1 and CBFbeta, and 3) Runx1 knockout mouse embryos using either cDNA or Affymetrix microarrays. We observe that our datasets (lists of differentially expressed genes) significantly correlate with published microarray data from sporadic AML patients with mutations in either RUNX1 or its cofactor, CBFbeta. A number of biological processes were identified among the differentially expressed genes and functional assays suggest that heterozygous RUNX1 point mutations in patients with FPD-AML impair cell proliferation, microtubule dynamics and possibly genetic stability. In addition, analysis of the regulatory regions of the differentially expressed genes has for the first time systematically identified numerous potential novel RUNX1 target genes. CONCLUSION: This work is the first large-scale study attempting to identify the genetic networks regulated by RUNX1, a master regulator in the development of the hematopoietic system and leukemia. The biological pathways and target genes controlled by RUNX1 will have considerable importance in disease progression in both familial and sporadic leukemia as well as therapeutic implications

Role of ATH5 in the specification of retinal ganglion cells and expression and cell compartmentalization of EFEMP1, a protein associated with Malattia Leventinese

ABSTRACT : The development of the retina is a very complex process, occurring through the progressive restriction of cell fates, from pluripotent cell populations to complex tissues and organs. In all vertebrate species analyzed so far, retinal differentiation starts with the generation of retinal ganglion cells (RGC)s. One of the documented key essential events in the specification of RGCs is the expression of ATHS, an atonal homolog encoding a bHLH transcription factor. Despite the putative role of master regulator of RGC differentiation, the mechanism of integrating its functions into a coherent program underlying the production of this subclass of retinal neurons has not yet been elucidated. By using chromatin immunoprecipitation combined with microarray (ChIP-on-chip) we have screened for ATH5 direct targets in the developing chick retina at two consecutive periods: E3.5 (stage HH22) and E6 (stage HH30), covering the stages of progenitor proliferation, neuroepithelium patterning, RGC specification, cell cycle exit and early neuronal differentiation. In parallel, complementary analysis with Affymetrix expression microarrays was conducted. We compared RGCs versus retina to see if the targets correspond to genes preferentially expressed in RGCs. We also precociously overexpressed ATH5 in the retina of individual embryo, and contralateral retina vas used as a control. Our integrated approach allowed us to establish a compendium of ATH5-targets and enabled us to position ATH5 in the transcription network underlying neurogenesis in the retina. Malattia Leventinese (ML) is an autosomal, dominant retinal dystrophy characterized by extracellular, amorphous deposits known as drusen, between the retinal pigment epithelium (RPE) and Bruch's membrane. On the genetic level, it has been associated with a single missense mutation (R345W) in a widely expressed gene with unknown function called EFEMP1. We determined expression patterns of the EFEMP1 gene in normal and ML human retinas. Our data shown that the upregulation of EFEMP1 is not specific to ML eye, except for the region of the ciliary body. We also analyzed the cell compartmentalization of different versions of the protein (both wild type and mutant). Our studies indicate that both abnormal expression of the EFEMP1 gene and mutation and accumulation of EFEMP 1 protein (inside or outside the cells) might contribute to the ML pathology. Résumé : 1er partie : L'ontogenèse de la rétine est un processus complexe au cours duquel des cellules progénitrices sont engagée, par vagues successives, dans des lignées où elles vont d'abord être déterminées puis vont se différencier pour finalement construire un tissu rétinien composé de cinq classes de neurones (les photorécepteurs, les cellules horizontales, bipolaires, amacrines et ganglionnaires) et d'une seule de cellules gliales (les cellules de Muller). Chez tous les vertébrés, la neurogenèse rétinienne est d'abord marquée par la production des cellules ganglionnaires (RGCs). La production de cette classe de neurone est liée à l'expression du gène ATH5 qui est un homologue du gène atonal chez la Drosophile et qui code pour un facteur de transcription de la famille des protéines basic Helix-Loop-Helix (bHLH). Malgré le rôle central que joue ATH5 dans la production des RGCs, le mécanisme qui intègre la fonction de cette protéine dans le programme de détermination neuronale et ceci en relation avec le développement de la rétine n'est pas encore élucidé. Grâce à une technologie qui permet de combiner la sélection de fragments de chromatine liant ATH5 et la recherche de séquences grâce à des puces d'ADN non-codants (ChIP-on-chip), nous avons recherché des cibles potentielles de la protéine ATH5 dans la rétine en développement. Nous avons conduit cette recherche à deux stades de développement de manière à englober la phase de prolifération cellulaire, la détermination des RGCs, la sortie du cycle cellulaire ainsi que les premières étapes de la différentiation de ces neurones. Des expériences complémentaires nous ont permis de définir les patrons d'expression des gènes sélectionnés ainsi que l'activité promotrice des éléments de régulation identifiés lors de notre criblage. Ces approches expérimentales diverses et complémentaires nous ont permis de répertorier des gènes cibles de la protéine ATH5 et d'établir ainsi des liens fonctionnels entre des voies métaboliques dont nous ne soupçonnions pas jusqu'alors qu'elles puissent être associées à la production d'une classe de neurones centraux. 2ème partie : Malattia Leventinese (ML) est une maladie génétique qui engendre une dystrophie de la rétine. Elle se caractérise par l'accumulation de dépôt amorphe entre l'épithélium pigmentaire et la membrane de Bruch et connu sous le nom de drusen. Cette maladie est liée à une simple mutation non-sens (R345W) dans un gène dénommé EFEMP1 qui est exprimé dans de nombreux tissus mais dont la fonction reste mal définie. Une étude détaillée de l'expression de ce gène dans des rétines humaines a révélé une expression à un niveau élevé du gène EFEMP1 dans divers tissus de l'oeil ML mais également dans des yeux contrôles. Alors que l'accumulation d'ARN messager EFEMP1 dans les cellules de l'épithélium pigmentaire n'est pas spécifique à ML, l'expression de ce gène dans le corps cilié n'a été observée que dans l'oeil ML. Nous avons également comparé la sécrétion de la protéine sauvage avec celle porteuse de la mutation. En résumé, notre étude révèle que le niveau élevé d'expression du gène EFEMP1 ainsi que l'accumulation de la protéine dans certains compartiments cellulaires pourraient contribuer au développement de pathologies rétiniennes liées à ML.

Genetics for clinicians: from candidate genes to whole genome scans (technological advances).

Human genetics has progressed at an unprecedented pace during the past 10 years. DNA microarrays currently allow screening of the entire human genome with high level of coverage and we are now entering the era of high-throughput sequencing. These remarkable technical advances are influencing the way medical research is conducted and have boosted our understanding of the structure of the human genome as well as of disease biology. In this context, it is crucial for clinicians to understand the main concepts and limitations of modern genetics. This review will describe key concepts in genetics, including the different types of genetic markers in the human genome, review current methods to detect DNA variation, describe major online public databases in genetics, explain key concepts in statistical genetics and finally present commonly used study designs in clinical and epidemiological research. This review will therefore concentrate on human genetic variation analysis.

Bgee: integrating and comparing heterogeneous transcriptome data among species

Gene expression patterns are a key feature in understanding gene function, notably in development. Comparing gene expression patterns between animals is a major step in the study of gene function as well as of animal evolution. It also provides a link between genes and phenotypes. Thus we have developed Bgee, a database designed to compare expression patterns between animals, by implementing ontologies describing anatomies and developmental stages of species, and then designing homology relationships between anatomies and comparison criteria between developmental stages. To define homology relationships between anatomical features we have developed the software Homolonto, which uses a modified ontology alignment approach to propose homology relationships between ontologies. Bgee then uses these aligned ontologies, onto which heterogeneous expression data types are mapped. These already include microarrays and ESTs.

Bioinformática: base de datos de matrices de expresión génica para su análisis vía web

El IBB ha desarrollado un servidor de aplicaciones: http://revolutionresearch.uab.es para el análisis de microarrays. Estas microarrays las obtienen y las suben a la base de datos local los usuarios de la aplicación. En la presente memoria se detalla el proceso realizado para automatizar la subida de microarrays públicas a la base de datos local. Estas microarrays se obtendrán del NCBI. El proceso de descarga de microarrays se realizará cada dos meses y estará sincronizado con un proceso de descarga de genes marcadores de microarrays del NCBI. En la memoria también se describen las fases realizadas para crear la interfaz web para gestionar estas microarrays públicas y las modificaciones realizadas sobre el aplicativo web para permitir realizar análisis con estas microarrays.

Pseudogenes in the ENCODE regions: consensus annotation, analysis of transcription, and evolution

Arising from either retrotransposition or genomic duplication of functional genes, pseudogenes are “genomic fossils” valuable for exploring the dynamics and evolution of genes and genomes. Pseudogene identification is an important problem in computational genomics, and is also critical for obtaining an accurate picture of a genome’s structure and function. However, no consensus computational scheme for defining and detecting pseudogenes has been developed thus far. As part of the ENCyclopedia Of DNA Elements (ENCODE) project, we have compared several distinct pseudogene annotation strategies and found that different approaches and parameters often resulted in rather distinct sets of pseudogenes. We subsequently developed a consensus approach for annotating pseudogenes (derived from protein coding genes) in the ENCODE regions, resulting in 201 pseudogenes, two-thirds of which originated from retrotransposition. A survey of orthologs for these pseudogenes in 28 vertebrate genomes showed that a significant fraction (∼80%) of the processed pseudogenes are primate-specific sequences, highlighting the increasing retrotransposition activity in primates. Analysis of sequence conservation and variation also demonstrated that most pseudogenes evolve neutrally, and processed pseudogenes appear to have lost their coding potential immediately or soon after their emergence. In order to explore the functional implication of pseudogene prevalence, we have extensively examined the transcriptional activity of the ENCODE pseudogenes. We performed systematic series of pseudogene-specific RACE analyses. These, together with complementary evidence derived from tiling microarrays and high throughput sequencing, demonstrated that at least a fifth of the 201 pseudogenes are transcribed in one or more cell lines or tissues.

The DART classification of unannotated transcription within the ENCODE regions: associating transcription with known and novel loci

For the ∼1% of the human genome in the ENCODE regions, only about half of the transcriptionally active regions (TARs) identified with tiling microarrays correspond to annotated exons. Here we categorize this large amount of “unannotated transcription.” We use a number of disparate features to classify the 6988 novel TARs—array expression profiles across cell lines and conditions, sequence composition, phylogenetic profiles (presence/absence of syntenic conservation across 17 species), and locations relative to genes. In the classification, we first filter out TARs with unusual sequence composition and those likely resulting from cross-hybridization. We then associate some of those remaining with proximal exons having correlated expression profiles. Finally, we cluster unclassified TARs into putative novel loci, based on similar expression and phylogenetic profiles. To encapsulate our classification, we construct a Database of Active Regions and Tools (DART.gersteinlab.org). DART has special facilities for rapidly handling and comparing many sets of TARs and their heterogeneous features, synchronizing across builds, and interfacing with other resources. Overall, we find that ∼14% of the novel TARs can be associated with known genes, while ∼21% can be clustered into ∼200 novel loci. We observe that TARs associated with genes are enriched in the potential to form structural RNAs and many novel TAR clusters are associated with nearby promoters. To benchmark our classification, we design a set of experiments for testing the connectivity of novel TARs. Overall, we find that 18 of the 46 connections tested validate by RT-PCR and four of five sequenced PCR products confirm connectivity unambiguously.

Hypothalamus transcriptome profile suggests an anorexia-cachexia syndrome in the anx/anx mouse model

The anx/anx mouse displays poor appetite and lean appearance and is considered a good model for the study of anorexia nervosa. To identify new genes involved in feeding behavior and body weight regulation we performed an expression profiling in the hypothalamus of the anx/anx mice. Using commercial microarrays we detected 156 differentially expressed genes and validated 92 of those using TaqMan low-density arrays. The expression of a set of 87 candidate genes selected based on literature evidences was also quantified by TaqMan low-density arrays. Our results showed enrichment in deregulated genes involved in cell death, cell morphology and cancer as well as an alteration of several signaling circuits involved in energy balance including neuropeptide Y and melanocortin signaling. The expression profile along with the phenotype led us to conclude that anx/anx mice resemble the anorexia-cachexia syndrome typically observed in cancer, infection with human immunodeficiency virus or chronic diseases, rather than starvation, and that anx/anx mice could be considered a good model for the treatment and investigation of this condition.

The quest for a message: budding yeast, a model organism to study the control of pre-mRNA splicing

Removal of introns during pre-mRNA splicing is a critical process in gene expression, and understanding its control at both single-gene and genomic levels is one of the great challenges in Biology. Splicing takes place in a dynamic, large ribonucleoprotein complex known as the spliceosome. Combining Genetics and Biochemistry, Saccharomyces cerevisiae provides insights into its mechanisms, including its regulation by RNA-protein interactions. Recent genome-wide analyses indicate that regulated splicing is broad and biologically relevant even in organisms with a relatively simple intronic structure, such as yeast. Furthermore, the possibility of coordination in splicing regulation at genomic level is becoming clear in this model organism. This should provide a valuable system to approach the complex problem of the role of regulated splicing in genomic expression.

Identification of molecular apocrine breast tumours by microarray analysis.

Previous microarray studies on breast cancer identified multiple tumour classes, of which the most prominent, named luminal and basal, differ in expression of the oestrogen receptor alpha gene (ER). We report here the identification of a group of breast tumours with increased androgen signalling and a 'molecular apocrine' gene expression profile. Tumour samples from 49 patients with large operable or locally advanced breast cancers were tested on Affymetrix U133A gene expression microarrays. Principal components analysis and hierarchical clustering split the tumours into three groups: basal, luminal and a group we call molecular apocrine. All of the molecular apocrine tumours have strong apocrine features on histological examination (P=0.0002). The molecular apocrine group is androgen receptor (AR) positive and contains all of the ER-negative tumours outside the basal group. Kolmogorov-Smirnov testing indicates that oestrogen signalling is most active in the luminal group, and androgen signalling is most active in the molecular apocrine group. ERBB2 amplification is commoner in the molecular apocrine than the other groups. Genes that best split the three groups were identified by Wilcoxon test. Correlation of the average expression profile of these genes in our data with the expression profile of individual tumours in four published breast cancer studies suggest that molecular apocrine tumours represent 8-14% of tumours in these studies. Our data show that it is possible with microarray data to divide mammary tumour cells into three groups based on steroid receptor activity: luminal (ER+ AR+), basal (ER- AR-) and molecular apocrine (ER- AR+).

Predicting cancer involvement of genes from heterogeneous data

Background: Systematic approaches for identifying proteins involved in different types of cancer are needed. Experimental techniques such as microarrays are being used to characterize cancer, but validating their results can be a laborious task. Computational approaches are used to prioritize between genes putatively involved in cancer, usually based on further analyzing experimental data. Results: We implemented a systematic method using the PIANA software that predicts cancer involvement of genes by integrating heterogeneous datasets. Specifically, we produced lists of genes likely to be involved in cancer by relying on: (i) protein-protein interactions; (ii) differential expression data; and (iii) structural and functional properties of cancer genes. The integrative approach that combines multiple sources of data obtained positive predictive values ranging from 23% (on a list of 811 genes) to 73% (on a list of 22 genes), outperforming the use of any of the data sources alone. We analyze a list of 20 cancer gene predictions, finding that most of them have been recently linked to cancer in literature. Conclusion: Our approach to identifying and prioritizing candidate cancer genes can be used to produce lists of genes likely to be involved in cancer. Our results suggest that differential expression studies yielding high numbers of candidate cancer genes can be filtered using protein interaction networks.

Differential Gene Expression between Fall- and Spring-Run Chinook Salmon Assessed by Long Serial Analysis of Gene Expression

Of all Pacific salmonids, Chinook salmon Oncorhynchus tshawytscha display the greatest variability in return times to freshwater. The molecular mechanisms of these differential return times have not been well described. Current methods, such as long serial analysis of gene expression (LongSAGE) and microarrays, allow gene expression to be analyzed for thousands of genes simultaneously. To investigate whether differential gene expression is observed between fall- and spring-run Chinook salmon from California's Central Valley, LongSAGE libraries were constructed. Three libraries containing between 25,512 and 29,372 sequenced tags (21 base pairs/tag) were generated using messenger RNA from the brains of adult Chinook salmon returning in fall and spring and from one ocean-caught Chinook salmon. Tags were annotated to genes using complementary DNA libraries from Atlantic salmon Salmo salar and rainbow trout O. mykiss. Differentially expressed genes, as estimated by differences in the number of sequence tags, were found in all pairwise comparisons of libraries (freshwater versus saltwater = 40 genes; fall versus spring = 11 genes: and spawning versus nonspawning = 51 genes). The gene for ependymin, an extracellular glycoprotein involved in behavioral plasticity in fish, exhibited the most differential expression among the three groupings. Reverse transcription polymerase chain reaction analysis verified the differential expression of ependymin between the fall- and spring-run samples. These LongSAGE libraries, the first reported for Chinook salmon, provide a window of the transcriptional changes during Chinook salmon return migration to freshwater and spawning and increase the amount of expressed sequence data.

Compositions in life science data

The aim of this talk is to convince the reader that there are a lot of interesting statisticalproblems in presentday life science data analysis which seem ultimately connected withcompositional statistics.Key words: SAGE, cDNA microarrays, (1D-)NMR, virus quasispecies

Mononuclear cell transcriptome response after sustained virgin olive oil consumption in humans: an exploratory nutrigenomics study

Virgin olive oil (VOO) is considered to be one of the main components responsible for the health benefits of the Mediterranean diet, particularly against atherosclerosis where peripheral blood mononuclear cells (PBMNCs) play a crucial role in atherosclerosis development and progression. The objective of this article was to identify the PBMNC genes that respond to VOO consumption in order to ascertain the molecular mechanisms underlying the beneficial action of VOO in the prevention of atherosclerosis. Gene expression profiles of PBMNCs from healthy individuals were examined in pooled RNA samples by microarrays after 3 weeks of moderate and regular consumption of VOO, as the main fat source in a diet controlled for antioxidant content. Gene expression was verified by qPCR. The response to VOO consumption was confirmed for individual samples (n = 10) by qPCR for 10 upregulated genes (ADAM17, ALDH1A1, BIRC1, ERCC5, LIAS, OGT, PPARBP, TNFSF10, USP48, and XRCC5). Their putative role in the molecular mechanisms involved in atherosclerosis development and progression is discussed, focusing on a possible relation with VOO consumption. Our data support the hypothesis that 3 weeks of nutritional intervention with VOO supplementation, at doses common in the Mediterranean diet, can alter the expression of genes related to atherosclerosis development and progression.

VARIABILITY OF PHENOTYPIC, PROTEOMIC AND TRANSCRIPTOMIC EXPRESSION OF STAPHYLOCOCCUS AUREUS SURFACE ADHESINS

Staphylococcus aureus is a highly successful pathogen responsible of a wide variety of diseases, from minor skin infection to life-threatening sepsis or infective endocarditis, as well as food poisoning and toxic shock syndrome. This heterogeneity of infections and the ability of S. aureus to develop antibiotic-resistance to virtually any available drugs reflect its extraordinary capacity to adapt and survive in a great variety of environments. The pathogenesis of S. aureus infection involves a wide range of cell wall-associated adhesins and extracellular toxins that promote host colonization and invasion. In addition, S. aureus is extremely well equipped with regulatory systems that sense environmental conditions and respond by fine tuning the expression of metabolic and virulence determinants. Surface adhesins referred to MSCRAMMs - for Microbial Surface Component Recognizing Adherence Matrix Molecules - mediate binding to the host extracellular matrix or serum components, including fibrinogen, fibronectin, collagen and elastin, and promote tissue colonization and invasion. Major MSCRAMMs include a family of surface-attached proteins covalently bound to the cell wall peptidoglycan via a conserved LPXTG motif. Genomic analyses indicate that S. aureus contain up to 22 LPXTG surface proteins, which could potentially act individually or in synergy to promote infection. In the first part of this study we determined the range of adherence phenotypes to fibrinogen and fibronectin among 30 carriage isolates of S. aureus and compared it to the adherence phenotypes of 30 infective endocarditis and 30 blood culture isolates. Overall there were great variations in in vitro adherence, but no differences were observed between carriage and infection strains. We further determined the relation between in vitro adherence and in vivo infectivity in a rat model of experimental endocarditis, using 4 isolates that displayed either extremely low or high adherence phenotypes. Unexpectedly, no differences were observed between the in vivo infectivity of isolates that were poorly and highly adherent in vitro. We concluded that the natural variability of in vitro adherence to fibrinogen and fibronectin did not correlate with in vivo infectivity, and thus that pathogenic differences between various strains might only be expressed in in vivo conditions, but not in vitro. Therefore, considering the importance of adhesins expression for infection, direct measurement of those adhesins present on the bacterial surface were made by proteomic approach. 5 In the second series of experiments we assessed the physical presence of the LPXTG species at the staphylococcal surface, as measured at various time points during growth in different culture media. S. aureus Newman was grown in either tryptic soy broth (TSB) or in Roswell Park Memorial Institute (RPMI) culture medium, and samples were removed from early exponential growth phase to late stationary phase. Experiments were performed with mutants in the global accessory-gene regulator (agr), surface protein A (Spa) and clumping factor A (ClfA). Peptides of surface proteins were recovered by "trypsin-shaving" of live bacteria, and semi-quantitative proteomic analysis was performed by tandem liquid-chromatography and mass-spectrometry (LC-MS). We also determined in parallel the mRNA expression by microarrays analysis, as well as the phenotypic adherence of the bacteria to fibrinogen in vitro. The surface proteome was highly complex and contained numerous proteins theoretically not belonging to the bacterial envelope, including ribosomal proteins and metabolic enzymes. Sixteen of the 21 known LPXTG species were detected, but were differentially expressed. As expected, 9 known agr-regulated proteins (e.g. including Spa, FnBPA, ClfA, IsdA, IsdB, SasH, SasD, SasG and FmtB) increased up to the late exponential growth phase, and were abrogated in agr-negative mutants. However, only Spa and SasH modified their proteomic and mRNA profiles in parallel in the parent and its agr negative mutant, while all other LPXTG proteins modified their proteomic profiles independently of their mRNA. Moreover, ClfA became highly transcribed and active in in vitro fibrinogen adherence tests during late growth (24h), whereas it remained poorly detected by proteomics. Differential expression was also detected in iron-rich TSB versus iron-poor RPMI. Proteins from the iron-regulated surface determinant (isd) system, including IsdA, IsdB and IsdH were barely expressed in iron-rich TSB, whereas they increased their expression by >10 time in iron-poor RPMI. We conclude that semi-quantitative proteomic analysis of specific protein species is feasible in S. aureus and that proteomic, transcriptomic and adherence phenotypes demonstrated differential profiles in S. aureus. Furthermore, peptide signatures released by trypsin shaving suggested differential protein domain exposures in various environments, which might be relevant for antiadhesins vaccines. A comprehensive understanding of the S. aureus physiology should integrate all these approaches.