Systems biology via redescription and ontologies (I): finding phase changes with applications to malaria temporal data

Biological systems are complex and often composed of many subtly interacting components. Furthermore, such systems evolve through time and, as the underlying biology executes its genetic program, the relationships between components change and undergo dynamic reorganization. Characterizing these relationships precisely is a challenging task, but one that must be undertaken if we are to understand these systems in sufficient detail. One set of tools that may prove useful are the formal principles of model building and checking, which could allow the biologist to frame these inherently temporal questions in a sufficiently rigorous framework. In response to these challenges, GOALIE (Gene ontology algorithmic logic and information extractor) was developed and has been successfully employed in the analysis of high throughput biological data (e.g. time-course gene-expression microarray data and neural spike train recordings). The method has applications to a wide variety of temporal data, indeed any data for which there exist ontological descriptions. This paper describes the algorithms behind GOALIE and its use in the study of the Intraerythrocytic Developmental Cycle (IDC) of Plasmodium falciparum, the parasite responsible for a deadly form of chloroquine resistant malaria. We focus in particular on the problem of finding phase changes, times of reorganization of transcriptional control

Functional Categories Associated with Clusters of Genes That Are Co-Expressed across the NCI-60 Cancer Cell Lines

The NCI-60 is a panel of 60 diverse human cancer cell lines used by the U.S. National Cancer Institute to screen compounds for anticancer activity. In the current study, gene expression levels from five platforms were integrated to yield a single composite transcriptome profile. The comprehensive and reliable nature of that dataset allows us to study gene co-expression across cancer cell lines.Hierarchical clustering revealed numerous clusters of genes in which the genes co-vary across the NCI-60. To determine functional categorization associated with each cluster, we used the Gene Ontology (GO) Consortium database and the GoMiner tool. GO maps genes to hierarchically-organized biological process categories. GoMiner can leverage GO to perform ontological analyses of gene expression studies, generating a list of significant functional categories.GoMiner analysis revealed many clusters of coregulated genes that are associated with functional groupings of GO biological process categories. Notably, those categories arising from coherent co-expression groupings reflect cancer-related themes such as adhesion, cell migration, RNA splicing, immune response and signal transduction. Thus, these clusters demonstrate transcriptional coregulation of functionally-related genes

RedundancyMiner: De-replication of redundant GO categories in microarray and proteomics analysis

<p>Abstract</p> <p>Background</p> <p>The Gene Ontology (GO) Consortium organizes genes into hierarchical categories based on biological process, molecular function and subcellular localization. Tools such as GoMiner can leverage GO to perform ontological analysis of microarray and proteomics studies, typically generating a list of significant functional categories. Two or more of the categories are often redundant, in the sense that identical or nearly-identical sets of genes map to the categories. The redundancy might typically inflate the report of significant categories by a factor of three-fold, create an illusion of an overly long list of significant categories, and obscure the relevant biological interpretation.</p> <p>Results</p> <p>We now introduce a new resource, RedundancyMiner, that de-replicates the redundant and nearly-redundant GO categories that had been determined by first running GoMiner. The main algorithm of RedundancyMiner, MultiClust, performs a novel form of cluster analysis in which a GO category might belong to several category clusters. Each category cluster follows a "complete linkage" paradigm. The metric is a similarity measure that captures the overlap in gene mapping between pairs of categories.</p> <p>Conclusions</p> <p>RedundancyMiner effectively eliminated redundancies from a set of GO categories. For illustration, we have applied it to the clarification of the results arising from two current studies: (1) assessment of the gene expression profiles obtained by laser capture microdissection (LCM) of serial cryosections of the retina at the site of final optic fissure closure in the mouse embryos at specific embryonic stages, and (2) analysis of a conceptual data set obtained by examining a list of genes deemed to be "kinetochore" genes.</p

Gitools: Analysis and Visualisation of Genomic Data Using Interactive Heat-Maps

