On the simultaneous use of clinical and microarray expression data in the cluster analysis of tissue samples

This paper considers a model-based approach to the clustering of tissue samples of a very large number of genes from microarray experiments. It is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. Frequently in practice, there are also clinical data available on those cases on which the tissue samples have been obtained. Here we investigate how to use the clinical data in conjunction with the microarray gene expression data to cluster the tissue samples. We propose two mixture model-based approaches in which the number of components in the mixture model corresponds to the number of clusters to be imposed on the tissue samples. One approach specifies the components of the mixture model to be the conditional distributions of the microarray data given the clinical data with the mixing proportions also conditioned on the latter data. Another takes the components of the mixture model to represent the joint distributions of the clinical and microarray data. The approaches are demonstrated on some breast cancer data, as studied recently in van't Veer et al. (2002).

Tissue-specific gene expression in soybean (Glycine max) detected by cDNA microarray analysis

We have constructed cDNA microarrays for soybean (Glycine max L. Merrill), containing approximately 4,100 Unigene ESTs derived from axenic roots, to evaluate their application and utility for functional genomics of organ differentiation in legumes. We assessed microarray technology by conducting studies to evaluate the accuracy of microarray data and have found them to be both reliable and reproducible in repeat hybridisations. Several ESTs showed high levels (>50 fold) of differential expression in either root or shoot tissue of soybean. A small number of physiologically interesting, and differentially expressed sequences found by microarray analysis were verified by both quantitative real-time RT-PCR and Northern blot analysis. There was a linear correlation (r(2) = 0.99, over 5 orders of magnitude) between microarray and quantitative real-time RT-PCR data. Microarray analysis of soybean has enormous potential not only for the discovery of new genes involved in tissue differentiation and function, but also to study the expression of previously characterised genes, gene networks and gene interactions in wild-type, mutant or transgenic; plants.

Effects of A1 adenosine receptor overexpression on normoxic and post-ischemic gene expression

Objectives: To identify potential molecular genetic determinants of cardiovascular ischemic tolerance in wild-type and transgenic hearts overexpressing A(1) adenosine receptors (A(1)ARs). Methods: cDNA microarrays were used to explore expression of 1824 genes ill wild-type hearts and ischemia-tolerant mouse hearts overexpressing A(1)ARs. Results: Overexpression of A(1)ARs reduced post-ischemic contractile dysfunction, limited arrhythmogenesis, and reduced necrosis by similar to80% in hearts subjected to 30 min global ischemia 60 mill reperfusion. Cardioprotection was abrogated by acute A(1)AR antagonism, and only a small number (19) of genes were modified by A(1)AR overexpression in normoxic hearts. Ischemia-reperfusion significantly altered expression of 75 genes in wild-type hearts (14 induced, 61 down-regulated), including genes for metabolic enzymes, structural/motility proteins, cell signaling proteins, defense/growth proteins, and regulators of transcription and translation. A(1)AR overexpression reversed the majority of gene down-regulation whereas gene induction was generally unaltered. Additionally, genes involved in cell defence, signaling and gene expression were selectively modified by ischemia in transgenic hearts (33 induced, 10 down-regulated), possibly contributing to the protected phenotype. Real-time PCR verified changes in nine selected genes, revealing concordance with array data. Transcription of the A(1)AR gene was also modestly reduced post-ischemia, consistent with impaired functional sensitivity to A(1)AR stimulation Conclusions: Data are presented regarding the early post-ischemic gene profile of intact heart. Reduced A(1)AR transcription is observed which may contribute to poor outcome from ischemia. A(1)AR overexpression selectively modifies post-ischemic gene expression, potentially contributing to ischemic-tolerance. (C) 2003 European Society of Cardiology. Published by Elsevier Science B.V. All rights reserved.

Engineering mRNA translation initiation to enhance transient gene expression in Chinese hamster ovary cells

