Weighted clustering of called array CGH data

Genetic Profiling of Colorectal Carcinomas of Patients with Primary Sclerosing Cholangitis and Inflammatory Bowel Disease

BACKGROUND

Patients with primary sclerosing cholangitis (PSC) and inflammatory bowel disease (IBD) run a 10-fold increased risk of developing colorectal cancer (CRC) compared to patients with IBD only. The aim of this study was to perform an extensive screen of known carcinogenic genomic alterations in patients with PSC-IBD, and to investigate whether such changes occur already in nondysplastic mucosa.

METHODS

Archival cancer tissue and nondysplastic mucosa from resection specimens of 19 patients with PSC-IBD-CRC were characterized, determining DNA copy-number variations, microsatellite instability (MSI), mutations on 48 cancer genes, and CpG island methylator phenotype (CIMP). Genetic profiles were compared with 2 published cohorts of IBD-associated CRC (IBD-CRC; n = 11) and sporadic CRC (s-CRC; n = 100).

RESULTS

Patterns of chromosomal aberrations in PSC-IBD-CRC were similar to those observed in IBD-CRC and s-CRC, MSI occurred only once. Mutation frequencies were comparable between the groups, except for mutations in KRAS, which were less frequent in PSC-IBD-CRC (5%) versus IBD-CRC (38%) and s-CRC (31%; P = .034), and in APC, which were less frequent in PSC-IBD-CRC (5%) and IBD-CRC (0%) versus s-CRC (50%; P < .001). Cases of PSC-IBD-CRC were frequently CIMP positive (44%), at similar levels to cases of s-CRC (34%; P = .574) but less frequent than in cases with IBD-CRC (90%; P = .037). Similar copy number aberrations and mutations were present in matched cancers and adjacent mucosa in 5/15 and 7/11 patients, respectively.

CONCLUSIONS

The excess risk of CRC in patients with PSC-IBD was not explained by copy number aberrations, mutations, MSI, nor CIMP status, in cancer tissue, nor in adjacent mucosa. These findings set the stage for further exome-wide and epigenetic studies.

A pipeline for copy number profiling of single circulating tumor cells to assess intra-patient tumor heterogeneity

Intrapatient tumour heterogeneity is likely a major determinant of clinical outcome in cancer patients. To assess heterogeneity in a minimally invasive manner, methods to perform single circulating tumour cell (CTC) genomics at high resolution are necessary. However, due to the rarity of CTCs, development of such methods is challenging. Here, we developed a modular single CTC analysis pipeline to assess intrapatient heterogeneity by copy number (CN) profiling. To optimize this pipeline, spike‐in experiments using MCF‐7 breast cancer cells were performed. The VyCAP puncher system was used to isolate single cells. The quality of whole genome amplification (WGA) products generated by REPLI‐g and Ampli1™ methods, as well as the results from the Illumina Truseq and the Ampli1™ LowPass library preparation techniques, was compared. Moreover, a bioinformatic pipeline was designed to generate CN profiles from single CTCs. The optimal combination of Ampli1™ WGA and Illumina Truseq library preparation was successfully validated on patient‐derived CTCs. In conclusion, we developed a novel modular pipeline to isolate single CTCs and subsequently generate detailed patient‐derived CN profiles that allow assessment of intrapatient heterogeneity in future studies.

Genetic analysis of single disseminated tumor cells in the lymph nodes and bone marrow of patients with head and neck squamous cell carcinoma

Considering the limited information on the biology and molecular characteristics of disseminated tumor cells (DTCs) in head and neck squamous cell carcinoma (HNSCC), we examined the genomic alterations in DTCs from HNSCCs and their potential clinical relevance. To analyze both the lymphatic and hematogenous routes of tumor cell dissemination, we investigated samples from lymph nodes (LNs) and bone marrow (BM) of 49 patients using immunofluorescence double staining for epithelial cells expressing cytokeratin 18 (KRT18) and/or epithelial cell adhesion molecules (EpCAM, CD326). The identified marker‐positive cells were isolated by micromanipulation followed by single‐cell whole‐genome amplification and metaphase‐based comparative genomic hybridization (mCGH) to determine genome‐wide copy number alterations. The findings were correlated with clinical parameters and follow‐up data. We detected chromosomal aberrations in KRT18‐ and EpCAM‐positive cells from both compartments; BM‐derived cells showed a significantly higher percentage of aberrant genome (PAG) per cell than cells detected in LNs. No significant association was found between DTC data and clinical follow‐up. Genomic profiling of BM‐DTCs revealed genomic alterations typical for HNSCC, suggesting hematogenous dissemination of subclones around the time of surgery. In contrast, DTC data in LNs revealed that several marker‐positive cells were not of malignant origin, indicating the presence of epithelial glandular inclusions in parts of the processed neck LN samples. Therefore, DTC detection of LNs in the neck based only on epithelial markers is not advisable and requires detection of chromosomal instability (CIN), gene mutations, or additional markers, which have yet to be identified. Nevertheless, our investigation paves the way for larger studies to focus on HNSCC BM‐DTCs with high‐resolution methods to gain deeper insights into the biology of hematogenous metastasis in this cancer.

Identification of Deregulated Pathways, Key Regulators, and Novel miRNA-mRNA Interactions in HPV-Mediated Transformation

Next to a persistent infection with high-risk human papillomavirus (HPV), molecular changes are required for the development of cervical cancer. To identify which molecular alterations drive carcinogenesis, we performed a comprehensive and longitudinal molecular characterization of HPV-transformed keratinocyte cell lines. Comparative genomic hybridization, mRNA, and miRNA expression analysis of four HPV-containing keratinocyte cell lines at eight different time points was performed. Data was analyzed using unsupervised hierarchical clustering, integrated longitudinal expression analysis, and pathway enrichment analysis. Biological relevance of identified key regulatory genes was evaluated in vitro and dual-luciferase assays were used to confirm predicted miRNA-mRNA interactions. We show that the acquisition of anchorage independence of HPV-containing keratinocyte cell lines is particularly associated with copy number alterations. Approximately one third of differentially expressed mRNAs and miRNAs was directly attributable to copy number alterations. Focal adhesion, TGF-beta signaling, and mTOR signaling pathways were enriched among these genes. PITX2 was identified as key regulator of TGF-beta signaling and inhibited cell growth in vitro, most likely by inducing cell cycle arrest and apoptosis. Predicted miRNA-mRNA interactions miR-221-3p_BRWD3, miR-221-3p_FOS, and miR-138-5p_PLXNB2 were confirmed in vitro. Integrated longitudinal analysis of our HPV-induced carcinogenesis model pinpointed relevant interconnected molecular changes and crucial signaling pathways in HPV-mediated transformation.

IBD-Associated Dysplastic Lesions Show More Chromosomal Instability Than Sporadic Adenomas

BACKGROUND

