A model-based scan statistic for identifying extreme chromosomal regions of gene expression in human tumors

MOTIVATION: The analysis of gene expression data in its chromosomal context has been a recent development in cancer research. However, currently available methods fail to account for variation in the distance between genes, gene density and genomic features (e.g. GC content) in identifying increased or decreased chromosomal regions of gene expression. RESULTS: We have developed a model-based scan statistic that accounts for these aspects of the complex landscape of the human genome in the identification of extreme chromosomal regions of gene expression. This method may be applied to gene expression data regardless of the microarray platform used to generate it. To demonstrate the accuracy and utility of this method, we applied it to a breast cancer gene expression dataset and tested its ability to predict regions containing medium-to-high level DNA amplification (DNA ratio values >2). A classifier was developed from the scan statistic results that had a 10-fold cross-validated classification rate of 93% and a positive predictive value of 88%. This result strongly suggests that the model-based scan statistic and the expression characteristics of an increased chromosomal region of gene expression can be used to accurately predict chromosomal regions containing amplified genes. AVAILABILITY: Functions in the R-language are available from the author upon request. CONTACT: fcouples@umich.edu.     

利用相互染色体涂色技术分析两株人肺腺癌细胞系(A549和GLC-82)的核型特征

肺癌是全世界癌症死亡中的一个主要的原因。除吸烟外,一些肺癌患者的发病与氡气污染相关。该研究采用包括染色体分选、正向和反向染色体涂色技术,分析了两株肺腺癌细胞系A549和GLC-82的核型特征。A549和GLC-82细胞系都属于非小细胞肺癌细胞系,但诱因不同,后者来源于一个长期生活在氡气污染环境肺癌病人的癌组织。染色体涂色结果表明,这两株肺癌细胞系发生了复杂的染色体重排。在A549和 GLC-82细胞系中,除正常染色体拷贝数变化外,还分别存在13条和24条畸变染色体。约一半的畸变染色体是通过非相互易位形成的,其余的畸变染色体则是通过一些正常染色体的片段缺失或重复而产生的。尽管这两株肺癌细胞系都没有共同的畸变染色体, 但它们似共享两个染色体易位断裂点：HSA8q24和12q14。     

Integrated genome-wide gene expression map and high-resolution analysis of aberrant chromosomal regions in squamous cell lung cancer

The recognition of recurrent aberrant regions in cancer is important to the discovery of candidate cancer related genes. Here we first constructed a genome-wide gene expression map of squamous lung carcinoma from the Stanford Microarray Database. High-resolution detection of aberrant chromosomal regions was performed by using moving-median method. 84% (27 of 32) of our results were consistent with the previous studies of comparative genomic hybridization or loss of heterozygosity. One overrepresented region in Xq28 was newly discovered to be related to squamous cell lung carcinoma. These observations could be of great interest for further studies.     

A Mechanism of Gene Amplification Driven by Small DNA Fragments

DNA amplification is a molecular process that increases the copy number of a chromosomal tract and often causes elevated expression of the amplified gene(s). Although gene amplification is frequently observed in cancer and other degenerative disorders, the molecular mechanisms involved in the process of DNA copy number increase remain largely unknown. We hypothesized that small DNA fragments could be the trigger of DNA amplification events. Following our findings that small fragments of DNA in the form of DNA oligonucleotides can be highly recombinogenic, we have developed a system in the yeast Saccharomyces cerevisiae to capture events of chromosomal DNA amplification initiated by small DNA fragments. Here we demonstrate that small DNAs can amplify a chromosomal region, generating either tandem duplications or acentric extrachromosomal DNA circles. Small fragment-driven DNA amplification (SFDA) occurs with a frequency that increases with the length of homology between the small DNAs and the target chromosomal regions. SFDA events are triggered even by small single-stranded molecules with as little as 20-nt homology with the genomic target. A double-strand break (DSB) external to the chromosomal amplicon region stimulates the amplification event up to a factor of 20 and favors formation of extrachromosomal circles. SFDA is dependent on Rad52 and Rad59, partially dependent on Rad1, Rad10, and Pol32, and independent of Rad51, suggesting a single-strand annealing mechanism. Our results reveal a novel molecular model for gene amplification, in which small DNA fragments drive DNA amplification and define the boundaries of the amplicon region. As DNA fragments are frequently found both inside cells and in the extracellular environment, such as the serum of patients with cancer or other degenerative disorders, we propose that SFDA may be a common mechanism for DNA amplification in cancer cells, as well as a more general cause of DNA copy number variation in nature.     

