SCEPTRANS: an online tool for analyzing periodic transcription in yeast

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/     

Transcriptional Analysis of the Candida albicans Cell Cycle

We have examined the periodic expression of genes through the cell cycle in cultures of the human pathogenic fungus Candida albicans synchronized by mating pheromone treatment. Close to 500 genes show increased expression during the G1, S, G2, or M transitions of the C. albicans cell cycle. Comparisons of these C. albicans periodic genes with those already found in the budding and fission yeasts and in human cells reveal that of 2200 groups of homologous genes, close to 600 show periodicity in at least one organism, but only 11 are periodic in all four species. Overall, the C. albicans regulatory circuit most closely resembles that of Saccharomyces cerevisiae but contains a simplified structure. Although the majority of the C. albicans periodically regulated genes have homologues in the budding yeast, 20% (100 genes), most of which peak during the G1/S or M/G1 transitions, are unique to the pathogenic yeast.     

Evolution and regulation of cellular periodic processes: a role for paralogues

Several cyclic processes take place within a single organism. For example, the cell cycle is coordinated with the 24 h diurnal rhythm in animals and plants, and with the 40 min ultradian rhythm in budding yeast. To examine the evolution of periodic gene expression during these processes, we performed the first systematic comparison in three organisms (Homo sapiens, Arabidopsis thaliana and Saccharomyces cerevisiae) by using public microarray data. We observed that although diurnal‐regulated and ultradian‐regulated genes are not generally cell‐cycle‐regulated, they tend to have cell‐cycle‐regulated paralogues. Thus, diverged temporal expression of paralogues seems to facilitate cellular orchestration under different periodic stimuli. Lineage‐specific functional repertoires of periodic‐associated paralogues imply that this mode of regulation might have evolved independently in several organisms.     

The latent process decomposition of cDNA microarray data sets