Patients with longstanding inflammatory bowel disease (IBD; ie, ulcerative colitis and Crohn's disease) have an increased risk of colorectal cancer (CRC). Due to ongoing inflammation, IBD-associated dysplastic lesions can develop. These lesions have an increased risk to progress to cancer compared with sporadic adenomas, which are also found in these patients. Differentiating between these 2 types of dysplasia remains challenging, both clinically and histologically, while treatment strategies may differ. Therefore, the aim of this study was to investigate molecular alterations associated with colorectal dysplasia to cancer progression in IBD and evaluate to what extent these alterations differ from sporadic adenomas.

METHODS

DNA copy number aberrations and mutation analyses of 48 genes were performed by next-generation sequencing in 43 IBD-associated dysplastic lesions, 30 of which were dysplastic and 13 of which were cancers. Results were compared with existing DNA copy number and mutation data from 118 sporadic adenomas and 24 sporadic cancers.

RESULTS

Inflammatory bowel disease-associated dysplastic lesions harbor patterns of DNA copy number aberrations comparable to carcinomas, which are rare in sporadic adenomas. TP53 mutation was the most frequent mutation observed in IBD-associated dysplastic lesions and in cancers. FBXW7 was mutated significantly more often in IBD-associated dysplastic lesions than in sporadic adenomas.

CONCLUSIONS

Inflammatory bowel disease-associated dysplastic lesions show more DNA copy number aberrations than sporadic adenomas. TP53 and FBXW7 mutations appear to be involved in the development of IBD-associated dysplastic lesions and cancer. These findings indicate that IBD-associated dysplastic lesions are more genomically unstable, possibly reflecting a faster progression toward cancer.

Molecular heterogeneity in human papillomavirus-dependent and -independent vulvar carcinogenesis

Vulvar squamous cell carcinoma (VSCC) and precancerous vulvar intraepithelial neoplasia (VIN) can develop through human papillomavirus (HPV)-dependent and -independent pathways, indicating a heterogeneous disease. Only a minority of VIN progress, but current clinicopathological classifications are insufficient to predict the cancer risk. Here we analyzed copy number alterations (CNA) to assess the molecular heterogeneity of vulvar lesions in relation to HPV and cancer risk. HPV-status and CNA by means of whole-genome next-generation shallow-sequencing were assessed in VSCC and VIN. The latter included VIN of women with associated VSCC (VINVSCC ) and women who did not develop VSCC during follow-up (VINnoVSCC ). HPV-testing resulted in 41 HPV-positive (16 VINVSCC , 14 VINnoVSCC , and 11 VSCC) and 24 HPV-negative (11 VINVSCC and 13 VSCC) lesions. HPV-positive and -negative VSCC showed a partially overlapping pattern of recurrent CNA, including frequent gains of 3q and 8q. In contrast, amplification of 11q13/cyclinD1 was exclusively found in HPV-negative lesions. HPV-negative VINVSCC had less CNA than HPV-negative VSCC (P = .009), but shared chromosome 8 alterations. HPV-positive VINnoVSCC had less CNA than VINVSCC (P = .022). Interestingly, 1pq gain was detected in 81% of HPV-positive VINVSCC and only in 21% of VINnoVSCC (P = .001). In conclusion, HPV-dependent and -independent vulvar carcinogenesis is characterized by distinct CNA patterns at the VIN stage, while more comparable patterns are present at the cancer stage. Cancer risk in VIN seems to be reflected by the extent of CNA, in particular chromosome 1 gain in HPV-positive cases.

Genomic profiling of stage II and III colon cancers reveals APC mutations to be associated with survival in stage III colon cancer patients

Tumor profiling of DNA alterations, i.e. gene point mutations, somatic copy number aberrations (CNAs) and structural variants (SVs), improves insight into the molecular pathology of cancer and clinical outcome. Here, associations between genomic aberrations and disease recurrence in stage II and III colon cancers were investigated. A series of 114 stage II and III microsatellite stable colon cancer samples were analyzed by high-resolution array-comparative genomic hybridization (array-CGH) to detect CNAs and CNA-associated chromosomal breakpoints (SVs). For 60 of these samples mutation status of APC, TP53, KRAS, PIK3CA, FBXW7, SMAD4, BRAF and NRAS was determined using targeted massive parallel sequencing. Loss of chromosome 18q12.1-18q12.2 occurred more frequently in tumors that relapsed than in relapse-free tumors (p < 0.001; FDR = 0.13). In total, 267 genes were recurrently affected by SVs (FDR < 0.1). CNAs and SVs were not associated with disease-free survival (DFS). Mutations in APC and TP53 were associated with increased CNAs. APC mutations were associated with poor prognosis in (5-fluorouracil treated) stage III colon cancers (p = 0.005; HR = 4.1), an effect that was further enhanced by mutations in MAPK pathway (KRAS, NRAS, BRAF) genes. We conclude that among multiple genomic alterations in CRC, strongest associations with clinical outcome were observed for common mutations in APC.

Fast Bayesian Inference of Copy Number Variants using Hidden Markov Models with Wavelet Compression

By integrating Haar wavelets with Hidden Markov Models, we achieve drastically reduced running times for Bayesian inference using Forward-Backward Gibbs sampling. We show that this improves detection of genomic copy number variants (CNV) in array CGH experiments compared to the state-of-the-art, including standard Gibbs sampling. The method concentrates computational effort on chromosomal segments which are difficult to call, by dynamically and adaptively recomputing consecutive blocks of observations likely to share a copy number. This makes routine diagnostic use and re-analysis of legacy data collections feasible; to this end, we also propose an effective automatic prior. An open source software implementation of our method is available at http://schlieplab.org/Software/HaMMLET/ (DOI: 10.5281/zenodo.46262). This paper was selected for oral presentation at RECOMB 2016, and an abstract is published in the conference proceedings.

A Systems Biology Interpretation of Array Comparative Genomic Hybridization (aCGH) Data through Phylogenetics

Array Comparative Genomic Hybridization (aCGH) is a rapid screening technique to detect gene deletions and duplications, providing an overview of chromosomal aberrations throughout the entire genome of a tumor, without the need for cell culturing. However, the heterogeneity of aCGH data obfuscates existing methods of data analysis. Analysis of aCGH data from a systems biology perspective or in the context of total aberrations is largely absent in the published literature. We present here a novel alternative to the functional analysis of aCGH data using the phylogenetic paradigm that is well-suited to high dimensional datasets of heterogeneous nature, but has not been widely adapted to aCGH data. Maximum parsimony phylogenetic analysis sorts out genetic data through the simplest presentation of the data on a cladogram, a graphical evolutionary tree, thus providing a powerful and efficient method for aCGH data analysis. For example, the cladogram models the multiphasic changes in the cancer genome and identifies shared early mutations in the disease progression, providing a simple yet powerful means of aCGH data interpretation. As such, applying maximum parsimony phylogenetic analysis to aCGH results allows for the differentiation between drivers and passenger genes aberrations in cancer specimens. In addition to offering a novel methodology to analyze aCGH results, we present here a crucial software suite that we wrote to carry out the analysis. In a broader context, we wish to underscore that phylogenetic analysis of aCGH data is a non-parametric method that circumvents the pitfalls and frustrations of standard analytical techniques that rely on parametric statistics. Organizing the data in a cladogram as explained in this research article provides insights into the disease common aberrations, as well as the disease subtypes and their shared aberrations (the synapomorphies) of each subtype. Hence, we report the method and make the software suite publicly and freely available at http://software.phylomcs.com so that researchers can test alternative and innovative approaches to the analysis of aCGH data.

HCsnip: An R Package for Semi-supervised Snipping of the Hierarchical Clustering Tree

Hierarchical clustering (HC) is one of the most frequently used methods in computational biology in the analysis of high-dimensional genomics data. Given a data set, HC outputs a binary tree leaves of which are the data points and internal nodes represent clusters of various sizes. Normally, a fixed-height cut on the HC tree is chosen, and each contiguous branch of data points below that height is considered as a separate cluster. However, the fixed-height branch cut may not be ideal in situations where one expects a complicated tree structure with nested clusters. Furthermore, due to lack of utilization of related background information in selecting the cutoff, induced clusters are often difficult to interpret. This paper describes a novel procedure that aims to automatically extract meaningful clusters from the HC tree in a semi-supervised way. The procedure is implemented in the R package HCsnip available from Bioconductor. Rather than cutting the HC tree at a fixed-height, HCsnip probes the various way of snipping, possibly at variable heights, to tease out hidden clusters ensconced deep down in the tree. The cluster extraction process utilizes, along with the data set from which the HC tree is derived, commonly available background information. Consequently, the extracted clusters are highly reproducible and robust against various sources of variations that "haunted" high-dimensional genomics data. Since the clustering process is guided by the background information, clusters are easy to interpret. Unlike existing packages, no constraint is placed on the data type on which clustering is desired. Particularly, the package accepts patient follow-up data for guiding the cluster extraction process. To our knowledge, HCsnip is the first package that is able to decomposes the HC tree into clusters with piecewise snipping under the guidance of patient time-to-event information. Our implementation of the semi-supervised HC tree snipping framework is generic, and can be combined with other algorithms that operate on detected clusters.

Semi-supervised adaptive-height snipping of the hierarchical clustering tree

BACKGROUND

In genomics, hierarchical clustering (HC) is a popular method for grouping similar samples based on a distance measure. HC algorithms do not actually create clusters, but compute a hierarchical representation of the data set. Usually, a fixed height on the HC tree is used, and each contiguous branch of samples below that height is considered a separate cluster. Due to the fixed-height cutting, those clusters may not unravel significant functional coherence hidden deeper in the tree. Besides that, most existing approaches do not make use of available clinical information to guide cluster extraction from the HC. Thus, the identified subgroups may be difficult to interpret in relation to that information.

RESULTS

We develop a novel framework for decomposing the HC tree into clusters by semi-supervised piecewise snipping. The framework, called guided piecewise snipping, utilizes both molecular data and clinical information to decompose the HC tree into clusters. It cuts the given HC tree at variable heights to find a partition (a set of non-overlapping clusters) which does not only represent a structure deemed to underlie the data from which HC tree is derived, but is also maximally consistent with the supplied clinical data. Moreover, the approach does not require the user to specify the number of clusters prior to the analysis. Extensive results on simulated and multiple medical data sets show that our approach consistently produces more meaningful clusters than the standard fixed-height cut and/or non-guided approaches.

CONCLUSIONS

The guided piecewise snipping approach features several novelties and advantages over existing approaches. The proposed algorithm is generic, and can be combined with other algorithms that operate on detected clusters. This approach represents an advancement in several regards: (1) a piecewise tree snipping framework that efficiently extracts clusters by snipping the HC tree possibly at variable heights while preserving the HC tree structure; (2) a flexible implementation allowing a variety of data types for both building and snipping the HC tree, including patient follow-up data like survival as auxiliary information. The data sets and R code are provided as supplementary files. The proposed method is available from Bioconductor as the R-package HCsnip.

Genomic landscape of metastatic colorectal cancer

Response to drug therapy in individual colorectal cancer (CRC) patients is associated with tumour biology. Here we describe the genomic landscape of tumour samples of a homogeneous well-annotated series of patients with metastatic CRC (mCRC) of two phase III clinical trials, CAIRO and CAIRO2. DNA copy number aberrations of 349 patients are determined. Within three treatment arms, 194 chromosomal subregions are associated with progression-free survival (PFS; uncorrected single-test P-values <0.005). These subregions are filtered for effect on messenger RNA expression, using an independent data set from The Cancer Genome Atlas which returned 171 genes. Three chromosomal regions are associated with a significant difference in PFS between treatment arms with or without irinotecan. One of these regions, 6q16.1–q21, correlates in vitro with sensitivity to SN-38, the active metabolite of irinotecan. This genomic landscape of mCRC reveals a number of DNA copy number aberrations associated with response to drug therapy.

Response to drug therapy in colorectal cancer (CRC) patients has been associated with tumour heterogeneity. Here the authors analyse DNA copy number aberrations in primary tumours from CRC patients and identify genetic variants that influence drug response.

Chromosome 20 aberrations at the diploid-aneuploid transition in sporadic colorectal cancer

DNA aneuploid sublines in sporadic colorectal cancers (CRCs) are quite frequent (about 85%) and likely the consequence of chromosomal instability and DNA copy number aberrations (CNAs). In order to gain insight into the mechanisms of the diploid-aneuploid transition in CRCs, we compared the CNA status in both diploid and aneuploid sublines. We used fresh/frozen material from 17 aneuploid CRCs, which was separated into 17 DNA diploid and 17 aneuploid sublines using enrichment of the epithelial component by multiparameter flow cytometry and sorting. CNA status of both sublines was obtained by array comparative genomic hybridization. The DNA diploid sublines from the aneuploid CRCs showed already CNAs, in particular, gains at 20 p and 20 q. The same aberrations were detected at increased frequencies in the corresponding DNA aneuploid sublines. Moreover, the very frequent gains/losses of chromosomes 4, 7, 8, 13, 15, and 18 in the DNA aneuploid sublines were absent or rare in the DNA diploid sublines from the same sporadic aneuploid CRCs. The comparison of the DNA diploid and aneuploid sublines from aneuploid CRCs suggests that 20 p and 20 q gains may play a role in the diploid-aneuploid transition. The 20 q chromosomal arm appears of particular interest since it harbors several genes implicated in chromosomal instability.

Deletion at 6q24.2-26 predicts longer survival of high-grade serous epithelial ovarian cancer patients

Standard treatments for advanced high-grade serous ovarian carcinomas (HGSOCs) show significant side-effects and provide only short-term survival benefits due to disease recurrence. Thus, identification of novel prognostic and predictive biomarkers is urgently needed. We have used 42 paraffin-embedded HGSOCs, to evaluate the utility of DNA copy number alterations, as potential predictors of clinical outcome. Copy number-based unsupervised clustering stratified HGSOCs into two clusters of different immunohistopathological features and survival outcome (HR = 0.15, 95%CI = 0.03-0.81; Padj = 0.03). We found that loss at 6q24.2-26 was significantly associated with the cluster of longer survival independently from other confounding factors (HR = 0.06, 95%CI = 0.01-0.43, Padj = 0.005). The prognostic value of this deletion was validated in two independent series, one consisting of 36 HGSOCs analyzed by fluorescent in situ hybridization (P = 0.04) and another comprised of 411 HGSOCs from the Cancer Genome Atlas study (TCGA) (HR = 0.67, 95%CI = 0.48-0.93, Padj = 0.019). In addition, we confirmed the association of low expression of the genes from the region with longer survival in 799 HGSOCs (HR = 0.74, 95%CI = 0.61-0.90, log-rank P = 0.002) and 675 high-FIGO stage HGSOCs (HR = 0.76, 95%CI = 0.61-0.96, log-rank P = 0.02) available from the online tool KM-plotter. Finally, by integrating copy number, RNAseq and survival data of 296 HGSOCs from TCGA we propose a few candidate genes that can potentially explain the association. Altogether our findings indicate that the 6q24.2-26 deletion is an independent marker of favorable outcome in HGSOCs with potential clinical value as it can be analyzed by FISH on tumor sections and guide the selection of patients towards more conservative therapeutic strategies in order to reduce side-effects and improve quality of life.

Spatial and temporal evolution of distal 10q deletion, a prognostically unfavorable event in diffuse low-grade gliomas

Background

The disease course of patients with diffuse low-grade glioma is notoriously unpredictable. Temporal and spatially distinct samples may provide insight into the evolution of clinically relevant copy number aberrations (CNAs). The purpose of this study is to identify CNAs that are indicative of aggressive tumor behavior and can thereby complement the prognostically favorable 1p/19q co-deletion.

Results

Genome-wide, 50 base pair single-end sequencing was performed to detect CNAs in a clinically well-characterized cohort of 98 formalin-fixed paraffin-embedded low-grade gliomas. CNAs are correlated with overall survival as an endpoint. Seventy-five additional samples from spatially distinct regions and paired recurrent tumors of the discovery cohort were analyzed to interrogate the intratumoral heterogeneity and spatial evolution. Loss of 10q25.2-qter is a frequent subclonal event and significantly correlates with an unfavorable prognosis. A significant correlation is furthermore observed in a validation set of 126 and confirmation set of 184 patients. Loss of 10q25.2-qter arises in a longitudinal manner in paired recurrent tumor specimens, whereas the prognostically favorable 1p/19q co-deletion is the only CNA that is stable across spatial regions and recurrent tumors.

Conclusions

CNAs in low-grade gliomas display extensive intratumoral heterogeneity. Distal loss of 10q is a late onset event and a marker for reduced overall survival in low-grade glioma patients. Intratumoral heterogeneity and higher frequencies of distal 10q loss in recurrences suggest this event is involved in outgrowth to the recurrent tumor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0471-6) contains supplementary material, which is available to authorized users.

