VAMP: visualization and analysis of array-CGH, transcriptome and other molecular profiles

Genome-wide analysis of cutaneous T-cell lymphomas identifies three clinically relevant classes

This study was undertaken to identify recurrent genetic alterations of the three main types of cutaneous T-cell lymphomas (CTCLs): mycosis fungoides (MF), Sézary syndrome (SS), and cutaneous anaplastic large-cell lymphoma (CALCL). Using array-based comparative genomic hybridization, the molecular cytogenetic profiles of 72 samples obtained from 58 patients with CTCL corresponding to 24 transformed MF (T-MF), 16 SS, and 18 CALCLs were determined. T-MF was characterized by gains of 1q25-31, 7p22-11.2, 7q21, 7q31, and 17q12, and losses of 9p21, 10p11.2, and 10q26. SS exhibited gains of 8q23-24.3 and 17q23-24, as well as losses of 9p21, 10p12-11.2, 10q22-24, 10q25-26, and 17p13-q11.1. Finally, CALCL exhibited 6q27 and 13q34 losses. Such imbalances were statistically associated with one CTCL subtype. Unsupervised hierarchical clustering defined three categories of clinical relevance: (1) CALCL apart from epidermotropic-CTCL, (2) an SS-only category, and (3) a mixed category with T-MF and SS cases, with both primary and secondary SS cases. In rare cases, the genetic classification did not correspond to the inclusion diagnosis, possibly reflecting the association of two diseases in the same patient or initial misdiagnosis according to follow-up. Finally, different samples in the same patient clustered together, showing reproducibility of such a classifier.

Extrachromosomal amplification mechanisms in a glioma with amplified sequences from multiple chromosome loci

Accumulation of extrachromosomal DNA molecules (double minute) is often responsible for gene amplification in cancers, but the mechanisms leading to their formation are still largely unknown. By using quantitative PCR, chromosome walking, in situ hybridization on metaphase chromosomes and whole genome analysis, we studied a glioma containing four extrachromosomally amplified loci (7p11, 1q32.1, 5p15 and 9p2). Complex extrachromosomal DNA molecules were formed by the fusion of several syntenic or non-syntenic DNA fragments from 7p11, 5p15 to 9p2. Fragments ranged from a few base pairs to megabase pairs. Scars of the amplification process remained at the original locus in the form of deletions or chromosome rearrangements. Chromosome fragmentation, due to replication stress, could explain this complex situation. In contrast, at 1q32.1, the initial extrachromosomal DNA molecule resulted from the circularization of a single fragment associated with an intrachromosomal deletion including, but larger than, the amplified sequence. The nature of the sequences involved in these rearrangements suggests that a V(D)J-like illegitimate recombination contributes to its formation.

Integrative oncogenomic analysis of microarray data in hematologic malignancies

During the last decade, gene expression microarrays and array-based comparative genomic hybridization (array-CGH) have unraveled the complexity of human tumor genomes more precisely and comprehensively than ever before. More recently, the simultaneous assessment of global changes in messenger RNA (mRNA) expression and in DNA copy number through "integrative oncogenomic" analyses has allowed researchers the access to results uncovered through the analysis of one-dimensional data sets, thus accelerating cancer gene discovery. In this chapter, we discuss the major contributions of DNA microarrays to the study of hematological malignancies, focusing on the integrative oncogenomic approaches that correlate genomic and transcriptomic data. We also present the basic aspects of these methodologies and their present and future application in clinical oncology.

CHESS (CgHExpreSS): a comprehensive analysis tool for the analysis of genomic alterations and their effects on the expression profile of the genome

Background

Genomic alterations frequently occur in many cancer patients and play important mechanistic roles in the pathogenesis of cancer. Furthermore, they can modify the expression level of genes due to altered copy number in the corresponding region of the chromosome. An accumulating body of evidence supports the possibility that strong genome-wide correlation exists between DNA content and gene expression. Therefore, more comprehensive analysis is needed to quantify the relationship between genomic alteration and gene expression. A well-designed bioinformatics tool is essential to perform this kind of integrative analysis. A few programs have already been introduced for integrative analysis. However, there are many limitations in their performance of comprehensive integrated analysis using published software because of limitations in implemented algorithms and visualization modules.

Results

To address this issue, we have implemented the Java-based program CHESS to allow integrative analysis of two experimental data sets: genomic alteration and genome-wide expression profile. CHESS is composed of a genomic alteration analysis module and an integrative analysis module. The genomic alteration analysis module detects genomic alteration by applying a threshold based method or SW-ARRAY algorithm and investigates whether the detected alteration is phenotype specific or not. On the other hand, the integrative analysis module measures the genomic alteration's influence on gene expression. It is divided into two separate parts. The first part calculates overall correlation between comparative genomic hybridization ratio and gene expression level by applying following three statistical methods: simple linear regression, Spearman rank correlation and Pearson's correlation. In the second part, CHESS detects the genes that are differentially expressed according to the genomic alteration pattern with three alternative statistical approaches: Student's t-test, Fisher's exact test and Chi square test. By successive operations of two modules, users can clarify how gene expression levels are affected by the phenotype specific genomic alterations. As CHESS was developed in both Java application and web environments, it can be run on a web browser or a local machine. It also supports all experimental platforms if a properly formatted text file is provided to include the chromosomal position of probes and their gene identifiers.

Conclusions

CHESS is a user-friendly tool for investigating disease specific genomic alterations and quantitative relationships between those genomic alterations and genome-wide gene expression profiling.

Genome-wide profiling using single-nucleotide polymorphism arrays identifies novel chromosomal imbalances in pediatric glioblastomas

