EXPath tool-a system for comprehensively analyzing regulatory pathways and coexpression networks from high-throughput transcriptome data.

Next generation sequencing (NGS) has become the mainstream approach for monitoring gene expression levels in parallel with various experimental treatments. Unfortunately, there is no systematical webserver to comprehensively perform further analysis based on the huge amount of preliminary data that is obtained after finishing the process of gene annotation. Therefore, a user-friendly and effective system is required to mine important genes and regulatory pathways under specific conditions from high-throughput transcriptome data. EXPath Tool (available at: http://expathtool.itps.ncku.edu.tw/) was developed for the pathway annotation and comparative analysis of user-customized gene expression profiles derived from microarray or NGS platforms under various conditions to infer metabolic pathways for all organisms in the KEGG database. EXPath Tool contains several functions: access the gene expression patterns and the candidates of co-expression genes; dissect differentially expressed genes (DEGs) between two conditions (DEGs search), functional grouping with pathway and GO (Pathway/GO enrichment analysis), and correlation networks (co-expression analysis), and view the expression patterns of genes involved in specific pathways to infer the effects of the treatment. Additionally, the effectively of EXPath Tool has been performed by a case study on IAA-responsive genes. The results demonstrated that critical hub genes under IAA treatment could be efficiently identified.

1. Introduction

In recent years, with the development of next-generation sequencing (NGS) and single molecule real time sequencing (SMRT sequencing), the cost per megabase (Mb; a million bases) has been decreasing. In particular, transcriptomic sequencing has advanced the study of gene regulatory pathways in non-model organisms. In general, assembly, annotation, gene expression calculation and systematic function analysis are the four major processes after transcriptome sequencing has been finished by biotech companies. At the step of assembly, the raw reads with low quality were removed before assembling to transcripts (contigs). Numerous software for sequencing assembly have been developed such as Trinity, Oases, SOAPdenovo-Trans, Trans-Abyss and Bridger. A Linux environment is required for all of these tools. After transcripts have been assembled, the annotation is critical to define their potential function. Several tools could be applied for annotation, such as Blast2GO, BLAST, Trinotate, MAPLE and TRAPID. As well as, many database can be utilized to annotate the transcripts with regard to NCBI non-redundancy proteins (NCBI NR), UniProt Swiss-Prot/TrEMBL, Gene Ontology (GO) and KEGG pathways. Following the annotation, it is important to identify the differentially expressed transcripts when comparing different samples. Therefore, the expression value of each transcript (contig) were calculated and normalized by RPKM (reads per million per kilo base), FPKM (fragments per million per kilo base) and TPM (transcript per million) depends on sequencing methods (i.e. single end reads or paired end reads).

Finally, systematic function analysis is the most important, but trouble issue for biologists. A useful bioinformatics resource is crucial to retrieve valuable information and experimental candidates from numerous transcripts. Recently, several open source software packages (e.g. Bioconductor (based on R script) and Graphical User Interface (GUI)) and commercial software (e.g. Ingenuity IPA, CLC Workbench) have been developed for functional analysis of high-throughput transcriptomic data. However, the use of commercial software might be restricted by research budgets, while open source software might be difficult for biologists to use if they do not have a programming background. Although R script is very powerful to analyze various biological data, a series courses of R language is needed before using those open resource. For some biologists, this might be a slow remedy that cannot meet an urgency. Consequently, the processing of these huge sequencing data usually becomes a bottleneck for biologists, due to the requirement of the efficient hardware and ability of programming language.

A database named EXPath was developed by our group, which contains the experimental microarray data of Arabidopsis thaliana, Oryza sativa and Zea mays, and provides comparative expression analysis inferring metabolic pathways. It is very useful with regard to dissecting the collected microarray datasets in terms of the pathways, differential expression, co-expression and pathway/GO enrichment analysis. Unfortunately, such a valuable database cannot be used for other high-throughput expression data, and at present, the system is limited to the datasets collected in the EXPath database. In order to provide more diverse transcriptome data, EXPath Tool is developed in this work to comprehensively analyze users’ transcriptome datasets. Several analysis functions were not only advanced but also added in EXPath Tool. The online interface of EXPath Tool is especially designed for biologists, and after uploading an annotation and expression profile, five main functions (Gene Search, Pathway Search, DEGs Search, Pathways/GO Enrichment and Co-expression analysis) are provided for users to dissect their data from multiple perspectives. The concept and construction of EXPath Tool is illustrated in Supplementary Fig. S1.

2. Materials and methods

2.1. Data preprocess

Due to the heavy server loading of sequence assembly, annotation and calculation of transcript expression values, several data preprocess should be finished by users before upload their data to the EXPath Tool. Usually, these processes are a part of service of sequencing company. However, the step-by-step processing flow and required tools for data preprocess were also suggested in our system. Additionally, numerous protein sequences from NCBI RefSeq and UniprotKB for the annotation of different organisms are provided in EXPath Tool. After data preprocess, the files of annotation and expression profiles could be uploaded for further analysis in the EXPath Tool. The system flow of data preprocess is shown in Supplementary Fig. S2. It briefly describes the analysis steps for sequence assembly, annotation and expression data calculation. The details about the utility of those tools and requirements are provided in the guide page of the EXPath Tool (http://expathtool.itps.ncku.edu.tw/guide_for_custom/guide/create_project_preprocess.php).

2.2. Processing of uploaded dataset in the EXPath Tool

After users create a personal account, they can login into the EXPath Tool system and upload a set of annotation and expression profiles (call a ‘Project’). In the EXPath Tool system, the annotation files of the pathway and Gene Ontology (GO) are required, as well as the gene (transcript) expression data. More additional annotations from other databases, such as Pfam and PANTHER are encouraged. After examination of the format of the uploaded files, all data are processed individually. Figure 1 illustrates the data processing flow of uploaded files. For the expression profiles, the samples are grouped into conditions by the column name and the average expression values of each gene (transcript) under each condition are calculated. The calculation results are applied for the function of differentially expressed genes (DEGs) search (DEGs search) and co-expression analysis. For the annotation files of the pathways, if there is no assigned KEGG ortholog (KO) ID, the system automatically converts the annotated ID (i.e. NCBI GI, NCBI RefSeq, UniprotKB accession number etc.) to KO ID using the UniProt ID mapping table. The KO ID could be used to connect the KEGG pathway map, and also applied to gene annotation, pathway searches and pathway enrichment analysis. Finally, the GO term file is applied to gene annotation and GO enrichment analysis.

Figure 1

The analysis flowchart of data processing in EXPath Tool system.

2.3. Algorithm and construction

2.3.1. Pathway search

A pathway describes a series links of chemical reactions. Several critical pathways are affected by environmental changes due to various expression values of numerous genes (transcripts) under different conditions. In ‘Pathway Search’, after setting a specific pathway and the conditions of interest, the annotated pathway map will be displayed (Fig. 2a). In addition, the annotation and expression level of the KEGG ortholog can be accessed by clicking the KEGG ortholog block (Fig. 2b). The expression of the KEGG ortholog is calculated by simply averaging the expression levels of the associated genes in each condition. Furthermore, the absolute and relative expression values of various transcripts of the KEGG ortholog are graphically illustrated (Fig. 2c and d).

Figure 2

The web interface of pathway search result. (a) The genes annotated with the KEGG enzyme, the KEGG enzyme block is filled with green rectangle. A user selected gene is marked with red rectangle. The average (b), absolute (c) and relative (d) expression level of KEGG orthologs (including every transcripts) will be displayed when user clicks the green rectangle in figure (a).

2.3.2. DEGs search

To determine genes that are differentially expressed under distinct conditions, the t-test statistical method was applied by using function t.test() of the R package in EXPath Tool. Users can specify the treatment and control, set the cut-off values of fold change and P value, and the DEGs results will then be determined.

2.3.3. Co-expression analysis

Co-expressed genes are a group of genes with similar expression profiles under specific conditions. Theoretically, they may be controlled by specific transcription factors and involved in identical pathways. To determine a group of genes which are co-expressed, Pearson’s correlation coefficient and Spearman’s rank methods are applied by using function cor() of the R package in EXPath Tool. There are two analysis modes in the co-expression analysis:

(A) Correlation of a single gene

This function is identifying positive/negative correlation gene groups of a query gene. The expression pattern of co-expression gene groups is illustrated based on the z-score transformation (Supplementary Fig. S3): z denotes the z-score value, while x, and represent the raw intensity, mean and standard deviation of the gene expression levels, respectively.

(B) Correlation network

This function is identifying the correlation network of a group of query genes. The algorithm of the correlation coefficient and cut-off value should be set first. The correlation networks are constructed by Cytoscape Web based on the setting parameters (Supplementary Fig. S4).

2.3.4. KEGG pathway and GO enrichment analysis

Enrichment analysis helps researchers to find the important pathways or functional GO terms from a group of query genes. To determine the significant pathways or GO terms, the cumulative probability (p-value) of the hypergeometric distribution, which is applied by dhyper() and phyper() of the R package, are calculated to evaluate the KEGG/GO enrichment of a group of input genes. The formula is as follows:

N and M denote the number of background genes and total genes involved in specific KEGG pathways or GO terms, respectively. i is the number of input genes (defined as n) which belong to those KEGG pathways or GO terms. An example of the results is presented in Supplementary Fig. S5 in the Supplementary data.

2.3.5. Project management

After logging into EXPath Tool, users can import their data via ‘Create New Project’, and all information can be modified via project management. In order to carry out different sets of experiments, numerous projects are allowed in EXPath Tool. After finishing all the requirements for the system, the functions of Gene Search, Pathway Search, DEGs Search, Pathways/GO Enrichment and Co-expression analysis can be carried out by clicking ‘ENTER PROJECT’. It is easy to access differentially expressed genes (transcripts), critical pathways, and coexpression networks under a specific treatment by using EXPath Tool. The detailed introduction and guidelines for each function can be retrieved from the main page.

3. Results and discussion

3.1. A case study: co-expression analysis of IAA-responsive genes

EXPath tool provides a platform to assess co-expressions among a set of genes or transcripts by using users’ expression data. EXPath tool offers not only the visualization of correlation networks but also the enriched GO/Pathway of co-expressed genes. The following case study demonstrates the application of EXPath tool.

As is well known, auxins are important plant hormones involved in a wide range of cellular processes such as hypocotyl elongation, organ development and gravitropism. Nemhauser et al. revealed that 51 Arabidopsis genes were vital IAA-responsive genes by using Affymetrix GeneChips as gene expression measurement (Supplementary Table S1). Among these genes, 49 genes showed significant up-regulation under IAA treatments. The other two genes (AT1G72430 and AT4G00880) were down-regulated. All genes were related to signaling and transcription. There is a series of microarray expression data (GSE39384) retrieved from AtGeneExpress were uploaded to demonstrate the usage of EXPath Tool. The samples included three time points of IAA and mock treatment. First, 51 genes were utilized to construct correlation networks. Figure 3 demonstrates that AT1G72430 and AT4G00880 were identified the negative co-expressions (marked in green), which consisted with the expression profile in the previous study. Moreover, two IAA genes, AT1G15580 (IAA5) and AT3G23030 (IAA2) showed the highest connection with other IAA-responsive genes. Afterwards, all these genes were sent to the pathways/GO enrichment analysis. As expected, the most significant pathways and GO terms were plant hormone signal transduction, response to auxin, and auxin-activated signaling pathway (Supplementary Tables S2 and S3). This case described above reveals that EXPath Tool is an effective webserver for users to explore an important regulatory mechanism from their own high-throughput transcriptome data.

Figure 3

The correlation networks of IAA-responsive genes. Two down-regulated genes (AT1G72430 and AT4G00880) show the negative co-expressions with other genes.

3.2.The utility and significance of EXPath Tool

There are several useful webservers or software packages that have been developed to help biologists to systematically analyze their high-throughput transcriptome data. Supplementary Table S4 compares the EXPath Tool and other similar resources. BioWardrob and T-REx conduct differentially expressed analysis and co-expression analysis according to the expression levels of genes. BioWardrob can conduct additional analyses (i.e., enrichment analysis) by executing R-scripts. However, code-editing is required for this, and R-script is usually difficult for a biologist with only wet lab experience. S-MART focuses on the detection of alternative splicing and transcript isoforms, and differentially expressed genes also can be identified. Although FunRich, GeneProf and Microscope can conduct GO enrichment analysis, differential expression analysis and network construction, visual pathway and co-expression analysis are unavailable in these systems. Degust not only displays the normalized expression level of each condition, but also combines the expression level and KEGG pathway map. Unfortunately, Degust is only supplied for RNA-seq data. Additionally, wapRNA is a good tool to analyze GO terms and KEGG pathways of the DEGs. However, the significance of GO and pathways could not be easily identified in the output results, neither can the number of related transcripts. Nevertheless, all of the tools mentioned above are useful for bioinformatics scientists, while some are inefficient for biologists. EXPath Tool is constructed especially for biologists, and the web interface is very user-friendly and effective. Comprehensive analysis functions for transcriptome data (including microarray and RNA-seq) such as Gene Search, Pathway Search, DEGs Search, Pathways/GO Enrichment and Co-expression analysis are freely provided. It is thus easy to identify critical candidate genes/pathways for further experiments. This efficient tool is now available at http://expathtool.itps.ncku.edu.tw/.

Availability

The EXPath Tool webserver is freely available at http://expathtool.itps.ncku.edu.tw/.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.

Click here for additional data file.