Integration of DNA copy number alterations and transcriptional expression analysis in human gastric cancer

Background

Genomic instability with frequent DNA copy number alterations is one of the key hallmarks of carcinogenesis. The chromosomal regions with frequent DNA copy number gain and loss in human gastric cancer are still poorly defined. It remains unknown how the DNA copy number variations contributes to the changes of gene expression profiles, especially on the global level.

Principal Findings

We analyzed DNA copy number alterations in 64 human gastric cancer samples and 8 gastric cancer cell lines using bacterial artificial chromosome (BAC) arrays based comparative genomic hybridization (aCGH). Statistical analysis was applied to correlate previously published gene expression data obtained from cDNA microarrays with corresponding DNA copy number variation data to identify candidate oncogenes and tumor suppressor genes. We found that gastric cancer samples showed recurrent DNA copy number variations, including gains at 5p, 8q, 20p, 20q, and losses at 4q, 9p, 18q, 21q. The most frequent regions of amplification were 20q12 (7/72), 20q12–20q13.1 (12/72), 20q13.1–20q13.2 (11/72) and 20q13.2–20q13.3 (6/72). The most frequent deleted region was 9p21 (8/72). Correlating gene expression array data with aCGH identified 321 candidate oncogenes, which were overexpressed and showed frequent DNA copy number gains; and 12 candidate tumor suppressor genes which were down-regulated and showed frequent DNA copy number losses in human gastric cancers. Three networks of significantly expressed genes in gastric cancer samples were identified by ingenuity pathway analysis.

Conclusions

This study provides insight into DNA copy number variations and their contribution to altered gene expression profiles during human gastric cancer development. It provides novel candidate driver oncogenes or tumor suppressor genes for human gastric cancer, useful pathway maps for the future understanding of the molecular pathogenesis of this malignancy, and the construction of new therapeutic targets.     

Differential gene alteration among hepatoma cell lines demonstrated by cDNA microarray-based comparative genomic hybridization

We assayed chromosomal abnormalities in hepatoma cell lines using the microarray-based comparative genomic hybridization (array-CGH) method and investigated the relationship between genomic copy number alterations and expression profiles in these hepatoma cell lines. We modified a cDNA array-CGH assay to compare genomic DNAs from seven hepatoma cell lines, as well as DNA from two non-hepatoma cell lines and from normal cells. The mRNA expression of each sample was assayed in parallel by cDNA microarray. We identified small amplified or deleted chromosomal regions, as well as alterations in DNA copy number not previously described. We predominantly found alterations of apoptosis-related genes in Hep3B and HepG2, cell adhesion and receptor molecules in HLE, and cytokine-related genes in PLC/PRF/5. About 40% of the genes showing amplification or loss showed altered levels of mRNA (p < 0.05). Hierarchical clustering analysis showed that the expression of these genes allows differentiation between alpha-fetoprotein (AFP)-producing and AFP-negative cell lines. cDNA array-CGH is a sensitive method that can be used to detect alterations in genomic copy number in tumor cells. Differences in DNA copy alterations between AFP-producing and AFP-negative cells may lead to differential gene expression and may be related to the phenotype of these cells.     

Detection of low level genomic alterations by comparative genomic hybridization based on cDNA micro-arrays

MOTIVATION: The accumulation of genomic alterations is an important process in tumor formation and progression. Comparative genomic hybridization performed on cDNA arrays (cDNA aCGH) is a common method to investigate the genomic alterations on a genome-wide scale. However, when detecting low-level DNA copy number changes this technology requires the use of noise reduction strategies due to a low signal to noise ratio. RESULTS: Currently a running average smoothing filter is the most frequently used noise reduction strategy. We analyzed this strategy theoretically and experimentally and found that it is not sensitive to very low level genomic alterations. The presence of systematic errors in the data is one of the main reasons for this failure. We developed a novel algorithm which efficiently reduces systematic noise and allows for the detection of low-level genomic alterations. The algorithm is based on comparison of the biological relevant data to data from so-called self-self hybridizations, additional experiments which contain no biological information but contain systematic errors. We find that with our algorithm the effective resolution for +/-1 DNA copy number changes is about 2 Mb. For copy number changes larger than three the effective resolution is on the level of single genes.     

Comparative genomics and experimental promoter analysis reveal functional liver-specific elements in mammalian hepatic lipase genes

Background

Array-based comparative genomic hybridization (aCGH) is a high-throughput method for measuring genome-wide DNA copy number changes. Current aCGH methods have limited resolution, sensitivity and reproducibility. Microarrays for aCGH are available only for a few organisms and combination of aCGH data with expression data is cumbersome.

