Functional assessment of time course microarray data

FUNNEL-GSEA: FUNctioNal ELastic-net regression in time-course gene set enrichment analysis

Motivation

Gene set enrichment analyses (GSEAs) are widely used in genomic research to identify underlying biological mechanisms (defined by the gene sets), such as Gene Ontology terms and molecular pathways. There are two caveats in the currently available methods: (i) they are typically designed for group comparisons or regression analyses, which do not utilize temporal information efficiently in time-series of transcriptomics measurements; and (ii) genes overlapping in multiple molecular pathways are considered multiple times in hypothesis testing.

Results

We propose an inferential framework for GSEA based on functional data analysis, which utilizes the temporal information based on functional principal component analysis, and disentangles the effects of overlapping genes by a functional extension of the elastic-net regression. Furthermore, the hypothesis testing for the gene sets is performed by an extension of Mann-Whitney U test which is based on weighted rank sums computed from correlated observations. By using both simulated datasets and a large-scale time-course gene expression data on human influenza infection, we demonstrate that our method has uniformly better receiver operating characteristic curves, and identifies more pathways relevant to immune-response to human influenza infection than the competing approaches.

Availability and Implementation

The methods are implemented in R package FUNNEL, freely and publicly available at: https://github.com/yunzhang813/FUNNEL-GSEA-R-Package.

Supplementary information

Supplementary data are available at Bioinformatics online.

Neurospora crassa transcriptomics reveals oxidative stress and plasma membrane homeostasis biology genes as key targets in response to chitosan

Chitosan is a natural polymer with antimicrobial activity. Chitosan causes plasma membrane permeabilization and induction of intracellular reactive oxygen species (ROS) in Neurospora crassa. We have determined the transcriptional profile of N. crassa to chitosan and identified the main gene targets involved in the cellular response to this compound. Global network analyses showed membrane, transport and oxidoreductase activity as key nodes affected by chitosan. Activation of oxidative metabolism indicates the importance of ROS and cell energy together with plasma membrane homeostasis in N. crassa response to chitosan. Deletion strain analysis of chitosan susceptibility pointed NCU03639 encoding a class 3 lipase, involved in plasma membrane repair by lipid replacement, and NCU04537 a MFS monosaccharide transporter related to assimilation of simple sugars, as main gene targets of chitosan. NCU10521, a glutathione S-transferase-4 involved in the generation of reducing power for scavenging intracellular ROS is also a determinant chitosan gene target. Ca(2+) increased tolerance to chitosan in N. crassa. Growth of NCU10610 (fig 1 domain) and SYT1 (a synaptotagmin) deletion strains was significantly increased by Ca(2+) in the presence of chitosan. Both genes play a determinant role in N. crassa membrane homeostasis. Our results are of paramount importance for developing chitosan as an antifungal.

Time-Course Gene Set Analysis for Longitudinal Gene Expression Data

Gene set analysis methods, which consider predefined groups of genes in the analysis of genomic data, have been successfully applied for analyzing gene expression data in cross-sectional studies. The time-course gene set analysis (TcGSA) introduced here is an extension of gene set analysis to longitudinal data. The proposed method relies on random effects modeling with maximum likelihood estimates. It allows to use all available repeated measurements while dealing with unbalanced data due to missing at random (MAR) measurements. TcGSA is a hypothesis driven method that identifies a priori defined gene sets with significant expression variations over time, taking into account the potential heterogeneity of expression within gene sets. When biological conditions are compared, the method indicates if the time patterns of gene sets significantly differ according to these conditions. The interest of the method is illustrated by its application to two real life datasets: an HIV therapeutic vaccine trial (DALIA-1 trial), and data from a recent study on influenza and pneumococcal vaccines. In the DALIA-1 trial TcGSA revealed a significant change in gene expression over time within 69 gene sets during vaccination, while a standard univariate individual gene analysis corrected for multiple testing as well as a standard a Gene Set Enrichment Analysis (GSEA) for time series both failed to detect any significant pattern change over time. When applied to the second illustrative data set, TcGSA allowed the identification of 4 gene sets finally found to be linked with the influenza vaccine too although they were found to be associated to the pneumococcal vaccine only in previous analyses. In our simulation study TcGSA exhibits good statistical properties, and an increased power compared to other approaches for analyzing time-course expression patterns of gene sets. The method is made available for the community through an R package.

