Multi-objective clustering ensemble for gene expression data analysis

In this paper, we present an algorithm for cluster analysis that integrates aspects from cluster ensemble and multi-objective clustering. The algorithm is based on a Pareto-based multi-objective genetic algorithm, with a special crossover operator, which uses clustering validation measures as objective functions. The algorithm proposed can deal with data sets presenting different types of clusters, without the need of expertise in cluster analysis. its result is a concise set of partitions representing alternative trade-offs among the objective functions. We compare the results obtained with our algorithm, in the context of gene expression data sets, to those achieved with multi-objective Clustering with automatic K-determination (MOCK). the algorithm most closely related to ours. (C) 2009 Elsevier B.V. All rights reserved.

Genomic deletions at 1p and 14q are associated with an abnormal cDNA microarray gene expression pattern in meningiomas but not in schwannomas

The molecular pathology of meningiomas and shwannomas involve the inactivation of the NF2 gene to generate grade I tumors. Genomic losses at 1p and 14q are observed in both neoplasms, although more frequently in meningiomas. The inactivation of unidentified genes located in these regions appears associated with tumor progression in meningiomas, but no clues to its molecular/clinical meaning are available in schwannomas. Recent microarray gene expression studies have demonstrated the existence of molecular subgroups in both entities. In the present study, we correlated the presence of genomic deletions at 1p, 14q, and 22q with the expression patterns of 96 tumor-related genes obtained by cDNA low-density microarrays in a series of 65 tumors including 42 meningiomas and 23 schwannomas. Two expression pattern groups were identified by cDNA mycroarray analysis when compared to the expression pattern in normal control RNA in both meningiomas and schwannomas, each one with patterns similar and different from the normal control. Meningioma and schwannoma subgroups differed in the expression of 38 and 16 genes, respectively. Using MLPA and microsatellites, we identified genomic losses at 1p, 14q, and 22q at nonrandom frequencies (12.5-69%) in meningiomas and schwannomas. Losses at 22q were almost equally frequent in both molecular expression subgroups in both neoplasms. However, deletions at 1p and 14q accumulated in meningiomas with a gene expression pattern different from the normal pattern, whereas the inverse situation occurred in schwannomas. Those anomalies characterized the schwannomas with expression pattern similar to the normal control. These findings suggest that deletions at 1p and 14q enhance the development of an abnormal tumor-related gene expression pattern in meningiomas, but this fact is not corroborated in schwannomas. (C) 2010 Elsevier Inc. All rights reserved.

Pre-processing for noise detection in gene expression classification data

Due to the imprecise nature of biological experiments, biological data is often characterized by the presence of redundant and noisy data. This may be due to errors that occurred during data collection, such as contaminations in laboratorial samples. It is the case of gene expression data, where the equipments and tools currently used frequently produce noisy biological data. Machine Learning algorithms have been successfully used in gene expression data analysis. Although many Machine Learning algorithms can deal with noise, detecting and removing noisy instances from the training data set can help the induction of the target hypothesis. This paper evaluates the use of distance-based pre-processing techniques for noise detection in gene expression data classification problems. This evaluation analyzes the effectiveness of the techniques investigated in removing noisy data, measured by the accuracy obtained by different Machine Learning classifiers over the pre-processed data.

Xylella fastidiosa gene expression analysis by DNA microarrays

Xylella fastidiosa genome sequencing has generated valuable data by identifying genes acting either on metabolic pathways or in associated pathogenicity and virulence. Based on available information on these genes, new strategies for studying their expression patterns, such as microarray technology, were employed. A total of 2,600 primer pairs were synthesized and then used to generate fragments using the PCR technique. The arrays were hybridized against cDNAs labeled during reverse transcription reactions and which were obtained from bacteria grown under two different conditions (liquid XDM2 and liquid BCYE). All data were statistically analyzed to verify which genes were differentially expressed. In addition to exploring conditions for X. fastidiosa genome-wide transcriptome analysis, the present work observed the differential expression of several classes of genes (energy, protein, amino acid and nucleotide metabolism, transport, degradation of substances, toxins and hypothetical proteins, among others). The understanding of expressed genes in these two different media will be useful in comprehending the metabolic characteristics of X. fastidiosa, and in evaluating how important certain genes are for the functioning and survival of these bacteria in plants.

A reliable measure of similarity based on dependency for short time series: an application to gene expression networks

