Comparative modular analysis of gene expression in vertebrate organs.

ABSTRACT: BACKGROUND: The degree of conservation of gene expression between homologous organs largely remains an open question. Several recent studies reported some evidence in favor of such conservation. Most studies compute organs' similarity across all orthologous genes, whereas the expression level of many genes are not informative about organ specificity. RESULTS: Here, we use a modularization algorithm to overcome this limitation through the identification of inter-species co-modules of organs and genes. We identify such co-modules using mouse and human microarray expression data. They are functionally coherent both in terms of genes and of organs from both organisms. We show that a large proportion of genes belonging to the same co-module are orthologous between mouse and human. Moreover, their zebrafish orthologs also tend to be expressed in the corresponding homologous organs. Notable exceptions to the general pattern of conservation are the testis and the olfactory bulb. Interestingly, some co-modules consist of single organs, while others combine several functionally related organs. For instance, amygdala, cerebral cortex, hypothalamus and spinal cord form a clearly discernible unit of expression, both in mouse and human. CONCLUSIONS: Our study provides a new framework for comparative analysis which will be applicable also to other sets of large-scale phenotypic data collected across different species.

MAR elements regulate the probability of epigenetic switching between active and inactive gene expression.

Gene expression often cycles between active and inactive states in eukaryotes, yielding variable or noisy gene expression in the short-term, while slow epigenetic changes may lead to silencing or variegated expression. Understanding how cells control these effects will be of paramount importance to construct biological systems with predictable behaviours. Here we find that a human matrix attachment region (MAR) genetic element controls the stability and heritability of gene expression in cell populations. Mathematical modeling indicated that the MAR controls the probability of long-term transitions between active and inactive expression, thus reducing silencing effects and increasing the reactivation of silent genes. Single-cell short-terms assays revealed persistent expression and reduced expression noise in MAR-driven genes, while stochastic burst of expression occurred without this genetic element. The MAR thus confers a more deterministic behavior to an otherwise stochastic process, providing a means towards more reliable expression of engineered genetic systems.

Hepatitis B Virus e Antigen Physically Associates With Receptor-Interacting Serine/Threonine Protein Kinase 2 and Regulates IL-6 Gene Expression.

We previously reported that hepatitis B virus (HBV) e antigen (HBeAg) inhibits production of interleukin 6 by suppressing NF-κB activation. NF-κB is known to be activated through receptor-interacting serine/threonine protein kinase 2 (RIPK2), and we examined the mechanisms of interleukin 6 regulation by HBeAg. HBeAg inhibits RIPK2 expression and interacts with RIPK2, which may represent 2 mechanisms through which HBeAg blocks nucleotide-binding oligomerization domain-containing protein 1 ligand-induced NF-κB activation in HepG2 cells. Our findings identified novel molecular mechanisms whereby HBeAg modulates intracellular signaling pathways by targeting RIPK2, supporting the concept that HBeAg could impair both innate and adaptive immune responses to promote chronic HBV infection.

ExpressionView--an interactive viewer for modules identified in gene expression data.

SUMMARY: ExpressionView is an R package that provides an interactive graphical environment to explore transcription modules identified in gene expression data. A sophisticated ordering algorithm is used to present the modules with the expression in a visually appealing layout that provides an intuitive summary of the results. From this overview, the user can select individual modules and access biologically relevant metadata associated with them. AVAILABILITY: http://www.unil.ch/cbg/ExpressionView. Screenshots, tutorials and sample data sets can be found on the ExpressionView web site.

MicroRNAs shape circadian hepatic gene expression on a transcriptome-wide scale.

