- Open Access
GenCLiP: a software program for clustering gene lists by literature profiling and constructing gene co-occurrence networks related to custom keywords
BMC Bioinformaticsvolume 9, Article number: 308 (2008)
Biomedical researchers often want to explore pathogenesis and pathways regulated by abnormally expressed genes, such as those identified by microarray analyses. Literature mining is an important way to assist in this task. Many literature mining tools are now available. However, few of them allows the user to make manual adjustments to zero in on what he/she wants to know in particular.
We present our software program, GenCLiP (Gene Cluster with Literature Profiles), which is based on the methods presented by Chaussabel and Sher (Genome Biol 2002, 3(10):RESEARCH0055) that search gene lists to identify functional clusters of genes based on up-to-date literature profiling. Four features were added to this previously described method: the ability to 1) manually curate keywords extracted from the literature, 2) search genes and gene co-occurrence networks related to custom keywords, 3) compare analyzed gene results with negative and positive controls generated by GenCLiP, and 4) calculate probabilities that the resulting genes and gene networks are randomly related. In this paper, we show with a set of differentially expressed genes between keloids and normal control, how implementation of functions in GenCLiP successfully identified keywords related to the pathogenesis of keloids and unknown gene pathways involved in the pathogenesis of keloids.
With regard to the identification of disease-susceptibility genes, GenCLiP allows one to quickly acquire a primary pathogenesis profile and identify pathways involving abnormally expressed genes not previously associated with the disease.
Biomedical researchers often want to explore pathogenesis and pathways regulated by abnormally expressed genes, such as those identified by microarray analyses. Literature mining is an important way to assist in this task . Many literature mining tools or methods have been developed [2–14] that can extract keywords or gene networks from the literature to functionally group genes and visualize their relationships. However, to our knowledge none of these tools are capable of automatically and easily constructing gene networks among the analyzed genes based on specified keywords, an aspect that would be helpful for investigating disease-associated signaling pathways. Here we report our software program, GenCLiP (Gene Cluster with Literature Profiles), which is based on the two methods provided by Chaussabel and Sher  and Jenssen et al. . GenCLiP can cluster functionally related genes based on up-to-date literature profiling and identify gene networks based on specified keywords.
The method of Chaussabel and Sher  can be used to analyze gene lists in order to cluster the genes based on up-to-date literature associations. This method first defines terms that occur frequently and exclusively in at least two analyzed genes' up-to-date related literature as keywords, then clusters the analyzed genes based on keyword occurrences. We have added a new feature to this method that allows the user to manually curate keywords through a series of operations including adding/removing, term weighting, defining synonyms, or defining singular/plural forms. Thus, the processed keywords will be more related to the pathogenesis of the specified disease.
The method of Jenssen et al.  shows how gene networks can be constructed based on simple co-occurrence. We have added the capacity to limit the extraction of analyzed gene co-occurrences from the literature to only those found in literature that contain certain keywords. This enables the user to readily identify pathways that are specifically associated with the pathogenesis of a particular disease by selecting suitable keywords.
To verify the literature profiling cluster result and the networks related to specified keywords, we have added two features to GenCLiP. One is that GenCLiP can generate two groups of control genes: 1) by using the full gene set to generate a group of negative control genes randomly, and 2) by using the full gene set and the curated keywords to search the database (e.g. PubMed  or Entrez Gene ) for all related genes to generate a group of positive control genes randomly. The other feature is that GenCLiP can use PubMed search results of the full gene set to perform 10,000 or more random simulations to calculate the probability that a set of randomly picked genes contain the same or more number of related genes (or gene pairs) as the analyzed genes do.
The workflow of GenCLiP is shown in Figure 1. First, a group of positive control genes and a group of negative control genes are generated based on the imported gene list. Then the literature pertaining to each gene of the three groups of genes is retrieved from PubMed. After that, the keywords related to each group of genes are auto-extracted from the literature. The keywords can be manually curated. Then, each group of genes is clustered with the keywords, and the gene co-occurrence networks can be constructed among each group of genes based on certain keywords. After that, the cluster results and the gene networks should be compared among the three groups of genes, and the user can select the keywords that are more related to the positive control genes and the analyzed genes compared to the negative control genes to construct gene networks. Once the analyzed genes and the positive control genes are found to contain more genes (or more complex gene networks) related to certain keywords than the negative control genes, 10,000 or more random simulations are done to decide whether it occurs randomly. Thus, an inference can be obtained for further experiment verification.
Generation of controls
To generate the negative control genes, the full gene set from which the analyzed genes are derived is used to generate a group of genes randomly. To generate the positive control genes, the full gene set and certain keywords are first used to search the database (e.g. PubMed or Entrez Gene) for all known genes related to the specified keywords. Then, the known-related gene set is used to generate a group of genes randomly. The analyzed genes, the positive control genes, and the negative control genes should have the same number of genes. And the average number of literature per gene for each of them should be comparative. We arbitrarily set that the average number of literature per gene for each group of randomly picked genes should be between 75 and 125% of that for the analyzed genes or as close as possible to that for the analyzed genes. A file recorded the literature number for each of all human genes has been generated and is updated at approximately 6-month intervals. This file is used to calculate the average number of literature per gene for each of the three groups of genes.
To retrieve literature pertaining to each gene of the three groups of genes, the NCBI EUtilities , ESearch, and EFetch web services, are used to access the PubMed database  for the description . The user can decide whether to provide the gene symbol directly, or provide the gene ID (HUGO , Entrez , or Unigene ). If an alternative gene ID is available, it can be converted to the appropriate input form using MatchMiner [20, 21]. Each gene's literature is saved in a text file with the gene's official symbol as the file name. To solve the ambiguity of gene names [22, 23], including synonyms (different names for the same gene) and homonyms (different genes or unrelated concepts with the same name), a human gene thesaurus that collected all the aliases for each gene name from the HUGO Nomenclature Committee database  and the Entrez Gene database  is used. The specificity of each gene name in the thesaurus has been improved by several methods for the description [2, 3, 5, 6]. The thesaurus is discussed in detail in additional files [see Additional file 1].
Auto-extraction of keywords
Auto-extraction of keywords is performed for the description . Briefly, terms are first extracted from literature titles and abstracts, and their occurrences (number of literature containing a given term divided by the number of total literature) for each gene are calculated. The terms are then filtered systematically using several criteria. (i) Terms with a baseline occurrence (the average occurrence of a set of 250 randomly picked genes, which has been proven unbiased ) of more than 5% are eliminated. (ii) Term-occurrence values for each gene are compared to the baseline. The difference cut-off between gene term occurrence and baseline occurrence is set as: cut-off = t + (k/n) (where n is the number of abstracts for a given gene, t and k are constants, with defaults t = 15% and k = 1.5). It is noted that if one sets a lower value of t and k, then more terms will pass the filter, but the noise will increase . (iii) Only the terms that pass through the filter for at least two of the analyzed genes are retained. These retained terms are considered keywords.
Manual curation of keywords
Keywords can be manually curated. The user can remove unrelated keywords and add relevant keywords (single terms or phrases). The user sets the weight for certain keywords that are perceived more important than others. The user defines certain keywords as one synonym entity. The user also determines which keywords have singular/plural forms.
Clustering analysis is performed for the description . Briefly, occurrences of all keywords for each gene pass through the following processes: (i) The occurrence of each keyword in its singular and plural forms is averaged into one unique occurrence; (ii) Each occurrence is multiplied by its weight; (iii) Occurrences of synonyms are averaged into one unique occurrence, and each synonym entity is represented by a keyword. An array file is then generated and used to do clustering analysis with the average linkage hierarchical clustering algorithm for the description . This file can also be used for clustering analysis with publicly available software, such as Cluster 3.0  and SpotFire (Göteborg, Sweden).
Gene co-occurrences are searched from the literature that contains certain keywords. The Neato program in the WinGraphviz software  is used to create a two-dimensional layout.
Random simulation is performed with two steps. First, each gene of the full gene set is used to search PubMed for whether its literature mentions certain keywords, and the resulting PubMed IDs are recorded. Second, for each simulation, the same number of genes as the number of the analyzed genes are randomly picked from the full gene set, and the number of genes (and then gene pairs, i.e. two genes sharing the same PubMed ID) related to the specified keywords are counted. The average number of literature per gene for the randomly picked genes should be comparative with that for the analyzed genes (for details, see the "Generation of controls" section). After 10,000 or more random simulations, if the distribution of the number of related genes (or gene pairs) is similar to the expected normal distribution and the probability that a set of randomly picked genes contain the same or more number of related genes (or gene pairs) as the analyzed genes do is less than 0.05 (i.e. P < 0.05), then it can be inferred that the gene relatedness is not random. [27, 28].
The literature containing certain genes and keywords can be searched and displayed with the genes and the keywords coded by different colors.
GenCLiP was used to analyze a list of 247 differentially expressed genes [see Additional file 2] between keloids and normal skin derived from a sample gene-expression dataset generated by a microarray (CSC-GE, Shenzhen Chipscreen Biosciences Ltd., China) . A keloid is a type of scar that results in an overgrowth of tissue at the site of a healed skin injury .
The list of 247 genes were transformed into 234 unique official symbols (13 genes were removed because 10 did not have an official symbol and the other 3 were repeated genes). Two hundred thirty-two of the 234 genes had related literature with an average of 848 literature per gene. The full gene set [see Additional file 3] of the microarray from which the analyzed genes had been derived was used to generate a list of 232 genes randomly [see Additional file 4] with an average of 887 literature per gene. And the full gene set of all human genes was used to search PubMed for all known genes related to the term "keloids". A group of 232 positive control genes [see Additional file 5] with an average of 2,451 literature per gene was then generated from the full known keloids-related gene set. The literature pertaining to each gene of the three groups of genes was then retrieved to the local machine with one text file per gene and one file folder per group of genes. After that, the keywords were auto-extracted from the literature with 502, 505, and 591 keywords for the analyzed genes, the positive control genes, and the negative control genes, respectively [see Additional files 6, 7 and 8] (Table 1).
The three groups of genes and their keywords were used to do cluster analysis. The cluster results showed that each of the three keyword lists was divided into many sub-groups (data not shown). Some of these sub-groups, such as sub-group containing keywords "fibroblast" or "collagen", are directly related to the biology of keloids. However, most of these sub-groups, such as sub-groups containing keywords '15', 'kda', and 'mg', are obviously un-related to the biology of keloids [see also Additional files 6, 7 and 8]. The cluster results were compared with each other and showed that there were more keywords related to the biology of keloids in the group of analyzed genes and the group of positive control genes than that in the group of negative control genes (data not shown).
The 502 keywords of the analyzed genes were manually curated further. First, 476 of the 502 keywords that did not seem directly related to keloid biology were manually removed and only 26 keywords remained (Table 2). Furthermore, some of the terms such as 'keloid' and 'scar' that are closely related to keloid pathogenesis were not included as keywords due to low occurrence frequency. Sixteen of them were manually added as keywords (Table 2, marked by asterisks), yielding 41 total keywords. Since some keywords such as 'keloid' and 'scar' might be more related to keloid pathogenesis than the other keywords, higher weights were given to them (Table 2). Therefore, after cluster analysis, genes related to them will be easier to group together. For example, to group keloid-related genes, the weight of the keyword 'keloid' was set from low to high until they can be grouped together (Table 2). Many of the 41 keywords are synonyms, such as 'hypoxia' and 'hypoxic', so the 41 keywords were reset as 16 distinct keyword entities (Table 2). Lastly, some of the 41 keywords might appear in the literature in singular and plural forms (e.g., 'keloid' and 'keloids'), and these should be considered as one keyword.
The 41 curated keywords were used to re-cluster the three groups of genes (Figure 2 and Table 1). The cluster results showed that the positive control genes (Figures 2A, 2D) and the analyzed genes (Figures 2B, 2E) were more highly related to the keywords, whereas the negative control genes (Figures 2C, 2F) were only related to some of the keywords to different degrees. Comparison of the cluster results showed that the keywords 'scar' and 'collagen' were specific to the analyzed genes and the positive control genes, and in a relative sense the keywords 'fibroblast', 'hypoxia', 'angiogenesis', and 'invasion' were more related to the positive control genes and the analyzed genes compared to the negative control genes. These results indicated that the curated keywords and the analyzed genes were indeed related to keloid disease, and the pathogenesis of the analyzed keloid patients might be related to the abnormalities of scars, collagen, fibroblasts, hypoxia, angiogenesis and invasion.
The co-occurrence networks were constructed among each of the three groups of genes based on the curated keyword "keloid" (Table 1). As a result, there were 25 known keloid-related genes reported among the analyzed genes, whereas two pairs of them (COL14A1 and TNC, COL1A2 and PEPD) were co-occurred in the literature that mentioned the curated keyword 'keloid' (data not shown). To determine whether the 25 known keloid-related genes were identified randomly, first, each gene of the full gene set of the microarray was used to search PubMed with the curated keyword "keloid". Then, 100,000 random simulations were done. As a result, the distribution of the number of keloid-related genes was found to be similar to the normal distribution, and the probability that a set of 232 randomly picked genes contained more than 24 keloid-related genes was P = 0.00003 (data not shown). Thus, the 232 analyzed genes are related to keloid disease.
Exposure of fibroblasts to hypoxia is involved in keloid pathogenesis . Thus, because only a few pairs of genes co-occurred in the literature that mentioned the keyword 'keloid', and because the above clustering analysis showed that the pathogenesis of the analyzed keloid patients might be related to the abnormalities of fibroblasts and hypoxia, gene networks among the three groups of genes were further constructed based on the keywords 'hypoxia' and 'fibroblast' (Table 1 and Figure 3), where genes known-related to the keyword 'keloid' were indicated too. The results showed that there were 31 of the 232 analyzed genes related to the keywords 'hypoxia' and 'fibroblast', while 19 of them (containing 20 pairs of genes) formed gene co-occurrence networks (Figure 3B). The probability that a set of randomly picked genes contain more than 30 genes (or 19 gene pairs) related to the keywords 'hypoxia' and 'fibroblast' was then calculated as P = 0.0228 (or P = 0.01589), and the distribution of the number of related gene and gene pairs was similar to the normal distribution (Figure 3D). Thus, the keywords 'fibroblast' and 'hypoxia' are related to keloid disease and the analyzed genes are engaged in the same networks related to the keywords 'hypoxia' and 'fibroblasts'.
Since the above analysis showed that the two keloid-related genes, HIF1A and MMP2, co-occurred in the literature that mentioned the keywords 'hypoxia' and 'fibroblast' but not in the literature that mentioned the keyword 'keloid', HIF1A and MMP2 might be involved in an unknown pathway that is abnormal in keloids. The literature that mentioned the keywords 'hypoxia' and 'fibroblast', and also the genes HIF1A and MMP2 was thus further searched and displayed. There were three literature obtained (data not shown). Two of the three literature showed that the expression of matrix metallopeptidase 2 (encoded by MMP2), which is involved in the breakdown of extracellular matrix in normal physiological processes, is up-regulated by hypoxia through over-expression of hypoxia-inducible factor HIF1A in fibroblasts [30, 31]. Thus, the pathway (MMP2 up-regulated by HIF1A) may be involved in the expansion of keloids and merit further experimental study.
Here we show that using GenCLiP, the researcher can gain the pathogenesis profiling of a specified disease and unknown pathways involved in the pathogenesis of the disease from a set of related genes. However, it should be noted that some of the pathways inferred by GenCLiP may already be well established, despite genes never having been mentioned together in any abstract . The user can check them in the curated databases (refer to Pathguide: the pathway resource list ) or read the full-text of related papers. Thus, it will be up to the user to decide whether some of the inferred pathways are trivial and thus unworthy of future investigation.
A obvious advantage of GenCLiP over other tools is that all of the associations between the analyzed human genes and the keywords outputted by GenCLiP are up-to-date, because they are calculated de novo. However, the cost is that the processes such as literature retrieval and keywords auto-extraction are time consuming. Generally, the cycle for GenCLiP analysis of a set of experimental data will be several days, especially in developing countries with slow internet speeds. While many of the existing publicly available literature mining tools will provide results to the user immediately. Because they use pre-calculated associations between all human genes and keywords, which are saved in a database that thus allows the user to search the database and obtain the results immediately.
Some of the existing publicly available tools such as PubGene  and Ali Baba  can also construct gene networks based on keywords. However, the essential differences between GenCLiP and these tools are that: (i) GenCLiP constructs gene networks from the list of analyzed genes, whereas PubGene and Ali Baba construct gene networks from all human genes and only allow the user to query sub-networks containing one or several of the analyzed genes – making the resulting networks quite unrelated to the list of analyzed genes; (ii) for keywords, GenCLiP can use any of the terms or phrases present in abstracts, whereas PubGene and Ali Baba can only use the terms from dictionaries such as MESH, Gene Ontology, and diseases, which may sometimes limit the capability of the user to construct gene networks of interest; and (iii) GenCLiP can use keywords in combination to construct gene networks, whereas PubGene and Ali Baba can only use one keyword to construct a gene network, which limits the capability of constructing gene networks that are highly specific.
In the future, we will develop GenCLiP as a web-based tool. And we will add several features to the next version of GenCLiP, such that: (i) the resulting gene networks will also show the expression value of each node (gene); and (ii) integrating GenCLiP with other data mining tools that can mine pathways from a large set of high-throughput data (such as microarray gene expression data) and the curated pathway databases (such as databases listed in Pathguide ) will form a platform [33–36] that can explore the pathogenesis of a specified disease more comprehensively and powerfully.
In this paper we present an program GenCLiP, a literature mining tool that can cluster a list of genes with keywords that are auto-extracted from their up-to-date related literature and then manually curated by the user. GenCLiP can also generate a group of negative control genes and a group of positive control genes for comparing the cluster results with the analyzed genes to filter out un-specific keywords. GenCLiP allows users to search gene and gene co-occurrence networks related to certain keywords among each of the three groups of genes, and decide whether their associations are random by using a random simulation. Further, GenCLiP can display literature mentioning specified genes and keywords for manual verification of their associations. All these features empower GenCLiP to interpret disease pathogenesis and find novel genes or pathways for further research.
Availability and requirements
GenCLiP is freely available under the BSD Open Source license for download from http://www.genclip.com/ or as additional files to this manuscript [see Additional file 9]. It runs on Windows platform. It is noted that a license is needed to include source code from the GenCLiP Package in commercial software projects.
Jensen LJ, Saric J, Bork P: Literature mining for the biologist: from information retrieval to biological discovery. Nat Rev Genet 2006, 7(2):119–129. 10.1038/nrg1768
Chaussabel D, Sher A: Mining microarray expression data by literature profiling. Genome Biol 2002, 3(10):RESEARCH0055. 10.1186/gb-2002-3-10-research0055
Jenssen TK, Laegreid A, Komorowski J, Hovig E: A literature network of human genes for high-throughput analysis of gene expression. Nat Genet 2001, 28(1):21–28. 10.1038/88213
Hoffmann R, Valencia A: A gene network for navigating the literature. Nat Genet 2004, 36(7):664. 10.1038/ng0704-664
Alako BT, Veldhoven A, van Baal S, Jelier R, Verhoeven S, Rullmann T, Polman J, Jenster G: CoPub Mapper: mining MEDLINE based on search term co-publication. BMC Bioinformatics 2005, 6(1):51. 10.1186/1471-2105-6-51
Rubinstein R, Simon I: MILANO--custom annotation of microarray results using automatic literature searches. BMC Bioinformatics 2005, 6(1):12. 10.1186/1471-2105-6-12
Plake C, Schiemann T, Pankalla M, Hakenberg J, Leser U: AliBaba: PubMed as a graph. Bioinformatics 2006, 22(19):2444–2445. 10.1093/bioinformatics/btl408
Li S, Wu L, Zhang Z: Constructing biological networks through combined literature mining and microarray analysis: a LMMA approach. Bioinformatics 2006, 22(17):2143–2150. 10.1093/bioinformatics/btl363
Daraselia N, Yuryev A, Egorov S, Mazo I, Ispolatov I: Automatic extraction of gene ontology annotation and its correlation with clusters in protein networks. BMC Bioinformatics 2007, 8(1):243. 10.1186/1471-2105-8-243
Gajendran VK, Lin JR, Fyhrie DP: An application of bioinformatics and text mining to the discovery of novel genes related to bone biology. Bone 2007, 40(5):1378–1388. 10.1016/j.bone.2006.12.067
Hammamieh R, Chakraborty N, Wang Y, Laing M, Liu Z, Mulligan J, Jett M: GeneCite: A Stand-alone Open Source Tool for High-Throughput Literature and Pathway Mining. Omics 2007, 11(2):143–151. 10.1089/omi.2007.4322
Rodriguez-Penagos C, Salgado H, Martinez-Flores I, Collado-Vides J: Automatic reconstruction of a bacterial regulatory network using Natural Language Processing. BMC Bioinformatics 2007, 8: 293. 10.1186/1471-2105-8-293
Kim H, Park H, Drake BL: Extracting unrecognized gene relationships from the biomedical literature via matrix factorizations. BMC Bioinformatics 2007, 8 Suppl 9: S6. 10.1186/1471-2105-8-S9-S6
Bundschus M, Dejori M, Stetter M, Tresp V, Kriegel HP: Extraction of semantic biomedical relations from text using conditional random fields. BMC Bioinformatics 2008, 9(1):207. 10.1186/1471-2105-9-207
HUGO Nomenclature Committee[http://www.genenames.org/]
Bussey KJ, Kane D, Sunshine M, Narasimhan S, Nishizuka S, Reinhold WC, Zeeberg B, Ajay W, Weinstein JN: MatchMiner: a tool for batch navigation among gene and gene product identifiers. Genome Biol 2003, 4(4):R27. 10.1186/gb-2003-4-4-r27
Fundel K, Zimmer R: Gene and protein nomenclature in public databases. BMC Bioinformatics 2006, 7: 372. 10.1186/1471-2105-7-372
Tsai RT, Wu SH, Chou WC, Lin YC, He D, Hsiang J, Sung TY, Hsu WL: Various criteria in the evaluation of biomedical named entity recognition. BMC Bioinformatics 2006, 7: 92. 10.1186/1471-2105-7-92
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A 1998, 95(25):14863–14868. 10.1073/pnas.95.25.14863
Huang Z, Desper R, Schaffer AA, Yin Z, Li X, Yao K: Construction of tree models for pathogenesis of nasopharyngeal carcinoma. Genes Chromosomes Cancer 2004, 40(4):307–315. 10.1002/gcc.20036
Schaffer AA, Simon R, Desper R, Richter J, Sauter G: Tree models for dependent copy number changes in bladder cancer. Int J Oncol 2001, 18(2):349–354.
Hu ZF, Gao JH, Li W, Song YB, Li CL: [Differential gene expression profile of keloids: a study with cDNA microarray]. Nan Fang Yi Ke Da Xue Xue Bao 2006, 26(3):308–312.
Karakiulakis G, Papakonstantinou E, Aletras AJ, Tamm M, Roth M: Cell type-specific effect of hypoxia and platelet-derived growth factor-BB on extracellular matrix turnover and its consequences for lung remodeling. J Biol Chem 2007, 282(2):908–915. 10.1074/jbc.M602178200
Ide T, Kitajima Y, Miyoshi A, Ohtsuka T, Mitsuno M, Ohtaka K, Koga Y, Miyazaki K: Tumor-stromal cell interaction under hypoxia increases the invasiveness of pancreatic cancer cells through the hepatocyte growth factor/c-Met pathway. Int J Cancer 2006, 119(12):2750–2759. 10.1002/ijc.22178
Suderman M, Hallett M: Tools for visually exploring biological networks. Bioinformatics 2007, 23(20):2651–2659. 10.1093/bioinformatics/btm401
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 13(11):2498–2504. 10.1101/gr.1239303
Cerami EG, Bader GD, Gross BE, Sander C: cPath: open source software for collecting, storing, and querying biological pathways. BMC Bioinformatics 2006, 7: 497. 10.1186/1471-2105-7-497
Lee D, Kim S, Kim Y: BioCAD: an information fusion platform for bio-network inference and analysis. BMC Bioinformatics 2007, 8 Suppl 9: S2. 10.1186/1471-2105-8-S9-S2
This work was supported by grants from the National Natural Sciences Foundation of China (No. 30400496 and No. 90612020). We thank Dr. Feng Huang for technical support and Mr. Hongchun Gu for developing the website. We also thank Dr. Divakaran Jeevan for critical reading of the manuscript and helpful suggestions. We thank two anonymous referees for numerous excellent suggestions that led to substantial improvement in the manuscript.
Z–XH was involved in the development of the method and responsible for development of software as well as for writing the paper. H–YT was involved in the development of method and software. Z–FH was responsible for the generation of microarray data and was involved in the analysis of data. Y–BZ and JZ were involved in the analysis of data and the development of method. K–TY was responsible for the development of the method. All authors read and approved the manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.