Array Comparative Genomic Hybridization in Sarcomas

Over the past years, much research on sarcomas based on low-resolution cytogenetic and molecular cytogenetic methods has been published, leading to the identification of genetic abnormalities partially underlying the tumourigenesis. Continued progress in the identification of genetic events such as copy number aberrations relies upon adapting the rapidly evolving high-resolution microarray technology, which will eventually provide novel insights into sarcoma biology, and targets for both diagnostics and drug development. The aim of this Thesis was to characterize DNA copy number changes that are involved in the pathogenesis of soft tissue leiomyosarcoma (LMS), dermatofibrosarcoma protuberans (DFSP), osteosarcoma (OS), malignant fibrous histiocytoma (MFH), and uterine leiomyosarcoma (ULMS) by applying fine resolution array comparative genomic hybridization (aCGH) technology. Both low- and high-grade LMS tumours showed distinct copy number patterns, in addition to sharing two minimal common regions of gains and losses. Small aberrations were detected by aCGH, which were beyond the resolution of chromosomal comparative genomic hybridization (cCGH). DFSP tumours analysed by aCGH showed gains in 17q, 22q, and 21 additional gained regions, but only one region (22q) with copy number loss. Recurrent amplicons identified in OS by aCGH were 12q11-q15, 8q, 6p12-p21, and 17p. Amplicons 12q and 17p were further characterized in detail. The amplicon at 17p was characterized by aCGH in low- and high-grade LMS, OS, and MFH. In all but one case this amplicon, with minimal common regions of gains at 17p11-p12, started with the distal loss of 17p13-pter. OS and high-grade LMS were grouped together as they showed a complex pattern of copy number gains and amplifications at 17p, whereas MFH and low-grade LMS showed a continuous pattern of copy number gains and amplification at 17p. In addition to the commonly gained and lost regions identified in ULMS by aCGH, various biological processes affected by these copy number changes were also indicated by pathway analysis. The three novel findings obtained in this work were: characterization of amplicon 17p in low- and high-grade LMS and MFH, profiles of DNA copy number changes in LMS, and detection of various pathways affected by copy number changes in ULMS. These studies have not been undertaken previously by aCGH technology, thus this Thesis adds new information regarding DNA copy number changes in sarcomas. In conclusion, the aCGH technique used in this Thesis has provided new insights into the genetics of sarcomas by detecting the precise regions affected by copy number changes and some potential candidate target genes within those regions, which had not been uncovered by previously applied low resolution techniques.

Functional genes and gene array analysis as tools for monitoring hydrocarbon biodegradation

Bioremediation, which is the exploitation of the intrinsic ability of environmental microbes to degrade and remove harmful compounds from nature, is considered to be an environmentally sustainable and cost-effective means for environmental clean-up. However, a comprehensive understanding of the biodegradation potential of microbial communities and their response to decontamination measures is required for the effective management of bioremediation processes. In this thesis, the potential to use hydrocarbon-degradative genes as indicators of aerobic hydrocarbon biodegradation was investigated. Small-scale functional gene macro- and microarrays targeting aliphatic, monoaromatic and low molecular weight polyaromatic hydrocarbon biodegradation were developed in order to simultaneously monitor the biodegradation of mixtures of hydrocarbons. The validity of the array analysis in monitoring hydrocarbon biodegradation was evaluated in microcosm studies and field-scale bioremediation processes by comparing the hybridization signal intensities to hydrocarbon mineralization, real-time polymerase chain reaction (PCR), dot blot hybridization and both chemical and microbiological monitoring data. The results obtained by real-time PCR, dot blot hybridization and gene array analysis were in good agreement with hydrocarbon biodegradation in laboratory-scale microcosms. Mineralization of several hydrocarbons could be monitored simultaneously using gene array analysis. In the field-scale bioremediation processes, the detection and enumeration of hydrocarbon-degradative genes provided important additional information for process optimization and design. In creosote-contaminated groundwater, gene array analysis demonstrated that the aerobic biodegradation potential that was present at the site, but restrained under the oxygen-limited conditions, could be successfully stimulated with aeration and nutrient infiltration. During ex situ bioremediation of diesel oil- and lubrication oil-contaminated soil, the functional gene array analysis revealed inefficient hydrocarbon biodegradation, caused by poor aeration during composting. The functional gene array specifically detected upper and lower biodegradation pathways required for complete mineralization of hydrocarbons. Bacteria representing 1 % of the microbial community could be detected without prior PCR amplification. Molecular biological monitoring methods based on functional genes provide powerful tools for the development of more efficient remediation processes. The parallel detection of several functional genes using functional gene array analysis is an especially promising tool for monitoring the biodegradation of mixtures of hydrocarbons.

