I Introduction
Technological research has been revolutionized using high power processing techniques and has highly increased chunks of data available. Specifically, gene expression data have transformed research in biological field by its potential to supervise RNA concentration change in large number of genes concurrently available [1].
Clustering techniques are applied for analysis of gene expression matrix which creates subset of genes which display same expression behavior [2]. Conventional clustering algorithms are applied on complete dataset, testing each gene under all conditions. Clustering algorithms does not always give best results because most of the gene patterns exist under a subset of experimental conditions. Hence, clustering mechanism should be shifted to methods which can discover local patterns from gene expression data.
Biclustering is able to find local patterns by finding gene patterns under a subset of experimental conditions. Still, both clustering and biclustering are not enough when gene expression micro array data is analyzed and focus is drawn on effect of time on gene’s behavior. These type of longitudinal experiments allow thorough analysis of molecular processes where time plays a crucial role. Cell cycles, evolution of diseases and development at molecular level are some of the examples [3]. So, specific tools must be used for analysis of genes which are analyzed under specific conditions and time points. Hence, triclustering emerges as an effective tool for the same.
Triclustering algorithms find genes with similar expression in a subset of conditions along a time segment. Coherent tricluster is defined as set of genes that expresses same numerical value in certain time and condition or same behavior despite of exactly same numerical value i.e positively or negatively correlated changes in expression values of genes. Both positive and negative coherent clusters and regulation reaction among genes are useful for identifying effective phenotypes.
Triclustering using evolutionary computation (genetic algorithm) is presented in this paper which searches for group of genes that express similar patterns across both condition and time points i.e in three dimensional space. Most of the clustering and biclustering approaches define similarity of elements on the basis of distance measures
[4] but these measures are not accurate to find similarities in genes as correlations might still exist in some genes inspite of having different magnitude levels. Thus, two different evaluation measures are combined: one being mean square residue (MSR) and other being least square approximation (LSL). MSR is used for finding coherent triclusters, which is a 3D form of MSR derived from biclustering distance measuring technique used for discovering similar gene patterns. LSL is least square approximation (LSL) which finds triclusters with coherent behavior by calculating distance between slopes of least square lines of the tricluster discovered.In 2005, Zhao and Zaki [5] proposed algorithm named triCluster which extracts three dimensional pattern from gene expression data. Set of measures are also defined to calculate goodness of triclusters. In 2006, [6], a generalized form of triclustering (gtriCluster) was defined which focused on finding more coherent triclusters which are resistant to noise. In [7], the authors labeled the triclustering problem in association with its NPcompleteness. They proposed a filterlabeled stream with a parallel approach, hence, showing greater improvement in computational cost. Later in [8], coherent triclusters found regulatory relationships among genes. This algorithm was applied on both real and synthetic datasets. In [9], an algorithm was proposed to find time delayed clusters. It discovered cycle time of gene expression which is necessary for forming a network of gene regulation.
In [10], a new approach named LagMiner was introduced which finds timelagged three dimensional clusters. This algorithm was able to remove constraint of coherence, size of subspace, number of genes, regulation and length of time period. All the above algorithms were implemented on real and synthetic dataset. Triclustering is also implemented using evolutionary computation. Specifically, the multiobjective algorithm optimizes the most conflicting requirements of a tricluster. This algorithm was implemented on real data set [11].
In [12], a strategy to find triclusters in real valued data was introduced. Major concern in discovering tricluster was to find subspaces with large number of elements from gene expression matrix but this should not be the only concern in triclustering approach. So, an approach was proposed [13]
which aimed to find triclusters with low variance and could mine quantitative data. This algorithm was tested for synthetic data set and on crossspecies genome data set, thus being a major advancement in this topic.
A triclustering approach to find temporal dependent association rules in micro array data was introduced in [14]. The rules generated were able to present regulated relations among genes.
Remaining paper is arranged as : section ii@ discusses the basic concepts of triclustering, section iii@ explains the EA based algorithm implemented to derive triclusters. The experiment details, numerical simulation results are detailed in section iv@. Section v@ concludes the paper.
Ii Paradigm of Triclustering
Triclustering appear as an expansion of biclustering which focuses on evolution of genes at different time points under certain experimental conditions. So, from a gene expression matrix EM, which consist of genes , conditions , and time points , triclustering is defined as a procedure that finds triclusters from EM, where a tricluster T is defined as T =G x C x T, where , , and , ie, a subset of genes from gene expression matrix EM that carries information about the behavior of some genes under conditions C at times T.
Two quality measures are defined below which are crucial to get good quality triclusters. These definitions are defined in [15] .

