M-CGH: Analysing microarray-based CGH experiments

BACKGROUND: Microarray-based comparative genomic hybridisation (array CGH) is a technique by which variation in relative copy numbers between two genomes can be analysed by competitive hybridisation to DNA microarrays. This technology has most commonly been used to detect chromosomal amplifications and deletions in cancer. Dedicated tools are needed to analyse the results of such experiments, which include appropriate visualisation, and to take into consideration the physical relation in the genome between the probes on the array. RESULTS: M-CGH is a MATLAB toolbox with a graphical user interface designed specifically for the analysis of array CGH experiments, with multiple approaches to ratio normalization. Specifically, the distributions of three classes of DNA copy numbers (gains, normal and losses) can be estimated using a maximum likelihood method. Amplicon boundaries are computed by either the fuzzy K-nearest neighbour method or a wavelet approach. The program also allows linking each genomic clone with the corresponding genomic information in the Ensembl database . CONCLUSIONS: M-CGH, which encompasses the basic tools needed for analysing array CGH experiments, is freely available for academics , and does not require any other MATLAB toolbox

DNA copy number changes in high-grade malignant peripheral nerve sheath tumors by array CGH

<p>Abstract</p> <p>Background</p> <p>Malignant peripheral nerve sheath tumors (MPNSTs) are rare and highly aggressive soft tissue tumors showing complex chromosomal aberrations. In order to identify recurrent chromosomal regions of gain and loss, and thereby novel gene targets of potential importance for MPNST development and/or progression, we have analyzed DNA copy number changes in seven high-grade MPNSTs using microarray-based comparative genomic hybridization (array CGH).</p> <p>Results</p> <p>Considerable more gains than losses were observed, and the most frequent minimal recurrent regions of gain included 1q24.1-q24.2, 1q24.3-q25.1, 8p23.1-p12, 9q34.11-q34.13 and 17q23.2-q25.3, all gained in five of seven samples. The 17q23.2-q25.3 region was gained in all five patients with poor outcome and not in the two patients with disease-free survival. cDNA microarray analysis and quantitative real-time reverse transcription PCR were used to investigate expression of genes located within these regions. The gene lysyl oxidase-like 2 (<it>LOXL2</it>) was identified as a candidate target for the 8p23.1-p12 gain. Within 17q, the genes topoisomerase II-α (<it>TOP2A</it>), ets variant gene 4 (E1A enhancer binding protein, <it>E1AF</it>) (<it>ETV4</it>) and baculoviral IAP repeat-containing 5 (survivin) (<it>BIRC5</it>) showed increased expression in all samples compared to two benign tumors. Increased expression of these genes has previously been associated with poor survival in other malignancies, and for <it>TOP2A</it>, in MPNSTs as well. In addition, we have analyzed the expression of five micro RNAs located within the 17q23.2-q25.3 region, but none of them showed high expression levels compared to the benign tumors.</p> <p>Conclusion</p> <p>Our study shows the potential of using DNA copy number changes obtained by array CGH to predict the prognosis of MPNST patients. Although no clear correlations between the expression level and patient outcome were observed, the genes <it>TOP2A</it>, <it>ETV4 </it>and <it>BIRC5 </it>are interesting candidate targets for the 17q gain associated with poor survival.</p

Clinical and molecular implications of NAB2-STAT6 fusion variants in solitary fibrous tumour

Solitary fibrous tumour (SFT) is a mesenchymal neoplasm characterised by pathognomonic NAB2-STAT6 gene fusions. The clinical implications and prognostic value of different fusion variants has not been clarified. In the current study, we explore the clinicopathological, prognostic and molecular differences between tumours with different fusions. Thirty-nine patients with localised, extrameningeal SFT were included, of whom 20 developed distant recurrence and 19 were without recurrence after long term follow-up. Capture-based RNA sequencing identified 12 breakpoint variants, which were categorised into two groups based on the STAT6 domain composition in the predicted chimeric proteins. Twenty-one of 34 (62%) sequenced tumours had fusions with most of the STAT6 domains intact and were classified as STAT6-Full. Thirteen tumours (38%) contained only the transactivation domain of STAT6 and were classified as STAT6-TAD. Tumours with STAT6-TAD fusions had a higher mitotic count (p=0.016) and were associated with inferior recurrence-free interval (p=0.004) and overall survival (p=0.012). Estimated 10-year recurrence-free survival was 25% for patients with STAT6-TAD tumours compared to 78% for the STAT6-Full group. Distinct transcriptional signatures between the fusion groups were identified, including higher expression of FGF2 in the STAT6-TAD group and IGF2, EGR2, PDGFRB, STAT6 and several extracellular matrix genes in STAT6-Full tumours. In summary, we demonstrate that NAB2-STAT6 fusion variants are associated with distinct clinicopathological and molecular characteristics and have prognostic significance in extrameningeal SFT.publishedVersio

