aCGH local copy number aberrations associated with overall copy number genomic instability in colorectal cancer: coordinate involvement of the regions including BCR and ABL

Identification of copy number alterations in colon cancer from analysis of amplicon-based next generation sequencing data

The objective of this study was to determine the feasibility to detect copy number alterations in colon cancer samples using Next Generation Sequencing data and to elucidate the association between copy number alterations in specific genes and the development of cancer in different colon segments. We report the successful detection of somatic changes in gene copy number in 37 colon cancer patients by analysis of sequencing data through Amplicon CNA Algorithm. Overall, we have found a total of 748 significant copy number alterations in 230 significant genes, of which 143 showed CN losses and 87 showed CN gains. Validation of results was performed on 20 representative genes by quantitative qPCR and/or immunostaining. By this analysis, we have identified 4 genes that were subjected to copy number alterations in tumors arising in all colon segments (defined "common genes") and the presence of copy number alterations in 14 genes that were significantly associated to one specific site (defined "site-associated genes"). Finally, copy number alterations in ASXL1, TSC1 and IL7R turned out to be clinically relevant since the loss of TSC1 and IL7R was associated with advanced stages and/or reduced survival whereas copy number gain of ASXL1 was associated with good prognosis.

Significant evidence of linkage for a gene predisposing to colorectal cancer and multiple primary cancers on 22q11

Objectives:

The genetic basis of colorectal cancer (CRC) is not completely specified. Part of the difficulty in mapping predisposition genes for CRC may be because of phenotypic heterogeneity. Using data from a population genealogy of Utah record linked to a statewide cancer registry, we identified a subset of CRC cases that exhibited familial clustering in excess of that expected for all CRC cases in general, which may represent a genetically homogeneous subset of CRC.

Methods:

Using a new familial aggregation method referred to as the subset genealogic index of familiality (subsetGIF), combined with detailed information from a statewide tumor registry, we identified a subset of CRC cases that exhibited excess familial clustering above that expected for CRC: CRC cases who had at least one other primary tumor at a different site. A genome-wide linkage analysis was performed on a set of high-risk CRC pedigrees that included multiple CRC cases with additional primaries to identify evidence for predisposition loci.

Results:

A total of 13 high-risk CRC pedigrees with multiple CRC cases with other primary cancers were identified. Linkage analysis identified one pedigree with a significant linkage signal at 22q11 (LOD (logarithm (base 10) of odds)=3.39).

Conclusions:

A predisposition gene or variant for CRC that also predisposes to other primary cancers likely resides on chromosome 22q11. The ability to use statewide population genealogy and tumor registry data was critical to identify an informative subset of CRC cases that is possibly more genetically homogeneous than CRC in general, and may have improved statistical power for predisposition locus identification in this study.

A comprehensive characterization of genome-wide copy number aberrations in colorectal cancer reveals novel oncogenes and patterns of alterations

To develop a comprehensive overview of copy number aberrations (CNAs) in stage-II/III colorectal cancer (CRC), we characterized 302 tumors from the PETACC-3 clinical trial. Microsatellite-stable (MSS) samples (n = 269) had 66 minimal common CNA regions, with frequent gains on 20 q (72.5%), 7 (41.8%), 8 q (33.1%) and 13 q (51.0%) and losses on 18 (58.6%), 4 q (26%) and 21 q (21.6%). MSS tumors have significantly more CNAs than microsatellite-instable (MSI) tumors: within the MSI tumors a novel deletion of the tumor suppressor WWOX at 16 q23.1 was identified (p<0.01). Focal aberrations identified by the GISTIC method confirmed amplifications of oncogenes including EGFR, ERBB2, CCND1, MET, and MYC, and deletions of tumor suppressors including TP53, APC, and SMAD4, and gene expression was highly concordant with copy number aberration for these genes. Novel amplicons included putative oncogenes such as WNK1 and HNF4A, which also showed high concordance between copy number and expression. Survival analysis associated a specific patient segment featured by chromosome 20 q gains to an improved overall survival, which might be due to higher expression of genes such as EEF1B2 and PTK6. The CNA clustering also grouped tumors characterized by a poor prognosis BRAF-mutant-like signature derived from mRNA data from this cohort. We further revealed non-random correlation between CNAs among unlinked loci, including positive correlation between 20 q gain and 8 q gain, and 20 q gain and chromosome 18 loss, consistent with co-selection of these CNAs. These results reinforce the non-random nature of somatic CNAs in stage-II/III CRC and highlight loci and genes that may play an important role in driving the development and outcome of this disease.

