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.

Searching for molecular markers in head and neck squamous cell carcinomas (HNSCC) by statistical and bioinformatic analysis of larynx-derived SAGE libraries

Background: Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancies in humans. The average 5-year survival rate is one of the lowest among aggressive cancers, showing no significant improvement in recent years. When detected early, HNSCC has a good prognosis, but most patients present metastatic disease at the time of diagnosis, which significantly reduces survival rate. Despite extensive research, no molecular markers are currently available for diagnostic or prognostic purposes. Methods: Aiming to identify differentially-expressed genes involved in laryngeal squamous cell carcinoma (LSCC) development and progression, we generated individual Serial Analysis of Gene Expression (SAGE) libraries from a metastatic and non-metastatic larynx carcinoma, as well as from a normal larynx mucosa sample. Approximately 54,000 unique tags were sequenced in three libraries. Results: Statistical data analysis identified a subset of 1,216 differentially expressed tags between tumor and normal libraries, and 894 differentially expressed tags between metastatic and non-metastatic carcinomas. Three genes displaying differential regulation, one down-regulated (KRT31) and two up-regulated (BST2, MFAP2), as well as one with a non-significant differential expression pattern (GNA15) in our SAGE data were selected for real-time polymerase chain reaction (PCR) in a set of HNSCC samples. Consistent with our statistical analysis, quantitative PCR confirmed the upregulation of BST2 and MFAP2 and the downregulation of KRT31 when samples of HNSCC were compared to tumor-free surgical margins. As expected, GNA15 presented a non-significant differential expression pattern when tumor samples were compared to normal tissues. Conclusion: To the best of our knowledge, this is the first study reporting SAGE data in head and neck squamous cell tumors. Statistical analysis was effective in identifying differentially expressed genes reportedly involved in cancer development. The differential expression of a subset of genes was confirmed in additional larynx carcinoma samples and in carcinomas from a distinct head and neck subsite. This result suggests the existence of potential common biomarkers for prognosis and targeted-therapy development in this heterogeneous type of tumor.

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.

Chicken skeletal muscle-associated macroarray for gene discovery

Macro- and microarrays are well-established technologies to determine gene functions through repeated measurements of transcript abundance. We constructed a chicken skeletal muscle-associated array based on a muscle-specific EST database, which was used to generate a tissue expression dataset of similar to 4500 chicken genes across 5 adult tissues (skeletal muscle, heart, liver, brain, and skin). Only a small number of ESTs were sufficiently well characterized by BLAST searches to determine their probable cellular functions. Evidence of a particular tissue-characteristic expression can be considered an indication that the transcript is likely to be functionally significant. The skeletal muscle macroarray platform was first used to search for evidence of tissue-specific expression, focusing on the biological function of genes/transcripts, since gene expression profiles generated across tissues were found to be reliable and consistent. Hierarchical clustering analysis revealed consistent clustering among genes assigned to 'developmental growth', such as the ontology genes and germ layers. Accuracy of the expression data was supported by comparing information from known transcripts and tissue from which the transcript was derived with macroarray data. Hybridization assays resulted in consistent tissue expression profile, which will be useful to dissect tissue-regulatory networks and to predict functions of novel genes identified after extensive sequencing of the genomes of model organisms. Screening our skeletal-muscle platform using 5 chicken adult tissues allowed us identifying 43 'tissue-specific' transcripts, and 112 co-expressed uncharacterized transcripts with 62 putative motifs. This platform also represents an important tool for functional investigation of novel genes; to determine expression pattern according to developmental stages; to evaluate differences in muscular growth potential between chicken lines, and to identify tissue-specific genes.

Origins of the Xylella fastidiosa Prophage-Like Regions and Their Impact in Genome Differentiation