Background: Microarray techniques have become an important tool to the investigation of genetic relationships and the assignment of different phenotypes. Since microarrays are still very expensive, most of the experiments are performed with small samples. This paper introduces a method to quantify dependency between data series composed of few sample points. The method is used to construct gene co-expression subnetworks of highly significant edges. Results: The results shown here are for an adapted subset of a Saccharomyces cerevisiae gene expression data set with low temporal resolution and poor statistics. The method reveals common transcription factors with a high confidence level and allows the construction of subnetworks with high biological relevance that reveals characteristic features of the processes driving the organism adaptations to specific environmental conditions. Conclusion: Our method allows a reliable and sophisticated analysis of microarray data even under severe constraints. The utilization of systems biology improves the biologists ability to elucidate the mechanisms underlying celular processes and to formulate new hypotheses.

Comprehensive gene expression profiling in lungs of mice infected with Mycobacterium tuberculosis following DNAhsp65 immunotherapy

Background The continued increase in tuberculosis (TB) rates and the appearance of extremely resistant Mycobacterium tuberculosis strains (XDR-TB) worldwide are some of the great problems of public health. In this context, DNA immunotherapy has been proposed as an effective alternative that could circumvent the limitations of conventional drugs. Nonetheless, the molecular events underlying these therapeutic effects are poorly understood. Methods We characterized the transcriptional signature of lungs from mice infected with M. tuberculosis and treated with heat shock protein 65 as a genetic vaccine (DNAhsp65) combining microarray and real-time polymerase chain reaction analysis. The gene expression data were correlated with the histopathological analysis of lungs. Results The differential modulation of a high number of genes allowed us to distinguish DNAhsp65-treated from nontreated animals (saline and vector-injected mice). Functional analysis of this group of genes suggests that DNAhsp65 therapy could not only boost the T helper (Th)1 immune response, but also could inhibit Th2 cytokines and regulate the intensity of inflammation through fine tuning of gene expression of various genes, including those of interleukin-17, lymphotoxin A, tumour necrosis factor-cl, interleukin-6, transforming growth factor-beta, inducible nitric oxide synthase and Foxp3. In addition, a large number of genes and expressed sequence tags previously unrelated to DNA-therapy were identified. All these findings were well correlated with the histopathological lesions presented in the lungs. Conclusions The effects of DNA therapy are reflected in gene expression modulation; therefore, the genes identified as differentially expressed could be considered as transcriptional biomarkers of DNAhsp65 immunotherapy against TB. The data have important implications for achieving a better understanding of gene-based therapies. Copyright (C) 2008 John Wiley & Sons, Ltd.

Sequential analysis of global gene expression profiles in immature and in vitro matured bovine oocytes: potential molecular markers of oocyte maturation

Background: Without intensive selection, the majority of bovine oocytes submitted to in vitro embryo production (IVP) fail to develop to the blastocyst stage. This is attributed partly to their maturation status and competences. Using the Affymetrix GeneChip Bovine Genome Array, global mRNA expression analysis of immature (GV) and in vitro matured (IVM) bovine oocytes was carried out to characterize the transcriptome of bovine oocytes and then use a variety of approaches to determine whether the observed transcriptional changes during IVM was real or an artifact of the techniques used during analysis. Results: 8489 transcripts were detected across the two oocyte groups, of which similar to 25.0% (2117 transcripts) were differentially expressed (p < 0.001); corresponding to 589 over-expressed and 1528 under-expressed transcripts in the IVM oocytes compared to their immature counterparts. Over expression of transcripts by IVM oocytes is particularly interesting, therefore, a variety of approaches were employed to determine whether the observed transcriptional changes during IVM were real or an artifact of the techniques used during analysis, including the analysis of transcript abundance in oocytes in vitro matured in the presence of a-amanitin. Subsets of the differentially expressed genes were also validated by quantitative real-time PCR (qPCR) and the gene expression data was classified according to gene ontology and pathway enrichment. Numerous cell cycle linked (CDC2, CDK5, CDK8, HSPA2, MAPK14, TXNL4B), molecular transport (STX5, STX17, SEC22A, SEC22B), and differentiation (NACA) related genes were found to be among the several over-expressed transcripts in GV oocytes compared to the matured counterparts, while ANXA1, PLAU, STC1and LUM were among the over-expressed genes after oocyte maturation. Conclusion: Using sequential experiments, we have shown and confirmed transcriptional changes during oocyte maturation. This dataset provides a unique reference resource for studies concerned with the molecular mechanisms controlling oocyte meiotic maturation in cattle, addresses the existing conflicting issue of transcription during meiotic maturation and contributes to the global goal of improving assisted reproductive technology.

Effects of the antimicrobial peptide gomesin on the global gene expression profile, virulence and biofilm formation of Xylella fastidiosa

