Integrating methylome and transcriptome signatures expands the molecular classification of the pituitary tumors

PURPOSE

To explore pituitary tumors by methylome and transcriptome signatures in a heterogeneous ethnic population.

DESIGN

Retrospective cross-sectional study.

PATIENTS AND METHODS

Clinicopathological features, methylome and transcriptome were evaluated in pituitary tumors from 77 patients (61% women, age: 12-72 years) followed due to functioning (FPT: GH-secreting n = 18, ACTH-secreting n = 14) and non-functioning pituitary tumors (NFPT, n = 45) at Ribeirao Preto Medical School, University of Sao Paulo.

RESULTS

Unsupervised hierarchical clustering analysis (UHCA) of methylome (n = 77) and transcriptome (n = 65 out of 77) revealed three clusters each: one enriched by FPT, other by NFPT, and another by ACTH-secreting and NFPT. Comparison between each omics-derived clusters identified 3,568 and 5,994 differentially methylated and expressed genes, respectively, which were associated with each other, with tumor clinical presentation, and with 2017 and 2022 WHO classifications. UHCA considering 11 transcripts related to pituitary development/differentiation also supported three clusters: POU1F1-driven somatotroph, TBX19-driven corticotroph, and NR5A1-driven gonadotroph adenomas, with rare exceptions (NR5A1 expressed in few GH-secreting and corticotroph-silent adenomas; POU1F1 in few ACTH-secreting adenomas; and TBX19 in few NFPTs).

CONCLUSIONS

This large heterogenic ethnic Brazilian cohort confirms that integrated methylome and transcriptome signatures classify FPT and NFPT, which are associated with clinical presentation and tumor invasiveness. Moreover, the cluster NFPT/ACTH-secreting adenomas raises interest regarding tumor heterogeneity, supporting the challenge raised by the 2017 and 2022 WHO definition regarding the discrepancy, in rare cases, between clinical presentation and pituitary lineage markers. Finally, making our data publicly available enables further studies to validate genes/pathways involved in pituitary tumor pathogenesis and prognosis.

Environmental Influences Measured by Epigenetic Clock and Vulnerability Components at Birth Impact Clinical ASD Heterogeneity

Although Autism Spectrum Disorders (ASD) is recognized as being heavily influenced by genetic factors, the role of epigenetic and environmental factors is still being established. This study aimed to identify ASD vulnerability components based on familial history and intrauterine environmental stress exposure, explore possible vulnerability subgroups, access DNA methylation age acceleration (AA) as a proxy of stress exposure during life, and evaluate the association of ASD vulnerability components and AA to phenotypic severity measures. Principal Component Analysis (PCA) was used to search the vulnerability components from 67 mothers of autistic children. We found that PC1 had a higher correlation with psychosocial stress (maternal stress, maternal education, and social class), and PC2 had a higher correlation with biological factors (psychiatric family history and gestational complications). Comparing the methylome between above and below PC1 average subgroups we found 11,879 statistically significant differentially methylated probes (DMPs, p < 0.05). DMPs CpG sites were enriched in variably methylated regions (VMRs), most showing environmental and genetic influences. Hypermethylated probes presented higher rates in different regulatory regions associated with functional SNPs, indicating that the subgroups may have different affected regulatory regions and their liability to disease explained by common variations. Vulnerability components score moderated by epigenetic clock AA was associated with Vineland Total score (p = 0.0036, adjR² = 0.31), suggesting risk factors with stress burden can influence ASD phenotype.

Halobacterium salinarum and Haloferax volcanii Comparative Transcriptomics Reveals Conserved Transcriptional Processing Sites

Post-transcriptional processing of messenger RNA is an important regulatory strategy that allows relatively fast responses to changes in environmental conditions. In halophile systems biology, the protein perspective of this problem (i.e., ribonucleases which implement the cleavages) is generally more studied than the RNA perspective (i.e., processing sites). In the present in silico work, we mapped genome-wide transcriptional processing sites (TPS) in two halophilic model organisms, Halobacterium salinarum NRC-1 and Haloferax volcanii DS2. TPS were established by reanalysis of publicly available differential RNA-seq (dRNA-seq) data, searching for non-primary (monophosphorylated RNAs) enrichment. We found 2093 TPS in 43% of H. salinarum genes and 3515 TPS in 49% of H. volcanii chromosomal genes. Of the 244 conserved TPS sites found, the majority were located around start and stop codons of orthologous genes. Specific genes are highlighted when discussing antisense, ribosome and insertion sequence associated TPS. Examples include the cell division gene ftsZ2, whose differential processing signal along growth was detected and correlated with post-transcriptional regulation, and biogenesis of sense overlapping transcripts associated with IS200/IS605. We hereby present the comparative, transcriptomics-based processing site maps with a companion browsing interface.

SAT-LB36 Distinct Vitamin D Receptor DNA Methylation Profiles Are Associated With the Outcome of Pediatric Patients With Adrenocortical Tumors

Pediatric adrenocortical tumors (pACT) are rare, display complex genomic background and lack robust prognostic markers. Very recently, distinct genomic methylation profiles of pACT were associated with prognosis. The vitamin D receptor (VDR) was shown to be underexpressed in ACT, especially in carcinomas (ACC). In adult ACC, VDR inactivation by methylation was demonstrated. On the other hand, VDR activation was shown to inhibit ACC proliferation in vitro and in vivo.

Aim: To evaluate VDR DNA methylation profile and its clinical and prognostic significance in pediatric ACT.

Methods: Genomic DNA methylation from 57 pACTs [40 girls; median age: 2.1 (0.2-16.4) years] was assessed using Infininium Methylation EPIC BeadChip Array. Unsupervised hierarchical clustering analysis (Ward method, R Stats Package) was performed considering the M-values of the 49 probes targeting the whole extension of VDR gene contained in the array. Clinical, histopathological and molecular features, as well as pACT VDR mRNA levels (qPCR) and nuclear immunoreactivity (IHC) were used for association analysis.