Xylella fastidiosa is a Gram negative plant pathogen causing many economically important diseases, and analyses of completely sequenced X. fastidiosa genome strains allowed the identification of many prophage-like elements and possibly phage remnants, accounting for up to 15% of the genome composition. To better evaluate the recent evolution of the X. fastidiosa chromosome backbone among distinct pathovars, the number and location of prophage-like regions on two finished genomes (9a5c and Temecula1), and in two candidate molecules (Ann1 and Dixon) were assessed. Based on comparative best bidirectional hit analyses, the majority (51%) of the predicted genes in the X. fastidiosa prophage-like regions are related to structural phage genes belonging to the Siphoviridae family. Electron micrograph reveals the existence of putative viral particles with similar morphology to lambda phages in the bacterial cell in planta. Moreover, analysis of microarray data indicates that 9a5c strain cultivated under stress conditions presents enhanced expression of phage anti-repressor genes, suggesting switches from lysogenic to lytic cycle of phages under stress-induced situations. Furthermore, virulence-associated proteins and toxins are found within these prophage-like elements, thus suggesting an important role in host adaptation. Finally, clustering analyses of phage integrase genes based on multiple alignment patterns reveal they group in five lineages, all possessing a tyrosine recombinase catalytic domain, and phylogenetically close to other integrases found in phages that are genetic mosaics and able to perform generalized and specialized transduction. Integration sites and tRNA association is also evidenced. In summary, we present comparative and experimental evidence supporting the association and contribution of phage activity on the differentiation of Xylella genomes.

Electrochemical Evidence of Strong Electronic Interaction of PtRu on Carbon Nanotubes with High Density of Defects

We show that carbon nanotubes (CNTs) with high density of defects can present a strong electronic interaction with nanoparticles of Pt-Ru with average particle size of 3.5 +/- 0.8 nm. Depending on the Pt-Ru loading on the CNTs, CO and methanol oxidation reactions suggest there is a charge transfer between Pt-Ru that in turn provokes a decrease in the electronic interaction taking place between Ru and Pt in the PtRu alloy. The CO stripping potentials were observed at about 0.65 and 0.5 V for Pt-Ru/CNT electrodes with Pt-Ru loadings of 10 and 20, and 30 wt %, respectively. (C) 2008 The Electrochemical Society. [DOI: 10.1149/1.2990222] All rights reserved.

Composite electrode of carbon nanotubes and vitreous carbon for electron field emission

In this work, the electron field emission behaviour of electrodes formed by carbon nanotubes (CNTs) grown onto monolithic vitreous carbon (VCarbon) substrates with microcavities is presented. Scanning electron microscopy was used to characterize the microstructure of the films. Tungsten probes, stainless steel sphere, and phosphor electrodes were employed in the electron field emission study. The CNT/VCarbon composite represents a route to inexpensive excellent large area electron emission cathodes with fields as low as 2.1 V mu m(-1). In preliminary lifetime tests for a period of about 24 h at an emission current of about 4 mA cm(-2), there is an onset degradation of the emission current of about 28%, which then stabilizes. Electron emission images of the composites show the cavity of the samples act as separate emission sites and predominantly control the emission process. The emission of CNTs/VCarbon was found to be stable for several hours. (c) 2008 American Institute of Physics.

Sugarcane genes associated with sucrose content

Background -: Sucrose content is a highly desirable trait in sugarcane as the worldwide demand for cost-effective biofuels surges. Sugarcane cultivars differ in their capacity to accumulate sucrose and breeding programs routinely perform crosses to identify genotypes able to produce more sucrose. Sucrose content in the mature internodes reach around 20% of the culms dry weight. Genotypes in the populations reflect their genetic program and may display contrasting growth, development, and physiology, all of which affect carbohydrate metabolism. Few studies have profiled gene expression related to sugarcane's sugar content. The identification of signal transduction components and transcription factors that might regulate sugar accumulation is highly desirable if we are to improve this characteristic of sugarcane plants. Results -: We have evaluated thirty genotypes that have different Brix (sugar) levels and identified genes differentially expressed in internodes using cDNA microarrays. These genes were compared to existing gene expression data for sugarcane plants subjected to diverse stress and hormone treatments. The comparisons revealed a strong overlap between the drought and sucrose-content datasets and a limited overlap with ABA signaling. Genes associated with sucrose content were extensively validated by qRT-PCR, which highlighted several protein kinases and transcription factors that are likely to be regulators of sucrose accumulation. The data also indicate that aquaporins, as well as lignin biosynthesis and cell wall metabolism genes, are strongly related to sucrose accumulation. Moreover, sucrose-associated genes were shown to be directly responsive to short term sucrose stimuli, confirming their role in sugar-related pathways. Conclusion -: Gene expression analysis of sugarcane populations contrasting for sucrose content indicated a possible overlap with drought and cell wall metabolism processes and suggested signaling and transcriptional regulators to be used as molecular markers in breeding programs. Transgenic research is necessary to further clarify the role of the genes and define targets useful for sugarcane improvement programs based on transgenic plants.

