Frequent Alteration of the Tumor Suppressor Gene APC in Sporadic Canine Colorectal Tumors

Sporadic canine colorectal cancers (CRCs) should make excellent models for studying the corresponding human cancers. To molecularly characterize canine CRC, we investigated exonic sequence mutations of adenomatous polyposis coli (APC), the best known tumor suppressor gene of human CRC, in 23 sporadic canine colorectal tumors, including 8 adenomas and 15 adenocarcinomas, via exon-resequencing analysis. As a comparison, we also performed the same sequencing analysis on 10 other genes, either located at human 5q22 (the same locus as APC) or 18q21 (also frequently altered in human CRC), or known to play a role in human carcinogenesis. We noted that APC was the most significantly mutated gene in both canine adenomas and adenocarcinomas among the 11 genes examined. Significantly, we detected large deletions of ≥10 bases, many clustered near the mutation cluster region, as well as single or two base deletions in ∼70% canine tumors of both subtypes. These observations indicate that like in the human, APC is also frequently altered in sporadic colorectal tumors in the dog and its alteration is an early event in canine colorectal tumorigenesis. Our study provides further evidence demonstrating the molecular similarity in pathogenesis between sporadic human and canine CRCs. This work, along with our previous copy number abnormality study, supports that sporadic canine CRCs are valid models of human CRCs at the molecular level.     

Zyxin as a potential cancer prognostic marker promotes the proliferation and metastasis of colorectal cancer cells

Colorectal cancer (CRC) is a leading cause of cancer death. This study was conducted to investigate the functions and mechanisms of Zyxin (ZYX) in CRC. Multiomics analysis associated ZYX with CRC metastasis. ZYX expression levels were increased in human CRC tissues and related to shorter recurrence-free survival. Knockdown of ZYX expression resulted in inhibition of cell growth, invasion, and migration in vitro and in vivo. Comprehensive analysis of gene microarray analysis showed that ZYX may activate the pathway of NUPR1 and JNK, inhibit CST5, regulate focal adhesion (FA), and affect epithelial–mesenchymal transition in CRC cells. Results of gene microarray and membrane protein isobaric tags with relative and absolute quantitation labeling mass spectrometry found ten differentially expressed genes, which were associated with ZYX activity. Furthermore, real-time polymerase chain reaction was used to validate the expression patterns of selected genes in the integrative analysis. Taken together, our findings provide the first evidence that decreased expression level of ZYX impairs CRC cell proliferation and metastasis probably via the FA pathway.     

A robust immune-related lncRNA signature for the prognosis of human colorectal cancer

Background: Colorectal cancer (CRC) is one of the most prevalent malignant cancers worldwide. Immune-related long non-coding RNAs (IRlncRNAs) are proved to be essential in the development and progression of carcinoma. The purpose of the present study was to develop and validate a prognostic IRlncRNA signature for CRC patients.Methods: Gene expression profiles of CRC samples were downloaded from The Cancer Genome Atlas (TCGA) database. Immune-related genes were obtained from the ImmPort database and were used to identify IRlncRNA by correlation analysis. Through LASSO Cox regression analyses, a prognostic signature was constructed. Functional enrichment analysis was performed by gene set enrichment analysis (GSEA). TIMER2.0 web server and tumor immune dysfunction and exclusion (TIDE) algorithm were employed to analyze the association between our model and tumor-infiltrating immune cells and immunotherapy response. The expression levels of IRlncRNAs in cell lines were detected by quantitative real-time PCR (qPCR).Results: A 9-IRlncRNA signature was developed by a LASSO Cox proportional regression model. Based on the signature, CRC patients were divided into high- and low-risk groups with different prognoses. GSEA results indicated that patients in high-risk group were associated with cancer-related pathways. In addition, patients in low-risk group were found to have more infiltration of anti-tumor immune cells and might show a favorable response to immunotherapy. Finally, the result of qPCR revealed that most IRlncRNAs were differently expressed between normal and tumor cell lines.Conclusion: The constructed 9-IRlncRNA signature has potential to predict the prognosis of CRC patients and may be helpful to guide personalized immunotherapy.     