In the xylem vessels of susceptible hosts, such as citrus trees, Xylella fastidiosa forms biofilm-like colonies that can block water transport, which appears to correlate to disease symptoms. Besides aiding host colonization, bacterial biofilms play an important role in resistance against antimicrobial agents, for instance antimicrobial peptides (AMPs). Here, we show that gomesin, a potent AMP from a tarantula spider, modulates X. fastidiosa gene expression profile upon 60 min of treatment with a sublethal concentration. DNA microarray hybridizations revealed that among the upregulated coding sequences, some are related to biofilm production. In addition, we show that the biofilm formed by gomesin-treated bacteria is thicker than that formed by nontreated cells or cells exposed to streptomycin. We have also observed that the treatment of X. fastidiosa with a sublethal concentration of gomesin before inoculation in tobacco plants correlates with a reduction in foliar symptoms, an effect possibly due to the trapping of bacterial cells to fewer xylem vessels, given the enhancement in biofilm production. These results warrant further investigation of how X. fastidiosa would respond to the AMPs produced by citrus endophytes and by the insect vector, leading to a better understanding of the mechanism of action of these molecules on bacterial virulence.

The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Gene expression profile of the plant pathogen Fusarium graminearum under the antagonistic effect of Pantoea agglomerans

The pathogenic fungus Fusarium graminearum is an ongoing threat to agriculture, causing losses in grain yield and quality in diverse crops. Substantial progress has been made in the identification of genes involved in the suppression of phytopathogens by antagonistic microorganisms; however, limited information regarding responses of plant pathogens to these biocontrol agents is available. Gene expression analysis was used to identify differentially expressed transcripts of the fungal plant pathogen F. graminearum under antagonistic effect of the bacterium Pantoea agglomerans. A macroarray was constructed, using 1014 transcripts from an F. graminearum cDNA library. Probes consisted of the cDNA of F. graminearum grown in the presence and in the absence of P. agglomerans. Twenty-nine genes were either up (19) or down (10) regulated during interaction with the antagonist bacterium. Genes encoding proteins associated with fungal defense and/or virulence or with nutritional and oxidative stress responses were induced. The repressed genes coded for a zinc finger protein associated with cell division, proteins containing cellular signaling domains, respiratory chain proteins, and chaperone-type proteins. These data give molecular and biochemical evidence of response of F. graminearum to an antagonist and could help develop effective biocontrol procedures for pathogenic plant fungi.

Global gene expression under nitrogen starvation in Xylella fastidiosa: contribution of the sigma(54) regulon

Background: Xylella fastidiosa, a Gram-negative fastidious bacterium, grows in the xylem of several plants causing diseases such as citrus variegated chlorosis. As the xylem sap contains low concentrations of amino acids and other compounds, X. fastidiosa needs to cope with nitrogen limitation in its natural habitat. Results: In this work, we performed a whole-genome microarray analysis of the X. fastidiosa nitrogen starvation response. A time course experiment (2, 8 and 12 hours) of cultures grown in defined medium under nitrogen starvation revealed many differentially expressed genes, such as those related to transport, nitrogen assimilation, amino acid biosynthesis, transcriptional regulation, and many genes encoding hypothetical proteins. In addition, a decrease in the expression levels of many genes involved in carbon metabolism and energy generation pathways was also observed. Comparison of gene expression profiles between the wild type strain and the rpoN null mutant allowed the identification of genes directly or indirectly induced by nitrogen starvation in a sigma(54)-dependent manner. A more complete picture of the sigma(54) regulon was achieved by combining the transcriptome data with an in silico search for potential sigma(54)-dependent promoters, using a position weight matrix approach. One of these sigma(54)-predicted binding sites, located upstream of the glnA gene (encoding glutamine synthetase), was validated by primer extension assays, confirming that this gene has a sigma(54)-dependent promoter. Conclusions: Together, these results show that nitrogen starvation causes intense changes in the X. fastidiosa transcriptome and some of these differentially expressed genes belong to the sigma(54) regulon.

Differential gene expression between the biotrophic-like and saprotrophic mycelia of the witches` broom pathogen Moniliophthora perniciosa

Moniliophthora perniciosa is a hemibiotrophic fungus that causes witches` broom disease (WBD) in cacao. Marked dimorphism characterizes this fungus, showing a monokaryotic or biotrophic phase that causes disease symptoms and a later dikaryotic or saprotrophic phase. A combined strategy of DNA microarray, expressed sequence tag, and real-time reverse-transcriptase polymerase chain reaction analyses was employed to analyze differences between these two fungal stages in vitro. In all, 1,131 putative genes were hybridized with cDNA from different phases, resulting in 189 differentially expressed genes, and 4,595 reads were clusterized, producing 1,534 unigenes. The analysis of these genes, which represent approximately 21% of the total genes, indicates that the biotrophic-like phase undergoes carbon and nitrogen catabollite repression that correlates to the expression of phytopathogenicity genes. Moreover, downregulation of mitochondrial oxidative phosphorylation and the presence of a putative ngr1 of Saccharomyces cerevisiae could help explain its lower growth rate. In contrast, the saprotrophic mycelium expresses genes related to the metabolism of hexoses, ammonia, and oxidative phosphorylation, which could explain its faster growth. Antifungal toxins were upregulated and could prevent the colonization by competing fungi. This work significantly contributes to our understanding of the molecular mechanisms of WBD and, to our knowledge, is the first to analyze differential gene expression of the different phases of a hemibiotrophic fungus.