To increase transient expression of recombinant proteins in Chinese hamster ovary cells, we have engineered their protein synthetic capacity by directed manipulation of mRNA translation initiation. To control this process we constructed a nonphosphorylatable Ser51Ala site-directed mutant of eIF2, a subunit of the trimeric eIF2 complex that is implicated in regulation of the global rate of mRNA translation initiation in eukaryotic cells. Phosphorylation of eIF2 by protein kinases inhibits eIF2 activity and is known to increase as cells perceive a range of stress conditions. Using single-and dual-gene plasmids introduced into CHO cells by electroporation, we found that transient expression of the eIF2 Ser51Ala mutant with firefly luciferase resulted in a 3-fold increase in reporter activity, relative to cells transfected with reporter only. This effect was maintained in transfected cells for at least 48 h after transfection. Expression of the wild-type eIF2 protein had no such effect. Elevated luciferase activity was associated with a reduction in the level of eIF2 phosphorylation in cells transfected with the mutant eIF2 construct. Transfection of CHO cells with the luciferase-only construct resulted in a marked decrease in the global rate of protein synthesis in the whole cell population 6 h post-transfection. However, expression of the mutant Ser51Ala or wild-type eIF2 proteins restored the rate of protein synthesis in transfected cells to a level equivalent to or exceeding that of control cells. Associated with this, entry of plasmid DNA into cells during electroporation was visualized by confocal microscopy using a rhodamine-labeled plasmid construct expressing green fluorescent protein. Six hours after transfection, plasmid DNA was present in all cells, albeit to a variable extent. These data suggest that entry of naked DNA into the cell itself functions to inhibit protein synthesis by signaling mechanisms affecting control of mRNA translation by eIF2. This work therefore forms the basis of a rational strategy to generically up-regulate transient expression of recombinant proteins by simultaneous host cell engineering.

Genetic study of alcoholism and novel gene expression in the alcoholic brain

Alcohol dependence may result from neuroadaptation involving alteration of gene expression after long-term alcohol exposure. The systematic study of gene expression profiles of the human alcoholic brain was initiated using the method of polymerase chain reaction (PCR)-differential display and was followed by DNA microarray. To date, more than 100 alcohol-responsive genes have been identified from the frontal cortex, motor cortex and nucleus accumbens of the human brain. These genes have a wide range of functions in the brain and indicate diverse actions of alcohol on neuronal function. This review discusses the current information on the genetic basis of alcoholism and the induction and characterization of these alcohol-responsive genes.

Correlating gene expression with larval competence and the effect of age and parentage on metamorphosis in the tropical abalone Haliotis asinina

The non-geniculate crustose coralline alga (CCA) Mastophora pacifica can induce the metamorphosis of competent Haliotis asinina (Vetigastropoda) larvae. The ability to respond to this natural cue varies considerably with larval age, with a higher proportion of older larvae (e.g. 90 h) able to metamorphose in response to M. pacifica than younger larvae (e.g. 66 h). Here we document the variation in time to acquisition of competence within a larval age class. For example, after 18 h of exposure to M. pacifica, approximately 15 and 36% of 84 and 90-h-old H. asinina larvae had initiated metamorphosis, respectively. This age-dependent response to M. pacifica is also observed when different aged larvae are exposed to CCA for varying periods. A higher proportion of older larvae require shorter periods of exposure to CCA than younger larvae in order to initiate metamorphosis. In this experiment, as in the previous, a small proportion of young larvae were able to respond to brief periods of CCA exposure, suggesting that they had developed the same state of competency as the majority of their older counterparts. Comparisons of the proportions of larvae undergoing metamorphosis between families reveals that parentage also has a significant (P < 0.05) affect on whether an individual will initiate metamorphosis at a given age. These familial differences are more pronounced when younger, largely pre-competent larvae (i.e. 66 h old) are exposed to M. pacifica, with proportions of larvae undergoing metamorphosis differing by as much as 10 fold between families. As these data suggest that variation in the rate of development of the competent state has a genetic basis, and as a first step towards identifying the molecular basis to this variation, we have identified numerous genes that are differentially expressed later in larval development using a differential display approach. Spatial expression analysis of these genes suggests that they may be directly involved in the acquisition of competence, or may play a functional role in the postlarva following metamorphosis.

Regulation of Amh during sex determination in chickens: Sox gene expression in male and female gonads

During mammalian sexual development, the SOX9 transcription factor up-regulates expression of the gene encoding anti-Mullerian hormone (AMH), but in chickens, Sox9 gene expression reportedly occurs after the onset of Amh expression. Here, we examined expression of the related gene Sox8 in chicken embryonic gonads during the sex-determining period. We found that cSox8 is expressed at similar levels in both sexes at embryonic day 6 and 7, and only at the anterior tip of the gonad, suggesting that SOX8 is not responsible for the sex-specific increase in cAmh gene expression at these stages. We also found that several other chicken Sox genes (cSox3, cSox4 and cSox11) are expressed in embryonic gonads, but at similar levels in both sexes. Our data suggest that the molecular mechanisms involved in the regulation of Amh genes of mouse and chicken are not conserved, despite similar patterns of Amh expression in both species.