Available data on genetic events in pediatric grade IV astrocytomas (glioblastoma [pGBM]) are scarce. This has traditionally been a major impediment in understanding the pathogenesis of this tumor and in developing ways for more effective management. Our aim is to chart DNA copy number aberrations (CNAs) and get insight into genetic pathways involved in pGBM. Using the Illumina Infinium Human-1 bead-chip-array (100K single-nucleotide polymorphisms [SNPs]), we genotyped 18 pediatric and 6 adult GBMs. Results were compared to BAC-array profiles harvested on 16 of the same pGBM, to an independent data set of 9 pediatric high-grade astrocytomas (HGAs) analyzed on Affymetrix 250K-SNP arrays, and to existing data sets on HGAs. CNAs were additionally validated by real-time qPCR in a set of genes in pGBM. Our results identify with nonrandom clustering of CNAs in several novel, previously not reported, genomic regions, suggesting that alterations in tumor suppressors and genes involved in the regulation of RNA processing and the cell cycle are major events in the pathogenesis of pGBM. Most regions were distinct from CNAs in aGBMs and show an unexpectedly low frequency of genetic amplification and homozygous deletions and a high frequency of loss of heterozygosity for a high-grade I rapidly dividing tumor. This first, complete, high-resolution profiling of the tumor cell genome fills an important gap in studies on pGBM. It ultimately guides the mapping of oncogenic networks unique to pGBM, identification of the related therapeutic predictors and targets, and development of more effective therapies. It further shows that, despite commonalities in a few CNAs, pGBM and aGBMs are two different diseases.

Detection of recurrent copy number alterations in the genome: taking among-subject heterogeneity seriously

BACKGROUND

Alterations in the number of copies of genomic DNA that are common or recurrent among diseased individuals are likely to contain disease-critical genes. Unfortunately, defining common or recurrent copy number alteration (CNA) regions remains a challenge. Moreover, the heterogeneous nature of many diseases requires that we search for common or recurrent CNA regions that affect only some subsets of the samples (without knowledge of the regions and subsets affected), but this is neglected by most methods.

RESULTS

We have developed two methods to define recurrent CNA regions from aCGH data. Our methods are unique and qualitatively different from existing approaches: they detect regions over both the complete set of arrays and alterations that are common only to some subsets of the samples (i.e., alterations that might characterize previously unknown groups); they use probabilities of alteration as input and return probabilities of being a common region, thus allowing researchers to modify thresholds as needed; the two parameters of the methods have an immediate, straightforward, biological interpretation. Using data from previous studies, we show that we can detect patterns that other methods miss and that researchers can modify, as needed, thresholds of immediate interpretability and develop custom statistics to answer specific research questions.

CONCLUSION

These methods represent a qualitative advance in the location of recurrent CNA regions, highlight the relevance of population heterogeneity for definitions of recurrence, and can facilitate the clustering of samples with respect to patterns of CNA. Ultimately, the methods developed can become important tools in the search for genomic regions harboring disease-critical genes.

Highly dynamic and sex-specific expression of microRNAs during early ES cell differentiation

Embryonic stem (ES) cells are pluripotent cells derived from the inner cell mass of the mammalian blastocyst. Cellular differentiation entails loss of pluripotency and gain of lineage-specific characteristics. However, the molecular controls that govern the differentiation process remain poorly understood. We have characterized small RNA expression profiles in differentiating ES cells as a model for early mammalian development. High-throughput 454 pyro-sequencing was performed on 19-30 nt RNAs isolated from undifferentiated male and female ES cells, as well as day 2 and 5 differentiating derivatives. A discrete subset of microRNAs (miRNAs) largely dominated the small RNA repertoire, and the dynamics of their accumulation could be readily used to discriminate pluripotency from early differentiation events. Unsupervised partitioning around meloids (PAM) analysis revealed that differentiating ES cell miRNAs can be divided into three expression clusters with highly contrasted accumulation patterns. PAM analysis afforded an unprecedented level of definition in the temporal fluctuations of individual members of several miRNA genomic clusters. Notably, this unravelled highly complex post-transcriptional regulations of the key pluripotency miR-290 locus, and helped identify miR-293 as a clear outlier within this cluster. Accordingly, the miR-293 seed sequence and its predicted cellular targets differed drastically from those of the other abundant cluster members, suggesting that previous conclusions drawn from whole miR-290 over-expression need to be reconsidered. Our analysis in ES cells also uncovered a striking male-specific enrichment of the miR-302 family, which share the same seed sequence with most miR-290 family members. Accordingly, a miR-302 representative was strongly enriched in embryonic germ cells derived from primordial germ cells of male but not female mouse embryos. Identifying the chromatin remodelling and E2F-dependent transcription repressors Ari4a and Arid4b as additional targets of miR-302 and miR-290 supports and possibly expands a model integrating possible overlapping functions of the two miRNA families in mouse cell totipotency during early development. This study demonstrates that small RNA sampling throughout early ES cell differentiation enables the definition of statistically significant expression patterns for most cellular miRNAs. We have further shown that the transience of some of these miRNA patterns provides highly discriminative markers of particular ES cell states during their differentiation, an approach that might be broadly applicable to the study of early mammalian development.

RJaCGH: Bayesian analysis of aCGH arrays for detecting copy number changes and recurrent regions

SUMMARY

