BACKGROUND
it is reported that some splicing isoform-specific transcriptional regulations are related to disease  <cit> .
to find disease specific transcriptional regulations, detection of disease specific splice variations is the first step.
however, conventional microarrays that produce gene-level information are not suitable for this purpose.
on the other hand, affymetrix human exon  <dig>  st array can measure exon-level expression profiles that are suitable to find differentially expressed exons in genome-wide scale.
affymetrix exon array can measure the transcript levels of more than  <dig> , <dig> exons with  <dig>  transcripts by about  <dig> , <dig> probes.
we have developed a supercomputer-based web service named exonminer to analyze exon array datasets for detecting genes that are spliced into different isoforms in two types of cells in comparison, e.g.
normal and cancer cells.
there are some noncommercial standalone applications for analyzing exon array data: igb  <cit>  is an application for visualizing exon array data and exact  <cit>  and affymetrix expression console  <cit>  are mainly focusing on normalizing exon array data.
also, bioconductor  <cit>   focuses on annotation as well as normalization.
the advantage of exonmap is that users can use other statistical tools implemented on r. these are well organized applications, however, these applications focus on data normalizations and we need to use other software for further analysis.
since exonminer is, however, an all-in-one web service on a supercomputer system, users can analyze more than  <dig>  transcripts spotted on exon array by data normalization, two-way anova analysis, visualization of the results, and detection of exon-level biomarkers.
based on our experiments, which used colon cancer exon array data that contains  <dig> exon arrays, on various situations of our system usages, the minimal computational time is four hours and the longest was finished in one day.
we also observed that the average computational time of colon cancer example is about eight hours.
we have implemented exonminer on our super computer system  in human genome center, institute of medical science, university of tokyo and created gui to use the all analysis tools of exonminer easily.
an illustrative example of colon cancer exon array data analysis  <cit>  is shown in the web site.
exonminer has five advantages:  a statistical analysis framework,  analysis for all transcripts completed,  effective visualization with heatmap and barplot images,  sophisticated and easy-to-use web interface, and  useful hyperlinks to major public databases, e.g.
pubmed and netaffx.
as shown in latter sections, the method implemented in exonminer requires more computational time than other software, due to the nonparametric test based on bootstrapping.
for example, we need to repeat bootstrap sampling more than  <dig>  times for computing accurate p-values of statistical tests finding aberrant splice variations, it requires  <dig>  times computation of usual statistical test of anova with gaussian error model.
therefore, we need high-performance parallel computing on super computer system.
also, more advanced methods implemented on exonminer in future possibly requires more computational resources, therefore, the use of super computer system can give flexible computational basis and is suitable for our purpose.
data normalization
before performing statistical analysis, we apply normalization method to raw exon array data.
exonminer can remove a bias related to gc-content in each probe.
the probes are categorized according to their gc-contents and gc-content specific bias will be removed from the probes in each category.
exonminer uses two types of control for data normalization: one is the median value for each gc category and the other is based on antigenomic background probes.
the antigenomic background probes are also categorized into gc categories and we compute their median values.
the median value of the probe intensities in each gc category will be transformed by subtracting corresponding control value.
in case that user chooses the median values of gc categories for control, the median of probe intensities in a gc category will be equal to one.
two-way analysis of variance
concept and model
for using exonminer to detect aberrant splice variations, user needs to prepare at least two exon array data from a pair of cells.
for example, in our illustrative example, one exon array is prepared for measuring exon profiles in colon cancer cell and the other exon array is used for normal cell.
in this case, we can find aberrant splice variants in colon cancer by comparing with normal cells.
in this purpose, we use two-way analysis of variance .
suppose that a gene  is composed of the m exonic regions , and that xijk is the background corrected probe intensity for the k th probe  on the i th exon  of a transcript, i.e.
this transcript has m exonic regions and each exonic region is spanned by nij probes.
here the index j denotes the type of cells, e.g.
j =  <dig> denotes normal cell and j =  <dig> for cancer cell.
if we observe xijk ≈ c for any i, j and k, the transcript does not show any transcriptional changes and splicing variations across cell types .
if we observe that xi1k ≈ c <dig> and xi2k ≈ c <dig>  for any i and k, it indicates that this transcript was differentially expressed between two cells and this information is equivalent to usual microarray expression data like cdna microarray, genechip and so on.
on the other hand, xijk ≈ c <dig> and xi'jk ≈ c <dig> for any j, k and i ≠ i' hold where c <dig> ≠ c <dig> and c <dig>  c <dig> >  <dig>  it means that this transcript has splice variations but these splice variations are commonly occurred between cell types.
finally, if we observe that two cells show different splice patterns, we define them aberrant splice variations.
we will capture this information by two-way anova model.
for anova in exon array data analysis, see also  <cit> .
for detecting transcripts that show aberrant splice variations, we use two-way anova model defined by
xijk = μ + αi + βj + γij + δijk,
where αi, βj and γij are parameters, εijk denotes the observational noise having zero mean and variance σ <dig>  and μ represents an overall mean of the probe intensities.
the parameter αi represents the baseline intensities in the i th exonic region , this parameter captures exon effect.
the parameters βj  capture difference in the overall means between two cells, this difference is called overall gene effect.
the γij s represent interaction effects for each combination of m exons and cell types, which is called effect of specific splice variations.
the effects of these parameters are shown in figure  <dig>  a given statistical evidence that one or more γij s are different with the others suggests that alternative splicing is present in a particular cell, but absent in the other.
we should note that midas  <cit>  is a similar method that uses anova model to analyze exon array data, but midas uses exon-level summarized data, while our model uses probe-level data.
also nonparametric test based on bootstrap method can be considered our advantage.
statistical tests for detecting alternative splicing, differentially expression, and aberrant splice variations
the estimates of γij s could capture presence of aberrant splice variations.
by the anova model, the probe fluctuations are decomposed into three orthogonal effects, i.e.
exon effect , overall gene effect  and effect of specific splice variations .
the statistical significance of each effect can be evaluated by the following three tests:
test  <dig> :
h0: αi =  <dig> for all i.
ha: αi ≠  <dig> at least one i.
test  <dig> :
h0: β <dig> = β2
ha: β <dig> ≠ β2
test  <dig> :
h0: γij =  <dig> for all i and j.
ha: γij ≠  <dig> for one or more pairs of .
here h <dig> and ha represent null and alternative hypotheses, respectively.
repeating these hypothesis tests for all transcript clusters, one can obtain the statistical evidences of aberrant splice variations which are scored by the computed p-values from test  <dig>  in exonminer, in addition to the usual f-test for test of parameter significance, the permutation method that is a nonparametric approach is developed to calculate the null distribution of the f-statistics; fexon, fgene and fsas, for assessing significance of exon effect, gene effect and effect of specific splice variations, respectively.
in order to evaluate the null distributions, we first generate a permutation set of samples by bootstrapping n = Σijnij samples from xijk s. repeating this process b times, we can approximately evaluate the null distribution of each f* with the q permutation statistics f0∗, q =  <dig>  ⋯, q .
note that * can be replaced by exon, gene and sas.
subsequently, the p-value for a given test statistic f* = f* obtained from the original data set is calculated by
p∗=#{q:f0∗≥f∗}q
for each effect.
in exonminer, users can choose parametric or nonparametric test for assessing significance of each parameter.
meta analysis
to detect aberrant splice variants as biomarkers, we need to check whether the detected aberrant splice variants are common in the targeted disease or not.
in this purpose, we establish a statistical testing procedure based on meta-analysis  <cit> .
suppose that we have g pair of exon array datasets, i.e.
normal and tumor exon expression data are measured from g patients.
by performing the whole transcript analysis based on two-way anova to g paired exon array datasets, one obtains a set of p-values for each effect, e.g.
effect of specific splice variations , p1g,⋯,prg, across patients, g =  <dig>  ⋯, g. here the total number of transcripts analyzed is denoted by r. intuitively, a transcript having a small p-value is strongly associated with tumor formation.
however, it is possible that some observed aberrant splice variants could be caused by the inter-individual differences of the analyzed samples.
our goal is to discover the "universal biomarkers", i.e.
aberrant splice variations which are shared by most individuals with a particular diagnostic category.
following this direction, we develop the statistical technique within the framework of meta-analysis based on the normal inversion method.
let xijkg denote the observed probe intensities of the k th probe which spans the i th exonic region for normal cell  or target cell  isolated from the g th individual.
we assume that the probe intensities xijkg of each individual can be modeled by the two-way anova defined by
xijkg=μg+αig+βjg+γijg+εijkg,
for g =  <dig>  ⋯, g. given these models, the statistical hypothesis testing of each effect, for example, effect of specific splice variations, is formulated by
test  <dig> :
h0: γijg =  <dig> for all i, j and k.
ha: γijg ≠  <dig> for one or more tuple .
in order to assess the h <dig>  we propose use of the normal inversion metric as a test statistic.
suppose that we have a set of p-values, ph <dig> ⋯,phg, for occurrence of the aberrant splice variations in the h th transcript cluster.
the method first converts these p-values into the z-scores as zhg=Φ− <dig>  where Φ- <dig> is the inversion of the standard normal cumulative distribution function, and then computed integrated z-score as
zh=∑g=1gzhgg.
the significance of ha can be assessed based on the integrated p-value which is computed by transforming the z-score with the standard normal cumulative distribution function Φ as
phintegrated=1−Φ.
we would like to show an actual example of meta-anlaysis.
in yoshida et al.
<cit> , colon cancer exon array dataset was analyzed by primary version of exonminer.
in this anlaysis, based on the test  <dig> of anova, gene muc <dig>  has p-values for ten patients:
ph <dig> =  <dig> ; ph <dig> =  <dig> ; ph <dig> =  <dig> ; ph <dig> =  <dig> ; ph <dig> =  <dig> ;
ph <dig> =  <dig> ; ph <dig> =  <dig> ; ph <dig> =  <dig> ; ph <dig> =  <dig> ; ph <dig> =  <dig> .
these p-values are transformed into z-scores as:
zh <dig> =  <dig> ; zh <dig> =  <dig> ; zh <dig> =  <dig> ; zh <dig> = - <dig> ; zh <dig> = - <dig> ;
zh <dig> =  <dig> ; zh <dig> = - <dig> ; zh <dig> =  <dig> ; zh <dig> =  <dig> ; zh <dig> =  <dig> .
the integrated z-score is  <dig>  and the integrated p-value is obtained as  <dig>  × 10- <dig>
in the colon cancer example, we compute q-values from integrated p-values of meta-analysis, the list of the genes identified as having aberrant splice variations including exon skipping/retaining has 10% false discovery rate  that corresponds to q-value <  <dig> .
in the above muc <dig>  the q-value is  <dig>  and it is determined as significant.
the computation of q-value is shown in yoshida et al.
<cit> .
by using exon array data with exonminer, it is possible to detect alternative splicing like exon skipping/retaining, alternative usage of donor and acceptor splice sites and so on.
however, since exon array does not have junction probes, custom array with junction probes or pcr method are needed for further analysis of detecting exact patterns of splice isoforms.
implementation
data upload
the users are required to upload their exon array data.
we prepared an ftp service for data upload.
a reason for choosing ftp service for our system is that a large dataset can easily be uploaded.
to increase the security level, we prepare one time account and password for ftp service.
note that one time account and password are different from the pair for login account and password of exonminer.
statistical analysis engine
exonminer performs anova for each transcript.
to test the significance of each effect in anova described in previous section, we implemented two types of tests: one is based on gaussian noise model and it performs f-test, the other is based on nonparametric approach using bootstrap method.
in the nonparametric approach, we need to compute test statistics repeatedly and it needs enormous computation.
therefore we implemented the anova program by fortran and optimized for high performance computing described in the latter section.
visualization engine
the information of exon expression pattern for each transcript needs to be shown visually.
we have developed two types of image generators and can make heatmap and barplot images optimized for exon array data.
these images are generated by using r. the graphics library is originally developed.
database
for the management of user information and probe annotation information, we use mysql database server.
for constructing a highly secure system, user login information is encrypted and stored in mysql database.
by keeping probe annotation information into mysql database, users are not necessary to explore other databases.
thus exonminer is an all-in-one web service.
high performance computing on supercomputers
since anova for the full set of transcripts needs high performance computing, we perform each anova computation in parallel on our supercomputer system.
our supercomputer system has eight sun fire  <dig> k and at most  <dig> cpus can be used for parallel statistical computation by using sun grid engine.
web interface
in exonminer, php scripts deal with connections between front end users and our supercomputer system and dynamically generate images by executing visualization engine described above based on user input.
php scripts generate html web pages with a uniformed style that increases usability.
results and discussions
overview of exonminer
create user account
ftp for data upload
for the upload of your data, you need to use ftp.
for using ftp service in exonminer, user needs to get one time password and account for ftp.
note that the account of ftp is different from login account.
using the one time password, the user can upload cel  files archived by zip via ftp.
exonminer supports cel files as text format , usual cel files are, however, binary format .
to convert binary cel files  in text format , "cel file conversion tool" provided by affymetrix inc. is available  <cit> .
analysis options
user should fill up all of the analysis options.
then user will start the analysis.
user must select all  –  options in figure  <dig> to start a statistical analysis by two-way anova and meta-analysis.
description: you can add a brief description of your analysis.
it may be convenient that you put a name of this analysis to organize your analyses.
select probe levels: you can select the level of expression information in exon array.
transcript level: for transcripts, there are three levels, core, extended and full transcripts, according to their information quality based on their information sources.
like transcript level, user can choose probe level and exon resolution.
select gff: you can select chromosomes.
transcripts on the selected chromosomes will be analyzed.
this selection can reduce computational time.
select which cel file is a patient or a control: user adds the outcome information to each cel file you have uploaded by ftp.
preprocessing data : user selects the type of normalization method.
gc-content: the median values in the same gc-content probe groups are used as control values.
antigenomic background: the median values in the same gc-content antigenomic background probes are used as control values.
preprocessing data : it is a possible case that probes with high gc-content work as noise.
so you can remove such probes.
in default, the probes with  <dig> or more gc-content are removed.
if you want to use the all probes for analysis, you choose  <dig> as the cut-off.
analysis type : user selects the analysis type from the following three types – don't analyze: exonminer does not perform anova.
only visualization and sequence information are available.
parametric analysis: gaussian distribution is assumed as the noise model.
nonparametric analysis: exonminer does not assume any distributions for the noise model.
bootstrap test will be applied for computing p-values.
analysis type : exonminer ignores probesets or exon clusters with small number of probes for stabilizing the results of anova.
you can choose this cut-off by this option.
nonparametric analysis options: the number of bootstraps in nonparametric anova is specified by this option.
visualization of the
RESULTS
setting the all options, user can start the analysis.
when the analysis is completed, exonminer sends an e-mail to the user to announce that the calculation is finished.
after that, the user can view result pages of the analysis with heatmaps, barplots, sequence information and calculated p-values of two-way anova and results of meta-analysis.
a screen shot of exonminer is given in figure  <dig>  in this figure, you can see the results of lgr <dig>  lgr <dig> is one of the most significant genes in colon cancer exon arrays reported by yoshida et al.
<cit> .
the colon cancer exon array data are provided by affymetrix.
we can reach the information for each transcript by either gene symbol or transcription cluster id.
the heatmap  represents the exon profiles of lgr <dig>  the user can download the heatmap image as bitmap or postscript file.
sequence information  for the transcript is shown with hyperlinks to the external web sites, entrez  <cit>  and netaffx  <cit> .
the table  shows calculated anova p-values.
user can view the barplot image of normalized exon expression for a pair of cells from the view hyperlinks.
the p-values for parameters calculated in meta-analysis are shown in the bottom table.
the user can download results in one excel file.
instead of the heatmap image, exonminer can produce barplot images.
figure  <dig> is a barplot image for lgr <dig>  a barplot image has three bar-graphs.
red bar-graph shows probe intensities in exon array of colon cancer cell and green bar-graph shows probe intensities in exon array of normal cell.
we show the bars with lower intensities in dark color.
if the color of the bar on a dark bar is red, the cell type of the dark bar is normal  and vice versa.
by using dark bar-graph, the users easily find the differences of exon expressions between two cells.
for example, from figure  <dig>  we can find that the exon expression levels of colon cancer cell are higher than those of normal cell in many exonic regions.
availability and requirements
• project name: exonminer
• project home page:
• anonymous accounts :
• operating systems: any os
• programming language: php, r, fortran, perl, ruby, mysql
CONCLUSIONS
exonminer is an all-in-one web service well suited for analysis of exon array data.
since it does not require any installation of software except for internet browsers, what all users need to do is to access the exonminer url .
exonminer can perform not only visualization of exon array data, but also can perform data normalization and user customized statistical analysis that is hard to run on a single computer.
with the support of supercomputers in human genome center, institute of medical science, university of tokyo, users can analyze full dataset of exon array data within hours with results of meta-analysis that finds aberrant splice variants as biomarkers.
authors' contributions
kn, as and mn designed exonminer and kn implemented.
kn and as prepared the figures.
ry and si developed statistical analysis in exonminer.
sm supervised the project.
kn wrote the manuscript.

 it is reported that some splicing isoform-specific transcriptional regulations are related to disease. to find disease specific transcriptional regulations, detection of disease specific splice variations is the first step. however, conventional microarrays that produce gene-level information are not suitable for this purpose. on the other hand, affymetrix human exon st array can measure exon-level expression profiles that are suitable to find differentially expressed exons in genome-wide scale. affymetrix exon array can measure the transcript levels of more than, exons with transcripts by about, probes. we have developed a supercomputer-based web service named exonminer to analyze exon array datasets for detecting genes that are spliced into different isoforms in two types of cells in comparison, e.g. normal and cancer cells. there are some noncommercial standalone applications for analyzing exon array data: igb is an application for visualizing exon array data and exact and affymetrix expression console are mainly focusing on normalizing exon array data. also, bioconductor focuses on annotation as well as normalization. the advantage of exonmap is that users can use other statistical tools implemented on r. these are well organized applications, however, these applications focus on data normalizations and we need to use other software for further analysis. since exonminer is, however, an all-in-one web service on a supercomputer system, users can analyze more than transcripts spotted on exon array by data normalization, two-way anova analysis, visualization of the results, and detection of exon-level biomarkers. based on our experiments, which used colon cancer exon array data that contains exon arrays, on various situations of our system usages, the minimal computational time is four hours and the longest was finished in one day. we also observed that the average computational time of colon cancer example is about eight hours. we have implemented exonminer on our super computer system in human genome center, institute of medical science, university of tokyo and created gui to use the all analysis tools of exonminer easily. an illustrative example of colon cancer exon array data analysis is shown in the web site. exonminer has five advantages: a statistical analysis framework, analysis for all transcripts completed, effective visualization with heatmap and barplot images, sophisticated and easy-to-use web interface, and useful hyperlinks to major public databases, e.g. pubmed and netaffx. as shown in latter sections, the method implemented in exonminer requires more computational time than other software, due to the nonparametric test based on bootstrapping. for example, we need to repeat bootstrap sampling more than times for computing accurate p-values of statistical tests finding aberrant splice variations, it requires times computation of usual statistical test of anova with gaussian error model. therefore, we need high-performance parallel computing on super computer system. also, more advanced methods implemented on exonminer in future possibly requires more computational resources, therefore, the use of super computer system can give flexible computational basis and is suitable for our purpose. data normalization before performing statistical analysis, we apply normalization method to raw exon array data. exonminer can remove a bias related to gc-content in each probe. the probes are categorized according to their gc-contents and gc-content specific bias will be removed from the probes in each category. exonminer uses two types of control for data normalization: one is the median value for each gc category and the other is based on antigenomic background probes. the antigenomic background probes are also categorized into gc categories and we compute their median values. the median value of the probe intensities in each gc category will be transformed by subtracting corresponding control value. in case that user chooses the median values of gc categories for control, the median of probe intensities in a gc category will be equal to one. two-way analysis of variance concept and model for using exonminer to detect aberrant splice variations, user needs to prepare at least two exon array data from a pair of cells. for example, in our illustrative example, one exon array is prepared for measuring exon profiles in colon cancer cell and the other exon array is used for normal cell. in this case, we can find aberrant splice variants in colon cancer by comparing with normal cells. in this purpose, we use two-way analysis of variance. suppose that a gene is composed of the m exonic regions, and that xijk is the background corrected probe intensity for the k th probe on the i th exon of a transcript, i.e. this transcript has m exonic regions and each exonic region is spanned by nij probes. here the index j denotes the type of cells, e.g. j = denotes normal cell and j = for cancer cell. if we observe xijk ≈ c for any i, j and k, the transcript does not show any transcriptional changes and splicing variations across cell types. if we observe that xi1k ≈ c and xi2k ≈ c for any i and k, it indicates that this transcript was differentially expressed between two cells and this information is equivalent to usual microarray expression data like cdna microarray, genechip and so on. on the other hand, xijk ≈ c and xi'jk ≈ c for any j, k and i ≠ i' hold where c ≠ c and c c > it means that this transcript has splice variations but these splice variations are commonly occurred between cell types. finally, if we observe that two cells show different splice patterns, we define them aberrant splice variations. we will capture this information by two-way anova model. for anova in exon array data analysis, see also. for detecting transcripts that show aberrant splice variations, we use two-way anova model defined by xijk = μ + αi + βj + γij + δijk, where αi, βj and γij are parameters, εijk denotes the observational noise having zero mean and variance σ and μ represents an overall mean of the probe intensities. the parameter αi represents the baseline intensities in the i th exonic region, this parameter captures exon effect. the parameters βj capture difference in the overall means between two cells, this difference is called overall gene effect. the γij s represent interaction effects for each combination of m exons and cell types, which is called effect of specific splice variations. the effects of these parameters are shown in figure a given statistical evidence that one or more γij s are different with the others suggests that alternative splicing is present in a particular cell, but absent in the other. we should note that midas is a similar method that uses anova model to analyze exon array data, but midas uses exon-level summarized data, while our model uses probe-level data. also nonparametric test based on bootstrap method can be considered our advantage. statistical tests for detecting alternative splicing, differentially expression, and aberrant splice variations the estimates of γij s could capture presence of aberrant splice variations. by the anova model, the probe fluctuations are decomposed into three orthogonal effects, i.e. exon effect, overall gene effect and effect of specific splice variations. the statistical significance of each effect can be evaluated by the following three tests: test: h0: αi = for all i. ha: αi ≠ at least one i. test: h0: β = β2 ha: β ≠ β2 test: h0: γij = for all i and j. ha: γij ≠ for one or more pairs of. here h and ha represent null and alternative hypotheses, respectively. repeating these hypothesis tests for all transcript clusters, one can obtain the statistical evidences of aberrant splice variations which are scored by the computed p-values from test in exonminer, in addition to the usual f-test for test of parameter significance, the permutation method that is a nonparametric approach is developed to calculate the null distribution of the f-statistics; fexon, fgene and fsas, for assessing significance of exon effect, gene effect and effect of specific splice variations, respectively. in order to evaluate the null distributions, we first generate a permutation set of samples by bootstrapping n = Σijnij samples from xijk s. repeating this process b times, we can approximately evaluate the null distribution of each f* with the q permutation statistics f0∗, q = ⋯, q. note that * can be replaced by exon, gene and sas. subsequently, the p-value for a given test statistic f* = f* obtained from the original data set is calculated by p∗=#{q:f0∗≥f∗}q for each effect. in exonminer, users can choose parametric or nonparametric test for assessing significance of each parameter. meta analysis to detect aberrant splice variants as biomarkers, we need to check whether the detected aberrant splice variants are common in the targeted disease or not. in this purpose, we establish a statistical testing procedure based on meta-analysis. suppose that we have g pair of exon array datasets, i.e. normal and tumor exon expression data are measured from g patients. by performing the whole transcript analysis based on two-way anova to g paired exon array datasets, one obtains a set of p-values for each effect, e.g. effect of specific splice variations, p1g,⋯,prg, across patients, g = ⋯, g. here the total number of transcripts analyzed is denoted by r. intuitively, a transcript having a small p-value is strongly associated with tumor formation. however, it is possible that some observed aberrant splice variants could be caused by the inter-individual differences of the analyzed samples. our goal is to discover the "universal biomarkers", i.e. aberrant splice variations which are shared by most individuals with a particular diagnostic category. following this direction, we develop the statistical technique within the framework of meta-analysis based on the normal inversion method. let xijkg denote the observed probe intensities of the k th probe which spans the i th exonic region for normal cell or target cell isolated from the g th individual. we assume that the probe intensities xijkg of each individual can be modeled by the two-way anova defined by xijkg=μg+αig+βjg+γijg+εijkg, for g = ⋯, g. given these models, the statistical hypothesis testing of each effect, for example, effect of specific splice variations, is formulated by test: h0: γijg = for all i, j and k. ha: γijg ≠ for one or more tuple. in order to assess the h we propose use of the normal inversion metric as a test statistic. suppose that we have a set of p-values, ph ⋯,phg, for occurrence of the aberrant splice variations in the h th transcript cluster. the method first converts these p-values into the z-scores as zhg=Φ− where Φ- is the inversion of the standard normal cumulative distribution function, and then computed integrated z-score as zh=∑g=1gzhgg. the significance of ha can be assessed based on the integrated p-value which is computed by transforming the z-score with the standard normal cumulative distribution function Φ as phintegrated=1−Φ. we would like to show an actual example of meta-anlaysis. in yoshida et al., colon cancer exon array dataset was analyzed by primary version of exonminer. in this anlaysis, based on the test of anova, gene muc has p-values for ten patients: ph =; ph =; ph =; ph =; ph =; ph =; ph =; ph =; ph =; ph =. these p-values are transformed into z-scores as: zh =; zh =; zh =; zh = -; zh = -; zh =; zh = -; zh =; zh =; zh =. the integrated z-score is and the integrated p-value is obtained as × 10- in the colon cancer example, we compute q-values from integrated p-values of meta-analysis, the list of the genes identified as having aberrant splice variations including exon skipping/retaining has 10% false discovery rate that corresponds to q-value <. in the above muc the q-value is and it is determined as significant. the computation of q-value is shown in yoshida et al.. by using exon array data with exonminer, it is possible to detect alternative splicing like exon skipping/retaining, alternative usage of donor and acceptor splice sites and so on. however, since exon array does not have junction probes, custom array with junction probes or pcr method are needed for further analysis of detecting exact patterns of splice isoforms. implementation data upload the users are required to upload their exon array data. we prepared an ftp service for data upload. a reason for choosing ftp service for our system is that a large dataset can easily be uploaded. to increase the security level, we prepare one time account and password for ftp service. note that one time account and password are different from the pair for login account and password of exonminer. statistical analysis engine exonminer performs anova for each transcript. to test the significance of each effect in anova described in previous section, we implemented two types of tests: one is based on gaussian noise model and it performs f-test, the other is based on nonparametric approach using bootstrap method. in the nonparametric approach, we need to compute test statistics repeatedly and it needs enormous computation. therefore we implemented the anova program by fortran and optimized for high performance computing described in the latter section. visualization engine the information of exon expression pattern for each transcript needs to be shown visually. we have developed two types of image generators and can make heatmap and barplot images optimized for exon array data. these images are generated by using r. the graphics library is originally developed. database for the management of user information and probe annotation information, we use mysql database server. for constructing a highly secure system, user login information is encrypted and stored in mysql database. by keeping probe annotation information into mysql database, users are not necessary to explore other databases. thus exonminer is an all-in-one web service. high performance computing on supercomputers since anova for the full set of transcripts needs high performance computing, we perform each anova computation in parallel on our supercomputer system. our supercomputer system has eight sun fire k and at most cpus can be used for parallel statistical computation by using sun grid engine. web interface in exonminer, php scripts deal with connections between front end users and our supercomputer system and dynamically generate images by executing visualization engine described above based on user input. php scripts generate html web pages with a uniformed style that increases usability. results and discussions overview of exonminer create user account ftp for data upload for the upload of your data, you need to use ftp. for using ftp service in exonminer, user needs to get one time password and account for ftp. note that the account of ftp is different from login account. using the one time password, the user can upload cel files archived by zip via ftp. exonminer supports cel files as text format, usual cel files are, however, binary format. to convert binary cel files in text format, "cel file conversion tool" provided by affymetrix inc. is available. analysis options user should fill up all of the analysis options. then user will start the analysis. user must select all – options in figure to start a statistical analysis by two-way anova and meta-analysis. description: you can add a brief description of your analysis. it may be convenient that you put a name of this analysis to organize your analyses. select probe levels: you can select the level of expression information in exon array. transcript level: for transcripts, there are three levels, core, extended and full transcripts, according to their information quality based on their information sources. like transcript level, user can choose probe level and exon resolution. select gff: you can select chromosomes. transcripts on the selected chromosomes will be analyzed. this selection can reduce computational time. select which cel file is a patient or a control: user adds the outcome information to each cel file you have uploaded by ftp. preprocessing data: user selects the type of normalization method. gc-content: the median values in the same gc-content probe groups are used as control values. antigenomic background: the median values in the same gc-content antigenomic background probes are used as control values. preprocessing data: it is a possible case that probes with high gc-content work as noise. so you can remove such probes. in default, the probes with or more gc-content are removed. if you want to use the all probes for analysis, you choose as the cut-off. analysis type: user selects the analysis type from the following three types – don't analyze: exonminer does not perform anova. only visualization and sequence information are available. parametric analysis: gaussian distribution is assumed as the noise model. nonparametric analysis: exonminer does not assume any distributions for the noise model. bootstrap test will be applied for computing p-values. analysis type: exonminer ignores probesets or exon clusters with small number of probes for stabilizing the results of anova. you can choose this cut-off by this option. nonparametric analysis options: the number of bootstraps in nonparametric anova is specified by this option. visualization of the  setting the all options, user can start the analysis. when the analysis is completed, exonminer sends an e-mail to the user to announce that the calculation is finished. after that, the user can view result pages of the analysis with heatmaps, barplots, sequence information and calculated p-values of two-way anova and results of meta-analysis. a screen shot of exonminer is given in figure in this figure, you can see the results of lgr lgr is one of the most significant genes in colon cancer exon arrays reported by yoshida et al.. the colon cancer exon array data are provided by affymetrix. we can reach the information for each transcript by either gene symbol or transcription cluster id. the heatmap represents the exon profiles of lgr the user can download the heatmap image as bitmap or postscript file. sequence information for the transcript is shown with hyperlinks to the external web sites, entrez and netaffx. the table shows calculated anova p-values. user can view the barplot image of normalized exon expression for a pair of cells from the view hyperlinks. the p-values for parameters calculated in meta-analysis are shown in the bottom table. the user can download results in one excel file. instead of the heatmap image, exonminer can produce barplot images. figure is a barplot image for lgr a barplot image has three bar-graphs. red bar-graph shows probe intensities in exon array of colon cancer cell and green bar-graph shows probe intensities in exon array of normal cell. we show the bars with lower intensities in dark color. if the color of the bar on a dark bar is red, the cell type of the dark bar is normal and vice versa. by using dark bar-graph, the users easily find the differences of exon expressions between two cells. for example, from figure we can find that the exon expression levels of colon cancer cell are higher than those of normal cell in many exonic regions. availability and requirements • project name: exonminer • project home page: • anonymous accounts: • operating systems: any os • programming language: php, r, fortran, perl, ruby, mysql  exonminer is an all-in-one web service well suited for analysis of exon array data. since it does not require any installation of software except for internet browsers, what all users need to do is to access the exonminer url. exonminer can perform not only visualization of exon array data, but also can perform data normalization and user customized statistical analysis that is hard to run on a single computer. with the support of supercomputers in human genome center, institute of medical science, university of tokyo, users can analyze full dataset of exon array data within hours with results of meta-analysis that finds aberrant splice variants as biomarkers. authors' contributions kn, as and mn designed exonminer and kn implemented. kn and as prepared the figures. ry and si developed statistical analysis in exonminer. sm supervised the project. kn wrote the manuscript. 
