Forkhead box E1, frequently downregulted by promoter methylation, inhibits colorectal cancer cell growth and migration

Forkhead box E1 (FOXE1), also known as thyroid transcription factor 2 (TTF-2), belongs to a large family of forkhead transcription factors. It plays important roles in embryogenesis, cell growth, and differentiation. Cancer-specific FOXE1 hypermethylation events have been identified in several cancers. However, the expression and function of FOXE1 in the tumorigenesis of colorectal cancer remain still unknown. In this study, we examined FOXE1 expression and methylation in normal colon mucosa, colorectal cancer (CRC) cell lines, and primary tumors by immunohistochemistry, semi-quantitative RT-PCR, methylation-specific PCR, and bisulfite genomic sequencing. We found that FOXE1 was frequently methylated and silenced in CRC cell lines and was downregulated in CRC tissues compared with paired adjacent non-tumor tissues. Meanwhile, low FOXE1 expression was significantly correlated with lymph node metastasis and advanced TNM stages, indicating its potential as a tumor marker. Subsequently, we established colon cancer cell lines with stable FOXE1 expression to observe the biological effect on colorectal cancer, including cell growth, migration, actin cytoskeleton, and growth of human colorectal xenografts in nude mice. Ectopic expression of FOXE1 could suppress tumor cell growth and migration and affect the organization of the actin cytoskeleton together with suppressing tumorigenicity in vivo. FOXE1 methylation was frequently seen in association with a complete absence of or downregulated gene expression, and FOXE1 plays a suppressive role in the development and progression of colorectal cancer.

Multifactor dimensionality reduction method identifies novel SNP interactions in the WNT protein interaction networks that are associated with recurrence risk in colorectal cancer

Background

Interactions among genetic variants are rarely studied but may explain a part of the variability in patient outcomes.

Objectives

In this study, we aimed to identify 1 to 3 way interactions among SNPs from five Wnt protein interaction networks that predict the 5-year recurrence risk in a cohort of stage I-III colorectal cancer patients.

Methods

423 patients recruited to the Newfoundland Familial Colorectal Cancer Registry were included. Five Wnt family member proteins (Wnt1, Wnt2, Wnt5a, Wnt5b, and Wnt11) were selected. The BioGRID database was used to identify the proteins interacting with each of these proteins. Genotypes of the SNPs located in the interaction network genes were retrieved from a genome-wide SNP genotype data previously obtained in the patient cohort. The GMDR 0.9 program was utilized to examine 1-, 2-, and 3-SNP interactions using a 5-fold cross validation step. Top GMDR 0.9 models were assessed by permutation testing and, if significant, prognostic associations were verified by multivariable logistic regression models.

Results

GMDR 0.9 has identified novel 1, 2, and 3-way SNP interactions associated with 5-year recurrence risk in colorectal cancer. Nine of these interactions were multi loci interactions (2-way or 3-way). Identified interaction models were able to distinguish patients based on their 5-year recurrence-free status in multivariable regression models. The significance of interactions was the highest in the 3-SNP models. Several of the identified SNPs were eQTLs, indicating potential biological roles of the genes they were associated with in colorectal cancer recurrence.

Conclusions

We identified novel interacting genetic variants that associate with 5-year recurrence risk in colorectal cancer. A significant portion of the genes identified were previously linked to colorectal cancer pathogenesis or progression. These variants and genes are of interest for future functional and prognostic studies. Our results provide further evidence for the utility of GMDR models in identifying novel prognostic biomarkers and the biological importance of the Wnt pathways in colorectal cancer.

Examining SNP-SNP interactions and risk of clinical outcomes in colorectal cancer using multifactor dimensionality reduction based methods

Background: SNP interactions may explain the variable outcome risk among colorectal cancer patients. Examining SNP interactions is challenging, especially with large datasets. Multifactor Dimensionality Reduction (MDR)-based programs may address this problem.