The genetic basis for inactivation of Wnt pathway in human osteosarcoma

BACKGROUND

Osteosarcoma is a highly genetically unstable tumor with poor prognosis. We performed microarray-based comparative genomic hybridization (aCGH), transcriptome sequencing (RNA-seq), and pathway analysis to gain a systemic view of the pathway alterations of osteosarcoma.

METHODS

aCGH experiments were carried out on 10 fresh osteosarcoma samples. The output data (Gene Expression Omnibus Series accession number GSE19180) were pooled with published aCGH raw data (GSE9654) to determine recurrent copy number changes. These were analyzed using Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis to identify altered pathways in osteosarcoma. Transcriptome sequencing of six osteosarcomas was performed to detect the expression profile of Wnt signaling pathway genes. Protein expression of WNT1, β-catenin, c-myc, and cyclin D1 in the Wnt pathway was detected by immunohistochemistry (IHC) in an independent group of 46 osteosarcoma samples.

RESULTS

KEGG pathway analysis identified frequent deletions of Wnt and other Wnt signaling pathway genes. At the mRNA level, transcriptome sequencing found reduced levels of mRNA expression of Wnt signaling pathway transcripts. While WNT1 protein expression was detected by IHC in 69.6% (32/46) of the osteosarcomas, no β-catenin protein was detected in the nucleus. β-catenin protein expression was, however, detected in the membrane and cytoplasm of 69.6% (32/46) of the osteosarcomas. c-myc protein expression was detected in only 47.8% (22/46) and cyclin D1 protein expression in 52.2% (24/46) of osteosarcoma samples. Kaplan-Meier survival analysis showed that WNT1-negative patients had a trend towards longer disease free survival than WNT1-positive patients. Interestingly, in WNT1-negative patients, those who were also cyclin D1-negative had significantly longer disease free survival than cyclin D1-positive patients. However, there was no significant association between any of the investigated proteins and overall survival of human osteosarcoma patients.

