Functional analysis: evaluation of response intensities--tailoring ANOVA for lists of expression subsets.

Background: Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms. Results: As several methods assume different null hypotheses, we propose to reformulate the main question that biologists seek to answer. To determine which genesets are associated with expression values that differ between two experiments, we focused on three ad hoc criteria: expression levels, the direction of individual gene expression changes (up or down regulation), and correlations between genes. We introduce the FAERI methodology, tailored from a two-way ANOVA to examine these criteria. The significance of the results was evaluated according to the self-contained null hypothesis, using label sampling or by inferring the null distribution from normally distributed random data. Evaluations performed on simulated data revealed that FAERI outperforms currently available methods for each type of set tested. We then applied the FAERI method to analyze three real-world datasets on hypoxia response. FAERI was able to detect more genesets than other methodologies, and the genesets selected were coherent with current knowledge of cellular response to hypoxia. Moreover, the genesets selected by FAERI were confirmed when the analysis was repeated on two additional related datasets. Conclusions: The expression values of genesets are associated with several biological effects. The underlying mathematical structure of the genesets allows for analysis of data from several genes at the same time. Focusing on expression levels, the direction of the expression changes, and correlations, we showed that two-step data reduction allowed us to significantly improve the performance of geneset analysis using a modified two-way ANOVA procedure, and to detect genesets that current methods fail to detect.

Up-regulation of the levels of expression and function of a constitutively active mutant of the hamster alpha1B-adrenoceptor by ligands that act as inverse agonists.

The alpha1-adrenergic agonist phenylephrine stimulated phospholipase D (PLD) activity in Rat 1 fibroblasts transfected to express either the wild-type hamster alpha1B-adrenoceptor or a constitutively active mutant (CAM) form of this receptor. The EC50 for agonist stimulation of PLD activity was substantially lower at the CAM receptor than at the wild-type receptor as previously noted for phenylephrine stimulation of phosphoinositidase C activity. Sustained treatment of cells expressing the CAM alpha1B-adrenoceptor with phentolamine resulted in a marked up-regulation in levels of this receptor with half-maximal effects produced within 24 h and with an EC50 of approx. 40 nM. Such an up-regulation could be produced with a range of other ligands generally viewed as alpha1-adrenoceptor antagonists but equivalent treatment of cells expressing the wild-type alpha1B-adrenoceptor was unable to mimic these effects. After sustained treatment of the CAM alpha1B-adrenoceptor expressing cells with phentolamine, basal PLD activity was increased and phenylephrine was now able to stimulate PLD activity to greater levels than in vehicle-treated CAM alpha1B-adrenoceptor-expressing cells. The EC50 for phenylephrine stimulation of PLD activity was not altered, however, by phentolamine pretreatment and the associated up-regulation of the receptor. After phentolamine-induced up-regulation of basal PLD activity, a range of alpha1-antagonists were shown to possess the characteristics of inverse agonists of the CAM alpha1B-adrenoceptor as they were able to substantially decrease the elevated basal PLD activity.

Post-transcriptional down-regulation of expression of transcription factor NF1 by Ha-ras oncogene.

NF1 is a family of polypeptides that binds to discrete DNA motifs and plays varying roles in the regulation of gene expression. These polypeptides are also thought to mediate the expression of differentiation-specific markers such as adipocyte and mammary cell type-specific genes. The expression of a number of cellular differentiation-specific markers is down-regulated during neoplastic transformation. We therefore investigated whether oncogenic transformation interferes with the action of NF1. Stable transfection of activated Ha-ras into a number of murine cells correlated with a down-regulation of the expression of the NF1 genes NF1/CTF and NF1/X. The down-regulation was not at the transcriptional level but at the level of stability of the NF1 mRNAs. The level of the DNA binding activity of the NF1 proteins was also reduced in Ha-v-ras-transformed cells, and the expression of a gene that depends on this family of transcription factors was specifically repressed. These results demonstrate that an activated Ha-ras-induced pathway destabilizes the half-life of mRNAs encoding specific members in the NF1 family of transcription factors, which leads to a decrease in NF1-dependent gene expression.