Three Dimensional Mean Squared Residue: This definition is a modification of MSR used in bicluster quality measurement to three dimensional MSR which evaluates the similarity of triclusters which contain subsets of genes, conditions, and time points. This measure is denoted as .
(1) where is measured as
(2) where
: Mean of all gene values under conditions at a time t
: Mean of all gene values across time in certain condition c
: Mean of gene g in time under conditions
: Mean of all gene values in a particular condition and time
: Mean of values of gene at time under conditions
: Mean of a gene under a condition in all times
: Mean of all values in tricluster.A graphical depiction of terms used in eq. (1) can be understood from Fig. 1. It can be seen that calculates similarity of the tricluster on the basis of difference of each element of gene expression (Fig. 1(a)), mean of all conditions at all times for a gene g (Fig. 1(b)), mean of all genes at all times for a condition c (Fig. 1(c)), mean of all genes under all conditions at time t (Fig. 1(d)) with the mean of a condition c and a time t under all genes (Fig. 1(e)), mean of a gene g and a time t under all conditions (Fig. 1(f)), mean of a gene g and a condition c under all times (Fig. 1(g)), and the mean of all values in TC (Fig. 1(h)). If the value of is closer to zero, then the tricluster is more homogeneous. is also able to find negatively correlated genes.

Least Square Approximation: This term is defined by equation:
(3) It measures similarity of the least squares approximation of points in each graphic of the three views which depicts a tricluster. Firstly, for every time coordinate, conditions are represented on xaxis, expression values on yaxis and genes are outlined ( in (3)). Secondly, for every condition coordinate, times on xaxis, expression values on yaxis and genes are outlined ( in (3)). Thirdly, for every condition coordinate, genes on xaxis, expression values on yaxis and times are outlined ( in (3)). A representation of this can be seen in Fig. 2. All elements of numerator in Eq. (3) have values indicated in equation group (4, 5) in common.
(4) (5) where summarizes all genes of the individuals undergoing evaluation, gives squared summation of all genes, gives summation of all times and gives squared summation of all times. In first view, distance among all least square approximations are calculated from Eq. (6) ().
(6) (7) (8) (9) where TD in (7) represent measures, and are sum of each expression value of each time point, every combination of particular gene and conditions multiplied by genes in (8) and for every particular time, sum of all expression values at that time and every combination of particular gene and condition as mentioned in (9).
Similarly, is defined by Eq. (10) which gives distance between every least square approximation generated in second view.
(10) (11) (12) (13) where CD in (11) represents measures, and gives values for every particular condition, sum of every expression value of that condition, every combination of particular time and gene multiplied by genes in (12) and for every particular condition, sum of all expression values at that paticular condition and every combination of particular time and gene as mentioned in (13).
Eq. (14) gives which calculate distance between every least square approximation generated in third view.
(14) (15) (16) (17) where GD represents measures (14) and and are for every particular condition, sum of every expression value of that condition, every combination of particular gene and time multiplied by time (15) and for every particular condition, sum of all expression values at that particular condition and every combination of particular gene and time as mentioned in (16).
Iii Evolutionary Algorithm for TriClustering
TriEA is implemented on the basis of genetic algorithm. Following steps are involved in this evolutionary process: initial population is generated in initialization step which minimizes overlapping with already discovered triclusters, another is evaluation step where quality of each individual is calculated using a fitness function; a selection step decides which individuals should be allowed to survive to next generation, crossover exchanges genetic material between pairs of individuals to create offspring and mutation flips a particular individual bit to make sure genetic variability is maintained of future generations.
Iiia Chromosome Encoding
A potential tricluster solution TC is represented by each individual of population. It contains genetic information which is operated upon by the genetic operators. An individual chromosome consist of three sequences: one series of genes G, conditions C, and time points T. Fig. 3 illustrates the chromosome generation of triclusters.
These chromosomes are formed on the basis of expression matrix,
(18) 
where A is number of genes in the expression data, for all genes and . Similarly:
(19) 
where B is number of conditions in expression data, for all genes, and . Finally, T represents different time stamps:
(20) 
where C is number of gene samples measured over time, .
So, the chromosome bit string in the population is formed in following sequence
(21)  
Triclusters are represented in the form of binary strings of X+Y+Z length, X being genes, Y being conditions and Z being time points in expression matrix. If value of bit in individual is 1, it indicates respective gene, condition or time point have a place in that tricluster.
Let binary string be of 15 bits (5 for genes, 5 for conditions and 5 for time points). The genotypic representation of chromosome will be:
101101000111001
IiiB Initial Population Generation
The population is initialized randomly. A subset of genes, conditions and times are randomly generated, and tricluster TC = T x C x G are assigned. Remaining individuals are also randomly generated but keeping into consideration that individuals should be non overlapping with already generated individuals. In order to achieve this, every time a new individual is initialized, its indices of genes, conditions and times are stored. After that new individuals are initialized from random subset that did not appear in any previous individual.
IiiC Operation of Crossover
Parent and individuals are combined to create two offsprings and
. A probability of crossover (
) is associated with this operation. Genetic material of parents is combined by a random crosspoint in the chromosome and mixing the coordinates to form offsprings.Formally, let and be two parents at some iteration i. Resultant offsprings are
(22) 
where, f randomly selects subset of (C, G, T) from parents with crossover probability ().
The procedure can be formalized as follows. Let = and = be two sequences of chromosome, where and m, n length of chromosome.
Let P be number randomly generated between 1 and min of (m, n). Two new offsprings are formed as follows:
(23) 
(24) 
IiiD Operation of Mutation
Depending on the probability of mutation (), any chromosome can be mutated. Standard mutation is applied which flips the value of a single bit randomly in chromosome.
IiiE Selection Operator
In the reproduction stage, offsprings are created by making them compete with other individuals for the placement in next generation. Elitism is applied which means only the best chromosomes should be allowed to survive to next generation. Tournament selection operator is applied with size of 2.
IiiF Proposed Fitness Equation
Fitness of every individual helps algorithm to find best solutions which are carried forward in further generations. Two quality measures are combined in TriEA algorithm, one being 3D Mean Square Residue measure (MSR), referred as . The second one is .
Weights term
(25) 
where , and represents weights for genes, conditions and times in the tricluster. Higher value of weights indicate that triclusters with many components are found by triEA.
Distinction term
(26) 
where (Coordinate Distinction Number of g) : Genes which are not present in tricluster solution being evaluated.
(Coordinate Distinction Number of c) : Conditions which are not present in tricluster solution being evaluated.
(Coordinate Distinction Number of t) : Time points which are absent in tricluster solution being evaluated.
, , are weights assigned at prior to genes, conditions and time points. Distinction calculates difference between chromosome under evaluation and triclusters already discovered.
If the value of weights is increased, nonoverlapping tricluster solutions can be discovered.
Fitness Equation
Following equation defines this function:
(27) 
This is a minimization equation which consist of four factors , Weight facto Distinction factor, particular MSR, LSL factors.
Iv Numerical Simulation Results and Discussion
The algorithm TriEA is implemented using R i386 3.3.2 and executed on system of 32 bit operating system and Intel(R)Core(TM)i74790 CPU @ 360GHz processor with 4 GB of RAM. TriEA is applied to the yeast (Saccharomyces cerevisiae) cell data [16]. The objective of this project is to find genes whose mRNA levels are controlled by cell cycle processes. When TriEA is applied to dataset of such form, meaningful sequences in cell cycle can be discovered. In this experiment, 6179 genes are examined under 6 conditions named clb2, cln3, cdc15, cdc28, pheromone and elutriation. Gene profiles are captured at 2 time points for cln3, 18 for pheromone, 2 for clb2, 24 for cdc15, 17 for cdc28 and 14 for elutriation. For applying TriEA, conditions with only 2 time points are not considered and first 14 time points of pheromone, elutriation, cdc15, cdc28 are considered in dataset. Hence, the expression matrix consist of 6179 genes, 4 conditions and 14 time points. 200 genes are taken to carry out this experiment due to its increased search space and high computation cost. The original expression data is preprocessed using maxmin normalization method such that all the values of matrix lies between 0 to l.
(28) 
Parameter values for triEA are given in TABLE i@. Values that are missing in dataset are replaced with randomly generated values between 0 and 1.
Parameter Name  TriEA 