CONCLUSIONS

Frequent deletions of Wnt and other Wnt signaling pathway genes suggest that the Wnt signaling pathway is genetically inactivated in human osteosarcoma.

Stability-based comparison of class discovery methods for DNA copy number profiles

MOTIVATION

Array-CGH can be used to determine DNA copy number, imbalances in which are a fundamental factor in the genesis and progression of tumors. The discovery of classes with similar patterns of array-CGH profiles therefore adds to our understanding of cancer and the treatment of patients. Various input data representations for array-CGH, dissimilarity measures between tumor samples and clustering algorithms may be used for this purpose. The choice between procedures is often difficult. An evaluation procedure is therefore required to select the best class discovery method (combination of one input data representation, one dissimilarity measure and one clustering algorithm) for array-CGH. Robustness of the resulting classes is a common requirement, but no stability-based comparison of class discovery methods for array-CGH profiles has ever been reported.

RESULTS

We applied several class discovery methods and evaluated the stability of their solutions, with a modified version of Bertoni's [Formula: see text]-based test [1]. Our version relaxes the assumption of independency required by original Bertoni's [Formula: see text]-based test. We conclude that Minimal Regions of alteration (a concept introduced by [2]) for input data representation, sim [3] or agree [4] for dissimilarity measure and the use of average group distance in the clustering algorithm produce the most robust classes of array-CGH profiles.

AVAILABILITY

The software is available from http://bioinfo.curie.fr/projects/cgh-clustering. It has also been partly integrated into "Visualization and analysis of array-CGH"(VAMP)[5]. The data sets used are publicly available from ACTuDB [6].

DNA copy number aberrations in endobronchial lesions: a validated predictor for cancer

We recently identified a DNA copy number aberration (CNA)-based classifier, including changes at 3p26.3-p11.1, 3q26.2-29, and 6p25.3-24.3, as a risk predictor for cancer in individuals presenting with endobronchial squamous metaplasia. The current study was set out to validate the prediction accuracy of this classifier in an independent series of endobronchial squamous metaplastic and dysplastic lesions. The study included 36 high-risk subjects who had endobronchial lesions of various histological grades that were identified and biopsied by autofluorescence bronchoscopy and were subjected to arrayCGH in a nested case-control design. Of the 36 patients, 12 had a carcinoma in situ or invasive carcinoma at the same site at follow-up (median 11 months, range 4-24), while 24 controls remained cancer free (78 months, range 21-142). The previously defined CNA-based classifier demonstrated 92% (95% CI 77% to 98%) accuracy for cancer (in situ) prediction. All nine subjects with CNA-based classifier-positive endobronchial lesions at baseline experienced cancer outcome, whereas all 24 controls and 3 cases were classified as being low risk. In conclusion, CNAs prove to be a highly accurate biomarker for assessing the progression risk of endobronchial squamous metaplastic and dysplastic lesions. This classifier could assist in selecting subjects with endobronchial lesions who might benefit from more aggressive therapeutic intervention or surveillance.

DNA copy number profiling reveals extensive genomic loss in hereditary BRCA1 and BRCA2 ovarian carcinomas

Background:

Few studies have attempted to characterise genomic changes occurring in hereditary epithelial ovarian carcinomas (EOCs) and inconsistent results have been obtained. Given the relevance of DNA copy number alterations in ovarian oncogenesis and growing clinical implications of the BRCA-gene status, we aimed to characterise the genomic profiles of hereditary and sporadic ovarian tumours.

Methods:

High-resolution array Comparative Genomic Hybridisation profiling of 53 familial (21 BRCA1, 6 BRCA2 and 26 non-BRCA1/2) and 15 sporadic tumours in combination with supervised and unsupervised analysis was used to define common and/or specific copy number features.

Results:

Unsupervised hierarchical clustering did not stratify tumours according to their familial or sporadic condition or to their BRCA1/2 mutation status. Common recurrent changes, spanning genes potentially fundamental for ovarian carcinogenesis, regardless of BRCA mutations, and several candidate subtype-specific events were defined. Despite similarities, greater contribution of losses was revealed to be a hallmark of BRCA1 and BRCA2 tumours.

Conclusion:

Somatic alterations occurring in the development of familial EOCs do not differ substantially from the ones occurring in sporadic carcinomas. However, some specific features like extensive genomic loss observed in BRCA1/2 tumours may be of clinical relevance helping to identify BRCA-related patients likely to respond to PARP inhibitors.

Aggressiveness of human melanoma xenograft models is promoted by aneuploidy-driven gene expression deregulation