Developmental regulation of expression of schizophrenia susceptibility genes in the primate hippocampal formation

The hippocampal formation is essential for normal memory function and is implicated in many neurodevelopmental, neurodegenerative and neuropsychiatric disorders. In particular, abnormalities in hippocampal structure and function have been identified in schizophrenic subjects. Schizophrenia has a strong polygenic component, but the role of numerous susceptibility genes in normal brain development and function has yet to be investigated. Here we described the expression of schizophrenia susceptibility genes in distinct regions of the monkey hippocampal formation during early postnatal development. We found that, as compared with other genes, schizophrenia susceptibility genes exhibit a differential regulation of expression in the dentate gyrus, CA3 and CA1, over the course of postnatal development. A number of these genes involved in synaptic transmission and dendritic morphology exhibit a developmental decrease of expression in CA3. Abnormal CA3 synaptic organization observed in schizophrenics might be related to some specific symptoms, such as loosening of association. Interestingly, changes in gene expression in CA3 might occur at a time possibly corresponding to the late appearance of the first clinical symptoms. We also found earlier changes in expression of schizophrenia susceptibility genes in CA1, which might be linked to prodromal psychotic symptoms. A number of schizophrenia susceptibility genes including APOE, BDNF, MTHFR and SLC6A4 are involved in other disorders, and thus likely contribute to nonspecific changes in hippocampal structure and function that must be combined with the dysregulation of other genes in order to lead to schizophrenia pathogenesis.

Correcting for the bias due to expression specificity improves the estimation of constrained evolution of expression between mouse and human.

MOTIVATION: Comparative analyses of gene expression data from different species have become an important component of the study of molecular evolution. Thus methods are needed to estimate evolutionary distances between expression profiles, as well as a neutral reference to estimate selective pressure. Divergence between expression profiles of homologous genes is often calculated with Pearson's or Euclidean distance. Neutral divergence is usually inferred from randomized data. Despite being widely used, neither of these two steps has been well studied. Here, we analyze these methods formally and on real data, highlight their limitations and propose improvements. RESULTS: It has been demonstrated that Pearson's distance, in contrast to Euclidean distance, leads to underestimation of the expression similarity between homologous genes with a conserved uniform pattern of expression. Here, we first extend this study to genes with conserved, but specific pattern of expression. Surprisingly, we find that both Pearson's and Euclidean distances used as a measure of expression similarity between genes depend on the expression specificity of those genes. We also show that the Euclidean distance depends strongly on data normalization. Next, we show that the randomization procedure that is widely used to estimate the rate of neutral evolution is biased when broadly expressed genes are abundant in the data. To overcome this problem, we propose a novel randomization procedure that is unbiased with respect to expression profiles present in the datasets. Applying our method to the mouse and human gene expression data suggests significant gene expression conservation between these species. CONTACT: marc.robinson-rechavi@unil.ch; sven.bergmann@unil.ch SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Expression analysis of the AtPHO1 gene family in Arabidopsis thaliana and the involvement of AtPHO1 in the regulation of stomatal movements