Array-based gene expression, CGH and tissue data defines a 12q24 gain in neuroblastic tumors with prognostic implication

Neuroblastoma has successfully served as a model system for the identification of neuroectoderm-derived oncogenes. However, in spite of various efforts, only a few clinically useful prognostic markers have been found. Here, we present a framework, which integrates DNA, RNA and tissue data to identify and prioritize genetic events that represent clinically relevant new therapeutic targets and prognostic biomarkers for neuroblastoma.

Molecular Genetic Background of Tumours in Hereditary Leiomyomatosis and Renal Cell Cancer Syndrome

Hereditary Leiomyomatosis and Renal Cell Cancer (HLRCC) is a hereditary tumour predisposition syndrome. Its phenotype includes benign cutaneous and uterine leiomyomas (CLM, ULM) with high penetrance and rarer renal cell cancer (RCC), most commonly of papillary type 2 subtype. Over 130 HLRCC families have been identified world-wide but the RCC phenotype seems to concentrate in families from Finland and North America for unknown reasons. HLRCC is caused by heterozygous germline mutations in the fumarate hydratase (FH) gene. FH encodes the enzyme fumarase from mitochondrial citric acid cycle. Fumarase enzyme activity or type or site of the FH mutation are unassociated with disease phenotype. The strongest evidence for tumourigenesis mechanism in HLRCC supports a hypoxia inducible factor driven process called pseudohypoxia resulting from accumulation of the fumarase substrate fumarate. In this study, to assess the importance of gene- or exon-level deletions or amplifications of FH in patients with HLRCC-associated phenotypes, multiplex ligation-dependent probe amplification (MLPA) method was used. One novel FH mutation, deletion of exon 1, was found in a Swedish male patient with an evident HLRCC phenotype with CLM, RCC, and a family history of ULM and RCC. Six other patients with CLM and 12 patients with only RCC or uterine leiomyosarcoma (ULMS) remained FH mutation-negative. These results suggest that copy number aberrations of FH or its exons are an infrequent cause of HLRCC and that only co-occurrence of benign tumour types justifies FH-mutation screening in RCC or ULMS patients. Determination of the genomic profile of 11 HLRCC-associated RCCs from Finnish patients was performed by array comparative genomic hybridization. The most common copy number aberrations were gains of 2, 7, and 17 and losses of 13q12.3-q21.1, 14, 18, and X. When compared to aberrations of sporadic papillary RCCs, HLRCC-associated RCCs harboured a distinct DNA copy number profile and lacked many of the changes characterizing the sporadic RCCs. The findings suggest a divergent molecular pathway for tumourigenesis of papillary RCCs in HLRCC. In order to find a genetic modifier of RCC risk in HLRCC, genome-wide linkage and identical by descent (IBD) analysis studies were performed in Finnish HLRCC families with microsatellite marker mapping and SNP-array platforms. The linkage analysis identified only one locus of interest, the FH gene locus in 1q43, but no mutations were found in the genes of the region. IBD analysis yielded no convincing haplotypes shared by RCC patients. Although these results do not exclude the existence of a genetic modifier for RCC risk in HLRCC, they emphasize the role of FH mutations in the malignant tumourigenesis of HLRCC. To study the benign tumours in HLRCC, genome-wide DNA copy number and gene expression profiles of sporadic and HLRCC ULMs were defined with modern SNP- and gene-expression array platforms. The gene expression array suggests novel genes involved in FH-deficient ULM tumourigenesis and novel genes with putative roles in propagation of sporadic ULM. Both the gene expression and copy number profiles of HLRCC ULMs differed from those of sporadic ULMs indicating distinct molecular basis of the FH-deficient HLRCC tumours.

