Meta-analysis of gene expression profiles in breast cancer: toward a unified understanding of breast cancer subtyping and prognosis signatures.

INTRODUCTION: Breast cancer subtyping and prognosis have been studied extensively by gene expression profiling, resulting in disparate signatures with little overlap in their constituent genes. Although a previous study demonstrated a prognostic concordance among gene expression signatures, it was limited to only one dataset and did not fully elucidate how the different genes were related to one another nor did it examine the contribution of well-known biological processes of breast cancer tumorigenesis to their prognostic performance. METHOD: To address the above issues and to further validate these initial findings, we performed the largest meta-analysis of publicly available breast cancer gene expression and clinical data, which are comprised of 2,833 breast tumors. Gene coexpression modules of three key biological processes in breast cancer (namely, proliferation, estrogen receptor [ER], and HER2 signaling) were used to dissect the role of constituent genes of nine prognostic signatures. RESULTS: Using a meta-analytical approach, we consolidated the signatures associated with ER signaling, ERBB2 amplification, and proliferation. Previously published expression-based nomenclature of breast cancer 'intrinsic' subtypes can be mapped to the three modules, namely, the ER-/HER2- (basal-like), the HER2+ (HER2-like), and the low- and high-proliferation ER+/HER2- subtypes (luminal A and B). We showed that all nine prognostic signatures exhibited a similar prognostic performance in the entire dataset. Their prognostic abilities are due mostly to the detection of proliferation activity. Although ER- status (basal-like) and ERBB2+ expression status correspond to bad outcome, they seem to act through elevated expression of proliferation genes and thus contain only indirect information about prognosis. Clinical variables measuring the extent of tumor progression, such as tumor size and nodal status, still add independent prognostic information to proliferation genes. CONCLUSION: This meta-analysis unifies various results of previous gene expression studies in breast cancer. It reveals connections between traditional prognostic factors, expression-based subtyping, and prognostic signatures, highlighting the important role of proliferation in breast cancer prognosis.

Na(+)-K(+)-ATPase gene expression during in vitro development of rat fetal forebrain.

The existence of at least three isoforms of Na(+)-K(+)-ATPase in adult brain tissues [alpha 1, kidney type; alpha 2 [or alpha(+)]; alpha 3] suggests that these genes might be regulated in a cell-specific and time-dependent manner during development. We have studied this question in serum-free aggregating cell cultures of mechanically dissociated rat fetal telencephalon. At the protein level, the relative rate of synthesis of the pool of alpha 1-, alpha 2-, and alpha 3-subunits increased approximately twofold over 15 days of culture, leading to a marked increase in the immunochemical pool of alpha-subunits as measured by a panspecific polyclonal antibody. Concomitantly, Na(+)-K(+)-ATPase enzyme-specific activity increased three- (lower forebrain) to sixfold (upper forebrain). The transcripts of all three alpha-isoforms and beta-subunit were detected in vitro in similar proportion to the level observed in vivo. alpha 3-mRNA (3.7 kb) was more abundant than alpha 1 (3.7 kb) or alpha 2 (5.3 and 3.4 kb). Cytosine arabinoside (0.4 microM) and cholera toxin (0.1 microM) were used to selectively eliminate glial cells or neurons, respectively. It was found that alpha 2-mRNA is predominantly transcribed in glial cell cultures, whereas alpha 3- and beta 1-mRNA (2.7, 2.3, and 1.8 kb) are predominant in neuronal cultures.

" Arabidopsis Thaliana " response to insect feeding : new components controlling defense gene expression and plant resistance

