File Dump Untuk Receiver Tanaka

 Posted admin

MySQL now supports a new collation, utf8mb4_0900_as_ci, for the utf8mb4 Unicode character set. This collation is accent sensitive and case insensitive. It is similar to the default utf8mb4 collation (utf8mb4_0900_ai_ci) except that the default collation is accent insensitive.

Firmware Tanaka T22 HD Jurasic silahkan download disini; Untuk software Hex Editor firmware receiver download disini Tanberg v2 tanaka. Receiver Tanaka T22 Hd Jurassic Baru www. Tanaka HD T22 Jurasic. Hex Editor firmware receiver download disini Tanberg v2 tanaka. SKYBOX F1+ Tandberg Powervu (by Dedy Syaputra) Download Software Tanaka T22 Jurassic HD Full Biss Key PowerVu Tanberg dan juga Tanaka T21 Download Software Firmware. Cara Upgrade Software Tanaka T-21 New Samurai Terbaru - Free download as PDF File (.pdf), Text File (.txt) or read online for free. Description: Berikut ini adalah cara untuk upgrade receiver satelit digital tanaka t-21 samurai atau t-22. TANAKA Salah satu pabrikan yang terkenal.Banyak produk-produk yang telah dikeluarkan seperti antena tv,boster,parabola,dan receiver. Dengan perkembangan teknologi sekarang ini khsususnya bagi penggemar dunia satelit parabola bahwa saat ini banyak sekali merk atau produk baru bermunculan dengan fitur terbaunya.

MySQL also now supports a new Japanese collation, utf8mb4_ja_0900_as_cs_ks, for the utf8mb4 Unicode character set. This collation is like utf8mb4_ja_0900_as_cs in that it is accent sensitive and case sensitive, but utf8mb4_ja_0900_as_cs_ks is also kana sensitive and distinguishes Katakana characters from Hiragana characters. utf8mb4_ja_0900_as_cs treats Katakana and Hiragana characters as equal for sorting. Applications that require a Japanese collation but not kana sensitivity may use utf8mb4_ja_0900_as_cs for better sort performance. utf8mb4_ja_0900_as_cs uses three weight levels for sorting; utf8mb4_ja_0900_as_cs_ks uses four.

utf8mb4_ja_0900_as_cs_ks is the first collation to use the _ks collation suffix. Japanese collations without this suffix are not kana sensitive.

For more information, see Unicode Character Sets.

Published online 2017 Jun 5. doi: 10.1186/s12864-017-3786-3
PMID: 28583129
Author informationArticle notesCopyright and License informationDisclaimer

Associated Data

Supplementary Materials
Additional file 1: ID correspondence table and information on the SRA Runs. The correspondences among Entrez Gene ID, Kegg Gene ID, and Ensemble Gene ID are shown in Table 1. Information (Run ID, Experiment ID, Study ID, sample tissue and cultivar) on the SRA Runs used to construct the gene expression matrix is shown in Table 2. (XLSX 1015 kb)
GUID: EA4D2E2A-9647-4F11-A124-F74B1958BC4A
Additional file 2: Dendrogram of the RNA-Seq samples. The dendrogram of the clustering analysis of the RNA-Seq samples. The first color bar indicates the sample tissue, and the second one indicates the cultivar. (PDF 65 kb)
GUID: 6360152A-4E48-45D8-AE80-9CB7CBEDFA90
Additional file 3: Distribution of the observed ES. The frequency distribution of ES for each KEGG pathway, which was obtained from the GSEA of all ranked gene lists. (PDF 372 kb)
GUID: F6666D28-4403-4501-A9A2-20768947F656
Additional file 4: Evaluation of the ORA of the top 1000–2500 ranked genes. ROC curves drawn from the ORA of the top (A) 1000, (B) 1500, (C) 2000, and (D) 2500 ranked genes. (PDF 132 kb)
GUID: 680C10DB-8923-4136-ACEF-CDDFA15DA979
Additional file 5: Effect of internal correlation in pathways. The average of −log10p-values (y-axis) was plotted against the internal correlation (x-axis), and the correlation coefficient between x and y, represented as r, was calculated. The ORA of the top (A) 1000, (B) 1500, (C) 2000, and (D) 2500 ranked genes. (PDF 205 kb)
GUID: 25D9ECA1-0C68-4F3E-8BF5-422E674283A0
Additional file 6: Evaluation of the ranked lists of KEGG pathways. ROC curves drawn from the ORA of the top (A) 1000, (B) 1500, (C) 2000, and (D) 2500 ranked genes. The evaluation is based on the rank order of KEGG pathways for each gene. (PDF 119 kb)
GUID: CF30CB83-ED18-4C57-A5A9-CF14A25D93E5
Additional file 7: Distribution of the p-score. The frequency distribution of p-score for each KEGG pathway. (PDF 352 kb)
GUID: 13A8810B-5138-411D-83A8-279EDFBD98ED
Additional file 8: KEGG pathway mapping of the “Carotenoid biosynthesis” pathway genes. The query gene (CRTISO gene) itself is colored purple. The intensity of the red and blue colors reflects the degree of positive and negative correlations, respectively. (PNG 61 kb)
GUID: 5B553465-41DE-4445-BF2F-A127605F28B8
Additional file 9: KEGG pathway mapping of the “Peroxisome” pathway genes. The intensity of the red and blue colors reflects the degree of positive and negative correlations with the query gene (CRTISO gene), respectively. (PNG 33 kb)
GUID: C9A0A3EE-2162-4DF0-A292-943465F7E285
Data Availability Statement