Résumé Le transfert du phosphate des racines vers les feuilles s'effectue par la voie du xylème. Il a été précédemment démontré que la protéine AtPHO1 était indispensable au transfert du phosphate dans les vaisseaux du xylème des racines chez la plante modèle Arabidopsis thaliana. Le séquençage et l'annotation du génome d'Arabidopsis ont permis d'identifier dix séquences présentant un niveau de similarité significatif avec le gène AtPHO1 et constituant une nouvelle famille de gène appelé la famille de AtPHO1. Basée sur une étude moléculaire et génétique, cette thèse apporte des éléments de réponse pour déterminer le rôle des membres de ia famille de AtPHO1 chez Arabidopsis, inconnue à ce jour. Dans un premier temps, une analyse bioinformatique des séquences protéiques des membres de la famille de AtPHO1 a révélé la présence dans leur région N-terminale d'un domaine nommé SPX. Ce dernier est conservé parmi de nombreuses protéines impliquées dans l'homéostasie du phosphate chez la levure, renforçant ainsi l'hypothèse que les membres de la famille de AtPHO1 auraient comme AtPHO1 un rôle dans l'équilibre du phosphate dans la plante. En parallèle, la localisation tissulaire de l'expression des gènes AtPHO dans Arabidopsis a été identifiée par l'analyse de plantes transgéniques exprimant le gène rapporteur uidA sous le contrôle des promoteurs respectifs des gènes AtPHO. Un profil d'expression de chaque gène AtPHO au cours du développement de la plante a été obtenu. Une expression prédominante au niveau des tissus vasculaires des racines, des feuilles, des tiges et des fleurs a été observée, suggérant que les gènes AtPHO pourraient avoir des fonctions redondantes au niveau du transfert de phosphate dans le cylindre vasculaire de ces différents organes. Toutefois, plusieurs régions promotrices des gènes AtPHO contrôlent également un profil d'expression GUS non-vasculaire, indiquant un rôle putatif des gènes AtPHO dans l'acquisition ou le recyclage de phosphate dans la plante. Dans un deuxième temps, l'analyse de l'expression des gènes AtPHO durant une carence en phosphate a établi que seule l'expression des gènes AtPHO1, AtPHO1; H1 et AtPHO1; H10 est régulée par cette carence. Une étude approfondie de leur expression en réponse à des traitements affectant l'homéostasie du phosphate dans la plante a ensuite démontré leur régulation par différentes voies de signalisation. Ensuite, une analyse détaillée de la régulation de l'expression du gène AtPHO1; H1O dans des feuilles d'Arabidopsis blessées ou déshydratées a révélé que ce gène constitue le premìer gène marqueur d'une nouvelle voie de signalisation induite par l'OPDA, pas par le JA et dépendante de la protéine COI1. Ces résultats démontrent pour la première fois que l'OPDA et le JA peuvent activer différents gènes via des voies de signalisation dépendantes de COI1. Enfin, cette thèse révèle l'identification d'un nouveau rôle de la protéine AtPHO1 dans la régulation de l'action de l'ABA au cours des processus de fermeture stomatique et de germination des graines chez Arabidopsis. Bien que les fonctions exactes des protéines AtPHO restent à être déterminées, ce travail de thèse suggère leur implication dans la propagation de différents signaux dans la plante via la modulation du potentiel membranaire et/ou l'affectation de la composition en ions des cellules comme le font de nombreux transporteurs ou régulateur du transport d'ions. Summary Phosphate is transferred from the roots to the shoot via the xylem. The requirement for AtPHO1 protein to transfer phosphate to the xylem vessels of the root has been previously demonstrated in Arabidopsis thaliana. The sequencing and the annotation of the Arabidopsis genome had allowed the identification of ten sequences that show a significant level of similarity with the AtPHO1 gene. These 10 genes, of unknown functions, constitute a new gene family called the AtPHO1 gene family. Based on a molecular and genetics study, this thesis reveals some information needed to understand the role of the AtPHO1 family members in the plant Arabidopsis. First, a bioinformatics study revealed that the AtPHO sequences contained, in the N-terminal hydrophilic region, a motif called SPX and conserved among multiple proteins involved in phosphate homeostasis in yeast. This finding reinforces the hypothesis that all AtPHO1 family members have, as AtPHO1, a role in phosphate homeostasis. In parallel, we identified the pattern of expression of AtPHO genes in Arabidopsis via analysis of transgenic plants expressing the uidA reporter gene under the control of respective AtPHO promoter regions. The results exhibit a predominant expression of AtPHO genes in vascular tissues of all organs of the plant, implying that these AtPHO genes could have redundant functions in the transfer of phosphate to the vascular cylinder of various organs. The GUS expression pattern for several AtPHO promoter regions was also detected in non-vascular tissue indicating a broad role of AtPHO genes in the acquisition or in the recycling of phosphate in the plant. In a second step, the analysis of the expression of AtPHO genes during phosphate starvation established that only the expression of the AtPHO1, AtPHO1; H1 and AtPHO1; H10 genes were regulated by Pi starvation. Interestingly, different signalling pathways appeared to regulate these three genes during various treatments affecting Pi homeostasis in the plant. The third chapter presents a detailed analysis of the signalling pathways regulating the expression of the AtPHO1; H10 gene in Arabidopsis leaves during wound and dehydrated stresses. Surprisingly, the expression of AtPHO1; H10 was found to be regulated by OPDA (the precursor of JA) but not by JA itself and via the COI1 protein (the central regulator of the JA signalling pathway). These results demonstrated for the first time that OPDA and JA could activate distinct genes via COI1-dependent pathways. Finally, this thesis presents the identification of a novel role of the AtPHO1 protein in the regulation of ABA action in Arabidopsis guard cells and during seed germination. Although the exact role and function of AtPHO1 still need to be determined, these last findings suggest that AtPHO1 and by extension other AtPHO proteins could mediate the propagation of various signals in the plant by modulating the membrane potential and/or by affecting cellular ion composition, as it is the case for many ion transporters or regulators of ion transport.