A considerable proportion of mammalian gene expression undergoes circadian oscillations. Post-transcriptional mechanisms likely make important contributions to mRNA abundance rhythms. We have investigated how microRNAs (miRNAs) contribute to core clock and clock-controlled gene expression using mice in which miRNA biogenesis can be inactivated in the liver. While the hepatic core clock was surprisingly resilient to miRNA loss, whole transcriptome sequencing uncovered widespread effects on clock output gene expression. Cyclic transcription paired with miRNA-mediated regulation was thus identified as a frequent phenomenon that affected up to 30% of the rhythmic transcriptome and served to post-transcriptionally adjust the phases and amplitudes of rhythmic mRNA accumulation. However, only few mRNA rhythms were actually generated by miRNAs. Overall, our study suggests that miRNAs function to adapt clock-driven gene expression to tissue-specific requirements. Finally, we pinpoint several miRNAs predicted to act as modulators of rhythmic transcripts, and identify rhythmic pathways particularly prone to miRNA regulation.DOI: http://dx.doi.org/10.7554/eLife.02510.001.

Gene expression signatures delineate biological and prognostic subgroups in peripheral T-cell lymphoma

Peripheral T-cell lymphoma (PTCL) encompasses a heterogeneous group of neoplasms with generally poor clinical outcome. Currently 50% of PTCL cases are not classifiable: PTCL-not otherwise specified (NOS). Gene-expression profiles on 372 PTCL cases were analyzed and robust molecular classifiers and oncogenic pathways that reflect the pathobiology of tumor cells and their microenvironment were identified for major PTCL-entities, including 114 angioimmunoblastic T-cell lymphoma (AITL), 31 anaplastic lymphoma kinase (ALK)-positive and 48 ALK-negative anaplastic large cell lymphoma, 14 adult T-cell leukemia/lymphoma and 44 extranodal NK/T-cell lymphoma that were further separated into NK-cell and gdT-cell lymphomas. Thirty-seven percent of morphologically diagnosed PTCL-NOS cases were reclassified into other specific subtypes by molecular signatures. Reexamination, immunohistochemistry, and IDH2 mutation analysis in reclassified cases supported the validity of the reclassification. Two major molecular subgroups can be identified in the remaining PTCL-NOS cases characterized by high expression of either GATA3 (33%; 40/121) or TBX21 (49%; 59/121). The GATA3 subgroup was significantly associated with poor overall survival (P = .01). High expression of cytotoxic gene-signature within the TBX21 subgroup also showed poor clinical outcome (P = .05). In AITL, high expression of several signatures associated with the tumor microenvironment was significantly associated with outcome. A combined prognostic score was predictive of survival in an independent cohort (P = .004).

Pattern and clinical significance of cancer-testis gene expression in head and neck squamous cell carcinoma

Les cancer-testis antigènes appartiennent à la famille des antigènes tumoraux spécifiques. Ils ont montré un pouvoir immunogène chez les patients porteurs de différents cancers. En effet, ils stimulent sélectivement les lymphocytes cytotoxiques, et leur expression spécifique dans les tissus tumoraux en fait une cible idéale pour une vaccination antitumorale. Le but de cette étude est d'identifier l'expression de certains de ces antigènes, d'analyser leur valeur pronostique et de déterminer la meilleure cible antigénique pour permettre une immunothérapie spécifique dans les carcinomes épidermoïdes des voies aérodigestives supérieures. Le profil et le taux d'expression de 12 cancer-testis antigènes (MAGE-A1, MAGE-A3, MAGE-A4, MAGEA10, MAGE-C2, NY-ESO-1, LAGE-1, SSX-2, SSX-4, BAGE, GAGE-1/2, GAGE-3/4) et de 3 autres antigènes tumoraux spécifiques (PRAME, HERV-K-MEL, NA-17A) ont été évalués par RT-PCR sur 57 échantillons de cancers ORL primaires. Les paramètres tumoraux et cliniques ont été prospectivement collectés afin de corréler ces données avec le résultat de nos investigations immunobiologiques. Quatre-vingt-huit pour cent des tumeurs expriment au moins 1 antigène. Une co-expression de 3 gènes ou plus est détectée chez 59% des patients. MAGE-A4 (60%), MAGE-A3 (51%), PRAME (49%) et HERV-K-MEL (42%) sont les gènes le plus fréquemment exprimés. Ils sont totalement absents des muqueuses saines avoisinantes. La présence de MAGE-A et NY-ESO-1 à la surface des cellules a été vérifiée par immunohistochimie. Nos analyses statistiques ont permis d'identifier une diminution de la survie liée au cancer chez les patients porteurs d'une tumeur exprimant de multiples cancer-testis antigènes et notamment MAGE-A4 dont l'expression indépendante d'autres éléments cliniques s'associe statistiquement à un taux de survie diminué. Nos résultats ont permis d'identifier un rôle pronostique de l'expression des gènes associés aux tumeurs dont l'expression est apparemment liée à un phénotype de malignité plus élevé. Cette constatation, corroborée par l'identification parallèle d'un infiltrat lymphocytaire spécifique confirme l'utilité potentielle de certains cancer-testis antigènes comme cible pour une immunothérapie ciblée dans les carcinomes des voies aérodigestives supérieures