Gene expression profiles in murine hematopoietic stem cells revisited: Analysis of cDNA libraries reveals high levels of translational and metabolic activities

Gene expression studies from hematopoietic stem cell (HSC) populations purified to variable degrees have defined a set of sternness genes. Unexpectedly, results also hinted toward a HSC chromatin poised in a wide-open state. With the aim of providing a robust tool for further studies into the molecular biology of HSCs, the studies herein describe the construction and comparative molecular analysis of A-phage cDNA libraries from highly purified HSCs that retained their long-term repopulating activities (long-term HSCs [LT-HSCs]) and from short-term repopulating HSCs that were largely depleted of these activities. Microarray analysis of the libraries confirmed the previous results but also revealed an unforeseen preferential expression of translation- and metabolism-associated genes in the LT-HSCs. Therefore, these data indicate that HSCs are quiescent only in regard of proliferative activities but are in a state of readiness to provide the metabolic and translational activities required after induction of proliferation and exit from the HSC pool.

Over-expression of Eph and ephrin genes in advanced ovarian cancer: ephrin gene expression correlates with shortened survival

Background: Increased expression of Eph receptor tyrosine kinases and their ephrin ligands has been implicated in tumor progression in a number of malignancies. This report describes aberrant expression of these genes in ovarian cancer, the commonest cause of death amongst gynaecological malignancies. Methods: Eph and ephrin expression was determined using quantitative real time RT-PCR. Correlation of gene expression was measured using Spearman's rho statistic. Survival was analysed using log-rank analysis and ( was visualised by) Kaplan-Meier survival curves. Results: Greater than 10 fold over-expression of EphA1 and a more modest over-expression of EphA2 were observed in partially overlapping subsets of tumors. Over-expression of EphA1 strongly correlated ( r = 0.801; p < 0.01) with the high affinity ligand ephrin A1. A similar trend was observed between EphA2 and ephrin A1 ( r = 0.387; p = 0.06). A striking correlation of both ephrin A1 and ephrin A5 expression with poor survival ( r = - 0.470; p = 0.02 and r = - 0.562; p < 0.01) was observed. Intriguingly, there was no correlation between survival and other clinical parameters or Eph expression. Conclusion: These data imply that increased levels of ephrins A1 and A5 in the presence of high expression of Ephs A1 and A2 lead to a more aggressive tumor phenotype. The known functions of Eph/ephrin signalling in cell de-adhesion and movement may explain the observed correlation of ephrin expression with poor prognosis.

A Mixture model with random-effects components for clustering correlated gene-expression profiles

Motivation: The clustering of gene profiles across some experimental conditions of interest contributes significantly to the elucidation of unknown gene function, the validation of gene discoveries and the interpretation of biological processes. However, this clustering problem is not straightforward as the profiles of the genes are not all independently distributed and the expression levels may have been obtained from an experimental design involving replicated arrays. Ignoring the dependence between the gene profiles and the structure of the replicated data can result in important sources of variability in the experiments being overlooked in the analysis, with the consequent possibility of misleading inferences being made. We propose a random-effects model that provides a unified approach to the clustering of genes with correlated expression levels measured in a wide variety of experimental situations. Our model is an extension of the normal mixture model to account for the correlations between the gene profiles and to enable covariate information to be incorporated into the clustering process. Hence the model is applicable to longitudinal studies with or without replication, for example, time-course experiments by using time as a covariate, and to cross-sectional experiments by using categorical covariates to represent the different experimental classes. Results: We show that our random-effects model can be fitted by maximum likelihood via the EM algorithm for which the E(expectation) and M(maximization) steps can be implemented in closed form. Hence our model can be fitted deterministically without the need for time-consuming Monte Carlo approximations. The effectiveness of our model-based procedure for the clustering of correlated gene profiles is demonstrated on three real datasets, representing typical microarray experimental designs, covering time-course, repeated-measurement and cross-sectional data. In these examples, relevant clusters of the genes are obtained, which are supported by existing gene-function annotation. A synthetic dataset is considered too.

Cell cycle alterations in biopsied olfactory neuroepithelium in schizophrenia and bipolar I disorder using cell culture and gene expression analyses