The SOCS3-independent expression of IDO2 supports the homeostatic generation of T regulatory cells by human dendritic cells.

Dendritic cells (DCs) are professional APCs that have a role in the initiation of adaptive immune responses and tolerance. Among the tolerogenic mechanisms, the expression of the enzyme IDO1 represents an effective tool to generate T regulatory cells. In humans, different DC subsets express IDO1, but less is known about the IDO1-related enzyme IDO2. In this study, we found a different pattern of expression and regulation between IDO1 and IDO2 in human circulating DCs. At the protein level, IDO1 is expressed only in circulating myeloid DCs (mDCs) and is modulated by PGE2, whereas IDO2 is expressed in both mDCs and plasmacytoid DCs and is not modulated by PGE2. In healthy subjects, IDO1 expression requires the presence of PGE2 and needs continuous transcription and translation, whereas IDO2 expression is constitutive, independent from suppressor of cytokine signaling 3 activity. Conversely, in patients suffering from inflammatory arthritis, circulating DCs express both IDO1 and IDO2. At the functional level, both mDCs and plasmacytoid DCs generate T regulatory cells through an IDO1/IDO2-dependent mechanism. We conclude that, in humans, whereas IDO1 provides an additional mechanism of tolerance induced by proinflammatory mediators, IDO2 is stably expressed in steady-state conditions and may contribute to the homeostatic tolerogenic capacity of DCs.

Unexpectedly late expression of intracellular CD3epsilon and TCR gammadelta proteins during adult thymus development.

During adult thymus development immature CD4(-)CD8(-) [double-negative (DN)] precursor cells pass through four phenotypically distinct stages defined by expression of CD44 and CD25: CD44(hi)CD25(-) (DN1), CD44(hi)CD25(+) (DN2), CD44(lo)CD25(+) (DN3) and CD44(lo)CD25(-) (DN4). Although it is well established that the TCR beta, gamma and delta genes are rearranged and expressed in association with the CD3 components in DN thymocytes, the precise timing of expression of the TCR and CD3 proteins has not been determined. In this report we have utilized a sensitive intracellular (ic) staining technique to analyze the expression of ic CD3epsilon, TCR beta and TCR gammadelta proteins in immature DN subsets. As expected from previous studies of TCR beta rearrangement and mRNA expression, icTCR beta(+) cells were first detected in the DN3 subset and their proportion increased thereafter. Surprisingly, however, both icCD3epsilon(+) and icTCR gammadelta(+) cells were detected at later stages of development than was predicted by molecular studies. In particular icCD3epsilon protein expression coincided with the transition from the DN2 to DN3 stage of development, whereas icTCR gammadelta protein expression was only detected in a minor subset of DN4 cells. The implications of these findings for alphabeta lineage divergence will be discussed.