Objectives: 1) To compare two MDR-based programs for their utility; and 2) to apply these programs to sets of MMP and VEGF-family gene SNPs in order to examine their interactions in relation to colorectal cancer survival outcomes.

Methods: This study applied two data reduction methods, Cox-MDR and GMDR 0.9, to study one to three way SNP interactions. Both programs were run using a 5-fold cross validation step and the top models were verified by permutation testing. Prognostic associations of the SNP interactions were verified using multivariable regression methods. Eight datasets, including SNPs from MMP family genes (n = 201) and seven sets of VEGF-family interaction networks (n = 1,517 SNPs) were examined.

Results: ∼90 million potential interactions were examined. Analyses in the MMP and VEGF gene family datasets found several novel 1- to 3-way SNP interactions. These interactions were able to distinguish between the patients with different outcome risks (regression p-values 0.03–2.2E-09). The strongest association was detected for a 3-way interaction including CHRM3.rs665159_EPN1.rs6509955_PTGER3.rs1327460 variants.

Conclusion: Our work demonstrates the utility of data reduction methods while identifying potential prognostic markers in colorectal cancer.

A comprehensive analysis of SNPs and CNVs identifies novel markers associated with disease outcomes in colorectal cancer

We aimed to examine the associations of a genome-wide set of single nucleotide polymorphisms (SNPs) and 254 copy number variations (CNVs) and/or insertion/deletions (INDELs) with clinical outcomes in colorectal cancer patients (n = 505). We also aimed to investigate whether their associations changed (e.g., appeared, diminished) over time. Multivariable Cox proportional hazards and piece-wise Cox regression models were used to examine the associations. The Cancer Genome Atlas (TCGA) datasets were used for replication purposes and to examine the gene expression differences between tumor and nontumor tissue samples. A common SNP (WBP11-rs7314075) was associated with disease-specific survival with P-value of 3.2 × 10-8 . Association of this region with disease-specific survival was also detected in the TCGA patient cohort. Two expression quantitative trait loci (eQTLs) were identified in this locus that were implicated in the regulation of ERP27 expression. Interestingly, expression levels of ERP27 and WBP11 were significantly different between colorectal tumors and nontumor tissues. Three SNPs predicted the risk of recurrent disease only after 5 years postdiagnosis. Overall, our study identified novel variants, one of which also showed an association in the TCGA dataset, but no CNVs/INDELs, that associated with outcomes in colorectal cancer. Three SNPs were candidate predictors of long-term recurrence/metastasis risk.

Computer-aided biomarker discovery for precision medicine: data resources, models and applications

Biomarkers are a class of measurable and evaluable indicators with the potential to predict disease initiation and progression. In contrast to disease-associated factors, biomarkers hold the promise to capture the changeable signatures of biological states. With methodological advances, computer-aided biomarker discovery has now become a burgeoning paradigm in the field of biomedical science. In recent years, the 'big data' term has accumulated for the systematical investigation of complex biological phenomena and promoted the flourishing of computational methods for systems-level biomarker screening. Compared with routine wet-lab experiments, bioinformatics approaches are more efficient to decode disease pathogenesis under a holistic framework, which is propitious to identify biomarkers ranging from single molecules to molecular networks for disease diagnosis, prognosis and therapy. In this review, the concept and characteristics of typical biomarker types, e.g. single molecular biomarkers, module/network biomarkers, cross-level biomarkers, etc., are explicated on the guidance of systems biology. Then, publicly available data resources together with some well-constructed biomarker databases and knowledge bases are introduced. Biomarker identification models using mathematical, network and machine learning theories are sequentially discussed. Based on network substructural and functional evidences, a novel bioinformatics model is particularly highlighted for microRNA biomarker discovery. This article aims to give deep insights into the advantages and challenges of current computational approaches for biomarker detection, and to light up the future wisdom toward precision medicine and nation-wide healthcare.

No associations of a set of SNPs in the Vascular Endothelial Growth Factor (VEGF) and Matrix Metalloproteinase (MMP) genes with survival of colorectal cancer patients