Identification of candidate predisposing copy number variants in familial and early-onset colorectal cancer patients

In the majority of colorectal cancers (CRCs) under clinical suspicion for a hereditary cause, the disease-causing genetic factors are still to be discovered. To identify such genetic factors we stringently selected a discovery cohort of 41 CRC index patients with microsatellite-stable tumors. All patients were below 40 years of age at diagnosis and/or exhibited an overt family history. We employed genome-wide copy number profiling using high-resolution SNP arrays on germline DNA, which resulted in the identification of novel copy number variants (CNVs) in six patients (15%) encompassing, among others, the cadherin gene CDH18, the bone morphogenetic protein antagonist family gene GREM1, and the breakpoint cluster region gene BCR. In addition, two genomic deletions were encountered encompassing two microRNA genes, hsa-mir-491/KIAA1797 and hsa-mir-646/AK309218. None of these CNVs has previously been reported in relation to CRC predisposition in humans, nor were they encountered in large control cohorts (>1,600 unaffected individuals). Since several of these newly identified candidate genes may be functionally linked to CRC development, our results illustrate the potential of this approach for the identification of novel candidate genes involved in CRC predisposition.

Genetic and expression analysis of MET, MACC1, and HGF in metastatic colorectal cancer: response to met inhibition in patient xenografts and pathologic correlations

PURPOSE

We determined the gene copy numbers for MET, for its transcriptional activator MACC1 and for its ligand hepatocyte growth factor (HGF) in liver metastases from colorectal carcinoma (mCRC). We correlated copy numbers with mRNA levels and explored whether gain and/or overexpression of MET and MACC1 predict response to anti-Met therapies. Finally, we assessed whether their genomic or transcriptional deregulation correlates with pathologic and molecular parameters of aggressive disease.

EXPERIMENTAL DESIGN

One hundred three mCRCs were analyzed. Copy numbers and mRNA were determined by quantitative PCR (qPCR). Thirty nine samples were implanted and expanded in NOD (nonobese diabetic)/SCID (severe combined immunodeficient) mice to generate cohorts that were treated with the Met inhibitor JNJ-38877605. In silico analysis of MACC1 targets relied on genome-wide mapping of promoter regions and on expression data from two CRC datasets.

RESULTS

No focal, high-grade amplifications of MET, MACC1, or HGF were detected. Chromosome 7 polysomy and gain of the p-arm were observed in 21% and 8% of cases, respectively, and significantly correlated with higher expression of both Met and MACC1. Met inhibition in patient-derived xenografts did not modify tumor growth. Copy number gain and overexpression of MACC1 correlated with unfavorable pathologic features better than overexpression of Met. Bioinformatic analysis of putative MACC1 targets identified elements besides Met, whose overexpression cosegregated with aggressive forms of colorectal cancer.

CONCLUSIONS

Experiments in patient-derived xenografts suggest that mCRCs do not rely on Met genomic gain and/or overexpression for growth. On the basis of pathologic correlations and bioinformatic analysis, MACC1 could contribute to CRC progression through mechanisms other than or additional to Met transcriptional upregulation.

Clinical omics analysis of colorectal cancer incorporating copy number aberrations and gene expression data

Background:

Colorectal cancer (CRC) is one of the most frequently occurring cancers in Japan, and thus a wide range of methods have been deployed to study the molecular mechanisms of CRC. In this study, we performed a comprehensive analysis of CRC, incorporating copy number aberration (CRC) and gene expression data. For the last four years, we have been collecting data from CRC cases and organizing the information as an “omics” study by integrating many kinds of analysis into a single comprehensive investigation.

In our previous studies, we had experienced difficulty in finding genes related to CRC, as we observed higher noise levels in the expression data than in the data for other cancers.

Because chromosomal aberrations are often observed in CRC, here, we have performed a combination of CNA analysis and expression analysis in order to identify some new genes responsible for CRC.