Results: Hierarchical clustering identified three clusters of pACT. Methylated VDR-targeted probes (M-values different from 0; n=37) composed the VDR methylation profile, which differed significantly between the clusters [M-values: C1=1.77 (1.1-1.9) (low), C2=2.15 (1.7-2.7) (intermediate), and C3=2.65 (2.2-3.1) (high); p<0.0001]. The C1 cluster comprised a set of patients with favorable outcome (n=18), who were younger (p=0.035), did not present metastasis at diagnosis (IPACTR stage IV) or after surgery, nor were diagnosed with carcinomas (Wieneke criteria >=4), were not carriers of somatic Beta-catenin activating mutations, or died. Although cluster C2 patients (n=21) presented intermediary disease features, only 2 patients died and the overall outcome was positive. Instead, the C3 cluster concentrated patients (n=18) with non-localized/metastatic disease (IPACTR stages I/II vs. III/IV; p=0.004), post-surgical metastasis/recurrence (p=0.009), and patients who needed adjuvant chemotherapy (p=0.005). Moreover, C3 patients had lower overall and disease-free survival rates (log-rank: p=0.001 and p=0.014, respectively). VDR methylation was not associated with sex, clinical presentation, P53 mutations, nor with tumor VDR mRNA expression or nuclear immunoreactivity.

Conclusions: Three VDR methylation profiles were associated with distinct pACT clinical features and outcome. High VDR methylation was associated with worst outcome. Fully functioning VDR may play a beneficial role against pediatric adrenocortical tumorigenesis. This finding highlights the potential of targeting VDR as an adjuvant therapeutic target.

Lineage relationship between prostate adenocarcinoma and small cell carcinoma

BACKGROUND

Prostate cancer displays different morphologies which, in turn, affect patient outcome. This fact prompted questions about the lineage relationship between differentiated, more treatable prostate adenocarcinoma and poorly differentiated, less treatable non-adenocarcinoma including small cell carcinoma, and the molecular mechanism underlying prostate cancer differentiation.

METHODS

Newly available non-adenocarcinoma/small cell carcinoma PDX LuCaP lines were analyzed for expression of stem cell transcription factors (scTF) LIN28A, NANOG, POU5F1, SOX2, which are responsible for reprogramming or de-differentiation. cDNA of these genes were cloned from small cell carcinoma LuCaP 145.1 into expression vectors to determine if they could function in reprogramming.

RESULTS

Expression of scTF was detected in small cell carcinoma LuCaP 93, 145.1, 145.2, and non-adenocarcinoma LuCaP 173.1, 173.2A. Transfection of scTF from LuCaP 145.1 altered the gene expression of prostate non-small cell carcinoma cells, as well as fibroblasts. The resultant cells grew in stem-like colonies. Of note was a 10-fold lower expression of B2M in the transfected cells. Low B2M was also characteristic of LuCaP 145.1. Conversely, B2M was increased when stem cells were induced to differentiate.

CONCLUSIONS

This work suggested a pathway in the emergence of non-adenocarcinoma/small cell carcinoma from adenocarcinoma through activation of scTF genes that produced cancer de-differentiation.

Internal RNAs overlapping coding sequences can drive the production of alternative proteins in archaea

Prokaryotic genomes show a high level of information compaction often with different molecules transcribed from the same locus. Although antisense RNAs have been relatively well studied, RNAs in the same strand, internal RNAs (intraRNAs), are still poorly understood. The question of how common is the translation of overlapping reading frames remains open. We address this question in the model archaeon Halobacterium salinarum. In the present work we used differential RNA-seq (dRNA-seq) in H. salinarum NRC-1 to locate intraRNA signals in subsets of internal transcription start sites (iTSS) and establish the open reading frames associated to them (intraORFs). Using C-terminally flagged proteins, we experimentally observed isoforms accurately predicted by intraRNA translation for kef1, acs3 and orc4 genes. We also recovered from the literature and mass spectrometry databases several instances of protein isoforms consistent with intraRNA translation such as the gas vesicle protein gene gvpC1. We found evidence for intraRNAs in horizontally transferred genes such as the chaperone dnaK and the aerobic respiration related cydA in both H. salinarum and Escherichia coli. Also, intraRNA translation evidence in H. salinarum, E. coli and yeast of a universal elongation factor (aEF-2, fusA and eEF-2) suggests that this is an ancient phenomenon present in all domains of life.

Conversion of Prostate Adenocarcinoma to Small Cell Carcinoma-Like by Reprogramming

The lineage relationship between prostate adenocarcinoma and small cell carcinoma was studied by using the LuCaP family of xenografts established from primary neoplasm to metastasis. Expression of four stem cell transcription factor (TF) genes, LIN28A, NANOG, POU5F1, SOX2, were analyzed in the LuCaP lines. These genes, when force expressed in differentiated cells, can reprogram the recipients into stem-like induced pluripotent stem (iPS) cells. Most LuCaP lines expressed POU5F1, while LuCaP 145.1, representative of small cell carcinoma, expressed all four. Through transcriptome database query, many small cell carcinoma genes were also found in stem cells. To test the hypothesis that prostate cancer progression from "differentiated" adenocarcinoma to "undifferentiated" small cell carcinoma could involve re-expression of stem cell genes, the four TF genes were transduced via lentiviral vectors into five adenocarcinoma LuCaP lines-70CR, 73CR, 86.2, 92, 105CR-as done in iPS cell reprogramming. The resultant cells from these five transductions displayed a morphology of small size and dark appearing unlike the parentals. Transcriptome analysis of LuCaP 70CR* ("*" to denote transfected progeny) revealed a unique gene expression close to that of LuCaP 145.1. In a prostate principal components analysis space based on cell-type transcriptomes, the different LuCaP transcriptome datapoints were aligned to suggest a possible ordered sequence of expression changes from the differentiated luminal-like adenocarcinoma cell types to the less differentiated, more stem-like small cell carcinoma types, and LuCaP 70CR*. Prostate cancer progression can thus be molecularly characterized by loss of differentiation with re-expression of stem cell genes. J. Cell. Physiol. 231: 2040-2047, 2016. © 2016 Wiley Periodicals, Inc.