Distinct expression pattern of microtubule-associated protein-2 in human oligodendrogliomas and glial precursor cells.

Microtubule-associated protein 2 (MAP2), a protein linked to the neuronal cytoskeleton in the mature central nervous system (CNS), has recently been identified in glial precursors indicating a potential role during glial development. In the present study, we systematically analyzed the expression of MAP2 in a series of 237 human neuroepithelial tumors including paraffin-embedded specimens and tumor tissue microarrays from oligodendrogliomas, mixed gliomas, astrocytomas, glioblastomas, ependymomas, as well as dysembryoplastic neuroepithelial tumors (DNT), and central neurocytomas. In addition, MAP2-immunoreactive precursor cells were studied in the developing human brain. Three monoclonal antibodies generated against MAP2A-B or MAP2A-D isoforms were used. Variable immunoreactivity for MAP2 could be observed in all gliomas with the exception of ependymomas. Oligodendrogliomas exhibited a consistently strong and distinct pattern of expression characterized by perinuclear cytoplasmic staining without significant process labeling. Tumor cells with immunoreactive bi- or multi-polar processes were mostly encountered in astroglial neoplasms, whereas the small cell component in neurocytomas and DNT was not labeled. These features render MAP2 immunoreactivity a helpful diagnostic tool for the distinction of oligodendrogliomas and other neuroepithelial neoplasms. RT-PCR, Western blot analysis, and in situ hybridization confirmed the expression of MAP2A-C (including the novel MAP2+ 13 transcript) in both oligodendrogliomas and astrocytomas. Double fluorescent laser scanning microscopy showed that GFAP and MAP2 labeled different tumor cell populations. In embryonic human brains, MAP2-immunoreactive glial precursor cells were identified within the subventricular or intermediate zones. These precursors exhibit morphology closely resembling the immunolabeled neoplastic cells observed in glial tumors. Our findings demonstrate MAP2 expression in astrocytic and oligodendroglial neoplasms. The distinct pattern of immunoreactivity in oligodendrogliomas may be useful as a diagnostic tool. Since MAP2 expression occurs transiently in migrating immature glial cells, our findings are in line with an assumed origin of diffuse gliomas from glial precursors.

Mini-review: Specificity and expression of CIITA, the master regulator of MHC class II genes.

The class II transactivator (CIITA) has been referred to as the "master control factor" for the expression of MHC class II (MHCII) genes. As our knowledge on the specificity and function of CIITA grows, it is becoming increasingly evident that this sobriquet is entirely justified. First, despite extensive investigations, the major target genes of CIITA remain those implicated in the presentation of antigenic peptides by MHCII molecules. Although other putative target genes have been reported, the contribution of CIITA to their expression remains indirect, controversial or comparatively minor relative to its decisive role as a regulator of MHCII and related genes. Second, the most important parameter dictating MHCII expression is by far the expression pattern of the gene encoding CIITA (MHC2TA). The vast majority of signals that activate or repress MHCII expression under physiological and pathological situations converge on one or more of the three alternative promoters that drive transcription of the MHC2TA gene. In short, with respect to its specificity and its exquisitely controlled pattern of expression, CIITA is by a long stretch the single most important transcription factor for the regulation of genes required for MHCII-restricted antigen-presentation.

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

