Genes with relevance for early to late progression of colon carcinoma based on combined genomic and transcriptomic information from the same patients

Novel mouse model recapitulates genome and transcriptome alterations in human colorectal carcinomas

Human colorectal carcinomas are defined by a nonrandom distribution of genomic imbalances that are characteristic for this disease. Often, these imbalances affect entire chromosomes. Understanding the role of these aneuploidies for carcinogenesis is of utmost importance. Currently, established transgenic mice do not recapitulate the pathognonomic genome aberration profile of human colorectal carcinomas. We have developed a novel model based on the spontaneous transformation of murine colon epithelial cells. During this process, cells progress through stages of pre-immortalization, immortalization and, finally, transformation, and result in tumors when injected into immunocompromised mice. We analyzed our model for genome and transcriptome alterations using ArrayCGH, spectral karyotyping (SKY), and array based gene expression profiling. ArrayCGH revealed a recurrent pattern of genomic imbalances. These results were confirmed by SKY. Comparing these imbalances with orthologous maps of human chromosomes revealed a remarkable overlap. We observed focal deletions of the tumor suppressor genes Trp53 and Cdkn2a/p16. High-level focal genomic amplification included the locus harboring the oncogene Mdm2, which was confirmed by FISH in the form of double minute chromosomes. Array-based global gene expression revealed distinct differences between the sequential steps of spontaneous transformation. Gene expression changes showed significant similarities with human colorectal carcinomas. Pathways most prominently affected included genes involved in chromosomal instability and in epithelial to mesenchymal transition. Our novel mouse model therefore recapitulates the most prominent genome and transcriptome alterations in human colorectal cancer, and might serve as a valuable tool for understanding the dynamic process of tumorigenesis, and for preclinical drug testing. © 2016 Wiley Periodicals, Inc.

DNA copy number alterations, gene expression changes and disease-free survival in patients with colorectal cancer: a 10 year follow-up

BACKGROUND

DNA copy number alterations (CNAs) and gene expression changes have amply been encountered in colorectal cancers (CRCs), but the extent at which CNAs affect gene expression, as well as their relevance for tumor development, are still poorly defined. Here we aimed at assessing the clinical relevance of these parameters in a 10 year follow-up study.

METHODS

Tumors and normal adjacent colon mucosa, obtained at primary surgery from 21 CRC patients, were subjected to (i) high-resolution array CGH (a-CGH) for the detection of CNAs and (ii) microarray-based transcriptome profiling for the detection of gene expression (GE) changes. Correlations between these genomic and transcriptomic changes and their associations with clinical and histopathological parameters were assessed with the aim to identify molecular signatures associated with disease-free survival of the CRC patients during a 10 year follow-up.

RESULTS

DNA copy number gains were frequently detected in chromosomes 7, 8q, 13, 19, 20q and X, whereas DNA copy number losses were frequently detected in chromosomes 1p, 4, 8p, 15, 17p, 18, 19 and 22q. None of these alterations were observed in all samples. In addition, we found that 2,498 genes were up- and that 1,094 genes were down-regulated in the tumor samples compared to their corresponding normal mucosa (p < 0.01). The expression of 65 genes was found to be significantly associated with prognosis (p < 0.01). Specifically, we found that up-regulation of the IL17RA, IGF2BP2 and ABCC2 genes, and of genes acting in the mTOR and cytokine receptor pathways, were strongly associated with a poor survival. Subsequent integrated analyses revealed that increased expression levels of the MMP9, BMP7, UBE2C, I-CAM, NOTCH3, NOTCH1, PTGES2, HMGB1 and ERBB3 genes were associated with copy number gains, whereas decreased expression levels of the MUC1, E2F2, HRAS and SIRT3 genes were associated with copy number losses. Pathways related to cell cycle progression, eicosanoid metabolism, and TGF-β and apoptosis signaling, were found to be most significantly affected.

CONCLUSIONS

Our results suggest that CNAs in CRC tumor tissues are associated with concomitant changes in the expression of cancer-related genes. In other genes epigenetic mechanism may be at work. Up-regulation of the IL17RA, IGF2BP2 and ABCC2 genes, and of genes acting in the mTOR and cytokine receptor pathways, appear to be associated with a poor survival. These alterations may, in addition to Dukes' staging, be employed as new prognostic biomarkers for the prediction of clinical outcome in CRC patients.