Co-expressed Pathways DataBase for Tomato is available at http://cox-path-db.kazusa.or.jp/tomato.

Abstract

Background

Gene co-expression, the similarity of gene expression profiles under various experimental conditions, has been used as an indicator of functional relationships between genes, and many co-expression databases have been developed for predicting gene functions. These databases usually provide users with a co-expression network and a list of strongly co-expressed genes for a query gene. Several of these databases also provide functional information on a set of strongly co-expressed genes (i.e., provide biological processes and pathways that are enriched in these strongly co-expressed genes), which is generally analyzed via over-representation analysis (ORA). A limitation of this approach may be that users can predict gene functions only based on the strongly co-expressed genes.

Results

In this study, we developed a new co-expression database that enables users to predict the function of tomato genes from the results of functional enrichment analyses of co-expressed genes while considering the genes that are not strongly co-expressed. To achieve this, we used the ORA approach with several thresholds to select co-expressed genes, and performed gene set enrichment analysis (GSEA) applied to a ranked list of genes ordered by the co-expression degree. We found that internal correlation in pathways affected the significance levels of the enrichment analyses. Therefore, we introduced a new measure for evaluating the relationship between the gene and pathway, termed the percentile (p)-score, which enables users to predict functionally relevant pathways without being affected by the internal correlation in pathways. In addition, we evaluated our approaches using receiver operating characteristic curves, which concluded that the p-score could improve the performance of the ORA.

Conclusions

We developed a new database, named Co-expressed Pathways DataBase for Tomato, which is available at http://cox-path-db.kazusa.or.jp/tomato. The database allows users to predict pathways that are relevant to a query gene, which would help to infer gene functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3786-3) contains supplementary material, which is available to authorized users.

Keywords: Co-expression database, Pathway, Over-representation analysis, Gene set enrichment analysis, Percentile-score

Background

Gene co-expression, the similarity of gene expression profiles under various experimental conditions, has been used as an indicator of functional relationships between genes [1], and many databases using co-expression analysis have been developed for plant research, e.g., ATTED-II [2], ALCOdb [3], AraNet v2 [4], RiceNet v2 [5], PlaNet [6], PODC [7], CoP [8], VTCdb [9], and TFGD [10]. These databases provide users with a co-expression network and a list of strongly co-expressed genes for a query gene, which has successfully contributed to the characterization of many genes [11–14].

To further facilitate the prediction of gene functions, several of these databases have also provided functional information on a set of strongly co-expressed genes with a query gene, i.e., providing biological processes and pathways that are enriched in strongly co-expressed genes [2, 4, 8, 9]. These genes are usually analyzed by over-representation analysis (ORA), which can identify biological processes and pathways enriched in the set of selected genes of interest and help to extract biological meanings, and therefore, has been used to facilitate interpretation of gene expression data [15, 16]. However, a limitation of ORA is that the results are highly dependent on the cutoff used in selecting a set of genes of interest and ignore the effect of the remaining genes [15, 16].

Gene set enrichment analysis (GSEA) has been developed to overcome the limitation of ORA [17]. Unlike ORA, GSEA can assess, without selecting genes of interest, whether biological processes and pathways are enriched at the top of a ranked list of genes ordered by the degree of differential expression [17, 18]. This enables GSEA to be performed without being dependent on the cutoff used to select differentially expressed genes. GSEA may also be effective in the case of the co-expression analysis. Namely, GSEA may be applied to a ranked list of genes ordered by the degree of co-expression, which would enable the examination of gene–pathway relationships without being dependent on the threshold used to determine strongly co-expressed genes. However, there is no co-expression database that uses GSEA in this way. Currently, only strongly co-expressed genes are considered, and therefore, users cannot predict gene functions from other co-expressed genes.

Solanum lycopersicum (tomato) is a major crop worldwide and a model system for fruit development [19]. Elucidating the metabolic functions of individual tomato genes will facilitate rational design of metabolic engineering and breeding. Tomato fruit metabolites have been intensively studied [20]. For example, the biosynthesis mechanism of lycopene, the red pigment in tomato fruits, has been well-characterized both in vitro and in vivo [21], and its consumption is reported to be associated with lowered risks of cancer and cardiovascular disease [22].