Sense overlapping transcripts in IS1341-type transposase genes are functional non-coding RNAs in archaea

The existence of sense overlapping transcripts that share regulatory and coding information in the same genomic sequence shows an additional level of prokaryotic gene expression complexity. Here we report the discovery of ncRNAs associated with IS1341-type transposase (tnpB) genes, at the 3'-end of such elements, with examples in archaea and bacteria. Focusing on the model haloarchaeon Halobacterium salinarum NRC-1, we show the existence of sense overlapping transcripts (sotRNAs) for all its IS1341-type transposases. Publicly available transcriptome compendium show condition-dependent differential regulation between sotRNAs and their cognate genes. These sotRNAs allowed us to find a UUCA tetraloop motif that is present in other archaea (ncRNA family HgcC) and in a H. salinarum intergenic ncRNA derived from a palindrome associated transposable elements (PATE). Overexpression of one sotRNA and the PATE-derived RNA harboring the tetraloop motif improved H. salinarum growth, indicating that these ncRNAs are functional.

Processing of voided urine for prostate cancer RNA biomarker analysis

BACKGROUND

Voided urine samples have been shown to contain cells released from prostate tumors. Could good quality RNA from cells in urine be obtained from every donor for multimarker analysis? In addition, could urine donation be as simple as possible, a practical consideration for a lab test, without involving a prostate massage (as indicated for PCA3 testing), which precludes frequent collection; needing it done at a specific time of day (e.g., first or second urine); and requiring prompt processing of samples in clinics with limited molecular biology capability?

METHODS

Collected urine samples were pelleted, and the RNA isolated was processed for cDNA synthesis and in vitro transcription to generate amplified sense aRNA. The resultant aRNA was rigorously analyzed for possible introduced changes. DMSO was used as a cell preservative for frozen storage of urine samples.

RESULTS

Good quality aRNA was obtained for over 100 samples collected at two different institutions. The process of RNA amplification removed co-isolated DNA in some samples, which did not affect RNA amplification. Amplification did not amplify genes that were absent and produce other expression alterations. The sense aRNA could be used to generate urinary transcriptomes specific to individual patients. No chaotropic agents for RNA preservation were added to the urine samples so that the supernatant could be used for analysis of secreted protein biomarkers. The time of donation was not important since patients were seen during the entire day. DMSO was an effective cell preservative for freezing urine.

CONCLUSIONS

Urinary RNA can be readily isolated and amplified for prostate cancer biomarker analysis. Individual patients had unique set of transcripts derived from their tumor.

A Methodology for the Development of RESTful Semantic Web Services for Gene Expression Analysis

Gene expression studies are generally performed through multi-step analysis processes, which require the integrated use of a number of analysis tools. In order to facilitate tool/data integration, an increasing number of analysis tools have been developed as or adapted to semantic web services. In recent years, some approaches have been defined for the development and semantic annotation of web services created from legacy software tools, but these approaches still present many limitations. In addition, to the best of our knowledge, no suitable approach has been defined for the functional genomics domain. Therefore, this paper aims at defining an integrated methodology for the implementation of RESTful semantic web services created from gene expression analysis tools and the semantic annotation of such services. We have applied our methodology to the development of a number of services to support the analysis of different types of gene expression data, including microarray and RNASeq. All developed services are publicly available in the Gene Expression Analysis Services (GEAS) Repository at http://dcm.ffclrp.usp.br/lssb/geas. Additionally, we have used a number of the developed services to create different integrated analysis scenarios to reproduce parts of two gene expression studies documented in the literature. The first study involves the analysis of one-color microarray data obtained from multiple sclerosis patients and healthy donors. The second study comprises the analysis of RNA-Seq data obtained from melanoma cells to investigate the role of the remodeller BRG1 in the proliferation and morphology of these cells. Our methodology provides concrete guidelines and technical details in order to facilitate the systematic development of semantic web services. Moreover, it encourages the development and reuse of these services for the creation of semantically integrated solutions for gene expression analysis.

Reprogramming of prostate cancer cells--technical challenges

Prostate cancer progression is characterized by tumor dedifferentiation. Cancer cells of less differentiated tumors have a gene expression/transcriptome more similar to that of stem cells. In dedifferentiation, cancer cells may follow a specific program of gene expression changes to a stem-like state. In order to treat cancer effectively, the stem-like cancer cells and cancer differentiation pathway need to be identified and studied. Due to the very low abundance of stem-like cancer cells, their isolation from fresh human tumors is technically challenging. Induced pluripotent stem cell technology can reprogram differentiated cells into stem-like, and this may be a tool to generate sufficient stem-like cancer cells.

A metabolomic protocol for plant systematics by matrix-assisted laser-desorption/ionization time-of flight mass spectrometry

Matrix-assisted laser desorption/ionization time-of flight mass spectrometry (MALDI-TOF MS) has been widely used for the identification and classification of microorganisms based on their proteomic fingerprints. However, the use of MALDI-TOF MS in plant research has been very limited. In the present study, a first protocol is proposed for metabolic fingerprinting by MALDI-TOF MS using three different MALDI matrices with subsequent multivariate data analysis by in-house algorithms implemented in the R environment for the taxonomic classification of plants from different genera, families and orders. By merging the data acquired with different matrices, different ionization modes and using careful algorithms and parameter selection, we demonstrate that a close taxonomic classification can be achieved based on plant metabolic fingerprints, with 92% similarity to the taxonomic classifications found in literature. The present work therefore highlights the great potential of applying MALDI-TOF MS for the taxonomic classification of plants and, furthermore, provides a preliminary foundation for future research.