ProbFAST: Probabilistic Functional Analysis System Tool

Background: The post-genomic era has brought new challenges regarding the understanding of the organization and function of the human genome. Many of these challenges are centered on the meaning of differential gene regulation under distinct biological conditions and can be performed by analyzing the Multiple Differential Expression (MDE) of genes associated with normal and abnormal biological processes. Currently MDE analyses are limited to usual methods of differential expression initially designed for paired analysis. Results: We proposed a web platform named ProbFAST for MDE analysis which uses Bayesian inference to identify key genes that are intuitively prioritized by means of probabilities. A simulated study revealed that our method gives a better performance when compared to other approaches and when applied to public expression data, we demonstrated its flexibility to obtain relevant genes biologically associated with normal and abnormal biological processes. Conclusions: ProbFAST is a free accessible web-based application that enables MDE analysis on a global scale. It offers an efficient methodological approach for MDE analysis of a set of genes that are turned on and off related to functional information during the evolution of a tumor or tissue differentiation. ProbFAST server can be accessed at http://gdm.fmrp.usp.br/probfast.

Gene expression down-regulation in CD90(+) prostate tumor-associated stromal cells involves potential organ-specific genes

Background: The prostate stroma is a key mediator of epithelial differentiation and development, and potentially plays a role in the initiation and progression of prostate cancer. The tumor-associated stroma is marked by increased expression of CD90/THYI. Isolation and characterization of these stromal cells could provide valuable insight into the biology of the tumor microenvironment. Methods: Prostate CD90(+) stromal fibromuscular cells from tumor specimens were isolated by cell-sorting and analyzed by DNA microarray. Dataset analysis was used to compare gene expression between histologically normal and tumor-associated stromal cells. For comparison, stromal cells were also isolated and analyzed from the urinary bladder. Results: The tumor-associated stromal cells were found to have decreased expression of genes involved in smooth muscle differentiation, and those detected in prostate but not bladder. Other differential expression between the stromal cell types included that of the CXC-chemokine genes. Conclusion: CD90(+) prostate tumor-associated stromal cells differed from their normal counterpart in expression of multiple genes, some of which are potentially involved in organ development.

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.

Gene Expression Profiles in Radiation Workers Occupationally Exposed to Ionizing Radiation