Unexpected frequency of genomic alterations in histologically normal colonic tissue from colon cancer patients

As shown by genomic studies, colorectal cancer (CRC) is a highly heterogeneous disease, where copy number alterations (CNAs) may greatly vary among different patients. To explore whether CNAs may be present also in histologically normal tissues from patients affected by CRC, we performed CGH + SNP Microarray on 15 paired tumoral and normal samples. Here, we report for the first time the occurrence of CNAs as a common feature of the histologically normal tissue from CRC patients, particularly CNAs affecting different oncogenes and tumor-suppressor genes, including some not previously reported in CRC and others known as being involved in tumor progression. Moreover, from the comparison of normal vs paired tumoral tissue, we were able to identify three groups: samples with an increased number of CNAs in tumoral vs normal tissue, samples with a similar number of CNAs in both tissues, and samples with a decrease of CNAs in tumoral vs normal tissue, which may be likely due to a selection of the cell population within the tumor. In conclusion, our approach allowed us to uncover for the first time an unexpected frequency of genetic alteration in normal tissue, suggesting that tumorigenic genetic lesions are already present in histologically normal colonic tissue and that the use in array comparative genomic hybridization (CGH) studies of normal samples as reference for the paired tumors can lead to misrepresented genomic data, which may be incomplete or limited, especially if used for the research of target molecules for personalized therapy and for the possible correlation with clinical outcome.

Alterations in Tumor DNA Are Related to Short Postoperative Survival in Patients Resected for Pancreatic Carcinoma Aimed at Cure

OBJECTIVES

Pancreatic ductal adenocarcinomas (PDACs) are found in more than 85% of patients with pancreatic cancer and with 5-year survival of less than 10%. Effective treatment may be radical surgery, which is hampered by rapid relapse. Therefore, our aim was to compare DNA sequence alterations in patients with short and long survival to evaluate if confirmed DNA alterations predict short postoperative survival.

METHODS

DNA was extracted from tumor tissue from 59 PDAC patients, analyzed for KRAS mutations, and hybridized to 180 K CGH + SNP microarrays and 450 K methylation arrays. Analyses were based on postoperative survival where less than 12 months was considered to be short survival and more than 18 months was considered long survival.

RESULTS

Ninety-three percent of the patients had KRAS mutations in tumor DNA. Great heterogeneity of whole genome DNA sequence alterations were observed among chromosomes within the patient materials. Specific DNA sequence alterations did not directly predict postoperative survival, although short survivors had significantly more and larger DNA amplifications (P < 0.006). Amplifications on chromosome 11 and 21 and deletions on chromosome 2 predicted short postoperative survival (P < 0.03). DNA methylation was not related to survival.

CONCLUSIONS

Highly variable genetic differences among DNA regions in PDAC tumors were demonstrated. Postoperative short survival was related to tumor sequence DNA alterations on chromosome 2, 11, and 21.

Multi-Scale Genomic, Transcriptomic and Proteomic Analysis of Colorectal Cancer Cell Lines to Identify Novel Biomarkers

Selecting colorectal cancer (CRC) patients likely to respond to therapy remains a clinical challenge. The objectives of this study were to establish which genes were differentially expressed with respect to treatment sensitivity and relate this to copy number in a panel of 15 CRC cell lines. Copy number variations of the identified genes were assessed in a cohort of CRCs. IC50's were measured for 5-fluorouracil, oxaliplatin, and BEZ-235, a PI3K/mTOR inhibitor. Cell lines were profiled using array comparative genomic hybridisation, Illumina gene expression analysis, reverse phase protein arrays, and targeted sequencing of KRAS hotspot mutations. Frequent gains were observed at 2p, 3q, 5p, 7p, 7q, 8q, 12p, 13q, 14q, and 17q and losses at 2q, 3p, 5q, 8p, 9p, 9q, 14q, 18q, and 20p. Frequently gained regions contained EGFR, PIK3CA, MYC, SMO, TRIB1, FZD1, and BRCA2, while frequently lost regions contained FHIT and MACROD2. TRIB1 was selected for further study. Gene enrichment analysis showed that differentially expressed genes with respect to treatment response were involved in Wnt signalling, EGF receptor signalling, apoptosis, cell cycle, and angiogenesis. Stepwise integration of copy number and gene expression data yielded 47 candidate genes that were significantly correlated. PDCD6 was differentially expressed in all three treatment responses. Tissue microarrays were constructed for a cohort of 118 CRC patients and TRIB1 and MYC amplifications were measured using fluorescence in situ hybridisation. TRIB1 and MYC were amplified in 14.5% and 7.4% of the cohort, respectively, and these amplifications were significantly correlated (p≤0.0001). TRIB1 protein expression in the patient cohort was significantly correlated with pERK, Akt, and Caspase 3 expression. In conclusion, a set of candidate predictive biomarkers for 5-fluorouracil, oxaliplatin, and BEZ235 are described that warrant further study. Amplification of the putative oncogene TRIB1 has been described for the first time in a cohort of CRC patients.