Intuitive visualization of data and results is very important in genomics, especially when many conditions are to be analyzed and compared. Heat-maps have proven very useful for the representation of biological data. Here we present Gitools (http://www.gitools.org), an open-source tool to perform analyses and visualize data and results as interactive heat-maps. Gitools contains data import systems from several sources (i.e. IntOGen, Biomart, KEGG, Gene Ontology), which facilitate the integration of novel data with previous knowledge

Gene expression changes associated with Barrett's esophagus and Barrett's-associated adenocarcinoma cell lines after acid or bile salt exposure

<p>Abstract</p> <p>Background</p> <p>Esophageal reflux and Barrett's esophagus represent two major risk factors for the development of esophageal adenocarcinoma. Previous studies have shown that brief exposure of the Barrett's-associated adenocarcinoma cell line, SEG-1, or primary cultures of Barrett's esophageal tissues to acid or bile results in changes consistent with cell proliferation. In this study, we determined whether similar exposure to acid or bile salts results in gene expression changes that provide insights into malignant transformation.</p> <p>Methods</p> <p>Using previously published methods, Barrett's-associated esophageal adenocarcinoma cell lines and primary cultures of Barrett's esophageal tissue were exposed to short pulses of acid or bile salts followed by incubation in culture media at pH 7.4. A genome-wide assessment of gene expression was then determined for the samples using cDNA microarrays. Subsequent analysis evaluated for statistical differences in gene expression with and without treatment.</p> <p>Results</p> <p>The SEG-1 cell line showed changes in gene expression that was dependent on the length of exposure to pH 3.5. Further analysis using the Gene Ontology, however, showed that representation by genes associated with cell proliferation is not enhanced by acid exposure. The changes in gene expression also did not involve genes known to be differentially expressed in esophageal adenocarcinoma. Similar experiments using short-term primary cultures of Barrett's esophagus also did not result in detectable changes in gene expression with either acid or bile salt exposure.</p> <p>Conclusion</p> <p>Short-term exposure of esophageal adenocarcinoma SEG-1 cells or primary cultures of Barrett's esophagus does not result in gene expression changes that are consistent with enhanced cell proliferation. Thus other model systems are needed that may reflect the impact of acid and bile salt exposure on the esophagus <it>in vivo</it>.</p

Glucocorticoids with different chemical structures but similar glucocorticoid receptor potency regulate subsets of common and unique genes in human trabecular meshwork cells

<p>Abstract</p> <p>Background</p> <p>In addition to their well-documented ocular therapeutic effects, glucocorticoids (GCs) can cause sight-threatening side-effects including ocular hypertension presumably via morphological and biochemical changes in trabecular meshwork (TM) cells. In the present study, we directly compared the glucocorticoid receptor (GR) potency for dexamethasone (DEX), fluocinolone acetonide (FA) and triamcinolone acetonide (TA), examined the expression of known GRα and GRβ isoforms, and used gene expression microarrays to compare the effects of DEX, FA, and TA on the complete transcriptome in two primary human TM cell lines.</p> <p>Methods</p> <p>GR binding affinity for DEX, FA, and TA was measured by a cell-free competitive radio-labeled GR binding assay. GR-mediated transcriptional activity was assessed using the GeneBLAzer beta-lactamase reporter gene assay. Levels of GRα and GRβ isoforms were assessed by Western blot. Total RNA was extracted from TM 86 and TM 93 cells treated with 1 μM DEX, FA, or TA for 24 hr and used for microarray gene expression analysis. The microarray experiments were repeated three times. Differentially expressed genes were identified by Rosetta Resolver Gene Expression Analysis System.</p> <p>Results</p> <p>The GR binding affinity (IC₅₀) for DEX, FA, and TA was 5.4, 2.0, and 1.5 nM, respectively. These values are similar to the GR transactivation EC₅₀of 3.0, 0.7, and 1.5 nM for DEX, FA, and TA, respectively. All four GRα translational isoforms (A-D) were expressed in TM 86 and TM 93 total cell lysates, however, the C and D isoforms were more highly expressed relative to A and B. All four GRβ isoforms (A-D) were also detected in TM cells, although GRβ-D isoform expression was lower compared to that of the A, B, or C isoforms. Microarray analysis revealed 1,968 and 1,150 genes commonly regulated by DEX, FA, and TA in TM 86 and TM 93, respectively. These genes included RGC32, OCA2, ANGPTL7, MYOC, FKBP5, SAA1 and ZBTB16. In addition, each GC specifically regulated a unique set of genes in both TM cell lines. Using Ingenuity Pathway Analysis (IPA) software, analysis of the data from TM 86 cells showed that DEX significantly regulated transcripts associated with RNA post-transcriptional modifications, whereas FA and TA modulated genes involved in lipid metabolism and cell morphology, respectively. In TM 93 cells, DEX significantly regulated genes implicated in histone methylation, whereas FA and TA altered genes associated with cell cycle and cell adhesion, respectively.</p> <p>Conclusion</p> <p>Human trabecular meshwork cells in culture express all known GRα and GRβ translational isoforms, and GCs with similar potency but subtly different chemical structure are capable of regulating common and unique gene subsets and presumably biologic responses in these cells. These GC structure-dependent effects appear to be TM cell-lineage dependent.</p

Modulation of gene expression in heart and liver of hibernating black bears (Ursus americanus)

<p>Abstract</p> <p>Background</p> <p>Hibernation is an adaptive strategy to survive in highly seasonal or unpredictable environments. The molecular and genetic basis of hibernation physiology in mammals has only recently been studied using large scale genomic approaches. We analyzed gene expression in the American black bear, <it>Ursus americanus</it>, using a custom 12,800 cDNA probe microarray to detect differences in expression that occur in heart and liver during winter hibernation in comparison to summer active animals.</p> <p>Results</p> <p>We identified 245 genes in heart and 319 genes in liver that were differentially expressed between winter and summer. The expression of 24 genes was significantly elevated during hibernation in both heart and liver. These genes are mostly involved in lipid catabolism and protein biosynthesis and include RNA binding protein motif 3 (<it>Rbm3</it>), which enhances protein synthesis at mildly hypothermic temperatures. Elevated expression of protein biosynthesis genes suggests induction of translation that may be related to adaptive mechanisms reducing cardiac and muscle atrophies over extended periods of low metabolism and immobility during hibernation in bears. Coordinated reduction of transcription of genes involved in amino acid catabolism suggests redirection of amino acids from catabolic pathways to protein biosynthesis. We identify common for black bears and small mammalian hibernators transcriptional changes in the liver that include induction of genes responsible for fatty acid β oxidation and carbohydrate synthesis and depression of genes involved in lipid biosynthesis, carbohydrate catabolism, cellular respiration and detoxification pathways.</p> <p>Conclusions</p> <p>Our findings show that modulation of gene expression during winter hibernation represents molecular mechanism of adaptation to extreme environments.</p

Analysis of Gene Regulatory Networks in the Mammalian Circadian Rhythm

Circadian rhythm is fundamental in regulating a wide range of cellular, metabolic, physiological, and behavioral activities in mammals. Although a small number of key circadian genes have been identified through extensive molecular and genetic studies in the past, the existence of other key circadian genes and how they drive the genomewide circadian oscillation of gene expression in different tissues still remains unknown. Here we try to address these questions by integrating all available circadian microarray data in mammals. We identified 41 common circadian genes that showed circadian oscillation in a wide range of mouse tissues with a remarkable consistency of circadian phases across tissues. Comparisons across mouse, rat, rhesus macaque, and human showed that the circadian phases of known key circadian genes were delayed for 4–5 hours in rat compared to mouse and 8–12 hours in macaque and human compared to mouse. A systematic gene regulatory network for the mouse circadian rhythm was constructed after incorporating promoter analysis and transcription factor knockout or mutant microarray data. We observed the significant association of cis-regulatory elements: EBOX, DBOX, RRE, and HSE with the different phases of circadian oscillating genes. The analysis of the network structure revealed the paths through which light, food, and heat can entrain the circadian clock and identified that NR3C1 and FKBP/HSP90 complexes are central to the control of circadian genes through diverse environmental signals. Our study improves our understanding of the structure, design principle, and evolution of gene regulatory networks involved in the mammalian circadian rhythm