A gene expression signature that distinguishes desmoid tumours from nodular fasciitis.

Nodular fasciitis (NF) is a rapidly growing cellular mass composed of fibroblasts/myofibroblasts, usually localized in subcutaneous tissues, that typically undergoes fibrosis and almost never recurs. Desmoid tumours (DTs) are rare forms of fibroblastic/myofibroblastic growth that arise in deep soft tissues, display a propensity for local infiltration and recurrence, but fail to metastasize. Given that both entities are primarily fibroblastic/myofibroblastic lesions with overlapping histological features, their gene expression profiles were compared to identify differentially expressed genes that may provide not only potential diagnostic markers, but also clues as to the pathogenesis of each disorder. Differentially expressed transcripts (89 clones displaying increased expression in DTs and 246 clones displaying increased expression in NF) included genes encoding several receptor and non-receptor tyrosine kinases (EPHB3, PTPRF, GNAZ, SYK, LYN, EPHA4, BIRC3), transcription factors (TWIST1, PITX2, EYA2, OAS1, MITF, TCF20), and members of the Wnt signalling pathway (AXIN2, WISP1, SFRP). Remarkably, almost one-quarter of the differentially expressed genes encode proteins associated with inflammation and tissue remodelling, including members of the interferon (IFN), tumour necrosis factor (TNF), and transforming growth factor beta (TGF-beta) signalling pathways as well as metalloproteinases (MMP1, 9, 13, 23), urokinase plasminogen activator (PLAU), and cathepsins. The observations provide the first comparative molecular characterization of desmoid tumours and nodular fasciitis and suggest that selected tyrosine kinases, transcription factors, and members of the Wnt, TGF-beta, IFN, and TNF signalling pathways may be implicated in influencing and distinguishing their fate.

Single cell gene expression profiling reveals qualitatively distinct CD8 T cells elicited by different gene-based vaccines

The CD8 T cell response generatedby gene-based vaccines is importantfor protective immunity againstmany infectious diseases but its complexityis incompletely understood.Here, we report that different vaccinesencoding HIV Env elicit qualitativelydistinct CD8 T cells that wereidentified by patterns of gene expressionin individual cells. Three alternativeprime-boost vector combinationsstimulated antigen-specific CD8 Tcell populations of similar magnitudeand function by intracellular cytokinestaining; however, single cell geneexpression profiling enabled the discriminationof distinct CM and EMCD8 cells elicited by the three vaccines.Two previously unrecognizedCD8 T cell subsets have been definedby their coexpression of Eomes,Cxcr3 and Ccr7; or Klrk1, Klrg1 andCcr5 in CM and EM cells respectively.

eVOC: a controlled vocabulary for unifying gene expression data.