This study was performed as part of the Clinical Omics Database Project at Tokyo Medical and Dental University. The purpose of this study was to investigate the mechanism of genetic instability in CRC by this combination of expression analysis and CNA, and to establish a new method for the diagnosis and treatment of CRC.

Materials and methods:

Comprehensive gene expression analysis was performed on 79 CRC cases using an Affymetrix Gene Chip, and comprehensive CNA analysis was performed using an Affymetrix DNA Sty array. To avoid the contamination of cancer tissue with normal cells, laser micro-dissection was performed before DNA/RNA extraction. Data analysis was performed using original software written in the R language.

Result:

We observed a high percentage of CNA in colorectal cancer, including copy number gains at 7, 8q, 13 and 20q, and copy number losses at 8p, 17p and 18. Gene expression analysis provided many candidates for CRC-related genes, but their association with CRC did not reach the level of statistical significance. The combination of CNA and gene expression analysis, together with the clinical information, suggested UGT2B28, LOC440995, CXCL6, SULT1B1, RALBP1, TYMS, RAB12, RNMT, ARHGDIB, S1000A2, ABHD2, OIT3 and ABHD12 as genes that are possibly associated with CRC. Some of these genes have already been reported as being related to CRC. TYMS has been reported as being associated with resistance to the anti-cancer drug 5-fluorouracil, and we observed a copy number increase for this gene. RALBP1, ARHGDIB and S100A2 have been reported as oncogenes, and we observed copy number increases in each. ARHGDIB has been reported as a metastasis-related gene, and our data also showed copy number increases of this gene in cases with metastasis.

Conclusion:

The combination of CNA analysis and gene expression analysis was a more effective method for finding genes associated with the clinicopathological classification of CRC than either analysis alone. Using this combination of methods, we were able to detect genes that have already been associated with CRC. We also identified additional candidate genes that may be new markers or targets for this form of cancer.

The evolutionary mechanism of cancer

Identification of the general molecular mechanism of cancer is the Holy Grail of cancer research. Since cancer is believed to be caused by a sequential accumulation of cancer gene mutations, the identification, characterization, and targeting of common genetic alterations and their defined pathways have dominated the field for decades. Despite the impressive data accumulated from studies of gene mutations, epigenetic dysregulation, and pathway alterations, an overwhelming amount of diverse molecular information has offered limited understanding of the general mechanisms of cancer. To solve this paradox, the newly established genome theory is introduced here describing how somatic cells evolve within individual patients. The evolutionary mechanism of cancer is characterized using only three key components of somatic cell evolution that include increased system dynamics induced by stress, elevated genetic and epigenetic heterogeneity, and genome alteration mediated natural selection. Cancer progression represents a macro-evolutionary process where karyotype change or genome replacement plays the key dominant role. Furthermore, the recently identified relationship between the evolutionary mechanism and a large number of diverse individual molecular mechanisms is discussed. The total sum of all the individual molecular mechanisms is equal to the evolutionary mechanism of cancer. Individual molecular mechanisms including all the molecular mechanisms described to date are stochastically selected and unpredictable and are therefore clinically impractical. Recognizing the fundamental importance of the underlying basis of the evolutionary mechanism of cancer mandates the development of new strategies in cancer research.

Intratumoural cytogenetic heterogeneity of sporadic colorectal carcinomas suggests several pathways to liver metastasis

