Chromosomal patterns of gene expression from microarray data: methodology, validation and clinical relevance in gliomas
- Federico E Turkheimer^{1}Email author,
- Federico Roncaroli^{2},
- Benoit Hennuy^{3},
- Christian Herens^{3},
- Minh Nguyen^{4},
- Didier Martin^{4},
- Annick Evrard^{3},
- Vincent Bours^{3},
- Jacques Boniver^{5} and
- Manuel Deprez^{5}
https://doi.org/10.1186/1471-2105-7-526
© Turkheimer et al; licensee BioMed Central Ltd. 2006
Received: 31 October 2006
Accepted: 01 December 2006
Published: 01 December 2006
Abstract
Background
Expression microarrays represent a powerful technique for the simultaneous investigation of thousands of genes. The evidence that genes are not randomly distributed in the genome and that their coordinated expression depends on their position on chromosomes has highlighted the need for mathematical approaches to exploit this dependency for the analysis of expression data-sets.
Results
We have devised a novel mathematical technique (CHROMOWAVE) based on the Haar wavelet transform and applied it to a dataset obtained with the Affymetrix^{®} HG-U133_Plus_2 array in 27 gliomas. CHROMOWAVE generated multi-chromosomal pattern featuring low expression in chromosomes 1p, 4, 9q, 13, 18, and 19q. This pattern was not only statistically robust but also clinically relevant as it was predictive of favourable outcome. This finding was replicated on a data-set independently acquired by another laboratory. FISH analysis indicated that monosomy 1p and 19q was a frequent feature of tumours displaying the CHROMOWAVE pattern but that allelic loss on chromosomes 4, 9q, 13 and 18 was much less common.
Conclusion
The ability to detect expression changes of spatially related genes and to map their position on chromosomes makes CHROMOWAVE a valuable screening method for the identification and display of regional gene expression changes of clinical relevance. In this study, FISH data showed that monosomy was frequently associated with diffuse low gene expression on chromosome 1p and 19q but not on chromosomes 4, 9q, 13 and 18. Comparative genomic hybridisation, allelic polymorphism analysis and methylation studies are in progress in order to identify the various mechanisms involved in this multi-chromosomal expression pattern.
Background
Genes are not randomly distributed and their coordinated expression is regulated by their position along chromosomes [1–3]. New mathematical approaches are therefore needed for the analysis of expression microarrays in order to identify variations in expression of spatially related genes and map them along chromosomes. The relationship between changes in DNA copy number and variations of mRNA expression has been previously investigated (for example [4–7]) but only a few studies have used microarrays to examine large chromosomal abnormalities [7–12].
Here we propose a novel mathematical model based on single value decomposition (SVD) and Haar wavelets, named CHROMOWAVE, that detect variations in expression of spatially related gene and visualise them on chromosomes. Wavelets are a recently introduced mathematical tool for the treatment of signals with "non-stationary behaviour" [13] (e.g a hammer blow, a plane flyover noise etc.). The counterpart of the wavelet transform is the Fourier transform that achieves optimal encoding of periodic signals. The use of wavelets for data encoding, transmission and compression is now pervasive in many fields including analysis of gene sequences and functional genomics data [14]. Application to microarrays has been proposed for the analysis of light signals of microarrays plates [15–17] or to de-noise microarray time-series [18]. Only three studies applied wavelets to explore the variation in expression of gene clusters and identify their position on chromosomes [19–21]. Allen et al. [19] first adopted the wavelet transform and used smooth wavelets to study periodical patterns of mRNA expression elicited by different promoters in the E. coli genome. Using a supervised statistical approach, we introduced the Haar wavelet analysis for the detection of chromosomal patterns of expression in neurodegenerative diseases [20]. The same methodology was validated by Hsu et al. [22] to denoise array-based comparative genomic hybridization (array-CGH) data. Aggarwal et al [21] combined wavelets with an empirical supervised approach to analyze chromosomal expression in a set of tumour cell-lines and matched the extracted clusters with abnormal karyotypes. Their technique was limited as it allowed the analysis of only one cell-line at a time. These two studies demonstrated that wavelets have the ability to identify the areas on chromosomes where genes show similar and coherent levels of expression.
Clinico Pathological Data
Case | Gender | Age | Tumour location | Karnovski score | Histology | Treatment | Follow up (weeks) | Status |
---|---|---|---|---|---|---|---|---|
1 – O1 | M | 27 | right frontal | 80 | Oligodendroglioma grade II | Incomplete surgical resection | 884 | alive and well |
2 – O20 | F | 69 | left temporal | 80 | Oligodendroglioma grade II | subtotal surgical resection, radiotherapy and chemotherapy | 99 | alive |
3 – O4 | F | 50 | right frontal | 80 | Oligodendroglioma grade II | complete surgical resection, radiotherapy | 216 | alive and well |
4 – O23 | F | 59 | left frontal | 80 | Oligodendroglioma grade II | subtotal surgical resection and radiotherapy | 135 | alive and well |
5 – O19 | M | 24 | right frontal | 90 | Mixed oligoastrocytoma grade II | complete surgical resection, radiotherapy | 306 | alive and well |
6 – O15 | M | 56 | left frontal | 60 | Mixed oligoastrocytoma grade II | complete surgical resection, radiotherapy, chemotherapy | 35 | dead of tumour recurrence |
7 – O10 | F | 30 | right frontal | 70 | Mixed oligoastrocytoma grade II | complete surgical resection and radiotherapy | 198 | alive and well |
8 – O7 | M | 32 | left frontal | 80 | Mixed oligoastrocytoma grade II | complete surgical resection, radiotherapy | 643 | alive and well |
9 – O8 | M | 38 | left temporal | 90 | Mixed oligoastrocytoma grade II | subtotal surgical resection, radiotherapy, chemotherapy | 50 | dead of tumour recurrence |
10 – O6 | F | 47 | left temporal | 80 | Mixed oligoastrocytoma grade II | complete surgical resection | 143 | alive and well |
11 – O16 | F | 31 | right parietal | 80 | Astrocytoma grade II | complete surgical resection | 105 | alive and well after 1 recurrence |
12 – AD6 | M | 40 | left temporo-occipital | 80 | Astrocytoma grade II | complete surgical resection | 169 | alive and well |
13 – AD9 | M | 38 | left frontal | 90 | Astrocytoma grade II | complete surgical resection | 136 | recurrence |
14 – AD10 | M | 70 | left temporal | 80 | Astrocytoma grade II | complete surgical resection and radiotherapy | 90 | recurrence (PET Scan) |
15 – AD11 | M | 25 | left temporal | 80 | Astrocytoma grade II | subtotal surgical resection | 137 | recurrence (PET scan) |
16 – AD12 | M | 34 | left temporal | 80 | Astrocytoma grade II | complete surgical resection | 25 | alive and well |
17 – O2 | F | 26 | right frontal | 80 | Anaplastic Oligodendroglioma | complete surgical resection, chemotherapy | 310 | dead of tumour recurrence |
18 – O3 | M | 21 | right frontal | 90 | Anaplastic Oligodendroglioma | subtotal surgical resection and chemotherapy | 867 | recurrence |
19 – O17 | F | 24 | left temporal | 80 | Anaplastic Oligodendroglioma | complete surgical resection and chemotherapy | 578 | alive and well |
20 – O18 | F | 45 | right parietal | 80 | Anaplastic Oligodendroglioma | subtotal surgical resection and chemotherapy | 105 | dead of tumour recurrence |
21 – O9 | F | 46 | left parietal | 80 | Anaplastic oligodendroglioma | complete surgical resection and chemotherapy | 624 | dead of tumour recurrence |
22 – O24 | M | 57 | left fronto-temporal | 80 | Anaplastic Oligoastrocytoma | complete surgical resection and radiotherapy | 15 | alive and well |
23 – O12 | F | 27 | left temporal | 60 | Anaplastic Astrocytoma | complete surgical resection and radiotherapy | 6 | dead of cerebral venous thrombosis |
24 – AA3 | F | 25 | right frontal | 80 | Anaplastic Astrocytoma | complete surgical resection, chemotherapy | 307 | dead of tumour recurrence |
25 – AA6 | F | 36 | left fronto-temporal | 80 | Anaplastic Astrocytoma | subtotal surgical resection and radiotherapy | 447 | alive and well |
26 – AA5 | M | 44 | right temporal | 90 | Anaplastic Astrocytoma | subtotal surgical resection and chemotherapy | 74 | alive and well |
Results
Chromowave and FISH results for the glioma data-set.
Case | Diagnosis | FISH Load Chr1p36 | Chromowave Load Chr1 (+ = loss) | FISH Load Chr4c | Chromowave Load Chr4 (+ = loss) | FISH Load Chr9q34 | Chromowave Load Chr9 (+ = loss) | FISH Load Chr13q14 | Chromowave Load Chr13 (+ = loss) | FISH Load Chr15c | Chromowave Load Chr15 (+ = loss) | FISH Load Chr18c | Chromowave Load Chr18 (+ = loss) | FISH Load Chr19q13 | Chromowave Load Chr19 (+ = loss) | Gender | Chromowave Load ChrY (+ = loss) |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
1 | OII (O1) | Mono | 0.21 | 0.03 | 0.06 | 0.08 | 0.06 | 0.02 | Mono | 0.37 | M | -0.20 | |||||
2 | OII (O20) | Mono | 0.17 | -0.15 | 0.25 | -0.22 | -0.18 | 0.40 | Mono | 0.20 | F | 0.22 | |||||
3 | OII (O4) | Mono | 0.09 | 0.00 | -0.01 | -0.01 | -0.16 | 0.03 | 0.03 | F | 0.14 | ||||||
4 | OII (O23) | 0.04 | Mono | 0.00 | -0.19 | -0.18 | -0.22 | -0.17 | Mono | 0.02 | F | 0.17 | |||||
5 | OAII (O19) | 0.09 | -0.22 | -0.18 | -0.32 | -0.14 | -0.19 | Mono | 0.16 | M | -0.28 | ||||||
6 | OAII (O15) | -0.17 | -0.15 | 0.42 | -0.19 | -0.12 | -0.08 | -0.51 | M | -0.22 | |||||||
7 | OII (O10) | 0.25 | 0.15 | 0.14 | Mono | 0.26 | 0.03 | 0.14 | 0.2 | F | 0.18 | ||||||
8 | OAII (O7) | -0.04 | 0.06 | 0.07 | Mono | 0.30 | Mono | 0.11 | -0.02 | 0.01 | M | -0.24 | |||||
9 | OAII (O8) | -0.19 | -0.12 | Mono | 0.08 | Mono | 0.03 | -0.24 | -0.19 | -0.22 | M | -0.16 | |||||
10 | OAII (O6) | -0.06 | NA | 0.20 | -0.01 | 0.04 | 0.1 | 0.05 | 0.1 | F | 0.16 | ||||||
11 | OII (O16) | -0.08 | -0.12 | 0.12 | -0.08 | -0.01 | 0.06 | 0.05 | F | 0.23 | |||||||
13 | AII (AD6) | -0.08 | -0.01 | -0.12 | -0.12 | 0.03 | -0.02 | -0.12 | M | -0.20 | |||||||
14 | AII (AD9) | -0.12 | -0.04 | -0.07 | -0.04 | 0.01 | -0.07 | 0.07 | M | -0.22 | |||||||
15 | AII (AD10) | -0.15 | -0.12 | -0.23 | -0.22 | -0.1 | -0.17 | -0.13 | M | -0.20 | |||||||
16 | AII (AD11) | Mono | 0. 1 | 0.12 | 0.12 | 0.12 | Mono | 0.27 | 0.09 | 0.01 | M | -0.19 | |||||
17 | AII (AD12) | -0.11 | -0.09 | -0.01 | -0.14 | -0.01 | -0.13 | -0.12 | M | -0.21 | |||||||
18 | OIII (O2) | Mono | 0.09 | 0.28 | Mono | -0.23 | Mono | 0.27 | -0.22 | 0.00 | Mono | 0.03 | F | 0.20 | |||
19 | OIII (O3) | Mono | 0.33 | Mono | 0.64 | 0.14 | Mono | 0.48 | Mono | 0.05 | Mono | 0.68 | Mono | 0.06 | M | 0.01 | |
20 | OIII (O17) | Mono | 0.27 | -0.05 | 0.17 | Mono | -0.10 | Mono | -0.08 | -0.07 | Mono | 0.4 | F | 0.18 | |||
21 | OIII (O18) | Mono | -0.16 | -0.15 | -0.13 | -0.17 | Mono | 0.22 | -0.16 | -0.36 | F | 0.13 | |||||
22 | OIII (O9) | Mono | 0.31 | 0.19 | 0.23 | 0.11 | Mono | 0.7 | 0.09 | Mono | 0.02 | F | 0.19 | ||||
23 | OAIII (O24) | -0.15 | -0.14 | -0.17 | -0.14 | Mono | 0.03 | -0.17 | -0.19 | M | -0.22 | ||||||
24 | OIII (O12) | -0.17 | 0.29 | -0.04 | 0.30 | 0.25 | 0.19 | 0.03 | F | 0.13 | |||||||
25 | AIII (AA3) | -0.15 | -0.10 | -0.48 | 0.02 | -0.11 | -0.26 | -0.1 | F | 0.25 | |||||||
26 | AIII (AA6) | -0.07 | -0.04 | 0.32 | 0.15 | -0.11 | 0.09 | -0.03 | F | 0.15 | |||||||
27 | AIII (AA5) | -0.05 | -0.30 | -0.10 | -0.10 | -0.08 | -0.04 | 0.09 | M | -0.18 |
We then compared FISH measurements to the individual chromosomal expressions extracted by CHROMOWAVE. When applied to one chromosome at a time, the SVD extracted as main components (> 70% of the total variability) chromosome wide diffuse signals on chromosomes 4, 9, 13, 15 and 18 and diffuse homogeneous expression on the chromosomal arms 1p and 19q. Case loadings for these patterns are displayed in Table 2. Note that the case-loadings in Table 2 reflect the normalization to the average that is performed by the SVD. Therefore, the more positive the case loading is the more the pattern is expressed (in this case the bigger the loss), the more negative the loading the less the pattern is expressed. When we compared the FISH measurements with these case-loadings we observed significant association for chromosome 1p (Pearson correlation R = 0.522, p = 0.005) and 19q (Pearson correlation R = 0.392, p = 0.043) indicating that loss of genetic material was the main cause of the reduction of expression seen by CHROMOWAVE.
Interestingly, one oligodendroglioma (case 010, see figure 2) clearly demonstrated low expression 1p/19q with CHROMOWAVE but no structural loss could be detected by FISH (Table 2). The reduced mRNA expression detected by CHROMOWAVE on chromosomes 4, 9q, 13, 15 and 18 was less commonly associated with chromosomal alterations seen by FISH (Table 2). In particular, six tumours showed low expression of 9q without alteration detected by FISH and among the seven lesions demonstrating low expression in chromosome 18, only one showed monosomy.
The discrepancy between FISH counts and expression data raises several hypotheses. First, this may be inherent to the FISH method where the probe targets only a short DNA sequence on the chromosome and is not informative of possible large losses of genetic material in regions flanking the probe target. The chromosomal areas targeted by the FISH probes are shown on Figure 1. Second, alternative genetic and epigenetic mechanisms can cause expression changes in adjacent genes in the absence of chromosomal loss such as translocation, uniparental disomy or methylation/acetylation silencing, all frequently reported in cancer. Third, hyperploidy which is frequently seen in malignant gliomas may also account for some of these observations.
When we tested the case loadings identified by CHROMOWAVE with outcome (tumour recurrence and patient survival) using Cox regression, we found that the pattern in Figure 1 was significantly predictive of favourable outcome (p = 0.007). By testing each of the chromosomes and their various combinations we then observed that the covariation of chromosomes 1p, 13 and 18had the strongest correlation with survival (p = 0.002).
In contrast, the major gene expression pattern obtained with the same SVD analysis but without Haar wavelet transformation did not correlate with survival (p = 0.802) suggesting that the distribution of gene expression changes on chromosomes is more relevant to tumour behaviour than their raw variations of amplitude.
Sensitivity to Individual Cases
In order to verify the stability of the pattern and its dependency upon single cases, we performed a jack-knife test to verify that association remained significant with the exclusion of single cases. The genome-wide SVD analysis was repeated after removal of one case (27 iterations), the chromosomal pattern was extracted and stored, and the correlation of the resulting case loadings with survival re-calculated using Cox regression, This process generated 27 p-values. The 95% confidence interval of the empirical p-value distribution obtained was calculated by normal approximation of the log-transformed p-values. Throughout the 27 permutations, the global patterns recovered were indistinguishable from that obtained with all cases, all the 27 patterns were significantly associated with outcome (p < 0.05) and resulting p-values were tightly distributed (p-values 95% confidence interval was [0.0372–0.0022]).
Sensitivity to Chromosome Biology: Chromosome Y and Gender Sensitivity
Sensitivity to Denoising Parameters
Sensitivity/Specificity Analysis
Although the technique presented here is unsupervised, its sensitivity/specificity can be evaluated through simulation studies by assuming the signal distribution known. The main problem with simulations in this context is the faithful generation of the noise covariance structure of chromosomal expression that is unknown. To recover the noise covariance we have selected Chr1 that has evident Chr1p loss pattern in this data-set. We have removed from the data-set this specific monosomy by zeroing the first singular component. The remaining singular components had Morgera's Complexity ~1 [24] indicating that all that was left was noise. We have therefore built a simulation by adding a telomeric pattern of various intensities (when signal was 0 the specificity was obtained) and of varying spatial dimensions (500 Kb, 1.5 MB, the whole "petit" arm and the whole chromosome) to 13 of the 27 arrays. In the second simulation we maintained a Chr1p loss pattern but varied the number of arrays to which the pattern was added.
Clinical Reproducibility: Application to Freije et al. (2004) dataset
Final validation was performed by assessing the reproducibility of the association between chromosomes 1p, 13 and 18 and outcome in a published set of arrays of comparable tumour types. The only comparable data-set that was publicly available at the time of writing was published by Freije et al. [23]. Microarray files (Affymetrix HG U133A and B oligonucleotide arrays) and clinical data were downloaded from the authors' website [25]. We excluded glioblastomas from the set and we examined the 25 arrays that included 10 anaplastic oligodendrogliomas, 8 anaplastic astrocytomas and 7 mixed anaplastic oligo-astrocytomas. Microarray files were pre-processed and normalized using the same procedure adopted for our data and then entered into CHROMOWAVE. Chromosomal expressions were obtained by application of the SVD to each chromosome independently. As in our data-set, diffuse patterns of expression were found in chromosomes 1p, 4, 9, 13, 18 and 19q. Case loadings for chromosomes 1p, 13 and 18 were then entered into a Cox regression model to test their association with survival that resulted being remarkably strong (p = 0.0028) and similar to what found in our data-set.
Discussion
CHROMOWAVE allows the unsupervised identification of clusters of adjacent genes with homogeneous changes of expression and their mapping on chromosomes, resulting in the display of multi-chromosomal gene expression patterns. Here we have demonstrated that these patterns are reliable, reproducible and statistically robust but that they are also clinically relevant. In this application, the SVD was the method of choice for statistical analysis in wavelet space: however this approach is amenable to treatment with any other unsupervised technique (independent component analysis, clustering etc.).
Low grade and anaplastic diffuse gliomas represent an interesting model to explore the variation of expression of spatially related genes rather than their individual expression. Characteristic genetic and epigenetic alterations have been found that predict favourable outcome in some tumour histological subtypes while being non informative or even carrying a poor prognosis in others. For instance, the extensively investigated allelic loss of chromosomes 1p/19q or the hypermethylation of the promoter of the drug-resistance gene O6-methylguanine-DNA methyltranferase are associated with chemosensitivity and therefore longer survival in a subset of oligodendrogliomas [26–30] but have unclear prognostic value or even correlate with more aggressive behaviour in astrocytomas and mixed oligoastrocytomas [28]. Moreover, no strong candidate tumour suppressor/promoter gene has been singled out yet on these frequently lost chromosomal segments [31]. A few studies have used expression microarrays to investigate diffuse gliomas [9, 23, 32–41] but none explored the expression change of gene clusters that could be relevant to tumour progression with regards to their distribution on chromosomes.
In our microarray dataset, derived for a group of 27 WHO grade II and III gliomas, CHROMOWAVE generated a multi-chromosomal pattern of variation that correlated with outcome. The pattern included diffuse losses on 1p and 19q but also on 4, 9q, 13 and 18. Among these, changes on 1p, 13 and 18 had the strongest correlation with survival. This finding was replicated on a comparable set of microarray data previously published by another laboratory [23].
Notably, while the main pattern of chromosomal variation extracted by the SVD correlated with survival, the main pattern of RNA variation extracted by the SVD from the raw data (no spatial transform applied) did not. This suggests a major role of chromosomal RNA modulation in tumour behaviour as opposed to the variation of the bulk of gene expression.
FISH studies suggested that, in our dataset, low expression on chromosomes 1p and 19q were most often the consequence of large allelic loss in these regions, an alteration commonly seen in oligodendrogliomas. However, FISH counts could not explain the diffusely reduced expression on 4, 9q, 13 and 18 raising several hypotheses: genetic loss occurring in regions flanking the FISH probe targets, genetic changes that do not result in gene loss detectable by FISH such as translocation and uniparental disomy, or epigenetic alterations such as methylation-based gene silencing. Clearly, extensive ancillary studies are needed to determine the various mechanisms underlying the CHROMOWAVE gene expression patterns.
Comparative genomic hybridisation (CGH), allelic polymorphism analysis and methylation studies are currently in progress in our laboratory. Whatever the causative mechanisms, the finding that large gene expression changes in chromosomes 4, 9q, 13 and 18 occur frequently in grade II and III diffuse gliomas and that they bear prognostic information is novel.
Conclusion
In conclusion, we propose a new mathematical model that has proved powerful in our dataset for detecting and mapping to chromosomes biologically meaningful gene expression changes. The possibility of visualising changes of spatially related genes and their position on chromosomes make CHROMOWAVE a valuable screening method to explore microarray datasets. The mechanisms contributing to these expression patterns are probably multiple and complex. Additional studies combining FISH, CGH/aCGH, allelic polymorphism and methylation analysis are clearly needed and should target those chromosomal areas identified by CHROMOWAVE as supporting clinically relevant gene expression changes.
Methods
We studied a dataset generated with Affymetrix U133_Plus_2 arrays (Affymetrix, Santa Clara, CA) in 27 low grade and anaplastic diffuse gliomas (clinicopathologic features are summarised in Table 1) and 11 samples of normal brain obtained in course of surgery for intractable epilepsy. Tissues were collected under the approved guidelines of the Ethics Committee of the Faculty of Medicine, University of Liège, Belgium and all patients gave informed consent for their participation in this study.
RNA extraction, target preparation and microarray hybridisation
Total RNA was extracted from cryostatic sections using the Qiagen RNeasy kit (Qiagen, Chatsworth, CA). The integrity of the RNA was confirmed with the Agilent Bioanalyser using the RNA 6000 Nano kit (Agilent). We used the GeneChip^{®} Expression 3' Amplification One-Cycle Target Labeling kit (Affymetrix, Santa Clara, CA) to label the RNA following the manufacturer protocol. The cRNA was hybridized to Affymetrix Human U133_Plus_2 arrays according to the manufacturer protocol. Briefly, double-stranded cDNA was synthesized routinely from five micrograms of total RNA primed with a poly-(dT) -T7 oligonucleotide. The cDNA was used in an in vitro transcription reaction (IVT) in the presence of T7 RNA polymerase and biotin-labelled modified nucleotides during 16 hours at 37°C. Biotinylated cRNA was purified and then fragmented (35–200 nucleotides), together with hybridization controls and hybridized to the microarrays for 16 h at 45°C. Using the Fluidics Station (Affymetrix), the biotin-labeled cRNA was revealed by successive reactions with streptavidin R-phycoerythrin conjugate, biotinylated antistreptavidine antibody and streptavidin R-phycoerythrin conjugate. The arrays were finally scanned in an Affymetrix/Hewlett-Packard GeneChip Scanner 3000
Preliminary data analysis
Preliminary data analysis was conducted using the software of the Affymetrix microarray suite (MAS, version 5.0) following the statistical procedures described in the Affymetrix: Statistical Algorithms Detection Document [42]). MAS produced an expression value plus an index parameter indicating positive or negative detection (present call index) for each of the 54,675 probe sets on the chip (settings used were standard for the U133_Plus_2 array: alpha1 = 0.05, alpha2 = 0.065, Tau = 0.015, TGT = 100). Statistical analysis and post-processing were performed using an in-house software (CHROMOWAVE) written in MATLAB 6.5 (The Mathworks Inc., Natick MA, USA). Individual arrays were normalized to the background by dividing intensities by the median value of those genes presented with positive detection. Expression values where then log2 transformed.
Mapping Target Sequence Values to Chromosomal Location
Expression values were mapped to their corresponding chromosomal location and then sorted within each vector using genome alignment information. Information on the physical location of each gene and the respective genome alignment information for each target sequence on the HG-U133_Plus_2 chip were obtained from the Affymetrix website [43].
Haar Wavelet Analysis of Chromosomal Expression
Gene expression values were analysed through CHROMOWAVE that uses the positional information of genes and statistical analysis to extract chromosomal pattern of gene expression. CHROMOWAVE applies the wavelet transform (WT) to the spatial distribution of the array probes and converts the original expression values in wavelet coefficients that are functions of the expression of adjacent genes. Wavelet coefficients are then filtered so that only those with high signal-to-noise ratio and/or representing probes with close genomic distance are retained. The application of the inverse WT produces a noise-free pattern of chromosomal gene expression. When multiple arrays are used, wavelet coefficients can be used in statistical analysis in the same fashion as gene expression values either with supervised methods of analysis (t-tests, ANOVA, discriminant analysis, etc.) or unsupervised (clustering techniques, independent component analysis etc.). The WT algorithm has been described elsewhere (see for example [44]) and will be only summarized here. The traditional WT scheme is limited by the decimation step that may "miss" relevant signal elements, particularly when noise levels are high as in this application. For this reason, CHROMOWAVE adopts the "cycle spinning" WT that has greater complexity but enjoys the translation-invariant property [45]. Chromosomal patterns of gene expression are not expected to be smooth but to have well defined boundaries. Therefore the WT transform adopts the simplest wavelet, the classic Haar wavelet [46], that allows a constant piece-wise approximation of the RNA profile. As a result, the WT can be described as follows. Firstly, gene expression values obtained from microarray measurements were sorted according to their chromosomal location as previously described. Then the coefficients of the first level of the WT were calculated as the difference of expression between two adjacent probes. The wavelet coefficients of the second level were obtained as the difference between the mean of pairs of adjacent probes. The coefficients of the next levels were then calculated as the differences between the means of P adjacent probes where P increases as a power of 2 (P = 1, 2, 4, 8, ...). The WT is an orthogonal operator and therefore the noise level is identical on the original raw data and at all WT levels. In contrast, clusters of genes with similar level of expression produce a WT coefficient that increases with the resolution level. In other words, genes whose individual expression is undetectable because they are below the noise level are detected through the WT when clustered together because their combined energy condenses into a greater wavelet coefficient.
Note that standard WT is usually applied to equally spaced data, which is not the case here because genes are not equally spaced on chromosomes. Sardy et al. [47] have showed that a wavelet estimator based on the Haar wavelet transform provides an estimate that is at least as good as that recovered by any other Haar wavelet implementation adapted to the unequally spaced case.
Unsupervised Analysis
Unsupervised analysis was performed using the Singular Value Decomposition algorithm [48] (SVD). SVD was applied to the set of Haar wavelet coefficients [49]. This produced a number of patterns of chromosomal expression equal to the number of the arrays. The contribution of each case to each pattern was calculated as a single number, the "case loading." Case loadings quantify the amount of the pattern expressed by each array. Case loadings were then used for further statistical analysis. Following sections describe the method in detail.:
Notation and Haar Wavelet Transform
Let C(h,k) be the matrix for the k = 1, 2, .., M arrays forming the experiment containing on the positions h = 1, 2, ..., 2^{n} the ordered probes. Note that, for algorithmic reasons, the number of probes must be a power of 2 and, if this is not the case, the matrix must be zero-padded. Application of the Haar WT to each column generates n levels of 2^{n} wavelet coefficients that are serially stored in the matrix C_{W}(i,k) where i = 1, 2, ..., n2^{n}.
Individual vs. Global Chromosomal Analysis
Matrix C_{W} so far contains the Haar wavelet decomposition for M arrays and one chromosome only. This allows the analysis of each chromosome independently. This option was used in this work for the comparison between individual chromosomal expression and structural changes detected by FISH. However, the core of the work was the analysis of the entire genome simultaneously to detect the combination of chromosomal patterns characteristic of this data-set and, possibly, of gliomas at large. The analysis can be extended to the whole genome by serially adding to the rows of C_{W} the Haar wavelet coefficients of all the other chromosomes. The following section is valid for both types of analysis and notation will simply refer to C_{W} irrespective of whether it contains the WT of a single chromosome or of all chromosomes.
SVD in Wavelet Space
Calculation of SVD for C_{W} cannot be solved directly because of the large size of the covariance matrix [C_{W}*${\text{C}}_{\text{W}}^{\text{T}}$] that has size n2^{n} × n2^{n}. This problem can be circumvented by using instead the covariance Cov(C_{W}) = ${\text{C}}_{\text{W}}^{\text{T}}$ * C_{W} as follows [48].
Matrix C_{W} is firstly normalized by removing the row means and then the SVD is applied to produce the decomposition:
Cov(C_{W}) = ${\text{C}}_{\text{W}}^{\text{T}}$ * C_{W} = V_{W}S_{W}${\text{V}}_{\text{W}}^{\text{T}}$ (1)
V_{W} is an MxM matrix. Each column [V_{W}]_{i}, or "singular vector", contains one of the M directions of maximal change for the arrays. These directions are orthogonal to each other and each is made of M coefficients representing the contribution of each array to that particular direction. We label these coefficients as "case loadings."
Case loadings can be used in any type of statistical analysis. For example, they can be entered in a bivariate correlation to test the association of the corresponding chromosomal pattern with clinical parameters by associating the load of each array with the correspondent external measure (say, survival of that particular patient).
S_{W} is an MxM diagonal matrix with diagonal elements S_{ii} = S_{i}. The M diagonal elements S_{i} are the singular values of C_{W} and, without loss of generality, it can be assumed that they are ordered in decreasing order so that S_{1} ≥ S_{2} ≥ ... ≥ S_{M}. The fraction of total variability in the expression data-set explained by any individual column [V_{W}]_{i} can be calculated as
f[V_{W}]_{i} = S_{i}/Σ (S_{i}) (2)
This means that the first singular vector explains the greatest amount of data variability; the second singular vector contains the direction of change with the second greatest variance and so forth. It is expected that the first singular vectors contain variability due to "real" (biological) signal, while noise contributions will be contained in the last ones [48].
The Haar wavelet patterns corresponding to the singular vectors can be recovered by projecting the matrix C_{W} on the rotated axis V_{W}.
U_{W} = C_{W} * V (3)
Similarly as before, the first columns [U_{W}]_{i}, corresponding to the first singular vectors, should mostly contain true signal. However a further reduction in the noise can be achieved by removing from matrix U_{W} all coefficients smaller than a suitable threshold.
In CHROMOWAVE the filtering procedure incorporates both the noise reduction and the introduction in the model of inter-probe genomic distances according to the procedure described in the next section.
Inter-Probes Distance and Noise Penalization
Previous wavelet models of chromosomal gene distribution [20, 21] were based on the simplifying assumption that the relation between adjacent probes, if present, does not depend on the absolute physical distance (as measured in base-pairs) but only on contiguity. This assumption may introduce inaccuracies, particularly for local processes involving a very small number of genes because wavelet coefficients (level 2 and above) pool together the expression of probes or groups of probes that may have quite varying distances among them.
In CHROMOWAVE, the likelihood of a Haar wavelet is made proportional to the distance between probes or probe-groups that the wavelet coefficient represents. Inter-probe likelihood was modeled in CROMOWAVE by adding a penalty function to the de-noising procedure.
In wavelet analysis, de-noising is achieved by suppression of all wavelet coefficients that are below an appropriate threshold dependent on the noise levels in the data. This operation requires an estimate of the noise variance of the data. The variances of gene expression measured with microarrays are usually heterogeneous, e.g, vary from gene to gene. However CHROMOWAVE aims to detect gene clusters only; therefore individual gene expressions which correspond to the first level of the WT are of no interest and are suppressed. The other levels of the WT are all generated by pooling together 2,4 8, ..., gene expression values and their variances are, therefore, more homogeneous. Besides, since the WT is an orthogonal operator, all Haar wavelet levels have approximately the same variance that can then be calculated from the robust estimator:
MAD denotes median absolute deviation from 0 and the factor 0.6745 is chosen for calibration with the normal distribution [50].
Spatial modeling and de-noising is therefore achieved by suppressing all coefficients in the matrix U_{W} that are under the threshold:
P(w) is a penalty of the form:
P(w) = 1-G(ln(d),μ,ν) (6)
In (6), d is the genomic distance between genes or gene groups represented by the Haar wavelet coefficient w and G is the Gaussian cumulative distribution with mean μ and standard deviation ν. This form of the penalty is justified by the fact that inter-wavelet genomic distribution, that we pre-calculated using the available information on the location of the probes, was Gaussian-like for all wavelet levels. Parameters μ and ν were directly obtained by the genomic-alignment information for the HG-U133_Plus_2 chip.
Confoundings Due to Aneuploidy and Errors in Normalisation
When applied to the whole chromosome set, the combination of SVD analysis and wavelets has the additional practical utility of identifying errors in the normalisation of data (if linear). With CHROMOWAVE, the application of an inefficient normalization procedure results in a genome-wide constant chromosomal pattern of expression that the SVD identifies and removes from the data. Besides, note that global RNA changes due to aneuploidy also result in the same genome-wide diffuse pattern that can be seemingly removed from the overall data-variability.
Chromosomal Pattern Reconstruction
Threshold (5) suppresses all those coefficients that are unlikely to be signal because of their relative height compared to noise and/or because they contain probes/probes clusters that are far apart. All those coefficients surviving the threshold (5) instead are likely contributions to the true signal and are passed through the inverse WT to produce the M filtered patterns of chromosomal variations C_{F}(h,k) where h = 1, 2, ..., 2^{n} and k = 1, 2, ..., M.
Single Profile Generation (Turkheimer et al., 2004 revisited)
CHROMOWAVE allows also the extraction of the differential profile of expression between a single case and a control group (supervised analysis). This application was described previously [20] and is just summarized here. For each chromosome, the differential profile is defined as:
d C(i) = C(i) - $\overline{\text{C}}$(i) (7)
C(i) contains on the positions i = 1, 2, ..., 2^{n} the expressions of the ordered probes for the single case of interest. As before, the number of probes must be a power of 2 and, if this is not the case, the matrix must be zero-padded.
d C_{W}(i,1) is the finest resolution level of the wavelet transform.
Upon application of the inverse wavelet transform to the filtered matrix d C_{W}(i,j) one obtains a de-noised approximation of the individual pattern d C(i).
We used CHROMOWAVE to extract patterns of chromosomal expression for individual tumour cases by contrasting its microarray measurement with the average expression of a normal database of 11 normal brain samples. RNA extracted from these cases was hybridized to Affymetrix U133_Plus_2 arrays and data processed as described previously.
Fluorescence in situ hybridization (FISH)
The 27 tumours were all also studied with FISH. Dual-colour assays were performed on 8-μm-thick cryostatic sections from the same tissue blocks as those used for microarray experiments. We tested the six chromosomes that revealed the major changes with CHROMOWAVE. Loss or gain on chromosomes 1 and 19 were detected with LSI^{®}1p36/LSI 1q25 and LSI 19q13/LSI 19p13 dual-color probe sets (Vysis, Inc., Downers Grove, IL, USA). Chromosome 9q was studied with the ABL (9q34) probe of the LSI ABL/BCR ES probe system (Vysis). For chromosomes 4, 13, 15 and 18, we used Vysis probes CEP4, LSI13 (440 kb including the RB gene in 13q14), CEP15 and CEP18. Samples were processed according to the manufacturer's protocol. Results were evaluated with an Olympus BX51 fluorescence microscope (Olympus, Melville, NY, USA) equipped with the appropriate fluorescence filters. At least 100 cells were examined for all signals and the mean signal numbers were recorded. Frozen sections of normal brain were analysed to establish a reference FISH copy number. Upper and Lower normality thresholds were calculated as mean +/- 2 standard deviations (SD). Tumour samples with mean signal below the lower threshold were reported as showing monosomy. Correlation between FISH mean copy numbers and CHROMOWAVE loadings were calculated using Pearson product moment correlation coefficient.
Availability and Requirements
The microarray data used in this work are deposited in the GEO database (GEO Submission GSE2817). Software available on request, free for academic users.
Project name: CHROMOWAVE
Project home page: http://www1.imperial.ac.uk/medicine/people/federico.turkheimer
Operating system(s): Platform independent;
Programming language: Matlab 7 (R14)
Licence: GNU
Non-academics: licence needed.
Declarations
Acknowledgements
The authors wish to express their gratitude to Dr Nelson S.F. for making his microarray dataset available for testing. This work was supported by the Télévie 2002 (grant n°7.4580.02), the Fonds National de la Recherche Scientifique (grant 1.5.286.05) and by a grant from the BTRC-Way Ahead Charity.
Authors’ Affiliations
References
- Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals. Nat Genet 2003, 33 Suppl: 245–254. 10.1038/ng1089View ArticlePubMedGoogle Scholar
- Vogelstein B, Kinzler KW: Cancer genes and the pathways they control. Nat Med 2004, 10: 789–799. 10.1038/nm1087View ArticlePubMedGoogle Scholar
- Sproul D, Gilbert N, Bickmore WA: The role of chromatin structure in regulating the expression of clustered genes. Nat Rev Genet 2005, 6: 775–781. 10.1038/nrg1688View ArticlePubMedGoogle Scholar
- Pollack JR, Sorlie T, Perou CM, Rees CA, Jeffrey SS, Lonning PE, Tibshirani R, Botstein D, Borresen-Dale AL, Brown PO: Microarray analysis reveals a major direct role of DNA copy number alteration in the transcriptional program of human breast tumors. Proc Natl Acad Sci U S A 2002, 99: 12963–12968. 10.1073/pnas.162471999PubMed CentralView ArticlePubMedGoogle Scholar
- Myers CL, Dunham MJ, Kung SY, Troyanskaya OG: Accurate detection of aneuploidies in array CGH and gene expression microarray data. Bioinformatics 2004, 20: 3533–3543. 10.1093/bioinformatics/bth440View ArticlePubMedGoogle Scholar
- Masayesva BG, Ha P, Garrett-Mayer E, Pilkington T, Mao R, Pevsner J, Speed T, Benoit N, Moon CS, Sidransky D, Westra WH, Califano J: Gene expression alterations over large chromosomal regions in cancers include multiple genes unrelated to malignant progression. Proc Natl Acad Sci U S A 2004, 101: 8715–8720. 10.1073/pnas.0400027101PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou Y, Luoh SM, Zhang Y, Watanabe C, Wu TD, Ostland M, Wood WI, Zhang Z: Genome-wide identification of chromosomal regions of increased tumor expression by transcriptome analysis. Cancer Res 2003, 63: 5781–5784.PubMedGoogle Scholar
- FitzPatrick DR, Ramsay J, McGill NI, Shade M, Carothers AD, Hastie ND: Transcriptome analysis of human autosomal trisomy. Hum Mol Genet 2002, 11: 3249–3256. 10.1093/hmg/11.26.3249View ArticlePubMedGoogle Scholar
- Mukasa A, Ueki K, Matsumoto S, Tsutsumi S, Nishikawa R, Fujimaki T, Asai A, Kirino T, Aburatani H: Distinction in gene expression profiles of oligodendrogliomas with and without allelic loss of 1p. Oncogene 2002, 21: 3961–3968. 10.1038/sj.onc.1205495View ArticlePubMedGoogle Scholar
- Kano M, Nishimura K, Ishikawa S, Tsutsumi S, Hirota K, Hirose M, Aburatani H: Expression imbalance map: a new visualization method for detection of mRNA expression imbalance regions. Physiol Genomics 2003, 13: 31–46.View ArticlePubMedGoogle Scholar
- Husing J, Zeschnigk M, Boes T, Jockel KH: Combining DNA expression with positional information to detect functional silencing of chromosomal regions. Bioinformatics 2003, 19: 2335–2342. 10.1093/bioinformatics/btg314View ArticlePubMedGoogle Scholar
- Zhou X, Cole SW, Hu S, Wong DT: Detection of DNA copy number abnormality by microarray expression analysis. Hum Genet 2004, 114: 464–467. 10.1007/s00439-004-1087-9View ArticlePubMedGoogle Scholar
- Mallat SG: A theory of multiresolution signal decomposition: the wavelet representation. IEEE Transactions on Pattern Analysis and Machine Intelligence 1989, 11: 673–693. 10.1109/34.192463View ArticleGoogle Scholar
- Lio P: Wavelets in bioinformatics and computational biology: state of art and perspectives. Bioinformatics 2003, 19: 2–9. 10.1093/bioinformatics/19.1.2View ArticlePubMedGoogle Scholar
- Wang J, Ma JZ, Li MD: Normalization of cDNA microarray data using wavelet regressions. Comb Chem High Throughput Screen 2004, 7: 783–791. 10.2174/1386207043328274View ArticlePubMedGoogle Scholar
- Wang XH, Istepanian RS, Song YH: Microarray image enhancement by denoising using stationary wavelet transform. IEEE Trans Nanobioscience 2003, 2: 184–189. 10.1109/TNB.2003.816225View ArticlePubMedGoogle Scholar
- Wang J, Meza-Zepeda LA, Kresse SH, Myklebost O: M-CGH: analysing microarray-based CGH experiments. BMC Bioinformatics 2004, 5: 74. 10.1186/1471-2105-5-74PubMed CentralView ArticlePubMedGoogle Scholar
- Klevecz RR: Dynamic architecture of the yeast cell cycle uncovered by wavelet decomposition of expression microarray data. Funct Integr Genomics 2000, 1: 186–192. 10.1007/s101420000027View ArticlePubMedGoogle Scholar
- Allen TE, Herrgard MJ, Liu M, Qiu Y, Glasner JD, Blattner FR, Palsson BO: Genome-scale analysis of the uses of the Escherichia coli genome: model-driven analysis of heterogeneous data sets. J Bacteriol 2003, 185: 6392–6399. 10.1128/JB.185.21.6392-6399.2003PubMed CentralView ArticlePubMedGoogle Scholar
- Turkheimer FE, Duke DC, Moran LB, Graeber MB: Wavelet analysis of Gene Expression. Volume Proceedings of the 2004 IEEE International Symposium on Biomedical Imaging. Arlington VA; 2004:1183–1186.Google Scholar
- Aggarwal A, Leong SH, Lee C, Kon OL, Tan P: Wavelet transformations of tumor expression profiles reveals a pervasive genome-wide imprinting of aneuploidy on the cancer transcriptome. Cancer Res 2005, 65: 186–194. 10.1158/0008-5472.CAN-05-1036View ArticlePubMedGoogle Scholar
- Hsu L, Self SG, Grove D, Randolph T, Wang K, Delrow JJ, Loo L, Porter P: Denoising array-based comparative genomic hybridization data using wavelets. Biostatistics 2005, 6: 211–226. 10.1093/biostatistics/kxi004View ArticlePubMedGoogle Scholar
- Freije WA, Castro-Vargas FE, Fang Z, Horvath S, Cloughesy T, Liau LM, Mischel PS, Nelson SF: Gene expression profiling of gliomas strongly predicts survival. Cancer Res 2004, 64: 6503–6510. 10.1158/0008-5472.CAN-04-0452View ArticlePubMedGoogle Scholar
- Morgera SD: Information theoretic complexity and relation to pattern recognition. IEEE Transactions in Systems Man and Cybernetics 1985, 15: 608–619.View ArticleGoogle Scholar
- Nelson DF: Freije Cancer Research.2004. [http://sumo.genetics.ucla.edu/~snelson/PublicDATASETS/Freije_CancerResearch_2004]Google Scholar
- Fuller CE, Perry A: Molecular diagnostics in central nervous system tumors. Adv Anat Pathol 2005, 12: 180–194. 10.1097/01.pap.0000175117.47918.f7View ArticlePubMedGoogle Scholar
- Mischel PS, Cloughesy TF, Nelson SF: DNA-microarray analysis of brain cancer: molecular classification for therapy. Nat Rev Neurosci 2004, 5: 782–792. 10.1038/nrn1518View ArticlePubMedGoogle Scholar
- McDonald JM, See SJ, Tremont IW, Colman H, Gilbert MR, Groves M, Burger PC, Louis DN, Giannini C, Fuller G, Passe S, Blair H, Jenkins RB, Yang H, Ledoux A, Aaron J, Tipnis U, Zhang W, Hess K, Aldape K: The prognostic impact of histology and 1p/19q status in anaplastic oligodendroglial tumors. Cancer 2005, 104: 1468–1477. 10.1002/cncr.21338View ArticlePubMedGoogle Scholar
- Papagikos MA, Shaw EG, Stieber VW: Lessons learned from randomised clinical trials in adult low grade glioma. Lancet Oncol 2005, 6: 240–244. 10.1016/S1470-2045(05)70095-4View ArticlePubMedGoogle Scholar
- Nutt CL, Noble M, Chambers AF, Cairncross JG: Differential expression of drug resistance genes and chemosensitivity in glial cell lineages correlate with differential response of oligodendrogliomas and astrocytomas to chemotherapy. Cancer Res 2000, 60: 4812–4818.PubMedGoogle Scholar
- Reifenberger G, Louis DN: Oligodendroglioma: toward molecular definitions in diagnostic neuro-oncology. J Neuropathol Exp Neurol 2003, 62: 111–126.PubMedGoogle Scholar
- Mukasa A, Ueki K, Ge X, Ishikawa S, Ide T, Fujimaki T, Nishikawa R, Asai A, Kirino T, Aburatani H: Selective expression of a subset of neuronal genes in oligodendroglioma with chromosome 1p loss. Brain Pathol 2004, 14: 34–42.View ArticlePubMedGoogle Scholar
- Rickman DS, Bobek MP, Misek DE, Kuick R, Blaivas M, Kurnit DM, Taylor J, Hanash SM: Distinctive molecular profiles of high-grade and low-grade gliomas based on oligonucleotide microarray analysis. Cancer Res 2001, 61: 6885–6891.PubMedGoogle Scholar
- Watson MA, Perry A, Budhjara V, Hicks C, Shannon WD, Rich KM: Gene expression profiling with oligonucleotide microarrays distinguishes World Health Organization grade of oligodendrogliomas. Cancer Res 2001, 61: 1825–1829.PubMedGoogle Scholar
- van den Boom J, Wolter M, Kuick R, Misek DE, Youkilis AS, Wechsler DS, Sommer C, Reifenberger G, Hanash SM: Characterization of gene expression profiles associated with glioma progression using oligonucleotide-based microarray analysis and real-time reverse transcription-polymerase chain reaction. Am J Pathol 2003, 163: 1033–1043.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoelzinger DB, Mariani L, Weis J, Woyke T, Berens TJ, McDonough WS, Sloan A, Coons SW, Berens ME: Gene expression profile of glioblastoma multiforme invasive phenotype points to new therapeutic targets. Neoplasia 2005, 7: 7–16. 10.1593/neo.04535PubMed CentralView ArticlePubMedGoogle Scholar
- Shai R, Shi T, Kremen TJ, Horvath S, Liau LM, Cloughesy TF, Mischel PS, Nelson SF: Gene expression profiling identifies molecular subtypes of gliomas. Oncogene 2003, 22: 4918–4923. 10.1038/sj.onc.1206753View ArticlePubMedGoogle Scholar
- Sallinen SL, Sallinen PK, Haapasalo HK, Helin HJ, Helen PT, Schraml P, Kallioniemi OP, Kononen J: Identification of differentially expressed genes in human gliomas by DNA microarray and tissue chip techniques. Cancer Res 2000, 60: 6617–6622.PubMedGoogle Scholar
- Godard S, Getz G, Delorenzi M, Farmer P, Kobayashi H, Desbaillets I, Nozaki M, Diserens AC, Hamou MF, Dietrich PY, Regli L, Janzer RC, Bucher P, Stupp R, de Tribolet N, Domany E, Hegi ME: Classification of human astrocytic gliomas on the basis of gene expression: a correlated group of genes with angiogenic activity emerges as a strong predictor of subtypes. Cancer Res 2003, 63: 6613–6625.PubMedGoogle Scholar
- Nutt CL, Mani DR, Betensky RA, Tamayo P, Cairncross JG, Ladd C, Pohl U, Hartmann C, McLaughlin ME, Batchelor TT, Black PM, von Deimling A, Pomeroy SL, Golub TR, Louis DN: Gene expression-based classification of malignant gliomas correlates better with survival than histological classification. Cancer Res 2003, 63: 1602–1607.PubMedGoogle Scholar
- Mischel PS, Shai R, Shi T, Horvath S, Lu KV, Choe G, Seligson D, Kremen TJ, Palotie A, Liau LM, Cloughesy TF, Nelson SF: Identification of molecular subtypes of glioblastoma by gene expression profiling. Oncogene 2003, 22: 2361–2373. 10.1038/sj.onc.1206344View ArticlePubMedGoogle Scholar
- Affymetrix: Statistical Algorithms Detection Document 2002.Google Scholar
- HG-U133_Plus_2 Data Sheet 2004.Google Scholar
- Mallat SG: A Wavelet Tour of Signal Processing. 2nd edition. San Diego, Academic Press; 1999.Google Scholar
- Coifman RR, Donoho DL: Translation-invariant denoising. Lecture Notes in Statistics 1995, 103: 125–150.View ArticleGoogle Scholar
- Haar A: Zur Theorie der orthogonalen Funktionensysteme. Annals of Mathematics 1910, 69: 331–371. 10.1007/BF01456326View ArticleGoogle Scholar
- Sardi S, Percival D, Bruce A, Gao HY, Stuetzle W: Wavelet shrinkage for unequally spaced data. Statistics and Computing 1999, 9: 65–75. 10.1023/A:1008818328241View ArticleGoogle Scholar
- Alter O, Brown PO, Botstein D: Singular value decomposition for genome-wide expression data processing and modeling. Proc Natl Acad Sci U S A 2000, 97: 10101–10106. 10.1073/pnas.97.18.10101PubMed CentralView ArticlePubMedGoogle Scholar
- Hou ZJ: Adaptive singular value decomposition in wavelet domain for image denoising. Pattern Recognition 2003, 36: 1747–1763. 10.1016/S0031-3203(02)00323-0View ArticleGoogle Scholar
- Donoho DL: De-noising by soft thresholding. IEEE Transactions on Information Theory 1995, 41: 613–627. 10.1109/18.382009View ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.