Comparative expressed sequence hybridization to chromosomes for tumor classification and identification of genomic regions of differential gene expression

Massive expansion of sex-specific SNPs, transposon-related elements, and neocentromere formation shape the young W-chromosome from the mosquitofish Gambusia affinis

BACKGROUND

The Western mosquitofish, Gambusia affinis, is a model for sex chromosome organization and evolution of female heterogamety. We previously identified a G. affinis female-specific marker, orthologous to the aminomethyl transferase (amt) gene of the related platyfish (Xiphophorus maculatus). Here, we have analyzed the structure and differentiation of the G. affinis W-chromosome, using a cytogenomics and bioinformatics approach.

RESULTS

The long arm of the G. affinis W-chromosome (Wq) is highly enriched in dispersed repetitive sequences, but neither heterochromatic nor epigenetically silenced by hypermethylation. In line with this, Wq sequences are highly transcribed, including an active nucleolus organizing region (NOR). Female-specific SNPs and evolutionary young transposable elements were highly enriched and dispersed along the W-chromosome long arm, suggesting constrained recombination. Wq copy number expanded elements also include female-specific transcribed sequences from the amt locus with homology to TE. Collectively, the G. affinis W-chromosome is actively differentiating by sex-specific copy number expansion of transcribed TE-related elements, but not (yet) by extensive sequence divergence or gene decay.

CONCLUSIONS

The G. affinis W-chromosome exhibits characteristic genomic properties of an evolutionary young sex chromosome. Strikingly, the observed sex-specific changes in the genomic landscape are confined to the W long arm, which is separated from the rest of the W-chromosome by a neocentromere acquired during sex chromosome evolution and may thus have become functionally insulated. In contrast, W short arm sequences were apparently shielded from repeat-driven differentiation, retained Z-chromosome like genomic features, and may have preserved pseudo-autosomal properties.

Fluorescence In Situ Hybridization and Rehybridization Using Bacterial Artificial Chromosome Probes

Fluorescence in situ hybridization (FISH) method enables in situ genetic analysis of both metaphase and interphase cells from different types of material, including cell lines, cell smears, and fresh and paraffin-embedded tissue. Despite the growing number of commercially available FISH probes, still for large number of gene loci or chromosomal regions commercial probes are not available. Here we describe a simple method for generating FISH probes using bacterial artificial chromosomes (BAC). Due to genome-wide coverage of BAC clones, there are almost unlimited possibilities for the analysis of any genomic regions using BAC FISH probes.

Random frog: an efficient reversible jump Markov Chain Monte Carlo-like approach for variable selection with applications to gene selection and disease classification