Much has been learned about the chromosomal abnormalities of colorectal carcinomas but the cytogenetic relationship between the neoplastic clones present in primary versus metastatic tumour samples remains unclear. We analyse the frequency of abnormalities for 47 chromosome regions using the interphase fluorescence in situ hybridization technique in a group of 48 tumours, including 24 primary colorectal tumours and 24 paired liver metastases. All tumours showed complex karyotypes with numerical/structural abnormalities for seven or more different chromosomes/chromosome regions both in the primary tumours and in their paired metastases. Chromosome 8 was the most frequently altered (22/24 primary tumours), consistently showing del(8p22) and/or gains/amplification of 8q24, followed by abnormalities of the entire chromosome 7 (21/24 primary tumours) and of chromosomes 17p and 20q (20/24 primary tumours). Simultaneous staining for multiple chromosome probes revealed the presence of two or more tumour cell clones in 23/24 cases (46/48 tumour samples). Interestingly, the liver metastases typically contained tumour cell clones similar to those found in the primary tumours, suggesting the absence of selective selection of specific tumour clones. Despite this, additional chromosomal abnormalities were detected in 23/24 metastatic tumours, which preferentially consisted of del(17p13) and gains/amplification of 11q13 and 20q13; moreover, compared to primary tumours, metastases showed an increased number of abnormalities of chromosomes 1p, 7q, 8q, 13q, and 18q, and new chromosomal abnormalities involving chromosomes 6, 10q23, 14q32, 15q22, and 19q13. Owing to the high frequency of numerical abnormalities of the entire chromosome 7 and loss and/or gain/amplification of specific regions of chromosome 8, eg del(8p22) and/or gains/amplification of 8q24 in primary colorectal tumours with associated metastases, it is suggested that their assessment at diagnosis could be of great clinical utility for the identification of colorectal cancer patients at higher risk of developing liver metastases.

Identification and analysis of specific chromosomal region adjacent to exogenous Dhfr-amplified region in Chinese hamster ovary cell genome

Chinese hamster ovary (CHO) cells are widely used for the stable production of recombinant proteins. Gene amplification techniques are frequently used to improve of protein production, and the dihydrofolate reductase (DHFR) gene amplification system is most widely used in the CHO cell line. We previously constructed a CHO genomic bacterial artificial chromosome (BAC) library from a mouse Dhfr-amplified CHO DR1000L-4N cell line and one BAC clone (Cg0031N14) containing the CHO genomic DNA sequence adjacent to Dhfr was selected. To identify the specific chromosomal region adjacent to the exogenous Dhfr-amplified region in the CHO cell genome, we performed further screening of BAC clones to obtain other Dhfr-amplified regions in the CHO genome. From the screening by high-density replica filter hybridization using a digoxigenin-labeled pSV2-dhfr/hGM-CSF probe, we obtained 8 new BAC clones containing a Dhfr-amplified region. To define the structures of the 8 BAC clones, Southern blot analysis, BAC end sequencing and fluorescence in situ hybridization (FISH) were performed. These results revealed that all the selected BAC clones contained a large palindrome structure with a small inverted repeat in the junction region. This suggests that the obtained amplicon structure in the Dhfr-amplified region in the CHO genome plays an important role in exogenous gene amplification.

Genetic and epigenetic heterogeneity in cancer: a genome-centric perspective

Genetic and epigenetic heterogeneity (the main form of non-genetic heterogeneity) are key elements in cancer progression and drug resistance, as they provide needed population diversity, complexity, and robustness. Despite drastically increased evidence of multiple levels of heterogeneity in cancer, the general approach has been to eliminate the "noise" of heterogeneity to establish genetic and epigenetic patterns. In particular, the appreciation of new types of epigenetic regulation like non-coding RNA, have led to the hope of solving the mystery of cancer that the current genetic theories seem to be unable to achieve. In this mini-review, we have briefly analyzed a number of mis-conceptions regarding cancer heterogeneity, followed by the re-evaluation of cancer heterogeneity within a framework of the genome-centric concept of evolution. The analysis of the relationship between gene, epigenetic and genome level heterogeneity, and the challenges of measuring heterogeneity among multiple levels have been discussed. Further, we propose that measuring genome level heterogeneity represents an effective strategy in the study of cancer and other types of complex diseases, as emphasis on the pattern of system evolution rather than specific pathways provides a global and synthetic approach. Compared to the degree of heterogeneity, individual molecular pathways will have limited predictability during stochastic cancer evolution where genome dynamics (reflected by karyotypic heterogeneity) will dominate.

Family cancer syndromes: inherited deficiencies in systems for the maintenance of genomic integrity

Familial cancer syndromes have revealed important fundamental features regarding how all cancers arise through destabilization of the genome, such that somatic evolution can select for the disruption of critical cellular coordinating and regulatory features. The authors examine those cellular genes and systems whose normal role is to preserve genomic integrity and relate them to the genetic foundations of heritable cancers. By examining how these cellular systems normally function, how family cancer genes are able to affect the process of tumor progression can be learned. In so doing, a clearer picture of how sporadic cancers arise is additionally gained.