Gene set analysis methods use prior biological knowledge to analyze gene expression data. This prior knowledge takes the form of predefined groups of genes, linked through their biological function. Gene set analysis methods have been successfully applied in transversal studies, their results being more sensitive and interpretable than those of methods investigating genomic data one gene at a time. The time-course gene set analysis (TcGSA) introduced here is an extension of such gene set analysis to longitudinal data. This method identifies a priori defined groups of genes whose expression is not stable over time, taking into account the potential heterogeneity between patients and between genes. When biological conditions are compared, it identifies the gene sets that have different expression dynamics according to these conditions. Data from 2 studies are analyzed: data from an HIV therapeutic vaccine trial, and data from a recent study on influenza and pneumococcal vaccines. In both cases, TcGSA provided new insights compared to standard approaches thanks to an increased sensitivity compared to other approaches. Those results highlight the benefits of the TcGSA method for analyzing gene expression dynamics.

Next maSigPro: updating maSigPro bioconductor package for RNA-seq time series

Motivation: The widespread adoption of RNA-seq to quantitatively measure gene expression has increased the scope of sequencing experimental designs to include time-course experiments. maSigPro is an R package specifically suited for the analysis of time-course gene expression data, which was developed originally for microarrays and hence was limited in its application to count data.

Results: We have updated maSigPro to support RNA-seq time series analysis by introducing generalized linear models in the algorithm to support the modeling of count data while maintaining the traditional functionalities of the package. We show a good performance of the maSigPro-GLM method in several simulated time-course scenarios and in a real experimental dataset.

