cancers in many tissues are heterogeneous , and the efficacy of therapeutic interventions depends on the specific subtype of the malignancy . hence , early and accurate identification of the cancer subtype is critical in designing an effective personalized therapy .
current methods for assessment rely on microscopic examinations of the malignant tissue for previously established histopathological abnormalities .
unfortunately , such features may not be apparent during early stages of the disease and moreover , differentiating between abnormalities in distinct cancer subtypes can be challenging .
recent advances in high - throughput genomics offer an exciting new alternative for early and reliable cancer prognosis .
mutations that underlie a malignancy modify the levels of many genes within a cell ; the goal of gene expression profiling is to define a signature for each cancer subtype through statistically significant up-/downregulation of a panel of genes .
the national institutes for health , through the cancer genome atlas ( tcga ) [ 1 , 2 ] , will aid this effort by establishing large sets of genomic data on human cancers in at least 20 tissues [ 38 ] .
the premise behind tcga is that statistically significant changes in gene expression levels due to malignant mutations can be placed in a few groups associated with subtypes , and that unsupervised ( or semisupervised ) clustering algorithms can be used to uncover these partitions .
this is illustrated through a schematic  malignancy  that can be partitioned using the expression levels of two genes . in this
schematic , each patient sample is represented by a point on a plane , see figure 1 .
the basic observation is that , while the patient samples are distributed over a broad range , there are  pockets  of high concentration , which can be identified using traditional clustering methods . the members of the pocket define the  core  samples of a cancer subtype . however , the presence of significant levels of noise in genomic data makes the partitioning a nontrivial task .
the variability is due to both the subject dependence of the expression levels and to imperfections in microarray technology .
unfortunately , costs associated with microarray experiments prohibit the use of a large set of replicates to reduce the effective error rates .
one focus of tgca has been glioblastoma multiforme ( gbm ) , the most common and aggressive form of brain cancer [ 9 , 10 ] .
tcga provides the expression levels of 11861 genes in 200 gbm and 2 normal brain samples [ 1 , 2 ] .
reference   identifies 1740 genes with consistent expression across affymetrix huex , affymetrix u133a ( affymetrix , santa clara , ca , usa ) , and agilent 244k common genomic hybridization arrays ( agilent technologies , santa clara , usa ) , to be used for the subgrouping .
they search for common partitions under sampling of genes and patients . the resulting consensus clustering   yields four robust clusters whose class boundaries are statistically significant .
173  core  representatives of the four groups were identified , and an 840-gene signature was determined on the basis of lowest cross - validation and prediction error .
this genomic partitioning of the 173 core samples was found to be consistent with the grouping into the four known subtypes classical , mesenchymal , proneural , and neural of glioblastoma . in this work ,
we introduce a new algorithm for gene expression profiling , which is illustrated through an application to gbm .
this approach avoids several difficulties associated with clustering algorithms commonly adopted to partition large sets of genomic data .
it provides robust partitions of patients and identifies a compact set of genes used to distinguish the clusters .
in particular , for gbm , the clusters of patients are sharply defined leading us to establish the quintessential genetic profile for the subtypes .
further , the newly identified set of genes will allow tumors of new patients to be subtyped quickly after diagnosis , which perhaps could lead to a more personalized and successful treatment regimen .
most of these genes have not been previously implicated in brain cancer or cancer in general ; they may be unidentified members of the cancer network .
unfortunately , significant structural variation in tumor samples renders it impossible to determine subtype by histological methods .
thus , we are unable to have independent corroboration of our results and are left to compare our calculation with the clustering analysis of the tgca data .
our algorithm requires the panel of genes used for clustering to be predetermined . due to stochasticity
the choice of an excessively large set of genes is likely to be counter productive as well , at least in part due to contamination from the many genes that are irrelevant for the comparison .
a large gene set will cause a another well - known difficulty , namely , that of partitioning a ( relatively ) small number of objects in a high dimensional space .
thus , it is important to preselect a sufficient yet compact set of genes for clustering .
this issue is resolved using a statistical invariant and persistent homology [ 1315 ] , one of the most exciting recent developments in topology .
this yields the appropriate number of clusters and provides information on the genetic proximity of the subtypes .
additionally , refinements to the panel of genes used in the partitioning are possible through repeated application of the method .
there is no natural  genomic distance  between two samples , although the correlation distance ( defined below ) is used in most studies .
hence , it is desirable to use a topological approach ( such as persistent homology that is adapted here ) for partitioning , rather than one which relies on computation of metrics such as eigenvalues .
this is likely the reason why hierarchical clustering   typically provides a more robust partitioning than spectral clustering [ 1720 ] or those based on principal components analysis .
empirical cumulative distribution functions , which quantify the proximity of a given clustering to a 
perfect  grouping ( where the membership of each object is unambiguous ) is one of the more successful approaches to address the problem in genomic data .
next , we wish to introduce a prediction model to assign the cancer subtype . to this end
, we randomly select 20 (  test  samples ) of the 202 patient samples which are only to be used for validation .
( the algorithm is  robust  in the sense that the patient subgroups obtained from the algorithm do not depend on the number of genes selected for the panel . )
one partition contains 60 samples , 27 of which labeled classical by tcga and 28 of which labeled mesenchymal .
its partitions contain 18 ( 17 classical ) and 14 ( all mesenchymal ) patient samples .
thus , our method identifies the three primary clusters of gbm . in our calculation , the neural group was not found to be a single , coherent cluster .
interestingly , core nodes of the biological network underlying gbm such as tp53 , pten , rb1 , nf1 , pik3r1 , and pik3ca [ 1 , 22 ] are not in this list .
finally we use the 59-gene set to assign the cancer subtypes in the 20 test samples , comparing the conclusions with those of the clustering calculation  http://tcga-data.nci.nih.gov/docs/publications/gbm_exp/  .
the assessment of the subtype of each test sample is made using the mean of its correlation distance from the corresponding core patient samples .
five of these were not assigned a group by either our calculation or that of tcga .
four samples , which were labeled neural , were not provided an assignment by our algorithm .
a gene can only be useful in differentiating between subtypes of cancer if its expression levels in the subtypes lie on statistically distinct distributions [ 1 , 23 ] .
a necessary ( but not sufficient ) condition is the bimodality of the corresponding distribution for all patient samples .
( note that , at this stage , we do not know the mean values , variances , or the membership of the putative subgroups .
hence standard statistical tools such as discriminant analysis [ 24 , 25 ] can not be used to decide if and how patient samples are to be partitioned . )
we propose the use of the nondimensionalized standard deviation(1)=var  ( x)[|x[x]| ] , 
 to select such genes . here