Expression data contribute significantly to the biological value of the sequenced human genome, providing extensive information about gene structure and the pattern of gene expression. ESTs, together with SAGE libraries and microarray experiment information, provide a broad and rich view of the transcriptome. However, it is difficult to perform large-scale expression mining of the data generated by these diverse experimental approaches. Not only is the data stored in disparate locations, but there is frequent ambiguity in the meaning of terms used to describe the source of the material used in the experiment. Untangling semantic differences between the data provided by different resources is therefore largely reliant on the domain knowledge of a human expert. We present here eVOC, a system which associates labelled target cDNAs for microarray experiments, or cDNA libraries and their associated transcripts with controlled terms in a set of hierarchical vocabularies. eVOC consists of four orthogonal controlled vocabularies suitable for describing the domains of human gene expression data including Anatomical System, Cell Type, Pathology and Developmental Stage. We have curated and annotated 7016 cDNA libraries represented in dbEST, as well as 104 SAGE libraries,with expression information,and provide this as an integrated, public resource that allows the linking of transcripts and libraries with expression terms. Both the vocabularies and the vocabulary-annotated libraries can be retrieved from http://www.sanbi.ac.za/evoc/. Several groups are involved in developing this resource with the aim of unifying transcript expression information.

Evolution under monogamy feminizes gene expression in Drosophila melanogaster.

Many genes have evolved sexually dimorphic expression as a consequence of divergent selection on males and females. However, because the sexes share a genome, the extent to which evolution can shape gene expression independently in each sex is controversial. Here, we use experimental evolution to reveal suboptimal sex-specific expression for much of the genome. By enforcing a monogamous mating system in populations of Drosophila melanogaster for over 100 generations, we eliminated major components of selection on males: female choice and male-male competition. If gene expression is subject to sexually antagonistic selection, relaxed selection on males should cause evolution towards female optima. Monogamous males and females show this pattern of feminization in both the whole-body and head transcriptomes. Genes with male-biased expression patterns evolved decreased expression under monogamy, while genes with female-biased expression evolved increased expression, relative to polygamous populations. Our results demonstrate persistent and widespread evolutionary tension between male and female adaptation.

Reactive electrophile species activate defense gene expression in Arabidopsis.

Compounds containing alpha,beta-unsaturated carbonyl groups are increasingly implicated as potent regulators of gene expression; some are powerful cytotoxins known to accumulate at the site of lesion formation in host-pathogen interactions. We used a robust measurement of photosynthetic efficiency to quantify the toxicity of a variety of lipid derivatives in Arabidopsis leaves. Small alpha,beta-unsaturated carbonyl compounds (e.g. acrolein and methyl vinyl ketone) were highly active and proved to be potent stimulators of expression of the pathogenesis-related gene HEL (PR4). These small volatile electrophiles were far more active than larger alkenal homologs like 2(E)-hexenal, and activated HEL expression in a manner independent of salicylate, ethylene, and jasmonate production/perception. Electrophile treatment massively increased the levels of unesterified cyclopentenone jasmonates, which themselves are electrophiles. Patterns of gene expression in response to electrophile treatment and in response to avirulent bacteria were compared, which revealed strikingly similar transcript profiles. The results broaden the range of known biologic effects of reactive electrophile species to include the activation of a pathogenesis-related gene (HEL) and genes involved in metabolism. Electrophiles can act as mediators of both genetic and biochemical effects on core defense signal transduction.

Regional and cellular gene expression changes in human Huntington's disease brain.