Copy number alterations and allelic ratio in relation to recurrence of rectal cancer

Background

In rectal cancer, total mesorectal excision surgery combined with preoperative (chemo)radiotherapy reduces local recurrence rates but does not improve overall patient survival, a result that may be due to the harmful side effects and/or co-morbidity of preoperative treatment. New biomarkers are needed to facilitate identification of rectal cancer patients at high risk for local recurrent disease. This would allow for preoperative (chemo)radiotherapy to be restricted to high-risk patients, thereby reducing overtreatment and allowing personalized treatment protocols. We analyzed genome-wide DNA copy number (CN) and allelic alterations in 112 tumors from preoperatively untreated rectal cancer patients. Sixty-six patients with local and/or distant recurrent disease were compared to matched controls without recurrence. Results were validated in a second cohort of tumors from 95 matched rectal cancer patients. Additionally, we performed a meta-analysis that included 42 studies reporting on CN alterations in colorectal cancer and compared results to our own data.

Results

The genomic profiles in our study were comparable to other rectal cancer studies. Results of the meta-analysis supported the hypothesis that colon cancer and rectal cancer may be distinct disease entities. In our discovery patient study cohort, allelic retention of chromosome 7 was significantly associated with local recurrent disease. Data from the validation cohort were supportive, albeit not statistically significant, of this finding.

Conclusions

We showed that retention of heterozygosity on chromosome 7 may be associated with local recurrence in rectal cancer. Further research is warranted to elucidate the mechanisms and effect of retention of chromosome 7 on the development of local recurrent disease in rectal cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1550-0) contains supplementary material, which is available to authorized users.

Effects of altered expression and activity levels of CK1δ and ɛ on tumor growth and survival of colorectal cancer patients

Colorectal cancer (CRC) is the fourth leading cause of cancer related death worldwide due to high apoptotic resistance and metastatic potential. Because mutations as well as deregulation of CK1 isoforms contribute to tumor development and tumor progression, CK1 has become an interesting drug target. In this study we show that CK1 isoforms are differently expressed in colon tumor cell lines and that growth of these cell lines can be inhibited by CK1-specific inhibitors. Furthermore, expression of CK1δ and ɛ is changed in colorectal tumors compared to normal bowel epithelium, and high CK1ɛ expression levels significantly correlate with prolonged patients' survival. In addition to changes in CK1δ and ɛ expression, mutations within exon 3 of CK1δ were detected in colorectal tumors. These mutations influence ATP binding resulting in changes in kinetic parameters of CK1δ. Overexpression of these mutants in HT29 cells alters their ability to grow anchorage independently. Consistent with these results, these CK1δ mutants lead to differences in proliferation rate and tumor size in xenografts due to changes in gene expression, especially in genes involved in regulation of cell proliferation, cell cycle, and apoptosis. In summary, our results provide evidence that changes in the expression levels of CK1 isoforms in colorectal tumors correlate with patients' survival. Furthermore, CK1 mutants affect growth and proliferation of tumor cells and induce tumor growth in xenografts, leading to the assumption that CK1 isoforms provide interesting targets for the development of novel effective therapeutic concepts to treat colorectal cancer.

Somatic copy-neutral loss of heterozygosity and copy number abnormalities in Malaysian sporadic colorectal carcinoma patients