[x ] represents the expectation value of a random variable x and var  ( x ) = [x ]  [x ] its variance .
the definition is inspired by kurtosis , which is a scale - invariant measure of the  width  of a distribution .
unfortunately , it is difficult to estimate the fourth moment reliably from 182 noisy points ( gene expression levels ) .
 is larger for distributions with broad tails : for gaussian and biexponential distributions its values are /2 and 2 , respectively .
its value for a finite - range characteristic function ( which has no tails ) is 2/3 .
variation of  as a function of the distance between two gaussian distributions is shown in figure 2(a ) , demonstrating that bimodality of a distribution can be inferred from the value of .   the gene expression profiles for the 202 patient samples are given at  http://tcga - data.nci.nih.gov / docs / publications / gbm_exp/.  the table 
unifiedscaled.txt  contains mean expression levels ( from three replicates ) for 11861 genes of each patient sample .
figure 2(b ) shows the values of  for the 11861 genes in increasing order .
( for each gene , we have discarded the extreme outliers , which are more than 8 standard deviations away from the mean . ) a very small fraction of genes have abnormally small values of  , which we associate with multimodal distributions .
the eight smallest values of  are found for genes eif1ay , rps4y1 , ddx3y , uty , usp9y , jarid1d , zfy , and nlgn4y .
the homology analysis ( see below ) shows that these eight genes and cyorf15b belong to one group , that is , the correlation distance between these genes is small .
however , each gene in this group belongs in the y - chromosome , and hence the associated partitioning is that between male and female subjects .
, we select n = 60 that have the smallest values of . the final partitioning of patients is independent of n for choices between 50 and 120 , and persistent homology fails to find appropriate partitions when n is outside this range .
not all of the 60 genes selected using  may be useful for subdividing patient samples .
our second refinement is to search for groups of genes whose expression levels change coherently between patient samples .
this coherence can , for example , result from the genes ' membership in a subnetwork .
for this task , we represent a gene gn by its expression levels xn   in all m = 182 patient samples ; the 60 genes form a  point cloud .  as is traditional in clustering algorithms , the proximity d(xn , xm ) between two points xn and xm is chosen to be the correlation distance : 
 ( 2)d(xn , xm)1c(xn , xm)=1[xnxm][x][xm]var  [ xn]var  [ xm ] , 
 where c(xn , xm ) is the pearson correlation function . we have now converted the search for critical genes to clustering a point cloud .