In this study, we aimed to investigate the associations of genetic variations within select genes functioning in angiogenesis, lymph‐angiogenesis, and metastasis pathways and the risk of outcome in colorectal cancer patients. We followed a two‐stage analysis: First, 381 polymorphisms from 30 genes (eight Vascular Endothelial Growth Factor (VEGF) and 22 Matrix Metalloproteinase [MMP] genes) were investigated in the discovery cohort (n = 505). Then, 16 polymorphisms with the lowest P‐value in this analysis were investigated in a separate replication cohort (n = 247). Genotypes were obtained using the Illumina^® HumanOmni‐1‐Quad (discovery cohort) and Sequenom MassArray^® (replication cohort) platforms. The primary outcome measure was overall survival (OS). Kaplan–Meier, univariate and multivariable Cox regression methods were used to test the associations between genotypes and OS. Four SNPs (rs12365082, rs11225389, rs11225388, and rs2846707) had the univariate analysis P < 0.05 in both the discovery and replication cohorts. These SNPs are in linkage disequilibrium with each other to varying extent and are located in the MMP8 and MMP27 genes. In the multivariable analysis adjusting for age, stage, and microsatellite instability status, three of these SNPs (rs12365082, rs11225389, rs11225388) were independent predictors of OS (P < 0.05) in the discovery cohort. However, the same analysis in the replication cohort did not yield statistically significant results. Overall, while the genetic variations in the VEGF and MMP genes are attractive candidates as prognostic markers, our study showed no evidence of associations of a large set of SNPs in these genes and overall survival of colorectal cancer patients in our study.

A Survival Association Study of 102 Polymorphisms Previously Associated with Survival Outcomes in Colorectal Cancer

Several published studies identified associations of a number of polymorphisms with a variety of survival outcomes in colorectal cancer. In this study, we aimed to explore 102 previously reported common genetic polymorphisms and their associations with overall survival (OS) and disease-free survival (DFS) in a colorectal cancer patient cohort from Newfoundland (n = 505). Genotypes were obtained using a genomewide SNP genotyping platform. For each polymorphism, the best possible genetic model was estimated for both overall survival and disease-free survival using a previously published approach. These SNPs were then analyzed under their genetic models by Cox regression method. Correction for multiple comparisons was performed by the False Discovery Rate (FDR) method. Univariate analysis results showed that RRM1-rs12806698, IFNGR1-rs1327474, DDX20-rs197412, and PTGS2-rs5275 polymorphisms were nominally associated with OS or DFS (p < 0.01). In stage-adjusted analysis, the nominal associations of DDX20-rs197412, PTGS2-rs5275, and HSPA5-rs391957 with DFS were detected. However, after FDR correction none of these polymorphisms remained significantly associated with the survival outcomes. We conclude that polymorphisms investigated in this study are not associated with OS or DFS in our colorectal cancer patient cohort.

Examining the polymorphisms in the hypoxia pathway genes in relation to outcome in colorectal cancer

Introduction

Colorectal cancer is a common malignancy. Identification of genetic prognostic markers may help prognostic estimations in colorectal cancer. Genes that regulate response to hypoxia and other genes that are regulated under the hypoxic conditions have been shown to play roles in cancer progression. In this study, we hypothesized that genetic variations in the hypoxia pathway genes were associated with the risk of outcome in colorectal cancer patients.

Methods

This study was performed in two phases. In the first phase, 49 SNPs from six hypoxia pathway genes (HIF1A, HIF1B, HIF2A, LOX, MIF and CXCL12) in 272 colorectal cancer patients were analyzed. In the second phase, 77 SNPs from seven hypoxia pathway genes (HIF1A, HIF1B, HIF2A, HIF2B, HIF3A, LOX and CXCL12) were analyzed in an additional cohort of 535 patients. Kaplan Meier, Cox univariate and multivariable regression analyses were performed to analyze the relationship between the SNPs and overall survival (OS), disease free survival (DFS) or disease specific survival (DSS). Since this was a hypothesis-generating study, no correction for multiple testing was applied.