Genomic microarrays in chromosomal analysis of leukemia

Chromosomal alterations in leukemia have been shown to have prognostic and predictive significance and are also important minimal residual disease (MRD) markers in the follow-up of leukemia patients. Although specific oncogenes and tumor suppressors have been discovered in some of the chromosomal alterations, the role and target genes of many alterations in leukemia remain unknown. In addition, a number of leukemia patients have a normal karyotype by standard cytogenetics, but have variability in clinical course and are often molecularly heterogeneous. Cytogenetic methods traditionally used in leukemia analysis and diagnostics; G-banding, various fluorescence in situ hybridization (FISH) techniques, and chromosomal comparative genomic hybridization (cCGH), have enormously increased knowledge about the leukemia genome, but have limitations in resolution or in genomic coverage. In the last decade, the development of microarray comparative genomic hybridization (array-CGH, aCGH) for DNA copy number analysis and the SNP microarray (SNP-array) method for simultaneous copy number and loss of heterozygosity (LOH) analysis has enabled investigation of chromosomal and gene alterations genome-wide with high resolution and high throughput. In these studies, genetic alterations were analyzed in acute myeloid leukemia (AML) and chronic lymphocytic leukemia (CLL). The aim was to screen and characterize genomic alterations that could play role in leukemia pathogenesis by using aCGH and SNP-arrays. One of the most important goals was to screen cryptic alterations in karyotypically normal leukemia patients. In addition, chromosomal changes were evaluated to narrow the target regions, to find new markers, and to obtain tumor suppressor and oncogene candidates. The work presented here shows the capability of aCGH to detect submicroscopic copy number alterations in leukemia, with information about breakpoints and genes involved in the alterations, and that genome-wide microarray analyses with aCGH and SNP-array are advantageous methods in the research and diagnosis of leukemia. The most important findings were the cryptic changes detected with aCGH in karyotypically normal AML and CLL, characterization of amplified genes in 11q marker chromosomes, detection of deletion-based mechanisms of MLL-ARHGEF12 fusion gene formation, and detection of LOH without copy number alteration in karyotypically normal AML. These alterations harbor candidate oncogenes and tumor suppressors for further studies.

Molecular Alterations in Asbestos-Related Lung Cancer