Evaluation of high-resolution microarray platforms for genomic profiling of bone tumours

<p>Abstract</p> <p>Background</p> <p>Several high-density oligonucleotide microarray platforms are available for genome-wide single nucleotide polymorphism (SNP) detection and microarray-based comparative genomic hybridisation (array CGH), which may be used to detect copy number aberrations in human tumours. As part of the EuroBoNeT network of excellence for research on bone tumours (eurobonet.eu), we have evaluated four different commercial high-resolution microarray platforms in order to identify the most appropriate technology for mapping DNA copy number aberrations in such tumours.</p> <p>Findings</p> <p>DNA from two different cytogenetically well-characterized bone sarcoma cell lines, representing a simple and a complex karyotype, respectively, was tested in duplicate on four high-resolution microarray platforms; Affymetrix Genome-Wide Human SNP Array 6.0, Agilent Human Genome CGH 244A, Illumina HumanExon510s-duo and Nimblegen HG18 CGH 385 k WG tiling v1.0. The data was analysed using the platform-specific analysis software, as well as a platform-independent analysis algorithm. DNA copy number was measured at six specific chromosomes or chromosomal regions, and compared with the expected ratio based on available cytogenetic information. All platforms performed well in terms of reproducibility and were able to delimit and score small amplifications and deletions at similar resolution, but Agilent microarrays showed better linearity and dynamic range. The platform-specific analysis software provided with each platform identified in general correct copy numbers, whereas using a platform-independent analysis algorithm, correct copy numbers were determined mainly for Agilent and Affymetrix microarrays.</p> <p>Conclusions</p> <p>All platforms performed reasonably well, but Agilent microarrays showed better dynamic range, and like Affymetrix microarrays performed well with the platform-independent analysis software, implying more robust data. Bone tumours like osteosarcomas are heterogeneous tumours with complex karyotypes that may be difficult to interpret, and it is of importance to be able to well separate the copy number levels and detect copy number changes in subpopulations. Taking all this into consideration, the Agilent and Affymetrix microarray platforms were found to be a better choice for mapping DNA copy numbers in bone tumours, the latter having the advantage of also providing heterozygosity information.</p

Optimality of broken extremals

In this paper we analyse the optimality of broken Pontryagin extremal for an n-dimensional affine control system with a control parameter, taking values in a k- dimensional closed ball. We prove the optimality of broken normal extremals when n = 3 and the controllable vector fields form a contact distribution, and when the Lie algebra of the controllable fields is locally orthogonal to the singular locus and the drift does not belong to it. Moreover, if k = 2, we show the optimality of any broken extremal even abnormal when the controllable fields do not form a contact distribution in the point of singularity.Comment: arXiv admin note: text overlap with arXiv:1610.0675

The interplay between neoantigens and immune cells in sarcomas treated with checkpoint inhibition

IntroductionSarcomas are comprised of diverse bone and connective tissue tumors with few effective therapeutic options for locally advanced unresectable and/or metastatic disease. Recent advances in immunotherapy, in particular immune checkpoint inhibition (ICI), have shown promising outcomes in several cancer indications. Unfortunately, ICI therapy has provided only modest clinical responses and seems moderately effective in a subset of the diverse subtypes.MethodsTo explore the immune parameters governing ICI therapy resistance or immune escape, we performed whole exome sequencing (WES) on tumors and their matched normal blood, in addition to RNA-seq from tumors of 31 sarcoma patients treated with pembrolizumab. We used advanced computational methods to investigate key immune properties, such as neoantigens and immune cell composition in the tumor microenvironment (TME).ResultsA multifactorial analysis suggested that expression of high quality neoantigens in the context of specific immune cells in the TME are key prognostic markers of progression-free survival (PFS). The presence of several types of immune cells, including T cells, B cells and macrophages, in the TME were associated with improved PFS. Importantly, we also found the presence of both CD8+ T cells and neoantigens together was associated with improved survival compared to the presence of CD8+ T cells or neoantigens alone. Interestingly, this trend was not identified with the combined presence of CD8+ T cells and TMB; suggesting that a combined CD8+ T cell and neoantigen effect on PFS was important.DiscussionThe outcome of this study may inform future trials that may lead to improved outcomes for sarcoma patients treated with ICI

Mapping and characterization of the amplicon near APOA2 in 1q23 in human sarcomas by FISH and array CGH