The pattern of gene amplification is determined by the chromosomal location of hairpin-capped breaks

DNA palindromes often colocalize in cancer cells with chromosomal regions that are predisposed to gene amplification. The molecular mechanisms by which palindromes can cause gene amplification are largely unknown. Using yeast as a model system, we found that hairpin-capped double-strand breaks (DSBs) occurring at the location of human Alu-quasipalindromes lead to the formation of intrachromosomal amplicons with large inverted repeats (equivalent to homogeneously staining regions in mammalian chromosomes) or extrachromosomal palindromic molecules (equivalent to double minutes [DM] in mammalian cells). We demonstrate that the specific outcomes of gene amplification depend on the applied selection, the nature of the break, and the chromosomal location of the amplified gene relative to the site of the hairpin-capped DSB. The rules for the palindrome-dependent pathway of gene amplification defined in yeast may operate during the formation of amplicons in human tumors.     

Lack of an Additive Effect between the Deletions of Klf5 and Nkx3-1 in Mouse Prostatic Tumorigenesis

Prostate cancer is one of the most common malignancies.The development and progression of prostate cancer are driven by a series of genetic and epigenetic events including gene amplification that activates oncogenes and chromosomal deletion that inactivates tumor suppressor genes.Whereas gene amplification occurs in human prostate cancer,gene deletion is more common,and a large number of chromosomal regions have been identified to have frequent deletion in prostate cancer,suggesting that tumor suppressor inactivation is more common than oncogene activation in prostatic carcinogenesis (Knuutila et al.,1998,1999;Dong,2001).Among the most frequently deleted chromosomal regions in prostate cancer,target genes such as NKX3-1 from 8p21,PTENfrom 10q23 andATBF1 from 16q22 have been identified by different approaches (He et al.,1997;Li et al.,1997;Sun et al.,2005),and deletion of these genes in mouse prostates has been demonstrated to induce and/or promote prostatic carcinogenesis.For example,knockout of Nkx3-1 in mice induces hyperplasia and dysplasia (Bhatia-Gaur et al.,1999;Abdulkadir et al.,2002) and promotes prostatic tumorigenesis (Abate-Shen et al.,2003),while knockout of Pten alone causes prostatic neoplasia (Wang et al.,2003).Therefore,gene deletion plays a causal role in prostatic carcinogenesis (Dong,2001).     

Etiology of specific molecular alterations in human malignancies

Cancer results from multiple genomic changes that affect DNA and its gene expression. The DNA sequences may be gained, lost or amplified, or translocated into different parts of the genome to form a fusion gene with oncogenic properties. The occurrence of specific chromosomal aberrations may be restricted to only one cancer type and it may be considered a primary carcinogenic event. Furthermore, the aberration profiles may be used to cluster tumors with similar origins. A variety of techniques exist for the detection of specific chromosomal and gene expression changes. However, the etiology of these molecular alterations remains unclear. Here we discuss the roles of Helicobacter pylori and asbestos burden as carcinogens that cause gastric cancer, mesothelioma and lung cancer.     

整合分析基因表达与拷贝数变异识别癌症的驱动基因及调控子miRNAs

