DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly.

Detection of DNA copy number aberrations by shallow whole-genome sequencing (WGS) faces many challenges, including lack of completion and errors in the human reference genome, repetitive sequences, polymorphisms, variable sample quality, and biases in the sequencing procedures. Formalin-fixed paraffin-embedded (FFPE) archival material, the analysis of which is important for studies of cancer, presents particular analytical difficulties due to degradation of the DNA and frequent lack of matched reference samples. We present a robust, cost-effective WGS method for DNA copy number analysis that addresses these challenges more successfully than currently available procedures. In practice, very useful profiles can be obtained with ∼0.1× genome coverage. We improve on previous methods by first implementing a combined correction for sequence mappability and GC content, and second, by applying this procedure to sequence data from the 1000 Genomes Project in order to develop a blacklist of problematic genome regions. A small subset of these blacklisted regions was previously identified by ENCODE, but the vast majority are novel unappreciated problematic regions. Our procedures are implemented in a pipeline called QDNAseq. We have analyzed over 1000 samples, most of which were obtained from the fixed tissue archives of more than 25 institutions. We demonstrate that for most samples our sequencing and analysis procedures yield genome profiles with noise levels near the statistical limit imposed by read counting. The described procedures also provide better correction of artifacts introduced by low DNA quality than prior approaches and better copy number data than high-resolution microarrays at a substantially lower cost.

Driver Gene and Novel Mutations in Asbestos-Exposed Lung Adenocarcinoma and Malignant Mesothelioma Detected by Exome Sequencing.

BACKGROUND: Asbestos is a carcinogen linked to malignant mesothelioma (MM) and lung cancer. Some gene aberrations related to asbestos exposure are recognized, but many associated mutations remain obscure. We performed exome sequencing to determine the association of previously known mutations (driver gene mutations) with asbestos and to identify novel mutations related to asbestos exposure in lung adenocarcinoma (LAC) and MM. METHODS: Exome sequencing was performed on DNA from 47 tumor tissues of MM (21) and LAC (26) patients, 27 of whom had been asbestos-exposed (18 MM, 9 LAC). In addition, 9 normal lung/blood samples of LAC were sequenced. Novel mutations identified from exome data were validated by amplicon-based deep sequencing. Driver gene mutations in BRAF, EGFR, ERBB2, HRAS, KRAS, MET, NRAS, PIK3CA, STK11, and ephrin receptor genes (EPHA1-8, 10 and EPHB1-4, 6) were studied for both LAC and MM, and in BAP1, CUL1, CDKN2A, and NF2 for MM. RESULTS: In asbestos-exposed MM patients, previously non-described NF2 frameshift mutation (one) and BAP1 mutations (four) were detected. Exome data mining revealed some genes potentially associated with asbestos exposure, such as MRPL1 and SDK1. BAP1 and COPG1 mutations were seen exclusively in MM. Pathogenic KRAS mutations were common in LAC patients (42 %), both in non-exposed (n = 5) and exposed patients (n = 6). Pathogenic BRAF mutations were found in two LACs. CONCLUSION: BAP1 mutations occurred in asbestos-exposed MM. MRPL1, SDK1, SEMA5B, and INPP4A could possibly serve as candidate genes for alterations associated with asbestos exposure. KRAS mutations in LAC were not associated with asbestos exposure.

Landscape of chromosomal copy number aberrations in gangliogliomas and dysembryoplastic neuroepithelial tumours.

AIM: Gangliogliomas (GGs) and dysembryoplastic neuroepithelial tumours (DNTs) represent the most common histological entities within the spectrum of glioneuronal tumours (GNTs). The wide variability of morphological features complicates histological classification, including discrimination from prognostically distinct diffuse low-grade astrocytomas (AIIs). This study was performed to increase our understanding of these tumours. METHODS: We studied chromosomal copy number aberrations (CNAs) by genome-wide sequencing in a large cohort of GNTs and linked these to comprehensive histological analysis and clinical characteristics. One hundred fourteen GNTs were studied: 50 GGs and 64 DNTs. Also, a data set of CNAs from 38 diffuse AIIs was included. RESULTS: The most frequent CNAs in both GGs and DNTs were gains at chromosomes 5 and 7, often concurrent, and gain at chromosome 6. None of the CNAs was linked to histological subtype, immunohistochemical features or to clinical characteristics. Comparison of AIIs and diffuse GNTs revealed that gain at whole chromosome 5 is only observed in GNTs. CNA patterns indicative of chromothripsis were detected in three GNTs. CONCLUSION: We conclude that GNTs with diverse morphologies share molecular features, and our findings support the need to improve classification and differential diagnosis of tumour entities within the spectrum of GNTs, as well as their distinction from other gliomas.