We present a new computational technique (a software implementation, data sets, and supplementary information are available at http://www.enm.bris.ac.uk/lpd/) which enables the probabilistic analysis of cDNA microarray data and we demonstrate its effectiveness in identifying features of biomedical importance. A hierarchical Bayesian model, called Latent Process Decomposition (LPD), is introduced in which each sample in the data set is represented as a combinatorial mixture over a finite set of latent processes, which are expected to correspond to biological processes. Parameters in the model are estimated using efficient variational methods. This type of probabilistic model is most appropriate for the interpretation of measurement data generated by cDNA microarray technology. For determining informative substructure in such data sets, the proposed model has several important advantages over the standard use of dendrograms. First, the ability to objectively assess the optimal number of sample clusters. Second, the ability to represent samples and gene expression levels using a common set of latent variables (dendrograms cluster samples and gene expression values separately which amounts to two distinct reduced space representations). Third, in constrast to standard cluster models, observations are not assigned to a single cluster and, thus, for example, gene expression levels are modeled via combinations of the latent processes identified by the algorithm. We show this new method compares favorably with alternative cluster analysis methods. To illustrate its potential, we apply the proposed technique to several microarray data sets for cancer. For these data sets it successfully decomposes the data into known subtypes and indicates possible further taxonomic subdivision in addition to highlighting, in a wholly unsupervised manner, the importance of certain genes which are known to be medically significant. To illustrate its wider applicability, we also illustrate its performance on a microarray data set for yeast.     

A Stochastic Model Correctly Predicts Changes in Budding Yeast Cell Cycle Dynamics upon Periodic Expression of CLN2

In this study, we focus on a recent stochastic budding yeast cell cycle model. First, we estimate the model parameters using extensive data sets: phenotypes of 110 genetic strains, single cell statistics of wild type and cln3 strains. Optimization of stochastic model parameters is achieved by an automated algorithm we recently used for a deterministic cell cycle model. Next, in order to test the predictive ability of the stochastic model, we focus on a recent experimental study in which forced periodic expression of CLN2 cyclin (driven by MET3 promoter in cln3 background) has been used to synchronize budding yeast cell colonies. We demonstrate that the model correctly predicts the experimentally observed synchronization levels and cell cycle statistics of mother and daughter cells under various experimental conditions (numerical data that is not enforced in parameter optimization), in addition to correctly predicting the qualitative changes in size control due to forced CLN2 expression. Our model also generates a novel prediction: under frequent CLN2 expression pulses, G1 phase duration is bimodal among small-born cells. These cells originate from daughters with extended budded periods due to size control during the budded period. This novel prediction and the experimental trends captured by the model illustrate the interplay between cell cycle dynamics, synchronization of cell colonies, and size control in budding yeast.     

The clustering of regression models method with applications in gene expression data

Identification of differentially expressed genes and clustering of genes are two important and complementary objectives addressed with gene expression data. For the differential expression question, many "per-gene" analytic methods have been proposed. These methods can generally be characterized as using a regression function to independently model the observations for each gene; various adjustments for multiplicity are then used to interpret the statistical significance of these per-gene regression models over the collection of genes analyzed. Motivated by this common structure of per-gene models, we proposed a new model-based clustering method--the clustering of regression models method, which groups genes that share a similar relationship to the covariate(s). This method provides a unified approach for a family of clustering procedures and can be applied for data collected with various experimental designs. In addition, when combined with per-gene methods for assessing differential expression that employ the same regression modeling structure, an integrated framework for the analysis of microarray data is obtained. The proposed methodology was applied to two microarray data sets, one from a breast cancer study and the other from a yeast cell cycle study.     

Genome-wide analysis of gene expression profiles associated with cell cycle transitions in growing organs of Arabidopsis

Organ growth results from the progression of component cells through subsequent phases of proliferation and expansion before reaching maturity. We combined kinematic analysis, flowcytometry, and microarray analysis to characterize cell cycle regulation during the growth process of leaves 1 and 2 of Arabidopsis (Arabidopsis thaliana). Kinematic analysis showed that the epidermis proliferates until day 12; thereafter, cells expand until day 19 when leaves reach maturity. Flowcytometry revealed that endoreduplication occurs from the time cell division rates decline until the end of cell expansion. Analysis of 10 time points with a 6k-cDNA microarray showed that transitions between the growth stages were closely reflected in the mRNA expression data. Subsequent genome-wide microarray analysis on the three main stages allowed us to categorize known cell cycle genes into three major classes: constitutively expressed, proliferative, and inhibitory. Comparison with published expression data obtained from root zones corresponding to similar developmental stages and from synchronized cell cultures supported this categorization and enabled us to identify a high confidence set of 131 proliferation genes. Most of those had an M phase-dependent expression pattern and, in addition to many known cell cycle-related genes, there were at least 90 that were unknown or previously not associated with proliferation.     

Tumor classification by partial least squares using microarray gene expression data

MOTIVATION: One important application of gene expression microarray data is classification of samples into categories, such as the type of tumor. The use of microarrays allows simultaneous monitoring of thousands of genes expressions per sample. This ability to measure gene expression en masse has resulted in data with the number of variables p(genes) far exceeding the number of samples N. Standard statistical methodologies in classification and prediction do not work well or even at all when N < p. Modification of existing statistical methodologies or development of new methodologies is needed for the analysis of microarray data. RESULTS: We propose a novel analysis procedure for classifying (predicting) human tumor samples based on microarray gene expressions. This procedure involves dimension reduction using Partial Least Squares (PLS) and classification using Logistic Discrimination (LD) and Quadratic Discriminant Analysis (QDA). We compare PLS to the well known dimension reduction method of Principal Components Analysis (PCA). Under many circumstances PLS proves superior; we illustrate a condition when PCA particularly fails to predict well relative to PLS. The proposed methods were applied to five different microarray data sets involving various human tumor samples: (1) normal versus ovarian tumor; (2) Acute Myeloid Leukemia (AML) versus Acute Lymphoblastic Leukemia (ALL); (3) Diffuse Large B-cell Lymphoma (DLBCLL) versus B-cell Chronic Lymphocytic Leukemia (BCLL); (4) normal versus colon tumor; and (5) Non-Small-Cell-Lung-Carcinoma (NSCLC) versus renal samples. Stability of classification results and methods were further assessed by re-randomization studies.     

Extracting binary signals from microarray time-course data

This article presents a new method for analyzing microarray time courses by identifying genes that undergo abrupt transitions in expression level, and the time at which the transitions occur. The algorithm matches the sequence of expression levels for each gene against temporal patterns having one or two transitions between two expression levels. The algorithm reports a P-value for the matching pattern of each gene, and a global false discovery rate can also be computed. After matching, genes can be sorted by the direction and time of transitions. Genes can be partitioned into sets based on the direction and time of change for further analysis, such as comparison with Gene Ontology annotations or binding site motifs. The method is evaluated on simulated and actual time-course data. On microarray data for budding yeast, it is shown that the groups of genes that change in similar ways and at similar times have significant and relevant Gene Ontology annotations.     

Meta-analysis of gene expression data: a predictor-based approach

MOTIVATION: With the increasing availability of cancer microarray data sets there is a growing need for integrative computational methods that evaluate multiple independent microarray data sets investigating a common theme or disorder. Meta-analysis techniques are designed to overcome the low sample size typical to microarray experiments and yield more valid and informative results than each experiment separately. RESULTS: We propose a new meta-analysis technique that aims at finding a set of classifying genes, whose expression level may be used to answering the classification question in hand. Specifically, we apply our method to two independent lung cancer microarray data sets and identify a joint core subset of genes which putatively play an important role in tumor genesis of the lung. The robustness of the identified joint core set is demonstrated on a third unseen lung cancer data set, where it leads to successful classification using very few top-ranked genes. Identifying such a set of genes is of significant importance when searching for biologically meaningful biomarkers. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.