In this study, we developed a new database that allows users to predict the function of tomato genes from the results of functional enrichment analyses of co-expressed genes. Our developed database provides, for each tomato gene, a ranked list of pathways in which higher-ranked pathways are more likely related to each gene. To create the ranked pathway list, we performed ORA with several thresholds to select co-expressed genes, and applied GSEA to a ranked list of genes ordered by the co-expression degree. This approach enables users to predict pathways that are relevant to the gene of interest while considering the genes that are not strongly co-expressed. In addition, we introduced a new measure for evaluating the relationship between the gene and pathway, which improved the prediction of functionally relevant pathways.

Construction and content

We constructed a database, named Co-expressed Pathways DataBase for Tomato (CoxPathDB) [23], which aims to help users infer relevant pathways to a query gene and assist to predict its gene functions. In this section, we describe the procedural steps taken to construct the database and to evaluate our approach.

Untuk

Creation of the gene–gene correlation matrix

RNA-Seq data from tomato plants generated on the Illumina HiSeq or MiSeq platforms were downloaded from the DDBJ Sequence Read Archive (SRA) database [24]. The 1,234 downloaded SRA files were converted to FASTQ format using the fastq-dump utility of the SRA toolkit [25].

To remove low-quality reads and adapter sequences, the reads were trimmed using Trimmomatic version 0.36 [26] with the following parameters: ILLUMINACLIP:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:20:20 MINLEN:50. Then, the reads were used to estimate gene expression levels by using kallisto version 0.43.0 [27] and the tomato cDNA sequences obtained from the RefSeq database [28]. In the case of single-end reads, the average fragment length was set to 200 bp. NCBI Entrez Gene IDs were converted to Ensembl Gene IDs by using BioMart [29] and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database [30] (Additional file 1), and the genes whose IDs could not be converted were removed from the analysis. We filtered out low-quality SRA data (total estimated counts < 1 million), and then performed manual curation (e.g., removed small RNA-Seq data annotated as RNA-Seq data). Consequently, 790 SRA Runs were selected for further analysis (Additional file 1).

The expression values (transcripts per million) were quantile-normalized using the preprocessCore package in the R statistical software [31], and were log2-transformed after adding pseudo-count of 4. The 790 SRA Runs were clustered based on their gene expression profiles by the unweighted pair-group method using arithmetic averages (Additional file 2). They were clustered largely according to the sample tissues, suggesting the validity of the gene expression matrix. Then, the gene–gene correlation matrix was calculated with the gene expression matrix; correlations between gene expression profiles were calculated using the Pearson’s correlation coefficient. The gene expression matrix and the correlation matrix can be downloaded from the CoxPathDB webpage [23].

Creation of the ranked gene lists

For each tomato gene, we created a ranked list of genes based on the values of correlation coefficients in the correlation matrix; all genes except for each target gene were ordered in decreasing order of correlation with the target gene. Consequently, 13,183 ranked gene lists were created.

ORA of co-expressed genes

For each ranked gene list, we selected the top 100, 500, 1000, 1500, 2000, 2500, and 3000 ranked genes and performed ORA via the Fisher’s exact test implemented in the SciPy Python library. In ORA and GSEA (described later), we used the KEGG pathways downloaded from the KEGG database [30], because they cover a wide variety of metabolic pathways and are less redundant than Gene Ontology terms. We omitted pathways containing more than 500 genes or less than 15 genes because they might be too general or meaningless.

GSEA of the ranked gene lists

We performed GSEA [17] for all 13,183 ranked gene lists using KEGG pathways. The sample permutation approach is not applicable for this analysis, and therefore, the gene permutation approach was used to obtain significance levels. To calculate the exact p-values, we used unweighted GSEA via dynamic programming [32], which is described as follows.

Given that a target gene is t and the ranked gene list for the target gene is Lt, the overall ranked gene list is represented as

1

where n is the total number of genes in the gene list and r∗ is the correlation coefficient of gene g∗ in the ranked gene list. We also assumed that the pathway to be tested is S, the number of genes in the pathway is m, and V is a vector where V(j) is the component corresponding to gene gj in the ranked gene list Lt. V(j) takes the value 1/m for the gene in pathway S and −1/(nm) for the gene not in pathway S. The enrichment score (ES), the test statistic of GSEA, for pathway S is calculated as

2
3
4

Namely, ES is the maximum deviation of the running sum statistic, j=1kV(j), from zero. The significance level is calculated depending on whether ES(S) is positive or negative. In the case of positive ES(S), the p-value for whether pathway S is enriched at the top of the ranked gene list Lt is computed by

5

where Pr means probability, ES+ represents positive ES, and ES+null is ES+ for a randomly generated pathway in which m genes are randomly distributed in the ranked gene list Lt. This probability can be calculated exactly by using dynamic programming [32]. If the p-value is small (e.g., p< 0.05), pathway S is significant. In the case of negative ES(S), the p-value for whether pathway S is enriched at the bottom of the ranked gene list can be computed similarly. However, in this study, the p-values of pathways with negative ES were set to 1 because we focused on detecting pathways enriched at the top.