Melanoma is a devastating skin cancer characterized by distinct biological subtypes. Besides frequent mutations in growth- and survival-promoting genes like BRAF and NRAS, melanomas additionally harbor complex non-random genomic alterations. Using an integrative approach, we have analysed genomic and gene expression changes in human melanoma cell lines (N=32) derived from primary tumors and various metastatic sites and investigated the relation to local growth aggressiveness as xenografts in immuno-compromised mice (N=22). Although the vast majority >90% of melanoma models harbored mutations in either BRAF or NRAS, significant differences in subcutaneous growth aggressiveness became obvious. Unsupervised clustering revealed that genomic alterations rather than gene expression data reflected this aggressive phenotype, while no association with histology, stage or metastatic site of the original melanoma was found. Genomic clustering allowed separation of melanoma models into two subgroups with differing local growth aggressiveness in vivo. Regarding genes expressed at significantly altered levels between these subgroups, a surprising correlation with the respective gene doses (>85% accordance) was found. Genes deregulated at the DNA and mRNA level included well-known cancer genes partly already linked to melanoma (RAS genes, PTEN, AURKA, MAPK inhibitors Sprouty/Spred), but also novel candidates like SIPA1 (a Rap1GAP). Pathway mining further supported deregulation of Rap1 signaling in the aggressive subgroup e.g. by additional repression of two Rap1GEFs. Accordingly, siRNA-mediated down-regulation of SIPA1 exerted significant effects on clonogenicity, adherence and migration in aggressive melanoma models. Together our data suggest that an aneuploidy-driven gene expression deregulation drives local aggressiveness in human melanoma.

Specific genomic regions are differentially affected by copy number alterations across distinct cancer types, in aggregated cytogenetic data

Background

Regional genomic copy number alterations (CNA) are observed in the vast majority of cancers. Besides specifically targeting well-known, canonical oncogenes, CNAs may also play more subtle roles in terms of modulating genetic potential and broad gene expression patterns of developing tumors. Any significant differences in the overall CNA patterns between different cancer types may thus point towards specific biological mechanisms acting in those cancers. In addition, differences among CNA profiles may prove valuable for cancer classifications beyond existing annotation systems.

Principal Findings

We have analyzed molecular-cytogenetic data from 25579 tumors samples, which were classified into 160 cancer types according to the International Classification of Disease (ICD) coding system. When correcting for differences in the overall CNA frequencies between cancer types, related cancers were often found to cluster together according to similarities in their CNA profiles. Based on a randomization approach, distance measures from the cluster dendrograms were used to identify those specific genomic regions that contributed significantly to this signal. This approach identified 43 non-neutral genomic regions whose propensity for the occurrence of copy number alterations varied with the type of cancer at hand. Only a subset of these identified loci overlapped with previously implied, highly recurrent (hot-spot) cytogenetic imbalance regions.

Conclusions

Thus, for many genomic regions, a simple null-hypothesis of independence between cancer type and relative copy number alteration frequency can be rejected. Since a subset of these regions display relatively low overall CNA frequencies, they may point towards second-tier genomic targets that are adaptively relevant but not necessarily essential for cancer development.

Genetic characterization of mesenchymal, clear cell, and dedifferentiated chondrosarcoma

Clear cell, mesenchymal, and dedifferentiated chondrosarcoma are rare, cartilaginous tumors with limited treatment options other than surgery. Conventional chondrosarcomas have been extensively studied at the genetic level, but for rare chondrosarcoma subtypes, this is merely restricted to case reports. Information on the genetics of rare chondrosarcomas may provide insight into the etiology of these specific disease subtypes and possible alternative treatment strategies. Therefore, the aim of this study was to genetically characterize this subset of rare tumors. Using array CGH, we gathered genomic information of 30 rare cartilaginous tumors. In addition, we constructed tissue microarrays with 2 mm cores of 23 clear cell, 23 mesenchymal, and 45 dedifferentiated chondrosarcomas, in triplicate. Using immunohistochemistry, we investigated expression of R132H IDH1, and p53 and retinoblastoma pathways. Results were verified and further investigated with a methylation assay and MLPA for CDKN2A/p16, and IDH1/2, and TP53 mutation analysis. Array-CGH showed numerous genomic alterations in all subtypes. However, only a limited number of recurrent alterations were detected, none of which seemed to be associated with the subtypes. The IDH1/2, p53, and retinoblastoma pathways were affected in 0, 9, and 95% of clear cell chondrosarcomas, in 0, 39, and 70% in mesenchymal chondrosarcomas, and in 50, 59, and 85% of dedifferentiated chondrosarcomas, respectively. Our results suggest an important role for the retinoblastoma pathway in all three rare chondrosarcoma subtypes investigated.

Genomic aberrations relate early and advanced stage ovarian cancer

Background

Because of the distinct clinical presentation of early and advanced stage ovarian cancer, we aim to clarify whether these disease entities are solely separated by time of diagnosis or whether they arise from distinct molecular events.

Methods

Sixteen early and sixteen advanced stage ovarian carcinomas, matched for histological subtype and differentiation grade, were included. Genomic aberrations were compared for each early and advanced stage ovarian cancer by array comparative genomic hybridization. To study how the aberrations correlate to the clinical characteristics of the tumors we clustered tumors based on the genomic aberrations.

Results

The genomic aberration patterns in advanced stage cancer equalled those in early stage, but were more frequent in advanced stage (p = 0.012). Unsupervised clustering based on genomic aberrations yielded two clusters that significantly discriminated early from advanced stage (p = 0.001), and that did differ significantly in survival (p = 0.002). These clusters however did give a more accurate prognosis than histological subtype or differentiation grade.

Conclusion

This study indicates that advanced stage ovarian cancer either progresses from early stage or from a common precursor lesion but that they do not arise from distinct carcinogenic molecular events. Furthermore, we show that array comparative genomic hybridization has the potential to identify clinically distinct patients.

Electronic supplementary material

The online version of this article (doi:10.1007/s13402-012-0077-5) contains supplementary material, which is available to authorized users.

Losses of chromosome 5q and 14q are associated with favorable clinical outcome of patients with gastric cancer

PURPOSE

To improve the clinical outcome of patients with gastric cancer, intensified combination strategies are currently in clinical development, including combinations of more extensive surgery, (neo)adjuvant chemotherapy, and radiotherapy. The present study used DNA copy number profiling to identify subgroups of patients with different clinical outcomes. We hypothesize that, by identification of subgroups, individual treatment strategies can be selected to improve clinical outcome and to reduce unnecessary treatment toxicity for patients with gastric cancer.

EXPERIMENTAL DESIGN

DNA from 206 gastric cancer patients was isolated and analyzed by genomewide array comparative genomic hybridization. DNA copy number profiles were correlated with lymph node status and patient survival. In addition, heat shock protein 90 (HSP90) expression was analyzed and correlated with survival in 230 gastric cancer patients.

RESULTS

Frequent (>20%) DNA copy number gains and losses were observed on several chromosomal regions. Losses on 5q11.2-q31.3 and 14q32.11-q32.33 (14% of patients) were correlated with good clinical outcome in univariate and multivariate analyses, with a median disease-free survival interval of 9.2 years. In addition, loss of expression of HSP90, located on chromosome 14q32.2, was correlated with better patient survival.

CONCLUSION

Genomewide DNA copy number profiling allowed the identification of a subgroup of gastric cancer patients, marked by losses on chromosomes 5q11.2-q31.3 and 14q32.11-q32.33 or low HSP90 protein expression, with an excellent clinical outcome after surgery alone. We hypothesize that this subgroup of patients most likely will not benefit from (neo)adjuvant systemic treatment and/or radiotherapy, whereas anti-HSP90 therapy may have clinical potential in patients with HSP90-expressing gastric cancer, pending validation in an independent dataset.

Chromosomal profiles of high-grade cervical intraepithelial neoplasia relate to duration of preceding high-risk human papillomavirus infection