Results

In phase I, one SNP (HIF2A rs11125070) was found to be associated with DFS in multivariable analysis; yet association of a proxy polymorphism (HIF2A rs4953342) was not detected in the phase II patient cohort. In phase II, associations of two SNPs (HIF2A rs4953352 and HIF2B rs12593988) were significant in both OS and DFS multivariable analyses. However, association of HIF2A rs4953352 was not replicated in the phase I cohort using a proxy SNP (HIF2A rs6706003).

Conclusion

Overall, our study did not find a convincing evidence of association of the investigated polymorphisms with the disease outcomes in colorectal cancer.

Mutation and expression analysis of the IDH1, IDH2, DNMT3A, and MYD88 genes in colorectal cancer

Colorectal cancer (CRC) is one of the leading causes of death around the world. Its genetic mechanism was intensively investigated in the past decades with findings of a number of canonical oncogenes and tumor-suppressor genes such as APC, KRAS, and TP53. Recent genome-wide association and sequencing studies have identified a series of promising oncogenes including IDH1, IDH2, DNMT3A, and MYD88 in hematologic malignancies. However, whether these genes are involved in CRC remains unknown. In this study, we screened the hotspot mutations of these four genes in 305 CRC samples from Han Chinese by direct sequencing. mRNA expression levels of these genes were quantified by quantitative real-time PCR (RT-qPCR) in paired cancerous and paracancerous tissues. Association analyses between mRNA expression levels and different cancerous stages were performed. Except for one patient harboring IDH1 mutation p.I99M, we identified no previously reported hotspot mutations in colorectal cancer tissues. mRNA expression levels of IDH1, DNMT3A, and MYD88, but not IDH2, were significantly decreased in the cancerous tissues comparing with the paired paracancerous normal tissues. Taken together, the hotspot mutations of IDH1, IDH2, DNMT3A, and MYD88 gene were absent in CRC. Aberrant mRNA expression of IDH1, DNMT3A, and MYD88 gene might be actively involved in the development of CRC.

Colorectal cancer risk and patients' survival: influence of polymorphisms in genes somatically mutated in colorectal tumors

PURPOSE

The first two studies aiming for the high-throughput identification of the somatic mutation spectrum of colorectal cancer (CRC) tumors were published in 2006 and 2007. Using exome sequencing, they described 69 and 140 candidate cancer genes (CAN genes), respectively. We hypothesized that germline variants in these genes may influence CRC risk, similar to APC, which is causing CRC through germline and somatic mutations.

METHODS

After excluding the well-established CRC genes APC, KRAS, TP53, and ABCA1, we analyzed 35 potentially functional single-nucleotide polymorphisms (SNPs) in 10 CAN genes (OBSCN, MLL3, PKHD1, SYNE1, ERCC6, FBXW7, EPHB6/TRPV6, ELAC1/SMAD4, EPHA3, and ADAMTSL3) using KBiosciences Competitive Allele-Specific PCR™ genotyping assays. In addition to CRC risk (1,399 CRC cases, 838 controls), we also considered the influence of the SNPs on patients' survival (406 cases).

RESULTS

In spite of the fact that our in silico analyses suggested functional relevance for the studied genes and SNPs, our data did not support a strong influence of the studied germline variants on CRC risk and survival. The strongest association with CRC risk and survival was found for MLL3 (rs6464211, OR 1.50, p = 0.002, dominant model; HR 2.12, p = 0.020, recessive model). Two SNPs in EPHB6/TRPV6 (dominant model) showed marginal associations with survival (rs4987622 HR 0.58 p = 0.028 and rs6947538 HR 0.64, p = 0.036, respectively).

CONCLUSION