Population Size  20 
Number of generation  100 
Probability of Crossover  0.95 
Probability of Mutation  0.50 
Weight for conditions  0.1 
Weight for genes  0.1 
Weight for time  0.1 
Genes from original dataset taken  200 
Threshold  1050 
Number of triclusters  20 
Tric No  Fitness Value  LSL  Weight  Distinction  MSR 

1  6246.74  19.74  1.0  0.0505  6228.04 
2  139141.44  492.01  7.7  0.0043  138657.13 
3  429.41  12.59  0.8  0.0280  417.64 
4  4003.45  662.39  1.7  0.0106  3342.77 
5  4120.64  576.88  1.6  0.0113  3545.37 
6  432.86  55.89  0.8  0.0280  377.79 
7  2837.64  123.84  1.0  0.0218  2714.82 
8  10885.81  819.77  2.5  0.0077  10068.54 
9  10763.06  834.15  1.6  0.0109  9930.51 
10  4533.08  745.21  1.6  0.0113  3789.48 
11  11241.41  171.32  1.4  0.0129  11071.50 
12  7714.68  352.65  2.2  0.0085  7364.23 
13  14278.17  693.94  1.8  0.0095  13586.042 
14  16013.35  134.89  1.2  0.0209  15879.68 
15  25446.38  125.77  3.4  0.0063  25324.015 
16  55523.39  737.98  4.9  0.0052  54790.31 
17  6966.46  49.03  2.1  0.0088  6919.50 
18  14581.88  225.15  2.5  0.0077  14359.23 
19  3589.76  255.23  1.5  0.0120  3336.03 
20  23074.08  873.73  1.9  0.0089  22202.25 
After the convergence of TriEA, 20 resultant triclusters were derived with minimum MSR and LSL score and it was observed all of them have least square approximation (LSL) score less than the threshold i.e 1050. The detailed results are shown in TABLE ii@. The average value of LSL for all triclusters obtained is 398.12. Tricluster obtained with least LSL value is 12.59. Also, average value of MSR for all triclusters obtained is 17545.24 which is quite large because this algorithm focus on triclusters with minimum LSL value but not least MSR value of triclusters. Least MSR value obtained is 417.64. Due to large MSR values of triclusters, overall fitness value also increases.
From Fig. 4 convergence of GA for 100 generations shows that triclusters with minimum MSR and LSL are preferred in evolution process resulting in to minimum fitness value of the tricluster. Thus, in order to find triclusters with minimum overall fitness value, triclusters with minimum MSR and LSL should be preferred in evolution process.
The results in TABLE ii@ are of TriEA for 200 genes, 4 conditions and 14 time points.
Fig. 4 shows convergence of GA towards minimum value for 100 generations.
V Conclusion and Future Work
For Triclustering, TriEA is implemented for 200 genes using two quality measures i.e MSR and LSL. It can be concluded from results that GA converges better when both the quality measures are combined to form one fitness function. TriEA yields good quality triclusters with least overlapping among triclusters discovered. But the tricluster quality can be improved more, if both MSR and LSL are minimized by algorithm simultaneously.
The microarray data is mined for longitudinal experiments using algorithm TriEA, but this algorithm can also be applied in other biological domains, for example gene expression data can be combined with gene regulatory data by methods of replacing dimensions of time by chipchip information depicting transcription factor gene interactions that provides information about gene regulatory networks. This algorithm can be applied on mining RNAseq data repositories.
References
 [1] P. O. Brown and D. Botstein, “Exploring the new world of the genome with dna microarrays,” Nature genetics, vol. 21, pp. 33–37, 1999.