High-grade cervical intraepithelial neoplasia (CIN2/3) represents a heterogeneous disease both with respect to clinical behavior and chromosomal aberrations detected. We hypothesized that the extent of chromosomal aberrations reflects the duration of their existence. Chromosomal profiles were determined of CIN3 of women with a known 5-year history of high-risk human papillomavirus virus (hrHPV) infection, in which duration of prior hrHPV infection was considered a proxy for duration of CIN3 existence. Eleven women had a <5 year preceding hrHPV infection (CIN3<5yrPHI) and 24 had a PHI lasting ≥5 years (CIN3≥5yrPHI). For comparison, six CIN3 adjacent to squamous cell carcinomas (CIN3-SCC), the corresponding SCCs, and six CIN1 were included. Unsupervised hierarchical clustering analysis of the chromosomal profiles revealed two clusters. One was characterized by a low number of chromosomal aberrations and included all CIN1, 81.8% of CIN3<5yrPHI and 33.3% of CIN3≥5yrPHI. Samples in the second cluster, displaying multiple aberrations, included 18.2% of CIN3<5yrPHI, 66.7% CIN3≥5yrPHI, all except one CIN3-SCC and all SCCs. The number of genomic aberrations increased according to lesion grade and also with longer duration of PHI. The increase in aberrations in CIN3≥5yrPHI compared to <5yrPHI was highly significant (p = 0.001), suggesting that CIN3≥5yrPHI represent more severe lesions. In conclusion, longer duration of preceding hrHPV infection is associated with an increased number of chromosomal aberrations. Hence, CIN3 with a longer duration of existence are likely more prone to have an increased short-term risk of cervical cancer.

Chipster: user-friendly analysis software for microarray and other high-throughput data

BACKGROUND

The growth of high-throughput technologies such as microarrays and next generation sequencing has been accompanied by active research in data analysis methodology, producing new analysis methods at a rapid pace. While most of the newly developed methods are freely available, their use requires substantial computational skills. In order to enable non-programming biologists to benefit from the method development in a timely manner, we have created the Chipster software.

RESULTS

