Gene Copy Number Profiles Research Articles

Pathway discovery in mantle cell lymphoma by integrated analysis of high-resolution gene expression and copy number profiling

AbstractThe genome of mantle cell lymphoma (MCL) is, in addition to the translocation t(11;14), characterized by a high number of secondary chromosomal gains and losses that probably account for the various survival times of MCL patients. We investigated 77 primary MCL tumors with available clinical information using high-resolution RNA expression and genomic profiling and applied our recently developed gene expression and dosage integrator algorithm to identify novel genes and pathways that may be of relevance for the pathobiology of MCL. We show that copy number neutral loss of heterozygosity is common in MCL and targets regions that are frequently affected by deletions. The molecular consequences of genomic copy number changes appear complex, even in genomic loci with identified tumor suppressors, such as the region 9p21 containing the CDKN2A locus. Moreover, the deregulation of novel genes, such as CUL4A, ING1, and MCPH1, may affect the 2 crucial pathogenetic mechanisms in MCL, the disturbance of the proliferation, and DNA damage response pathways. Deregulation of the Hippo pathway may have a pathogenetic role in MCL because decreased expression of its members MOBKL2A, MOBKL2B, and LATS2 was associated with inferior outcome, including an independent validation series of 32 MCLs.

Prognostic classification of patients with acute lymphoblastic leukemia by using gene copy number profiles identified from array-based comparative genomic hybridization data

The development of risk-adapted therapy has improved the treatment results of acute lymphoblastic leukemia (ALL) especially in children. However, more accurate risk classifiers are warranted. In this study we aimed at defining a prognostic classifier based on DNA copy number alterations of adolescent and young adult (AYA) (10–25 yrs) ALL patients ( n = 60) determined by microarray CGH and the relapse status of the patients. As a result of prognostic model identification procedure, we got a model of four genes: BAK1, CDKN2C, GSTM1, and MT1F, the copy number profile combinations of which differentiated AYA ALL patients at diagnosis depending on their risk of relapse. The performance of the model was poorer on other age groups. We suggest that this kind of approach produces models simple and accurate enough for potential use in ALL routine classification.

Genome-Wide Profiling of Archived Material from CALGB 9840 and 9342 for Paclitaxel (P) and Trastuzumab (T) Response Biomarkers Using Gene Expression and Copy Number Analysis.

Abstract Background:Emerging data suggest that RNA obtained from formalin-fixed, paraffin-embedded (FFPE) tissue can produce reliable gene expression profiles. Archived material from two taxane monotherapy studies, CALGB 9342 (comparison of 3 doses of P) and 9840 (weekly vs. q3 week P)were profiled for gene expression and DNA copy number.Methods: A total of 238 patients had primary tumor blocks available from a combined sample size of 680. DNA and RNA was extracted from 1.5mm punch cores and the Ambion Recover-All kit™. quality was measured by spectrophotometric analysis, Bioanalyzer RNA Integrity Number (RIN), and housekeeping genes (RPL13A and Actin). A custom DASL™ array containing 779 genes in two-fold redundancy was designed with genes selected to represent the PAM50 intrinsic subtypes, the Oncotype Dx Score, the Netherlands prognostic signature and the genes most frequently found on recurrent breast cancer amplicons. Several methods for identifying outliers were evaluated, including principal components analysis, pairwise correlations as well as the reproducibility of the platform based on replicate samples.Results: Adequate RNA was obtained from 237/238 of these cases which ranged in age from 12-18 years. Of these, 215/237DASL arrays passed further quality control measures. Adequate DNA for CGH was obtained from 227/238 samples. Analysis of PAM50 intrinsic subtypes showed an excess of basal-like tumors (30%) in the primary tumors of this metastatic cohort compared with expected frequency in an early stage population. Luminal A tumors were less frequent than expected (20%). Patients with basal-like tumors did far worse than other tumor types for both PFS on P (p=0.015) and OS (p=2.7X10-6), which persisted in multivariable analysis (p=0.0047), however the interaction term was not significant (Wald p=0.26). While basal-like tumors had similar PFS and OS on both weekly and q3 week P, luminal A tumors appear to achieve more benefit from weekly P (p= 0.0041).The HER2-enriched expression subtype had a similar prognosis to Luminal A and B tumors. This appeared to be due to the presence of T, as the addition of this agent improved PFS (p=0.026) and OS (p=2.0X10-4). Of note, some centrally confirmed HER2 FISH amplified tumors were classified into luminal A, B, and basal-like subtypes. These tumors have similar prognoses to the overall group, for example the basal-like and HER2 tumors had a poor prognosis despite T (p=0.00086). This suggests that HER2 FISH positive tumors may behave based on the underlying tumor subtype. Sawtooth genomes (45% vs 15%) were more frequent than predicted by an early stage tumor dataset as were simplex genomes (3% vs 24%). HER2 by FISH and CGH were highly concordant suggesting data on gene amplification from this platform is robust.Conclusions: Gene expression and copy number profiling of FFPE material from archived tumor blocks (>10 years) produces quality data for biomarker discovery in CALGB clinical trial datasets. These tools allow discovery of novel patterns of gene expression and genomic aberrations that are associated with differential response to P and T. Further studies using these platforms should be performed. Citation Information: Cancer Res 2009;69(24 Suppl):Abstract nr 4032.