we implemented exonminer on a supercomputer system in order to perform genome-wide analysis for more than transcripts in exon array data, which has the potential to reveal the aberrant splice variations in cancer cells as exon level biomarkers. exonminer is well suited for analysis of exon array data and does not require any installation of software except for internet browsers. exonminer can perform the following analyses: data normalization, statistical analysis based on two-way anova, finding transcripts with significantly different splice patterns, efficient visualization based on heatmaps and barplots, and meta-analysis to detect exon level biomarkers. users can analyze full dataset of exon array data within hours by high-level statistical analysis with sound theoretical basis that finds aberrant splice variants as biomarkers. some splicing isoform-specific transcriptional regulations are related to disease. therefore, detection of disease specific splice variations is the first step for finding disease specific transcriptional regulations. affymetrix human exon st array can measure exon-level expression profiles that are suitable to find differentially expressed exons in genome-wide scale. what all users need to do is to access the exonminer url. however, exon array produces massive datasets that are more than we can handle and analyze on personal computer. we have developed exonminer that is the first all-in-one web service for analysis of exon array data to detect transcripts that have significantly different splicing patterns in two cells, e.g.
 some splicing isoform-specific transcriptional regulations are related to disease. therefore, detection of disease specific splice variations is the first step for finding disease specific transcriptional regulations. affymetrix human exon st array can measure exon-level expression profiles that are suitable to find differentially expressed exons in genome-wide scale. however, exon array produces massive datasets that are more than we can handle and analyze on personal computer.  we have developed exonminer that is the first all-in-one web service for analysis of exon array data to detect transcripts that have significantly different splicing patterns in two cells, e.g. normal and cancer cells. exonminer can perform the following analyses: data normalization, statistical analysis based on two-way anova, finding transcripts with significantly different splice patterns, efficient visualization based on heatmaps and barplots, and meta-analysis to detect exon level biomarkers. we implemented exonminer on a supercomputer system in order to perform genome-wide analysis for more than transcripts in exon array data, which has the potential to reveal the aberrant splice variations in cancer cells as exon level biomarkers.  exonminer is well suited for analysis of exon array data and does not require any installation of software except for internet browsers. what all users need to do is to access the exonminer url. users can analyze full dataset of exon array data within hours by high-level statistical analysis with sound theoretical basis that finds aberrant splice variants as biomarkers. 