Several methods have been proposed to detect copy number changes and recurrent regions of copy number variation from aCGH, but few methods return probabilities of alteration explicitly, which are the direct answer to the question 'is this probe/region altered?' RJaCGH fits a Non-Homogeneous Hidden Markov model to the aCGH data using Markov Chain Monte Carlo with Reversible Jump, and returns the probability that each probe is gained or lost. Using these probabilites, recurrent regions (over sets of individuals) of copy number alteration can be found.

AVAILABILITY

RJaCGH is available as an R package from CRAN repositories (e.g. http://cran.r-project.org/web/packages).

Rhomboid domain containing 2 (RHBDD2): a novel cancer-related gene over-expressed in breast cancer

In the course of breast cancer global gene expression studies, we identified an uncharacterized gene known as RHBDD2 (Rhomboid domain containing 2) to be markedly over-expressed in primary tumors from patients with recurrent disease. In this study, we identified RHBDD2 mRNA and protein expression significantly elevated in breast carcinomas compared with normal breast samples as analyzed by SAGE (n=46) and immunohistochemistry (n=213). Interestingly, specimens displaying RHBDD2 over-expression were predominantly advanced stage III breast carcinomas (p=0.001). Western-blot, RT-PCR and cDNA sequencing analyses allowed us to identify two RHBDD2 alternatively spliced mRNA isoforms expressed in breast cancer cell lines. We further investigated the occurrence and frequency of gene amplification and over-expression affecting RHBDD2 in 131 breast samples. RHBDD2 gene amplification was detected in 21% of 98 invasive breast carcinomas analyzed. However, no RHBDD2 amplification was detected in normal breast tissues (n=17) or breast benign lesions (n=16) (p=0.014). Interestingly, siRNA-mediated silencing of RHBDD2 expression results in a decrease of MCF7 breast cancer cells proliferation compared with the corresponding controls (p=0.001). In addition, analysis of publicly available gene expression data showed a strong association between high RHBDD2 expression and decreased overall survival (p=0.0023), relapse-free survival (p=0.0013), and metastasis-free interval (p=0.006) in patients with primary ER-negative breast carcinomas. In conclusion, our findings suggest that RHBDD2 over-expression behaves as an indicator of poor prognosis and may play a role facilitating breast cancer progression.

Merkel cell carcinoma of the skin: pathological and molecular evidence for a causative role of MCV in oncogenesis

Merkel cell carcinoma (MCC), a skin tumour with neuroendocrine features, was recently found to be associated with a new type of human polyomavirus, called Merkel cell virus (MCV). We investigated the specificity of this association as well as a causal role of MCV in oncogenesis. DNA and RNA from ten cases of MCC were analysed using PCR and RT-PCR. DNA from 1241 specimens of a wide range of human tumours was also analysed. The DIPS technique was used to identify the integration locus of viral DNA sequences. Array CGH was performed to analyse structural alterations of the cell genome. MCV DNA sequences were found in all ten cases of MCC and in none of the 1241 specimens of other tumour types. Clonal integration of MCV into the host genome was seen in all MCC cases and was checked by FISH in one case. A recurrent pattern of conserved viral sequences which encompassed the replication origin, the small tumour (ST), and the 5' part of the large tumour (LT) antigen DNA sequences was observed. Both ST and LT viral sequences were found to be significantly expressed in all MCCs. Neither recurrent site of integration nor alteration of cellular genes located near the viral sequences was observed. The tight association of MCV with MCC, the clonal pattern of MCV integration, and the expression of the viral oncoproteins strongly support a causative role for MCV in the tumour process. This information will help the development of novel approaches for the assessment and therapy of MCC and biologically related tumours.

Overall genomic pattern is a predictor of outcome in neuroblastoma

PURPOSE

For a comprehensive overview of the genetic alterations of neuroblastoma, their association and clinical significance, we conducted a whole-genome DNA copy number analysis.

PATIENTS AND METHODS

A series of 493 neuroblastoma (NB) samples was investigated by array-based comparative genomic hybridization in two consecutive steps (224, then 269 patients).

RESULTS

Genomic analysis identified several types of profiles. Tumors presenting exclusively whole-chromosome copy number variations were associated with excellent survival. No disease-related death was observed in this group. In contrast, tumors with any type of segmental chromosome alterations characterized patients with a high risk of relapse. Patients with both numerical and segmental abnormalities clearly shared the higher risk of relapse of segmental-only patients. In a multivariate analysis, taking into account the genomic profile, but also previously described individual genetic and clinical markers with prognostic significance, the presence of segmental alterations with (HR, 7.3; 95% CI, 3.7 to 14.5; P < .001) or without MYCN amplification (HR, 4.5; 95% CI, 2.4 to 8.4; P < .001) was the strongest predictor of relapse; the other significant variables were age older than 18 months (HR, 1.8; 95% CI, 1.2 to 2.8; P = .004) and stage 4 (HR, 1.8; 95% CI, 1.2 to 2.7; P = .005). Finally, within tumors showing segmental alterations, stage 4, age, MYCN amplification, 1p and 11q deletions, and 1q gain were independent predictors of decreased overall survival.

CONCLUSION

The analysis of the overall genomic pattern, which probably unravels particular genomic instability mechanisms rather than the analysis of individual markers, is essential to predict relapse in NB patients. It adds critical prognostic information to conventional markers and should be included in future treatment stratification.

Hormone escape is associated with genomic instability in a human prostate cancer model

Lack of hormone dependency in prostate cancers is an irreversible event that occurs through generation of genomic instability induced by androgen deprivation. Indeed, the cytogenetic profile of hormone-dependent (HD) prostate cancer remains stable as long as it received a hormone supply, whereas the profile of hormone-independent (HID) variants acquired new and various alterations. This is demonstrated here using a HD xenografted model of a human prostate cancer, PAC120, transplanted for 11 years into male nude mice and 4 HID variants obtained by surgical castration. Cytogenetic analysis, done by karyotype, FISH, CGH and array-CGH, shows that PAC120 at early passage presents numerous chromosomal alterations. Very few additional alterations were found between the 5th and 47th passages, indicating the stability of the parental tumor. HID variants largely maintained the core of chromosomal alterations of PAC120 - losses at 6q, 7p, 12q, 15q and 17q sites. However, each HID variant displayed a number of new alterations, almost all being specific to each variant and very few shared by all. None of the HID had androgen receptor mutations. Our study indicates that hormone castration is responsible for genomic instability generating new cytogenetic abnormalities susceptible to alter the properties of cancer cell associated with tumor progression, such as increased cell survival and ability to metastasize.

Frequent PTEN genomic alterations and activated phosphatidylinositol 3-kinase pathway in basal-like breast cancer cells

INTRODUCTION

Basal-like carcinomas (BLCs) and human epidermal growth factor receptor 2 overexpressing (HER2+) carcinomas are the subgroups of breast cancers that have the most aggressive clinical behaviour. In contrast to HER2+ carcinomas, no targeted therapy is currently available for the treatment of patients with BLCs. In order to discover potential therapeutic targets, we aimed to discover deregulated signalling pathways in human BLCs.

METHODS

In this study, we focused on the oncogenic phosphatidylinositol 3-kinase (PI3K) pathway in 13 BLCs, and compared it with a control series of 11 hormonal receptor negative- and grade III-matched HER2+ carcinomas. The two tumour populations were first characterised by immunohistochemistry and gene expression. The PI3K pathway was then investigated by gene copy-number analysis, gene expression profiling and at a proteomic level using reverse-phase protein array technology and tissue microarray. The effects of the PI3K inhibition pathway on proliferation and apoptosis was further analysed in three human basal-like cell lines.

RESULTS

The PI3K pathway was found to be activated in BLCs and up-regulated compared with HER2+ tumours as shown by a significantly increased activation of the downstream targets Akt and mTOR (mammalian target of rapamycin). BLCs expressed significantly lower levels of the tumour suppressor PTEN and PTEN levels were significantly negatively correlated with Akt activity within that population. PTEN protein expression correlated significantly with PTEN DNA copy number and more importantly, reduced PTEN DNA copy numbers were observed specifically in BLCs. Similar to human samples, basal-like cell lines exhibited an activation of PI3K/Akt pathway and low/lack PTEN expression. Both PI3K and mTOR inhibitors led to basal-like cell growth arrest. However, apoptosis was specifically observed after PI3K inhibition.

CONCLUSIONS

These data provide insight into the molecular pathogenesis of BLCs and implicate the PTEN-dependent activated Akt signalling pathway as a potential therapeutic target for the management of patients with poor prognosis BLCs.

Characterization of amplicons in neuroblastoma: high-resolution mapping using DNA microarrays, relationship with outcome, and identification of overexpressed genes

Somatically acquired chromosomal imbalances are a key feature of neuroblastoma, a heterogeneous pediatric solid tumor. Among these alterations, genomic amplification targeting the MYCN oncogene and observed in about 25-30% of the cases, strongly correlates with advanced stage and poor outcome. In this work, we have used BAC and SNP arrays as well as gene expression arrays to characterize amplifications in neuroblastoma. Eighty-eight distinct BACs defining high-level amplification events were identified in 65 samples, including 43 tumors and 22 cell lines. Although the highest recurrence was observed on chromosome 2, clones on chromosomes 8, 12, 16, and 17 also revealed genomic amplification in several samples. A detailed analysis of the 2p22-2p25 MYCN containing region indicated highly complex patterns in a number of cases. Coamplifications involving MYCN and other regions were explored by FISH in three cell lines. High-resolution arrays then allowed us to further refine the mapping of 25 amplicons in 19 samples, either reducing the size of a single continuous amplicon or increasing the complexity by highlighting multiple noncontiguous regions of amplification. Combined analysis of gene expression profiling and array-CGH data indicated that 12 to 25% of the genes that are targeted by genomic amplification are actually over-expressed in tumor cells, several of them having already been implicated in cancer. Finally, our results suggest that the presence of amplicons localized outside of chromosome 2, in addition to MYCN amplification, may be linked to a particularly severe outcome in neuroblastoma patients.

SIGMA2: a system for the integrative genomic multi-dimensional analysis of cancer genomes, epigenomes, and transcriptomes

Background

High throughput microarray technologies have afforded the investigation of genomes, epigenomes, and transcriptomes at unprecedented resolution. However, software packages to handle, analyze, and visualize data from these multiple 'omics disciplines have not been adequately developed.

Results

Here, we present SIGMA², a system for the integrative genomic multi-dimensional analysis of cancer genomes, epigenomes, and transcriptomes. Multi-dimensional datasets can be simultaneously visualized and analyzed with respect to each dimension, allowing combinatorial integration of the different assays belonging to the different 'omics.

Conclusion

The identification of genes altered at multiple levels such as copy number, loss of heterozygosity (LOH), DNA methylation and the detection of consequential changes in gene expression can be concertedly performed, establishing SIGMA² as a novel tool to facilitate the high throughput systems biology analysis of cancer.

Oncogenomic analysis of mycosis fungoides reveals major differences with Sezary syndrome

Mycosis fungoides (MF), the most common cutaneous T-cell lymphoma, is a malignancy of mature, skin-homing T cells. Sézary syndrome (Sz) is often considered to represent a leukemic phase of MF. In this study, the pattern of numerical chromosomal alterations in MF tumor samples was defined using array-based comparative genomic hybridization (CGH); simultaneously, gene expression was analyzed using microarrays. Highly recurrent chromosomal alterations in MF include gain of 7q36, 7q21-7q22 and loss of 5q13 and 9p21. The pattern characteristic of MF differs markedly from chromosomal alterations observed in Sz. Integration of data from array-based CGH and gene-expression analysis yielded several candidate genes with potential relevance in the pathogenesis of MF. We confirmed that the FASTK and SKAP1 genes, residing in loci with recurrent gain, demonstrated increased expression. The RB1 and DLEU1 tumor suppressor genes showed diminished expression associated with loss. In addition, it was found that the presence of chromosomal alterations on 9p21, 8q24, and 1q21-1q22 was associated with poor prognosis in patients with MF. This study provides novel insight into genetic alterations underlying MF. Furthermore, our analysis uncovered genomic differences between MF and Sz, which suggest that the molecular pathogenesis and therefore therapeutic requirements of these cutaneous T-cell lymphomas may be distinct.

Genomic changes in progression of low-grade gliomas

Using a one-megabase BAC-based array comparative genomic hybridization technique (aCGH), we have investigated a series of 16 low-grade gliomas (LGGs) and their subsequent progression to higher-grade malignancies. The most frequent chromosome imbalances in primary tumors were gains of chromosomes 7q, 8q, and 22q, and losses of chromosomes 1p, 13q, and 19q. In tumor progression, gains of chromosomes 11q, 7q, 20q, and 21q, and losses of chromosomes 9p, including CDKN2A locus, 19q, 14q, 1p, and 6q were the most frequent genomic disequilibria. Progressive tumors were more imbalanced than primary tumors in terms of altered chromosomal arms (3.8 vs. 6.6 in mean abnormal chromosomal arm) and altered BACs (17 vs. 21%). Interestingly, putative novel candidate genes associated with glioma progression were identified, in particular DOCK8, PTPRD, CER1, TPHO, DHFR, MSH3, ETS1, ACACA, and CSE1L.

Integrated genomic and transcriptomic analysis of ductal carcinoma in situ of the breast

PURPOSE

To gain insight into genomic and transcriptomic subtypes of ductal carcinomas in situ of the breast (DCIS).

EXPERIMENTAL DESIGN

We did a combined phenotypic and genomic analysis of a series of 57 DCIS integrated with gene expression profile analysis for 26 of the 57 cases.

RESULTS

Thirty-two DCIS exhibited a luminal phenotype; 21 were ERBB2 positive, and 4 were ERBB2/estrogen receptor (ER) negative with 1 harboring a bona fide basal-like phenotype. Based on a CGH analysis, genomic types were identified in this series of DCIS with the 1q gain/16q loss combination observed in 3 luminal DCIS, the mixed amplifier pattern including all ERBB2, 12 luminal and 2 ERBB2(-)/ER(-) DCIS, and the complex copy number alteration profile encompassing 14 luminal and 1 ERBB2(-)/ER(-) DCIS. Eight cases (8 of 57; 14%) presented a TP53 mutation, all being amplifiers. Unsupervised analysis of gene expression profiles of 26 of the 57 DCIS showed that luminal and ERBB2-amplified, ER-negative cases clustered separately. We further investigated the effect of high and low copy number changes on gene expression. Strikingly, amplicons but also low copy number changes especially on 1q, 8q, and 16q in DCIS regulated the expression of a subset of genes in a very similar way to that recently described in invasive ductal carcinomas.

CONCLUSIONS

These combined approaches show that the molecular heterogeneity of breast ductal carcinomas exists already in in situ lesions and further indicate that DCIS and invasive ductal carcinomas share genomic alterations with a similar effect on gene expression profile.

MD-SeeGH: a platform for integrative analysis of multi-dimensional genomic data

Background

Recent advances in global genomic profiling methodologies have enabled multi-dimensional characterization of biological systems. Complete analysis of these genomic profiles require an in depth look at parallel profiles of segmental DNA copy number status, DNA methylation state, single nucleotide polymorphisms, as well as gene expression profiles. Due to the differences in data types it is difficult to conduct parallel analysis of multiple datasets from diverse platforms.

Results

To address this issue, we have developed an integrative genomic analysis platform MD-SeeGH, a software tool that allows users to rapidly and directly analyze genomic datasets spanning multiple genomic experiments. With MD-SeeGH, users have the flexibility to easily update datasets in accordance with new genomic builds, make a quality assessment of data using the filtering features, and identify genetic alterations within single or across multiple experiments. Multiple sample analysis in MD-SeeGH allows users to compare profiles from many experiments alongside tracks containing detailed localized gene information, microRNA, CpG islands, and copy number variations.

Conclusion

MD-SeeGH is a new platform for the integrative analysis of diverse microarray data, facilitating multiple profile analyses and group comparisons.

High-resolution array comparative genomic hybridization of single micrometastatic tumor cells

Only few selected cancer cells drive tumor progression and are responsible for therapy resistance. Their specific genomic characteristics, however, are largely unknown because high-resolution genome analysis is currently limited to DNA pooled from many cells. Here, we describe a protocol for array comparative genomic hybridization (array CGH), which enables the detection of DNA copy number changes in single cells. Combining a PCR-based whole genome amplification method with arrays of highly purified BAC clones we could accurately determine known chromosomal changes such as trisomy 21 in single leukocytes as well as complex genomic imbalances of single cell line cells. In single T47D cells aberrant regions as small as 1-2 Mb were identified in most cases when compared to non-amplified DNA from 10(6) cells. Most importantly, in single micrometastatic cancer cells isolated from bone marrow of breast cancer patients, we retrieved and confirmed amplifications as small as 4.4 and 5 Mb. Thus, high-resolution genome analysis of single metastatic precursor cells is now possible and may be used for the identification of novel therapy target genes.

BAC array CGH distinguishes mutually exclusive alterations that define clinicogenetic subtypes of gliomas

The pathological classification of gliomas constitutes a critical step of the clinical management of patients, yet it is frequently challenging. To assess the relationship between genetic abnormalities and clinicopathological characteristics, we have performed a genetic and clinical analysis of a series of gliomas. A total of 112 gliomas were analyzed by comparative genomic hybridization on a BAC array with a 1 megabase resolution. Altered regions were identified and correlation analysis enabled to retrieve significant associations and exclusions. Whole chromosomes (chrs) 1p and 19q losses with centromeric breakpoints and EGFR high level amplification were found to be mutually exclusive, permitting identification of 3 distinct, nonoverlapping groups of tumors with striking clinicopathological differences. Type A tumors with chrs 1p and 19q co-deletion exhibited an oligodendroglial phenotype and a longer patient survival. Type B tumors were characterized by EGFR amplification. They harbored a WHO high grade of malignancy and a short patient survival. Finally, type C tumors displayed none of the previous patterns but the presence of chr 7 gain, chr 9p deletion and/or chr 10 loss. It included astrocytic tumors in patients younger than in type B and whose prognosis was highly dependent upon the number of alterations. A multivariate analysis based on a Cox model shows that age, WHO grade and genomic type provide complementary prognostic informations. Finally, our results highlight the potential of a whole-genome analysis as an additional diagnostic in cases of unclear conventional genetic findings.

Chromosome 1p loss evaluation in anaplastic oligodendrogliomas

The chromosome (chr) 1p deletion is a favorable biomarker in oligodendroglial tumors and is even more powerful a marker when combined with chr 19q loss. As a result, the 1p deletion is taken into account more and more in clinical trials and the management of patients. However, the laboratory technique implemented for detection of this biomarker has been a topic of debate. To illustrate the usefulness of evaluating multiple loci, we here report two anaplastic oligodendrogliomas that were investigated using fluorescent in situ hybridization (FISH) and bacterial artificial chromosome (BAC)-array-based comparative genomic hybridization (aCGH). Indeed, segmental analysis using FISH, limited to chr 1p36 was unable to discriminate between complete and partial deletions of chrs 1p. However, complete and partial deletions of 1p are reported to have distinct clinical outcomes. Our results illustrate that aCGH (or other multiple loci technologies) provide complementary information to single locus technologies such as FISH because multiple loci technologies can evaluate the extent of the chr 1p deletion.

ITALICS: an algorithm for normalization and DNA copy number calling for Affymetrix SNP arrays

MOTIVATION

Affymetrix SNP arrays can be used to determine the DNA copy number measurement of 11 000-500 000 SNPs along the genome. Their high density facilitates the precise localization of genomic alterations and makes them a powerful tool for studies of cancers and copy number polymorphism. Like other microarray technologies it is influenced by non-relevant sources of variation, requiring correction. Moreover, the amplitude of variation induced by non-relevant effects is similar or greater than the biologically relevant effect (i.e. true copy number), making it difficult to estimate non-relevant effects accurately without including the biologically relevant effect.

RESULTS

We addressed this problem by developing ITALICS, a normalization method that estimates both biological and non-relevant effects in an alternate, iterative manner, accurately eliminating irrelevant effects. We compared our normalization method with other existing and available methods, and found that ITALICS outperformed these methods for several in-house datasets and one public dataset. These results were validated biologically by quantitative PCR.

AVAILABILITY

The R package ITALICS (ITerative and Alternative normaLIzation and Copy number calling for affymetrix Snp arrays) has been submitted to Bioconductor.

Dedifferentiated liposarcoma of the pleura mimicking a malignant solitary fibrous tumor and associated with dedifferentiated liposarcoma of the mediastinum: usefulness of cytogenetic and molecular genetic analyses

Dedifferentiated liposarcoma of the pleura is an extremely rare malignancy mimicking a variety of tumors, such as other sarcomas, mesothelioma, and malignant solitary fibrous tumor of the pleura. Liposarcoma of the pleura can be combined with mediastinal involvement, and in most cases it may be impossible to be certain where the primary tumor originated. In this report, we describe a very rare occurence of a dedifferentiated liposarcoma of the pleura in a 76-year-old woman associated with a distinct second dedifferentiated liposarcoma of the mediastinum. Histologically, the pleural tumor demonstrated spindle cells arranged in a fascicular pattern, whereas the mediastinal tumor was mostly adipocytic with small areas of spindle cells. Vimentin and protein S100 were focally expressed by the tumor cells. The differential diagnosis of the pleural mass included malignant solitary fibrous tumor. Cytogenetic analysis showed supernumerary ring chromosomes in the pleural tumor, as well as strong amplification of MDM2 and CDK4 genes in both tumors. Array comparative genomic hybridization showed amplifications of chromosome arms 6q, 12q, and 15q, shared by both tumors and strongly pointing to a common origin.

Characterization of breast cancer by array comparative genomic hybridization

Cancer progression is due to the accumulation of recurrent genomic alterations that induce growth advantage and clonal expansion. Most of these genomic changes can be detected using the array comparative genomic hybridization (CGH) technique. The accurate classification of these genomic alterations is expected to have an important impact on translational and basic research. Here we review recent advances in CGH technology used in the characterization of different features of breast cancer. First, we present bioinformatics methods that have been developed for the analysis of CGH arrays; next, we discuss the use of array CGH technology to classify tumor stages and to identify and stratify subgroups of patients with different prognoses and clinical behaviors. We finish our review with a discussion of how CGH arrays are being used to identify oncogenes, tumor suppressor genes, and breast cancer susceptibility genes.

High-resolution mapping of DNA breakpoints to define true recurrences among ipsilateral breast cancers

BACKGROUND

To distinguish new primary breast cancers from true recurrences, pangenomic analyses of DNA copy number alterations (CNAs) using single-nucleotide polymorphism arrays have proven useful.

METHODS

The pangenomic profiles of 22 pairs of primary breast carcinoma (ductal or lobular) and ipsilateral breast cancers from the same patients were analyzed. Hierarchical clustering was performed using CNAs and DNA breakpoint information. A partial identity score developed using DNA breakpoint information was used to quantify partial identities between two tumors. The nature of ipsilateral breast cancers (true recurrence vs new primary tumor) as defined using the clustering methods and the partial identity score was compared with that based on clinical characteristics. Metastasis-free survival was compared among patients with primary tumors and true recurrences as defined using the partial identity score and by clinical characteristics. All statistical tests were two-sided.

RESULTS

All methods agreed on the nature of ipsilateral breast cancers for 14 pairs of samples. For five pairs, the clinical definition disagreed with both clustering methods. For three pairs, the two clustering methods were discordant and the one using DNA breakpoints agreed with the clinical definition. The partial identity score confirmed the nature of ipsilateral breast cancers as defined by clustering of DNA breakpoints in 21 of 22 pairs. The difference in metastasis-free survival of patients with new primary tumors and those with true recurrences was not statistically significant when tumors were defined based on clinical and histologic characteristics (5-year metastasis-free survival: 76%, 95% confidence interval [CI] = 52% to 100% for new primary tumors and 38%, 95% CI = 17% to 83% for true recurrences; P = .18; new primary tumor vs true recurrence, hazard ratio = 2.8, 95% CI = 0.6 to 13.7), but the difference was statistically significant when tumors were defined using the partial identity score (5-year metastasis-free survival: 100% for new primary tumors and 29%, 95% CI = 11% to 78% for true recurrences; P = .01).

CONCLUSIONS

DNA breakpoint information more often agreed with the clinical determination than CNAs in this population. The partial identity score, which was calculated based on DNA breakpoints, allows statistical discrimination between new primary tumors and true recurrences that could outperform the clinical determination in terms of prognosis.

In neuroblastic tumours, Schwann cells do not harbour the genetic alterations of neuroblasts but may nevertheless share the same clonal origin

Neuroblastic tumours are composed of variable proportions of neuroblasts and Schwann cells. Whether both components share a common neoplastic origin is highly debated and discrepant results have been reported about the presence of tumour-related genetic alterations in Schwann cells. We have used X-methylation analysis and array-CGH to investigate contiguous Schwannian and neuroblastic areas in tumours with a nodular pattern. A skewed X inactivation was observed in four out of five stromal components. Interestingly, in these four cases, the X-inactivation profiles of the neuroblastic components were identical to the matched stromal areas. However, whereas all neuroblastic areas displayed chromosomal imbalances, no alteration was found in any Schwann cell components. Similarly, no alteration was observed in a series of 19 tumours with a single stroma-rich component, which occasionally exhibited a skewed X-inactivation pattern (3/17 informative tumours). Altogether, this indicates that most stroma-rich tumours display a polyclonal proliferation and that Schwann cells do not derive from neuroblasts. However, in tumours with both stroma-rich and -poor components, our results suggest that cells from both areas share a common progenitor.

An annotation infrastructure for the analysis and interpretation of Affymetrix exon array data

An annotation database (X:MAP) and BioConductor/R package (exonmap) have been developed to support fine-grained analysis of exon array data.

Affymetrix exon arrays contain probesets intended to target every known and predicted exon in the entire genome, posing significant challenges for high-throughput genome-wide data analysis. X:MAP , an annotation database, and exonmap , a BioConductor/R package, are designed to support fine-grained analysis of exon array data. The system supports the application of standard statistical techniques, prior to the use of genome scale annotation to provide gene-, transcript- and exon-level summaries and visualization tools.

Tumor genomic profiling and TP53 germline mutation analysis of first-degree relative familial gliomas

About 5% of gliomas occur in a familial context, which suggests a genetic origin, but the predisposing molecular factors remain unknown in most cases. A series of nine familial gliomas were characterized with 1-megabase resolution BAC array-based comparative genomic hybridization (aCGH) together with germline sequence analysis of TP53. This series was compared with a literature series of familial gliomas and a personal series of sporadic gliomas, analyzed by chromosome CGH and aCGH, respectively. No significant difference was noted between the three populations in terms of clinical characteristics, pathologic features, and the most frequent chromosomal alterations, including loss of 1p, 10p, 10q, 13q, and 19q, and gain of 7p, 7q, 16p, 18q, 19p, 19q, 20p, and 22q. However, a genomic region located in 6q was more frequently gained in our series of familial as compared to sporadic gliomas (P=0.028). A germline TP53 mutation was observed in 1/9 cases, which suggests Li-Fraumeni syndrome. Interestingly, the Pro allele in the codon 72 of TP53 was observed in 5/9 tumors. Although familial and sporadic gliomas share very similar cytogenetic quantitative patterns, aCGH is a promising technique for the detection of small genomic differences of potential significance.

ACTuDB, a new database for the integrated analysis of array-CGH and clinical data for tumors

In recent years, an increasing number of projects have investigated tumor genome structure, using microarray-based techniques like array comparative genomic hybridization (array-CGH) or single nucleotide polymorphism (SNP) arrays. The forthcoming studies have to integrate these former results and compare their findings to the existing sets of copy number data for validation. These sets also form the basis from which many comparative retrospective analyses can be carried out. Nevertheless, exploitation of this mass of data relies on a homogeneous preparation of copy number data, which will make it possible to compare them together, and their integration into a unified bioinformatics environment with ad hoc analysis tools and interfaces. To our knowledge, no such data integration has been proposed yet. Therefore the biologists and clinicians involved in cancer research urgently need such an integrative tool, which motivated us to undertake the construction of a database for array-CGH and other DNA copy number data for tumors (ACTuDB). When available, the associated clinical, transcriptome and loss of heterozygosity data were also integrated into ACTuDB. ACTuDB contains currently about 1500 genomic profiles for tumors and cell lines for the bladder, brain, breast, colon, liver, lymphoma, neuroblastoma, mouth and pancreas, together with data for replication timing experiments. The CGH array data were processed, using ad hoc algorithms (probe mapping, breakpoint detection, gain or loss status assignment and visualization) developed at Institut Curie. The database is available from http://bioinfo.curie.fr/actudb/ and can be browsed with a user-friendly interface. This database will be a useful resource for the genomic profiling of tumors, a field of highly active research. We invite research groups involved in tumor genome profiling to submit their data to ACTuDB.

The C-MYB locus is involved in chromosomal translocation and genomic duplications in human T-cell acute leukemia (T-ALL), the translocation defining a new T-ALL subtype in very young children

The C-Myb transcription factor is essential for hematopoiesis, including in the T-cell lineage. The C-Myb locus is a common site of retroviral insertional mutagenesis, however no recurrent genomic involvement has been reported in human malignancies. Here, we identified 2 types of genomic alterations involving the C-MYB locus at 6q23 in human T-cell acute leukemia (T-ALL). First, we found a reciprocal translocation, t(6;7)(q23;q34), that juxtaposed the TCRB and C-MYB loci (n = 6 cases). Second, a genome-wide copy-number analysis by array-based comparative genomic hybridization (array-CGH) identified short somatic duplications that include C-MYB (MYB(dup), n = 13 cases of 84 T-ALL, 15%). Expression analysis, including allele-specific approaches, showed stronger C-MYB expression in the MYB-rearranged cases compared with other T-ALLs, and a dramatically skewed C-MYB allele expression in the TCRB-MYB cases, which suggests that a translocation-driven deregulated expression may overcome a cellular attempt to down-regulate C-MYB. Strikingly, profiling of the T-ALLs by clinical, genomic, and large-scale gene expression analyses shows that the TCRB-MYB translocation defines a new T-ALL subtype associated with a very young age for T-cell leukemia (median, 2.2 years) and with a proliferation/mitosis expression signature. By contrast, the MYB(dup) alteration was associated with the previously defined T-ALL subtypes.

SIGMA: a system for integrative genomic microarray analysis of cancer genomes

BACKGROUND

The prevalence of high resolution profiling of genomes has created a need for the integrative analysis of information generated from multiple methodologies and platforms. Although the majority of data in the public domain are gene expression profiles, and expression analysis software are available, the increase of array CGH studies has enabled integration of high throughput genomic and gene expression datasets. However, tools for direct mining and analysis of array CGH data are limited. Hence, there is a great need for analytical and display software tailored to cross platform integrative analysis of cancer genomes.

RESULTS

We have created a user-friendly java application to facilitate sophisticated visualization and analysis such as cross-tumor and cross-platform comparisons. To demonstrate the utility of this software, we assembled array CGH data representing Affymetrix SNP chip, Stanford cDNA arrays and whole genome tiling path array platforms for cross comparison. This cancer genome database contains 267 profiles from commonly used cancer cell lines representing 14 different tissue types.

CONCLUSION

In this study we have developed an application for the visualization and analysis of data from high resolution array CGH platforms that can be adapted for analysis of multiple types of high throughput genomic datasets. Furthermore, we invite researchers using array CGH technology to deposit both their raw and processed data, as this will be a continually expanding database of cancer genomes. This publicly available resource, the System for Integrative Genomic Microarray Analysis (SIGMA) of cancer genomes, can be accessed at http://sigma.bccrc.ca.

Visualization of genomic aberrations using Affymetrix SNP arrays

MOTIVATION

DNA copy number aberrations are frequently found in different types of cancer. Recent developments of microarray-based approaches have broadened the knowledge on number and structure of such aberrations. High-density single nucleotide polymorphism (SNP) microarrays provide an extremely high resolution with up to 500,000 SNPs per genome. Owing to the enormous amount of data the detection of common aberrations in large datasets is a great challenge. We describe a novel open source software tool--IdeogramBrowser--which was specifically designed for use with the Affymetrix SNP arrays. It provides an interactive karyotypic visualization of multiple aberration profiles and direct links to GeneCards. Visualization of consensus regions together with gene representation allows the explorative assessment of the data.

AVAILABILITY

IdeogramBrowser and its source code are freely available under a creative commons license and can be obtained from http://www.informatik.uni-ulm.de/ni/staff/HKestler/ideo/. IdeogramBrowser is a platform independent Java application.