Meta-mining of Neuroblastoma and Neuroblast Gene Expression Profiles Reveals Candidate Therapeutic Compounds

Neuroblastoma is a heterogeneous childhood tumor with poor survival outcome for the aggressive type despite intensive multimodal therapies. In this study, we aimed to identify new treatment options for neuroblastoma based on integrative genomic analysis. The Connectivity Map is a database comprising expression profiles in response to known therapeutic compounds. This renders it a useful tool in the search for potential therapeutic compounds based on comparison of gene expression profiles of diseased cells and a database of profiles in response to known therapeutic compounds. We have used this strategy in the search for new therapeutic molecules for neuroblastoma based on data of an integrative meta-analysis of gene copy number and expression profiles from 146 primary neuroblastoma tumors and normal fetal neuroblasts. In a first step, a 132-gene classifier was established that discriminates three major genomic neuroblastoma subgroups, reflecting inherent differences in gene expression between these subgroups. Subsequently, we screened the Connectivity Map database using gene lists generated by comparing expression profiles of fetal adrenal neuroblasts and the genomic subgroups of neuroblastomas. A putative therapeutic effect was predicted for several compounds of which six were empirically tested. A significant reduction in cell viability was shown for five of these molecules: 17-allylamino-geldanamycin, monorden, fluphenazine, trichostatin, and rapamycin. This proof-of-principle study indicates that an integrative genomic meta-analysis approach with inclusion of neuroblast data enables the identification of promising compounds for treatment of children with neuroblastoma. Further studies are warranted to explore in detail the therapeutic potential of these compounds.

Array comparative genomic hybridization analysis of olfactory neuroblastoma

Olfactory neuroblastoma is an unusual neuroectodermal malignancy, which is thought to arise at the olfactory membrane of the sinonasal tract. Due to its rarity, little is understood regarding its molecular and cytogenetic abnormalities. The aim of the current study is to identify specific DNA copy number changes in olfactory neuroblastoma. Thirteen dissected tissue samples were analyzed using array comparative genomic hybridization. Our results show that gene copy number profiles of olfactory neuroblastoma samples are complex. The most frequent changes included gains at 7q11.22–q21.11, 9p13.3, 13q, 20p/q, and Xp/q, and losses at 2q31.1, 2q33.3, 2q37.1, 6q16.3, 6q21.33, 6q22.1, 22q11.23, 22q12.1, and Xp/q. Gains were more frequent than losses, and high-stage tumors showed more alterations than low-stage olfactory neuroblastoma. Frequent changes in high-stage tumors were gains at 13q14.2–q14.3, 13q31.1, and 20q11.21–q11.23, and loss of Xp21.1 (in 66% of cases). Gains at 5q35, 13q, and 20q, and losses at 2q31.1, 2q33.3, and 6q16–q22, were present in 50% of cases. The identified regions of gene copy number change have been implicated in a variety of tumors, especially carcinomas. In addition, our results indicate that gains in 20q and 13q may be important in the progression of this cancer, and that these regions possibly harbor genes with functional relevance in olfactory neuroblastoma.

311 POSTER Early discovery of genomic correlates of drug response by high-resolution gene copy number profiling

