Search

Abstract Background Xenografts have been shown to provide a suitable source of tumor tissue for molecular analysis in the absence of primary tumor material. We utilized ES xenograft series for integrated microarray analyses to identify novel biomarkers. Method Microarray technology (array comparative genomic hybridization (aCGH) and micro RNA arrays) was used to screen and identify copy number changes and differentially expressed miRNAs of 34 and 14 passages, respectively. Incubated cells used for xenografting (Passage 0) were considered to represent the primary tumor. Four important differentially expressed miRNAs (miR-31, miR-31*, miR-145, miR-106) were selected for further validation by real time polymerase chain reaction (RT-PCR). Integrated analysis of aCGH and miRNA data was performed on 14 xenograft passages by bioinformatic methods. Results The most frequent losses and gains of DNA copy number were detected at 9p21.3, 16q and at 8, 15, 17q21.32-qter, 1q21.1-qter, respectively. The presence of these alterations was consistent in all tumor passages. aCGH profiles of xenograft passages of each series resembled their corresponding primary tumors (passage 0). MiR-21, miR-31, miR-31*, miR-106b, miR-145, miR-150*, miR-371-5p, miR-557 and miR-598 showed recurrently altered expression. These miRNAS were predicted to regulate many ES-associated genes, such as genes of the IGF1 pathway, EWSR1, FLI1 and their fusion gene (EWS-FLI1). Twenty differentially expressed miRNAs were pinpointed in regions carrying altered copy numbers. Conclusion In the present study, ES xenografts were successfully applied for integrated microarray analyses. Our findings showed expression changes of miRNAs that were predicted to regulate many ES associated genes, such as IGF1 pathway genes, FLI1, EWSR1, and the EWS-FLI1 fusion genes.