Transcription start site associated RNAs (TSSaRNAs) are ubiquitous in all domains of life

A plethora of non-coding RNAs has been discovered using high-resolution transcriptomics tools, indicating that transcriptional and post-transcriptional regulation is much more complex than previously appreciated. Small RNAs associated with transcription start sites of annotated coding regions (TSSaRNAs) are pervasive in both eukaryotes and bacteria. Here, we provide evidence for existence of TSSaRNAs in several archaeal transcriptomes including: Halobacterium salinarum, Pyrococcus furiosus, Methanococcus maripaludis, and Sulfolobus solfataricus. We validated TSSaRNAs from the model archaeon Halobacterium salinarum NRC-1 by deep sequencing two independent small-RNA enriched (RNA-seq) and a primary-transcript enriched (dRNA-seq) strand-specific libraries. We identified 652 transcripts, of which 179 were shown to be primary transcripts (∼7% of the annotated genome). Distinct growth-associated expression patterns between TSSaRNAs and their cognate genes were observed, indicating a possible role in environmental responses that may result from RNA polymerase with varying pausing rhythms. This work shows that TSSaRNAs are ubiquitous across all domains of life.

ProbMetab: an R package for Bayesian probabilistic annotation of LC-MS-based metabolomics

Summary: We present ProbMetab, an R package that promotes substantial improvement in automatic probabilistic liquid chromatography–mass spectrometry-based metabolome annotation. The inference engine core is based on a Bayesian model implemented to (i) allow diverse source of experimental data and metadata to be systematically incorporated into the model with alternative ways to calculate the likelihood function and (ii) allow sensitive selection of biologically meaningful biochemical reaction databases as Dirichlet-categorical prior distribution. Additionally, to ensure result interpretation by system biologists, we display the annotation in a network where observed mass peaks are connected if their candidate metabolites are substrate/product of known biochemical reactions. This graph can be overlaid with other graph-based analysis, such as partial correlation networks, in a visualization scheme exported to Cytoscape, with web and stand-alone versions.

Availability and implementation: ProbMetab was implemented in a modular manner to fit together with established upstream (xcms, CAMERA, AStream, mzMatch.R, etc) and downstream R package tools (GeneNet, RCytoscape, DiffCorr, etc). ProbMetab, along with extensive documentation and case studies, is freely available under GNU license at: http://labpib.fmrp.usp.br/methods/probmetab/.

Contact:rvencio@usp.br

Supplementary information:Supplementary data are available at Bioinformatics online.

Semantic integration of gene expression analysis tools and data sources using software connectors

BACKGROUND

The study and analysis of gene expression measurements is the primary focus of functional genomics. Once expression data is available, biologists are faced with the task of extracting (new) knowledge associated to the underlying biological phenomenon. Most often, in order to perform this task, biologists execute a number of analysis activities on the available gene expression dataset rather than a single analysis activity. The integration of heterogeneous tools and data sources to create an integrated analysis environment represents a challenging and error-prone task. Semantic integration enables the assignment of unambiguous meanings to data shared among different applications in an integrated environment, allowing the exchange of data in a semantically consistent and meaningful way. This work aims at developing an ontology-based methodology for the semantic integration of gene expression analysis tools and data sources. The proposed methodology relies on software connectors to support not only the access to heterogeneous data sources but also the definition of transformation rules on exchanged data.

RESULTS

We have studied the different challenges involved in the integration of computer systems and the role software connectors play in this task. We have also studied a number of gene expression technologies, analysis tools and related ontologies in order to devise basic integration scenarios and propose a reference ontology for the gene expression domain. Then, we have defined a number of activities and associated guidelines to prescribe how the development of connectors should be carried out. Finally, we have applied the proposed methodology in the construction of three different integration scenarios involving the use of different tools for the analysis of different types of gene expression data.

CONCLUSIONS

The proposed methodology facilitates the development of connectors capable of semantically integrating different gene expression analysis tools and data sources. The methodology can be used in the development of connectors supporting both simple and nontrivial processing requirements, thus assuring accurate data exchange and information interpretation from exchanged data.

A UML profile for the OBO relation ontology

Background

Ontologies have increasingly been used in the biomedical domain, which has prompted the emergence of different initiatives to facilitate their development and integration. The Open Biological and Biomedical Ontologies (OBO) Foundry consortium provides a repository of life-science ontologies, which are developed according to a set of shared principles. This consortium has developed an ontology called OBO Relation Ontology aiming at standardizing the different types of biological entity classes and associated relationships. Since ontologies are primarily intended to be used by humans, the use of graphical notations for ontology development facilitates the capture, comprehension and communication of knowledge between its users. However, OBO Foundry ontologies are captured and represented basically using text-based notations. The Unified Modeling Language (UML) provides a standard and widely-used graphical notation for modeling computer systems. UML provides a well-defined set of modeling elements, which can be extended using a built-in extension mechanism named Profile. Thus, this work aims at developing a UML profile for the OBO Relation Ontology to provide a domain-specific set of modeling elements that can be used to create standard UML-based ontologies in the biomedical domain.

Results

We have studied the OBO Relation Ontology, the UML metamodel and the UML profiling mechanism. Based on these studies, we have proposed an extension to the UML metamodel in conformance with the OBO Relation Ontology and we have defined a profile that implements the extended metamodel. Finally, we have applied the proposed UML profile in the development of a number of fragments from different ontologies. Particularly, we have considered the Gene Ontology (GO), the PRotein Ontology (PRO) and the Xenopus Anatomy and Development Ontology (XAO).

Conclusions

The use of an established and well-known graphical language in the development of biomedical ontologies provides a more intuitive form of capturing and representing knowledge than using only text-based notations. The use of the profile requires the domain expert to reason about the underlying semantics of the concepts and relationships being modeled, which helps preventing the introduction of inconsistencies in an ontology under development and facilitates the identification and correction of errors in an already defined ontology.