The identification of disease-relevant genes represents a challenge in microarray-based disease diagnosis where the sample size is often limited. Among established methods, reversible jump Markov Chain Monte Carlo (RJMCMC) methods have proven to be quite promising for variable selection. However, the design and application of an RJMCMC algorithm requires, for example, special criteria for prior distributions. Also, the simulation from joint posterior distributions of models is computationally extensive, and may even be mathematically intractable. These disadvantages may limit the applications of RJMCMC algorithms. Therefore, the development of algorithms that possess the advantages of RJMCMC methods and are also efficient and easy to follow for selecting disease-associated genes is required. Here we report a RJMCMC-like method, called random frog that possesses the advantages of RJMCMC methods and is much easier to implement. Using the colon and the estrogen gene expression datasets, we show that random frog is effective in identifying discriminating genes. The top 2 ranked genes for colon and estrogen are Z50753, U00968, and Y10871_at, Z22536_at, respectively. (The source codes with GNU General Public License Version 2.0 are freely available to non-commercial users at: http://code.google.com/p/randomfrog/.).

Chemical proteomics identifies Nampt as the target of CB30865, an orphan cytotoxic compound

Drug discovery based on cellular phenotypes is impeded by the challenge of identifying the molecular target. To alleviate this problem, we developed a chemical proteomic process to identify cellular proteins that bind to small molecules. CB30865 is a potent (subnanomolar) and selective cytotoxic compound of previously unknown mechanism of action. By combining chemical proteomics with biochemical and cellular pharmacology we have determined that CB30865 cytotoxicity is due to subnanomolar inhibition of nicotinamide phosphoribosyltransferase (Nampt), an enzyme present in the NAD biosynthetic pathway. Cancer cells develop dependence on Nampt due to increased energy requirements and the elevated activity of NAD consuming enzymes such as sirtuins and mono and poly(ADP-ribose) polymerases (PARPs). These findings suggest new chemical starting points for Nampt inhibitors and further implicate this enzyme as a target in cancer.

The association of CCND1 overexpression and cisplatin resistance in testicular germ cell tumors and other cancers

Development of chemoresistance limits the clinical efficiency of platinum-based therapy. Although many resistance mechanisms have been demonstrated, genetic/molecular alterations responsible for drug resistance in the majority of clinical cases have not been identified. We analyzed three pairs of testicular germ cell tumor cell lines using Affymetrix expression microarrays and revealed a limited number of differentially expressed genes across the cell lines when comparing the parental and resistant cells. Among them, CCND1 was the most significantly differentially expressed gene. Analysis of testicular germ cell tumor clinical samples by quantitative reverse transcription PCR analysis revealed that overall expression of CCND1 was significantly higher in resistant cases compared with sensitive samples (P < 0.0001). We also found that CCND1 was dramatically overexpressed both in induced and intrinsically resistant samples of ovarian and prostate cancer. Finally combined CCND1 knockdown using small-interfering RNA and cisplatin treatment inhibited cell growth in vitro significantly more effectively than any of these single treatments. Therefore, deregulation of CCND1 may be a major cause of cisplatin resistance in testicular germ cell tumors and may also be implicated in ovarian and prostate cancers. CCND1 could be potentially used as a marker for treatment stratification and as a molecular target to improve the treatment of platinum-resistant tumors.

Comparative expressed sequence hybridization studies of t(11;18)(q21;q21)-positive and -negative gastric MALT lymphomas reveal both unique and overlapping gene programs

Among the genetic abnormalities reported to occur in MALT lymphomas, the translocation t(11;18)(q21;q21) is of particular interest because it is exclusively documented in MALT lymphomas, mainly with gastrointestinal location. It results in the creation of a fusion protein API2-MALT1 that activates the transcription factor NF-kappaB through enhanced IKK gamma polyubiquitination. Here, we apply the recently developed molecular technique termed comparative expressed sequence hybridization to identify differentially expressed chromosomal regions related to the pathogenesis of gastric MALT lymphomas. By comparing t(11;18)(q21;q21)-positive gastric MALT lymphomas to their t(11;18)(q21;q21)-negative counterparts, we found that the location of the MALT1 break point determines a difference in expression pattern within the t(11;18)(q21;q21)-positive group. Moreover, we could define a gastric MALT lymphoma signature, which most likely comprises the regions and genes with significance in the development of MALT lymphomas, by comparing both t(11;18)(q21;q21)-positive and -negative MALT lymphomas to normal lymphoid tissue. Finally, a significant imprint of the marginal zone signature, established by comparing microdissected, splenic B follicles with and without marginal zone, was evident in the expression profile of MALT lymphoma, further supporting a marginal zone origin for this type of B-cell non-Hodgkin's lymphoma.

Is there a difference in childhood T-cell acute lymphoblastic leukaemia and T-cell lymphoblastic lymphoma?

To distinguish the similarities or differences between T-cell acute lymphoblastic leukaemia (T-ALL) and T-cell lymphoblastic lymphoma (T-LBL), we retrospectively analyzed the clinical, immunophenotypic, cytogenetic, and molecular characteristics in 37 children diagnosed between December 1990 and December 2003. Comparative Expressed Sequence Hybridisation (CESH) was used to determine gene expressing profile in both diseases. Twenty two patients suffered from T-ALL and 15 patients were diagnosed as T-LBL. Immunophenotyping demonstrated a more immature phenotype in T-ALL and a more mature phenotype in T-LBL. Cytogenetic and molecular genetic aberrations were found in 82% of T-ALL compared with 73% of T-LBL. By CESH gene expression profiling, the investigated cases were segregated into two groups that largely corresponded with T-ALL and T-LBL. The clinical presentation and cytogenetic characteristics are largely similar for T-ALL and T-LBL supporting the concept that both represent a spectrum of one single disease. The differences that were found between both neoplasms, in particular in their phenotype and in their expression profile may suggest that most T-ALL derive from a T-cell progenitor of the bone marrow, while thymocytes represent the normal counterpart of T-LBL.

Translocations targeting CCND2, CCND3, and MYCN do occur in t(11;14)-negative mantle cell lymphomas

The genetics of t(11;14)(q13;q32)/cyclin D1-negative mantle cell lymphoma (MCL) is poorly understood. We report here 8 MCL cases lacking t(11;14) or variant CCND1 rearrangement that showed expression of cyclin D1 (2 cases), D2 (2 cases), and D3 (3 cases). One case was cyclin D negative. Cytogenetics and fluorescence in situ hybridization detected t(2;12)(p11;p13)/IGK-CCND2 in one of the cyclin D2-positive cases and t(6;14)(p21;q32)/IGH-CCND3 in one of the cyclin D3-positive cases. Moreover, we identified a novel cryptic t(2;14)(p24;q32) targeting MYCN in 2 blastoid MCLs: one negative for cyclin D and one expressing cyclin D3. Interestingly, both cases showed expression of cyclin E. Notably, all 3 blastoid MCLs showed a monoallelic deletion of RB1 associated with a lack of expression of RB1 protein and monoallelic loss of p16. In sum-mary, this study confirms frequent aberrant expression of cyclin D2 and D3 in t(11;14)-negative MCLs and shows a t(11;14)-independent expression of cy-clin D1 in 25% of present cases. Novel findings include cyclin E expression in 2 t(11;14)-negative MCLs characterized by a cryptic t(2;14)(p24;q32) and identification of MYCN as a new lymphoma oncogene associated with a blastoid MCL. Clinically important is a predisposition of t(11;14)-negative MCLs to the central nervous system involvement.

Genes proximal and distal to MYCN are highly expressed in human neuroblastoma as visualized by comparative expressed sequence hybridization

MYCN amplification is associated with poor prognosis in neuroblastoma disease. To improve our understanding of the influence of the MYCN amplicon and its corresponding expression, we investigated the 2p expression pattern of MYCN amplified (n = 13) and nonamplified (n = 4) cell lines and corresponding primary tumors (n = 3) using the comparative expressed sequence hybridization technique. All but one MYCN amplified cell line displayed overexpression at 2p. Expression peaks were observed frequently at 2pter and less frequently at 2p24 (MYCN locus), 2p23.3-23.2, and/or 2p23.1. Importantly, cell lines and two corresponding primary tumors displayed expression peaks at similar loci. No significant 2p24 expression level was observed for those cell lines displaying a low amplification rate (n = 3) by comparative genomic hybridization. Only the cell lines with an enhanced peak at 2p23.2-23.3 displayed coamplification of the ALK gene (2p23.2), reported to be associated with unfavorable prognosis. Finally, two of four cell lines without MYCN amplification, both derived from patients with poor outcome, also showed an expression peak at 2p23.2. These data indicate that, besides MYCN, other genes proximal and distal to MYCN are highly expressed in neuroblastoma. The prognostic significance of expression peaks at 2p23.2-23.3, independent of MYCN and ALK status, remains to be investigated.

Genomic copy number and expression patterns in testicular germ cell tumours

Testicular germ cell tumours of adults and adolescents (TGCT) include seminomas (SE) and nonseminomas (NS), with spermatocytic seminomas (SSE) representing a distinct entity in older men. SE and NS have gain of 12p material in all cases, whereas SSE are associated with overrepresentation of chromosome 9. Here, we compare at the chromosomal level, copy number imbalances with global expression changes, identified by comparative expressed sequence hybridisation analyses, in seven SE, one combined tumour, seven NS and seven cell lines. Positive correlations were found consistent with copy number as a main driver of expression change, despite reported differences in methylation status in SE and NS. Analysis of chromosomal copy number and expression data could not distinguish between SE and NS, in-keeping with a similar genetic pathogenesis. However, increased expression from 4q22, 5q23.2 and 9p21 distinguished SSE from SE and NS and decreased copy number and expression from 2q36–q37 and 6q24 was a specific feature of NS-derived cell lines. Our analysis also highlights 19 regions with both copy number and expression imbalances in greater than 40% of cases. Mining available expression array data identified genes from these regions as candidates for involvement in TGCT development. Supplementary data is available at http://www.crukdmf.icr.ac.uk/array/array.html.

Relationship between classic Hodgkin lymphoma and overlapping large cell lymphoma investigated by comparative expressed sequence hybridization expression profiling

There is a diagnostic grey zone between classic Hodgkin lymphoma (cHL) and some non-Hodgkin lymphoma (NHL), including primary mediastinal B cell lymphoma, diffuse large B cell lymphoma, and anaplastic large cell lymphoma. They all have some morphological and/or phenotypic features in common. To investigate this, we undertook an expression profiling study of these lymphomas using comparative expressed sequence hybridization. This technique detects chromosomal regions that are differentially expressed between a test and a reference tissue in a manner similar to comparative genomic hybridization, and is particularly suitable when the number of informative biopsies is limited. Using this approach, we identified a unique expression profile for all lymphoma types investigated. Unsupervised hierarchical cluster analysis of the acquired data showed that cHL separates from all investigated NHLs, including ALCL-like HL. Moreover, anaplastic lymphoma kinase (ALK)-negative ALCL clustered in a separate branch together with ALCL-like HL. Thus, analysing the neoplastic cells concurrently with their microenvironment, ALK-negative ALCL and ALCL-like HL seem to be related to each other, while cHL constitutes a separate lymphoma entity.

Optimization of comparative expressed sequence hybridization for genome-wide expression profiling at chromosome level

Comparative expressed sequence hybridization (CESH) has recently been developed for global expression profiling at chromosome level. To improve its specificity and sensitivity, we examined the effects of cDNA amplification and labeling methods on CESH profiles, using a gastric cancer cell line, Kato III, and compared the CESH profiles to cDNA microarray and reverse transcriptase-polymerase chain reaction (RT-PCR) data. CESH results were scarcely affected by the amplification process, either at the RNA level with T7 polymerase or at the cDNA level with degenerate oligonucleotide-primed PCR (DOP-PCR). The labeling method, however, did remarkably affect the CESH results; false positive shifts of the test/reference ratio (T/R) were not detected in self-matched CESH with pre-cDNA labeling and random priming labeling of cDNA but were consistently seen with DOP-PCR labeling in 11 chromosomes. The use of cDNA deriving from mRNA either with pre-cDNA or random priming labeling gave results of higher detection sensitivity for regions of up- or downregulated expression and higher concordance with the microarray and RT-PCR data in the corresponding regions than with conventional CESH. This modification of CESH with random priming labeling was found feasible by its application to Kato III cells with and without 5-aza-2'-deoxycytidine treatment; the regions identified as epigenetically silenced included genes that were reportedly methylated in Kato III.

Role for amplification and expression of glypican-5 in rhabdomyosarcoma

Overexpression of genes, through genomic amplification and other mechanisms, can critically affect the behavior of tumor cells. Genomic amplification of the 13q31-32 region is reported in many tumors, including rhabdomyosarcomas that are primarily pediatric sarcomas resembling developing skeletal muscle. The minimum overlapping region of amplification at 13q31-32 in rhabdomyosarcomas was defined as containing two genes: Glypican-5 (GPC5) encoding a cell surface proteoglycan and C13orf25 encompassing the miR-17-92 micro-RNA cluster. Genomic copy number and gene expression analyses of rhabdomyosarcomas indicated that GPC5 was the only gene consistently expressed and up-regulated in all cases with amplification. Constitutive overexpression and knockdown of GPC5 expression in rhabdomyosarcoma cell lines increased and decreased cell proliferation, respectively. A correlation between expression levels of nascent pre-rRNA and GPC5 (P = 0.001), but not a C13orf25 transcript containing miR-17-92, in primary samples supports an association of GPC5 with proliferative capacity in vivo. We show that GPC5 increases proliferation through potentiating the action of the growth factors fibroblast growth factor 2 (FGF2), hepatocyte growth factor (HGF), and Wnt1A. GPC5 enhanced the intracellular signaling of FGF2 and HGF and altered the cellular distribution of FGF2. The mesoderm-inducing effect of FGF2 and FGF4 in Xenopus blastocysts was also enhanced. Our data are consistent with a role of GPC5, in the context of sarcomagenesis, in enhancing FGF signaling that leads to mesodermal cell proliferation without induction of myogenic differentiation. Furthermore, the properties of GPC5 make it an attractive target for therapeutic intervention in rhabdomyosarcomas and other tumors that amplify and/or overexpress the gene.

Nascent pre-rRNA overexpression correlates with an adverse prognosis in alveolar rhabdomyosarcoma

Comparative expressed sequence hybridization (CESH) is an expression profiling technique which identifies chromosomal regions corresponding to differential gene expression. Here, we observe that various tumor samples including rhabdomyosarcoma show very prominent staining on the short arms of the acrocentric chromosomes suggesting an increase in expression of ribosomal RNA synthesized from the repetitive rDNA of the nucleolar organizer regions located on these chromosomes. Survival analysis showed a correlation with overexpression from this region and a poor prognosis in rhabdomyosarcoma. This phenomenon was studied in an extended set of rhabdomyosarcoma tumor samples using quantitative real-time reverse transcriptase-PCR to quantify levels of pre-rRNA (precursor ribosomal RNA). It was demonstrated first that the strong CESH signals did correspond to a marked increase in pre-rRNA expression and second that high pre-rRNA expression correlated with an adverse prognosis in alveolar subtype rhabdomyosarcoma. In addition, we demonstrate that pre-rRNA expression is significantly correlated with tumor stage. We conclude that measuring expression of pre-rRNA by real-time PCR is a useful prognostic marker in alveolar rhabdomyosarcoma. Furthermore, given that we have observed similar rDNA staining in all cancer types that we have studied by CESH, we propose that pre-rRNA overexpression is a general phenomenon in cancer and that our real-time PCR assay may be applicable as a prognostic marker in other tumor types.

A locally adaptive statistical procedure (LAP) to identify differentially expressed chromosomal regions

MOTIVATION

The systematic integration of expression profiles and other types of gene information, such as chromosomal localization, ontological annotations and sequence characteristics, still represents a challenge in the gene expression arena. In particular, the analysis of transcriptional data in context of the physical location of genes in a genome appears promising in detecting chromosomal regions with transcriptional imbalances often characterizing cancer.

RESULTS

A computational tool named locally adaptive statistical procedure (LAP), which incorporates transcriptional data and structural information for the identification of differentially expressed chromosomal regions, is described. LAP accounts for variations in the distance between genes and in gene density by smoothing standard statistics on gene position before testing the significance of their differential levels of gene expression. The procedure smooths parameters and computes p-values locally to account for the complex structure of the genome and to more precisely estimate the differential expression of chromosomal regions. The application of LAP to three independent sets of raw expression data allowed identifying differentially expressed regions that are directly involved in known chromosomal aberrations characteristic of tumors.

AVAILABILITY

Functions in R for implementing the LAP method are available at http://www.dpci.unipd.it/Bioeng/Publications/LAP.htm

Molecular characterisation of the t(1;15)(p22;q22) translocation in the prostate cancer cell line LNCaP

Although chromosome translocations are well-documented recurrent events in hematological malignancies and soft tissue sarcomas, their significance in carcinomas is less clear. We report here the molecular characterization of the reciprocal translocation t(1;15)(p22;q22) in the prostate carcinoma cell line, LNCaP. The chromosome 1 breakpoint was localized to a single BAC clone, RP11-290M5, by sequential FISH analysis of clones selected from the NCBI chromosome 1 map. This was further refined to a 580-bp region by Southern blot analysis. A 2.85-kb fragment spanning the der(1) breakpoint was amplified by long-range inverse PCR. The breakpoint on chromosome 1 was shown to lie between the CYR61 and the DDAH1 genes with the der(1) junctional sequence linking the CYR61 gene to the TSPAN3 (TM4SF8) gene on chromosome 15. Confirmatory PCR and FISH mapping of the der(15) showed loss of chromosome material proximal to the breakpoint on chromosome 15, containing the PSTPIP1 and RCN2 genes. On the available evidence we conclude that this translocation does not result in an in-frame gene fusion. Comparative expressed sequence hybridization (CESH) and comparative genomic hybridization (CGH) analysis, showed relative down-regulation of gene expression surrounding the breakpoint, but no gross change in genomic copy number. Real-time quantitative RT-PCR for genes around the breakpoint supported the CESH data. Therefore, here we may have revealed a gene down-regulation mechanism associated with a chromosome translocation, either through small deletion at the breakpoint or through another means of chromosome domain related gene regulation.

Genomic and Expression Profiling of Human Spermatocytic Seminomas: Primary Spermatocyte as Tumorigenic Precursor and DMRT1 as Candidate Chromosome 9 Gene

Spermatocytic seminomas are solid tumors found solely in the testis of predominantly elderly individuals. We investigated these tumors using a genome-wide analysis for structural and numerical chromosomal changes through conventional karyotyping, spectral karyotyping, and array comparative genomic hybridization using a 32 K genomic tiling-path resolution BAC platform (confirmed by in situ hybridization). Our panel of five spermatocytic seminomas showed a specific pattern of chromosomal imbalances, mainly numerical in nature (range, 3-24 per tumor). Gain of chromosome 9 was the only consistent anomaly, which in one case also involved amplification of the 9p21.3-pter region. Parallel chromosome level expression profiling as well as microarray expression analyses (Affymetrix U133 plus 2.0) was also done. Unsupervised cluster analysis showed that a profile containing transcriptional data on 373 genes (difference of > or = 3.0-fold) is suitable for distinguishing these tumors from seminomas/dysgerminomas. The diagnostic markers SSX2-4 and POU5F1 (OCT3/OCT4), previously identified by us, were among the top discriminatory genes, thereby validating the experimental set-up. In addition, novel discriminatory markers suitable for diagnostic purposes were identified, including Deleted in Azospermia (DAZ). Although the seminomas/dysgerminomas were characterized by expression of stem cell-specific genes (e.g., POU5F1, PROM1/CD133, and ZFP42), spermatocytic seminomas expressed multiple cancer testis antigens, including TSP50 and CTCFL (BORIS), as well as genes known to be expressed specifically during prophase meiosis I (TCFL5, CLGN, and LDHc). This is consistent with different cells of origin, the primordial germ cell and primary spermatocyte, respectively. Based on the region of amplification defined on 9p and the associated expression plus confirmatory immunohistochemistry, DMRT1 (a male-specific transcriptional regulator) was identified as a likely candidate gene for involvement in the development of spermatocytic seminomas.

Comparative expressed sequence hybridization reveals differential gene expression in morphological breast cancer subtypes

In this study, comparative expressed sequence hybridization (CESH) has been used to compare gene expression patterns in three morphologically different breast cancer subtypes: classic-type invasive lobular carcinoma (ILC), poorly differentiated ERBB2-negative invasive ductal carcinoma-not otherwise specified (IDC-NOS), and poorly differentiated ERBB2-positive IDC-NOS. CESH allows global detection of chromosomal regions with differential gene expression in a way similar to that of comparative genomic hybridization (CGH). Eight cases of each breast cancer subtype were included in the study. For each subtype, two pools of four cases each were constructed. CESH was used to compare both pools within the same morphological subtype, followed by a comparison of pools belonging to different subtypes. This revealed three chromosomal regions that were differentially expressed in ductal and lobular carcinomas, including relative overexpression at 8q13-q23 and 16q22, and relative underexpression at 8p21-p22. In addition, an expression signature characterized by relative overexpression at 3q24-q26.3, 14q23-31, 17q12, and 20q12-13 was identified for ERBB2-positive IDC-NOS. In summary, CESH analysis highlights chromosomal regions of differential gene expression that are associated with morphologically defined breast cancer subtypes and suggests that regions on chromosome 8 are of interest in the discrimination between ductal and lobular carcinomas. In addition, using CESH, it was possible to identify an ERBB2 expression signature, comprising four chromosomal regions with potential significance in the aggressive behaviour of ERBB2-positive IDC-NOS.

CGH-Profiler: Data mining based on genomic aberration profiles

Background

CGH-Profiler is a program that supports the analysis of genomic aberrations measured by Comparative Genomic Hybridisation (CGH). Comparative genomic hybridisation (CGH) is a well-established, molecular cytogenetic method that allows the detection of chromosomal imbalances in entire genomes. This technique is widely used in routine molecular diagnostics. Typically, chromosomal imbalances are described in a complex syntax based on the International Standard for Cytogenetic Nomenclature (ISCN). This semantic description of chromosomal imbalances hinders a large-scale statistical analysis across different experiments, e.g. for finding aberration patterns associated with a particular disease type or state.

Results

CGH-Profiler circumvents the semantic ISCN description by importing data from different CGH system vendors and by directly transferring the data into a table format that is readily accessible for subsequent statistical analysis. CGH-profiler comes with different consistency checks, calculates various statistics and automatically assigns a median copy number ratio to each chromosomal band.

Import of CGH profiles from different CGH system vendors is already supported; its extension to other systems can be readily achieved through Perl scripts.

CGH profiler can also be used to analyse comparative expressed sequence hybridisation (CESH) data. CESH reveals gene expression patterns according to chromosomal locations in a similar manner as CGH detects chromosomal imbalances.

Conclusion

CGH-Profiler is a useful tool for processing of CGH and CESH data.

Identification of candidate genes associated with salivary adenoid cystic carcinomas using combined comparative genomic hybridization and oligonucleotide microarray analyses

Adenoid cystic carcinoma (ACC) of the salivary gland often has a variable clinical course with a poor prognosis. To investigate DNA copy number aberrations associated with ACCs, we compared comparative genome hybridization data from ACCs (n = 6) with other types of salivary gland tumors such as adenocarcinomas (n = 3) and pleomorphic adenomas (n = 6). While 15 gain loci (1q32, 6p25, 6q21-q24, 7q11.2, 7q31, 10q11.2, 11p12-q12, 12q13, 12q14, 13q24, 16p13.3-13.2, 18p11.3, 18q23, 19q13.4, and Xq28) were detected, no DNA loss locus was evident. To examine the expression status of genes on the ACC-associated loci, transcriptional measurements of approximately 38000 human genes then were monitored using Affymetrix U133 Plus 2.0 GeneChips. A total of 4431 genes were found differentially expressed by at least two-fold between ACCs and normal salivary glands. Of them, 3162 genes were up-regulated and 1269 genes were down-regulated in ACCs. After obtaining locus information about the RNA transcripts from the Affymetrix database, we found 262 ACC-associated genes with increased expression on ACC-associated loci. To investigate functional network and gene ontology, the 262 genes were analyzed using Ingenuity Pathway Analysis Tool. The function with the highest P value was a cancer-related function (P = 2.52e-4 to 4.71e-2). In addition, we identified pituitary tumor-transforming gene 1 and transformation related protein 63 genes that were up-regulated by increasing DNA copy number and modulated expression of oncogenes. These results suggested that the combination of copy number and gene expression profiling provides an improved strategy for gene identification in salivary gland ACCs.

Diminished transcription of chromosome arm 4q is inversely related to local spreading of hepatocellular carcinoma

To identify transcriptional features associated with local spreading of hepatocellular carcinoma (HCC), regional expression biases were determined from the gene expression profiles of solitary HCC tumors and of HCC tumors associated with additional intrahepatic lesions and/or portal vein involvement. Regional expression biases are genetic intervals in which gene expression is coordinately changed. Often, but not always, a regional expression bias indicates an underlying chromosomal aberration. This study demonstrated that the presence or absence of a regional expression bias on chromosome arm 4q can predict the presence or absence of intrahepatic spread/portal vein involvement in 83% of cases analyzed (n = 40). These results suggest that the transcriptional state of 4q may be a marker for local spreading of HCC and that a more detailed genetic/epigenetic characterization of this locus may provide additional insights into HCC development.

Molecular cytogenetics in haematological malignancy: current technology and future prospects

Cytogenetics has played a pivotal role in haematological malignancy, both as an aid to diagnosis and in identifying recurrent chromosomal rearrangements, an essential prerequisite to identifying genes involved in leukaemia and lymphoma pathogenesis. In the late 1980s, a series of technologies based around fluorescence in situ hybridisation (FISH) revolutionised the field. Interphase FISH, multiplex-FISH (M-FISH, SKY) and comparative genomic hybridisation (CGH) have emerged as the most significant of these. More recently, microarray technologies have come to prominence. In the acute leukaemias, the finding of characteristic gene expression signatures corresponding to biological subgroups has heralded gene expression profiling as a possible future alternative to current cytogenetic and morphological methods for diagnosis. In the lymphomas, high-resolution array CGH has successfully identified new regions of deletion and amplification, providing the prospect of disease-specific arrays.

Processing of radical prostatectomy specimens for correlation of data from histopathological, molecular biological, and radiological studies: a new whole organ technique

AIMS

To develop a method of processing non-formalin fixed prostate specimens removed at radical prostatectomy to obtain fresh tissue for research and for correlating diagnostic and molecular results with preoperative imaging.

METHODS/RESULTS

The method involves a prostate slicing apparatus comprising a tissue slicer with a series of juxtaposed planar stainless steel blades linked to a support, and a cradle adapted to grip the tissue sample and receive the blades. The fresh prostate gland is held in the cradle and the blades are moved through the cradle slits to produce multiple 4 mm slices of the gland in a plane perpendicular to its posterior surface. One of the resulting slices is preserved in RNAlater. The areas comprising tumour and normal glands within this preserved slice can be identified by matching it to the haematoxylin and eosin stained sections of the adjacent slices that are formalin fixed and paraffin wax embedded. Intact RNA can be extracted from the identified tumour and normal glands within the RNAlater preserved slice. Preoperative imaging studies are acquired with the angulation of axial images chosen to be similar to the slicing axis, such that stained sections from the formalin fixed, paraffin wax embedded slices match their counterparts on imaging.

CONCLUSIONS

A novel method of sampling fresh prostate removed at radical prostatectomy that allows tissue samples to be used both for diagnosis and molecular analysis is described. This method also allows the integration of preoperative imaging data with histopathological and molecular data obtained from the prostate tissue slices.

Comparison of array-based comparative genomic hybridization with gene expression-based regional expression biases to identify genetic abnormalities in hepatocellular carcinoma

BACKGROUND

Regional expression biases (REBs) are genetic intervals where gene expression is coordinately changed. For example, if a region of the genome is amplified, often the majority of genes that map within the amplified region show increased expression when compared to genes located in cytogenetically normal regions. As such, REBs have the potential to act as surrogates for cytogenetic data traditionally obtained using molecular technologies such as comparative genomic hybridization. However as REBs are identified using transcriptional information, detection of REBs may also identify local transcriptional abnormalities produced by both genetic and epigenetic mechanisms.

RESULTS

REBs were identified from a set of hepatocellular carcinoma (HCC) gene expression profiles using a multiple span moving binomial test and compared to genetic abnormalities identified using array-based comparative genomic hybridization (aCGH). In the majority of cases, REBs overlapped genetic abnormalities as determined by aCGH. For example, both methods identified narrow regions of frequent amplification on chromosome 1p and narrow regions of frequent deletion on 17q. In a minority of cases, REBs were identified in regions not determined to be abnormal via other cytogenetic technologies. Specifically, expression biases reflective of cell proliferation were frequently identified on chromosome 6p21-23.

CONCLUSION

Identification of REBs using a multiple span moving binomial test produced reasonable approximations of underlying cytogenetic abnormalities. However, caution should be used when attributing REBs identified on chromosome 6p to cytogenetic events in rapidly proliferating cells.

Relationship between MYCN copy number and expression in rhabdomyosarcomas and correlation with adverse prognosis in the alveolar subtype

PURPOSE

Amplification of the transcription factor MYCN is an important molecular diagnostic tool in stratifying treatment for neuroblastoma. Increased copy number and overexpression of MYCN in the pediatric cancer rhabdomyosarcoma has been described in a number of small studies with conflicting conclusions about its association with clinicopathologic characteristics. We aimed to study the phenomenon in the largest series to date.

PATIENTS AND METHODS

Using quantitative polymerase chain reaction, we measured MYCN copy number and expression levels in rhabdomyosarcoma samples from 113 and 92 individuals with a confirmed diagnosis of rhabdomyosarcoma, respectively.

RESULTS

Increased copy number of MYCN was found to be a feature of both the embryonal and alveolar subtypes. The copy number and expression levels were significantly greater in the alveolar subtype, although the range of expression in both subtypes spanned several orders of magnitude. MYCN copy number showed a significant correlation with expression in the alveolar subtype; this relationship between copy number and expression could be modeled as a logarithmic function. It is notable that relatively high expression frequently occurred in embryonal rhabdomyosarcoma without high copy number and that low expression was found in some cases with high copy number. In patients with alveolar rhabdomyosarcoma, overexpression (greater than median) or gain of genomic copies of MYCN were significantly associated with adverse outcome.

CONCLUSION

MYCN deregulation is a feature of rhabdomyosarcoma tumorigenesis, defines groups of patients with a poor prognosis, and is a potential target for novel therapies.

Conventional and array-based comparative genomic hybridization analysis of nasopharyngeal carcinomas from the Mediterranean area

Nasopharyngeal carcinoma (NPC) occurs with a high incidence in Southeast Asia and to a lesser extent in the Mediterranean area, especially in Tunisia, Algeria, and Morocco. Cellular gene alterations combined with latent Epstein-Barr virus infection are thought to be essential for NPC oncogenesis. To date, chromosome analysis with comparative genomic hybridization (CGH) has been reported exclusively for NPCs from Southeast Asia. Although NPCs from the Mediterranean area have several distinct clinical and epidemiological features, CGH investigations have been lacking. Chromosome analysis was therefore undertaken on a series of NPC xenografts and biopsies derived from patients of Mediterranean origin. Four xenografts were investigated with a combination of conventional CGH, array-based CGH, and comparative expressed sequence hybridization. In addition, 23 fresh NPC biopsies were analyzed with conventional CGH. Data obtained from xenografts and fresh biopsies were consistent, except that amplification of genes at 18p was observed only in xenografts derived from metastatic tissues. Frequent gains associated with gene overexpression were detected at 1q25 approximately qter (64%) and 12p13 (50%). Losses were noticed mainly at 11q14 approximately q23 (50%), 13q12 approximately q31 (50%), 14q24 approximately q31 (43%), and 3p13 approximately p23 (43%). Comparison with previous reports suggests that Mediterranean NPCs have higher frequencies of gains at 1q and losses at 13q than their Asian counterparts.

Novel genes, possibly relevant for molecular diagnosis or therapy of human rhabdomyosarcoma, detected by genomic expression profiling

Transcriptional profiles of an alveolar rhabdomyosarcoma (RMS) and of a RMS cell line were reconstructed by a computational and statistical approach. Expression data of 29,963 genes in 11 adult human healthy tissues and in 37 tumour tissues were analysed for comparison. We identified 202 genes differentially expressed in at least one RMS sample, as compared with normal skeletal muscle. Among them, 107 resulted specifically overexpressed in RMS, but in no tumour affecting other tissues. Cluster analysis applied to expression data detected a series of genes presumably co-expressed with genes encoding known tumour markers and/or reportedly involved in genesis or development of rhabdomyosarcoma. This study succeeded in identifying a number of genes, which become candidates for in vitro study, thus facilitating discovery of novel tumour markers or targets for drug therapy.

Overexpression of genes on 16q associated with cisplatin resistance of testicular germ cell tumor cell lines

Testicular germ-cell tumors (TGCTs) show exquisite sensitivity to cisplatin-based chemotherapy, and therefore this is considered a good model system for studying the mechanism of chemotherapy resistance. Although the genetic alterations related to TGCT have been well studied, little is known about the genetic basis of chemotherapy resistance, which occurs in a small proportion of TGCTs. In this study, we investigated genomic and expression differences between three cisplatin-sensitive and their paired cisplatin-resistant lines using combined whole-genome screen approaches. Comparative genomic hybridization (CGH) analysis on chromosomes revealed genetic differences between the resistant and parent cell lines in each pair, but did not show any consistent chromosome changes in all three lines. Microarray CGH analysis generated some additional information of DNA copy number gains and losses including some important oncogenes, tumor-suppressor genes, and drug-resistance-related genes. However, no consistent genomic region changes were found in the three cell lines. Interestingly, when comparative expressed sequence hybridization, a technique for gene expression profiling along chromosomes, was applied, we discovered a consistently overexpressed chromosomal region in all three resistant lines compared with their parent lines. The minimum overlapping chromosomal region is at 16q22-23. Further definition of genes in this chromosomal region will aid our understanding of the mechanism of cisplatin resistance and may offer novel therapeutic targets.

Comparative expressed sequence hybridization studies of high-hyperdiploid childhood acute lymphoblastic leukemia

The functional consequences of a high-hyperdiploid karyotype, found in up to one-third of cases of acute lymphoblastic leukemia (ALL), are unknown. Using the technique of comparative expressed sequence hybridization (CESH), we sought to address the question of whether increased chromosome copies in hyperdiploid ALL lead to increased gene expression. Relative expression of hyperdiploid ALL blasts versus peripheral blood mononuclear cells was analyzed in 18 patients. Common regions of overexpression corresponding to the presence of tri-/tetrasomies included: Xp22.1-22.2, 4q28, 6q14-15, 6q24, 10p13, 14q23-24, 17q21, 18q12, and 21q21, identified in 28-89% of cases. However, increased expression without underlying trisomy occurred at 3p21.3, 7q11.2, 8p21, and 8q24.1 in 39-90% of cases. High expression at 7q11.2, the most consistent change detected, was confirmed by quantitative PCR. Poor correlation between the presence of tri-/tetrasomy and overexpression was observed for chromosomes 14 and 17. Two cases were reanalyzed versus (i) B cells, (ii) transformed B cells, and (iii) CD34+19+ cells (the putative counterpart of the leukemic cell). A reduction in the number of relatively overexpressed regions was observed with CD34+19+ cells. In particular, the peak at 7q11.2 disappeared, suggesting up-regulation of genes from this region in the early ontology of normal B-cell development. In conclusion, we have shown that tri-/tetrasomies in hyperdiploid ALL lead to an increase in the expression of associated sequences. The choice of a biologically relevant reference is crucial for data interpretation.

Detection of chromosomal regions showing differential gene expression in human skeletal muscle and in alveolar rhabdomyosarcoma

Background

Rhabdomyosarcoma is a relatively common tumour of the soft tissue, probably due to regulatory disruption of growth and differentiation of skeletal muscle stem cells. Identification of genes differentially expressed in normal skeletal muscle and in rhabdomyosarcoma may help in understanding mechanisms of tumour development, in discovering diagnostic and prognostic markers and in identifying novel targets for drug therapy.

Results

A Perl-code web client was developed to automatically obtain genome map positions of large sets of genes. The software, based on automatic search on Human Genome Browser by sequence alignment, only requires availability of a single transcribed sequence for each gene. In this way, we obtained tissue-specific chromosomal maps of genes expressed in rhabdomyosarcoma or skeletal muscle. Subsequently, Perl software was developed to calculate gene density along chromosomes, by using a sliding window. Thirty-three chromosomal regions harbouring genes mostly expressed in rhabdomyosarcoma were identified. Similarly, 48 chromosomal regions were detected including genes possibly related to function of differentiated skeletal muscle, but silenced in rhabdomyosarcoma.

Conclusion

In this study we developed a method and the associated software for the comparative analysis of genomic expression in tissues and we identified chromosomal segments showing differential gene expression in human skeletal muscle and in alveolar rhabdomyosarcoma, appearing as candidate regions for harbouring genes involved in origin of alveolar rhabdomyosarcoma representing possible targets for drug treatment and/or development of tumor markers.

Comparative expressed sequence hybridization studies of hairy cell leukemia show uniform expression profile and imprint of spleen signature

Comparative expressed sequence hybridization (CESH) to chromosomes is a recently introduced technique that identifies chromosomal regions corresponding to a differential gene expression. This technique is analogous to comparative genomic hybridization (CGH) that detects genomic imbalances. We applied CESH for the study of hairy cell leukemia (HCL), a disorder with a largely unknown expression profile. Twelve HCL cases with spleen involvement were investigated by CESH and CGH. While the latter analysis identified only a few nonrecurrent genomic imbalances, CESH showed a consistent expression profile in all HCL cases. In addition, pairing normal spleen with normal lymph node, a "spleen signature" was established by CESH. This signature most likely reflects the expression profile of spleen-specific components, such as the sinusoidal lining cells from the red pulp and the marginal zone B cells from the white pulp. Imprint of the spleen signature was found in the HCL expression profile, suggesting that HCL may originate from a particular B-cell subset present in these splenic components. Besides pairing HCL with normal lymph node and spleen, we identified an "HCL signature" comprising several chromosome regions with altered expression. The most significantly underexpressed regions include 3p24, 3p21, 3q13.3-q22, 4p16, 11q23, 14q22-q24, 15q21-q22, 15q24-q25, and 17q22-q24; and 13q31 and Xq13.3-q21 were the most significantly overexpressed. These regions possibly harbor genes related to the biology and the pathogenesis of HCL. Their identification warrants further molecular investigations.

Expression profiling targeting chromosomes for tumor classification and prediction of clinical behavior

Tumors are associated with altered or deregulated gene products that affect critical cellular functions. Here we assess the use of a global expression profiling technique that identifies chromosome regions corresponding to differential gene expression, termed comparative expressed sequence hybridization (CESH). CESH analysis was performed on a total of 104 tumors with a diagnosis of rhabdomyosarcoma, leiomyosarcoma, prostate cancer, and favorable-histology Wilms tumors. Through the use of the chromosome regions identified as variables, support vector machine analysis was applied to assess classification potential, and feature selection (recursive feature elimination) was used to identify the best discriminatory regions. We demonstrate that the CESH profiles have characteristic patterns in tumor groups and were also able to distinguish subgroups of rhabdomyosarcoma. The overall CESH profiles in favorable-histology Wilms tumors were found to correlate with subsequent clinical behavior. Classification by use of CESH profiles was shown to be similar in performance to previous microarray expression studies and highlighted regions for further investigation. We conclude that analysis of chromosomal expression profiles can group, subgroup, and even predict clinical behavior of tumors to a level of performance similar to that of microarray analysis. CESH is independent of selecting sequences for interrogation and is a simple, rapid, and widely accessible approach to identify clinically useful differential expression.

Expressive genomic hybridisation: gene expression profiling at the cytogenetic level

AIMS

To describe a cytogenetic technique suitable for the rapid assessment of global gene expression that is based on comparative genomic hybridisation (CGH), and to use it to understand the relation between genetic amplifications and gene expression.

METHODS

Whereas traditional CGH uses DNA as test and reference in hybridisations, expressive genomic hybridisation (EGH) uses globally amplified mRNA as test and normal DNA as reference. EGH is a rapid and powerful tool for localising and studying global gene expression profiles and correlating them with loci of genetic amplifications using traditional CGH.

RESULTS

EGH was used to correlate genetic amplifications detected by CGH with the expression profile of two independent cell lines-Colo320 and T47D. Although many amplifications resulted in overexpression, other amplifications were partially or completely silenced at the cytogenetic level.

CONCLUSION

This technique will assist in the analysis of overexpressed genes within amplicons and could resolve a controversial issue in cancer cytogenetics; namely, the relation between genetic amplifications and overexpression.

Expression profile of genes from 12p in testicular germ cell tumors of adolescents and adults associated with i(12p) and amplification at 12p11.2-p12.1

Gain of 12p material is invariably associated with testicular germ cell tumors (TGCTs) of adolescents and adults, most usually as an isochromosome 12p. We analyzed TGCTs with i(12p) using a global approach to expression profiling targeting chromosomes (comparative expressed sequence hybridization, CESH). This indicated overexpression of genes from 12p11.2-p12.1 relative to testis tissue and fibroblasts. The nonseminoma subtype showed higher levels of expression than seminomas. Notably, 12p11.2-p12.1 is amplified in about 10% of TGCTs and CESH analysis of such amplicon cases showed high levels of overexpression from this region. Microarray analysis, including cDNA clones representing most UniGene clusters from 12p11.2-p12.1, was applied to DNA and RNA from 5 TGCTs with amplification of 12p11.2-p12.1 and seven TGCTs with gain of the entire short arm of chromosome 12. Expression profiles were consistent with the CESH data and overexpression of EST595078, MRPS35 and LDHB at 12p11.2-p12.1 was detected in most TGCTs. High-level overexpression of BCAT1 was specific to nonseminomas and overexpression of genes such as CMAS, EKI1, KRAS2, SURB7 and various ESTs correlated with their amplification. Genes such as CCND2, GLU3, LRP6 and HPH1 at 12p13 were also overexpressed. The overexpressed sequences identified, particularly those in the region amplified, represent candidate genes for involvement in TGCT development.

Role of gain of 12p in germ cell tumour development

Within the human testis, three entities of germ cell tumours are distinguished: the teratomas and yolk sac tumors of newborn and infants, the seminomas and nonseminomas of adolescents and young adults, referred to as testicular germ cell tumours (TGCT), and the spermatocytic seminomas. Characteristic chromosomal anomalies have been reported for each group, supporting their distinct pathogenesis. TGCT are the most common cancer in young adult men. The initiating pathogenetic event of these tumours occurs during embryonal development, affecting a primordial germ cell or gonocyte. Despite this intra-uterine initiation, the tumour will only be clinically manifest after puberty, with carcinoma in situ (IS) as the precursor. All invasive TGCT, both seminomas and nonseminomas, as well as CIS cells are aneuploid. The only consistent (structural) chromosomal abnormalities in invasive TGCT are gains of the short arm of chromosome 12, mostly due to isochromosome (i(12p)) formation. This suggests that an increase in copy number of a gene(s) on 12p is associated with the development of a clinically manifest TGCT. Despite the numerous (positional) candidate gene approaches that have been undertaken thus far, identification of a causative gene(s) has been hampered by the fact that most 12p gains involve rather large genomic intervals, containing unmanageable numbers of candidate genes. Several years ago, we initiated a search for 12p candidate genes using TGCT with a restricted 12p-amplification, cytogenetically identified as 12p11.2-p12.1. This approach is mainly based on identification of candidate genes mapped within the shortest region of overlap of amplification (SROA). In this review, data will be presented, which support the model that gain of 12p-sequences is associated with suppression of apoptosis and Sertoli cell-independence of CIS cells. So far, DAD-R is one of the most likely candidate genes involved in this process, possibly via N-glycosylation. Preliminary results on high through-put DNA- and cDNA array analyses of 12p-sequences will be presented.

Chromosome and gene alterations in breast cancer as markers for diagnosis and prognosis as well as pathogenetic targets for therapy

Chromosomal abnormalities have been implicated in cancer development since the turn of the last century. Only during the past two decades, with advances in cytogenetics and molecular biology, has the genetic basis of neoplasia been firmly established, however, with chromosomal alterations being recognized as critical in the pathogenesis of human cancer. Recurrent chromosomal alterations provide cytological and molecular markers for the diagnosis and prognosis of disease. They also facilitate the identification of genes that are important in carcinogenesis and, ultimately, may lead to the development of targeted therapy. In breast cancer, the most prevalent malignancy among females, substantial progress has been achieved in identifying genes located at sites of recurrent chromosomal alterations and in profiling gene expression through the application of powerful cytogenetic and functional genomic techniques. Characterization of the molecular pathologic characteristics and gene-expression profiles of breast cancer should provide new clinical tools for the accurate diagnosis and prediction of prognosis as well as new targets for the development of therapeutic agents.

Epstein-barr virus-induced changes in B-lymphocyte gene expression

To elucidate the mechanisms by which Epstein-Barr virus (EBV) latency III gene expression transforms primary B lymphocytes to lymphoblastoid cell lines (LCLs), the associated alterations in cell gene expression were assessed by using 4,146 cellular cDNAs arrayed on nitrocellulose filters and real-time reverse transcription-PCR (RT-PCR). A total of 1,405 of the 4,146 cDNAs were detected using cDNA probes from poly(A)(+) RNA of IB4 LCLs, a non-EBV-infected Burkitt's lymphoma (BL) cell line, BL41, or EBV latency III-converted BL41 cells (BL41EBV). Thirty-eight RNAs were consistently twofold more abundant in the IB4 LCL and BL41EBV than in BL41 by microarray analysis. Ten of these are known to be EBV induced. A total of 23 of 28 newly identified EBV-induced genes were confirmed by real-time RT-PCR. In addition, nine newly identified genes and CD10 were EBV repressed. These EBV-regulated genes encode proteins involved in signal transduction, transcription, protein biosynthesis and degradation, and cell motility, shape, or adhesion. Seven of seven newly identified EBV-induced RNAs were more abundant in newly established LCLs than in resting B lymphocytes. Surveys of eight promoters of newly identified genes implicate NF-kappaB or PU.1 as potentially important mediators of EBV-induced effects through LMP1 or EBNA2, respectively. Thus, examination of the transcriptional effects of EBV infection can elucidate the molecular mechanisms by which EBV latency III alters B lymphocytes.

Chromosome 1q expression profiling and relapse in Wilms' tumour

Gene-expression patterns have proved useful for tumour classification and prediction of prognosis. We investigated a new gene-expression profiling method that targets individual chromosomes--comparative expressed sequence hybridisation (CESH)--to establish if any expression pattern was associated with relapse, and if this pattern could accurately predict outcome. We studied 18 cases of Wilms' tumour with favourable histological findings. Relative overexpression of genes on the long arm of chromosome 1 was shown in all tumours that subsequently relapsed but in none of those from patients in remission (odds ratio 357, 95% CI 6.4-19960). This finding shows the value of CESH analysis to rapidly screen for chromosomal regions that carry genes whose differential expression is important in establishment of tumour behaviour.