Calculation of the percentile (p)-scores

In addition to the p-value, we calculated the p-score from the GSEA results, which is defined as follows. Additional file 3 shows the observed ES distribution of each pathway, which were obtained from the GSEA of 13,183 ranked gene lists. The p-score of each pathway for a gene of interest was calculated using the observed ES distribution. To estimate the probability density function of the observed ES, we used kernel density estimation, implemented in R version 3.3.1 [31]. Based on the estimated probability density function, the p-score of each pathway for a gene of interest was calculated as

6

where ESe is the ES for the gene of interest. As described in Eq. (5), the p-value is derived from the ES distribution of randomly generated pathways, whereas the p-score is derived from the observed ES distribution that was obtained from the GSEA of all 13,183 ranked gene lists.

Evaluation of the ORA and GSEA results

We evaluated the ORA and GSEA results using receiver operating characteristic (ROC) curves. We classified gene–pathway pairs as “condition positive” or “condition negative” by using the relationship between genes and pathways in the KEGG database. If a gene was a member of a KEGG pathway, the gene–pathway pair was classified as condition positive; otherwise the gene–pathway pair was classified as condition negative. We calculated the true positive rate and false positive rate for each approach, drew ROC curves, and calculated the area under the curves (AUCs).

Results

Comparative analysis of ORA and GSEA