目的:通过整合分析基因的表达与拷贝数变异(CNV)识别癌症的驱动基因及调控子mi RNAs。方法:通过整合基因表达与CNV数据,分别计算了乳腺癌、结肠癌、肺癌、肾癌、膀胱癌、头颈癌六种癌症中mi RNAs的调控得分,提出了一个识别驱动基因和显著调控子mi RNAs的方法。结果:本文研究发现,CNV区域上编码的基因相比于非CNV区域上编码的基因更倾向于受mi RNAs调控。但是,癌相关CNV区域上的基因相比正常CNV区域上的基因更少受mi RNAs调控。本研究识别出了EXOSC4、ZNF7、BOP1等原癌基因,以及mi R-488、mi R-27a、mi R-454等在多种癌症中都起调控作用的调控子mi RNAs。结论:本文的方法为癌症研究带来了新的启发,这些具有调控扩增基因过表达作用的mi RNAs的发现,有助于我们更进一步了解癌基因表达的复杂调控机制,进而推动癌症的诊断、治疗和预后。     

Wild-type p53 restores cell cycle control and inhibits gene amplification in cells with mutant p53 alleles.

Loss of cell cycle control and acquisition of chromosomal rearrangements such as gene amplification often occur during tumor progression, suggesting that they may be correlated. We show here that the wild-type p53 allele is lost when fibroblasts from patients with the Li-Fraumeni syndrome (LFS) are passaged in vitro. Normal and LFS cells containing wild-type p53 arrested in G1 when challenged with the uridine biosynthesis inhibitor PALA and did not undergo PALA-selected gene amplification. The converse occurred in cells lacking wild-type p53 expression. Expression of wild-type p53 in transformants of immortal and tumor cells containing mutant p53 alleles restored G1 control and reduced the frequency of gene amplification to undetectable levels. These studies reveal that p53 contributes to a metabolically regulated G1 check-point, and they provide a model for understanding how abnormal cell cycle progression leads to the genetic rearrangements involved in tumor progression.     

From amplification to function: the case of the MDR1 gene.

This review describes the features of gene amplification associated with the selection of multidrug-resistant cell lines. Some of these lines carry multiple copies of the MDR1 gene that encodes P-glycoprotein, a broad specificity efflux pump. The MDR1 gene was initially identified as the common component of the amplicons found in multidrug-resistant cell lines selected with different drugs. Subsequent studies have established that increased MDR1 expression is sufficient for the multidrug-resistant phenotype. MDR1-containing amplicons may include a number of additional transcribed genes that do not appear to contribute to multidrug resistance. MDR1 amplification is associated with specific chromosomal changes and apparently non-random recombinational events. Increased expression of the MDR1 gene, however, does not necessarily require gene amplification. Although amplification of the MDR1 gene has not been found in clinical tumor samples, increased expression of this gene is commonly observed in different types of cancer and appears to be a significant marker of clinical drug resistance.     

Origin and functional significance of large-scale chromosomal imbalances in neuroblastoma

Neuroblastoma is characterized by numerous recurrent large-scale chromosomal imbalances and gene amplifications which are associated with poor clinical outcome. The most common include MYCN amplification, loss of 1p, 3p and 11q, and gain of 17q genomic regions. Two of these abnormalities, MYCN amplification and loss of 11q, define different genetic subtypes of the disease with vastly different global gene expression profiles. The progress towards the identification of the genes and genetic pathways that have been affected by these abnormalities is reviewed and high resolution mapping of the chromosomal breakpoint regions using oligonucleotide array CGH (oaCGH) is discussed. oaCGH analysis is proving useful for both defining minimal regions of overlap of imbalances, as well as providing information on the molecular mechanisms that generate the chromosomal imbalances. These high resolution analyses have also permitted the detection of micro-deletions in the tumors that further assist in identifying genes important for neuroblastoma pathogenesis.     

How a Replication Origin and Matrix Attachment Region Accelerate Gene Amplification under Replication Stress in Mammalian Cells