Background: Asbestos is a well known cancer-causing mineral fibre, which has a synergistic effect on lung cancer risk in combination with tobacco smoking. Several in vitro and in vivo experiments have demonstrated that asbestos can evoke chromosomal damage and cause alterations as well as gene expression changes. Lung tumours, in general, have very complex karyotypes with several recurrently gained and lost chromosomal regions and this has made it difficult to identify specific molecular changes related primarily to asbestos exposure. The main aim of these studies has been to characterize asbestos-related lung cancer at a molecular level. Methods: Samples from asbestos-exposed and non-exposed lung cancer patients were studied using array comparative genomic hybridization (aCGH) and fluorescent in situ hybridization (FISH) to detect copy number alterations (CNA) as well as microsatellite analysis to detect allelic imbalance (AI). In addition, asbestos-exposed cell lines were studied using gene expression microarrays. Results: Eighteen chromosomal regions showing differential copy number in the lung tumours of asbestos-exposed patients compared to those of non-exposed patients were identified. The most significant differences were detected at 2p21-p16.3, 5q35.3, 9q33.3-q34.11, 9q34.13-q34.3, 11p15.5, 14q11.2 and 19p13.1-p13.3 (p<0.005). The alterations at 2p and 9q were validated and characterized in detail using AI and FISH analysis in a larger study population. Furthermore, in vitro studies were performed to examine the early gene expression changes induced by asbestos in three different lung cell lines. The results revealed specific asbestos-associated gene expression profiles and biological processes as well as chromosomal regions enriched with genes believed to contribute to the common asbestos-related responses in the cell lines. Interestingly, the most significant region enriched with asbestos-response genes was identified at 2p22, close to the previously identified region showing asbestos-related CNA in lung tumours. Additionally, in this thesis, the dysregulated biological processes (Gene Ontology terms) detected in the cell line experiment were compared to dysregulated processes identified in patient samples in a later study (Ruosaari et al., 2008a). Commonly affected processes such as those related to protein ubiquitination, ion transport and surprisingly sensory perception of smell were identified. Conclusions: The identification of specific CNA and dysregulated biological processes shed some light on the underlying genes acting as mediators in asbestos-related lung carcinogenesis. It is postulated that the combination of several asbestos-specific molecular alterations could be used to develop a diagnostic method for the identification of asbestos-related lung cancer.

Sampled longest common prefix array

When augmented with the longest common prefix (LCP) array and some other structures, the suffix array can solve many string processing problems in optimal time and space. A compressed representation of the LCP array is also one of the main building blocks in many compressed suffix tree proposals. In this paper, we describe a new compressed LCP representation: the sampled LCP array. We show that when used with a compressed suffix array (CSA), the sampled LCP array often offers better time/space trade-offs than the existing alternatives. We also show how to construct the compressed representations of the LCP array directly from a CSA

Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line

Background: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. Methods: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. Results: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). Conclusion: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.

Micropillar Array Based Microchips for Electrospray Ionization Mass Spectrometry, Microreactors, and Liquid Chromatographic Separation

This dissertation deals with the design, fabrication, and applications of microscale electrospray ionization chips for mass spectrometry. The microchip consists of microchannel, which leads to a sharp electrospray tip. Microchannel contain micropillars that facilitate a powerful capillary action in the channels. The capillary action delivers the liquid sample to the electrospray tip, which sprays the liquid sample to gas phase ions that can be analyzed with mass spectrometry. The microchip uses a high voltage, which can be utilized as a valve between the microchip and mass spectrometry. The microchips can be used in various applications, such as for analyses of drugs, proteins, peptides, or metabolites. The microchip works without pumps for liquid transfer, is usable for rapid analyses, and is sensitive. The characteristics of performance of the single microchips are studied and a rotating multitip version of the microchips are designed and fabricated. It is possible to use the microchip also as a microreactor and reaction products can be detected online with mass spectrometry. This property can be utilized for protein identification for example. Proteins can be digested enzymatically on-chip and reaction products, which are in this case peptides, can be detected with mass spectrometry. Because reactions occur faster in a microscale due to shorter diffusion lengths, the amount of protein can be very low, which is a benefit of the method. The microchip is well suited to surface activated reactions because of a high surface-to-volume ratio due to a dense micropillar array. For example, titanium dioxide nanolayer on the micropillar array combined with UV radiation produces photocatalytic reactions which can be used for mimicking drug metabolism biotransformation reactions. Rapid mimicking with the microchip eases the detection of possibly toxic compounds in preclinical research and therefore could speed up the research of new drugs. A micropillar array chip can also be utilized in the fabrication of liquid chromatographic columns. Precisely ordered micropillar arrays offer a very homogenous column, where separation of compounds has been demonstrated by using both laser induced fluorescence and mass spectrometry. Because of small dimensions on the microchip, the integrated microchip based liquid chromatography electrospray microchip is especially well suited to low sample concentrations. Overall, this work demonstrates that the designed and fabricated silicon/glass three dimensionally sharp electrospray tip is unique and facilitates stable ion spray for mass spectrometry.