Results

We present a novel method of using commercial oligonucleotide expression microarrays for aCGH, enabling DNA copy number measurements and expression profiles to be combined using the same platform. This method yields aCGH data from genomic DNA without complexity reduction at a median resolution of approximately 17,500 base pairs. Due to the well-defined nature of oligonucleotide probes, DNA amplification and deletion can be defined at the level of individual genes and can easily be combined with gene expression data.

Conclusion

A novel method of gene resolution analysis of copy number variation (graCNV) yields high-resolution maps of DNA copy number changes and is applicable to a broad range of organisms for which commercial oligonucleotide expression microarrays are available. Due to the standardization of oligonucleotide microarrays, graCNV results can reliably be compared between laboratories and can easily be combined with gene expression data using the same platform.     

Array-comparative genomic hybridization reveals loss of SOCS6 is associated with poor prognosis in primary lung squamous cell carcinoma

Background

Primary tumor recurrence commonly occurs after surgical resection of lung squamous cell carcinoma (SCC). Little is known about the genes driving SCC recurrence.

Methods

We used array comparative genomic hybridization (aCGH) to identify genes affected by copy number alterations that may be involved in SCC recurrence. Training and test sets of resected primary lung SCC were assembled. aCGH was used to determine genomic copy number in a training set of 62 primary lung SCCs (28 with recurrence and 34 with no evidence of recurrence) and the altered copy number of candidate genes was confirmed by quantitative PCR (qPCR). An independent test set of 72 primary lung SCCs (20 with recurrence and 52 with no evidence of recurrence) was used for biological validation. mRNA expression of candidate genes was studied using qRT-PCR. Candidate gene promoter methylation was evaluated using methylation microarrays and Sequenom EpiTYPER analysis.

Results

18q22.3 loss was identified by aCGH as being significantly associated with recurrence (p = 0.038). Seven genes within 18q22.3 had aCGH copy number loss associated with recurrence but only SOCS6 copy number was both technically replicated by qPCR and biologically validated in the test set. SOCS6 copy number loss correlated with reduced mRNA expression in the study samples and in the samples with copy number loss, there was a trend for increased methylation, albeit non-significant. Overall survival was significantly poorer in patients with SOCS6 loss compared to patients without SOCS6 loss in both the training (30 vs. 43 months, p = 0.023) and test set (27 vs. 43 months, p = 0.010).

Conclusion

Reduced copy number and mRNA expression of SOCS6 are associated with disease recurrence in primary lung SCC and may be useful prognostic biomarkers.     

Expression of catalase and myeloperoxidase genes in hydrogen peroxide-resistant HL-60 cells.

We studied the expression of catalase and myeloperoxidase genes in the hydrogen peroxide-resistant variants of human myeloid leukemia HL-60 cells HP50-2 and HP100-1. Southern blot hybridization with catalase and myeloperoxidase cDNA probes indicated that the copy number of the catalase gene in HP50-2 and HP100-1 cells was two and eight times, respectively, higher than that in HL-60 cells, whereas the copy number of the myeloperoxidase gene was the same. The amplified catalase and c-myc genes in HP100-1 cells were not decreased by treatment of the cells with inhibitors of poly(ADP-Ribose) polymerase, such as nicotinamide and benzamide. RNA blot hybridization with cDNA probes indicated that the content of catalase mRNA in HP50-2 and HP100-1 cells was four and 16 times higher, respectively, than that in HL-60 cells. By contrast, the content of myeloperoxidase mRNA in HP50-2 and HP100-1 cells was only a few percent of that in HL-60 cells. Furthermore, fluorescent in situ hybridization of a catalase cDNA probe to chromosomes indicated that the catalase gene in HP100-1 was amplified in the p13 region of a derivative chromosome 11. These results indicate that the increased synthesis of catalase in these resistant cells is mainly due to increased expression of the catalase gene, and that the lack of myeloperoxidase synthesis in these cells is due to the absence of its mRNA.     

Inverted duplication pattern in anaphase bridges confirms the breakage-fusion-bridge (BFB) cycle model for 11q13 amplification