Copy number alterations and neoplasia-specific mutations in MELK, PDCD1LG2, TLN1, and PAX5 at 9p in different neoplasias.

Genetic alterations affecting 9p are commonly present in many cancer types and many cancer-related genes are located in this chromosomal region. We sequenced all of the genes located in a 32Mb region of 9p by targeted next generation sequencing (NGS) in 96 patients with different cancer types, including acute lymphoblastic leukemia, bone malignant fibrous histiocytoma/undifferentiated pleomorphic sarcoma, fibrosarcoma, Ewing's sarcoma, and lung carcinoma. Copy number alterations (CNA), and mutations were studied from the NGS data. We detected a deletion at the CDKN2A locus as being the most frequent genetic alteration in all cancer types. In addition to this locus, NGS also identified other small regions of copy number loss and gain. However, different cancer types did not reveal any statistically significant differences with regard to CNA frequency or type. Of the 191 genes within the target region, two novel recurrent mutations were found in the MELK and PDCD1LG2 genes. The most commonly mutated gene in sarcomas was TLN1 (8%) and PAX5 in ALL (9%). Mutations in PAX5, and RUSC2, were seen exclusively in ALL patients and those in KIAA1432, CA9, TLN1, and MELK only in sarcomas (MFH, FS, EFT). Thus using targeted NGS of the 9p region, in addition to commonly deleted CDKN2A locus, we were able to identify a number of small deletions and gains, as well as novel recurrent mutations in different cancer types. © 2014 Wiley Periodicals, Inc.

Targeted resequencing of 9p in acute lymphoblastic leukemia yields concordant results with array CGH and reveals novel genomic alterations.

Genetic alterations of the short arm of chromosome 9 are frequent in acute lymphoblastic leukemia. We performed targeted sequencing of 9p region in 35 adolescent and adult acute lymphoblastic leukemia patients and sought to investigate the sensitivity of detecting copy number alterations in comparison with array comparative genomic hybridization (aCGH), and besides, to detect novel genetic anomalies. We found a high concordance of copy number variations (CNVs) as detected by next generation sequencing (NGS) and aCGH. By both methodologies, the recurrent deletion at CDKN2A/B locus was identified, whereas NGS revealed additional, small regions of CNVs, seen more frequently in adult patients, while aCGH was better at detecting larger CNVs. Also, by NGS, we detected novel structural variations, novel SNVs and small insertion/deletion variants. Our results show that NGS, in addition to detecting mutations and other genetic aberrations, can be used to study CNVs.

High-resolution SNP array analysis of patients with developmental disorder and normal array CGH results.

BACKGROUND: Diagnostic analysis of patients with developmental disorders has improved over recent years largely due to the use of microarray technology. Array methods that facilitate copy number analysis have enabled the diagnosis of up to 20% more patients with previously normal karyotyping results. A substantial number of patients remain undiagnosed, however. METHODS AND RESULTS: Using the Genome-Wide Human SNP array 6.0, we analyzed 35 patients with a developmental disorder of unknown cause and normal array comparative genomic hybridization (array CGH) results, in order to characterize previously undefined genomic aberrations. We detected no seemingly pathogenic copy number aberrations. Most of the vast amount of data produced by the array was polymorphic and non-informative. Filtering of this data, based on copy number variant (CNV) population frequencies as well as phenotypically relevant genes, enabled pinpointing regions of allelic homozygosity that included candidate genes correlating to the phenotypic features in four patients, but results could not be confirmed. CONCLUSIONS: In this study, the use of an ultra high-resolution SNP array did not contribute to further diagnose patients with developmental disorders of unknown cause. The statistical power of these results is limited by the small size of the patient cohort, and interpretation of these negative results can only be applied to the patients studied here. We present the results of our study and the recurrence of clustered allelic homozygosity present in this material, as detected by the SNP 6.0 array.

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.

CGHpower: exploring sample size calculations for chromosomal copy number experiments.