Integration of miRNA and mRNA expression with DNA copy number of Ewing sarcoma cell lines

Ewing sarcoma is an aggressive and poorly differentiated malignancy of bone and soft tissue. It primarily affects children, adolescents, and young adults, with a slight male predominance. It is characterized by a translocation between chromosomes 11 and 22 resulting in the EWSR1-FLI1fusion transcription factor. The aim of this study is to identify putative Ewing sarcoma target genes through an integrative analysis of three microarray data sets. Array comparative genomic hybridization is used to measure changes in DNA copy number, and analyzed to detect common chromosomal aberrations. mRNA and miRNA microarrays are used to measure expression of protein-coding and miRNA genes, and these results integrated with the copy number data. Chromosomal aberrations typically contain also bystanders in addition to the driving tumor suppressor and oncogenes, and integration with expression helps to identify the true targets. Correlation between expression of miRNAs and their predicted target mRNAs is also evaluated to assess the results of post-transcriptional miRNA regulation on mRNA levels. The highest frequencies of copy number gains were identified in chromosome 8, 1q, and X. Losses were most frequent in 9p21.3, which also showed an enrichment of copy number breakpoints relative to the rest of the genome. Copy number losses in 9p21.3 were found have a statistically significant effect on the expression of MTAP, but not on CDKN2A, which is a known tumor-suppressor in the same locus. MTAP was also down-regulated in the Ewing sarcoma cell lines compared to mesenchymal stem cells. Genes exhibiting elevated expression in association with copy number gains and up-regulation compared to the reference samples included DCAF7, ENO2, MTCP1, andSTK40. Differentially expressed miRNAs were detected by comparing Ewing sarcoma cell lines against mesenchymal stem cells. 21 up-regulated and 32 down-regulated miRNAs were identified, includingmiR-145, which has been previously linked to Ewing sarcoma. The EWSR1-FLI1 fusion gene represses miR-145, which in turn targets FLI1 forming a mutually repressive feedback loop. In addition higher expression linked to copy number gains and compared to mesenchymal stem cells, STK40 was also found to be a target of four different miRNAs that were all down-regulated in Ewing sarcoma cell lines compared to the reference samples. SLCO5A1 was identified as the only up-regulated gene within a frequently gained region in chromosome 8. This region was gained in over 90 % of the cell lines, and also with a higher frequency than the neighboring regions. In addition, SLCO5A1 was found to be a target of three miRNAs that were down-regulated compared to the mesenchymal stem cells.

Genomic alterations in Myeloproliferative Neoplasms and Myelodysplastic Syndromes