Reprogramming of prostate cancer-associated stromal cells to embryonic stem-like

BACKGROUND

CD90(+) prostate cancer-associated (CP) stromal cells represent a diseased cell type found only in tumor tissue. They differ from their normal counterpart in gene expression and inductive signaling. Genetic reprogramming by induced pluripotent stem (iPS) cell technology can effectively change adult cells into stem-like cells through wholesale alteration of the gene expression program. This technology might be used to 'erase' the abnormal gene expression of diseased cells. The resultant iPS cells would no longer express the disease phenotype, and behave like stem cells.

METHODS

CP stromal cells, isolated from tumor tissue of a surgically resected prostate by anti-CD90-mediated sorting and cultured in vitro, were transfected with in vitro packaged lentiviral expression vectors containing stem cell transcription factor genes POU5F1, LIN28, NANOG, and SOX2.

RESULTS

Alkaline phosphatase-positive iPS cells were obtained in about 3 weeks post-transfection at a frequency of 10(-4) . Their colony morphology was indistinguishable from that of human embryonic stem (ES) cells. Transcriptome analysis showed a virtually complete match in gene expression between the iPS and ES cells.

CONCLUSIONS

Genes of CP stromal cells could be fully inactivated by genetic reprogramming. As a consequence, the disease phenotype was 'cured'.

Bladder expression of CD cell surface antigens and cell-type-specific transcriptomes

Many cell types have no known functional attributes. In the bladder and prostate, basal epithelial and stromal cells appear similar in cytomorphology and share several cell surface markers. Their total gene expression (transcriptome) should provide a clear measure of the extent to which they are alike functionally. Since urologic stromal cells are known to mediate organ-specific tissue formation, these cells in cancers might exhibit aberrant gene expression affecting their function. For transcriptomes, cluster designation (CD) antigens have been identified for cell sorting. The sorted cell populations can be analyzed by DNA microarrays. Various bladder cell types have unique complements of CD molecules. CD9(+) urothelial, CD104(+) basal and CD13(+) stromal cells of the lamina propria were therefore analyzed, as were CD9(+) cancer and CD13(+) cancer-associated stromal cells. The transcriptome datasets were compared by principal components analysis for relatedness between cell types; those with similarity in gene expression indicated similar function. Although bladder and prostate basal cells shared CD markers such as CD104, CD44 and CD49f, they differed in overall gene expression. Basal cells also lacked stem cell gene expression. The bladder luminal and stromal transcriptomes were distinct from their prostate counterparts. In bladder cancer, not only the urothelial but also the stromal cells showed gene expression alteration. The cancer process in both might thus involve defective stromal signaling. These cell-type transcriptomes provide a means to monitor in vitro models in which various CD-isolated cell types can be combined to study bladder differentiation and bladder tumor development based on cell-cell interaction.

Differential Inductive Signaling of CD90 Prostate Cancer-Associated Fibroblasts Compared to Normal Tissue Stromal Mesenchyme Cells

Prostate carcinomas are surrounded by a layer of stromal fibroblastic cells that are characterized by increased expression of CD90. These CD90⁺ cancer-associated stromal fibroblastic cells differ in gene expression from their normal counterpart, CD49a⁺CD90^lo stromal smooth muscle cells; and were postulated to represent a less differentiated cell type with altered inductive properties. CD90⁺ stromal cells were isolated from tumor tissue specimens and co-cultured with the pluripotent embryonal carcinoma cell line NCCIT in order to elucidate the impact of tumor-associated stroma on stem cells, and the ‘cancer stem cell.’ Transcriptome analysis identified a notable decreased induction of smooth muscle and prostate stromal genes such as PENK, BMP2 and ChGn compared to previously determined NCCIT response to normal prostate stromal cell induction. CD90⁺ stromal cell secreted factors induced an increased expression of CD90 and differential induction of genes involved in extracellular matrix remodeling and the RECK pathway in NCCIT. These results suggest that, compared to normal tissue stromal cells, signaling from cancer-associated stromal cells has a markedly different effect on stem cells as represented by NCCIT. Given that stromal cells are important in directing organ-specific differentiation, stromal cells in tumors appear to be defective in this function, which may contribute to abnormal differentiation found in diseases such as cancer.

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 relationship between prostate cancer cells of Gleason 3, 4 and normal epithelial cells as revealed by cell type-specific transcriptomes

BACKGROUND

Prostate cancer cells in primary tumors have been typed CD10-/CD13-/CD24hi/CD26+/CD38lo/CD44-/CD104-. This CD phenotype suggests a lineage relationship between cancer cells and luminal cells. The Gleason grade of tumors is a descriptive of tumor glandular differentiation. Higher Gleason scores are associated with treatment failure.

METHODS

CD26+ cancer cells were isolated from Gleason 3+3 (G3) and Gleason 4+4 (G4) tumors by cell sorting, and their gene expression or transcriptome was determined by Affymetrix DNA array analysis. Dataset analysis was used to determine gene expression similarities and differences between G3 and G4 as well as to prostate cancer cell lines and histologically normal prostate luminal cells.

RESULTS

The G3 and G4 transcriptomes were compared to those of prostatic cell types of non-cancer, which included luminal, basal, stromal fibromuscular, and endothelial. A principal components analysis of the various transcriptome datasets indicated a closer relationship between luminal and G3 than luminal and G4. Dataset comparison also showed that the cancer transcriptomes differed substantially from those of prostate cancer cell lines.

CONCLUSIONS

Genes differentially expressed in cancer are potential biomarkers for cancer detection, and those differentially expressed between G3 and G4 are potential biomarkers for disease stratification given that G4 cancer is associated with poor outcomes. Differentially expressed genes likely contribute to the prostate cancer phenotype and constitute the signatures of these particular cancer cell types.