BACKGROUND: Determining a suitable sample size is an important step in the planning of microarray experiments. Increasing the number of arrays gives more statistical power, but adds to the total cost of the experiment. Several approaches for sample size determination have been developed for expression array studies, but so far none has been proposed for array comparative genomic hybridization (aCGH). RESULTS: Here we explore power calculations for aCGH experiments comparing two groups. In a pilot experiment CGHpower estimates the biological diversity between groups and provides a statistical framework for estimating average power as a function of sample size. As the method requires pilot data, it can be used either in the planning stage of larger studies or in estimating the power achieved in past experiments. CONCLUSIONS: The proposed method relies on certain assumptions. According to our evaluation with public and simulated data sets, they do not always hold true. Violation of the assumptions typically leads to unreliable sample size estimates. Despite its limitations, this method is, at least to our knowledge, the only one currently available for performing sample size calculations in the context of aCGH. Moreover, the implementation of the method provides diagnostic plots that allow critical assessment of the assumptions on which it is based and hence on the feasibility and reliability of the sample size calculations in each case.The CGHpower web application and the program outputs from evaluation data sets can be freely accessed at http://www.cangem.org/cghpower/

CanGEM: mining gene copy number changes in cancer.

The use of genome-wide and high-throughput screening methods on large sample sizes is a well-grounded approach when studying a process as complex and heterogeneous as tumorigenesis. Gene copy number changes are one of the main mechanisms causing cancerous alterations in gene expression and can be detected using array comparative genomic hybridization (aCGH). Microarrays are well suited for the integrative systems biology approach, but none of the existing microarray databases is focusing on copy number changes. We present here CanGEM (Cancer GEnome Mine), which is a public, web-based database for storing quantitative microarray data and relevant metadata about the measurements and samples. CanGEM supports the MIAME standard and in addition, stores clinical information using standardized controlled vocabularies whenever possible. Microarray probes are re-annotated with their physical coordinates in the human genome and aCGH data is analyzed to yield gene-specific copy numbers. Users can build custom datasets by querying for specific clinical sample characteristics or copy number changes of individual genes. Aberration frequencies can be calculated for these datasets, and the data can be visualized on the human genome map with gene annotations. Furthermore, the original data files are available for more detailed analysis. The CanGEM database can be accessed at http://www.cangem.org/.

Integrated gene copy number and expression microarray analysis of gastric cancer highlights potential target genes.

We performed an integrated array comparative genomic hybridization (aCGH) and expression microarray analysis of 8 normal gastric tissues and 38 primary tumors, including 25 intestinal and 13 diffuse gastric adenocarcinomas to identify genes whose expression is deregulated in association with copy number alteration. Our aim was also to identify molecular genetic alterations that are specific to particular clinicopathological characteristics of gastric cancer. Distinct molecular genetic profiles were identified for intestinal and diffuse gastric cancers and for tumors obtained from 2 different locations of the stomach. Interestingly, the ERBB2 amplification and gains at 20q13.12-q13.33 almost exclusively discriminated intestinal cancers from the diffuse type. In addition, the 17q12-q25 gain was characteristic to cancers located in corpus and the 20q13.12-q13.13 gain was more common in the antrum. Statistical analysis was performed using integrated copy number and expression data to identify genes showing differential expression associated with a copy number alteration. Genes with the highest statistical significance included ERBB2, MUC1, GRB7, PPP1R1B and PPARBP with concomitant changes in copy number and expression. Immunohistochemical analysis of ERBB2 and MUC1 on a tissue microarray containing 78 independent gastric tissues showed statistically significant differences (p < 0.05 and <0.001) in immunopositivity in the intestinal (31 and 70%) and diffuse subtypes (14 and 41%), respectively. In conclusion, our results demonstrate that intestinal and diffuse type gastric cancers as well as cancers located in different sites of the stomach have distinct molecular profiles which may have clinical value.

Clinical value of chromosome arms 19q and 11p losses in low-grade gliomas.

BACKGROUND: Diffuse low-grade gliomas (LGGs) form a heterogeneous subgroup of gliomas in adults. Chromosome (chr) arms 1p/19q codeletion and IDH mutation have been shown to be closely associated with oligodendroglial phenotype and better prognosis. We sought to identify relevant biomarkers in non 1p/19q codeleted LGGs. METHODS: We characterized a retrospective series of 126 LGGs using genomic arrays, microsatellite analysis, IDH sequencing, MGMT promoter methylation assay, and p53 expression analysis. RESULTS: Our study confirms that 1p/19q codeletion, mutually exclusive with p53 overexpression, was associated with: (i) better prognosis, (ii) oligodendroglial phenotype, (iii) MGMT promoter methylation, and (iv) IDH mutation. Interestingly, 1p/19q codeleted tumors occur in older patients at diagnosis. Our study shows that non 1p/19q codeleted LGGs can be divided in 5 main genomic subgroups: (i) 11p loss, (ii) 19q loss (iii) 7 gain, (iv) 19 gain, and (v) unclassified. In non 1p/19q codeleted LGGs, we demonstrated that (i) 11p loss is associated with astrocytoma phenotype and has an independent negative prognostic value, and (ii) 19q loss diminished the favorable prognostic value of IDH mutation. Our findings were validated in an independent cohort of 98 LGGs. CONCLUSION: Novel genomic entities and biomarkers have been identified in non 1p/19q codeleted LGGs. Our findings may help to stratify non 1p/19q codeleted LGGs, facilitating future individualization of treatment. Further prospective studies are warranted to support our findings.

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.