BACKGROUND: Amplification of the q21-q23 region on chromosome 1 is frequently found in sarcomas and a variety of other solid tumours. Previous analyses of sarcomas have indicated the presence of at least two separate amplicons within this region, one located in 1q21 and one located near the apolipoprotein A-II (APOA2) gene in 1q23. In this study we have mapped and characterized the amplicon in 1q23 in more detail. RESULTS: We have used fluorescence in situ hybridisation (FISH) and microarray-based comparative genomic hybridisation (array CGH) to map and define the borders of the amplicon in 10 sarcomas. A subregion of approximately 800 kb was identified as the core of the amplicon. The amplification patterns of nine possible candidate target genes located to this subregion were determined by Southern blot analysis. The genes activating transcription factor 6 (ATF6) and dual specificity phosphatase 12 (DUSP12) showed the highest level of amplification, and they were also shown to be over-expressed by quantitative real-time reverse transcription PCR (RT-PCR). In general, the level of expression reflected the level of amplification in the different tumours. DUSP12 was expressed significantly higher than ATF6 in a subset of the tumours. In addition, two genes known to be transcriptionally activated by ATF6, glucose-regulated protein 78 kDa and -94 kDa (GRP78 and GRP94), were shown to be over-expressed in the tumours that showed over-expression of ATF6. CONCLUSION: ATF6 and DUSP12 seem to be the most likely candidate target genes for the 1q23 amplification in sarcomas. Both genes have possible roles in promoting cell growth, which makes them interesting candidate targets

The Interplay Between Neoantigens and Immune Cells in Sarcomas Treated With Checkpoint Inhibition

INTRODUCTION: Sarcomas are comprised of diverse bone and connective tissue tumors with few effective therapeutic options for locally advanced unresectable and/or metastatic disease. Recent advances in immunotherapy, in particular immune checkpoint inhibition (ICI), have shown promising outcomes in several cancer indications. Unfortunately, ICI therapy has provided only modest clinical responses and seems moderately effective in a subset of the diverse subtypes. METHODS: To explore the immune parameters governing ICI therapy resistance or immune escape, we performed whole exome sequencing (WES) on tumors and their matched normal blood, in addition to RNA-seq from tumors of 31 sarcoma patients treated with pembrolizumab. We used advanced computational methods to investigate key immune properties, such as neoantigens and immune cell composition in the tumor microenvironment (TME). RESULTS: A multifactorial analysis suggested that expression of high quality neoantigens in the context of specific immune cells in the TME are key prognostic markers of progression-free survival (PFS). The presence of several types of immune cells, including T cells, B cells and macrophages, in the TME were associated with improved PFS. Importantly, we also found the presence of both CD8+ T cells and neoantigens together was associated with improved survival compared to the presence of CD8+ T cells or neoantigens alone. Interestingly, this trend was not identified with the combined presence of CD8+ T cells and TMB; suggesting that a combined CD8+ T cell and neoantigen effect on PFS was important. DISCUSSION: The outcome of this study may inform future trials that may lead to improved outcomes for sarcoma patients treated with ICI

GeneCount: genome-wide calculation of absolute tumor DNA copy numbers from array comparative genomic hybridization data

Absolute tumor DNA copy numbers can currently be achieved only on a single gene basis by using fluorescence in situ hybridization (FISH). We present GeneCount, a method for genome-wide calculation of absolute copy numbers from clinical array comparative genomic hybridization data. The tumor cell fraction is reliably estimated in the model. Data consistent with FISH results are achieved. We demonstrate significant improvements over existing methods for exploring gene dosages and intratumor copy number heterogeneity in cancers

Genome-wide Map of Quantified Epigenetic Changes during In vitro Chondrogenic Differentiation of Primary Human Mesenchymal Stem Cells

Background: For safe clinical application of engineered cartilage made from mesenchymal stem cells (MSCs), molecular mechanisms for chondrogenic differentiation must be known in detail. Changes in gene expression and extracellular matrix synthesis have been extensively studied, but the epigenomic modifications underlying these changes have not been described. To this end we performed whole-genome chromatin immunoprecipitation and deep sequencing to quantify six histone modifications, reduced representation bisulphite sequencing to quantify DNA methylation and mRNA microarrays to quantify gene expression before and after 7 days of chondrogenic differentiation of MSCs in an alginate scaffold. To add to the clinical relevance of our observations, the study is based on primary bone marrow-derived MSCs from four donors, allowing us to investigate inter-individual variations. Results: We see two levels of relationship between epigenetic marking and gene expression. First, a large number of genes ontogenetically linked to MSC properties and the musculoskeletal system are epigenetically prepatterned by moderate changes in H3K4me3 and H3K9ac near transcription start sites. Most of these genes remain transcriptionally unaltered. Second, transcriptionally upregulated genes, more closely associated with chondrogenesis, are marked by H3K36me3 in gene bodies, highly increased H3K4me3 and H3K9ac on promoters and 5' end of genes, and increased H3K27ac and H3K4me1 marking in at least one enhancer region per upregulated gene. Within the 7-day time frame, changes in promoter DNA methylation do not correlate significantly with changes in gene expression. Inter-donor variability analysis shows high level of similarity between the donors for this data set. Conclusions: Histone modifications, rather than DNA methylation, provide the primary epigenetic control of early differentiation of MSCs towards the chondrogenic lineage.Stem Cell and Regenerative Biolog