Availability and implementation: The package is freely available under the LGPL license from the Bioconductor Web site (http://bioconductor.org).

Contact:mj.nueda@ua.es or aconesa@cipf.es

timeClip: pathway analysis for time course data without replicates

Background

Time-course gene expression experiments are useful tools for exploring biological processes. In this type of experiments, gene expression changes are monitored along time. Unfortunately, replication of time series is still costly and usually long time course do not have replicates. Many approaches have been proposed to deal with this data structure, but none of them in the field of pathway analysis. Pathway analyses have acquired great relevance for helping the interpretation of gene expression data. Several methods have been proposed to this aim: from the classical enrichment to the more complex topological analysis that gains power from the topology of the pathway. None of them were devised to identify temporal variations in time course data.

Results

Here we present timeClip, a topology based pathway analysis specifically tailored to long time series without replicates. timeClip combines dimension reduction techniques and graph decomposition theory to explore and identify the portion of pathways that is most time-dependent. In the first step, timeClip selects the time-dependent pathways; in the second step, the most time dependent portions of these pathways are highlighted. We used timeClip on simulated data and on a benchmark dataset regarding mouse muscle regeneration model. Our approach shows good performance on different simulated settings. On the real dataset, we identify 76 time-dependent pathways, most of which known to be involved in the regeneration process. Focusing on the 'mTOR signaling pathway' we highlight the timing of key processes of the muscle regeneration: from the early pathway activation through growth factor signals to the late burst of protein production needed for the fiber regeneration.

Conclusions

timeClip represents a new improvement in the field of time-dependent pathway analysis. It allows to isolate and dissect pathways characterized by time-dependent components. Furthermore, using timeClip on a mouse muscle regeneration dataset we were able to characterize the process of muscle fiber regeneration with its correct timing.

The common ground of genomics and systems biology

The rise of systems biology is intertwined with that of genomics, yet their primordial relationship to one another is ill-defined. We discuss how the growth of genomics provided a critical boost to the popularity of systems biology. We describe the parts of genomics that share common areas of interest with systems biology today in the areas of gene expression, network inference, chromatin state analysis, pathway analysis, personalized medicine, and upcoming areas of synergy as genomics continues to expand its scope across all biomedical fields.

Transcriptomic analysis highlights epigenetic and transcriptional regulation during zygotic embryo development of Pinus pinaster

BACKGROUND

It is during embryogenesis that the plant body plan is established and the meristems responsible for all post-embryonic growth are specified. The molecular mechanisms governing conifer embryogenesis are still largely unknown. Their elucidation may contribute valuable information to clarify if the distinct features of embryo development in angiosperms and gymnosperms result from differential gene regulation. To address this issue, we have performed the first transcriptomic analysis of zygotic embryo development in a conifer species (Pinus pinaster) focusing our study in particular on regulatory genes playing important roles during plant embryo development, namely epigenetic regulators and transcription factors.

RESULTS

Microarray analysis of P. pinaster zygotic embryogenesis was performed at five periods of embryo development from early developing to mature embryos. Our results show that most changes in transcript levels occurred in the first and the last embryo stage-to-stage transitions, namely early to pre-cotyledonary embryo and cotyledonary to mature embryo. An analysis of functional categories for genes that were differentially expressed through embryogenesis highlighted several epigenetic regulation mechanisms. While putative orthologs of transcripts associated with mechanisms that target transposable elements and repetitive sequences were strongly expressed in early embryogenesis, PRC2-mediated repression of genes seemed more relevant during late embryogenesis. On the other hand, functions related to sRNA pathways appeared differentially regulated across all stages of embryo development with a prevalence of miRNA functions in mid to late embryogenesis. Identification of putative transcription factor genes differentially regulated between consecutive embryo stages was strongly suggestive of the relevance of auxin responses and regulation of auxin carriers during early embryogenesis. Such responses could be involved in establishing embryo patterning. Later in development, transcripts with homology to genes acting on modulation of auxin flow and determination of adaxial-abaxial polarity were up-regulated, as were putative orthologs of genes required for meristem formation and function as well as establishment of organ boundaries. Comparative analysis with A. thaliana embryogenesis also highlighted genes involved in auxin-mediated responses, as well as epigenetic regulation, indicating highly correlated transcript profiles between the two species.

CONCLUSIONS

This is the first report of a time-course transcriptomic analysis of zygotic embryogenesis in a conifer. Taken together our results show that epigenetic regulation and transcriptional control related to auxin transport and response are critical during early to mid stages of pine embryogenesis and that important events during embryogenesis seem to be coordinated by putative orthologs of major developmental regulators in angiosperms.

Integrated analysis of drug-induced gene expression profiles predicts novel hERG inhibitors

Growing evidence suggests that drugs interact with diverse molecular targets mediating both therapeutic and toxic effects. Prediction of these complex interactions from chemical structures alone remains challenging, as compounds with different structures may possess similar toxicity profiles. In contrast, predictions based on systems-level measurements of drug effect may reveal pharmacologic similarities not evident from structure or known therapeutic indications. Here we utilized drug-induced transcriptional responses in the Connectivity Map (CMap) to discover such similarities among diverse antagonists of the human ether-à-go-go related (hERG) potassium channel, a common target of promiscuous inhibition by small molecules. Analysis of transcriptional profiles generated in three independent cell lines revealed clusters enriched for hERG inhibitors annotated using a database of experimental measurements (hERGcentral) and clinical indications. As a validation, we experimentally identified novel hERG inhibitors among the unannotated drugs in these enriched clusters, suggesting transcriptional responses may serve as predictive surrogates of cardiotoxicity complementing existing functional assays.

Clustering of High Throughput Gene Expression Data

High throughput biological data need to be processed, analyzed, and interpreted to address problems in life sciences. Bioinformatics, computational biology, and systems biology deal with biological problems using computational methods. Clustering is one of the methods used to gain insight into biological processes, particularly at the genomics level. Clearly, clustering can be used in many areas of biological data analysis. However, this paper presents a review of the current clustering algorithms designed especially for analyzing gene expression data. It is also intended to introduce one of the main problems in bioinformatics - clustering gene expression data - to the operations research community.

Microarray analysis of Etrog citron (Citrus medica L.) reveals changes in chloroplast, cell wall, peroxidase and symporter activities in response to viroid infection

Viroids are small (246-401 nucleotides), single-stranded, circular RNA molecules that infect several crop plants and can cause diseases of economic importance. Citrus are the hosts in which the largest number of viroids have been identified. Citrus exocortis viroid (CEVd), the causal agent of citrus exocortis disease, induces considerable losses in citrus crops. Changes in the gene expression profile during the early (pre-symptomatic) and late (post-symptomatic) stages of Etrog citron infected with CEVd were investigated using a citrus cDNA microarray. MaSigPro analysis was performed and, on the basis of gene expression profiles as a function of the time after infection, the differentially expressed genes were classified into five clusters. FatiScan analysis revealed significant enrichment of functional categories for each cluster, indicating that viroid infection triggers important changes in chloroplast, cell wall, peroxidase and symporter activities.

Transcriptome data analysis for cell culture processes

In the past decade, DNA microarrays have fundamentally changed the way we study complex biological systems. By measuring the expression levels of thousands of transcripts, the paradigm of studying organisms has shifted from focusing on the local phenomena of a few genes to surveying the whole genome. DNA microarrays are used in a variety of ways, from simple comparisons between two samples to more intricate time-series studies. With the large number of genes being studied, the dimensionality of the problem is inevitably high. The analysis of microarray data thus requires specific approaches. In the case of time-series microarray studies, data analysis is further complicated by the correlation between successive time points in a series.In this review, we survey the methodologies used in the analysis of static and time-series microarray data, covering data pre-processing, identification of differentially expressed genes, profile pattern recognition, pathway analysis, and network reconstruction. When available, examples of their use in mammalian cell cultures are presented.

Fortunella margarita transcriptional reprogramming triggered by Xanthomonas citri subsp. citri

BACKGROUND

Citrus canker disease caused by the bacterial pathogen Xanthomonas citri subsp. citri (Xcc) has become endemic in areas where high temperature, rain, humidity, and windy conditions provide a favourable environment for the dissemination of the bacterium. Xcc is pathogenic on many commercial citrus varieties but appears to elicit an incompatible reaction on the citrus relative Fortunella margarita Swing (kumquat), in the form of a very distinct delayed necrotic response. We have developed subtractive libraries enriched in sequences expressed in kumquat leaves during both early and late stages of the disease. The isolated differentially expressed transcripts were subsequently sequenced. Our results demonstrate how the use of microarray expression profiling can help assign roles to previously uncharacterized genes and elucidate plant pathogenesis-response related mechanisms. This can be considered to be a case study in a citrus relative where high throughput technologies were utilized to understand defence mechanisms in Fortunella and citrus at the molecular level.

RESULTS

cDNAs from sequenced kumquat libraries (ESTs) made from subtracted RNA populations, healthy vs. infected, were used to make this microarray. Of 2054 selected genes on a customized array, 317 were differentially expressed (P < 0.05) in Xcc challenged kumquat plants compared to mock-inoculated ones. This study identified components of the incompatible interaction such as reactive oxygen species (ROS) and programmed cell death (PCD). Common defence mechanisms and a number of resistance genes were also identified. In addition, there were a considerable number of differentially regulated genes that had no homologues in the databases. This could be an indication of either a specialized set of genes employed by kumquat in response to canker disease or new defence mechanisms in citrus.

CONCLUSION

Functional categorization of kumquat Xcc-responsive genes revealed an enhanced defence-related metabolism as well as a number of resistant response-specific genes in the kumquat transcriptome in response to Xcc inoculation. Gene expression profile(s) were analyzed to assemble a comprehensive and inclusive image of the molecular interaction in the kumquat/Xcc system. This was done in order to elucidate molecular mechanisms associated with the development of the hypersensitive response phenotype in kumquat leaves. These data will be used to perform comparisons among citrus species to evaluate means to enhance the host immune responses against bacterial diseases.

Large-scale transcriptional profiling and functional assays reveal important roles for Rho-GTPase signalling and SCL during haematopoietic differentiation of human embryonic stem cells

Understanding the transcriptional cues that direct differentiation of human embryonic stem cells (hESCs) and human-induced pluripotent stem cells to defined and functional cell types is essential for future clinical applications. In this study, we have compared transcriptional profiles of haematopoietic progenitors derived from hESCs at various developmental stages of a feeder- and serum-free differentiation method and show that the largest transcriptional changes occur during the first 4 days of differentiation. Data mining on the basis of molecular function revealed Rho-GTPase signalling as a key regulator of differentiation. Inhibition of this pathway resulted in a significant reduction in the numbers of emerging haematopoietic progenitors throughout the differentiation window, thereby uncovering a previously unappreciated role for Rho-GTPase signalling during human haematopoietic development. Our analysis indicated that SCL was the 11th most upregulated transcript during the first 4 days of the hESC differentiation process. Overexpression of SCL in hESCs promoted differentiation to meso-endodermal lineages, the emergence of haematopoietic and erythro-megakaryocytic progenitors and accelerated erythroid differentiation. Importantly, intrasplenic transplantation of SCL-overexpressing hESC-derived haematopoietic cells enhanced recovery from induced acute anaemia without significant cell engraftment, suggesting a paracrine-mediated effect.

Curve-based multivariate distance matrix regression analysis: application to genetic association analyses involving repeated measures

Most, if not all, human phenotypes exhibit a temporal, dosage-dependent, or age effect. Despite this fact, it is rare that data are collected over time or in sequence in relevant studies of the determinants of these phenotypes. The costs and organizational sophistication necessary to collect repeated measurements or longitudinal data for a given phenotype are clearly impediments to this, but greater efforts in this area are needed if insights into human phenotypic expression are to be obtained. Appropriate data analysis methods for genetic association studies involving repeated or longitudinal measures are also needed. We consider the use of longitudinal profiles obtained from fitted functions on repeated data collections from a set of individuals whose similarities are contrasted between sets of individuals with different genotypes to test hypotheses about genetic influences on time-dependent phenotype expression. The proposed approach can accommodate uncertainty of the fitted functions, as well as weighting factors across the time points, and is easily extended to a wide variety of complex analysis settings. We showcase the proposed approach with data from a clinical study investigating human blood vessel response to tyramine. We also compare the proposed approach with standard analytic procedures and investigate its robustness and power via simulation studies. The proposed approach is found to be quite flexible and performs either as well or better than traditional statistical methods.

Serial Expression Analysis: a web tool for the analysis of serial gene expression data

Serial transcriptomics experiments investigate the dynamics of gene expression changes associated with a quantitative variable such as time or dosage. The statistical analysis of these data implies the study of global and gene-specific expression trends, the identification of significant serial changes, the comparison of expression profiles and the assessment of transcriptional changes in terms of cellular processes. We have created the SEA (Serial Expression Analysis) suite to provide a complete web-based resource for the analysis of serial transcriptomics data. SEA offers five different algorithms based on univariate, multivariate and functional profiling strategies framed within a user-friendly interface and a project-oriented architecture to facilitate the analysis of serial gene expression data sets from different perspectives. SEA is available at sea.bioinfo.cipf.es.

The 20th anniversary of EMBnet: 20 years of bioinformatics for the Life Sciences community

The EMBnet Conference 2008, focusing on 'Leading Applications and Technologies in Bioinformatics', was organized by the European Molecular Biology network (EMBnet) to celebrate its 20th anniversary. Since its foundation in 1988, EMBnet has been working to promote collaborative development of bioinformatics services and tools to serve the European community of molecular biology laboratories. This conference was the first meeting organized by the network that was open to the international scientific community outside EMBnet. The conference covered a broad range of research topics in bioinformatics with a main focus on new achievements and trends in emerging technologies supporting genomics, transcriptomics and proteomics analyses such as high-throughput sequencing and data managing, text and data-mining, ontologies and Grid technologies. Papers selected for publication, in this supplement to BMC Bioinformatics, cover a broad range of the topics treated, providing also an overview of the main bioinformatics research fields that the EMBnet community is involved in.