We previously demonstrated that olfactory cultures front individuals with schizophrenia had increased cell proliferation compared to Cultures from healthy controls. The aims of this study were to (a) replicate this observation in a new group Of individuals with schizophrenia, (b) examine the specificity of these findings by including individuals with bipolar I disorder and (c) explore gene expression differences that may underlie cell cycle differences in these diseases. Compared to controls (n = 10), there was significantly more mitosis in schizophrenia patient cultures (it = 8) and significantly more cell death in the bipolar I disorder patient cultures (n=8). Microarray data showed alterations to the cell cycle and phosphatidylinositol signalling pathways in schizophrenia and bipolar I disorder, respectively. Whilst caution is required in the interpretation of the array results, the study provides evidence indicating that cell proliferation and cell death in olfactory neuroepithelial cultures is differentially altered in schizophrenia and bipolar disorder. (c) 2005 Elsevier B.V. All rights reserved.

Comparative gene expression in brain regions of human alcoholics

The mesocorticolimbic system is the reward centre of the brain and the major target for drugs of abuse including alcohol. Neuroadaptive changes in this region are thought to underlie the process of tolerance and dependence. Recently, several research groups have searched for alcohol-responsive genes using high-throughput microarrays and well-characterized human post-mortem material. Comparison of data from these studies of cortical regions highlights the differences in experimental approach and selection of cases. However, alcohol-responsive gene sets associated with transcription, oxidative stress and energy production were common to these studies. In marked contrast, alcohol-responsive genes in the nucleus accumbens and the ventral tegmental area are primarily associated with changes in neurotransmission and signal transduction. These data support the concept that, within cortical regions, changes in gene expression are associated with alcoholism-related pathology. In the dopaminergic tract of the mesocorticolimbic system, alcohol-responsive gene sets suggest long-term neuroplastic changes in synaptic transmission.

Transcriptome analysis of a cnidarian-dinoflagellate mutualism reveals complex modulation of host gene expression

Background: Cnidarian - dinoflagellate intracellular symbioses are one of the most important mutualisms in the marine environment. They form the trophic and structural foundation of coral reef ecosystems, and have played a key role in the evolutionary radiation and biodiversity of cnidarian species. Despite the prevalence of these symbioses, we still know very little about the molecular modulators that initiate, regulate, and maintain the interaction between these two different biological entities. In this study, we conducted a comparative host anemone transcriptome analysis using a cDNA microarray platform to identify genes involved in cnidarian - algal symbiosis. Results: We detected statistically significant differences in host gene expression profiles between sea anemones ( Anthopleura elegantissima) in a symbiotic and non-symbiotic state. The group of genes, whose expression is altered, is diverse, suggesting that the molecular regulation of the symbiosis is governed by changes in multiple cellular processes. In the context of cnidarian dinoflagellate symbioses, we discuss pivotal host gene expression changes involved in lipid metabolism, cell adhesion, cell proliferation, apoptosis, and oxidative stress. Conclusion: Our data do not support the existence of symbiosis- specific genes involved in controlling and regulating the symbiosis. Instead, it appears that the symbiosis is maintained by altering expression of existing genes involved in vital cellular processes. Specifically, the finding of key genes involved in cell cycle progression and apoptosis have led us to hypothesize that a suppression of apoptosis, together with a deregulation of the host cell cycle, create a platform that might be necessary for symbiont and/or symbiont-containing host cell survival. This first comprehensive molecular examination of the cnidarian - dinoflagellate associations provides critical insights into the maintenance and regulation of the symbiosis.

Real-Time Gene-Expression Extraction with Fast Classification (FC) Networks

Fast Classification (FC) networks were inspired by a biologically plausible mechanism for short term memory where learning occurs instantaneously. Both weights and the topology for an FC network are mapped directly from the training samples by using a prescriptive training scheme. Only two presentations of the training data are required to train an FC network. Compared with iterative learning algorithms such as Back-propagation (which may require many hundreds of presentations of the training data), the training of FC networks is extremely fast and learning convergence is always guaranteed. Thus FC networks may be suitable for applications where real-time classification is needed. In this paper, the FC networks are applied for the real-time extraction of gene expressions for Chlamydia microarray data. Both the classification performance and learning time of the FC networks are compared with the Multi-Layer Proceptron (MLP) networks and support-vector-machines (SVM) in the same classification task. The FC networks are shown to have extremely fast learning time and comparable classification accuracy.