Chipster (http://chipster.csc.fi/) brings a powerful collection of data analysis methods within the reach of bioscientists via its intuitive graphical user interface. Users can analyze and integrate different data types such as gene expression, miRNA and aCGH. The analysis functionality is complemented with rich interactive visualizations, allowing users to select datapoints and create new gene lists based on these selections. Importantly, users can save the performed analysis steps as reusable, automatic workflows, which can also be shared with other users. Being a versatile and easily extendable platform, Chipster can be used for microarray, proteomics and sequencing data. In this article we describe its comprehensive collection of analysis and visualization tools for microarray data using three case studies.

CONCLUSIONS

Chipster is a user-friendly analysis software for high-throughput data. Its intuitive graphical user interface enables biologists to access a powerful collection of data analysis and integration tools, and to visualize data interactively. Users can collaborate by sharing analysis sessions and workflows. Chipster is open source, and the server installation package is freely available.

Breast tumors from CHEK2 1100delC-mutation carriers: genomic landscape and clinical implications

INTRODUCTION

Checkpoint kinase 2 (CHEK2) is a moderate penetrance breast cancer risk gene, whose truncating mutation 1100delC increases the risk about twofold. We investigated gene copy-number aberrations and gene-expression profiles that are typical for breast tumors of CHEK2 1100delC-mutation carriers.

METHODS

In total, 126 breast tumor tissue specimens including 32 samples from patients carrying CHEK2 1100delC were studied in array-comparative genomic hybridization (aCGH) and gene-expression (GEX) experiments. After dimensionality reduction with CGHregions R package, CHEK2 1100delC-associated regions in the aCGH data were detected by the Wilcoxon rank-sum test. The linear model was fitted to GEX data with R package limma. Genes whose expression levels were associated with CHEK2 1100delC mutation were detected by the bayesian method.

RESULTS

We discovered four lost and three gained CHEK2 1100delC-related loci. These include losses of 1p13.3-31.3, 8p21.1-2, 8p23.1-2, and 17p12-13.1 as well as gains of 12q13.11-3, 16p13.3, and 19p13.3. Twenty-eight genes located on these regions showed differential expression between CHEK2 1100delC and other tumors, nominating them as candidates for CHEK2 1100delC-associated tumor-progression drivers. These included CLCA1 on 1p22 as well as CALCOCO1, SBEM, and LRP1 on 12q13. Altogether, 188 genes were differentially expressed between CHEK2 1100delC and other tumors. Of these, 144 had elevated and 44, reduced expression levels.Our results suggest the WNT pathway as a driver of tumorigenesis in breast tumors of CHEK2 1100delC-mutation carriers and a role for the olfactory receptor protein family in cancer progression. Differences in the expression of the 188 CHEK2 1100delC-associated genes divided breast tumor samples from three independent datasets into two groups that differed in their relapse-free survival time.

CONCLUSIONS

We have shown that copy-number aberrations of certain genomic regions are associated with CHEK2 mutation 1100delC. On these regions, we identified potential drivers of CHEK2 1100delC-associated tumorigenesis, whose role in cancer progression is worth investigating. Furthermore, poorer survival related to the CHEK2 1100delC gene-expression signature highlights pathways that are likely to have a role in the development of metastatic disease in carriers of the CHEK2 1100delC mutation.

DNA copy number alterations in endobronchial squamous metaplastic lesions predict lung cancer

RATIONALE

Autofluorescence bronchoscopy (AFB) is a valid strategy for detecting premalignant endobronchial lesions. However, no biomarker can reliably predict lung cancer risk of subjects with AFB-visualized premalignant lesions.

OBJECTIVES

The present study set out to identify AFB-visualized squamous metaplastic (SqM) lesions with malignant potential by DNA copy number profiling.

METHODS

Regular AFB examinations in 474 subjects at risk of lung cancer identified six subjects with SqM lesions at baseline, and carcinoma in situ or carcinoma (carcinoma in situ or greater) at the initial SqM site at follow-up bronchoscopy. These progressive SqM lesions were compared for immunostaining pattern and array comparative genomic hybridization-based chromosomal profiles with 23 SqM lesions of subjects who remained cancer-free. Specific DNA copy number alterations (CNAs) linked to cancer risk were identified and accuracy of CNAs to predict endobronchial cancer in this series was determined.

MEASUREMENTS AND MAIN RESULTS

At baseline, p53, p63, and Ki-67 immunostaining were not predictive for a differential clinical outcome of SqM lesions. The mean number of CNAs in baseline SqM of cases was significantly higher compared with control subjects (P < 0.01). Chromosomal regions significantly more frequently altered in SqM of cases were 3p26.3-p11.1, 3q26.2-q29, 9p13.3-p13.2, and 17p13.3-p11.2 (family-wise error rate <0.10). CNAs were specifically detected at the site of future cancer. In cases, baseline-detected CNAs persisted in subsequent biopsies taken from the initial site, and levels increased toward cancer progression. In this series, a model based on CNAs at 3p26.3-p11.1, 3q26.2-29, and 6p25.3-24.3 predicted cancer with 97% accuracy.

CONCLUSIONS

The data suggest that the presence of specific CNAs in SqM lesions predict endobronchial cancer.

Small bowel adenocarcinoma copy number profiles are more closely related to colorectal than to gastric cancers

BACKGROUND

Small bowel adenocarcinoma (SBA) is a rare cancer and consequently, the options for clinical trials are limited. As they are treated according to either a colorectal or a gastric cancer regimen and the molecular biology of a tumor is a pivotal determinant for therapy response, chromosomal copy number aberrations were compared with the colorectal and gastric adenocarcinomas.

MATERIALS AND METHODS

A total of 85 microsatellite stable (MSS) adenocarcinomas from the stomach, colorectum and small bowel were selected from existing array comparative genomic hybridization (aCGH) datasets. We compared the aCGH profiles of the three tumor sites by supervised analysis and hierarchical clustering.

RESULTS

Hierarchical clustering revealed substantial overlap of 27 SBA copy number profiles with matched colorectal adenocarcinomas but less overlap with profiles of gastric adenocarcinomas. DNA copy number aberrations located at chromosomes 1p36.3-p34.3, 4p15.3-q35.2, 9p24.3-p11.1, 13q13.2-q31.3 and 17p13.3-p13.2 were the strongest features discriminating SBAs and colorectal adenocarcinomas from gastric adenocarcinomas.

CONCLUSIONS

We show that MSS SBAs are more similar to colorectal than to gastric cancer, based on the 27 genome-wide DNA copy number profiles that are currently available. These molecular similarities provide added support for treatment of MSS small bowel cancers according to colorectal cancer regimens.

Genetic amplification of the vascular endothelial growth factor (VEGF) pathway genes, including VEGFA, in human osteosarcoma

BACKGROUND

Osteosarcoma is the most common primary tumor of bone. It is a highly vascular and extremely destructive malignancy that mainly affects children and young adults. The authors conducted microarray-based comparative genomic hybridization (aCGH) and pathway analyses to gain a systemic view of pathway alterations in the genetically altered genes.

METHODS

Recurrent amplified and deleted genes that were detected by aCGH were subjected to an analysis based on the Kyoto Encyclopedia of Genes and Genomes to identify the altered pathways. Among the enriched pathways, vascular endothelial growth factor (VEGF) pathway genes collectively were amplified, and alterations of this pathway were validated by fluorescence in situ hybridization (FISH) and immunohistochemistry analyses in 58 formalin-fixed, paraffin-embedded osteosarcoma archival tissues that had clinical follow-up information.

RESULTS

The pathway enrichment analyses of the aCGH data revealed that VEGF pathway genes, including the VEGFA gene itself, were amplified significantly in osteosarcoma. Genetic amplification of the VEGFA gene, both focally and in larger fragment, was validated by FISH analysis. It is noteworthy that amplification of the VEGFA gene and elevated expression of the VEGFA protein were associated significantly with microvascular density and adverse tumor-free survival in patients with osteosarcoma.

CONCLUSIONS

The authors report for the first time that VEGF pathway genes, including the VEGFA gene, are amplified in osteosarcoma. Amplification of the VEGFA gene is not only an important mechanism for elevated VEGFA protein expression but also is a poor prognostic factor for tumor-free survival. Combined classification of VEGFA gene amplification and positive VEGFA protein expression may provide a more accurate stratification method of selecting anti-VEGF therapy for patients with osteosarcoma.

Gastric cancers of Western European and African patients show different patterns of genomic instability

Background

Infection with H. pylori is important in the etiology of gastric cancer. Gastric cancer is infrequent in Africa, despite high frequencies of H. pylori infection, referred to as the African enigma. Variation in environmental and host factors influencing gastric cancer risk between different populations have been reported but little is known about the biological differences between gastric cancers from different geographic locations. We aim to study genomic instability patterns of gastric cancers obtained from patients from United Kingdom (UK) and South Africa (SA), in an attempt to support the African enigma hypothesis at the biological level.

Methods

DNA was isolated from 67 gastric adenocarcinomas, 33 UK patients, 9 Caucasian SA patients and 25 native SA patients. Microsatellite instability and chromosomal instability were analyzed by PCR and microarray comparative genomic hybridization, respectively. Data was analyzed by supervised univariate and multivariate analyses as well as unsupervised hierarchical cluster analysis.

Results

Tumors from Caucasian and native SA patients showed significantly more microsatellite instable tumors (p < 0.05). For the microsatellite stable tumors, geographical origin of the patients correlated with cluster membership, derived from unsupervised hierarchical cluster analysis (p = 0.001). Several chromosomal alterations showed significantly different frequencies in tumors from UK patients and native SA patients, but not between UK and Caucasian SA patients and between native and Caucasian SA patients.

Conclusions

Gastric cancers from SA and UK patients show differences in genetic instability patterns, indicating possible different biological mechanisms in patients from different geographical origin. This is of future clinical relevance for stratification of gastric cancer therapy.

Joint segmentation, calling, and normalization of multiple CGH profiles

The statistical analysis of array comparative genomic hybridization (CGH) data has now shifted to the joint assessment of copy number variations at the cohort level. Considering multiple profiles gives the opportunity to correct for systematic biases observed on single profiles, such as probe GC content or the so-called "wave effect." In this article, we extend the segmentation model developed in the univariate case to the joint analysis of multiple CGH profiles. Our contribution is multiple: we propose an integrated model to perform joint segmentation, normalization, and calling for multiple array CGH profiles. This model shows great flexibility, especially in the modeling of the wave effect that gives a likelihood framework to approaches proposed by others. We propose a new dynamic programming algorithm for break point positioning, as well as a model selection criterion based on a modified bayesian information criterion proposed in the univariate case. The performance of our method is assessed using simulated and real data sets. Our method is implemented in the R package cghseg.

Identification of recurrent regions of Copy-Number Variants across multiple individuals

BACKGROUND

Algorithms and software for CNV detection have been developed, but they detect the CNV regions sample-by-sample with individual-specific breakpoints, while common CNV regions are likely to occur at the same genomic locations across different individuals in a homogenous population. Current algorithms to detect common CNV regions do not account for the varying reliability of the individual CNVs, typically reported as confidence scores by SNP-based CNV detection algorithms. General methodologies for identifying these recurrent regions, especially those directed at SNP arrays, are still needed.

RESULTS

In this paper, we describe two new approaches for identifying common CNV regions based on (i) the frequency of occurrence of reliable CNVs, where reliability is determined by high confidence scores, and (ii) a weighted frequency of occurrence of CNVs, where the weights are determined by the confidence scores. In addition, motivated by the fact that we often observe partially overlapping CNV regions as a mixture of two or more distinct subregions, regions identified using the two approaches can be fine-tuned to smaller sub-regions using a clustering algorithm. We compared the performance of the methods with sequencing-based results in terms of discordance rates, rates of departure from Hardy-Weinberg equilibrium (HWE) and average frequency and size of the identified regions. The discordance rates as well as the rates of departure from HWE decrease when we select CNVs with higher confidence scores. We also performed comparisons with two previously published methods, STAC and GISTIC, and showed that the methods we consider are better at identifying low-frequency but high-confidence CNV regions.

CONCLUSIONS

The proposed methods for identifying common CNV regions in multiple individuals perform well compared to existing methods. The identified common regions can be used for downstream analyses such as group comparisons in association studies.

High-resolution array comparative genomic hybridization in sporadic and celiac disease-related small bowel adenocarcinomas

PURPOSE

The molecular pathogenesis of small intestinal adenocarcinomas is not well understood. Understanding the molecular characteristics of small bowel adenocarcinoma may lead to more effective patient treatment.

EXPERIMENTAL DESIGN

Forty-eight small bowel adenocarcinomas (33 non-celiac disease related and 15 celiac disease related) were characterized for chromosomal aberrations by high-resolution array comparative hybridization, microsatellite instability, and APC promoter methylation and mutation status. Findings were compared with clinicopathologic and survival data. Furthermore, molecular alterations were compared between celiac disease-related and non-celiac disease-related small bowel adenocarcinomas.

RESULTS

DNA copy number changes were observed in 77% small bowel adenocarcinomas. The most frequent DNA copy number changes found were gains on 5p15.33-5p12, 7p22.3-7q11.21, 7q21.2-7q21.3, 7q22.1-7q34, 7q36.1, 7q36.3, 8q11.21-8q24.3, 9q34.11-9q34.3, 13q11-13q34, 16p13.3, 16p11.2, 19q13.2, and 20p13-20q13.33, and losses on 4p13-4q35.2, 5q15-5q21.1, and 21p11.2-21q22.11. Seven highly amplified regions were identified on 6p21.1, 7q21.1, 8p23.1, 11p13, 16p11.2, 17q12-q21.1, and 19q13.2. Celiac disease-related and non-celiac disease-related small bowel adenocarcinomas displayed similar chromosomal aberrations. Promoter hypermethylation of the APC gene was found in 48% non-celiac disease-related and 73% celiac disease-related small bowel adenocarcinomas. No nonsense mutations were found. Thirty-three percent of non-celiac disease-related small bowel adenocarcinomas showed microsatellite instability, whereas 67% of celiac disease-related small bowel adenocarcinomas were microsatellite unstable.

CONCLUSIONS

Our study characterized chromosomal aberrations and amplifications involved in small bowel adenocarcinoma. At the chromosomal level, celiac disease-related and non-celiac disease-related small bowel adenocarcinomas did not differ. A defect in the mismatch repair pathways seems to be more common in celiac disease-related than in non-celiac disease-related small bowel adenocarcinomas. In contrast to colon and gastric cancers, no APC nonsense mutations were found in small bowel adenocarcinoma. However, APC promoter methylation seems to be a common event in celiac disease-related small bowel adenocarcinoma. Clin Cancer Res; 16(5); 1391-401.

Preprocessing and downstream analysis of microarray DNA copy number profiles

Analysis of DNA copy number profiles requires methods tailored to the specific nature of these data. The number of available data analysis methods has grown enormously in the last 5 years. We discuss the typical characteristics of DNA copy number data, as measured by microarray technology and review the extensive literature on preprocessing methods such as segmentation and calling. Subsequently, the focus narrows to applications of DNA copy number in cancer, in particular, several downstream analyses of multi-sample data sets such as testing, clustering and classification. Finally, we look ahead: what should we prepare for and which methodology-related topics may deserve attention in the near future?

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.

Accurate confidence aware clustering of array CGH tumor profiles

MOTIVATION

Chromosomal aberrations tend to be characteristic for given (sub)types of cancer. Such aberrations can be detected with array comparative genomic hybridization (aCGH). Clustering aCGH tumor profiles aids in identifying chromosomal regions of interest and provides useful diagnostic information on the cancer type. An important issue here is to what extent individual aCGH tumor profiles can be reliably assigned to clusters associated with a given cancer type.

RESULTS

We introduce a novel evolutionary fuzzy clustering (EFC) algorithm, which is able to deal with overlapping clusterings. Our method assesses these overlaps by using cluster membership degrees, which we use here as a confidence measure for individual samples to be assigned to a given tumor type. We first demonstrate the usefulness of our method using a synthetic aCGH dataset and subsequently show that EFC outperforms existing methods on four real datasets of aCGH tumor profiles involving four different cancer types. We also show that in general best performance is obtained using 1- Pearson correlation coefficient as a distance measure and that extra preprocessing steps, such as segmentation and calling, lead to decreased clustering performance.

AVAILABILITY

The source code of the program is available from http://ibi.vu.nl/programs/efcwww

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.

Copy number gain at 8q12.1-q22.1 is associated with a malignant tumor phenotype in salivary gland myoepitheliomas

Salivary gland myoepithelial tumors are relatively uncommon tumors with an unpredictable clinical course. More knowledge about their genetic profiles is necessary to identify novel predictors of disease. In this study, we subjected 27 primary tumors (15 myoepitheliomas and 12 myoepithelial carcinomas) to genome-wide microarray-based comparative genomic hybridization (array CGH). We set out to delineate known chromosomal aberrations in more detail and to unravel chromosomal differences between benign myoepitheliomas and myoepithelial carcinomas. Patterns of DNA copy number aberrations were analyzed by unsupervised hierarchical cluster analysis. Both benign and malignant tumors revealed a limited amount of chromosomal alterations (median of 5 and 7.5, respectively). In both tumor groups, high frequency gains (> or =20%) were found mainly at loci of growth factors and growth factor receptors (e.g., PDGF, FGF(R)s, and EGFR). In myoepitheliomas, high frequency losses (> or =20%) were detected at regions of proto-cadherins. Cluster analysis of the array CGH data identified three clusters. Differential copy numbers on chromosome arm 8q and chromosome 17 set the clusters apart. Cluster 1 contained a mixture of the two phenotypes (n = 10), cluster 2 included mostly benign tumors (n = 10), and cluster 3 only contained carcinomas (n = 7). Supervised analysis between malignant and benign tumors revealed a 36 Mbp-region at 8q being more frequently gained in malignant tumors (P = 0.007, FDR = 0.05). This is the first study investigating genomic differences between benign and malignant myoepithelial tumors of the salivary glands at a genomic level. Both unsupervised and supervised analysis of the genomic profiles revealed chromosome arm 8q to be involved in the malignant phenotype of salivary gland myoepitheliomas.

Genomic profiling identifies common HPV-associated chromosomal alterations in squamous cell carcinomas of cervix and head and neck

Background

It is well known that a persistent infection with high-risk human papillomavirus (hrHPV) is causally involved in the development of squamous cell carcinomas of the uterine cervix (CxSCCs) and a subset of SCCs of the head and neck (HNSCCs). The latter differ from hrHPV-negative HNSCCs at the clinical and molecular level.

Methods

To determine whether hrHPV-associated SCCs arising from different organs have specific chromosomal alterations in common, we compared genome-wide chromosomal profiles of 10 CxSCCs (all hrHPV-positive) with 12 hrHPV-positive HNSCCs and 30 hrHPV-negative HNSCCs. Potential organ-specific alterations and alterations shared by SCCs in general were investigated as well.

Results

Unsupervised hierarchical clustering resulted in one mainly hrHPV-positive and one mainly hrHPV-negative cluster. Interestingly, loss at 13q and gain at 20q were frequent in HPV-positive carcinomas of both origins, but uncommon in hrHPV-negative HNSCCs, indicating that these alterations are associated with hrHPV-mediated carcinogenesis. Within the group of hrHPV-positive carcinomas, HNSCCs more frequently showed gains of multiple regions at 8q whereas CxSCCs more often showed loss at 17p. Finally, gains at 3q24-29 and losses at 11q22.3-25 were frequent (>50%) in all sample groups.

Conclusion

In this study hrHPV-specific, organ-specific, and pan-SCC chromosomal alterations were identified. The existence of hrHPV-specific alterations in SCCs of different anatomical origin, suggests that these alterations are crucial for hrHPV-mediated carcinogenesis.

DNA copy number profiles of primary tumors as predictors of response to chemotherapy in advanced colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is biologically a heterogeneous disease, which may affect response to drug therapy. We investigated the correlation of genome-wide DNA copy number profiles of primary tumors with response to systemic chemotherapy in advanced CRC.

PATIENTS AND METHODS

DNA was isolated from formaldehyde-fixed paraffin-embedded primary tumors of 32 patients with advanced CRC, which were selected based on either a good response (n = 16) or a poor response (n = 16) to first-line combination therapy with capecitabine and irinotecan. High-resolution DNA copy number profiles were obtained by means of 30 K oligonucleotide-based array comparative genomic hybridization (aCGH).

RESULTS

Unsupervised hierarchical cluster analysis of the aCGH data revealed two clusters of 19 and 13 tumors, respectively, and cluster membership showed a significant correlation with response status (P < 0.03). The nonresponders had fewer chromosomal alterations compared with the responders, in particular less losses were found (P < 0.03). Most prominent differences between the two groups were losses of regions 18p11.32-q11.2 (P < 0.02) and 18q12.1-q23 (P < 0.03), which were more frequently observed in responders.

CONCLUSIONS

Differences in DNA copy number profiles of primary CRCs are associated with response to systemic combination chemotherapy with capecitabine and irinotecan. Responders overall had more chromosomal alterations, especially loss of chromosome 18.