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.

Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling

Background Diet plays a role on the development of the immune system, and polyunsaturated fatty acids can modulate the expression of a variety of genes. Human milk contains conjugated linoleic acid (CLA), a fatty acid that seems to contribute to immune development. Indeed, recent studies carried out in our group in suckling animals have shown that the immune function is enhanced after feeding them with an 80:20 isomer mix composed of c9,t11 and t10,c12 CLA. However, little work has been done on the effects of CLA on gene expression, and even less regarding immune system development in early life. Results The expression profile of mesenteric lymph nodes from animals supplemented with CLA during gestation and suckling through dam's milk (Group A) or by oral gavage (Group B), supplemented just during suckling (Group C) and control animals (Group D) was determined with the aid of the specific GeneChip® Rat Genome 230 2.0 (Affymettrix). Bioinformatics analyses were performed using the GeneSpring GX software package v10.0.2 and lead to the identification of 89 genes differentially expressed in all three dietary approaches. Generation of a biological association network evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network. Gene underexpression was confirmed by Real-Time RT-PCR. Conclusions Ctgf, Timp1, Gal and Syt1, among others, are genes modulated by CLA supplementation that may have a role on mucosal immune responses in early life.

Thermal evolution of gene expression profiles in Drosophila subobscura

Background: Despite its pervasiveness, the genetic basis of adaptation resulting in variation directly or indirectly related to temperature (climatic) gradients is poorly understood. By using 3-fold replicated laboratory thermal stocks covering much of the physiologically tolerable temperature range for the temperate (i.e., cold tolerant) species Drosophila subobscura we have assessed whole-genome transcriptional responses after three years of thermal adaptation, when the populations had already diverged for inversion frequencies, pre-adult life history components, and morphological traits. Total mRNA from each population was compared to a reference pool mRNA in a standard, highly replicated two-colour competitive hybridization experiment using cDNA microarrays.Results: A total of 306 (6.6%) cDNA clones were identified as 'differentially expressed' (following a false discovery rate correction) after contrasting the two furthest apart thermal selection regimes (i.e., 13°C vs . 22°C), also including four previously reported candidate genes for thermotolerance in Drosophila (Hsp26, Hsp68, Fst, and Treh). On the other hand, correlated patterns of gene expression were similar in cold- and warm-adapted populations. Analysis of functional categories defined by the Gene Ontology project point to an overrepresentation of genes involved in carbohydrate metabolism, nucleic acids metabolism and regulation of transcription among other categories. Although the location of differently expressed genes was approximately at random with respect to chromosomes, a physical mapping of 88 probes to the polytene chromosomes of D. subobscura has shown that a larger than expected number mapped inside inverted chromosomal segments.Conclusion: Our data suggest that a sizeable number of genes appear to be involved in thermal adaptation in Drosophila, with a substantial fraction implicated in metabolism. This apparently illustrates the formidable challenge to understanding the adaptive evolution of complex trait variation. Furthermore, some clustering of genes within inverted chromosomal sections was detected. Disentangling the effects of inversions will be obviously required in any future approach if we want to identify the relevant candidate genes.

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.

Semiquantitative RT-PCR measurement of gene expression in rat tissues including a correction for varying cell size and number

Background: Current methodology of gene expression analysis limits the possibilities of comparison between cells/tissues of organs in which cell size and/or number changes as a consequence of the study (e.g. starvation). A method relating the abundance of specific mRNA copies per cell may allow direct comparison or different organs and/or changing physiological conditions. Methods: With a number of selected genes, we analysed the relationship of the number of bases and the fluorescence recorded at a present level using cDNA standards. A lineal relationship was found between the final number of bases and the length of the transcript. The constants of this equation and those of the relationship between fluorescence and number of bases in cDNA were determined and a general equation linking the length of the transcript and the initial number of copies of mRNA was deduced for a given pre-established fluorescence setting. This allowed the calculation of the concentration of the corresponding mRNAs per g of tissue. The inclusion of tissue RNA and the DNA content per cell, allowed the calculation of the mRNA copies per cell. Results: The application of this procedure to six genes: Arbp, cyclophilin, ChREBP, T4 deiodinase 2, acetyl-CoA carboxylase 1 and IRS-1, in liver and retroperitoneal adipose tissue of food-restricted rats allowed precise measures of their changes irrespective of the shrinking of the tissue, the loss of cells or changes in cell size, factors that deeply complicate the comparison between changing tissue conditions. The percentage results obtained with the present methods were essentially the same obtained with the delta-delta procedure and with individual cDNA standard curve quantitative RT-PCR estimation. Conclusion: The method presented allows the comparison (i.e. as copies of mRNA per cell) between different genes and tissues, establishing the degree of abundance of the different molecular species tested.