Huntington's disease (HD) pathology is well understood at a histological level but a comprehensive molecular analysis of the effect of the disease in the human brain has not previously been available. To elucidate the molecular phenotype of HD on a genome-wide scale, we compared mRNA profiles from 44 human HD brains with those from 36 unaffected controls using microarray analysis. Four brain regions were analyzed: caudate nucleus, cerebellum, prefrontal association cortex [Brodmann's area 9 (BA9)] and motor cortex [Brodmann's area 4 (BA4)]. The greatest number and magnitude of differentially expressed mRNAs were detected in the caudate nucleus, followed by motor cortex, then cerebellum. Thus, the molecular phenotype of HD generally parallels established neuropathology. Surprisingly, no mRNA changes were detected in prefrontal association cortex, thereby revealing subtleties of pathology not previously disclosed by histological methods. To establish that the observed changes were not simply the result of cell loss, we examined mRNA levels in laser-capture microdissected neurons from Grade 1 HD caudate compared to control. These analyses confirmed changes in expression seen in tissue homogenates; we thus conclude that mRNA changes are not attributable to cell loss alone. These data from bona fide HD brains comprise an important reference for hypotheses related to HD and other neurodegenerative diseases.

Analysis of gene expression patterns in animals

During my PhD, my aim was to provide new tools to increase our capacity to analyse gene expression patterns, and to study on a large-scale basis the evolution of gene expression in animals. Gene expression patterns (when and where a gene is expressed) are a key feature in understanding gene function, notably in development. It appears clear now that the evolution of developmental processes and of phenotypes is shaped both by evolution at the coding sequence level, and at the gene expression level.Studying gene expression evolution in animals, with complex expression patterns over tissues and developmental time, is still challenging. No tools are available to routinely compare expression patterns between different species, with precision, and on a large-scale basis. Studies on gene expression evolution are therefore performed only on small genes datasets, or using imprecise descriptions of expression patterns.The aim of my PhD was thus to develop and use novel bioinformatics resources, to study the evolution of gene expression. To this end, I developed the database Bgee (Base for Gene Expression Evolution). The approach of Bgee is to transform heterogeneous expression data (ESTs, microarrays, and in-situ hybridizations) into present/absent calls, and to annotate them to standard representations of anatomy and development of different species (anatomical ontologies). An extensive mapping between anatomies of species is then developed based on hypothesis of homology. These precise annotations to anatomies, and this extensive mapping between species, are the major assets of Bgee, and have required the involvement of many co-workers over the years. My main personal contribution is the development and the management of both the Bgee database and the web-application.Bgee is now on its ninth release, and includes an important gene expression dataset for 5 species (human, mouse, drosophila, zebrafish, Xenopus), with the most data from mouse, human and zebrafish. Using these three species, I have conducted an analysis of gene expression evolution after duplication in vertebrates.Gene duplication is thought to be a major source of novelty in evolution, and to participate to speciation. It has been suggested that the evolution of gene expression patterns might participate in the retention of duplicate genes. I performed a large-scale comparison of expression patterns of hundreds of duplicated genes to their singleton ortholog in an outgroup, including both small and large-scale duplicates, in three vertebrate species (human, mouse and zebrafish), and using highly accurate descriptions of expression patterns. My results showed unexpectedly high rates of de novo acquisition of expression domains after duplication (neofunctionalization), at least as high or higher than rates of partitioning of expression domains (subfunctionalization). I found differences in the evolution of expression of small- and large-scale duplicates, with small-scale duplicates more prone to neofunctionalization. Duplicates with neofunctionalization seemed to evolve under more relaxed selective pressure on the coding sequence. Finally, even with abundant and precise expression data, the majority fate I recovered was neither neo- nor subfunctionalization of expression domains, suggesting a major role for other mechanisms in duplicate gene retention.

Plants respond to pathogen infection by enhancing the antifungal gene expression of root-associated bacteria.

Plant health and fitness widely depend on interactions with soil microorganisms. Some bacteria such as pseudomonads can inhibit pathogens by producing antibiotics, and controlling these bacteria could help improve plant fitness. In the present study, we tested whether plants induce changes in the antifungal activity of root-associated bacteria as a response to root pathogens. We grew barley plants in a split-root system with one side of the root system challenged by the pathogen Pythium ultimum and the other side inoculated with the biocontrol strain Pseudomonas fluorescens CHA0. We used reporter genes to follow the expression of ribosomal RNA indicative of the metabolic state and of the gene phlA, required for production of 2,4-diacetylphloroglucinol, a key component of antifungal activity. Infection increased the expression of the antifungal gene phlA. No contact with the pathogen was required, indicating that barley influenced gene expression by the bacteria in a systemic way. This effect relied on increased exudation of diffusible molecules increasing phlA expression, suggesting that communication with rhizosphere bacteria is part of the pathogen response of plants. Tripartite interactions among plants, pathogens, and bacteria appear as a novel determinant of plant response to root pathogens.