Myeloproliferative neoplasms (MPN) and myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic disorders whose etiology and molecular pathogenesis are poorly understood. During the past decade, enormous developments in microarray technology and bioinformatics methods have made it possible to mine novel molecular alterations in a large number of malignancies, including MPN and MDS, which has facilitated the detection of new prognostic, predictive and therapeutic biomarkers for disease stratification. By applying novel microarray techniques, we profiled copy number alterations and microRNA (miRNA) expression changes in bone marrow aspirate and blood samples. In addition, we set up and validated an miRNA expression test for bone marrow core biopsies in order to utilize the large archive material available in many laboratories. We also tested JAK2 mutation status and compare it with the in vitro growth pattern of hematologic progenitors cells. In the study focusing on 100 MPN cases, we detected a Janus kinase 2 (JAK2) mutation in 71 cases. We observed spontaneous erythroid colony growth in all mutation-positive cases in addition to nine mutation negative cases. Interestingly, seven JAK2V167F negative ET cases showed spontaneous megakaryocyte colony formation, one case of which also harbored a myeloproliferative leukemia virus oncogene (MPL) mutation. We studied copy number alterations in 35 MPN and 37 MDS cases by using oligonucleotide-based array comparative hybridization (array CGH). Only one essential thrombocythemia (ET) case presented copy number alterations in chromosomes 1q and 13q. In contrast, MDS cases were characterized by numerous novel cryptic chromosomal aberrations with the most common copy number losses at 5q21.3q33.1 and 7q22.1q33, while the most common copy number gain was trisomy 8. As for the study of the bone marrow core biopsy samples, we showed that even though these samples were embedded in paraffin and underwent decalcification, they were reliable sources of miRNA and suitable for array expression analysis. Further, when studying the miRNA expression profiles of the 19 MDS cases, we found that, compared to controls, two miRNAs (one human Epstein-Barr virus (miR-BART13) miRNA and one human (has-miR-671-5p) miRNA) were downregulated, whereas two other miRNAs (hsa-miR-720 and hsa-miR-21) were upregulated. However, we could find no correlation between copy number alterations and microRNA expression when integrating these two data. This thesis brings to light new information about genomic changes implicated in the development of MPN and MDS, and also underlines the power of applying genome-wide array screening techniques in neoplasias. Rapid advances in molecular techniques and the integration of different genomic data will enable the discovery of the biological contexts of many complex disorders, including myeloid neoplasias.

Molecular diagnosis of developmental disorders

ABSTRACT Idiopathic developmental disorders (DDs) affect ~1% of the population worldwide. This being a considerable amount, efforts are being made to elucidate the disease mechanisms. One or several genetic factors cause 30-40% of DDs, and only 10% are caused by environmental factors. The remaining 50% of DD patients go undiagnosed, mostly due to a lack of diagnostic techniques. The cause in most undiagnosed cases is though to be a genetic factor or a combination of genetic and environmental factors. Despite the surge of new technologies entering the market, their implementation into diagnostic laboratories is hampered by costs, lack of information about the expected diagnostic yield, and the wide range of selection. This study evaluates new microarray methods in diagnosing idiopathic DDs, providing information about their added diagnostic value. Study I analysed 150 patients by array comparative genomic hybridization (array CGH, 44K and 244K), with a subsequent 18% diagnostic yield. These results are supported by other studies, indicating an enourmous added diagnostic value of array CGH, compared with conventional cytogenetic analysis. Nevertheless, 80% of the patients remained undiagnosed in Study I. In an effort to diagnose more patients, in Study IV the resolution was increased from 8.9 Kb of the 244K CGH array to 0.7 Kb, by using a single-nucleotide polymorphism (SNP) array. However, no additional pathogenic changes were detected in the 35 patients assessed, and thus, for diagnostic purposes, an array platform with ca 9 Kb resolution appears adequate. The recent vast increase in reports of detected aberrations and associated phenotypes has enabled characterization of several new syndromes first based on a common aberration and thereafter by delineation of common clinical characteristics. In Study II, a familial deletion at 9q22.2q22.32 with variable penetrance was described. Despite several reports of aberrations in the adjacent area at 9q associated with Gorlin syndrome, the patients in this family had a unique phenotype and did not present with the syndrome. In Study III, a familial duplication of chromosome 6p22.2 was described. The duplication caused increased expression of an important enzyme of the γ-aminobutyric acid (GABA) degradation pathway, causing oxidative stress of the brain, and thus, very likely, the mild mental retardation of these patients. These two case studies attempted to pinpoint candidate genes and to resolve the pathogenic mechanism causing the clinical characteristics of the patients. Presenting rare genetic and clinical findings to the international science and medical community enables interpretation of similar findings in other patients. The added value of molecular karyotyping in patients with idiopathic DD is evident. As a first line of testing, arrays with a median resolution of at least 9 Kb should be considered and further characterization of detected aberrations undertaken when possible. Diagnostic whole-exome sequencing may be the best option for patients who remain undiagnosed after high-resolution array analysis.