We evaluated the ORA of the top 100, 500, 1000, 1500, 2000, 2500, and 3000 ranked genes and the GSEA of the ranked lists of genes ordered by the co-expression degree (Fig. (Fig.11 and Additional file 4). We expected that if the p-value was small, the gene–pathway pair would be related to each other. Fig. Fig.11 shows the ROC curves drawn from the ORA of the top 100, 500, and 3000 ranked genes and the GSEA of the ranked gene lists (see Additional file 4 for the ROC curves generated from the ORA of the top 1000–2500 ranked genes). These results demonstrated that the ORA of the top 500 ranked genes performed best, with the largest AUC value of 0.782. The AUC of GSEA was smaller than that of ORA, although GSEA is a threshold free approach.

Open in a separate window

Evaluation of the ORA and GSEA. ROC curves drawn from the ORA of the top (a) 100, (b) 500, and (c) 3000 ranked genes and from (d) the GSEA of the ranked gene list

Effect of internal correlation in pathways

We examined the effect of internal correlation in pathways (Fig. (Fig.22 and Additional file 5), which is reported to lead to an overestimation of the statistical significance of GSEA [18]. To calculate the internal correlation in each pathway, we averaged the correlation coefficients between all gene pairs in each pathway. Then, to examine the effect of the internal correlation, we averaged, for each pathway and approach, the −log10p-values of gene–pathway pairs where the gene was not a member of the pathway tested. We plotted the p-value averages against the internal correlation (Fig. (Fig.22 and Additional file 5).

Open in a separate window

Effect of internal correlation in pathways. The average of −log10p-values (y-axis) was plotted against the internal correlation (x-axis), and the correlation coefficient between x and y, represented as r, was calculated. The ORA of the top (a) 100, (b) 500, and (c) 3000 ranked genes and (d) the GSEA of the ranked gene list are presented

Figure Figure22 and Additional file 5 indicate that the internal correlation in pathways affected the significance levels of ORA and GSEA, i.e., even when a gene is not a member of a pathway tested, the higher the internal correlation in the pathway, the smaller the p-value tended to become. The effect of internal correlation was larger in GSEA than in ORA (for GSEA, r= 0.775, and for ORA of the top 100–3000 ranked genes, r= 0.383, 0.577, 0.630, 0.628, 0.613, 0.598, and 0.588), which explains the lower performance of GSEA.

Evaluation of the ranked lists of KEGG pathways

CoxPathDB provides, for each tomato gene, a ranked list of KEGG pathways, where higher ranked pathways are more likely related to each gene. Therefore, we evaluated whether the highly ranked pathways are actually related to each gene (Fig. (Fig.3;3; see Additional file 6 for the ROC curves generated from the ORA of the top 1000–2500 ranked genes). Figure Figure33 and Additional file 6 show that in addition to the evaluation of the p-value (Fig. (Fig.11 and Additional file 4), the ORA of the top 500 ranked genes had the second largest AUC value of 0.758.

Open in a separate window

Evaluation of the ranked lists of KEGG pathways. ROC curves drawn from the ORA of the top (a) 100, (b) 500, and (c) 3000 ranked genes and from (d) the GSEA of the ranked gene list. The evaluation is based on the rank order of KEGG pathways for each gene

Improvement in the ranked lists of KEGG pathways

The findings highlighted in Fig. Fig.22 may be biologically meaningful; the pathways with high internal correlation may tend to be related to many genes, and therefore, had small p-value averages (Fig. (Fig.2).2). However, if all pathways tested are ordered by increasing p-values, the pathway with high internal correlation tends to be ranked higher than that with low internal correlation, which prevents users from exploring a wide variety of gene–pathway relationships. For example, the “Photosynthesis” pathway had a high internal correlation (0.577) and small p-value averages (the −log10p-values were 3.82 for GSEA and 0.560, 1.80, 2.68, 2.92, 2.89, 2.80, and 2.68 for the ORA of the top 100–3000 ranked genes). This is consistent with a previous report stating that in the Arabidopsis gene co-expression network, genes involved in photosynthesis are strongly co-expressed and over-represented in the largest co-expression module [33]. This centrality of photosynthesis genes indicates that many genes are related to photosynthesis, and therefore, the Photosynthesis pathway is often ranked high, which may hinder the discovery of relationships among other pathways and genes.

Therefore, to compare pathways without being affected by the difference in their internal correlation, we calculated p-scores from the GSEA results, as described in the “Construction and content” section. The distribution of p-scores was similar among all pathways (Additional file 7), indicating that p-scores are not affected by internal correlation and suitable for examining diverse gene–pathway relationships.

Another merit of the p-score is that it can be calculated with respect to all gene–pathway pairs. Although the ORA of the top 500 ranked genes performed well overall, these genes often do not contain any genes from the pathways to be tested. In such cases, the significance levels of the pathways cannot be compared, which decreases the performance of ORA. The p-scores can be used to order such pathways. The performance of GSEA using the p-score itself is smaller than those of ORA except for the ORA of the top 100 genes (Fig. (Fig.44a).a). However, the combination of the ORA of the top 500 ranked genes and the GSEA using the p-score had the largest AUC value among the approaches we used (Fig. (Fig.44bb).

Open in a separate window

Evaluation of the GSEA using the p-score. ROC curves drawn from (a) the GSEA using the p-score and from (b) the analysis combining the ORA of the top 500 ranked genes and the GSEA using the p-score. The evaluation is based on the rank order of KEGG pathways for each gene

Utility and discussion

In this section, we give an example of the usage of CoxPathDB (Fig. (Fig.5).5). The search box on the CoxPathDB website is shown in Fig. Fig.55a.a. In this example, the word “CRTISO” is entered. The CRTISO enzyme catalyzes the isomerization of prolycopene to lycopene [21]. After clicking the “Submit” button, the search results are displayed on the search results page (Fig. (Fig.55b).b). The “Ensembl Gene Id” column provides the links for the ranked list of KEGG pathways for the query gene. In this example, the blue link, “Solyc10g081650.1,” provides the link for the ranked list of the CRTISO gene (Fig. (Fig.55bb).

Open in a separate window

Example of the usage of CoxPathDB. a The search box on CoxPathDB. The word “CRTISO” is entered. b The search results page for the query, CRTISO. c, d The page for the CRTISO gene displaying (c) the brief gene information and (d) the ranked list of KEGG pathways

The webpage for each gene (Figs. 55cc and andd)d) also provides brief gene information from the Entrez database [34], Ensembl Plants [29] and the KEGG database [30], and the KEGG pathways that the gene belongs to (Fig. (Fig.55c).c). Each id is a link to the external database where the detailed information is available. The KEGG pathways list is displayed just below the gene information (Fig. (Fig.55d).d). In the default setting, it provides pathways that were ordered primarily by the increasing p-values obtained from the ORA of the top 500 ranked genes, and then, by p-scores calculated from GSEA, because this performed best among the approaches we used (Figs. (Figs.33 and and44 and Additional file 6). The results for pathways with less than 15 genes are available by clicking the “Show pathways with less than 15 genes” button.

The “KEGG Pathway Mapping” column provides links to the KEGG database (Fig. (Fig.55d).d). In this example, the CRTISO gene is highly correlated with “Carotenoid biosynthesis,” the pathway to which the CRTISO gene belongs (the p-value [ORA Top 500] = 0.00732 and the p-score = 0.0166). When the corresponding link in the column is clicked, the genes in the Carotenoid biosynthesis pathway are mapped to the KEGG pathway database and colored in red or blue (Additional file 8); the CRTISO gene itself is colored in purple. The intensity of red and blue colors reflects the degree of positive and negative correlations, respectively. Green color means that the corresponding gene is present in the KEGG database but is not present in CoxPathDB, whereas white color means that the corresponding gene is not present in the KEGG database.

Figure Figure55dd shows that the CRTISO gene is also highly correlated with the “Aminoacyl-tRNA biosynthesis” pathway (the p-value [ORA Top 500] = 3.25E −7 and the p-score = 0.185). Although the p-value from the ORA of the top 500 ranked genes is small, the p-score is not very low. The p-score can compare pathways without being affected by the internal correlation in pathways (Additional file 7), and therefore, the reason for the small p-value may be the high internal correlation. On the other hand, the p-score of the Carotenoid biosynthesis pathway is relatively low (0.0166), indicating that its small p-value is not caused by high internal correlation. By checking the p-score, users can examine gene–pathway relationships while considering the context of pathways.

The “Peroxisome” pathway also exhibits a low p-value for the ORA of the top 500 ranked genes (0.0160) and a low p-score (0.0120). It has been proposed that antioxidative enzymes in peroxisomes may act as modulators of reactive oxygen species (ROS) signaling during pepper fruit maturation [35]. Solyc07g063430.2 (Entrez Gene ID: 101247444), encoding an MPV17 protein which may be involved in ROS metabolism, is strongly co-expressed with the CRTISO gene (Additional file 9), which may suggest that peroxisomal ROS generated via this protein modulates the lycopene biosynthesis during tomato fruit maturation.

Conclusions

In this study, we developed a database named Co-expressed Pathways DataBase for Tomato [23]. The database provides, for each tomato gene, a ranked list of KEGG pathways, where the higher-ranked pathways are more likely related to each gene. The p-score enables users to predict functionally relevant pathways without being affected by internal correlation in pathways.

Additional files

Additional file 1(1014K, xlsx)

ID correspondence table and information on the SRA Runs. The correspondences among Entrez Gene ID, Kegg Gene ID, and Ensemble Gene ID are shown in Table 1. Information (Run ID, Experiment ID, Study ID, sample tissue and cultivar) on the SRA Runs used to construct the gene expression matrix is shown in Table 2. (XLSX 1015 kb)

Sambandar thevaram in tamil pdf. Apr 7, 2016 - TamilBrahmins.com is an online forum for the Tamil Brahmin Community spread across. Here is the PDF of Thirugnana sambandar thevaram.

Additional file 2(64K, pdf)

Dendrogram of the RNA-Seq samples. The dendrogram of the clustering analysis of the RNA-Seq samples. The first color bar indicates the sample tissue, and the second one indicates the cultivar. (PDF 65 kb)

Additional file 3(371K, pdf)

Distribution of the observed ES. The frequency distribution of ES for each KEGG pathway, which was obtained from the GSEA of all ranked gene lists. (PDF 372 kb)

Additional file 4(131K, pdf)

Evaluation of the ORA of the top 1000–2500 ranked genes. ROC curves drawn from the ORA of the top (A) 1000, (B) 1500, (C) 2000, and (D) 2500 ranked genes. (PDF 132 kb)

Additional file 5(204K, pdf)

Effect of internal correlation in pathways. The average of −log10p-values (y-axis) was plotted against the internal correlation (x-axis), and the correlation coefficient between x and y, represented as r, was calculated. The ORA of the top (A) 1000, (B) 1500, (C) 2000, and (D) 2500 ranked genes. (PDF 205 kb)

Additional file 6(118K, pdf)

Evaluation of the ranked lists of KEGG pathways. ROC curves drawn from the ORA of the top (A) 1000, (B) 1500, (C) 2000, and (D) 2500 ranked genes. The evaluation is based on the rank order of KEGG pathways for each gene. (PDF 119 kb)

Additional file 7(352K, pdf)

Distribution of the p-score. The frequency distribution of p-score for each KEGG pathway. (PDF 352 kb)

Additional file 8(61K, png)

KEGG pathway mapping of the “Carotenoid biosynthesis” pathway genes. The query gene (CRTISO gene) itself is colored purple. The intensity of the red and blue colors reflects the degree of positive and negative correlations, respectively. (PNG 61 kb)

Additional file 9(33K, png)

KEGG pathway mapping of the “Peroxisome” pathway genes. The intensity of the red and blue colors reflects the degree of positive and negative correlations with the query gene (CRTISO gene), respectively. (PNG 33 kb)

Acknowledgements

Computations were partially performed on the NIG supercomputer at ROIS National Institute of Genetics.

Funding

This work was conducted as part of the Low Carbon Technology Research and Development program funded by the Ministry of the Environment, Japan.

Availability of data and materials

Co-expressed Pathways DataBase for Tomato is available at http://cox-path-db.kazusa.or.jp/tomato.

Authors’ contributions

File Dump Untuk Receiver Tanaka

TN designed the study, analyzed the data, developed the database, and wrote the paper. NS gave technical support, supervised the database development, and complemented the writing. TO, HO, and DS designed the study, and supervised and complemented the writing. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Dongle Wifi Untuk Receiver

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Abbreviations

AUCArea under the curve
CoxPathDBCo-expressed Pathways DataBase for Tomato
GSEAGene set enrichment analysis
KEGGKyoto Encyclopedia of Genes and Genomes
ORAOver-representation analysis
p-scorepercentile-score
ROCReceiver operating characteristic
ROSReactive oxygen species
SRASequence Read Archive

Footnotes

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3786-3) contains supplementary material, which is available to authorized users.

Contributor Information

Takafumi Narise, Email: pj.ro.asuzak@esiran.

Nozomu Sakurai, Email: pj.ro.asuzak@iarukas.

Takeshi Obayashi, Email: pj.ca.ukohot.iece@ihsayabo.

Hiroyuki Ohta, Email: pj.ca.hcetit.oib@athoh.

Daisuke Shibata, Email: pj.ro.asuzak@atabihs.

References

1. Usadel B, Obayashi T, Mutwil M, Giorgi FM, Bassel GW, Tanimoto M, Chow A, Steinhauser D, Persson S, Provart NJ. Co-expression tools for plant biology: opportunities for hypothesis generation and caveats. Plant Cell Environ. 2009;32(12):1633–51. doi: 10.1111/j.1365-3040.2009.02040.x. [PubMed] [CrossRef] [Google Scholar]
2. Aoki Y, Okamura Y, Tadaka S, Kinoshita K, Obayashi T. Atted-II in 2016: a plant coexpression database towards lineage-specific coexpression. Plant Cell Physiol. 2016;57(1):e5. doi: 10.1093/pcp/pcv165.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
3. Aoki Y, Okamura Y, Ohta H, Kinoshita K, Obayashi T. Alcodb: Gene coexpression database for microalgae. Plant Cell Physiol. 2016;57(1):e3. doi: 10.1093/pcp/pcv190.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
4. Lee T, Yang S, Kim E, Ko Y, Hwang S, Shin J, Shim JE, Shim H, Kim H, Kim C, et al. AraNet v2: an improved database of co-functional gene networks for the study of Arabidopsis thaliana and 27 other nonmodel plant species. Nucleic Acids Res. 2015;43(D1):D996–D1002. doi: 10.1093/nar/gku1053.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
5. Lee T, Oh T, Yang S, Shin J, Hwang S, Kim CY, Kim H, Shim H, Shim JE, Ronald PC, et al. RiceNet v2: an improved network prioritization server for rice genes. Nucleic Acids Res. 2015;43(W1):W122–W127. doi: 10.1093/nar/gkv253.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
6. Mutwil M, Klie S, Tohge T, Giorgi FM, Wilkins O, Campbell MM, Fernie AR, Usadel B, Nikoloski Z, Persson S. PlaNet: combined sequence and expression comparisons across plant networks derived from seven species. The Plant Cell. 2011;23(3):895–910. doi: 10.1105/tpc.111.083667.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
7. Ohyanagi H, Takano T, Terashima S, Kobayashi M, Kanno M, Morimoto K, Kanegae H, Sasaki Y, Saito M, Asano S, et al. Plant Omics Data Center: an integrated web repository for interspecies gene expression networks with NLP-based curation. Plant Cell Physiol. 2015;56(1):e9. doi: 10.1093/pcp/pcu188.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
8. Ogata Y, Suzuki H, Sakurai N, Shibata D. Cop: a database for characterizing co-expressed gene modules with biological information in plants. Bioinformatics. 2010;26(9):1267–8. doi: 10.1093/bioinformatics/btq121. [PubMed] [CrossRef] [Google Scholar]
9. Wong DC, Sweetman C, Drew DP, Ford CM. VTCdb: a gene co-expression database for the crop species Vitis vinifera (grapevine) BMC Genomics. 2013;14:882. doi: 10.1186/1471-2164-14-882.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
10. Fei Z, Joung JG, Tang X, Zheng Y, Huang M, Lee JM, McQuinn R, Tieman DM, Alba R, Klee HJ, et al. Tomato Functional Genomics Database: a comprehensive resource and analysis package for tomato functional genomics. Nucleic Acids Res. 2011;39(Suppl 1):D1156–D1163. doi: 10.1093/nar/gkq991.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
11. Hirai MY, Sugiyama K, Sawada Y, Tohge T, Obayashi T, Suzuki A, Araki R, Sakurai N, Suzuki H, Aoki K, et al. Omics-based identification of Arabidopsis Myb transcription factors regulating aliphatic glucosinolate biosynthesis. Proc Natl Acad Sci USA. 2007;104(15):6478–83. doi: 10.1073/pnas.0611629104.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
12. Okazaki Y, Shimojima M, Sawada Y, Toyooka K, Narisawa T, Mochida K, Tanaka H, Matsuda F, Hirai A, Hirai MY, et al. A chloroplastic UDP-glucose pyrophosphorylase from Arabidopsis is the committed enzyme for the first step of sulfolipid biosynthesis. The Plant Cell. 2009;21(3):892–909. doi: 10.1105/tpc.108.063925.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
13. Sugano SS, Shimada T, Imai Y, Okawa K, Tamai A, Mori M, Hara-Nishimura I. Stomagen positively regulates stomatal density in Arabidopsis. Nature. 2010;463(7278):241–4. doi: 10.1038/nature08682. [PubMed] [CrossRef] [Google Scholar]
14. Block A, Widhalm JR, Fatihi A, Cahoon RE, Wamboldt Y, Elowsky C, Mackenzie SA, Cahoon EB, Chapple C, Dudareva N, et al. The origin and biosynthesis of the benzenoid moiety of ubiquinone (coenzyme Q) in Arabidopsis. The Plant Cell. 2014;26(5):1938–48. doi: 10.1105/tpc.114.125807.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
15. Khatri P, Sirota M, Butte AJ. Ten years of pathway analysis: current approaches and outstanding challenges. PLoS Comput Biol. 2012;8(2):e1002375. doi: 10.1371/journal.pcbi.1002375.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
16. Irizarry RA, Wang C, Zhou Y, Speed TP. Gene set enrichment analysis made simple. Stat Methods Med Res. 2009;18(6):565–75. doi: 10.1177/0962280209351908.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
17. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al.Proc Natl Acad Sci USA. 2005; 102(43):15545–50. [PMC free article] [PubMed]
18. Nam D, Kim SY. Gene-set approach for expression pattern analysis. Brief Bioinform. 2008;9(3):189–97. doi: 10.1093/bib/bbn001. [PubMed] [CrossRef] [Google Scholar]
19. Tomato Genome Consortium: The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012; 485(7400):635–41. [PMC free article] [PubMed]
20. Seymour GB, Chapman NH, Chew BL, Rose JK. Regulation of ripening and opportunities for control in tomato and other fruits. Plant Biotechnol J. 2013;11(3):269–78. doi: 10.1111/j.1467-7652.2012.00738.x. [PubMed] [CrossRef] [Google Scholar]
21. Fantini E, Falcone G, Frusciante S, Giliberto L, Giuliano G. Dissection of tomato lycopene biosynthesis through virus-induced gene silencing. Plant Physiol. 2013;163(2):986–98. doi: 10.1104/pp.113.224733.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
22. Rao AV, Agarwal S. Role of antioxidant lycopene in cancer and heart disease. J Am Coll Nutr. 2000;19(5):563–9. doi: 10.1080/07315724.2000.10718953. [PubMed] [CrossRef] [Google Scholar]
23. Co-expressed Pathways DataBase for Tomato. http://cox-path-db.kazusa.or.jp/tomato/. [PMC free article] [PubMed]
24. Kodama Y, Mashima J, Kaminuma E, Gojobori T, Ogasawara O, Takagi T, Okubo K, Nakamura Y. The DNA Data Bank of Japan launches a new resource, the DDBJ Omics Archive of functional genomics experiments. Nucleic Acids Res. 2012;40(D1):D38–D42. doi: 10.1093/nar/gkr994.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
25. Leinonen R, Sugawara H, Shumway M. The sequence read archive. Nucleic Acids Res. 2011;39(Suppl 1):D19–D21. doi: 10.1093/nar/gkq1019.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
26. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. doi: 10.1093/bioinformatics/btu170.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
27. Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016;34(5):525–7. doi: 10.1038/nbt.3519. [PubMed] [CrossRef] [Google Scholar]
28. O’Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, Rajput B, Robbertse B, Smith-White B, Ako-Adjei D, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016;44(D1):D733–D745. doi: 10.1093/nar/gkv1189.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
29. Bolser D, Staines DM, Pritchard E, Kersey P. Ensembl Plants: integrating tools for visualizing, mining, and analyzing plant genomics data. Methods Mol Biol. 2016;1374:115–40. doi: 10.1007/978-1-4939-3167-5_6. [PubMed] [CrossRef] [Google Scholar]
30. Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44(D1):D457–D462. doi: 10.1093/nar/gkv1070.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
31. R Core Team. R: A Language and Environment for Statistical Computing. 2016. http://www.r-project.org/.
32. Keller A, Backes C, Lenhof HP. Computation of significance scores of unweighted gene set enrichment analyses. BMC Bioinforma. 2007;8:290. doi: 10.1186/1471-2105-8-290.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
33. Mao L, Van Hemert JL, Dash S, Dickerson JA. Arabidopsis gene co-expression network and its functional modules. BMC Bioinforma. 2009;10:346. doi: 10.1186/1471-2105-10-346.[PMC free article] [PubMed] [CrossRef] [Google Scholar]
34. NCBI Resource Coordinators: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2016; 44(D1):D7–D19. [PMC free article] [PubMed]
35. Mateos RM, León AM, Sandalio LM, Gómez M, del Río LA, Palma JM. Peroxisomes from pepper fruits (Capsicum annuum L.): purification, characterisation and antioxidant activity. J Plant Physiol. 2003;160(12):1507–16. doi: 10.1078/0176-1617-01008. [PubMed] [CrossRef] [Google Scholar]
Articles from BMC Genomics are provided here courtesy of BioMed Central