Homozygous deletions of cadherin genes in chondrosarcoma-an array comparative genomic hybridization study.

Chondrosarcoma is a malignant bone tumor that is often resistant to chemotherapy and radiotherapy. We applied high resolution oligonucleotide array comparative genomic hybridization to 46 tumor specimens from 44 patients with chondrosarcoma and identified several genes with potential importance for the development of chondrosarcoma. Several homozygous deletions were detected. The tumor suppressor genes CDKN2A and MTAP were each homozygously deleted in four of the cases, and the RB1 gene was homozygously deleted in one. Two homozygous deletions of MTAP did not affect CDKN2A. Deletions were also found to affect genes of the cadherin family, including CDH4 and CDH7, each of which had a targeted homozygous loss in one case, and CDH19, which had a targeted homozygous loss in two cases. Loss of the EXT1 and EXT2 genes was uncommon; EXT1 was homozygously deleted in none and EXT2 in two of the cases, and large heterozygous losses including EXT1 and/or EXT2 were seen in three cases. Targeted gains and amplifications affected the MYC, E2F3, CDK6, PDGFRA, KIT, and PDGFD genes in one case each. The data indicate that chondrosarcomas develop through a combination of genomic imbalances that often affect the RB1 signaling pathway. The inactivation of cadherin genes may also be critical in the pathogenesis of the tumor.

miRNA expression profiles in myelodysplastic syndromes reveal Epstein-Barr virus miR-BART13 dysregulation.

Recently, the microRNA (miRNA) signature has been used for better characterization and understanding of the pathogenesis of different malignancies, including myelodysplastic syndromes (MDS). MDS are a heterogeneous group of stem cell disorders in which the genetic and molecular defects are not well defined. In the present study, we applied array based miRNA profiling to study 19 bone marrow cell samples of de novo MDS compared with eight healthy individuals. In addition, integration of the miRNA profiling data with our previous array comparative genomic hybridization data, from the same cohort of patients, was performed. We observed up-regulation of hsa-miR-720 and hsa-miR-21, and down-regulation of hsa-miR-671-5p and one human virus miRNA (Epstein-Barr virus miR-BART13) in MDS samples compared with normal samples. In our study, the copy number alteration harboring miRNA was not affecting miRNA expression, but a distinct microRNA expression pattern was observed, not only in MDS compared with controls, but also between MDS entities.

New Insights into the Cellular Pathways Affected in Primary Uterine Leiomyosarcoma.

Resistance to chemotherapeutic agents and radiotherapy has kept surgery the primary treatment of uterine leiomyosarcoma (ULMS). In search of leads for potential therapeutic targets, array CGH (aCGH) was used to obtain a genomewide pattern of ULMS-specific genetic imbalances and to define the affected biological processes. Fine-resolution genomewide aCGH analysis was performed using customised 16K cDNA microarrays on 18 primary ULMS cases. Furthermore, patterns of DNA copy number changes were assessed for associations with clinical parameters, i.e., tumour grade, tumour size and patient status at last follow-up. Our aCGH results demonstrated extensive DNA copy number changes in all chromosomes. Of the 10,590 gene loci included in the analysis, 4,387 were found to be affected by DNA copy number gains and 4,518 by DNA copy number losses in at least one case. Further analyses revealed that 231 of these were commonly gained, and 265 lost in at least 20% of the cases. The gains affected loci at 1p, 1q, 2p, 3p, 6p, 8q, 10q and 18q, whereas losses were observed at 2q, 4q, 6p, 6q, 7p, 7q, 13q, 14p, 16q, 19p, Xp and Xq. Enrichment analysis of biological processes revealed the gained genes to be involved in the G1/S transition of mitotic cell cycle, co-translational protein targeting to membrane, actin filament polymerisation and positive regulation of cytokine biosynthesis, whereas the genes affected by losses were associated with DNA replication, chromatin modification, telomere maintenance, meiosis, mitosis and angiogenesis. These biological processes featured prominently two well-established tumour suppressors (BRCA2, EREG) and one proto-oncogene (GFI1). No statistically significant associations were found between the aberration patterns and clinical variables. Analysis of gene pathways using aCGH uncovered the biological networks involved in malignant progression of ULMS.