Although somatic mutations in the CAN genes have been related to the development and progression of various types of cancers in several next-generation sequencing or expression analyses, our study suggests that the studied potentially functional germline variants are not likely to affect CRC risk or survival.

Adult genetic risk screening

Recent advances in genetic analysis especially DNA sequencing technology open a new strategy for adult disease prevention by genetic screening. Physicians presently treat disease pathology with less emphasis on disease risk prevention/reduction. Genetic screening has reduced the incidence of untreatable childhood genetic diseases and improved the care of newborns. The opportunity exists to expand screening programs and reduce the incidence of adult onset diseases via genetic risk identification and disease intervention. This article outlines the approach, challenges, and benefits of such screening for adult genetic disease risks.

Genetic polymorphisms of the CASP8 gene promoter may not be associated with colorectal cancer in Han Chinese from southwest China

PURPOSE

Caspase 8 (CASP8) plays a critical role in the apoptotic pathway and aberrant regulation of this pathway causes many diseases including cancers. Genetic variants rs3834129 (CTTACT/-) and rs3769821 (T/C) in the promoter region of the CASP8 gene were documented to be associated with multiple solid cancers and non-Hodgkin's lymphoma (NHL), respectively, despite of some controversies. We aimed to discern potential association of these two variants and rs113686495 (CTGTCATT/-), as well as CASP8 mRNA and protein expression levels with colorectal cancer (CRC) in Han Chinese.

METHODS

We genotyped CASP8 genetic variants in 305 CRC patients and 342 healthy individuals from Kunming, Southwest China. Expression levels of CASP8 mRNA and protein were quantified in paired cancerous and paracancerous normal tissues by using real-time quantitative PCR and western blot, respectively. We compared the frequencies of alleles, genotypes, and haplotypes between the cases and controls. Correlation of CASP8 mRNA and protein expression levels in paired cancerous and paracancerous normal tissues from patients with different genotypes and clinical expression were also evaluated.

RESULTS

There was no association of the CASP8 genetic variants with CRC in our case-control study. The CASP8 gene mRNA expression levels in cancerous and paracancerous normal tissues were similar and there was no significant difference between subjects with different genotypes and clinical features. However, we found that CASP8 protein level was significantly lower in cancerous tissues than in paired paracancerous normal tissues.

CONCLUSIONS

Our results suggest that the three CASP8 genetic variants may not be associated with CRC risk in Han Chinese from southwest China. Aberrant CASP8 protein expression may play a role in the pathogenesis of CRC.

MTHFR Glu429Ala and ERCC5 His46His polymorphisms are associated with prognosis in colorectal cancer patients: analysis of two independent cohorts from Newfoundland

Introduction

In this study, 27 genetic polymorphisms that were previously reported to be associated with clinical outcomes in colorectal cancer patients were investigated in relation to overall survival (OS) and disease free survival (DFS) in colorectal cancer patients from Newfoundland.

Methods

The discovery and validation cohorts comprised of 532 and 252 patients, respectively. Genotypes of 27 polymorphisms were first obtained in the discovery cohort and survival analyses were performed assuming the co-dominant genetic model. Polymorphisms associated with disease outcomes in the discovery cohort were then investigated in the validation cohort.

Results

When adjusted for sex, age, tumor stage and microsatellite instability (MSI) status, four polymorphisms were independent predictors of OS in the discovery cohort MTHFR Glu429Ala (HR: 1.72, 95%CI: 1.04–2.84, p = 0.036), ERCC5 His46His (HR: 1.78, 95%CI: 1.15–2.76, p = 0.01), SERPINE1 −675indelG (HR: 0.52, 95%CI: 0.32–0.84, p = 0.008), and the homozygous deletion of GSTM1 gene (HR: 1.4, 95%CI: 1.03–1.92, p = 0.033). In the validation cohort, the MTHFR Glu429Ala polymorphism was associated with shorter OS (HR: 1.71, 95%CI: 1.18–2.49, p = 0.005), although with a different genotype than the discovery cohort (CC genotype in the discovery cohort and AC genotype in the validation cohort). When stratified based on treatment with 5-Fluorouracil (5-FU)-based regimens, this polymorphism was associated with reduced OS only in patients not treated with 5-FU. In the DFS analysis, when adjusted for other variables, the TT genotype of the ERCC5 His46His polymorphism was associated with shorter DFS in both cohorts (discovery cohort: HR: 1.54, 95%CI: 1.04–2.29, p = 0.032 and replication cohort: HR: 1.81, 95%CI: 1.11–2.94, p = 0.018).