Abstract: Plants cannot run away to escape attacking herbivores, but they defend themselves by producing anti-digestive proteins and toxic compounds (for example glucosinolates). The first goal of this thesis was to study changes in gene expression after insect attack using microarrays. The responses of Arabidopsis thaliana to feeding by the specialist Pieris rapae and the generalist Spodoptera liffora is were compared. We found that the transcript profiles after feeding by the two chewing insects were remarkably similar, although the generalist induced a slightly stronger response. The second goal was to evaluate the implication of the four signals jasmonic acid (JA), salicylic acid (SA), ethylene (ET), and abscisic acid (ABA) in the control of insect-regulated gene expression. Using signaling mutants, we observed that JA was the predominant signal and that ABA modulated defense gene expression. In contrast, SA and ET appeared to control slightly gene expression, but only after feeding by S. litforalis. The third goal was to establish whether plant responses are really effective against insects. In accordance with the transcript profile, both insects were affected by the JA-dependent defenses, as they performed better on the JA-insensitive mutant. S. littoralis also performed better on ABA-deficient mutants, providing evidence for the role of ABA in defense against insects. When testing indole or aliphatic glucosinolate deficient mutants, we found that they were also more susceptible to insect feeding, providing some of the first genetic evidence for the defensive role of glucosinolates in planta. Finally, a glutathione-deficient mutant, pad2-1, was also more susceptible to insect feeding and we could attribute this phenotype to a lowered accumulation of the major indole glucosinolate. In this thesis, we provide a comprehensive list of insect-regulated genes, including many transcription factors that constitute interesting candidate genes for the further study of insect-induced expression changes. Understanding how the plant responses to insects are regulated will provide tools for a better management of insect pest in the field. Résumé: Les plantes ne peuvent s'échapper pour fuir les insectes qui les attaquent, mais elles se défendent en produisant des protéines anti-digestives et des composés toxiques (par exemple des glucosinolates). Le premier but de cette thèse était d'étudier les changements de l'expression génétique lors d'attaque par des insectes en utilisant des puces à ADN. Nous avons comparé la réponse d'Arabidopsis thaliana à deux espèces d'insectes avec des habitudes alimentaires différentes : le spécialiste Pieris rapae et le généraliste Spodoptera littoralis. Nous avons trouvé que les profils de transcription après l'attaque par les deux insectes sont remarquablement similaires, bien que le généraliste induise une réponse légèrement plus forte. Le deuxième but était de déterminer l'implication de quatre signaux dans le contrôle de la réponse :l'acide jasmonique (JA), l'acide salicylique (SA), l'éthylène (ET), et l'acide abscissique (ABA). En utilisant de mutants de signalisation, nous avons montré que l'acide jasmonique était le signal prédominant et que l'acide abscissique modulait l'expression génétique. D'autre part, l'acide salicylique et l'éthylène contrôlent à un degré moindre l'expression génétique, mais seulement après l'attaque par S. littoralís. Le troisième but était d'établir si les réponses des plantes sont efficaces contre les insectes. En accord avec le profil de transcription, les deux espèces d'insectes se sont mieux développées sur un mutant insensible au JA, indiquant que les défenses contrôlées par ce signal sont cruciales pour la plante. De plus, les larves de S. littorales se sont mieux développées sur des mutants déficients en ABA, ce qui fournit une preuve du rôle de l'acide abscissique dans la défense contre les insectes. En testant des mutants déficients en glucosinolates de type indole ou aliphatique, nous avons trouvé qu'ils étaient plus sensibles aux insectes, démontrant ainsi le rôle défensif des glucosinolates in planta. Finalement, le mutant déficient en glutathion pad2-1 était aussi plus sensible à l'attaque des insectes, et nous avons pu attribuer ce phénotype à une plus faible augmentation d'un indole glucosinolate dans ce mutant. Dans cette thèse, nous avons mis en évidence un nombre important de gènes contrôlés par les insectes, comprenant de nombreux facteurs de transcription qui constituent des candidats intéressants pour`étudier plus en détail les changements d'expression génétique induits par les insectes. Une meilleure compréhension de la réponse des plantes contre l'attaque des insectes devrait nous permettre de développer de nouvelles stratégies pour mieux gérer les ravageurs des cultures.

The mammalian circadian timing system: from gene expression to physiology.

Many physiological processes in organisms from bacteria to man are rhythmic, and some of these are controlled by self-sustained oscillators that persist in the absence of external time cues. Circadian clocks are perhaps the best characterized biological oscillators and they exist in virtually all light-sensitive organisms. In mammals, they influence nearly all aspects of physiology and behavior, including sleep-wake cycles, cardiovascular activity, endocrinology, body temperature, renal activity, physiology of the gastro-intestinal tract, and hepatic metabolism. The master pacemaker is located in the suprachiasmatic nuclei, two small groups of neurons in the ventral part of the hypothalamus. However, most peripheral body cells contain self-sustained circadian oscillators with a molecular makeup similar to that of SCN (suprachiasmatic nucleus) neurons. This organization implies that the SCN must synchronize countless subsidiary oscillators in peripheral tissues, in order to coordinate cyclic physiology. In this review, we will discuss some recent studies on the structure and putative functions of the mammalian circadian timing system, but we will also point out some apparent inconsistencies in the currently publicized model for rhythm generation.

A comprehensive analysis of human gene expression profiles identifies stromal immunoglobulin κ C as a compatible prognostic marker in human solid tumors.

PURPOSE: Although the central role of the immune system for tumor prognosis is generally accepted, a single robust marker is not yet available. EXPERIMENTAL DESIGN: On the basis of receiver operating characteristic analyses, robust markers were identified from a 60-gene B cell-derived metagene and analyzed in gene expression profiles of 1,810 breast cancer; 1,056 non-small cell lung carcinoma (NSCLC); 513 colorectal; and 426 ovarian cancer patients. Protein and RNA levels were examined in paraffin-embedded tissue of 330 breast cancer patients. The cell types were identified with immunohistochemical costaining and confocal fluorescence microscopy. RESULTS: We identified immunoglobulin κ C (IGKC) which as a single marker is similarly predictive and prognostic as the entire B-cell metagene. IGKC was consistently associated with metastasis-free survival across different molecular subtypes in node-negative breast cancer (n = 965) and predicted response to anthracycline-based neoadjuvant chemotherapy (n = 845; P < 0.001). In addition, IGKC gene expression was prognostic in NSCLC and colorectal cancer. No association was observed in ovarian cancer. IGKC protein expression was significantly associated with survival in paraffin-embedded tissues of 330 breast cancer patients. Tumor-infiltrating plasma cells were identified as the source of IGKC expression. CONCLUSION: Our findings provide IGKC as a novel diagnostic marker for risk stratification in human cancer and support concepts to exploit the humoral immune response for anticancer therapy. It could be validated in several independent cohorts and carried out similarly well in RNA from fresh frozen as well as from paraffin tissue and on protein level by immunostaining.

Concordance among digital gene expression, microarrays, and qPCR when measuring differential expression of microRNAs.

Profiling microRNA (miRNA) expression is of widespread interest given the critical role of miRNAs in many cellular functions. Profiling can be achieved via hybridization-based (microarrays), sequencing-based, or amplification-based (quantitative reverse transcription-PCR, qPCR) technologies. Among these, microarrays face the significant challenge of accurately distinguishing between mature and immature miRNA forms, and different vendors have developed different methods to meet this challenge. Here we measure differential miRNA expression using the Affymetrix, Agilent, and Illumina microarray platforms, as well as qPCR (Applied Biosystems) and ultra high-throughput sequencing (Illumina). We show that the differential expression measurements are more divergent when the three types of microarrays are compared than when the Agilent microarray, qPCR, and sequencing technology measurements are compared, which exhibit a good overall concordance.

Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species.

Ants provide remarkable examples of equivalent genotypes developing into divergent and discrete phenotypes. Diploid eggs can develop either into queens, which specialize in reproduction, or workers, which participate in cooperative tasks such as building the nest, collecting food, and rearing the young. In contrast, the differentiation between males and females generally depends upon whether eggs are fertilized, with fertilized (diploid) eggs giving rise to females and unfertilized (haploid) eggs giving rise to males. To obtain a comprehensive picture of the relative contributions of gender (sex), caste, developmental stage, and species divergence to gene expression evolution, we investigated gene expression patterns in pupal and adult queens, workers, and males of two species of fire ants, Solenopsis invicta and S. richteri. Microarray hybridizations revealed that variation in gene expression profiles is influenced more by developmental stage than by caste membership, sex, or species identity. The second major contributor to variation in gene expression was the combination of sex and caste. Although workers and queens share equivalent diploid nuclear genomes, they have highly distinctive patterns of gene expression in both the pupal and the adult stages, as might be expected given their extraordinary level of phenotypic differentiation. Overall, the difference in the proportion of differentially expressed genes was greater between workers and males than between workers and queens or queens and males, consistent with the fact that workers and males share neither gender nor reproductive capability. Moreover, between-species comparisons revealed that the greatest difference in gene expression patterns occurred in adult workers, a finding consistent with the fact that adult workers most directly experience the distinct external environments characterizing the different habitats occupied by the two species. Thus, much of the evolution of gene expression in ants may occur in the worker caste, despite the fact that these individuals are largely or completely sterile. Analyses of gene expression evolution revealed a combination of positive selection and relaxation of stabilizing selection as important factors driving the evolution of such genes.

Detection of plant-modulated alterations in antifungal gene expression in Pseudomonas fluorescens CHA0 on roots by flow cytometry.

The biocontrol activity of the root-colonizing Pseudomonas fluorescens strain CHA0 is largely determined by the production of antifungal metabolites, especially 2,4-diacetylphloroglucinol. The expression of these metabolites depends on abiotic and biotic environmental factors, in particular, elements present in the rhizosphere. In this study, we have developed a new method for the in situ analysis of antifungal gene expression using flow cytometry combined with green fluorescent protein (GFP)-based reporter fusions to the phlA and prnA genes essential for the production of the antifungal compounds 2,4-diacetylphloroglucinol and pyrrolnitrin, respectively, in strain CHA0. Expression of phlA-gfp and prnA-gfp in CHA0 cells harvested from the rhizosphere of a set of plant species as well as from the roots of healthy, leaf pathogen-attacked, and physically stressed plants were analyzed using a FACSCalibur. After subtraction of background fluorescence emitted by plant-derived particles and CHA0 cells not carrying the gfp reporters, the average gene expression per bacterial cell could be calculated. Levels of phlA and prnA expression varied significantly in the rhizospheres of different plant species. Physical stress and leaf pathogen infection lowered phlA expression levels in the rhizosphere of cucumber. Our results demonstrate that the newly developed approach is suitable to monitor differences in levels of antifungal gene expression in response to various plant-derived factors. An advantage of the method is that it allows quantification of bacterial gene expression in rhizosphere populations at a single-cell level. To our best knowledge, this is the first study using flow cytometry for the in situ analysis of biocontrol gene expression in a plant-beneficial bacterium in the rhizosphere.

Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.

Hypermethylation-mediated regulation of CD44 gene expression in human neuroblastoma

The CD44 adhesion receptor is silenced in highly malignant neuroblastomas (NBs) with MYCN amplification. Because its functional expression is associated with decreased tumorigenic properties, CD44 behaves as a tumor suppressor gene in NB and other cancers. Given that the precise mechanisms responsible for CD44 silencing are not elucidated, we investigated whether CD44 expression could be regulated by DNA hypermethylation. The methylation status of CD44 gene promoter and exon 1 regions was analyzed in 12 NB cell lines and 21 clinical samples after bisulfite genomic modification, followed by PCR and single-strand conformation polymorphism analysis and genomic sequencing. The results showed that almost all CD44-negative cell lines displayed hypermethylation in both regions, whereas all CD44-expressing cell lines were unmethylated. These observations correlated with the ability to restore CD44 mRNA and protein expression by treatment of CD44-negative cells with the 5-aza-2'-deoxycytidine demethylating agent. In contrast, no CD44 gene hypermethylation could be detected in 21 NB clinical samples of different stages, irrespective of CD44 expression. Although our results suggest that aberrant methylation of promoter and exon 1 regions is involved in CD44 silencing in NB cell lines, they also indicate that methylation of unidentified regulatory sequences or methylation-independent mechanisms also control the expression of CD44 in primary NB tumors and cell lines. We therefore conclude that CD44 silencing is controlled by complex and tumor cell-specific processes, including gene hypermethylation. Further investigation of other mechanisms and genes involved in CD44 regulation will be needed before demethylation-mediated reactivation of the CD44 gene can be considered as therapeutic strategy for neuroblastoma and perhaps other related cancers.

Perinuclear Mlp proteins down-regulate gene expression in response to a defect in mRNA export

Résumé Une caractéristique des cellules eucaryotes est le confinement du matériel génétique (ADN/DNA) dans le noyau. Pour décoder cette information, un ARN messager (mRNA) est d'abord transcrit sous forme d'un ARN prémessager (pré-mRNA). Ce-dernier doit subir plusieurs étapes de maturation pour aboutir à une particule ribonucléoprotéique (mRNP) qui sera exportée vers le cytoplasme et traduite en protéine. La protéine de levure Mex67p et son homologue humain TAP sont des récepteurs d'export médiant la translocation du mRNP au travers des complexes du pore nucléaire (NPC). Mex67p/TAP ne se lient pas directement au mRNA, mais nécessitent la présence de protéines adaptatrices, telles que Yra1p et son homologue humain REF1. Afin d'identifier de nouveaux facteurs impliqués dans l'export des mRNPs ou de nouvelles fonctions pour Yra1p, nous avons effectué un crible génétique avec un mutant thermosensible de Yra1p, GFP-yra 1 -8. Ce mutant présente un défaut d'export des mRNAs et une diminution des niveaux de transcrits du gène rapporteur LacZ ainsi que de certains transcrits endogènes. Nous avons trouvé que la perte de Mlp2p, ou d'une protéine hautement similaire, Mlp1p, restaure la croissance du mutant GFP-yra1-8 à température restrictive. Mlp1p et Mlp2p sont des protéines nucléaires, dont l'homologue humain est TPR. Les Mlp (myosin¬like proteins) ainsi que TPR forment des structures filamenteuses ancrées aux NPC. Bien que la fonction des Mlp ne soit pas clairement définie, un rôle dans la biogenèse et la surveillance des mRNPs a été récemment proposé. Notre étude montre que la perte des Mlp, non seulement restaure la croissance de GFP-yra1-8, mais augmente aussi les niveaux des transcrits LacZ et facilite leur apparition dans le cytoplasme. Des expériences d'immunoprécipitations de la chromatine révèlent que Mlp2p diminue le taux de synthèse du transcrit LacZ dans GFP-yra1-8. Des analyses du transcriptome montrent que Mlp2p réduit aussi les niveaux d'une population de transcrits endogènes dans le mutant. Finalement, des localisations in situ suggèrent que la transcription du rapporteur LacZ a lieu à la périphérie du noyau, à proximité des Mlp. Ainsi, les protéines Mlp pourraient préférentiellement diminuer la transcription de gènes exprimés à la périphérie nucléaire. Nous montrons aussi que Yra1p interagit génétiquement avec Nab2p une protéine liée au mRNA et impliquée dans son export, mais non avec d'autres protéines également impliquées dans l'export des mRNAs. Les résultats obtenus soutiennent un modèle où les protéines Yra1p et Nab2p sont nécessaires à l'arrimage des mRNPs sur la plate-forme des Mlp. Si ces signaux manquent ou sont défectueux, les mRNPs ne peuvent pas poursuivre leur trajet vers le canal central du NPC. Ce bloc induirait par la suite une diminution de la transcription d'une population de gènes potentiellement localisée à la périphérie nucléaire. Dans son ensemble, cette étude suggère que les protéines Mlp établissent un lien entre la transcription de certains mRNAs et leur export au travers du pore nucléaire. Summary A hallmark of the eukaryotic cell is the packaging of DNA in the nucleus. To decode the genetic information, a messenger RNA (mRNA) is first synthesized as a pre-mRNA molecule, which undergoes different maturation steps resulting in an mRNP (messenger RNA ribonucleoprotein), which can be actively transported to the cytoplasm and translated into a protein. Yeast Mex67p and its human homologue TAP are export receptors mediating mRNP translocation through the nuclear pore complex (NPC). The recruitment of Mex67p/TAP to mRNA is mediated by mRNA export adaptors of the evolutionarily conserved REF (RNA and Export Factor binding) family: yeast Yra1p and human REF1. To uncover new functions of Yra1p or new factors implicated in mRNA export, we performed a genetic screen with a themiosensitive (ts) yra1 mutant, GFP-yra1-8. This mutant exhibits mRNA export defects and a decrease in the levels of LacZ reporter and certain endogenous transcripts. We found that the loss of Mlp2p, or the related Mlp1p protein, substantially rescues the growth defect of the GFP-yra1 -8 mutant. Mlp1p and M1p2p are large non-essential proteins, homologous to human TPR, proposed to form intra-nuclear filamentous structures anchored at the NPC. Their role is not clearly defined, but they have been implicated in mRNP biogenesis and surveillance. Our study shows that loss of Mlp proteins not only restores growth of GFP-yra1-8, but also rescues LacZ mRNA levels and increases their appearance in the cytoplasm. Chromatin immunoprecipitation and pulse chase experiments indicate that Mlp2p down-regulates LacZ mRNA synthesis in GFP-yra1-8. DNA micro- array analyses reveal that Mlp2p also reduces the levels of a subset of cellular transcripts in the yra1 mutant strain. In situ localizations suggest that LacZ transcription occurs at the nuclear periphery, in close proximity to Mlp proteins. Thus, Mlp proteins may preferentially down-regulate genes expressed at the nuclear periphery. Finally, we show that Yra1p genetically interacts with the shuttling mRNA-binding protein Nab2p and that loss of Mlp proteins rescues the growth defect of yra1 and nab2, but not other mRNA export mutants. The data support a model in which Nab2p and Yra1p are required for rnRNP docking to the Mlp platform. Lack of these signals prevents mRNPs from crossing the Mlp gate. This block may then negatively feed-back on the transcription of a subset of genes, potentially located at the nuclear envelope. Overall, this study suggests that perinuclear Mlp proteins establish a link between mRNA transcription and export.

Sensitive gene expression profiling of human T cell subsets reveals parallel post-thymic differentiation for CD4+ and CD8+ lineages.

The differentiation of CD4(+) or CD8(+) T cells following priming of naive cells is central in the establishment of the immune response against pathogens or tumors. However, our understanding of this complex process and the significance of the multiple subsets of differentiation remains controversial. Gene expression profiling has opened new directions of investigation in immunobiology. Nonetheless, the need for substantial amount of biological material often limits its application range. In this study, we have developed procedures to perform microarray analysis on amplified cDNA from low numbers of cells, including primary T lymphocytes, and applied this technology to the study of CD4 and CD8 lineage differentiation. Gene expression profiling was performed on samples of 1000 cells from 10 different subpopulations, defining the major stages of post-thymic CD4(+) or CD8(+) T cell differentiation. Surprisingly, our data revealed that while CD4(+) and CD8(+) T cell gene expression programs diverge at early stages of differentiation, they become increasingly similar as cells reach a late differentiation stage. This suggests that functional heterogeneity between Ag experienced CD4(+) and CD8(+) T cells is more likely to be located early during post-thymic differentiation, and that late stages of differentiation may represent a common end in the development of T-lymphocytes.

Light-regulated interactions with SPA proteins underlie cryptochrome-mediated gene expression.

Cryptochromes are a class of photosensory receptors that control important processes in animals and plants primarily by regulating gene expression. How photon absorption by cryptochromes leads to changes in gene expression has remained largely elusive. Three recent studies, including Lian and colleagues (pp. 1023-1028) and Liu and colleagues (pp. 1029-1034) in this issue of Genes & Development, demonstrate that the interaction of light-activated Arabidopsis cryptochromes with a class of regulatory components of E3 ubiquitin ligase complexes leads to environmentally controlled abundance of transcriptional regulators.