Conventional cytogenetic analyses and comparative genomic hybridization have revealed a complex and even chaotic nature of chromosomal aberrations in pleural malignant mesothelioma (MM). We set out to describe the complex gene copy number changes and screen for novel genetic aberrations using a high-density oligonucleotide microarray platform for comparative genomic hybridization (aCGH) of a series of 26 well-characterized MM tumor samples. The number of copy number changes varied from zero to 40 per sample. Gene copy number losses predominated over gains, and the most frequent region of loss was 9p21.3 (17/26 cases), the locus of CDKN2A and CDKN2B, both known to be commonly lost in MM. The most recurrent minimal regions of losses were 1p31.1--> p13.2, 3p22.1-->p14.2, 6q22.1, 9p21.3, 13cen-->q14.12, 14q22.1-->qter, and 22qcen-->q12.3. Previously unreported gains included 9p13.3, 7p22.3-->p22.2, 12q13.3, and 17q21.32-->qter. The results suggest that gene copy number losses are a major mechanism of MM carcinogenesis and reveal a recurrent pattern of copy number changes in MM.     

GATA6 activates Wnt signaling in pancreatic cancer by negatively regulating the Wnt antagonist Dickkopf-1

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal disease characterized by late diagnosis and treatment resistance. Recurrent genetic alterations in defined genes in association with perturbations of developmental cell signaling pathways have been associated with PDAC development and progression. Here, we show that GATA6 contributes to pancreatic carcinogenesis during the temporal progression of pancreatic intraepithelial neoplasia by virtue of Wnt pathway activation. GATA6 is recurrently amplified by both quantitative-PCR and fluorescent in-situ hybridization in human pancreatic intraepithelial neoplasia and in PDAC tissues, and GATA6 copy number is significantly correlated with overall patient survival. Forced overexpression of GATA6 in cancer cell lines enhanced cell proliferation and colony formation in soft agar in vitro and growth in vivo, as well as increased Wnt signaling. By contrast siRNA mediated knockdown of GATA6 led to corresponding decreases in these same parameters. The effects of GATA6 were found to be due to its ability to bind DNA, as forced overexpression of a DNA-binding mutant of GATA6 had no effects on cell growth in vitro or in vivo, nor did they affect Wnt signaling levels in these same cells. A microarray analysis revealed the Wnt antagonist Dickopf-1 (DKK1) as a dysregulated gene in association with GATA6 knockdown, and direct binding of GATA6 to the DKK1 promoter was confirmed by chromatin immunoprecipitation and electrophoretic mobility shift assays. Transient transfection of GATA6, but not mutant GATA6, into cancer cell lines led to decreased DKK1 mRNA expression and secretion of DKK1 protein into culture media. Forced overexpression of DKK1 antagonized the effects of GATA6 on Wnt signaling in pancreatic cancer cells. These findings illustrate that one mechanism by which GATA6 promotes pancreatic carcinogenesis is by virtue of its activation of canonical Wnt signaling via regulation of DKK1.     

Chemical Inhibition of Histone Deacetylases 1 and 2 Induces Fetal Hemoglobin through Activation of GATA2

Therapeutic intervention aimed at reactivation of fetal hemoglobin protein (HbF) is a promising approach for ameliorating sickle cell disease (SCD) and β-thalassemia. Previous studies showed genetic knockdown of histone deacetylase (HDAC) 1 or 2 is sufficient to induce HbF. Here we show that ACY-957, a selective chemical inhibitor of HDAC1 and 2 (HDAC1/2), elicits a dose and time dependent induction of γ-globin mRNA (HBG) and HbF in cultured primary cells derived from healthy individuals and sickle cell patients. Gene expression profiling of erythroid progenitors treated with ACY-957 identified global changes in gene expression that were significantly enriched in genes previously shown to be affected by HDAC1 or 2 knockdown. These genes included GATA2, which was induced greater than 3-fold. Lentiviral overexpression of GATA2 in primary erythroid progenitors increased HBG, and reduced adult β-globin mRNA (HBB). Furthermore, knockdown of GATA2 attenuated HBG induction by ACY-957. Chromatin immunoprecipitation and sequencing (ChIP-Seq) of primary erythroid progenitors demonstrated that HDAC1 and 2 occupancy was highly correlated throughout the GATA2 locus and that HDAC1/2 inhibition led to elevated histone acetylation at well-known GATA2 autoregulatory regions. The GATA2 protein itself also showed increased binding at these regions in response to ACY-957 treatment. These data show that chemical inhibition of HDAC1/2 induces HBG and suggest that this effect is mediated, at least in part, by histone acetylation-induced activation of the GATA2 gene.     

High-resolution oligonucleotide array-CGH pinpoints genes involved in cryptic losses in chronic lymphocytic leukemia