Conclusions

In this study, associations of the MTHFR Glu429Ala polymorphism with OS and the ERCC5 His46His polymorphism with DFS were identified in two colorectal cancer patient cohorts. Our results also suggest that the MTHFR Glu429Ala polymorphism may be an adverse prognostic marker in patients not treated with 5-FU.

A curated database of genetic markers from the angiogenesis/VEGF pathway and their relation to clinical outcome in human cancers

INTRODUCTION

Angiogenesis causes local growth, aggressiveness and metastasis in solid tumors, and thus, is almost always associated with poor prognosis and survival in cancer patients. Because of this clinical importance, several chemotherapeutic agents targeting angiogenesis have also been developed. Genes and genetic variations in angiogenesis/VEGF pathway thus may be correlated with clinical outcome in cancer patients.

MATERIAL AND METHODS

Here, we describe a manually curated public database, dbANGIO, which posts the results of studies testing the possible correlation of genetic variations (polymorphisms and mutations) from the angiogenesis/VEGF pathway with demographic features, clinicopathological features, treatment response and toxicity, and prognosis and survival-related endpoints in human cancers. The scientific findings are retrieved from PUBMED and posted in the dbANGIO website in a summarized form.

RESULTS AND CONCLUSION

As of September 2011, dbANGIO includes 362 entries from 83 research articles encompassing 154 unique genetic variations from 39 genes investigated in several solid and hematological cancers. By curating the literature findings and making them freely available to researchers, dbANGIO will expedite the research on genetic factors from the angiogenesis pathway and will assist in their utility in clinical management of cancer patients. dbANGIO is freely available for non-profit institutions at http://www.med.mun.ca/angio.

Polymorphisms of caspase 8 and caspase 9 gene and colorectal cancer susceptibility and prognosis

PURPOSE

Caspase-8 (CASP8) and caspase-9 (CASP9) play crucial roles in regulating apoptosis, and their functional polymorphisms may alter cancer risk. Our aim was to investigate the association between CASP8 and CASP9 gene polymorphisms and colorectal cancer (CRC) susceptibility.

METHODS

A case-control study at 402 CRC patients and 480 healthy controls was undertaken in order to investigate the association between the genotype and allelic frequencies of CASP8 -652 6N ins/del and CASP9 -1263 A>G polymorphisms and the CRC susceptibility. The polymerase chain reaction (PCR) restriction fragment length polymorphism method was used and the incidence of polymorphisms on messenger RNA (mRNA) expression levels was detected by quantitative reverse-transcriptase PCR in CRC tissues.

RESULTS

No statistical significant association was observed between CASP8 -652 6N ins/del polymorphism frequencies and CRC susceptibility. CASP9 -1263 G allele was observed to be significant associated with reduced risk of CRC. Homozygotes for the -1263 GG CASP9 genotype, and hetrozygotes for the -1263 AG genotype expressed 6.64- and 3.69-fold higher mRNA levels of Caspase-9, respectively compared to the -1263 AA genotype cases. No significant association was observed between CASP9 -1263 A>G polymorphism and tumor characteristics. The CASP9 -1263 GG genotype was associated with increased overall survival in CRC patients.

CONCLUSION

The CASP9 -1263 A>G polymorphism was observed to play a protective role in CRC predisposition, while the CASP9 -1263 GG genotype may confer a better prognosis at CRC patients.