Loss of heterozygosity (LOH) is frequently associated with the inactivation of tumor suppressor genes. 18q LOH has been frequently reported in colorectal cancer and lung cancer; however, allelic loss on 18q has not been investigated in renal cell carcinoma (RCC). We evaluated LOH on 18q using nine microsatellite markers in 126 with conventional RCC (cRCC). LOH was observed in more than one 18q microsatellite locus in 24 cRCC (19%). We found the highest frequency of LOH (13.5%) at 18q21.3, where the DCC gene is located. We also assessed the relationship between LOH frequency and patient clinical parameters. Patients with a family history of cancer had a significantly higher frequency of 18q LOH than those without such a history (P = 0.0017). No associations were found with other parameters, including gender, tumor grade, tumor stage, smoking status, and body mass index. The results suggest that inactivation of tumor suppressor genes at 18q21.3, including DCC and SMAD4 as candidates, may be involved in the tumorigenesis of some conventional RCCs.

Genomic imbalances in Schistosoma-associated and non– Schistosoma-associated bladder carcinoma. An array comparative genomic hybridization analysis

Carcinoma of the urinary bladder is the most common malignancy in many tropical and subtropical countries due to endemic infection by Schistosoma hematobium (bilharzia). In the current study, we performed a high-resolution analysis of gene copy number amplifications using array comparative genomic hybridization to compare DNA copy number changes in pools of Schistosoma-associated (SA) and non– Schistosoma-associated (NSA) bladder cancer (BC). Many DNA copy number changes were detected in all studies, with multiple gains and losses of genetic material. The most frequent alterations were gains on 5p15.2∼p15.33, 8q13.1, and 11q13, and losses on 8p21.3∼p22 and 22q13. Even when SA pools showed no Schistosoma-specific gene copy number profiling as compared to NSA pools, some genes seemed to be gained ( ELN on 7q11.23) and some lost ( PRKAG3 on 2q35 and PRDM6 on 5q23.2) in SA-SCC. The following genes were gained in all histopathologic categories: SRC (20q11.23), CEBPB (20q13.13), and GPR9 (Xq13.1). Our study did not provide clear evidence of differences in carcinogenesis of SA-BC and NSA-BC.

Gene expression and copy number profiling suggests the importance of allelic imbalance in 19p in asbestos-associated lung cancer

Asbestos is a pulmonary carcinogen known to give rise to DNA and chromosomal damage, but the exact carcinogenic mechanisms are still largely unknown. In this study, gene expression arrays were performed on lung tumor samples from 14 heavily asbestos-exposed and 14 non-exposed patients matched for other characteristics. Using a two-step statistical analysis, 47 genes were revealed that could differentiate the tumors of asbestos-exposed from those of non-exposed patients. To identify asbestos-associated regions with DNA copy number and expressional changes, the gene expression data were combined with comparative genomic hybridization microarray data. As a result, a combinatory profile of DNA copy number aberrations and expressional changes significantly associated with asbestos exposure was obtained. Asbestos-related areas were detected in 2p21-p16.3, 3p21.31, 5q35.2-q35.3, 16p13.3, 19p13.3-p13.1 and 22q12.3-q13.1. The most prominent of these, 19p13, was further characterized by microsatellite analysis in 62 patients for the differences in allelic imbalance (AI) between the two groups of lung tumors. 79% of the exposed and 45% of the non-exposed patients (P=0.008) were found to be carriers of AI in their lung tumors. In the exposed group, AI in 19p was prevalent regardless of the histological tumor type. In adenocarcinomas, AI in 19p appeared to occur independently of the asbestos exposure.

Gene copy number profiling of soft-tissue leiomyosarcomas by array-comparative genomic hybridization

