Bioinformatics Advance Access originally published online on March 30, 2007
Bioinformatics 2007 23(12):1559-1561; doi:10.1093/bioinformatics/btm126
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
SCEPTRANS: an online tool for analyzing periodic transcription in yeast
University of Texas Southwestern Medical Center, Department of Biochemistry, 5323 Harry Hines Blvd. Dallas, TX 75390-8816, USA
*To whom correspondence should be addressed.
 |
ABSTRACT |
|---|
 TOP ABSTRACT 1 INTRODUCTION2 IMPLEMENTATION3 DISCUSSIONACKNOWLEDGEMENTREFERENCES |
|---|
Summary: SCEPTRANS is designed for analysis of microarray timecourse data related to periodic phenomena in the budding yeast. The server allows for easy viewing of temporal profiles of multiple genes in a number of datasets. Additional functionality includes searching for coexpressed genes, periodicity and correlation analysis, integrating functional annotation and localization data as well as advanced operations on sets of genes.
Availability: Available online at http://sceptrans.org/
Contact: andrzej{at}work.swmed.edu
Supplementary information: Supplementary data are available at Bioinformatics online.
|
1 INTRODUCTION |
|---|
|
TOPABSTRACT1 INTRODUCTION2 IMPLEMENTATION3 DISCUSSIONACKNOWLEDGEMENTREFERENCES |
|---|
Periodic processes in the living organisms (e.g. cell cycle), involve precisely orchestrated transcription of hundreds of genes. Studying whole-genome time-course expression data provides insight into regulation and function of individual genes in context of cellular processes. A number of microarray datasets related to periodic phenomena in yeast are available (Table 1). The functionality of the existing servers (Pramila et al., 2006; Spellman et al., 1998, SGD) does not allow for precise comparison of expression profiles within a user-defined group of genes or for comprehensive comparison of expression patterns of a group of genes across conditions, both essential for in-depth analysis. The methods of assessing periodicity also vary substantially, from relying on visual examination (Cho et al., 1998) or assuming all genes in such a system are periodic (Klevecz et al., 2004) to various complex approaches (Pramila et al., 2006; Tu et al., 2005), making comparisons between datasets difficult for general users not willing to completely re-analyze each dataset. Here, we present an integrative resource applying the same set of tools to data from different experiments.
|
|
2 IMPLEMENTATION |
|---|
|
TOPABSTRACT1 INTRODUCTION2 IMPLEMENTATION3 DISCUSSIONACKNOWLEDGEMENTREFERENCES |
|---|
The interface is build on an apache/CGI server and consists of two functional modules: gene selection and analysis.
2.1 Gene selection
Gene selection is done by subsequent modification of the current set of genes, by applying the AND, OR, NOT and REPLACE set operators. Genes may be thus added or removed from the current set by one of the following criteria:
- Explicit list of gene names or identifiers (e.g. YNL300w, cln3)
- Regular expression search of gene name or description (e.g. to select autophagy-related genes, one can either search names for ‘ATG[0–9]’ or search descriptions for ‘autophag’ and then remove ‘non-autophag’)
- Correlation with a given gene in one of the experiments
- Peak-to-trough ratio in each experiment
- Periodicity score in one of the experiments
- Cellular localization (e.g. mitochondrion or spindle pole)
- Classified as periodic by different studies
- Functional annotation (general, e.g. metabolism, or detailed, e.g. nitrogen and sulfur metabolism)
2.2 Gene analysis
The gene analysis module allows the user to inspect expression profiles of chosen genes in selected datasets. The following analysis tools are available:
- Table of gene features. The table includes the common and systematic names of a gene, its cellular localization (Huh et al., 2003), functional annotation (Comprehensive Yeast Genome Database), description (Saccharomyces Genome Database), peak-to-trough transcription ratios, cell cycle regulated flags (Cho et al., 1998; Pramila et al., 2006; Rowicka et al., 2007; Spellman et al., 1998) and metabolic regulated flag (Tu et al., 2005).
- Plots of temporal expression profiles of selected genes. Profiles are shown in the linear scale, normalized to unit average or unit maximum. With a limited number of plots displayed at a time, the user can highlight each expression profile for clarity.
- Periodicity analysis. Two measures of periodicity are implemented: the magnitude of the Fourier mode and autocorrelation at different time shifts. The computed Fourier mode (Lomb, 1976; Scargle, 1982), is presented as a periodogram (i.e. using period length as the abscissa). Autocorrelations are computed on normalized profiles interpolated with cubic splines (Press et al., 1992). Periodograms or autocorellograms can be displayed graphically as an alternative to time series plots. Statistical significance of periodograms (Horne and Baliunas, 1986) can also be displayed.
- Correlation table. Correlation tables are available to distinguish co-regulation modules in the group of genes. To facilitate interpretation, correlation tables are displayed both as heat maps and numerical values of Pearson correlation coefficients. Genes may be sorted according to time of expression in each dataset, thus visually separating different expression patterns in the correlation table.
2.3 Limitations
Our implementation of the interface relies on the use of JavaScript, which has to be enabled in the user's browser. To save server-side CPU time and bandwidth, correlation tables will not be displayed for more than 175 transcripts; no more than 100 temporal profiles will be shown at a time (and maximum 50 can be highlighted), and the gene description table will be hidden if it contains more than 1000 lines. These maximum numbers of transcripts may be further reduced if more than one dataset is displayed at a time, or if the current load on the server is too high.
|
3 DISCUSSION |
|---|
|
TOPABSTRACT1 INTRODUCTION2 IMPLEMENTATION3 DISCUSSIONACKNOWLEDGEMENTREFERENCES |
|---|
To exemplify SCEPTRANS functionality, we show how it speeds up the functional analysis we have reported previously (Tu et al., 2005). For instance, lists of genes periodic during yeast metabolic cycle (YMC) with selected cellular localizations can be readily obtained using the ‘localized’ and ‘regulated’ boxes in the gene selection panel. This procedure shows immediately that the most significantly overrepresented localization among the genes found periodic during YMC is mitochondrion (430 out of 523 are periodic in YMC). Similarly, one can obtain that disproportionally large number of genes from the ‘energy’ and ‘metabolism’ functional categories are periodic during YMC (80 and 67% respectively, compared to 52% in the whole genome). These results support the view that the observed oscillations are of metabolic nature (Tu et al., 2005).
Another task much facilitated by our server is the analysis of expression of mitochondrial ribosomal genes (). Genes with this annotation can be selected by just three clicks of the mouse. In the YMC dataset, the expression pattern of one such transcript, PPE1, is significantly different than that of other mitoribosomal transcripts, raising a possibility that PPE1 is either misannotated, or that it is regulated in a different manner than other mitoribosomal transcripts. Indeed, Ppe1 is the only gene annotated as mitochondrial ribosomal subunit that lacks the regulatory motif characteristic of this group (Tu et al., 2005).
|
Another very useful tool is a correlation table: it allows us to see data structure more easily than time-course plots (Fig. 2).
|
In summary, SCEPTRANS provides an opportunity for general researchers to perform a comprehensive yet flexible analysis of various microarray time-course data. We expect that SCEPTRANS will facilitate substantially such analysis and will help to formulate and test novel research hypotheses.
|
ACKNOWLEDGEMENT |
|---|
|
TOPABSTRACT1 INTRODUCTION2 IMPLEMENTATION3 DISCUSSIONACKNOWLEDGEMENTREFERENCES |
|---|
This work was supported by NIH GM U54GM074942. We thank Maciej Puzio for technical assistance.
Conflict of Interest: none declared.
|
FOOTNOTES |
|---|
Associate Editor: Chris Stoeckert
Received on March 2, 2007; revised on March 22, 2007; accepted on March 24, 2007
|
REFERENCES |
|---|
|
TOPABSTRACT1 INTRODUCTION2 IMPLEMENTATION3 DISCUSSIONACKNOWLEDGEMENTREFERENCES |
|---|
Cho R, et al. A genome-wide transcriptional analysis of the mitotic cell cycle. Mol. Cell (1998) 2:65–73.[CrossRef][Web of Science][Medline]
Horne JH, Baliunas SL. Astrophysical J (1986) 302:757.[CrossRef]
Huh W, et al. Global analysis of protein localization in budding yeast. Nature (2003) 425:686.[CrossRef][Medline]
Klevecz R, et al. A genomewide oscillation in transcription gates DNA replication and cell cycle. PNAS (2004) 101:1200–1205.
Lomb N. Astrophy. Space Sci (1976) 39:447.[CrossRef]
Pramila T, et al. The Forkhead transcription factor Hcm1 regulates chromosome segregation genes and fills the Sphase gap in the transcriptional circuitry of the cell cycle. Genes Dev (2006) 20:2266.
Press WH, et al. Numerical Recipes in C. The Art of Scientific Computing. (1992) Cambridge University Press.
Rowicka M, et al. High-resolution timing of the cell-cycle regulated gene expression. (2007) Submitted.
Scargle JD. Astrophys. J (1982) 263:835.[CrossRef][Web of Science]
Spellman P, et al. Comprehensive identification of cell cycle-regulated genes of the yeast S. cerevisiae by microarray hybridization. Mol. Biol. Cell (1998) 9:3273.
Tu BP, et al. Logic of the yeast metabolic cycle: temporal compartmentalization of cellular processes. Science (2005) 310:1152.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  M. J. Brauer, C. Huttenhower, E. M. Airoldi, R. Rosenstein, J. C. Matese, D. Gresham, V. M. Boer, O. G. Troyanskaya, and D. Botstein Coordination of Growth Rate, Cell Cycle, Stress Response, and Metabolic Activity in Yeast Mol. Biol. Cell, January 1, 2008; 19(1): 352 - 367. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. Rowicka, A. Kudlicki, B. P. Tu, and Z. Otwinowski High-resolution timing of cell cycle-regulated gene expression PNAS, October 23, 2007; 104(43): 16892 - 16897. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||