this task could be implemented through algorithms such as hierarchical clustering [ 16 , 29 ] , spectral clustering [ 20 , 30 ] , or community structures .
spectral approaches are highly sensitive to the precise genes chosen for the analysis and to cutoff values , and hierarchical clustering does not provide an unambiguous estimate for the number of partitions .
consider a point cloud and a set of connections between a subset of pairs of points .
a k - simplex is a set of ( k + 1 ) points , each pair of whose vertices is connected .
for example , a 1 simplex is a line along with its two vertices , a 2 simplex is a triangle along with its edges and vertices , and a 3 simplex is a tetrahedran along with its faces , edges , and vertices .
the face of a k simplex is a subset which is a simplex of lower order .
for example , the vertices , edges , and triangular sides of a tetrahedran are its faces . a k  simplicial
complex is an object formed by gluing together faces of a set of simplices each of whose dimension is less than or equal to k ; the intersection of any two such simplices is required to be a face of both .
for example , an object formed by gluing a set of triangles on their edges or vertices , with perhaps some lines and isolated points , is a 2-simplicial complex .
the question is addressed by constructing a simplicial complex and searching for its disjoint components .
first , for a prespecified cutoff radius rc , connect points xm and xn if and only if d(xn , xm ) < rc . the object so formed is a simplicial complex . for sufficiently small rc
when rc is sufficiently large , all pairs of points are connected , and the simplicial complex is a single component .
simplicial complexes can be characterized using topological invariants ( i.e. , features that do not change under smooth distortions of space ) .
one of these invariants is the number of disjoint components , which is referred to as the zeroth betti number ( betti0 ) .
the first betti number , betti1 , is the number of connected components with a two - dimensional hole .
the second betti number ( betti2 ) is the number of surfaces enclosing a three - dimensional region .
the first three betti numbers for the crust of a bagel are 1 , 2 , and 1 .
betti numbers can provide a more comprehensive invariant characterization of the point cloud than possible through clustering .
( persistent homology has been used successfully to search for recurrent genomic instability in breast cancer [ 31 , 32 ] . )
a partition will be robust if the components ( of the point cloud ) are far apart from each other ; consequently , the number of components will remain unchanged for a broad range in rc .
a generalization of this statement forms the novel proposal in [ 13 , 14 ] , which associates robust characteristics of the point cloud with  persistent  topological invariants .
our computations are implemented using the package jplex ( available at  http://comptop.stanford.edu/programs/jplex/index.html  ) .
we note first that if r1 < r2 , the simpicial complex at cutoff r1 is contained in that at r2 .
the primary construction used in jplex is the  rips stream ,  which computes the simplices at cutoff radii rc in a prespecified range , and assigns a  filtration time  for each simplex , when it first appears .
 barcodes ,  such as those shown in figure 3 , show how the number of components evolves as rc is increased .
the abscissa is the cutoff radius and each horizontal line represents a component that maintains its identity as rc increases .
it should be noted that the membership of a component may increase with rc . in figure 3(a ) , we only show components that contain a minimum of nmin points .
this constraint is imposed by constructing a new rips stream which contains the original simplices ; the filtration time associated with each simplex in the new stream is defined to be the maximum of the original filtration time and the time when the component reaches nmin elements .
the top panel of figure 3 shows the results of the analysis for the 60 genes with nmin = 5 .
the membership of the groups is maximized at rc  0.56 , just prior to their combining into a single component .
( the persistence width of the component that begins after rc = 0.56 decreases as the number n of genes used for the analysis increases . ) genes in each component are found using the jplex routine  verticesineachcomponent .  at rc = 0.56 , the components contain 17 and 20 genes , respectively .
however , we find that identical partitioning of the patient samples can be achieved with a smaller number of genes . for the analysis here , we select the first ng = 15 genes to join each component .
this choice is made in order to include a sufficiently inclusive , yet not too large a set of genes for the analysis .
the two gene sets are 
 ( 3)c1={ugt8,mbp , c11orf9,mog , klk6,rp1135n6.1,rab33a , dcx , gpr17,cxorf1,lrrtm4,tmsl8,sn , ap91,atp10b , luzp2},c2={meox2,pipox , fgfr3,dlc1,pla2g5 ,  egfr , grik1,erbb2,ptrf , postn , baiap3 ,  pdlim4,kcnf1,eya4,spag4}. 
 genes in c1 lie on chromosomes 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 10 , 13 , 17 , 20 , and 21 .
genes meox2 , egfr , dcl1 , efgr , ptrf , and eya4 have the gene ontology ( go ) classification  nucleus .
 genes meox2 , dcl1 , erbb2 , ptrf , eya4 , spag4 have the classification 
cytoplasm ,  and fgfr3 , grik1 , kcnf1 have the classifications  membrane  and  integral to membrane .  genes in c2 belong to chromosomes 1 , 2 , 4 , 5 , 6 , 9 , 11 , 15 , and the x - chromosome . c11orf9 and klk6 have the gene ontology classification  nucleus ,  while dcx , klk6 , and tmsl8 have the classification 
cytoplasm ,  and ugt8 , c11orf9 , mog , gpr17 , lrrtm4 , and atp10b have the classification  membrane .
 the 182 patient samples are partitioned using the expression levels of the 30 genes in c1 and c2. each sample is represented by the expression levels y of these 30 genes . as before , the distance between two points ( patient samples ) ym and yn in  is defined as the correlation distance .
the barcode for the analysis , where we retain only groups of size 10 or larger , is shown in figure 3(b ) .
two groups are seen to persist over a significant range in the cutoff radius rc .
since the components grow with rc , we choose rc  0.34 , in order to maximize their membership .
they contains 84 and 76 patient samples , and the heat map [ 3335 ] is shown in figure 4(a ) .
the patient samples in each group are arranged in the order they appear in the component as rc is increased .
the fold inductions for genes in the left - most samples in each group are significantly closer than those joining the component later .
we propose to only include the most coherent patient samples in the  core  group associated with each category , see figure 1 . with this trimming
the new heat map , figure 4(b ) , provides a significantly sharper partitioning than those computed through many previous approaches [ 1 , 2 , 11 , 12 , 36 , 37 ] , although the membership in each group is smaller .
the first group of 60 contains 27 patient samples that have been assigned the classical subtype and 28 samples that have been assigned the mesenchymal subtype of gbm by the clustering analysis of tcga .
37 of the 44 patient samples in the second group have been assigned to the proneural subtype by tcga .
the rest of the samples in this group are unassigned with the exception of one sample labeled neural and one labeled mesenchymal .
one of our goals is to search for a panel of genes to differentiate between the subtypes of gbm . thus far
, we have used the genes in c1 and c2 to find two subgroups of the 182 patient sample , and ( using the heat map ) select the most tightly correlated members in each subgroup . however , when only the 60 and 44 core samples are considered , an alternative gene set may provide a sharper separation of the patient samples .
we thus repeat the computations of the last three subsections , using only the 104 patient samples in the two core groups .
the new gene groups are 
 ( 4)c1={ugt8,rp1135n6.1,dcx , rab33a , gpr17 ,  erbb3,sox10,pak7,il1rapl1,dnm3 ,  snap91,atp1a3,rundc3a , dusp26,ephb1},c2={elovl2,arsj , flj21963,pipox , nr2e1 ,  znf217,pla2g5,meox2,ephb4,dcl1 ,  postn , lama2,egfr , kcnf1,eya4}. 
 genes in c1 and c2 are spread among several chromosomes .
genes rp11 - 35n6.1 , ugt8 , gpr17 , il1rapl1 , snap91 , and ephb1 of c1 have the go classifications 
membrane  and/or  integral to membrane ,  while dcx , rab33a , and dnm3 have the classification  intracellular .
 genes nr2e1 , znf217 , meox2 , egfr , and eya4 of c2 have the go classification  nucleus  and arsj , pipox , pla2g5 , postn , lama2 , egfr , and kcnf1 have the classification  extracellular region .  the heat map , shown in figure 4(b ) , clearly indicates that genes in c1 and c2 provide a significantly better separation of the core members in the two groups .
we conducted the corresponding analysis for the validation data set combined from four public cohorts ( see  http://tcga-data.nci.nih.gov/docs/publications/gbm_exp/  ) using genes in c1 and c2 .
we again find a clear difference in fold changes of the genes between the two partitions .
in addition , the fractional sizes of the partitions are similar to those of figure 4(b ) .
( the assignments of the cancer subtype of validation samples are not given at the tcga site . )
the next step is to partition patient samples in group a , which we found to contain exclusively classical and mesenchymal samples .
the computations described above are repeated , however , only using 60 patient samples assigned to the partition .
the gene groups for the partition are 
 ( 5)c1(cm)={ascl1,mpped2,cspg5,bcan , dscam ,  dpf1,znf606,znp30,znf419,sgta ,  znf8,dgkb , egfr , phc1,blm},c2(cm)={tnfaip8,slc25a24,olfml2b , sft2d2 ,  kynu , dkk1,rras , gusb , heph , enpp4 ,  ndn , spag4,mtap , klhl9,col11a1}. 
 genes in c1 lie in chromosomes 1 , 3 , 7 , 11 , 12 , 15 , 19 , and 21 .
dpf1 , sgta , and the four zinc fingers belong to chromosome 19 ( 19q13 ) and have the go classification  intracellular .
 they along with ascl1 , blm , and phc1 have the go classification  nucleus .
 genes in c2 belong to chromosomes 1 , 2 , 5 , 6 , 7 , 10 , 15 , 19 , and the x - chromosome .
tnfaip8 , slc25a24 , kynu , spag4 , and mtap have the classification  cytoplasm .
 the first group of patient samples contains 18 members , 17 of which have been assigned to the classical subtype in .
the second group contains 14 samples , all of which have been assigned the mesenchymal subtype in .
the heat map for the bisection of the classical / mesenchymal group is given in figure 5 .
the complete list of core patient samples is given in supplementary materials ( see supplementary materials available online at doi:10.5402/2012/381023 ) .
a similar bisection of patient group b , which contains the proneural samples , was not possible .
fold inductions of genes belonging to c1 , c2 , c1 , and c2 can be used to predict the membership of a new sample in the three core groups .
the choice is based on the mean correlation distance between the patient sample and members of each core group .
the first step is to determine if the new sample belongs to the classical / mesenchymal or proneural core groups .
the determination is made using genes belonging to c1 and c2 and computing the mean correlation distance between the new sample and those in each partition .
however , we must first determine the  radius  ( for the correlation distance ) within which the membership for each group is assigned .
for every patient sample in the classical / mesenchymal core group , we compute its mean correlation distance from the remaining members of the group .
the largest of these values is defined to be the radius of the classical / mesenchymal subgroup of patient samples .
this choice guarantees that all core patient samples in the group are assigned correctly , and that those outside the core group are not assigned to it .
the computed cutoffs for the classical / mesenchymal and proneural core groups are rcm = 0.70 and rp = 0.42 , respectively .
a new patient sample whose mean correlation distance from members of the classical / mesenchymal core group is smaller than rcm is assigned to that group , and any sample closer than rp to the proneural subgroup is assigned to it .
next , we use genes in c1 and c2 to assign the membership of any patient sample in the classical / mesenchymal core group . the corresponding cutoff radii are rc = 0.63 and rm = 0.30 .
table 1 shows the mean correlation distances between the test samples and the core groups .
the last column of the table gives the assignment of the gbm subtype given by tcga ( see  http://tcga-data.nci.nih.gov/docs/publications/gbm_exp/  ) .
five of the test samples did not belong in the core subgroups identified in , and our calculation agrees . of the remaining 15 patient samples , 5 do not belong to any of our three core groups .
( in fact , 4 of these have been assigned to the neural subgroup . )
robust prognostics for many subtypes of cancer are yet to be discovered [ 23 , 32 ] .
high - throughput genomics offers a possible approach for early and reliable cancer prognosis [ 8 , 38 ] .
clustering algorithms require the definition of a genomic distance between patient samples . since there is no natural measure , topological approaches to partitioning are likely to prove more robust .
second , although techniques have been proposed to identify the optimal number of partitions of a set of objects , its value can be sensitive to the choice of the gene panel used for profiling .
it was used to partition the 202 gbm patient samples whose genetic profiles are given at the publicly accessible tcga site .
we identified a set of genes most useful for the partitioning and to determine core groups of samples for each subtype of gbm .
we argued that a necessary condition for the inclusion of a gene was that its fold inductions lie on a distribution with a ( relatively ) small value of the nondimensionalized standard deviation  , indicating a bimodal distribution .
we then used persistent homology to streamline the set further by only retaining genes that belong in groups ; genes within a group exhibit highly correlated variations of expression levels between patient samples .
many of them do not play a pivotal role in gene regulatory networks associated with gbm [ 14 , 22 , 23 , 39 ] and may be downstream nodes .
it would be interesting to determine if they can be implicated directly in the histopathological criteria used to define the malignancies .
next , we use the differential expression of genes in c1 and c2 to bisect the patient samples .
the two partitions contain 84 and 76 patient samples . using the heat map for guidance
, we selected 60 and 44 samples from the two partitions as its core members .
once the core patient samples in each group are known , it is possible to determine the most suitable gene set to differentiate between the partitions . the new gene groups c1 and c2 provide a significantly sharper heat map .
we find that the first core group contains patient samples of the classical and mesenchymal subtypes , while the second contains proneural samples .
( 1 ) patient subgroups derived from the analysis depend only weakly on the number of genes used for the persistent homology step ( n = 60 for the results reported here ) . for the gbm data , barcodes ( figure 3(a ) ) for a range of choices for n
between 50 and 120 are similar , with the component appearing at rc  0.6 lasting for a smaller range in cutoff as n increases .
values of n < 50 do not give uniform results , and those larger than 120 can not be justified from figure 2(b ) . ( 2 ) the patient groups are very weakly dependent on the number of genes ng chosen for the panels c1 and c2 ( ng = 15 for each set ) .
we could choose ng as small as 10 or use all genes in the sets c1 and c2 for the analysis with similar results .
a unique aspect of our approach is that each subgroup can now be further bisected , albeit using different gene sets .
specifically , the classical / mesenchymal group was subdivided using sets c1 and c2 , each containing 15 genes .
it is known that that classical patient samples exhibit a high level of expression of egfr [ 1 , 9 ] .
consistent with this , we find that egfr is over expressed in the classical / mesenchymal partition in figure 4 and in the classical group in figure 5 .
the proneural subtype of gbm is associated with high alterations in tp53 , pdgfra , and idh1 .
it is possible that , although there are large alterations in the levels of these genes , the mean values in the ensemble ( of proneural samples ) do not change significantly .
the high levels of mutations of nf1 of mesenchymal samples are not reflected in the gene sets as well .
our analysis suggests that patient samples in this group lie , not on a single partition , but on several small partitions .
we used the panel of 59 genes in c1 , c2 , c1 , and c2 to introduce a prediction model for the subtypes of gbm .
the test involved the computation of the mean correlation distance between a new sample and members of each core group .
fifteen out of the 20 test sample predictions agreed with the results from the clustering calculation .
four of the remaining samples were categorized as neural by the clustering method but were unassigned by ours since the neural group was not found to be a single cluster . of the genes in c1 , c2 , c1 , c2 , only egfr has been previously implicated in gbm . however , other gene groups may have functions related to cancer ( the functional annotations can be found in  genecards ,  the human gene compendium , weizmann institute of science ) .
one example is a group related to intercellular scaffolding : dcx translates to a microtubule - associated protein , dnm3 is involved in producing microtubule bundles , postn plays a role in extracellular matrix mineralization , rras regulates the actin cytoskeleton , and coll11a1 may play a role in fibrillogenesis .
another group is known to be expressed in brain function : rab33a is a member of the ras oncogene superfamily , gpr17 may mediate in brain damage , and nr2e1 may be required for brain development .
a third group of genes is associated with differentiation , development , and the cell cycle : they include edhb4 , ascl1 involved in early stages of development , mpped2 , which plays a role in the development of the nervous system , cspg4 , which functions as a growth and differentiation factor , bcan involved in terminal differentiation , meox2 , which plays a regulatory role in muscle cells , and ndn , which is a growth suppressor involved in cell cycle arrest .
the zinc fingers zfp30 , znf419 , znf8 , and znf606 are involved in transcriptional regulation .
we have introduced a novel algorithm for genomic subtyping of cancer samples . unlike prior clustering techniques , we predetermine a compact set of genes to be used in the analysis through the nondimensionalized standard deviation and persistent homology .
importantly , the gene set selected for the analysis depends on the group of patient samples to be analyzed .
hence , if needed , partitions obtained from the first step can be further bisected .
our approach not only provides a robust partitioning of patient samples , but can also be used to identify a panel of genes to be used as a biomarker .
the gene panels have been used to introduce a prediction model to determine if a new patient sample belongs to the core subgroups associated with classical , mesenchymal , and proneural subtypes of gbm .
 early and accurate diagnoses of cancer can significantly improve the design of personalized therapy and enhance the success of therapeutic interventions . 
 histopathological approaches , which rely on microscopic examinations of malignant tissue , are not conducive to timely diagnoses . 
 high throughput genomics offers a possible new classification of cancer subtypes . 
 unfortunately , most clustering algorithms have not been proven sufficiently robust . 
 we propose a novel approach that relies on the use of statistical invariants and persistent homology , one of the most exciting recent developments in topology . 
 it identifies a sufficient but compact set of genes for the analysis as well as a core group of tightly correlated patient samples for each subtype . 
 partitioning occurs hierarchically and allows for the identification of genetically similar subtypes . 
 we analyzed the gene expression profiles of 202 tumors of the brain cancer glioblastoma multiforme ( gbm ) given at the cancer genome atlas ( tcga ) site . 
 we identify core patient groups associated with the classical , mesenchymal , and proneural subtypes of gbm . in our analysis , 
 the neural subtype consists of several small groups rather than a single component . 
 a subtype prediction model is introduced which partitions tumors in a manner consistent with clustering algorithms but requires the genetic signature of only 59 genes . 