SVM-T-RFE: a novel gene selection algorithm for identifying metastasis-related genes in colorectal cancer using gene expression profiles

Although metastasis is the principal cause of death cause for colorectal cancer (CRC) patients, the molecular mechanisms underlying CRC metastasis are still not fully understood. In an attempt to identify metastasis-related genes in CRC, we obtained gene expression profiles of 55 early stage primary CRCs, 56 late stage primary CRCs, and 34 metastatic CRCs from the expression project in Oncology (http://www.intgen.org/expo/). We developed a novel gene selection algorithm (SVM-T-RFE), which extends support vector machine recursive feature elimination (SVM-RFE) algorithm by incorporating T-statistic. We achieved highest classification accuracy (100%) with smaller gene subsets (10 and 6, respectively), when classifying between early and late stage primary CRCs, as well as between metastatic CRCs and late stage primary CRCs. We also compared the performance of SVM-T-RFE and SVM-RFE gene selection algorithms on another large-scale CRC dataset and the five public microarray datasets. SVM-T-RFE bestowed SVM-RFE algorithm in identifying more differentially expressed genes, and achieving highest prediction accuracy using equal or smaller number of selected genes. A fraction of selected genes have been reported to be associated with CRC development or metastasis.     

Identification of novel TGF-β regulated genes with pro-migratory roles

Transforming growth factor-β (TGF-β) signaling plays fundamental roles in the development and homeostasis of somatic cells. Dysregulated TGF-β signaling contributes to cancer progression and relapse to therapies by inducing epithelial-to-mesenchymal transition (EMT), enriching cancer stem cells, and promoting immunosuppression. Although many TGF-β-regulated genes have been identified, only a few datasets were obtained by next-generation sequencing. In this study, we performed RNA-sequencing analysis of MCF10A cells and identified 1166 genes that were upregulated and 861 genes that were downregulated by TGF-β. Gene set enrichment analysis revealed that focal adhesion and metabolic pathways were the top enriched pathways of the up- and downregulated genes, respectively. Genes in these pathways also possess significant predictive value for renal cancers. Moreover, we confirmed that TGF-β induced expression of MICAL1 and 2, and the histone demethylase, KDM7A, and revealed their regulatory roles on TGF-β-induced cell migration. We also show a critical effect of KDM7A in regulating the acetylation of H3K27 on TGF-β-induced genes. In sum, this study identified novel effectors that mediate the pro-migratory role of TGF-β signaling, paving the way for future studies that investigate the function of MICAL family members in cancer and the novel epigenetic mechanisms downstream TGF-β signaling.     

Tyrosine phosphatases as a superfamily of tumor suppressors in colorectal cancer

Phosphorylation and dephosphorylation processes catalyzed by numerous kinases and phosphorylases are essential for cell homeostasis and may lead to disturbances in a variety of vital cellular pathways, such as cell proliferation and differentiation, and thus to complex diseases including cancer. As over 80 % of all oncogenes encode protein tyrosine kinases (PTKs), protein tyrosine phosphatases (PTPs), which can reverse the effects of tyrosine kinases, are very important tumor suppressors. Alterations in tyrosine kinase and phosphatase genes including point mutations, changes in epigenetic regulation, as well as chromosomal aberrations involving regions critical to these genes, are frequently observed in a variety of cancers. Colorectal cancer (CRC) is one of the most common cancers in humans. CRCs occur in a familial (about 15 % of all cases), hereditary (about 5%) and sporadic (almost 75-80 %) form. As genetic-environmental interrelations play an important role in the susceptibility to sporadic forms of CRCs, many studies are focused on genetic alterations in such tumors. Mutational analysis of the tyrosine phosphatome in CRCs has identified somatic mutations in PTPRG, PTPRT, PTPN3, PTPN13 and PTPN14. The majority of these mutations result in a loss of protein function. Also, alterations in the expression of these genes, such as decreased expression of PTPRR, PTPRO, PTPRG and PTPRD, mediated by epigenetic mechanisms have been observed in a variety of tumors. Since cancer is a social and global problem, there will be a growing number of studies on alterations in the candidate cancer genes, including protein kinases and phosphatases, to determine the origin, biology and potential pathways for targeted anticancer therapy.     

Genome-wide gene expression analysis suggests an important role of hypoxia in the pathogenesis of endemic osteochondropathy Kashin-Beck disease

Kashin-Beck Disease (KBD) is an endemic osteochondropathy, the pathogenesis of which remains unclear now. In this study, we compared gene expression profiles of articular cartilage derived respectively from KBD patients and normal controls. Total RNA were isolated, amplified, labeled and hybridized to Agilent human 1A 22 k whole genome microarray chip. qRT-PCR was conducted to validate our microarray data. We detected 57 up-regulated genes (ratios ≥2.0) and 24 down-regulated genes (ratios ≤0.5) in KBD cartilage. To further identify the key genes involved in the pathogenesis of KBD, Bayesian analysis of variance for microarrays (BAM) software was applied and identified 12 potential key genes with an average ratio 6.64, involved in apoptosis, metabolism, cytokine & growth factor and cytoskeleton & cell movement. Gene Set Enrichment Analysis (GSEA) software was used to identify differently expressed gene ontology categories and pathways. GSEA found that a set of apoptosis, hypoxia and mitochondrial function related gene ontology categories and pathways were significantly up-regulated in KBD compared to normal controls. Based on the results of this study, we suggest that chronic hypoxia-induced mitochondrial damage and apoptosis might play an important role in the pathogenesis of KBD. Our efforts may help to understand the pathogenesis of KBD as well as other osteoarthrosis with similar articular cartilage lesions.     

PAGE: Parametric Analysis of Gene Set Enrichment

Background  

Gene set enrichment analysis (GSEA) is a microarray data analysis method that uses predefined gene sets and ranks of genes to identify significant biological changes in microarray data sets. GSEA is especially useful when gene expression changes in a given microarray data set is minimal or moderate.     

Network cluster analysis of protein–protein interaction network identified biomarker for early onset colorectal cancer

Colorectal cancer (CRC) is a major cause of morbidity and mortality throughout the world. However, the genetic alterations and molecular mechanism of the early onset CRCs are not fully investigated. The present study aimed to characterize early onset CRC by analyzing its gene expression compared with normal controls and to identify network-based biomarkers of early onset CRC. The gene expression profiles of early onset CRC were downloaded from Gene Expression Omnibus and the differentially expressed genes (DEGs) in CRC patients were identified. Then, a protein–protein interaction (PPI) network was constructed and the clusters in PPI were analyzed by ClusterONE. Furthermore, the gene ontology functional analysis and pathway enrichment analysis were conducted to the modules in PPI network. A systems biology approach integrating microarray data and PPI was further applied to construct a PPI network in CRC. Total 631 DEGs were identified from the early onset CRC compared to healthy controls. These genes were found to be involved in several biological processes, including cell communication, cell proliferation, cell shape and apoptosis. Five functional modules which may play important roles in the initiation of early onset CRC were identified from the PPI network. Functional annotation revealed that these five modules were involved in the pathways of signal transduction, carcinogenesis and metastasis. The hub nodes of these five modules, CDC42, TEX11, QKI, CAV1 and FN1, may serve as the biomarkers of early onset CRC and have the potential to be targets for therapeutic intervention. However, further investigations are still needed to confirm our findings.     

A proteomics analysis of cell signaling alterations in colorectal cancer

To gain further insight into alterations in cellular pathways, tumor profiling, and marker discovery in colorectal cancer (CRC) we used a new antibody microarray specific for cell signaling. Soluble protein extracts were prepared from paired tumor/normal biopsies of 11 patients diagnosed with colorectal carcinoma at different stages; four liver carcinomas were used as a reference. Antibody microarray analysis identified 46 proteins that were differentially expressed between normal colorectal epithelium and adenocarcinoma. These proteins gave a specific signature for CRC, different from other tumors, as well as a panel of novel markers and potential targets for CRC. Twenty-four proteins were validated by using a specific colorectal cancer tissue microarray and immunoblotting analysis. Together with some previously well known deregulated proteins in CRC (beta-catenin, c-MYC, or p63), we found new potential markers preferentially expressed in CRC tumors: cytokeratin 13, calcineurin, CHK1, clathrin light chain, MAPK3, phospho-PTK2/focal adhesion kinase (Ser-910), and MDM2. CHK1 antibodies were particularly effective in discriminating between tumoral and normal mucosa in CRC. Moreover a global picture of alterations in signaling pathways in CRC was observed, including a significant up-regulation of different components of the epidermal growth factor receptor and Wnt/beta-catenin pathways and the down-regulation of p14(ARF). The experimental approach described here should be applicable to other pathologies and neoplastic processes.     

A High Degree of LINE-1 Hypomethylation Is a Unique Feature of Early-Onset Colorectal Cancer

Objective

Early-onset colorectal cancer (CRC) represents a clinically distinct form of CRC that is often associated with a poor prognosis. Methylation levels of genomic repeats such as LINE-1 elements have been recognized as independent factors for increased cancer-related mortality. The methylation status of LINE-1 elements in early-onset CRC has not been analyzed previously.

Design

We analyzed 343 CRC tissues and 32 normal colonic mucosa samples, including 2 independent cohorts of CRC diagnosed ≤50 years old (n = 188), a group of sporadic CRC >50 years (MSS n = 89; MSI n = 46), and a group of Lynch syndrome CRCs (n = 20). Tumor mismatch repair protein expression, microsatellite instability status, LINE-1 and MLH1 methylation, somatic BRAF V600E mutation, and germline MUTYH mutations were evaluated.

Results

Mean LINE-1 methylation levels (±SD) in the five study groups were early-onset CRC, 56.6% (8.6); sporadic MSI, 67.1% (5.5); sporadic MSS, 65.1% (6.3); Lynch syndrome, 66.3% (4.5) and normal mucosa, 76.5% (1.5). Early-onset CRC had significantly lower LINE-1 methylation than any other group (p<0.0001). Compared to patients with <65% LINE-1 methylation in tumors, those with ≥65% LINE-1 methylation had significantly better overall survival (p = 0.026, log rank test).

Conclusions

LINE-1 hypomethylation constitutes a potentially important feature of early-onset CRC, and suggests a distinct molecular subtype. Further studies are needed to assess the potential of LINE-1 methylation status as a prognostic biomarker for young people with CRC.     

PathAct: a novel method for pathway analysis using gene expression profiles

We developed PathAct, a novel method for pathway analysis to investigate the biological and clinical implications of the gene expression profiles. The advantage of PathAct in comparison with the conventional pathway analysis methods is that it can estimate pathway activity levels for individual patient quantitatively in the form of a pathway-by-sample matrix. This matrix can be used for further analysis such as hierarchical clustering and other analysis methods. To evaluate the feasibility of PathAct, comparison with frequently used gene-enrichment analysis methods was conducted using two public microarray datasets. The dataset #1 was that of breast cancer patients, and we investigated pathways associated with triple-negative breast cancer by PathAct, compared with those obtained by gene set enrichment analysis (GSEA). The dataset #2 was another breast cancer dataset with disease-free survival (DFS) of each patient. Contribution by each pathway to prognosis was investigated by our method as well as the Database for Annotation, Visualization and Integrated Discovery (DAVID) analysis. In the dataset #1, four out of the six pathways that satisfied p < 0.05 and FDR < 0.30 by GSEA were also included in those obtained by the PathAct method. For the dataset #2, two pathways (“Cell Cycle” and “DNA replication”) out of four pathways by PathAct were commonly identified by DAVID analysis. Thus, we confirmed a good degree of agreement among PathAct and conventional methods. Moreover, several applications of further statistical analyses such as hierarchical cluster analysis by pathway activity, correlation analysis and survival analysis between pathways were conducted.     

Commensal bacteria, redox stress, and colorectal cancer: mechanisms and models

The potential role for commensal bacteria in colorectal carcinogenesis is explored in this review. Most colorectal cancers (CRCs) occur sporadically and arise from the gradual accumulation of mutations in genes regulating cell growth and DNA repair. Genetic mutations followed by clonal selection result in the transformation of normal cells into malignant derivatives. Numerous toxicological effects of colonic bacteria have been reported. However, those recognized as damaging epithelial cell DNA are most easily reconciled with the currently understood genetic basis for sporadic CRC. Thus, we focus on mechanisms by which particular commensal bacteria may convert dietary procarcinogens into DNA damaging agents (e.g., ethanol and heterocyclic amines) or directly generate carcinogens (e.g., fecapentaenes). Although these and other metabolic activities have yet to be linked directly to sporadic CRC, several lines of investigation are reviewed to highlight difficulties and progress in the area. Particular focus is given to commensal bacteria that alter the epithelial redox environment, such as production of oxygen radicals by Enterococcus faecalis or production of hydrogen sulfide by sulfate-reducing bacteria (SRB). Super-oxide-producing E. faecalis has conclusively been shown to cause colonic epithelial cell DNA damage. Though SRB-derived hydrogen sulfide (H(2)S) has not been reported thus far to induce DNA damage or function as a carcinogen, recent data demonstrate that this reductant activates molecular pathways implicated in CRC. These observations combined with evidence that SRB carriage may be genetically encoded evoke a working model that incorporates multifactorial gene-environment interactions that appear to underlie the development of sporadic CRC.     

Comprehensive Suppression of All Apoptosis-Induced Proliferation Pathways as a Proposed Approach to Colorectal Cancer Prevention and Therapy

Mutations in the WNT/beta-catenin pathway are present in the majority of all sporadic colorectal cancers (CRCs), and histone deacetylase inhibitors induce apoptosis in CRC cells with such mutations. This apoptosis is counteracted by (1) the signaling heterogeneity of CRC cell populations, and (2) the survival pathways induced by mitogens secreted from apoptotic cells. The phenomena of signaling heterogeneity and apoptosis-induced survival constitute the immediate mechanisms of resistance to histone deacetylase inhibitors, and probably other chemotherapeutic agents. We explored the strategy of augmenting CRC cell death by inhibiting all survival pathways induced by the pro-apoptotic agent LBH589, a histone deacetylase inhibitor: AKT, JAK/STAT, and ERK signaling. The apoptosis-enhancing ability of a cocktail of synthetic inhibitors of proliferation was compared to the effects of the natural product propolis. We utilized colorectal adenoma, drug-sensitive and drug-resistant colorectal carcinoma cells to evaluate the apoptotic potential of the combination treatments. The results suggest that an effective approach to CRC combination therapy is to combine apoptosis-inducing drugs (e.g., histone deacetylase inhibitors, such as LBH589) with agents that suppress all compensatory survival pathways induced during apoptosis (such as the cocktail of inhibitors of apoptosis-associated proliferation). The same paradigm can be applied to a CRC prevention approach, as the apoptotic effect of butyrate, a diet-derived histone deacetylase inhibitor, is augmented by other dietary agents that modulate survival pathways (e.g., propolis and coffee extract). Thus, dietary supplements composed by fermentable fiber, propolis, and coffee extract may effectively counteract neoplastic growth in the colon.     

Molecular genetic analysis of 103 sporadic colorectal tumours in Czech patients

The Czech Republic has one of the highest incidences of colorectal cancer (CRC) in Europe. To evaluate whether sporadic CRCs in Czech patients have specific mutational profiles we analysed somatic genetic changes in known CRC genes (APC, KRAS, TP53, CTNNB1, MUTYH and BRAF, loss of heterozygosity (LOH) at the APC locus, microsatellite instability (MSI), and methylation of the MLH1 promoter) in 103 tumours from 102 individuals. The most frequently mutated gene was APC (68.9% of tumours), followed by KRAS (31.1%), TP53 (27.2%), BRAF (8.7%) and CTNNB1 (1.9%). Heterozygous germline MUTYH mutations in 2 patients were unlikely to contribute to the development of their CRCs. LOH at the APC locus was found in 34.3% of tumours, MSI in 24.3% and MLH1 methylation in 12.7%. Seven tumours (6.9%) were without any changes in the genes tested. The analysis yielded several findings possibly specific for the Czech cohort. Somatic APC mutations did not cluster in the mutation cluster region (MCR). Tumours with MSI but no MLH1 methylation showed earlier onset and more severe mutational profiles compared to MSI tumours with MLH1 methylation. TP53 mutations were predominantly located outside the hot spots, and transitions were underrepresented. Our analysis supports the observation that germline MUTYH mutations are rare in Czech individuals with sporadic CRCs. Our findings suggest the influence of specific ethnic genetic factors and/or lifestyle and dietary habits typical for the Czech population on the development of these cancers.     

肿瘤表达谱基因芯片数据的混合效应模型分析

目的:研究混合效应模型(Mixed Effects Model)在肿瘤表达谱基因芯片数据分析中的检验效能,并探讨其分析效果。方法:采用混合效应模型分析肿瘤实例基因芯片数据,并以基因集富集分析方法(GSEA)作为参照比较分析结果的有效性和科学性,探讨其检验效果。结果:通过混合效应模型和基因集富集分析(GSEA)两种方法对肿瘤基因芯片数据的分析和比较,两种方法筛选出共同的差异表达通路外,混合效应模型额外地筛选出来GSEA未能检验到的8条差异表达通路,且得到文献支持;混和效应模型筛选出的前10个差异表达通路中有6个已有生物学证明而基因集富集分析方法(GSEA)筛选出的前10个差异表达通路中仅有4个已有生物学证明。结论:混合效应模型作为top-down方法中的典型代表,其优势在于通过构建潜变量达到降维目的,可有效地减少多个复杂的变异来源从而保证了结果的准确性和科学性,其检验效能优于基因集富集分析方法(GSEA),是一种行之有效的筛选肿瘤基因芯片数据的分析方法。     

NCAM-induced focal adhesion assembly: a functional switch upon loss of E-cadherin

Loss of expression of the cell-cell adhesion molecule E-cadherin is a hallmark of epithelial-mesenchymal transition (EMT) in development and in the progression from epithelial tumours to invasive and metastatic cancers. Here, we demonstrate that the loss of E-cadherin function upregulates expression of the neuronal cell adhesion molecule (NCAM). Subsequently, a subset of NCAM translocates from fibroblast growth factor receptor (FGFR) complexes outside lipid rafts into lipid rafts where it stimulates the non-receptor tyrosine kinase p59(Fyn) leading to the phosphorylation and activation of focal adhesion kinase and the assembly of integrin-mediated focal adhesions. Ablation of NCAM expression during EMT inhibits focal adhesion assembly, cell spreading and EMT. Conversely, forced expression of NCAM induces epithelial cell delamination and migration, and high NCAM expression correlates with tumour invasion. These results establish a mechanistic link between the loss of E-cadherin expression, NCAM function, focal adhesion assembly and cell migration and invasion.