Ionizing radiation OR) imposes risks to human health and the environment. IR at low doses and low (lose rates has the potency to initiate carcinogenesis. Genotoxic environmental agents such as IR trigger a cascade of signal transduction pathways for cellular protection. In this study, using cDNA microarray technique, we monitored the gene expression profiles in lymphocytes derived from radiation-ex posed individuals (radiation workers). Physical dosimetry records on these patients indicated that the absorbed dose ranged from 0.696 to 39.088 mSv. Gene expression analysis revealed statistically significant transcriptional changes in a total of 78 genes (21 up-regulated and 57 clown-regulated) involved in several biological processes such as ubiquitin cycle (UHRF2 and PIAS1), DNA repair (LIG3, XPA, ERCC5, RAD52, DCLRE1C), cell cycle regulation/proliferation (RHOA, CABLES2, TGFB2, IL16), and stress response (GSTP1, PPP2R5A, DUSP22). Some of the genes that showed altered expression profiles in this study call be used as biomarkers for monitoring the chronic low level exposure in humans. Additionally, alterations in gene expression patterns observed in chronically exposed radiation workers reinforces the need for defining the effective radiation dose that causes immediate genetic damage as well as the long-term effects on genomic instability, including cancer.

A role for mammalian target of rapamycin (mTOR) pathway in non alcoholic steatohepatitis related-cirrhosis

Non-alcoholic fatty liver disease (NAFLD) encompasses the whole spectrum of steatosis, nonalcoholic steatohepatitis (NASH), and NASH-related cirrhosis (NASH/Cir). Although molecular advances have been made in this field, the pathogenesis of NAFLD is not completely understood. The gene expression profiling associated to NASH/Cir was assessed, in an attempt to better characterize the pathways involved in its etiopathogenesis. Methods: In the first step, we used cDNA microarray to evaluate the gene expression profiles in normal liver (n=3) and NASH/Cir samples (n=3) by GeneSifter (TM) analysis to identify differentially expressed genes and biological pathways. Second, tissue microarray was used to determine immunohistochemical expression of phosphorylated mTOR and 4E-BP1 in 11 normal liver samples, 10 NASH/Cir samples and in 37 samples of cirrhosis of other etiologies to further explore the involvement of the mTOR pathway evidenced by the gene expression analysis. Results: 138 and 106 genes were, respectively, up and down regulated in NASH/Cir in comparison to normal liver. Among the 9 pathways identified as significantly modulated in NASH/Cir, the participation of the mTOR pathway was confirmed, since expression of cytoplasmic and membrane phospho-mTOR were higher in NASH/Cir in comparison to cirrhosis of other etiologies and to normal liver. Conclusions: Recent findings have suggested a role for the cellular ""nutrient sensor"" mTOR in NAFLD and the present study corroborates the participation of this pathway in NASH/Cir. Phospho-mTOR evaluation might be of clinical utility as a potential marker for identification of NASH/Cir in cases mistakenly considered as cryptogenic cirrhosis owing to paucity of clinical data.

Endometrial claudin-4 and leukemia inhibitory factor are associated with assisted reproduction outcome

Background: Claudin-4 (CLDN4) is one of several proteins that act as molecular mediators of embryo implantation. Recently, we examined immunolabeling of leukemia inhibitory factor (LIF) in the endometrial tissue of 52 IVF patients, and found that LIF staining intensity was strongly correlated with successful pregnancy initiation. In the same set of patients, we have now examined endometrial CLDN4 expression, to see how expression intensity may vary with LIF. We examined CLDN4 in the luteal phase of the menstrual cycle, immediately preceding IVF treatment. Our aim was to compare expression of LIF and CLDN4 in the luteal phase, and document these patterns as putative biomarkers for pregnancy. Methods: Endometrial tissue was collected from women undergoing IVF. Endometrial biopsies were obtained during the luteal phase preceding IVF, and were then used for tissue microarray (TMA) immunolabeling of CLDN4. Previously published LIF expression data were then combined with CLDN4 expression data, to determine CLDN4/LIF expression patterns. Associations between successful pregnancy after IVF and combined CLDN4/LIF expression patterns were evaluated. Results: Four patterns of immunolabeling were observed in the endometrial samples: 16% showed weak CLDN4 and strong LIF (CLDN4(-)/LIF(+)); 20% showed strong CLDN4 and strong LIF (LIF(+)/CLDN4(+)); 28% showed strong CLDN4 and weak LIF (CLDN4(+)/LIF(-)); and 36% showed weak CLDN4 and weak LIF (CLDN4(-)/LIF(-)). Successful implantation after IVF was associated with CLDN4(-)/LIF(+)(p = 0.003). Patients showing this endometrial CLDN4(-)/LIF(+) immunolabeling were also 6 times more likely to achieve pregnancy than patients with endometrial CLDN4(+)/LIF(-) immunolabeling (p = 0.007). Conclusion: The combined immunolabeling expression of CLDN4(-)/LIF(+) in endometrial tissue is a potential biomarker for predicting successful pregnancy in IVF candidates.

Proteome analysis of schizophrenia patients Wernicke's area reveals an energy metabolism dysregulation

Background: Schizophrenia is likely to be a consequence of DNA alterations that, together with environmental factors, will lead to protein expression differences and the ultimate establishment of the illness. The superior temporal gyrus is implicated in schizophrenia and executes functions such as the processing of speech, language skills and sound processing. Methods: We performed an individual comparative proteome analysis using two-dimensional gel electrophoresis of 9 schizophrenia and 6 healthy control patients' left posterior superior temporal gyrus (Wernicke's area - BA22p) identifying by mass spectrometry several protein expression alterations that could be related to the disease. Results: Our analysis revealed 11 downregulated and 14 upregulated proteins, most of them related to energy metabolism. Whereas many of the identified proteins have been previously implicated in schizophrenia, such as fructose-bisphosphate aldolase C, creatine kinase and neuron-specific enolase, new putative disease markers were also identified such as dihydrolipoyl dehydrogenase, tropomyosin 3, breast cancer metastasis-suppressor 1, heterogeneous nuclear ribonucleoproteins C1/C2 and phosphate carrier protein, mitochondrial precursor. Besides, the differential expression of peroxiredoxin 6 (PRDX6) and glial fibrillary acidic protein (GFAP) were confirmed by western blot in schizophrenia prefrontal cortex. Conclusion: Our data supports a dysregulation of energy metabolism in schizophrenia as well as suggests new markers that may contribute to a better understanding of this complex disease.