Colorectal cancer is one of the most common cancers in many countries, including Malaysia. The accumulation of genomic alterations is an important feature of colorectal carcinogenesis. A better understanding of the molecular events underlying the stages of colorectal carcinogenesis might be helpful in the detection and management of the disease. We used a commercially available single-nucleotide polymorphism genotyping array to detect both copy number abnormalities (CNAs) and copy-neutral loss of heterozygosity (LOH) in sporadic colorectal carcinomas. Matched tumor and normal tissues of 13 colorectal carcinomas (Dukes' stages A-D) were analyzed using a 250K single nucleotide polymorphism array. An additional assay was performed to determine the microsatellite instability status by using the National Cancer Institute-recommended BAT-26 panel. In general, copy number gain (92.3%) was most common, followed by copy number loss (53.8%) and copy-neutral LOH (46.2%). Frequent CNAs of gains and losses were observed on chromosomes 7p, 8, 13q, 17p, 18q, and 20q, and copy-neutral LOH was observed on chromosomes 2, 6, 12, 13q, 14q, 17, 20p, 19q, and 22q. Even though genomic alterations are associated with colorectal cancer progression, our results showed that DNA CNAs and copy-neutral LOH do not reflect disease progression in at least 50% tumors. Copy-neutral LOH was observed in both early and advanced tumors, which favors the involvement of these genomic alterations in the early stages of tumor development.

Integrated analysis of genome-wide copy number alterations and gene expression in microsatellite stable, CpG island methylator phenotype-negative colon cancer

Microsatellite stable (MSS), CpG island methylator phenotype (CIMP)-negative colorectal tumors, the most prevalent molecular subtype of colorectal cancer, are associated with extensive copy number alteration (CNA) events and aneuploidy. We report on the identification of characteristic recurrent CNA (with frequency >25%) events and associated gene expression profiles for a total of 40 paired tumor and adjacent normal colon tissues using genome-wide microarrays. We observed recurrent CNAs, namely gains at 1q, 7p, 7q, 8p12-11, 8q, 12p13, 13q, 20p, 20q, Xp, and Xq and losses at 1p36, 1p31, 1p21, 4p15-12, 4q12-35, 5q21-22, 6q26, 8p, 14q, 15q11-12, 17p, 18p, 18q, 21q21-22, and 22q. Within these genomic regions we identified 356 genes with significant differential expression (P < 0.0001 and ±1.5-fold change) in the tumor compared to adjacent normal tissue. Gene ontology and pathway analyses indicated that many of these genes were involved in functional mechanisms that regulate cell cycle, cell death, and metabolism. An amplicon present in >70% of the tumor samples at 20q11-20q13 contained several cancer-related genes (AHCY, POFUT1, RPN2, TH1L, and PRPF6) that were upregulated and demonstrated a significant linear correlation (P < 0.05) for gene dosage and gene expression. Copy number loss at 8p, a CNA associated with adenocarcinoma and poor prognosis, was observed in >50% of the tumor samples and demonstrated a significant linear correlation for gene dosage and gene expression for two potential tumor suppressor genes, MTUS1 (8p22) and PPP2CB (8p12). The results from our integration analysis illustrate the complex relationship between genomic alterations and gene expression in colon cancer.

Genomic CGH-assessed structural DNA alterations in rectal carcinoma as related to local recurrence following primary operation for cure

Several factors determine overall outcome and possible local recurrence after curative surgery for rectal carcinoma. Surgical performance is usually believed to be the most pertinent factor, followed by adjuvant oncological treatment and tumor histopathology. However, chromosomal instability is common in colorectal cancer and tumor clones are assumed to differ in aggressiveness and potential of causing local recurrence. The aim of this study was, therefore, to evaluate if genetic alterations in primary rectal carcinoma are predictive of local recurrences. A large clinical database with linked bio-bank allowed for careful matching of two patient groups (R0) resected for rectal carcinoma. One group had developed early, isolated local recurrences and the other group seemed cured after 93 months follow-up. DNA from the primary tumors was analysed with array-CGH (comparative genomic hybridization) including 55,000 genomic probes. DNA from all primary tumors in both groups displayed previously reported and well-recognised DNA aberrations in colorectal carcinoma. Significant copy number gains were confirmed in the 4q31.1-31.22 region in DNA from tumors with subsequent local recurrence. Twenty-two affected genes in this region code for products with high relevance in tumor biology (p53 regulation, cell cycle activity, transcription). DNA from rectal carcinoma displayed well-known aberrations as described for colon carcinoma with no obvious prediction of local rectal recurrence. Gains in the 4q31.1-31.22 DNA region are highly potential for local recurrence despite R0 resection to be confirmed in larger patient materials.

Genome-wide screening for methylation-silenced genes in colorectal cancer

Identification of methylation-silenced genes in colorectal cancer (CRC) is of great importance. We employed oligonucleotide microarrays to identify differences in global gene expression of five CRC cell lines (HCT116, RKO, Colo320, SW480 and HT29) that were analyzed before and after treatment with 5-aza-2'-deoxycitidine. Selected candidates were subjected to methylation-specific PCR and real-time quantitative reverse transcription-PCR using 15 CRC cell lines and 23 paired tumor and normal samples from CRC patients. After 5-aza-2'-deoxycitidine treatment, 139 genes were re-expressed in all 5 CRC cell lines collectively with a fold change of more than 1.5 in at least one cell line. These genes include known methylated and silenced genes in CRC. After applying study selection criteria we identified 20 candidates. The GADD45B and THSD1 genes were selected for further analysis. Among 15 colon cancer cell lines, methylation was only identified in THSD1 (27%). THSD1 methylation was subsequently investigated in 23 colorectal tumors and methylation was detected in 9% of the analyzed samples; the observed promoter hypermethylation was cancer-specific. THSD1 mRNA down-regulation was observed in tumor tissues. This genome-wide screening led to the identification of genes putatively affected by methylation in CRC. The THSD1 gene may play a role in the tumorigenesis of CRC.

Cell cycle and aging, morphogenesis, and response to stimuli genes are individualized biomarkers of glioblastoma progression and survival

BACKGROUND

Glioblastoma is a complex multifactorial disorder that has swift and devastating consequences. Few genes have been consistently identified as prognostic biomarkers of glioblastoma survival. The goal of this study was to identify general and clinical-dependent biomarker genes and biological processes of three complementary events: lifetime, overall and progression-free glioblastoma survival.

METHODS

A novel analytical strategy was developed to identify general associations between the biomarkers and glioblastoma, and associations that depend on cohort groups, such as race, gender, and therapy. Gene network inference, cross-validation and functional analyses further supported the identified biomarkers.

RESULTS

A total of 61, 47 and 60 gene expression profiles were significantly associated with lifetime, overall, and progression-free survival, respectively. The vast majority of these genes have been previously reported to be associated with glioblastoma (35, 24, and 35 genes, respectively) or with other cancers (10, 19, and 15 genes, respectively) and the rest (16, 4, and 10 genes, respectively) are novel associations. Pik3r1, E2f3, Akr1c3, Csf1, Jag2, Plcg1, Rpl37a, Sod2, Topors, Hras, Mdm2, Camk2g, Fstl1, Il13ra1, Mtap and Tp53 were associated with multiple survival events.Most genes (from 90 to 96%) were associated with survival in a general or cohort-independent manner and thus the same trend is observed across all clinical levels studied. The most extreme associations between profiles and survival were observed for Syne1, Pdcd4, Ighg1, Tgfa, Pla2g7, and Paics. Several genes were found to have a cohort-dependent association with survival and these associations are the basis for individualized prognostic and gene-based therapies. C2, Egfr, Prkcb, Igf2bp3, and Gdf10 had gender-dependent associations; Sox10, Rps20, Rab31, and Vav3 had race-dependent associations; Chi3l1, Prkcb, Polr2d, and Apool had therapy-dependent associations. Biological processes associated glioblastoma survival included morphogenesis, cell cycle, aging, response to stimuli, and programmed cell death.

CONCLUSIONS

Known biomarkers of glioblastoma survival were confirmed, and new general and clinical-dependent gene profiles were uncovered. The comparison of biomarkers across glioblastoma phases and functional analyses offered insights into the role of genes. These findings support the development of more accurate and personalized prognostic tools and gene-based therapies that improve the survival and quality of life of individuals afflicted by glioblastoma multiforme.