Uniparentalism in sporadic colorectal cancer is independent of imprint status, and coordinate for chromosomes 14 and 18

Our previous allelotyping studies of 59 sporadic colorectal cancers revealed that loss of heterozygosity is most frequent for regions of chromosomes 14 and 18. Yet subsequent BAC microarray comparative genomic hybridization studies of the same tumor DNAs showed no corresponding pattern of copy number alteration for chromosome 14. To clarify this apparent discrepancy, we utilized hybridization to SNP microarrays; this revealed frequent uniparentalism for chromosome 14 and for chromosome 18. Based on the BAC array results combined with fluorescent in situ hybridization data, it was evident that uniparental disomy was occurring in many colorectal cancers as well as in additional chromosomes, and often coordinately involved chromosomes 14 and 18. Further studies examined the possibility that uniparentalism was directed towards the selection for imprinted genes, but no association with imprinting was observed.

Biological and functional analysis of statistically significant pathways deregulated in colon cancer by using gene expression profiles

Gene expression profiling offers a great opportunity for studying multi-factor diseases and for understanding the key role of genes in mechanisms which drive a normal cell to a cancer state. Single gene analysis is insufficient to describe the complex perturbations responsible for cancer onset, progression and invasion. A deeper understanding of the mechanisms of tumorigenesis can be reached focusing on deregulation of gene sets or pathways rather than on individual genes. We apply two known and statistically well founded methods for finding pathways and biological processes deregulated in pathological conditions by analyzing gene expression profiles. In particular, we measure the amount of deregulation and assess the statistical significance of predefined pathways belonging to a curated collection (Molecular Signature Database) in a colon cancer data set. We find that pathways strongly involved in different tumors are strictly connected with colon cancer. Moreover, our experimental results show that the study of complex diseases through pathway analysis is able to highlight genes weakly connected to the phenotype which may be difficult to detect by using classical univariate statistics. Our study shows the importance of using gene sets rather than single genes for understanding the main biological processes and pathways involved in colorectal cancer. Our analysis evidences that many of the genes involved in these pathways are strongly associated to colorectal tumorigenesis. In this new perspective, the focus shifts from finding differentially expressed genes to identifying biological processes, cellular functions and pathways perturbed in the phenotypic conditions by analyzing genes co-expressed in a given pathway as a whole, taking into account the possible interactions among them and, more importantly, the correlation of their expression with the phenotypical conditions.

Aberrant crypt foci as precursors in colorectal cancer progression

Colorectal cancer progression originates when accumulated genetic and epigenetic alterations cause genomic instability and a malignant phenotype. Subsequent molecular pathway deregulation leads to histopathologic changes that are clinically evident as aberrant crypt foci (ACF) and visualized by high-magnification chromoscopic colonoscopy. ACF are biomarkers of increased colorectal cancer risk, particularly those with dysplastic features. Genetic profiling using genomic instability, loss of heterozygosity, and methylation analysis has revealed a minority population of ACF genotypically analogous to cancer.

A probe-density-based analysis method for array CGH data: simulation, normalization and centralization

MOTIVATION

Genomic instability is one of the fundamental factors in tumorigenesis and tumor progression. Many studies have shown that copy-number abnormalities at the DNA level are important in the pathogenesis of cancer. Array comparative genomic hybridization (aCGH), developed based on expression microarray technology, can reveal the chromosomal aberrations in segmental copies at a high resolution. However, due to the nature of aCGH, many standard expression data processing tools, such as data normalization, often fail to yield satisfactory results.

RESULTS

We demonstrated a novel aCGH normalization algorithm, which provides an accurate aCGH data normalization by utilizing the dependency of neighboring probe measurements in aCGH experiments. To facilitate the study, we have developed a hidden Markov model (HMM) to simulate a series of aCGH experiments with random DNA copy number alterations that are used to validate the performance of our normalization. In addition, we applied the proposed normalization algorithm to an aCGH study of lung cancer cell lines. By using the proposed algorithm, data quality and the reliability of experimental results are significantly improved, and the distinct patterns of DNA copy number alternations are observed among those lung cancer cell lines.

SUPPLEMENTARY INFORMATION

Source codes and.gures may be found at http://ntumaps.cgm.ntu.edu.tw/aCGH_supplementary.