Reciprocal changes in gene expression profiles of cocultured breast epithelial cells and primary fibroblasts

The importance of epithelial-stroma interaction in normal breast development and tumor progression has been recognized. To identify genes that were regulated by these reciprocal interactions, we cocultured a nonmalignant (MCF10A) and a breast cancer derived (MDA-MB231) basal cell lines, with fibroblasts isolated from breast benign-disease adjacent tissues (NAF) or with carcinoma-associated fibroblasts (CAF), in a transwell system. Gene expression profiles of each coculture pair were compared with the correspondent monocultures, using a customized microarray. Contrariwise to large alterations in epithelial cells genomic profiles, fibroblasts were less affected. In MDA-MB231 highly represented genes downregulated by CAF derived factors coded for proteins important for the specificity of vectorial transport between ER and golgi, possibly affecting cell polarity whereas the response of MCF10A comprised an induction of genes coding for stress responsive proteins, representing a prosurvival effect. While NAF downregulated genes encoding proteins associated to glycolipid and fatty acid biosynthesis in MDA-MB231, potentially affecting membrane biogenesis, in MCF10A, genes critical for growth control and adhesion were altered. NAFs responded to coculture with MDA-MB231 by a decrease in the expression of genes induced by TGF beta 1 and associated to motility. However, there was little change in NAFs gene expression profile influenced by MCF10A. CAFs responded to the presence of both epithelial cells inducing genes implicated in cell proliferation. Our data indicate that interactions between breast fibroblasts and basal epithelial cells resulted in alterations in the genomic profiles of both cell types which may help to clarify some aspects of this heterotypic signaling. (C) 2009 UICC

Insights on PRAME and osteosarcoma by means of gene expression profiling

Osteosarcoma (OS) is the most frequent bone tumor in children and adolescents. Tumor antigens are encoded by genes that are expressed in many types of solid tumors but are silent in normal tissues, with the exception of placenta and male germ-line cells. It has been proposed that antigen tumors are potential tumor markers. The premise of this study is that the identification of novel OS-associated transcripts will lead to a better understanding of the events involved in OS pathogenesis and biology. We analyzed the expression of a panel of seven tumor antigens in OS samples to identify possible tumor markers. After selecting the tumor antigen expressed in most samples of the panel, gene expression profiling was used to identify osteosarcoma-associated molecular alterations. A microarray was employed because of its ability to accurately produce comprehensive expression profiles. PRAME was identified as the tumor antigen expressed in most OS samples; it was detected in 68% of the cases. Microarray results showed differences in expression for genes functioning in cell signaling and adhesion as well as extracellular matrix-related genes, implying that such tumors could indeed differ in regard to distinct patterns of tumorigenesis. The hypothesis inferred in this study was gathered mostly from available data concerning other kinds of tumors. There is circumstantial evidence that PRAME expression might be related to distinct patterns of tumorigenesis. Further investigation is needed to validate the differential expression of genes belonging to tumorigenesis-related pathways in PRAME-positive and PRAME-negative tumors.

Alterations in myocardial gene expression associated with experimental Trypanosoma cruzi infection

Chagas disease, characterized by acute myocarditis and chronic cardiomyopathy, is caused by infection with the protozoan parasite Trypanosoma cruzi. We sought to identify genes altered during the development of parasite-induced cardiomyopathy. Microarrays containing 27,400 sequence-verified mouse cDNAs were used to analyze global gene expression changes in the myocardium of a murine model of chagasic cardiomyopathy. Changes in gene expression were determined as the acute stage of infection developed into the chronic stage. This analysis was performed on the hearts of male CD-1 mice infected with trypomastigotes of T. cruzi (Brazil strain). At each interval we compared infected and uninfected mice and confirmed the microarray data with dye reversal. We identified eight distinct categories of mRNAs that were differentially regulated during infection and identified dysregulation of several key genes. These data may provide insight into the pathogenesis of chagasic cardiomyopathy and provide new targets for intervention. (c) 2008 Elsevier Inc. All rights reserved.