Leiomyosarcoma (LMS) is a rare malignant mesenchymal tumor of smooth muscle cells. Chromosomal aberrations in LMS have been studied, but the cytogenetic data that have been published so far are complex, limited, and incomplete. Here, we performed for the first time a high-resolution genome-wide array comparative genomic hybridization (CGH) analysis (aCGH) on a pool of 14 low- and high-grade LMS cases to obtain gene-level information about the amplified and deleted regions that may play a role in the development and progression of LMS. Our aCGH results indicated that 2,218 genes were involved in 25 altered chromosomal regions; 9 regions in low-grade LMS, 12 regions in high-grade LMS, and 4 minimal common regions shared by low- and high-grade LMS. The frequency of DNA copy number gains in high-grade LMS was threefold compared to low-grade LMS, whereas losses in low-grade LMS were almost twice as frequent as in high-grade LMS. Both low- and high-grade tumors shared two minimal common regions of gain (15q26∼qter and 17p13.1∼q11) and loss (6p12∼p21.3 and 13q14.3∼qter). Moreover, our findings indicated that low- and high-grade LMS and osteosarcoma share 12 genes located in the 17p amplicon. In conclusion, by using aCGH, we were able to define the precise location of the altered chromosomal areas and to identify putative tumor suppressor genes and oncogenes therein. The list of altered genes in the minimal common regions is available as Appendix 1 at our web site ( http://www.helsinki.fi/cmg/microarray_data).

High‐resolution gene copy number and expression profiling of human chromosome 22 in ovarian carcinomas

Previous low-resolution studies of chromosome 22 in ovarian carcinoma have suggested its involvement in the development of the disease. We report a high-resolution analysis of DNA copy number and gene expression of 22q in 18 ovarian carcinomas using a 22q-specific genomic microarray. We identified aberrations in 67% of the studied tumors, which displayed 3 distinct gene copy number profiles. The majority of the cases (11 of 18) demonstrated heterozygous terminal deletions of various sizes, the smallest of which was 3.5 Mb. The second profile, detected in 3 tumors, revealed the coexistence of heterozygous deletions and different patterns of low-copy-number gain that involved the proximal half of 22q. The latter finding has not been reported previously in ovarian carcinoma. One case displayed a continuous deletion encompassing the entire 22q, consistent with monosomy 22. Furthermore, we compared the results with the available data on these tumors by using cDNA microarrays to define the degree of correlation between abnormalities at the DNA level and variation in mRNA expression. By a comparison with the expression data, we were able to identify 21 deleted genes showing low mRNA levels and 12 amplified genes displaying elevated gene expression, several of which play roles in cell cycle control and the induction of apoptosis. Our results indicated significant correlation between DNA copy number aberrations and variation in mRNA expression. We also identified several regions and candidate genes on 22q that should be studied further to determine their role in the development of ovarian cancer.

DNA copy-number analysis of the 22q11 deletion-syndrome region using array-CGH with genomic and PCR-based targets

Deletions and duplications of genomic segments commonly cause developmental disorders. The resolution and efficiency in diagnosing such gene-dosage alterations can be drastically increased using microarray-based comparative genomic hybridization (array-CGH). However, array-CGH currently relies on spotting genomic clones as targets, which confers severe limitations to the approach including resolution of analysis and reliable gene-dosage assessment of regions with high content of redundant sequences. To improve the methodology for analysis, we compared the use of genomic clones, repeat-free pools of amplified genomic DNA and cDNAs (single and pooled) as targets on the array. For this purpose, we chose q11.2 locus on chromosome 22 as a testing ground. Microdeletions at 22q11 cause birth defects collectively described as the DiGeorge/velocardiofacial syndrome. The majority of patients present 3 Mb typical deletions. Here, we report the construction of a gene-dosage array, covering 6 Mb of 22q11 and including the typically deleted region. We hybridized DNA from six DiGeorge syndrome patients to the array, and show that as little as 11.5 kb non-redundant, repeat-free PCR-generated sequence can be used for reliable detection of hemizygous deletions. By extrapolation, this would allow analysis of the genome with an average resolution of 25 kb. In the case of cDNAs our results indicate that 3.5 kb sequence is necessary for accurate identification of haploid/diploid dosage alterations. Thus, for regions rich in redundant sequences and repeats, such as 22q11, a specifically tailored array-CGH approach is good for gene copy number profiling.

Characterization of genomically amplified segments using PCR: Optimizing relative-PCR for reliable and simple gene expression and gene copy analyses

Gene amplification is one of the mechanisms for oncogene activation in solid tumors. The size of the amplified regions may vary considerably among individual tumors, and more than one gene may be affected within the same amplicon. The main objective in analyzing genomic amplifications has therefore been to map the shortest region involved and to identify genes with increased expression as a result of the increased gene copy number. To facilitate such an analysis, we have developed simple polymerase chain reaction (PCR) procedures using the internal standards beta-actin (ACTB) and L1Hs for gene expression and gene copy number analyses, respectively. We used cDNA derived from pancreatic carcinoma cell lines, and genomic DNA extracted from the same cell lines, as templates in the gene expression and in the gene copy number analyses, respectively. To determine the optimal number of PCR cycles, dilution series of the templates were made. Furthermore, competing primers were used to adjust for differences in target sequence levels. We show that by these simple means it is possible to determine optimal conditions for expression analyses. In addition, the procedure was adapted for the analysis of gene copy number changes at the genomic level using L1Hs as the internal standard. This PCR method makes it possible to produce detailed gene copy number profiles of amplified genomic regions.