Abstract Background Only 29 cases of constitutional 9q22 deletions have been published and all have been sporadic. Most associate with Gorlin syndrome or nevoid basal cell carcinoma syndrome (NBCCS, MIM #109400) due to haploinsufficiency of the PTCH1 gene (MIM *601309). Methods and Results We report two mentally retarded female siblings and their cognitively normal father, all carrying a similar 5.3 Mb microdeletion at 9q22.2q22.32, detected by array CGH (244 K). The deletion does not involve the PTCH1 gene, but instead 30 other gene,s including the ROR2 gene (MIM *602337) which causing both brachydactyly type 1 (MIM #113000) and Robinow syndrome (MIM #268310), and the immunologically active SYK gene (MIM *600085). The deletion in the father was de novo and FISH analysis of blood lymphocytes did not suggest mosaicism. All three patients share similar mild dysmorphic features with downslanting palpebral fissures, narrow, high bridged nose with small nares, long, deeply grooved philtrum, ears with broad helix and uplifted lobuli, and small toenails. All have significant dysarthria and suffer from continuous middle ear and upper respiratory infections. The father also has a funnel chest and unilateral hypoplastic kidney but the daughters have no malformations. Conclusions This is the first report of a familial constitutional 9q22 deletion and the first deletion studied by array-CGH which does not involve the PTCH1 gene. The phenotype and penetrance are variable and the deletion found in the cognitively normal normal father poses a challenge in genetic counseling.

Abstract 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/

Abstract Background Ewing sarcoma family of tumors (ESFT), characterized by t(11;22)(q24;q12), is one of the most common tumors of bone in children and young adults. In addition to EWS/FLI1 gene fusion, copy number changes are known to be significant for the underlying neoplastic development of ESFT and for patient outcome. Our genome-wide high-resolution analysis aspired to pinpoint genomic regions of highest interest and possible target genes in these areas. Methods Array comparative genomic hybridization (CGH) and expression arrays were used to screen for copy number alterations and expression changes in ESFT patient samples. A total of 31 ESFT samples were analyzed by aCGH and in 16 patients DNA and RNA level data, created by expression arrays, was integrated. Time of the follow-up of these patients was 5–192 months. Clinical outcome was statistically evaluated by Kaplan-Meier/Logrank methods and RT-PCR was applied on 42 patient samples to study the gene of the highest interest. Results Copy number changes were detected in 87% of the cases. The most recurrent copy number changes were gains at 1q, 2, 8, and 12, and losses at 9p and 16q. Cumulative event free survival (ESFT) and overall survival (OS) were significantly better (P < 0.05) for primary tumors with three or less copy number changes than for tumors with higher number of copy number aberrations. In three samples copy number imbalances were detected in chromosomes 11 and 22 affecting the FLI1 and EWSR1 loci, suggesting that an unbalanced t(11;22) and subsequent duplication of the derivative chromosome harboring fusion gene is a common event in ESFT. Further, amplifications on chromosomes 20 and 22 seen in one patient sample suggest a novel translocation type between EWSR1 and an unidentified fusion partner at 20q. In total 20 novel ESFT associated putative oncogenes and tumor suppressor genes were found in the integration analysis of array CGH and expression data. Quantitative RT-PCR to study the expression levels of the most interesting gene, HDGF, confirmed that its expression was higher than in control samples. However, no association between HDGF expression and patient survival was observed. Conclusion We conclude that array CGH and integration analysis proved to be effective methods to identify chromosome regions and novel target genes involved in the tumorigenesis of ESFT.

Abstract Background Studies on asbestos-induced tumourigenesis have indicated the role of, e.g., reactive oxygen/nitrogen species, mitochondria, as well as NF-κB and MAPK signalling pathways. The exact molecular mechanisms contributing to asbestos-mediated carcinogenesis are, however, still to be characterized. Methods In this study, gene expression data analyses together with gene annotation data from the Gene Ontology (GO) database were utilized to identify pathways that are differentially regulated in lung and tumour tissues between asbestos-exposed and non-exposed lung cancer patients. Differentially regulated pathways were identified from gene expression data from 14 asbestos-exposed and 14 non-exposed lung cancer patients using custom-made software and Iterative Group Analysis (iGA). Western blotting was used to further characterize the findings, specifically to determine the protein levels of UBA1 and UBA7. Results Differences between asbestos-related and non-related lung tumours were detected in pathways associated with, e.g., ion transport, NF-κB signalling, DNA repair, as well as spliceosome and nucleosome complexes. A notable fraction of the pathways down-regulated in both normal and tumour tissue of the asbestos-exposed patients were related to protein ubiquitination, a versatile process regulating, for instance, DNA repair, cell cycle, and apoptosis, and thus being also a significant contributor of carcinogenesis. Even though UBA1 or UBA7, the early enzymes involved in protein ubiquitination and ubiquitin-like regulation of target proteins, did not underlie the exposure-related deregulation of ubiquitination, a difference was detected in the UBA1 and UBA7 levels between squamous cell carcinomas and respective normal lung tissue (p = 0.02 and p = 0.01) without regard to exposure status. Conclusion Our results indicate alterations in protein ubiquitination related both to cancer type and asbestos. We present for the first time pathway analysis results on asbestos-associated lung cancer, providing important insight into the most relevant targets for future research.

Abstract Background Urothelial carcinomas originate from the epithelial cells of the inner lining of the bladder and may appear as single or as multiple synchronous tumors. Patients with urothelial carcinomas frequently show recurrences after treatment making follow-up necessary. The leading hypothesis explaining the origin of meta- and synchronous tumors assumes a monoclonal origin. However, the genetic relationship among consecutive tumors has been shown to be complex in as much as the genetic evolution does not adhere to the chronological appearance of the metachronous tumors. Consequently, genetically less evolved tumors may appear chronologically later than genetically related but more evolved tumors. Methods Forty-nine meta- or synchronous urothelial tumors from 22 patients were analyzed using expression profiling, conventional CGH, LOH, and mutation analyses. Results We show by CGH that partial chromosomal losses in the initial tumors may not be present in the recurring tumors, by LOH that different haplotypes may be lost and that detected regions of LOH may be smaller in recurring tumors, and that mutations present in the initial tumor may not be present in the recurring ones. In contrast we show that despite apparent genomic differences, the recurrent and multiple bladder tumors from the same patients display remarkably similar expression profiles. Conclusion Our findings show that even though the vast majority of the analyzed meta- and synchronous tumors from the same patients are not likely to have originated directly from the preceding tumor they still show remarkably similar expressions profiles. The presented data suggests that an expression profile is established early in tumor development and that this profile is stable and maintained in recurring tumors.

Abstract Background DNA amplifications alter gene dosage in cancer genomes by multiplying the gene copy number. Amplifications are quintessential in a considerable number of advanced cancers of various anatomical locations. The aims of this study were to classify human cancers based on their amplification patterns, explore the biological and clinical fundamentals behind their amplification-pattern based classification, and understand the characteristics in human genomic architecture that associate with amplification mechanisms. Methods We applied a machine learning approach to model DNA copy number amplifications using a data set of binary amplification records at chromosome sub-band resolution from 4400 cases that represent 82 cancer types. Amplification data was fused with background data: clinical, histological and biological classifications, and cytogenetic annotations. Statistical hypothesis testing was used to mine associations between the data sets. Results Probabilistic clustering of each chromosome identified 111 amplification models and divided the cancer cases into clusters. The distribution of classification terms in the amplification-model based clustering of cancer cases revealed cancer classes that were associated with specific DNA copy number amplification models. Amplification patterns – finite or bounded descriptions of the ranges of the amplifications in the chromosome – were extracted from the clustered data and expressed according to the original cytogenetic nomenclature. This was achieved by maximal frequent itemset mining using the cluster-specific data sets. The boundaries of amplification patterns were shown to be enriched with fragile sites, telomeres, centromeres, and light chromosome bands. Conclusions Our results demonstrate that amplifications are non-random chromosomal changes and specifically selected in tumor tissue microenvironment. Furthermore, statistical evidence showed that specific chromosomal features co-localize with amplification breakpoints and link them in the amplification process.

Abstract Background Asbestos has been shown to cause chromosomal damage and DNA aberrations. Exposure to asbestos causes many lung diseases e.g. asbestosis, malignant mesothelioma, and lung cancer, but the disease-related processes are still largely unknown. We exposed the human cell lines A549, Beas-2B and Met5A to crocidolite asbestos and determined time-dependent gene expression profiles by using Affymetrix arrays. The hybridization data was analyzed by using an algorithm specifically designed for clustering of short time series expression data. A canonical correlation analysis was applied to identify correlations between the cell lines, and a Gene Ontology analysis method for the identification of enriched, differentially expressed biological processes. Results We recognized a large number of previously known as well as new potential asbestos-associated genes and biological processes, and identified chromosomal regions enriched with genes potentially contributing to common responses to asbestos in these cell lines. These include genes such as the thioredoxin domain containing gene (TXNDC) and the potential tumor suppressor, BCL2/adenovirus E1B 19kD-interacting protein gene (BNIP3L), GO-terms such as "positive regulation of I-kappaB kinase/NF-kappaB cascade" and "positive regulation of transcription, DNA-dependent", and chromosomal regions such as 2p22, 9p13, and 14q21. We present the complete data sets as Additional files. Conclusion This study identifies several interesting targets for further investigation in relation to asbestos-associated diseases.