The gene amplification plays a critical role in the malignant transformation of mammalian cells. The most widespread method for amplifying a target gene in cell culture is the use of methotrexate (Mtx) treatment to amplify dihydrofolate reductase (Dhfr). Whereas, we found that a plasmid bearing both a mammalian origin of replication (initiation region; IR) and a matrix attachment region (MAR) was spontaneously amplified in mammalian cells. In this study, we attempted to uncover the underlying mechanism by which the IR/MAR sequence might accelerate Mtx induced Dhfr amplification. The plasmid containing the IR/MAR was extrachromosomally amplified, and then integrated at multiple chromosomal locations within individual cells, increasing the likelihood that the plasmid might be inserted into a chromosomal environment that permits high expression and further amplification. Efficient amplification of this plasmid alleviated the genotoxicity of Mtx. Clone-based cytogenetic and sequence analysis revealed that the plasmid was amplified in a chromosomal context by breakage-fusion-bridge cycles operating either at the plasmid repeat or at the flanking fragile site activated by Mtx. This mechanism explains how a circular molecule bearing IR/MAR sequences of chromosomal origin might be amplified under replication stress, and also provides insight into gene amplification in human cancer.     

DNA profiling by array comparative genomic hybridization (CGH) of peripheral blood mononuclear cells (PBMC) and tumor tissue cell in non-small cell lung cancer (NSCLC)

Lung tumor cell DNA copy number alteration (CNA) was expected to display specific patterns such as a large-scale amplification or deletion of chromosomal arms, as previously published data have reported. Peripheral blood mononuclear cell (PBMC) CNA however, was expected to show normal variations in cancer patients as well as healthy individuals, and has thus been used as normal control DNA samples in various published studies. We performed array CGH to measure and compare genetic changes in terms of the CNA of PBMC samples as well as DNA isolated from tumor tissue samples, obtained from 24 non-small cell lung cancer patients. Contradictory to expectations, our studies showed that the PBMC CNA also showed chromosomal variant regions. The list included well-known tumor-associated NTRK1, FGF8, TP53, and TGFβ1 genes and potentially novel oncogenes such as THPO (3q27.1), JMJD1B, and EGR1 (5q31.2), which was investigated in this study. The results of this study highlighted the connection between PBMC and tumor cell genomic DNA in lung cancer patients. However, the application of these studies to cancer prognosis may pose a challenge due to the large amount of information contained in genetic predisposition and family history that has to be processed for useful downstream clinical applications.     

Regularities of karyotypic evolution during stepwise amplification of genes determining drug resistance.

Analysis of chromosomal alterations during stepwise development of mdr1, dhfr, or CAD gene amplifications in a large number of independently selected Djungarian hamster DM-15 and murine P388 sublines revealed typical patterns of karyotypic evolution, specific for multiplication of each of these genes in each cell type. Some principal similarities of karyotypic evolution were noted in at least two different systems. They include: (i) appearance at the first selection step of a new chromosomal arm bearing the resident gene copy followed at the next selection steps by the formation in these specific chromosomal arms of amplified DNA tandem arrays; (ii) translocations of amplified DNA from its initial site to other, also non-random, chromosomal sites; and (iii) emergence in the cell variants with high degrees of gene amplification of multiple extra-chromosomal elements. The most prominent distinctions among the systems were as follows: (i) different structures, evidently containing amplified DNAs, appeared at the initial steps of amplification of different genes--additional heterogeneously staining regions in specific chromosomal segments in the case of amplification of dhfr or CAD genes in DM-15 cells, and mini-chromosomes in the case of mdr1 gene amplification in both DM-15 and P380 cells; (ii) distinct patterns of location of the amplified mdr1 gene copies are characteristic of Djungarian hamster DM-15 and murine P388 cell derivatives after subsequent steps of selection--at the site of resident gene localization or in some other, also non-random, chromosomal sites in DM-15 sublines, and predominantly extra-chromosomal in P388 sublines. We propose that different mechanisms are responsible for the initial steps of amplification of dhfr and CAD genes on the one hand and the mdr1 gene on the other: non-equal sister-chromatid exchanges and autonomous replication of the extra-chromosomal elements. It seems, however, that both mechanisms may be involved in further rounds of amplification of each of these three genes.     

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.