[2]
C. RubioEscudero, F. MartínezÁlvarez, R. RomeroZaliz, and I. Zwir, “Classification of gene expression profiles: comparison of kmeans and expectation maximization algorithms,” in
Hybrid Intelligent Systems, 2008. HIS’08. Eighth International Conference on. IEEE, 2008, pp. 831–836.  [3] Z. BarJoseph, “Analyzing time series gene expression data,” Bioinformatics, vol. 20, no. 16, pp. 2493–2503, 2004.
 [4] R. Tibshirani, T. Hastie, M. Eisen, D. Ross, D. Botstein, P. Brown et al., “Clustering methods for the analysis of dna microarray data,” Dept. Statist., Stanford Univ., Stanford, CA, Tech. Rep, 1999.
 [5] L. Zhao and M. J. Zaki, “Tricluster: an effective algorithm for mining coherent clusters in 3d microarray data,” in Proceedings of the 2005 ACM SIGMOD international conference on Management of data. ACM, 2005, pp. 694–705.
 [6] H. Jiang, S. Zhou, J. Guan, and Y. Zheng, “gtricluster: a more general and effective 3d clustering algorithm for genesampletime microarray data,” in International Workshop on Data Mining for Biomedical Applications. Springer, 2006, pp. 48–59.
 [7] R. Braga Araújo, G. H. Trielli Ferreira, G. H. Orair, W. Meira, R. A. Celso Ferreira, D. Olavo Guedes Neto, and M. J. Zaki, “The partricluster algorithm for gene expression analysis,” International Journal of Parallel Programming, vol. 36, no. 2, pp. 226–249, 2008.
 [8] Y. Yin, Y. Zhao, B. Zhang, and G. Wang, “Mining timeshifting coregulation patterns from gene expression data,” in Advances in data and web management. Springer, 2007, pp. 62–73.
 [9] G. Wang, L. Yin, Y. Zhao, and K. Mao, “Efficiently mining timedelayed gene expression patterns,” IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics), vol. 40, no. 2, pp. 400–411, 2010.
 [10] X. Xu, Y. Lu, K.L. Tan, and A. K. Tung, “Finding timelagged 3d clusters,” in Data Engineering, 2009. ICDE’09. IEEE 25th International Conference on. IEEE, 2009, pp. 445–456.
 [11] J. Liu, Z. Li, X. Hu, and Y. Chen, “Multiobjective evolutionary algorithm for mining 3d clusters in genesampletime microarray data,” in Granular Computing, 2008. GrC 2008. IEEE International Conference on. IEEE, 2008, pp. 442–447.
 [12] K. Sim, Z. Aung, and V. Gopalkrishnan, “Discovering correlated subspace clusters in 3d continuousvalued data,” in Data Mining (ICDM), 2010 IEEE 10th International Conference on. IEEE, 2010, pp. 471–480.
 [13] Z. Hu and R. Bhatnagar, “Algorithm for discovering lowvariance 3clusters from realvalued datasets,” in Data Mining (ICDM), 2010 IEEE 10th International Conference on. IEEE, 2010, pp. 236–245.
 [14] Y.C. Liu, C.H. Lee, W.C. Chen, J. Shin, H.H. Hsu, and V. S. Tseng, “A novel method for mining temporally dependent association rules in threedimensional microarray datasets,” in Computer Symposium (ICS), 2010 International. IEEE, 2010, pp. 759–764.
 [15] D. GutiérrezAvilés, C. RubioEscudero, F. MartínezÁlvarez, and J. C. Riquelme, “Trigen: A genetic algorithm to mine triclusters in temporal gene expression data,” Neurocomputing, vol. 132, pp. 42–53, 2014.
 [16] P. T. Spellman, G. Sherlock, M. Q. Zhang, V. R. Iyer, K. Anders, M. B. Eisen, P. O. Brown, D. Botstein, and B. Futcher, “Comprehensive identification of cell cycle–regulated genes of the yeast saccharomyces cerevisiae by microarray hybridization,” Molecular biology of the cell, vol. 9, no. 12, pp. 3273–3297, 1998.
Comments
There are no comments yet.