Although recurrent chromosomal alterations occur in chronic lymphocytic leukemia (CLL), relatively few affected tumor suppressors and oncogenes have been implicated. To improve genetic characterization of CLL, we performed high-resolution gene copy number analysis of 20 CLL patients using oligonucleotide array comparative genomic hybridization (aCGH). The most recurrent losses were observed in 13q and 11q with variable sizes. The 11q losses varied between 7.44 Mb and 41.72 Mb in size and targeted ATM among others. Lost regions in 13q were generally smaller, spanning from 0.79 Mb to 29.33 Mb. The minimal common region (158 kb) in 13q14.3, which was also homozygously deleted in some cases, harbored five genes: TRIM13, KCNRG, DLEU2, DLEU1, and FAM10A4. Additionally, two micro-RNA genes (MIRN15A and MIRN16-1) locate to the region. New cryptic losses were detected in 1q23.2-->q23.3, 3p21.31, 16pter-->p13.3, 17p13.3-->p13.2, 17q25.3-->qter, and 22q11.22. In conclusion, our oligonucleotide aCGH study revealed novel aberrations and provided detailed genomic profiles of the altered regions.     

Copy Number Gains at 8q24 and 20q11-q13 in Gastric Cancer Are More Common in Intestinal-Type than Diffuse-Type

The present study was aimed at discovering DNA copy number alterations (CNAs) involved in the carcinogenesis of stomach and at understanding their clinicopathological significances in the Korean population. DNA copy numbers were analyzed using Agilent 244K or 400K array comparative genomic hybridization (aCGH) in fresh-frozen tumor and matched normal tissues from 40 gastric cancer patients. Some of the detected CNA regions were validated using multiplex ligation-dependent probe amplification (MLPA) in six of the 40 patients and customized Agilent 60K aCGH in an independent set of 48 gastric cancers. The mRNA levels of genes at common CNA regions were analyzed using quantitative real-time PCR. Copy number gains were more common than losses across the entire genome in tumor tissues compared to matched normal tissues. The mean number of alterations per case was 64 for gains and 40 for losses, and the median aberration length was 44016 bp for gains and 4732 bp for losses. Copy number gains were frequently detected at 7p22.1 (20%), 8q24.21 (27%–30%), 8q24.3 (22%–48%), 13q34 (20%–31%), and 20q11-q13 (25%–30%), and losses at 3p14.2 (43%), 4q35.2 (27%), 6q26 (23%), and 17p13.3 (20%–23%). CNAs at 7p22.1, 13q34, and 17p13.3 have not been reported in other populations. Most of the copy number losses were associated with down-regulation of mRNA levels, but the correlation between copy number gains and mRNA expression levels varied in a gene-dependent manner. In addition, copy number gains tended to occur more commonly in intestinal-type cancers than in diffuse-type cancers. In conclusion, the present study suggests that copy number gains at 8q24 and 20q11-q13 and losses at 3p14.2 may be common events in gastric cancer but CNAs at 7p22.1, 13q34, and 17p13.3 may be Korean-specific.     

Alternatively spliced transcripts of the Drosophila tramtrack gene encode zinc finger proteins with distinct DNA binding specificities.

A protein present in nuclear extracts of Drosophila embryos binds multiple sites in the promoter and genetically defined autoregulatory element of the pair-rule gene even-skipped (eve). We reported here the isolation of a cDNA encoding this binding activity, the sequence of which identifies it as the 69 kDa zinc finger tramtrack (ttk) protein. As ttk was previously implicated in controlling the expression of another pair-rule gene, fushi tarazu (ftz), our findings suggest that ttk plays a role in the regulation of at least two developmentally important genes. An additional ttk-related cDNA clone was isolated which gives rise to an 88 kDa protein with an alternative set of zinc fingers having a DNA binding specificity distinct from that of the 69 kDa protein. Both proteins were shown to be encoded by the ttk gene through alternative splicing, providing the first example of the use of this mechanism to generate related proteins with distinct DNA binding specificities. Whole mount in situ hybridization analysis revealed different patterns of embryonic expression of the two ttk mRNA isoforms.     

Analysis of fibrinogen genes in patients with congenital afibrinogenemia

Several cDNA clones coding for A alpha, B beta and gamma chains of fibrinogen have been isolated from a human liver cDNA library. They were selected by differential hybridization with probes raised against fractionated liver mRNA (positive probes) and muscle and albumin mRNA (negative probes), then firmly identified by positive hybridization selection. Three of these clones, encoding A alpha, B beta and gamma fibrinogen chain sequences, were further characterized by restriction mapping and used as probes to characterize fibrinogen mRNAs from adult and fetal liver and fibrinogen genes in normal individuals and two afibrinogenemic patients. The results indicate that there is a single copy of the fibrinogen genes which are present and grossly intact in afibrinogenemic DNA.