Temporal expression profiling of the effects of secreted factors from prostate stromal cells on embryonal carcinoma stem cells

BACKGROUND

There is a growing body of evidence indicating that epigenetic influences originating from stromal cells in the immediate microenvironment may play a role in carcinogenesis. Determining the molecular mechanisms involved in stromal-stem cell interaction could provide critical insight into prostate development and disease progression, particularly with regard to their relationship to and influence on the putative cancer stem cell.

METHODS

Prostate and bladder stromal cells prepared from tissue specimens were co-cultured with the pluripotent embryonal carcinoma cell line NCCIT. Transcriptome analysis was used to characterize NCCIT cell response to prostate or bladder signaling.

RESULTS

A systems approach demonstrated that prostate stromal cells were capable of inducing gene expression changes in NCCIT through secreted factors. Induction led to a loss of embryonic stem cell markers, with concurrent up-regulation of many genes characteristic of stromal mesenchyme cells as well as some of epithelial and cancer stem cells. Bladder stromal signaling produced gene expression changes different from those of prostate signaling.

CONCLUSIONS

This study indicates that paracrine stromal cell signaling can affect cancer stem cell response in an organ-specific manner and may provide insight for future development of treatment strategies such as differentiation therapy.

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.

Plasmodium vivax and the importance of the subtelomeric multigene vir superfamily

Plasmodium vivax is responsible for more than 100 million clinical cases yearly. Unlike P. falciparum, in which infected red blood cells cytoadhere via variant proteins, avoiding passage through the spleen, P.-vivax-infected reticulocytes seem not to cytoadhere. However, a variant subtelomeric multigene vir family has been identified in P. vivax. Thus, questions remain about how P. vivax circulates through the spleen and the role of Vir proteins. In this review, the importance of the vir multigene superfamily is reviewed in the light of the completion of the entire genome sequence of P. vivax and from data gathered from experimental infections in reticulocyte-prone non-lethal malaria parasites and natural P. vivax infections.

Multisite fMRI reproducibility of a motor task using identical MR systems

PURPOSE

To estimate between-scanner functional MRI (fMRI) reproducibility in a multisite study.

MATERIALS AND METHODS

A total of five identical 1.5T MR systems were used to repeatedly scan five subjects while performing a finger tapping task. A two-way (scanners, subjects) random effects analysis of variance (ANOVA) was used to estimate between-scanner and between-subject variability on two outcome variables: task-related mean blood oxygenation level dependent (BOLD) signal change and volume of activation within a motor region-of-interest (ROI).

RESULTS

Between-scanner variability of fMRI data accounted for a small proportion of the total variation in the BOLD signal change (8.34%, P = 0.114) and volume of activation (5.46%, P = 0.203). Between-subject variation accounted for more than half of the total variation for both measurements (57.17% and 54.46%, respectively, P < 0.01).

CONCLUSION

These results support the feasibility of multisite studies using identical scanner systems.

HTself: self-self based statistical test for low replication microarray studies