Cancer-testis (CT) antigens comprise families of tumor-associated antigens that are immunogenic in patients with various cancers. Their restricted expression makes them attractive targets for immunotherapy. The aim of this study was to determine the expression of several CT genes and evaluate their prognostic value in head and neck squamous cell carcinoma (HNSCC). The pattern and level of expression of 12 CT genes (MAGE-A1, MAGE-A3, MAGE-A4, MAGE-A10, MAGE-C2, NY-ESO-1, LAGE-1, SSX-2, SSX-4, BAGE, GAGE-1/2, GAGE-3/4) and the tumor-associated antigen encoding genes PRAME, HERV-K-MEL, and NA-17A were evaluated by RT-PCR in a panel of 57 primary HNSCC. Over 80% of the tumors expressed at least 1 CT gene. Coexpression of three or more genes was detected in 59% of the patients. MAGE-A4 (60%), MAGE-A3 (51%), PRAME (49%) and HERV-K-MEL (42%) were the most frequently expressed genes. Overall, the pattern of expression of CT genes indicated a coordinate regulation; however there was no correlation between expression of MAGE-A3/A4 and BORIS, a gene whose product has been implicated in CT gene activation. The presence of MAGE-A and NY-ESO-1 proteins was verified by immunohistochemistry. Analysis of the correlation between mRNA expression of CT genes with clinico-pathological characteristics and clinical outcome revealed that patients with tumors positive for MAGE-A4 or multiple CT gene expression had a poorer overall survival. Furthermore, MAGE-A4 mRNA positivity was prognostic of poor outcome independent of clinical parameters. These findings indicate that expression of CT genes is associated with a more malignant phenotype and suggest their usefulness as prognostic markers in HNSCC.

Early expression of the type III secretion system of Parachlamydia acanthamoebae during a replicative cycle within its natural host cell Acanthamoeba castellanii.

The type three secretion system (T3SS) operons of Chlamydiales bacteria are distributed in different clusters along their chromosomes and are conserved at both the level of sequence and genetic organization. A complete characterization of the temporal expression of multiple T3SS components at the transcriptional and protein levels has been performed in Parachlamydia acanthamoebae, replicating in its natural host cell Acanthamoeba castellanii. The T3SS components were classified in four different temporal clusters depending on their pattern of expression during the early, mid- and late phases of the infectious cycle. The putative T3SS transcription units predicted in Parachlamydia are similar to those described in Chlamydia trachomatis, suggesting that T3SS units of transcriptional expression are highly conserved among Chlamydiales bacteria. The maximal expression and activation of the T3SS of Parachlamydia occurred during the early to mid-phase of the infectious cycle corresponding to a critical phase during which the intracellular bacterium has (1) to evade and/or block the lytic pathway of the amoeba, (2) to differentiate from elementary bodies (EBs) to reticulate bodies (RBs), and (3) to modulate the maturation of its vacuole to create a replicative niche able to sustain efficient bacterial growth.

Gene expression variation and expression quantitative trait mapping of human chromosome 21 genes.

Inter-individual differences in gene expression are likely to account for an important fraction of phenotypic differences, including susceptibility to common disorders. Recent studies have shown extensive variation in gene expression levels in humans and other organisms, and that a fraction of this variation is under genetic control. We investigated the patterns of gene expression variation in a 25 Mb region of human chromosome 21, which has been associated with many Down syndrome (DS) phenotypes. Taqman real-time PCR was used to measure expression variation of 41 genes in lymphoblastoid cells of 40 unrelated individuals. For 25 genes found to be differentially expressed, additional analysis was performed in 10 CEPH families to determine heritabilities and map loci harboring regulatory variation. Seventy-six percent of the differentially expressed genes had significant heritabilities, and genomewide linkage analysis led to the identification of significant eQTLs for nine genes. Most eQTLs were in trans, with the best result (P=7.46 x 10(-8)) obtained for TMEM1 on chromosome 12q24.33. A cis-eQTL identified for CCT8 was validated by performing an association study in 60 individuals from the HapMap project. SNP rs965951 located within CCT8 was found to be significantly associated with its expression levels (P=2.5 x 10(-5)) confirming cis-regulatory variation. The results of our study provide a representative view of expression variation of chromosome 21 genes, identify loci involved in their regulation and suggest that genes, for which expression differences are significantly larger than 1.5-fold in control samples, are unlikely to be involved in DS-phenotypes present in all affected individuals.