Different statistical methods have been used to classify a gene as differentially expressed in microarray experiments. They usually require a number of experimental observations to be adequately applied. However, many microarray experiments are constrained to low replication designs for different reasons, from financial restrictions to scarcely available RNA samples. Although performed in a high-throughput framework, there are few experimental replicas for each gene to allow the use of traditional or state-of-art statistical methods. In this work, we present a web-based bioinformatics tool that deals with real-life problems concerning low replication experiments. It uses an empirically derived criterion to classify a gene as differentially expressed by combining two widely accepted ideas in microarray analysis: self-self experiments to derive intensity-dependent cutoffs and non-parametric estimation techniques. To help laboratories without a bioinformatics infrastructure, we implemented the tool in a user-friendly website (http://blasto.iq.usp.br/~rvencio/HTself).

ProbCD: enrichment analysis accounting for categorization uncertainty

Background

As in many other areas of science, systems biology makes extensive use of statistical association and significance estimates in contingency tables, a type of categorical data analysis known in this field as enrichment (also over-representation or enhancement) analysis. In spite of efforts to create probabilistic annotations, especially in the Gene Ontology context, or to deal with uncertainty in high throughput-based datasets, current enrichment methods largely ignore this probabilistic information since they are mainly based on variants of the Fisher Exact Test.

Results

We developed an open-source R-based software to deal with probabilistic categorical data analysis, ProbCD, that does not require a static contingency table. The contingency table for the enrichment problem is built using the expectation of a Bernoulli Scheme stochastic process given the categorization probabilities. An on-line interface was created to allow usage by non-programmers and is available at: .

Conclusion

We present an analysis framework and software tools to address the issue of uncertainty in categorical data analysis. In particular, concerning the enrichment analysis, ProbCD can accommodate: (i) the stochastic nature of the high-throughput experimental techniques and (ii) probabilistic gene annotation.

Comparative transcriptome analysis of Listeria monocytogenes strains of the two major lineages reveals differences in virulence, cell wall, and stress response

Listeria monocytogenes is a food-borne, opportunistic, bacterial pathogen causing a wide spectrum of diseases, including meningitis, septicemia, abortion, and gastroenteritis, in humans and animals. Among the 13 L. monocytogenes serovars described, human listeriosis is mostly associated with strains of serovars 4b, 1/2b, and 1/2a. Within the species L. monocytogenes, three phylogenetic lineages are described. Serovar 1/2a belongs to phylogenetic lineage I, while serovars 4b and 1/2b group in phylogenetic lineage II. To explore the role of gene expression in the adaptation of L. monocytogenes strains of these two major lineages to different environments, as well as in virulence, we performed whole-genome expression profiling of six L. monocytogenes isolates of serovars 4b, 1/2b, and 1/2a of distinct origins, using a newly constructed Listeria multigenome DNA array. Comparison of the global gene expression profiles revealed differences among strains. The expression profiles of two strains having distinct 50% lethal doses, as assessed in the mouse model, were further analyzed. Gene ontology term enrichment analysis of the differentially expressed genes identified differences in protein-, nucleic acid-, carbon metabolism-, and virulence-related gene expression. Comparison of the expression profiles of the core genomes of all strains revealed differences between the two lineages with respect to cell wall synthesis, the stress-related sigma B regulon and virulence-related genes. These findings suggest different patterns of interaction with host cells and the environment, key factors for host colonization and survival in the environment.

Global gene expression analysis of the heat shock response in the phytopathogen Xylella fastidiosa

Xylella fastidiosa is a phytopathogenic bacterium that is responsible for diseases in many economically important crops. Although different strains have been studied, little is known about X. fastidiosa stress responses. One of the better characterized stress responses in bacteria is the heat shock response, which induces the expression of specific genes to prevent protein misfolding and aggregation and to promote degradation of the irreversibly denatured polypeptides. To investigate X. fastidiosa genes involved in the heat shock response, we performed a whole-genome microarray analysis in a time course experiment. Globally, 261 genes were induced (9.7%) and 222 genes were repressed (8.3%). The expression profiles of the differentially expressed genes were grouped, and their expression patterns were validated by quantitative reverse transcription-PCR experiments. We determined the transcription start sites of six heat shock-inducible genes and analyzed their promoter regions, which allowed us to propose a putative consensus for sigma(32) promoters in Xylella and to suggest additional genes as putative members of this regulon. Besides the induction of classical heat shock protein genes, we observed the up-regulation of virulence-associated genes such as vapD and of genes for hemagglutinins, hemolysin, and xylan-degrading enzymes, which may indicate the importance of heat stress to bacterial pathogenesis. In addition, we observed the repression of genes related to fimbriae, aerobic respiration, and protein biosynthesis and the induction of genes related to the extracytoplasmic stress response and some phage-related genes, revealing the complex network of genes that work together in response to heat shock.

Differential transcription profiles in Trypanosoma cruzi associated with clinical forms of Chagas disease: Maxicircle NADH dehydrogenase subunit 7 gene truncation in asymptomatic patient isolates

The majority of individuals in the chronic phase of Chagas disease are asymptomatic (indeterminate form). Every year 2-3% of these individuals develop severe clinical manifestations (cardiac and digestive forms). In this study a Trypanosoma cruzi DNA microarray was used to compare the transcript profiles of six human isolates: three from asymptomatic and three from cardiac patients. Seven signals were expressed differentially between the two classes of isolates, including tryparedoxin, surface protease GP63, cyclophilin, some hypothetical proteins and the pre-edited maxicircle gene NADH dehydrogenase subunit 7 (ND7). The approximately 30-fold greater signal in cardiac strains for ND7 was the most pronounced of the group, and differential levels of pre-edited ND7 transcript confirmed the microarray analysis. The ND7 gene from asymptomatic isolates showed a deletion of 455bp from nt 222 to nt 677 relative to ND7 of the CL Brener reference strain. The ND7 gene structure correlated with disease manifestation for 20 isolates from clinically characterised, chronic phase patients. The ND7 lesion produces a truncated product that could impair the function of mitochondrial complex I. Possible links between the integrity of the electron transport chain and symptom presentation are discussed. We propose that ND7 and other genes of the pathway constitute valuable targets for PCR assays in the differential diagnosis of the infective T. cruzi strain. While this hypothesis requires validation by the examination of additional recent parasite isolates from patients with defined pathologies, the identification of specific molecular markers represents a promising advance in the association between parasite genetics and disease pathology.

BayBoots: a model-free Bayesian tool to identify class markers from gene expression data

One of the goals of gene expression experiments is the identification of differentially expressed genes among populations that could be used as markers. For this purpose, we implemented a model-free Bayesian approach in a user-friendly and freely available web-based tool called BayBoots. In spite of a common misunderstanding that Bayesian and model-free approaches are incompatible, we merged them in the BayBoots implementation using the Kernel density estimator and Rubin 's Bayesian Bootstrap. We used the Bayes error rate (BER) instead of the usual P values as an alternative statistical index to rank a class marker's discriminative potential, since it can be visualized by a simple graphical representation and has an intuitive interpretation. Subsequently, Bayesian Bootstrap was used to assess BER 's credibility. We tested BayBoots on microarray data to look for markers for Trypanosoma cruzi strains isolated from cardiac and asymptomatic patients. We found that the three most frequently used methods in microarray analysis: t-test, non-parametric Wilcoxon test and correlation methods, yielded several markers that were discarded by a time-consuming visual check. On the other hand, the BayBoots graphical output and ranking was able to automatically identify markers for which classification performance was consistent. BayBoots is available at: http://www.vision.ime.usp.br/~rvencio/BayBoots.

SpotWhatR: a user-friendly microarray data analysis system

SpotWhatR is a user-friendly microarray data analysis tool that runs under a widely and freely available R statistical language (http://www.r-project.org) for Windows and Linux operational systems. The aim of SpotWhatR is to help the researcher to analyze microarray data by providing basic tools for data visualization, normalization, determination of differentially expressed genes, summarization by Gene Ontology terms, and clustering analysis. SpotWhatR allows researchers who are not familiar with computational programming to choose the most suitable analysis for their microarray dataset. Along with well-known procedures used in microarray data analysis, we have introduced a stand-alone implementation of the HTself method, especially designed to find differentially expressed genes in low-replication contexts. This approach is more compatible with our local reality than the usual statistical methods. We provide several examples derived from the Blastocladiella emersonii and Xylella fastidiosa Microarray Projects. SpotWhatR is freely available at http://blasto.iq.usp.br/~tkoide/SpotWhatR, in English and Portuguese versions. In addition, the user can choose between "single experiment" and "batch processing" versions.

Gene discovery and expression profile analysis through sequencing of expressed sequence tags from different developmental stages of the chytridiomycete Blastocladiella emersonii

Blastocladiella emersonii is an aquatic fungus of the chytridiomycete class which diverged early from the fungal lineage and is notable for the morphogenetic processes which occur during its life cycle. Its particular taxonomic position makes this fungus an interesting system to be considered when investigating phylogenetic relationships and studying the biology of lower fungi. To contribute to the understanding of the complexity of the B. emersonii genome, we present here a survey of expressed sequence tags (ESTs) from various stages of the fungal development. Nearly 20,000 cDNA clones from 10 different libraries were partially sequenced from their 5' end, yielding 16,984 high-quality ESTs. These ESTs were assembled into 4,873 putative transcripts, of which 48% presented no matches with existing sequences in public databases. As a result of Gene Ontology (GO) project annotation, 1,680 ESTs (35%) were classified into biological processes of the GO structure, with transcription and RNA processing, protein biosynthesis, and transport as prevalent processes. We also report full-length sequences, useful for construction of molecular phylogenies, and several ESTs that showed high similarity with known proteins, some of which were not previously described in fungi. Furthermore, we analyzed the expression profile (digital Northern analysis) of each transcript throughout the life cycle of the fungus using Bayesian statistics. The in silico approach was validated by Northern blot analysis with good agreement between the two methodologies.

Whole-genome expression profiling of Xylella fastidiosa in response to growth on glucose

Xylella fastidiosa is the etiologic agent of diseases in a wide range of economically important crops including citrus variegated chlorosis, a major threat to the Brazilian citrus industry. The genomes of several strains of this phytopathogen have been completely sequenced enabling large-scale functional studies. In this work we used whole-genome DNA microarrays to investigate the transcription profile of X. fastidiosa grown in defined media with different glucose concentrations. Our analysis revealed that while transcripts related to fastidian gum production were unaffected, colicin-V-like and fimbria precursors were induced in high glucose medium. Based on these results, we suggest a model for colicin-defense mechanism in X. fastidiosa.

DNA microarrays for comparative genomics and analysis of gene expression in Trypanosoma cruzi

Trypanosoma cruzi presents high genetic diversity and parasite isolates show remarkable differences in biological parameters. In this study, we evaluated whether DNA microarrays containing CL Brener cDNAs can be used for comparative genomics and for the analysis of gene expression in T. cruzi. We constructed a prototype microarray with 710 expression sequence tags of CL Brener and 20 sequences of T. cruzi strains. These probes represent 665 unique genes. Results from four hybridisations with genomic DNA of Silvio (T. cruzi I) and CL Brener (hybrid genotype) identified 9.3% of the probes (68/730) differentially represented in the two genomes. Data from eight hybridisations with cDNA obtained from three independent parasite harvests of Silvio and CL Brener disclosed 84 sequences of 730 (11.5%) that showed statistical significant (P < or = 0.01) changes in expression (1.6-6.5-fold). Some of the array-identified sequences were confirmed by Southern and Northern blot analysis. Only 20% of the probes with increased expression in Silvio or CL Brener presented higher hybridisation with genomic DNA of either strain. Approximately 2.5% (18/730) and 9.0% (65/730) of the probes were differentially expressed (P < or = 0.01), respectively, in epimastigotes and metacyclic trypomastigotes of two T. cruzi II strains isolated from chronic chagasic patients. Microarrays identified several sequences for which differences in gene copy number and/or in the levels of RNA transcripts were previously demonstrated by different approaches. The data indicate that DNA microarrays are a useful tool for comparative studies between strains and provide further evidence for a high level of post-transcriptional regulation of RNA abundance in T. cruzi.

Using credibility intervals instead of hypothesis tests in SAGE analysis

MOTIVATION

Statistical methods usually used to perform Serial Analysis of Gene Expression (SAGE) analysis are based on hypothesis testing. They answer the biologist's question: 'what are the genes with differential expression greater than r with P-value smaller than P?'. Another useful and not yet explored question is: 'what is the uncertainty in differential expression ratio of a gene?'.

RESULTS

We have used Bayesian model for SAGE differential gene expression ratios as a more informative alternative to hypothesis tests since it provides credibility intervals.

AVAILABILITY

The model is implemented in R statistical language script and is available under GNU/GLP copyleft at supplemental web site.

SUPPLEMENTARY INFORMATION

http://www.ime.usp.br/~rvencio/SAGEci/

DNA microarray-based genome comparison of a pathogenic and a nonpathogenic strain of Xylella fastidiosa delineates genes important for bacterial virulence

Xylella fastidiosa is a phytopathogenic bacterium that causes serious diseases in a wide range of economically important crops. Despite extensive comparative analyses of genome sequences of Xylella pathogenic strains from different plant hosts, nonpathogenic strains have not been studied. In this report, we show that X. fastidiosa strain J1a12, associated with citrus variegated chlorosis (CVC), is nonpathogenic when injected into citrus and tobacco plants. Furthermore, a DNA microarray-based comparison of J1a12 with 9a5c, a CVC strain that is highly pathogenic and had its genome completely sequenced, revealed that 14 coding sequences of strain 9a5c are absent or highly divergent in strain J1a12. Among them, we found an arginase and a fimbrial adhesin precursor of type III pilus, which were confirmed to be absent in the nonpathogenic strain by PCR and DNA sequencing. The absence of arginase can be correlated to the inability of J1a12 to multiply in host plants. This enzyme has been recently shown to act as a bacterial survival mechanism by down-regulating host nitric oxide production. The lack of the adhesin precursor gene is in accordance with the less aggregated phenotype observed for J1a12 cells growing in vitro. Thus, the absence of both genes can be associated with the failure of the J1a12 strain to establish and spread in citrus and tobacco plants. These results provide the first detailed comparison between a nonpathogenic strain and a pathogenic strain of X. fastidiosa, constituting an important step towards understanding the molecular basis of the disease.