Transcriptome profile of human colorectal adenomas

Genotypes of cancer stem cells characterized by epithelial-to-mesenchymal transition and proliferation related functions

Cancer stem cells (CSCs), or cancer cells with stem cell-like properties, generally exhibit drug resistance and have highly potent cancer inducing capabilities. Genome-wide expression data collected at public repositories over the last few years provide excellent material for studies that can lead to insights concerning the molecular and functional characteristics of CSCs. Here, we conducted functional genomic studies of CSC based on fourteen PCA-screened high quality public CSC whole genome gene expression datasets and, as control, four high quality non-stem-like cancer cell and non-cancerous stem cell datasets from the Gene Expression Omnibus database. A total of 6,002 molecular signatures were taken from the Molecular Signatures Database and used to characterize the datasets, which, under two-way hierarchical clustering, formed three genotypes. Type 1, consisting of mainly glia CSCs, had significantly enhanced proliferation, and significantly suppressed epithelial-mesenchymal transition (EMT), related functions. Type 2, mainly breast CSCs, had significantly enhanced EMT, but not proliferation, related functions. Type 3, composed of ovarian, prostate, and colon CSCs, had significantly suppressed proliferation related functions and mixed expressions on EMT related functions.

Identification of Aging-Associated Gene Expression Signatures That Precede Intestinal Tumorigenesis

Aging-associated alterations of cellular functions have been implicated in various disorders including cancers. Due to difficulties in identifying aging cells in living tissues, most studies have focused on aging-associated changes in whole tissues or certain cell pools. Thus, it remains unclear what kinds of alterations accumulate in each cell during aging. While analyzing several mouse lines expressing fluorescent proteins (FPs), we found that expression of FPs is gradually silenced in the intestinal epithelium during aging in units of single crypt composed of clonal stem cell progeny. The cells with low FP expression retained the wild-type Apc allele and the tissues composed of them did not exhibit any histological abnormality. Notably, the silencing of FPs was also observed in intestinal adenomas and the surrounding normal mucosae of Apc-mutant mice, and mediated by DNA methylation of the upstream promoter. Our genome-wide analysis then showed that the silencing of FPs reflects specific gene expression alterations during aging, and that these alterations occur in not only mouse adenomas but also human sporadic and hereditary (familial adenomatous polyposis) adenomas. Importantly, pharmacological inhibition of DNA methylation, which suppresses adenoma development in Apc-mutant mice, reverted the aging-associated silencing of FPs and gene expression alterations. These results identify aging-associated gene expression signatures that are heterogeneously induced by DNA methylation and precede intestinal tumorigenesis triggered by Apc inactivation, and suggest that pharmacological inhibition of the signature genes could be a novel strategy for the prevention and treatment of intestinal tumors.

Exportin-5 Functions as an Oncogene and a Potential Therapeutic Target in Colorectal Cancer

Purpose: Dysregulated expression of miRNAs has emerged as a hallmark feature in human cancers. Exportin-5 (XPO5), a karyopherin family member, is a key protein responsible for transporting precursor miRNAs from the nucleus to the cytoplasm. Although XPO5 is one of the key regulators of miRNA biogenesis, its functional role and potential clinical significance in colorectal cancer remains unclear.Experimental Design: The expression levels of XPO5 were initially assessed in three genomic datasets, followed by determination and validation of the relationship between XPO5 expression and clinicopathologic features in two independent colorectal cancer patient cohorts. A functional characterization of XPO5 in colorectal cancer was examined by targeted gene silencing in colorectal cancer cell lines and a xenograft animal model.Results: XPO5 is upregulated, both at mRNA and protein levels, in colorectal cancers compared with normal tissues. High XPO5 expression is associated with worse clinicopathologic features and poor survival in colorectal cancer patient cohorts. The siRNA knockdown of XPO5 resulted in reduced cellular proliferation, attenuated invasion, induction of G₁-S cell-cycle arrest, and downregulation of key oncogenic miRNAs in colorectal cancer cells. These findings were confirmed in a xenograft animal model, wherein silencing of XPO5 resulted in the attenuation of tumor growth.Conclusions: XPO5 acts like an oncogene in colorectal cancer by regulating the expression of miRNAs and may be a potential therapeutic target in colorectal cancer. Clin Cancer Res; 23(5); 1312-22. ©2016 AACR.

Induction of KIAA1199/CEMIP is associated with colon cancer phenotype and poor patient survival

Genes induced in colon cancer provide novel candidate biomarkers of tumor phenotype and aggressiveness. We originally identified KIAA1199 (now officially called CEMIP) as a transcript highly induced in colon cancer: initially designating the transcript as Colon Cancer Secreted Protein 1. We molecularly characterized CEMIP expression both at the mRNA and protein level and found it is a secreted protein induced an average of 54-fold in colon cancer. Knockout of CEMIPreduced the ability of human colon cancer cells to form xenograft tumors in athymic mice. Tumors that did grow had increased deposition of hyaluronan, linking CEMIP participation in hyaluronan degradation to the modulation of tumor phenotype. We find CEMIP mRNA overexpression correlates with poorer patient survival. In stage III only (n = 31) or in combined stage II plus stage III colon cancer cases (n = 73), 5-year overall survival was significantly better (p = 0.004 and p = 0.0003, respectively) among patients with low CEMIP expressing tumors than those with high CEMIP expressing tumors. These results demonstrate that CEMIP directly facilitates colon tumor growth, and high CEMIP expression correlates with poor outcome in stage III and in stages II+III combined cohorts. We present CEMIP as a candidate prognostic marker for colon cancer and a potential therapeutic target.

Knockdown of CSE1L Gene in Colorectal Cancer Reduces Tumorigenesis in Vitro

Human cellular apoptosis susceptibility (chromosomal segregation 1-like, CSE1L) gene plays a role in nuclear-to-cytoplasm transport and chromosome segregation during mitosis, cellular proliferation, and apoptosis. CSE1L is involved in colon carcinogenesis. CSE1L gene expression was assessed with three data sets using Affymetrix U133 + gene chips on normal human colonic mucosa (NR), adenomas (ADs), and colorectal carcinoma (CRC). CSE1L protein expression in CRC, AD, and NR from the same patients was measured by immunohistochemistry using a tissue microarray. We evaluated CSE1L expression in CRC cells (HCT116, SW480, and HT29) and its biological functions. CSE1L mRNA was significantly increased in all AD and CRC compared with NR (P < 0.001 and P = 0.02, respectivly). We observed a change in CSE1L staining intensity and cellular localization by immunohistochemistry. CSE1L was significantly increased during the transition from AD to CRC when compared with NR in a CRC tissue microarray (P = 0.01 and P < 0.001). HCT116, SW480, and HT29 cells also expressed CSE1L protein. CSE1L knockdown by shRNA inhibited protein, resulting in decreased cell proliferation, reduced colony formation in soft agar, and induction of apoptosis. CSE1L protein is expressed early and across all stages of CRC development. shRNA knockdown of CSE1L was associated with inhibition of tumorigenesis in CRC cells. CSE1L may represent a potential target for treatment of CRC.

Exosomes in colorectal carcinoma formation: ALIX under the magnifying glass

Exosomes are small membrane vesicles that have important roles in transporting a great variety of bioactive molecules between epithelial compartment and their microenvironment during tumor formation including colorectal adenoma-carcinoma sequence. We tested the mRNA expression of the top 25 exosome-related markers based on ExoCharta database in healthy (n=49), adenoma (n=49) and colorectal carcinoma (n=49) patients using Affymetrix HGU133 Plus2.0 microarrays. Most related genes showed significantly elevated expression including PGK1, PKM, ANXA5, ENO1, HSP90AB1 and MSN during adenoma-carcinoma sequence. Surprisingly, the expression of ALIX (ALG 2-interacting protein X), involved in multivesicular body (MVB) and exosome formation, was significantly reduced in normal vs adenoma (P=5.02 × 10(-13)) and in normal vs colorectal carcinoma comparisons (P=1.51 × 10(-10)). ALIX also showed significant reduction (P<0.05) at the in situ protein level in the epithelial compartment of adenoma (n=35) and colorectal carcinoma (n=37) patients compared with 27 healthy individuals. Furthermore, significantly reduced ALIX protein levels were accompanied by their gradual transition from diffuse cytoplasmic expression to granular signals, which fell into the 0.6-2 μm diameter size range of MVBs. These ALIX-positive particles were seen in the tumor nests, including tumor-stroma border, which suggest their exosome function. MVB-like structures were also detected in tumor microenvironment including α-smooth muscle actin-positive stromal cells, budding off cancer cells in the tumor front as well as in cancer cells entrapped within lymphoid vessels. In conclusion, we determined the top aberrantly expressed exosome-associated markers and revealed the transition of diffuse ALIX protein signals into a MVB-like pattern during adenoma-carcinoma sequence. These tumor-associated particles seen both in the carcinoma and the surrounding microenvironment can potentially mediate epithelial-stromal interactions involved in the regulation of tumor growth, metastatic invasion and therapy response.

A genome-wide assessment of rare copy number variants in colorectal cancer

Colorectal cancer (CRC) is a complex disease with an estimated heritability of approximately 35%. However, known CRC-related common single nucleotide polymorphisms (SNPs) can only explain ~0.65% of the heritability. This “missing heritability” may be explained partially by rare copy number variants (CNVs). In this study, we performed a genome-wide scan using Illumina Human-Omni Express BeadChip, 694 sporadic CRC cases and 1641 controls were eventually included in our analysis after quality control. The global burden analysis revealed a 1.53-fold excess of rare CNVs in CRC cases compared with controls (P < 1 × 10⁻⁶), and the difference being more pronounced for genic rare CNVs and CNVs overlapped with coding regions (1.65-fold and 1.84-fold, respectively, both P < 1 × 10⁻⁶). Interestingly, both the cases in the lowest and middle tertile of age carried a higher burden of rare CNVs comparing to the highest tertile. Furthermore, 639 CNV-disrupted genes exclusive to CRC cases were found to be significantly enriched in gene ontology (GO) terms concerning nucleosome assembly and olfactory receptor activity. Our study was the first to evaluate the burden of rare CNVs in sporadic CRC and suggested that rare CNVs contributed to the missing heritability of CRC.

Feasibility of Unbiased RNA Profiling of Colorectal Tumors: A Proof of Principle

Despite recent advances in molecular profiling of colorectal cancer (CRC), as of yet this has not translated into an unbiased molecular liquid biopsy profile which can accurately screen for early CRC. In this study we depict the profile of early stage CRC as well as for advanced adenomas (AA) by combination of current molecular knowledge with microarray technology, using efficient circulating free plasma RNA purification from blood and RNA amplification technologies. We joined literature search with Affymetrix gene chip experimental procedure to draw the circulating free plasma RNA profile of colorectal cancer disease reflected in blood. The RNA panel was tested by two datasets comparing patients with CRC with healthy subjects and patients with AA to healthy subjects. For the CRC patient cohort (28 CRC cases vs. 41 healthy controls), the ROC analysis of the selected biomarker panel generated a sensitivity of 75% and a specificity of 93% for the detection of CRC using 8-gene classification model. For the AA patient cohort (28 subjects vs. 46 healthy controls), a sensitivity of 60% and a specificity of 87% were calculated using a 2-gene classification model. We have identified a panel of 8 plasma RNA markers as a preliminary panel for CRC detection and subset markers suitable for AA detection. Subjected to extensive clinical validation we suggest that this panel represents a feasible approach and a potential strategy for noninvasive early diagnosis, as a first-line screening test for asymptomatic, average-risk population before colonoscopy.

Neuromedin U is upregulated by Snail at early stages of EMT in HT29 colon cancer cells

BACKGROUND

The epithelial-mesenchymal transition (EMT) is considered a core process that facilitates the escape of cancer cells from the primary tumor site. The transcription factor Snail was identified as a key regulator of EMT; however, the cascade of regulatory events leading to metastasis remains unknown and new predictive markers of the process are awaited.

METHODS

Gene expressions were analysed using real-time PCR, protein level by Western immunoblotting and confocal imaging. The motility of the cells was examined using time-lapse microscopy. Affymetrix GeneChip Human Genome U133 Plus 2.0 analysis was performed to identify transcriptomic changes upon Snail. Snail silencing was performed using siRNA nucleofection. NMU detection was performed by ELISA.

RESULTS

HT29 cells overexpressing Snail showed changed morphology, functions and transcriptomic profile indicating EMT induction. Changes in expression of 324 genes previously correlated with cell motility were observed. Neuromedin U was the second highest upregulated gene in HT29-Snail cells. This increase was validated by real-time PCR. Additionally elevated NMU protein was detected by ELISA in cell media.

CONCLUSIONS

These results show that Snail in HT29 cells regulates early phenotype conversion towards an intermediate epithelial state. We provided the first evidence that neuromedin U is associated with Snail regulatory function of metastatic induction in colon cancer cells.

GENERAL SIGNIFICANCE

We described the global, early transcriptomic changes induced through Snail in HT29 colon cancer cells and suggested NMU involvement in this process.

Hypoxia promotes colon cancer dissemination through up-regulation of cell migration-inducing protein (CEMIP)

Hypoxic stress drives cancer progression by causing a transcriptional reprogramming. Recently, KIAA1199 was discovered to be a cell-migration inducing protein (renamed CEMIP) that is upregulated in human cancers. However, the mechanism of induction of CEMIP in cancer was hitherto unknown. Here we demonstrate that hypoxia induces CEMIP expression leading to enhanced cell migration. Immunohistochemistry of human colon cancer tissues revealed that CEMIP is upregulated in cancer cells located at the invasive front or in the submucosa. CEMIP localization inversely correlated with E-cadherin expression, which is characteristic of the epithelial-to-mesenchymal transition. Mechanistically, hypoxia-inducible-factor-2α (HIF-2α), but not HIF-1α binds directly to the hypoxia response element within the CEMIP promoter region resulting in increased CEMIP expression. Functional characterization reveals that CEMIP is a downstream effector of HIF-2α-mediated cell migration. Expression of CEMIP was demonstrated to negatively correlate with the expression of Jarid1A, a histone demethylase that removes methyl groups from H3K4me3 (an activation marker for transcription), resulting in altered gene repression. Low oxygen tension inhibits the function of Jarid1A, leading to increased presence of H3K4me3 within the CEMIP promoter. These results provide insight into the upregulation of CEMIP within cancer and can lead to novel treatment strategies targeting this cancer cell migration-promoting gene.

IgG Glycome in Colorectal Cancer

PURPOSE

Alternative glycosylation has significant structural and functional consequences on IgG and consequently also on cancer immunosurveillance. Because of technological limitations, the effects of highly heritable individual variations and the differences in the dynamics of changes in IgG glycosylation on colorectal cancer were never investigated before.

EXPERIMENTAL DESIGN

Using recently developed high-throughput UPLC technology for IgG glycosylation analysis, we analyzed IgG glycome composition in 760 patients with colorectal cancer and 538 matching controls. Effects of surgery were evaluated in 28 patients sampled before and three times after surgery. A predictive model was built using regularized logistic regression and evaluated using a 10-cross validation procedure. Furthermore, IgG glycome composition was analyzed in 39 plasma samples collected before initial diagnosis of colorectal cancer.

RESULTS

We have found that colorectal cancer associates with decrease in IgG galactosylation, IgG sialylation and increase in core-fucosylation of neutral glycans with concurrent decrease of core-fucosylation of sialylated glycans. Although a model based on age and sex did not show discriminative power (AUC = 0.499), the addition of glycan variables into the model considerably increased the discriminative power of the model (AUC = 0.755). However, none of these differences were significant in the small set of samples collected before the initial diagnosis.

CONCLUSIONS

Considering the functional relevance of IgG glycosylation for both tumor immunosurveillance and clinical efficacy of therapy with mAbs, individual variation in IgG glycosylation may turn out to be important for prediction of disease course or the choice of therapy, thus warranting further, more detailed studies of IgG glycosylation in colorectal cancer. Clin Cancer Res; 22(12); 3078-86. ©2016 AACR.

Distinct Expression Pattern of a Deafness Gene, KIAA1199, in a Primate Cochlea

Deafness is one of the most common types of congenital impairments, and at least half of the cases are caused by hereditary mutations. Mutations of the gene KIAA1199 are associated with progressive hearing loss. Its expression is abundant in human cochlea, but interestingly the spatial expression patterns are different between mouse and rat cochleae; the pattern in humans has not been fully investigated. We performed immunohistochemical analysis of a nonhuman primate, common marmoset (Callithrix jacchus), cochlea with a KIAA1199-specific antibody. In the common marmoset cochlea, KIAA1199 protein expression was more widespread than in rodents, with all epithelial cells, including hair cells, expressing KIAA1199. Our results suggest that the primate pattern of KIAA1199 expression is wider in comparison with rodents and may play an essential role in the maintenance of cochlear epithelial cells.

Aberrant DNA methylation of WNT pathway genes in the development and progression of CIMP-negative colorectal cancer

The WNT signaling pathway has an essential role in colorectal carcinogenesis and progression, which involves a cascade of genetic and epigenetic changes. We aimed to analyze DNA methylation affecting the WNT pathway genes in colorectal carcinogenesis in promoter and gene body regions using whole methylome analysis in 9 colorectal cancer, 15 adenoma, and 6 normal tumor adjacent tissue (NAT) samples by methyl capture sequencing. Functional methylation was confirmed on 5-aza-2'-deoxycytidine-treated colorectal cancer cell line datasets. In parallel with the DNA methylation analysis, mutations of WNT pathway genes (APC, β-catenin/CTNNB1) were analyzed by 454 sequencing on GS Junior platform. Most differentially methylated CpG sites were localized in gene body regions (95% of WNT pathway genes). In the promoter regions, 33 of the 160 analyzed WNT pathway genes were differentially methylated in colorectal cancer vs. normal, including hypermethylated AXIN2, CHP1, PRICKLE1, SFRP1, SFRP2, SOX17, and hypomethylated CACYBP, CTNNB1, MYC; 44 genes in adenoma vs. NAT; and 41 genes in colorectal cancer vs. adenoma comparisons. Hypermethylation of AXIN2, DKK1, VANGL1, and WNT5A gene promoters was higher, while those of SOX17, PRICKLE1, DAAM2, and MYC was lower in colon carcinoma compared to adenoma. Inverse correlation between expression and methylation was confirmed in 23 genes, including APC, CHP1, PRICKLE1, PSEN1, and SFRP1. Differential methylation affected both canonical and noncanonical WNT pathway genes in colorectal normal-adenoma-carcinoma sequence. Aberrant DNA methylation appears already in adenomas as an early event of colorectal carcinogenesis.

Low molecular weight protein tyrosine phosphatase (LMWPTP) upregulation mediates malignant potential in colorectal cancer

Phosphatases have long been regarded as tumor suppressors, however there is emerging evidence for a tumor initiating role for some phosphatases in several forms of cancer. Low Molecular Weight Protein Tyrosine Phosphatase (LMWPTP; acid phosphatase 1 [ACP1]) is an 18 kDa enzyme that influences the phosphorylation of signaling pathway mediators involved in cancer and is thus postulated to be a tumor-promoting enzyme, but neither unequivocal clinical evidence nor convincing mechanistic actions for a role of LMWPTP have been identified. In the present study, we show that LMWPTP expression is not only significantly increased in colorectal cancer (CRC), but also follows a step-wise increase in different levels of dysplasia. Chemical inhibition of LMWPTP significantly reduces CRC growth. Furthermore, downregulation of LMWPTP in CRC leads to a reduced migration ability in both 2D- and 3D-migration assays, and sensitizes tumor cells to the chemotherapeutic agent 5-FU. In conclusion, this study shows that LMWPTP is not only overexpressed in colorectal cancer, but it is correlated with the malignant potential of this cancer, suggesting that this phosphatase may act as a predictive biomaker of CRC stage and represents a rational novel target in the treatment of this disease.

SPATIAL: A System-level PAThway Impact AnaLysis approach

The goal of pathway analysis is to identify the pathways that are significantly impacted when a biological system is perturbed, e.g. by a disease or drug. Current methods treat pathways as independent entities. However, many signals are constantly sent from one pathway to another, essentially linking all pathways into a global, system-wide complex. In this work, we propose a set of three pathway analysis methods based on the impact analysis, that performs a system-level analysis by considering all signals between pathways, as well as their overlaps. Briefly, the global system is modeled in two ways: (i) considering the inter-pathway interaction exchange for each individual pathways, and (ii) combining all individual pathways to form a global, system-wide graph. The third analysis method is a hybrid of these two models. The new methods were compared with DAVID, GSEA, GSA, PathNet, Crosstalk and SPIA on 23 GEO data sets involving 19 tissues investigated in 12 conditions. The results show that both the ranking and the P-values of the target pathways are substantially improved when the analysis considers the system-wide dependencies and interactions between pathways.

Repositioning of a cyclin-dependent kinase inhibitor GW8510 as a ribonucleotide reductase M2 inhibitor to treat human colorectal cancer

Colorectal cancer (CRC) is the second leading cause of cancer-related death in males and females in the world. It is of immediate importance to develop novel therapeutics. Human ribonucleotide reductase (RRM1/RRM2) has an essential role in converting ribonucleoside diphosphate to 2'-deoxyribonucleoside diphosphate to maintain the homeostasis of nucleotide pools. RRM2 is a prognostic biomarker and predicts poor survival of CRC. In addition, increased RRM2 activity is associated with malignant transformation and tumor cell growth. Bioinformatics analyses show that RRM2 was overexpressed in CRC and might be an attractive target for treating CRC. Therefore, we attempted to search novel RRM2 inhibitors by using a gene expression signature-based approach, connectivity MAP (CMAP). The result predicted GW8510, a cyclin-dependent kinase inhibitor, as a potential RRM2 inhibitor. Western blot analysis indicated that GW8510 inhibited RRM2 expression through promoting its proteasomal degradation. In addition, GW8510 induced autophagic cell death. In addition, the sensitivities of CRC cells to GW8510 were associated with the levels of RRM2 and endogenous autophagic flux. Taken together, our study indicates that GW8510 could be a potential anti-CRC agent through targeting RRM2.

Downregulation of myosin VI reduced cell growth and increased apoptosis in human colorectal cancer

Colorectal cancer (CRC) is the third most commonly diagnosed cancer worldwide, with the mortality increasing steadily over the last decade. Myosin VI (MYO6) expression is found to be elevated in some types of human carcinoma cell types, suggesting that it may be a sensitive biomarker for the diagnosis and follow-up. In this study, we first used the Oncomine database to explore the expression of MYO6 in CRC tissues, and then constructed a plasmid of RNA interference targeting MYO6 gene. After transfection of lentivirus targeting MYO6 into SW1116 cells, cell viability and proliferation were measured with 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide and colony formation assay. Cell cycle distribution was assayed by flow cytometry and apoptosis was evaluated by Annexin V. MYO6 expression was detected by quantitative real-time polymerase chain reaction and western blot analysis. It was found that MYO6 mRNA was upregulated in CRC tissues using data mining of public Oncomine microarray datasets. Depletion of MYO6 significantly inhibited cell proliferation and colony formation. In addition, knockdown of MYO6 slightly arrested cell cycle in G0/G1 phase, but remarkably increased the proportion of the sub-G1 phase of cell with the increase of apoptotic cells. These results suggest that MYO6 may promote cell growth and may be used as a potential target for anticancer therapy of CRC.

Gene-set activity toolbox (GAT): A platform for microarray-based cancer diagnosis using an integrative gene-set analysis approach

Cancer is a complex disease that cannot be diagnosed reliably using only single gene expression analysis. Using gene-set analysis on high throughput gene expression profiling controlled by various environmental factors is a commonly adopted technique used by the cancer research community. This work develops a comprehensive gene expression analysis tool (gene-set activity toolbox: (GAT)) that is implemented with data retriever, traditional data pre-processing, several gene-set analysis methods, network visualization and data mining tools. The gene-set analysis methods are used to identify subsets of phenotype-relevant genes that will be used to build a classification model. To evaluate GAT performance, we performed a cross-dataset validation study on three common cancers namely colorectal, breast and lung cancers. The results show that GAT can be used to build a reasonable disease diagnostic model and the predicted markers have biological relevance. GAT can be accessed from http://gat.sit.kmutt.ac.th where GAT's java library for gene-set analysis, simple classification and a database with three cancer benchmark datasets can be downloaded.

SNPs in transporter and metabolizing genes as predictive markers for oxaliplatin treatment in colorectal cancer patients

Oxaliplatin is frequently used as part of a chemotherapeutic regimen with 5-fluorouracil in the treatment of colorectal cancer (CRC). The cellular availability of oxaliplatin is dependent on metabolic and transporter enzymes. Variants in genes encoding these enzymes may cause variation in response to oxaliplatin and could be potential predictive markers. Therefore, we used a two-step procedure to comprehensively investigate 1,444 single nucleotide polymorphisms (SNPs) from these pathways for their potential as predictive markers for oxaliplatin treatment, using 623 stage II-IV CRC patients (of whom 201 patients received oxaliplatin) from a German prospective patient cohort treated with adjuvant or palliative chemotherapy. First, all genes were screened using the global test that evaluated SNP*oxaliplatin interaction terms per gene. Second, one model was created by backward elimination on all SNP*oxaliplatin interactions of the selected genes. The statistical procedure was evaluated using bootstrap analyses. Nine genes differentially associated with overall survival according to oxaliplatin treatment (unadjusted p values < 0.05) were selected. Model selection resulted in the inclusion of 14 SNPs from eight genes (six transporter genes, ABCA9, ABCB11, ABCC10, ATP1A1, ATP1B2, ATP8B3, and two metabolism genes GSTM5, GRHPR), which significantly improved model fit. Using bootstrap analysis we show an improvement of the prediction error of 3.7% in patients treated with oxaliplatin. Several variants in genes involved in metabolism and transport could thus be potential predictive markers for oxaliplatin treatment in CRC patients. If confirmed, inclusion of these variants in a predictive test could identify patients who are more likely to benefit from treatment with oxaliplatin.

Overexpression of GRO-β is associated with an unfavorable outcome in colorectal cancer

Growth-related oncogene (GRO)-β, or chemokine (C-X-C motif) ligand 2 (CXCL2), is a member of the CXC chemokine family that may mediate various functions, including attracting neutrophils to sites of inflammation, and participating in tumorigenesis and progression. However, the expression of GRO-β in colorectal cancer (CRC) and the association with the clinical outcome of the disease remains poorly understood. In the present study, CXCL2 mRNA expression in CRC was analyzed using six independent datasets from the Oncomine microarray database. The immunohistochemical analysis of tissue microarrays (TMA) was used to characterize the expression of the GRO-β protein in CRC. The association between GRO-β expression and the clinicopathological features and prognosis of patients was determined by statistical analysis. The results indicated that GRO-β was highly expressed in CRC tissues, and that high GRO-β cytoplasmic expression was associated with the tumor location, extent of the primary tumor, and lymph node metastasis. Kaplan-Meier survival and Cox regression analysis revealed that high GRO-β expression was an independent indicator of poor prognosis for CRC patients. The results indicate that high GRO-β expression in CRC may correlate with an unfavorable outcome and facilitate cancer cell invasion and metastasis.

The high expressed serum soluble neural cell adhesion molecule, a high risk factor indicating hepatic encephalopathy in hepatocelular carcinoma patients

OBJECTIVE

To investigate whether the expression of serum soluble neural cell adhesion molecule (sNCAM) is associated with hepatic encephalopathy (HE) in hepatocelular carcinoma (HCC) patients.

MATERIALS AND METHODS

The Oncomine Cancer Microarray database was used to determine the clinical relevance of NCAM expression in different kinds of human cancers. Sera from 75 HCC cases enrolled in this study were assessed for expression of sNCAM by enzyme linked immunosorbent assay (ELISA).

RESULTS

Dependent on the Oncomine Cancer Microarray database analysis, NCAM was down regulated in 10 different kinds of cancer, like bladder cancer, brain and central nervous system cancer, while up-regulated in lung cancer, uterine corpus leiomyoma and sarcoma, compared to normal groups. Puzzlingly, NCAM expression demonstrated no significant difference between normal and HCC groups. However, we found by quantitative ELISA that the level of sNCAM in sera from HCC patients with HE (347.4±151.9 ng/ml) was significantly more up-regulated than that in HCC patients without HE (260.3±104.2 ng/ml), the p-value being 0.008. sNCAM may be an important risk factor of HE in HCC patients, the correlation coefficients was 0.278 (P< 0.05) on rank correlation analysis.

CONCLUSIONS

This study highlights that up-regulated level of serum sNCAM is associated with HE in HCC patients and suggests that the high expression can be used as an indicator.

The MAP3K ZAK, a novel modulator of ERK-dependent migration, is upregulated in colorectal cancer

Often described as a mediator of cell cycle arrest or as a pro-apoptotic factor in stressful conditions, the MAP3K ZAK (Sterile alpha motif and leucine zipper-containing kinase) has also been proven to positively regulate epidermal growth factor receptor (EGFR) and WNT signaling pathways, cancer cell proliferation and cellular neoplastic transformation. Here, we show that both isoforms of ZAK, ZAK-α and ZAK-β are key factors in cancer cell migration. While ZAK depletion reduced cell motility of HeLa and HCT116 cells, its overexpression triggered the activation of all three mitogen-activated protein kinases (MAPKs), extracellular signal-regulated kinase (ERK), c-JUN N-terminal kinase (JNK) and p38, as well as an increase in cell motion. On the contrary, the kinase-dead mutants, ZAK-α K45M and ZAK-β K45M, were not able to provoke such events, and instead exerted a dominant-negative effect on MAPK activation and cell migration. Pharmacological inhibition of ZAK by nilotinib, preventing ZAK-autophosphorylation and thereby auto-activation, led to the same results. Activated by epidermal growth factor (EGF), we further showed that ZAK constitutes an essential element of the EGF/ERK-dependent cell migration pathway. Using public transcriptomic databases and tissue microarrays, we finally established that, as strong factors of the EGFR signaling pathway, ZAK-α and/or ZAK-β transcripts and protein(s) are frequently upregulated in colorectal adenoma and carcinoma patients. Notably, gene set enrichment analysis disclosed a significant correlation between ZAK+ colorectal premalignant lesions and gene sets belonging to the MAPK/ERK and motility-related signaling pathways of the reactome database, strongly suggesting that ZAK induces such pro-tumoral reaction cascades in human cancers.

Comparative study on gene set and pathway topology-based enrichment methods

Background

Enrichment analysis is a popular approach to identify pathways or sets of genes which are significantly enriched in the context of differentially expressed genes. The traditional gene set enrichment approach considers a pathway as a simple gene list disregarding any knowledge of gene or protein interactions. In contrast, the new group of so called pathway topology-based methods integrates the topological structure of a pathway into the analysis.

Methods

We comparatively investigated gene set and pathway topology-based enrichment approaches, considering three gene set and four topological methods. These methods were compared in two extensive simulation studies and on a benchmark of 36 real datasets, providing the same pathway input data for all methods.

Results

In the benchmark data analysis both types of methods showed a comparable ability to detect enriched pathways. The first simulation study was conducted with KEGG pathways, which showed considerable gene overlaps between each other. In this study with original KEGG pathways, none of the topology-based methods outperformed the gene set approach. Therefore, a second simulation study was performed on non-overlapping pathways created by unique gene IDs. Here, methods accounting for pathway topology reached higher accuracy than the gene set methods, however their sensitivity was lower.

Conclusions

We conducted one of the first comprehensive comparative works on evaluating gene set against pathway topology-based enrichment methods. The topological methods showed better performance in the simulation scenarios with non-overlapping pathways, however, they were not conclusively better in the other scenarios. This suggests that simple gene set approach might be sufficient to detect an enriched pathway under realistic circumstances. Nevertheless, more extensive studies and further benchmark data are needed to systematically evaluate these methods and to assess what gain and cost pathway topology information introduces into enrichment analysis. Both types of methods for enrichment analysis require further improvements in order to deal with the problem of pathway overlaps.

Electronic supplementary material

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

Apc Restoration Promotes Cellular Differentiation and Reestablishes Crypt Homeostasis in Colorectal Cancer

The adenomatous polyposis coli (APC) tumor suppressor is mutated in the vast majority of human colorectal cancers (CRC) and leads to deregulated Wnt signaling. To determine whether Apc disruption is required for tumor maintenance, we developed a mouse model of CRC whereby Apc can be conditionally suppressed using a doxycycline-regulated shRNA. Apc suppression produces adenomas in both the small intestine and colon that, in the presence of Kras and p53 mutations, can progress to invasive carcinoma. In established tumors, Apc restoration drives rapid and widespread tumor-cell differentiation and sustained regression without relapse. Tumor regression is accompanied by the re-establishment of normal crypt-villus homeostasis, such that once aberrantly proliferating cells reacquire self-renewal and multi-lineage differentiation capability. Our study reveals that CRC cells can revert to functioning normal cells given appropriate signals and provide compelling in vivo validation of the Wnt pathway as a therapeutic target for treatment of CRC.

Diacylglycerol kinase-ζ regulates mTORC1 and lipogenic metabolism in cancer cells through SREBP-1

Diacylglycerol kinases (DGKs) transform diacylglycerol (DAG) into phosphatidic acid (PA), balancing the levels of these key metabolic and signaling lipids. We previously showed that PA derived from the DGKζ isoform promotes mammalian target of rapamycin complex 1 (mTORC1) activation. This function might be crucial for the growth and survival of cancer cells, especially for those resistant to the allosteric mTOR inhibitor rapamycin. How this positive function of DGKζ coordinates with DAG metabolism and signaling is unknown. In this study, we used a rapamycin-resistant colon cancer cell line as a model to address the role of DGKζ in tumor cells. We found that DGKζ predominated over other PA sources such as DGKα or phospholipase D to activate mTORC1, and that its activity was a component of the rapamycin-induced feedback loops. We show that the DGKζ DAG-consuming function is central to cell homeostasis, as DAG negatively regulates levels of the lipogenic transcription factor SREBP-1. Our findings suggest a model in which simultaneous regulation of DAG and PA levels by DGKζ is integrated with mTOR function to maintain tumor cell homeostasis; we provide new evidence of the crosstalk between mTOR and lipid metabolism that will be advantageous in the design of drug therapies.

Long Non-Coding RNAs: The Key Players in Glioma Pathogenesis

Long non-coding RNAs (LncRNAs) represent a novel class of RNAs with no functional protein-coding ability, yet it has become increasingly clear that interactions between lncRNAs with other molecules are responsible for important gene regulatory functions in various contexts. Given their relatively high expressions in the brain, lncRNAs are now thought to play important roles in normal brain development as well as diverse disease processes including gliomagenesis. Intriguingly, certain lncRNAs are closely associated with the initiation, differentiation, progression, recurrence and stem-like characteristics in glioma, and may therefore be exploited for the purposes of sub-classification, diagnosis and prognosis. LncRNAs may also serve as potential therapeutic targets as well as a novel biomarkers in the treatment of glioma. In this article, the functional aspects of lncRNAs, particularly within the central nervous system (CNS), will be briefly discussed, followed by highlights of the important roles of lncRNAs in mediating critical steps during glioma development. In addition, the key lncRNA players and their possible mechanistic pathways associated with gliomagenesis will be addressed.

Krüppel-like factor 9 (KLF9) prevents colorectal cancer through inhibition of interferon-related signaling

Expression of the transcription factor Krüppel-like factor 9 (KLF9) is frequently reduced in colorectal cancers, although a tumor suppressive role has not been established. To determine if KLF9 suppresses intestinal adenoma formation, we generated mice of distinct Klf9 genotypes in the background of the Apc (Min/+) mouse and compared their adenoma burdens at 16 weeks of age. While small intestine adenoma burden remained unchanged among Klf9 genotypes, male and female Apc(Min/+)/Klf9(-/-) and Apc(Min/+)/Klf9(+/-) mice exhibited significantly more colon adenomas than their Apc(Min/+)/Klf9(+/+) counterparts. Microarray analysis showed significant increases in the expression of interferon-induced genes in the colon mucosa of female Apc (Min/+)/Klf9(+/-) and Apc(Min/+)/Klf9(-/-) compared to Apc(Min/+)/Klf9(+/+) mice, prior to overt adenoma occurrence. Gene upregulation was confirmed by qPCR of colon mucosa and by siRNA knockdown of KLF9 in human HT29 colorectal cancer cells. Increases in expression of these genes were further augmented by supplementation with Interferon β1. Circulating levels of the cytokine, interferon-stimulated gene 15 (ISG15) were increased in Apc(Min/+)/Klf9(+/-) and Apc(Min/+)/Klf9(-/-) mice relative to Apc(Min/+)/Klf9(+/+). Additionally, colon mucosal levels of ISG15 were increased in Apc(Min/+)/Klf9(+/-) mice. Chromatin immunoprecipitation demonstrated KLF9 recruitment to the ISG15 promoter. Lastly, treatment with ISG15 suppressed apoptosis in HT29 cells, in the presence and absence of 5-fluorouracil (5FU). Results show KLF9 to be a haploinsufficient suppressor of colon tumorigenesis in Apc(Min/+) mice in part, by repression of ISG15 and the latter's antiapoptotic function.

SUMMARY

Krüppel-like factor 9 (KLF9) is a haploinsufficient tumor suppressor in the ApcMin/+ mouse colon by suppressing expression of ISG15, an apoptosis-inhibiting cytokine.

Endoplasmic reticulum ribosome-binding protein 1, RRBP1, promotes progression of colorectal cancer and predicts an unfavourable prognosis

BACKGROUND

Ribosome-binding protein 1 (RRBP1) has been implicated in the regulation of unfolded protein response, which is involved in almost every aspect of cancer development. We aimed to explore the significance of RRBP1 in the progression and prognosis of colorectal cancer (CRC).

METHODS

The study population consisted of 856 patients with stage I-III CRC from two hospitals. RRBP1 expression was examined by immunohistochemistry (IHC) in colorectal tissues. The correlation of RRBP1 expression and CRC occurrence was assessed in paired cancer-adjacent tissues. Factors contributing to prognosis were evaluated in a training-validation design with univariate and multivariate Cox analysis. Colorectal cancer aggressiveness caused by RRBP1 knockdown or overexpression was evaluated in CRC cells.

RESULTS

RRBP1 was aberrantly overexpressed in CRC. Compared with low-RRBP1 patients, high-RRBP1 patients had shorter disease-specific survival in the training (hazard ratio (HR), 2.423; 95% confidence interval (CI), 1.531-3.835) and validation cohorts (HR, 3.749; 95% CI, 2.166-6.448) in multivariate Cox analysis. High-RRBP1 independently predicted a shorter disease-free survival (HR, 4.821; 95% CI, 3.220-7.218) in the validation cohort. RRBP1 knockdown reduced the aggressiveness of CRC cells in vitro and inhibited the growth of CRC xenografts in vivo.

CONCLUSIONS

High RRBP1 expression facilitates CRC progression and predicts an unfavourable post-operative prognosis.

Meta-Analysis of Public Microarray Datasets Reveals Voltage-Gated Calcium Gene Signatures in Clinical Cancer Patients

Voltage-gated calcium channels (VGCCs) are well documented to play roles in cell proliferation, migration, and apoptosis; however, whether VGCCs regulate the onset and progression of cancer is still under investigation. The VGCC family consists of five members, which are L-type, N-type, T-type, R-type and P/Q type. To date, no holistic approach has been used to screen VGCC family genes in different types of cancer. We analyzed the transcript expression of VGCCs in clinical cancer tissue samples by accessing ONCOMINE (www.oncomine.org), a web-based microarray database, to perform a systematic analysis. Every member of the VGCCs was examined across 21 different types of cancer by comparing mRNA expression in cancer to that in normal tissue. A previous study showed that altered expression of mRNA in cancer tissue may play an oncogenic role and promote tumor development; therefore, in the present findings, we focus only on the overexpression of VGCCs in different types of cancer. This bioinformatics analysis revealed that different subtypes of VGCCs (CACNA1C, CACNA1D, CACNA1B, CACNA1G, and CACNA1I) are implicated in the development and progression of diverse types of cancer and show dramatic up-regulation in breast cancer. CACNA1F only showed high expression in testis cancer, whereas CACNA1A, CACNA1C, and CACNA1D were highly expressed in most types of cancer. The current analysis revealed that specific VGCCs likely play essential roles in specific types of cancer. Collectively, we identified several VGCC targets and classified them according to different cancer subtypes for prospective studies on the underlying carcinogenic mechanisms. The present findings suggest that VGCCs are possible targets for prospective investigation in cancer treatment.

Mobilization of epithelial mesenchymal transition genes distinguishes active from inactive lesional tissue in patients with ulcerative colitis

Ulcerative colitis (UC) is a chronic, relapsing and debilitating idiopathic inflammation, with variable and complex pathophysiologies. Our objective was to elucidate patterns of gene expression underlying the progression of UC disease. Single endoscopic pinch FFPE biopsies (n = 41) were sampled at both active and inactive stages at the same site in individual UC patients and compared with each other and with non-inflammatory bowel disease healthy controls. Gene expression results were validated by quantitative reverse transcriptase-PCR (QRT-PCR), and results at the protein level were validated by immunohistochemistry and western blot. Analysis of microarray results demonstrated that UC patients in remission display an intermediate gene expression phenotype between active UC patients and controls. It is clear that UC active site recovery does not revert fully back to a healthy control phenotype. Both UC active and inactive tissue displayed evidence, at both the gene expression and protein level, of a positive precancerous state as indicated by increases in the expression of Chitinase 3-Like-1, and the colorectal cancer metastasis marker MMP1. A key distinguishing feature between active and inactive UC, however, was the mobilization of marker genes and proteins for the Epithelial Mesenchymal Transition (EMT) pathway only in active UC. Analysis of the gene expression signatures associated with UC remission identified multiple pathways which appear to be permanently dysregulated in UC patients at formerly active sites in spite of clear histological recovery. Among these pathways, the EMT pathway was specifically up-regulated only in active UC emphasizing the potential for cancer progression in these patients.

KIAA1199 as a potential diagnostic biomarker of rheumatoid arthritis related to angiogenesis

Introduction

Our previous proteomic study on fibroblast-like synoviocytes (FLSs) derived from the synovial tissues found that the expression of KIAA1199 was higher in rheumatoid arthritis (RA) patients than in healthy controls. The aim of this study was to examine the biological function of KIAA1199 and evaluate its clinical diagnosis value in RA.

Methods

The over-expression of KIAA1199 was verified by quantitative real-time polymerase chain reaction (qPCR), Immunohistochemistry, Immunofluorescence and enzyme linked immunosorbent assay (ELISA) in inactive and active RA patients and healthy controls. The effect of KIAA1199 expression on FLSs proliferation, angiogenesis and related pathway were analyzed by MTT, cell migration, tube formation, chorioallantoic membrane (CAM) assay, qPCR and western-blotting after KIAA1199 knockdown and over-expression.

Results

The verification results show the up-regulation of KIAA1199 in RA patients at mRNA and protein level as compared to that in healthy controls. ELISA and receiver operator characteristic (ROC) analysis shows that KIAA1199 concentration in serum, synovial fluid and synovial tissues could be used as dependable biomarkers for the diagnosis of active RA, provided an area under roc curve (AUC) of 0.83, 0.92 and 0.92. Sensitivity and specificity, which were determined by cut-off points, reached 72% 84% and 80% in sensitivity and 80%, 93.3%, 93.3% in specificity, respectively. Moreover, KIAA1199 also enhance the proliferation and angiogenesis of synovial membrane, and KIAA1199/ PLXNB3/ SEMA5A/CTGF axis may be a newly found pathway enhancing cell proliferation and angiogenesis.

Conclusion

KIAA1199 may be a potential diagnostic biomarker of RA related to angiogenesis.

Electronic supplementary material

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

Regenerating gene 1B silencing inhibits colon cancer cell HCT116 proliferation and invasion

BACKGROUND

The human regenerating gene 1B (REG1B) is found to be frequently up-regulated in many types of human tumors. It is unclear whether REG1B expression may have therapeutic value in colorectal carcinoma. Additionally, how REG1B is associated with the clinical features of colorectal carcinoma is not known. To investigate the relationship between REG1B and colorectal cancer, we analyzed REG1B expression in clinical specimens and cell lines and the effect of down-regulation of REG1B by short hairpin RNA (shRNA) in HCT116 cells.

METHODS

Paraffin-embedded specimens from 30 pairs of colorectal cancer tissues and adjacent colon tissues were used to investigate the expression of REG1B by immunohistochemistry. We also examined whether REG1B itself may be related to cell proliferation, cell cycle arrest, apoptosis, migration and invasion in colon cancer HCT116 cells.

RESULTS

Our results showed that REG1B was highly expressed in colorectal carcinoma and was significantly associated with cell differentiation status. The results also illustrated that REG1B silencing with shRNA inhibited cell proliferation, migration and invasion but did not induce apoptosis. Furthermore, down-regulation of REG1B induces G1-phase cell cycle arrest in colon cancer cells.

CONCLUSIONS

Knockdown of REG1B can inhibit cell proliferation, migration and invasion. It may act by a mechanism regulating cell cycle progression. Thus, REG1B may be a novel candidate therapeutic target for colorectal cancer.

The Novel Gene CRNDE Encodes a Nuclear Peptide (CRNDEP) Which Is Overexpressed in Highly Proliferating Tissues

CRNDE, recently described as the lncRNA-coding gene, is overexpressed at RNA level in human malignancies. Its role in gametogenesis, cellular differentiation and pluripotency has been suggested as well. Herein, we aimed to verify our hypothesis that the CRNDE gene may encode a protein product, CRNDEP. By using bioinformatics methods, we identified the 84-amino acid ORF encoded by one of two CRNDE transcripts, previously described by our research team. This ORF was cloned into two expression vectors, subsequently utilized in localization studies in HeLa cells. We also developed a polyclonal antibody against CRNDEP. Its specificity was confirmed in immunohistochemical, cellular localization, Western blot and immunoprecipitation experiments, as well as by showing a statistically significant decrease of endogenous CRNDEP expression in the cells with transient shRNA-mediated knockdown of CRNDE. Endogenous CRNDEP localizes predominantly to the nucleus and its expression seems to be elevated in highly proliferating tissues, like the parabasal layer of the squamous epithelium, intestinal crypts or spermatocytes. After its artificial overexpression in HeLa cells, in a fusion with either the EGFP or DsRed Monomer fluorescent tag, CRNDEP seems to stimulate the formation of stress granules and localize to them. Although the exact role of CRNDEP is unknown, our preliminary results suggest that it may be involved in the regulation of the cell proliferation. Possibly, CRNDEP also participates in oxygen metabolism, considering our in silico results, and the correlation between its enforced overexpression and the formation of stress granules. This is the first report showing the existence of a peptide encoded by the CRNDE gene.

IL-33 activates tumor stroma to promote intestinal polyposis

Tumor epithelial cells develop within a microenvironment consisting of extracellular matrix, growth factors, and cytokines produced by nonepithelial stromal cells. In response to paracrine signals from tumor epithelia, stromal cells modify the microenvironment to promote tumor growth and metastasis. Here, we identify interleukin 33 (IL-33) as a regulator of tumor stromal cell activation and mediator of intestinal polyposis. In human colorectal cancer, IL-33 expression was induced in the tumor epithelium of adenomas and carcinomas, and expression of the IL-33 receptor, IL1RL1 (also referred to as IL1-R4 or ST2), localized predominantly to the stroma of adenoma and both the stroma and epithelium of carcinoma. Genetic and antibody abrogation of responsiveness to IL-33 in the Apc(Min/+) mouse model of intestinal tumorigenesis inhibited proliferation, induced apoptosis, and suppressed angiogenesis in adenomatous polyps, which reduced both tumor number and size. Similar to human adenomas, IL-33 expression localized to tumor epithelial cells and expression of IL1RL1 associated with two stromal cell types, subepithelial myofibroblasts and mast cells, in Apc(Min/+) polyps. In vitro, IL-33 stimulation of human subepithelial myofibroblasts induced the expression of extracellular matrix components and growth factors associated with intestinal tumor progression. IL-33 deficiency reduced mast cell accumulation in Apc(Min/+) polyps and suppressed the expression of mast cell-derived proteases and cytokines known to promote polyposis. Based on these findings, we propose that IL-33 derived from the tumor epithelium promotes polyposis through the coordinated activation of stromal cells and the formation of a protumorigenic microenvironment.

Can "normal" protein expression ranges be estimated with high-throughput proteomics?

Although biological science discovery often involves comparing conditions to a normal state, in proteomics little is actually known about normal. Two Human Proteome studies featured in Nature offer new insights into protein expression and an opportunity to assess how high-throughput proteomics measures normal protein ranges. We use data from these studies to estimate technical and biological variability in protein expression and compare them to other expression data sets from normal tissue. Results show that measured protein expression across same-tissue replicates vary by ±4- to 10-fold for most proteins. Coefficients of variation (CV) for protein expression measurements range from 62% to 117% across different tissue experiments; however, adjusting for technical variation reduced this variability by as much as 50%. In addition, the CV could also be reduced by limiting comparisons to proteins with at least 3 or more unique peptide identifications as the CV was on average 33% lower than for proteins with 2 or fewer peptide identifications. We also selected 13 housekeeping proteins and genes that were expressed across all tissues with low variability to determine their utility as a reference set for normalization and comparative purposes. These results present the first step toward estimating normal protein ranges by determining the variability in expression measurements through combining publicly available data. They support an approach that combines standard protocols with replicates of normal tissues to estimate normal protein ranges for large numbers of proteins and tissues. This would be a tremendous resource for normal cellular physiology and comparisons of proteomics studies.

BioMiner: Paving the Way for Personalized Medicine

Personalized medicine is promising a revolution for medicine and human biology in the 21st century. The scientific foundation for this revolution is accomplished by analyzing biological high-throughput data sets from genomics, transcriptomics, proteomics, and metabolomics. Currently, access to these data has been limited to either rather simple Web-based tools, which do not grant much insight or analysis by trained specialists, without firsthand involvement of the physician. Here, we present the novel Web-based tool “BioMiner,” which was developed within the scope of an international and interdisciplinary project (SYSTHER†) and gives access to a variety of high-throughput data sets. It provides the user with convenient tools to analyze complex cross-omics data sets and grants enhanced visualization abilities. BioMiner incorporates transcriptomic and cross-omics high-throughput data sets, with a focus on cancer. A public instance of BioMiner along with the database is available at http://systherDB.microdiscovery.de/, login and password: “systher”; a tutorial detailing the usage of BioMiner can be found in the Supplementary File.

Regulation of VDR Expression in Apc-Mutant Mice, Human Colon Cancers and Adenomas

One variable that may affect the ability of vitamin D to reduce colon cancer risk is the expression of its high-affinity receptor, VDR. Here, we show that vitamin D does not reduce tumor formation in Apc(Δ14/+) mice and that VDR expression is lost in the majority of the colon tumor cells. The extent of VDR loss corresponded inversely to the level of β-catenin nuclear localization and could be observed in early lesions composed of just a few crypts. Analysis of reported VDR regulators showed that the repressing class I histone deacetylases (HDAC) were significantly elevated in the tumors (up to 4-fold), whereas the VDR-activating retinoid X receptors (RXR) were downregulated (∼50%). Expression of the Slug repressor was also increased, but was found primarily in stromal cells. Analysis of epigenetically active compounds on colon cell lines and intestinal organoids showed that HDAC inhibitors were particularly adept at stimulating VDR expression. Treatment of tumor-bearing Apc(Δ14/+) mice with the HDAC inhibitor panobinostat increased VDR expression in the tumors and normal mucosa. The RXR agonist bexarotene failed to activate VDR expression, indicating that RXR ligands were not limiting. Analysis of human microarray data indicated that VDR mRNA is frequently downregulated in colon adenomas, which correlated positively with RXRA expression and inversely with HDAC 2 and 8 expression. Human adenomas showed variable VDR protein expression levels, both between and within individual lesions. Determining the mechanisms of VDR regulation in colon neoplasms may significantly enhance our ability to use vitamin D as a cancer prevention agent.

Tissue-specific patterns of gene expression in the epithelium and stroma of normal colon in healthy individuals in an aspirin intervention trial

Background

Regular aspirin use reduces colon adenoma and carcinoma incidence. UDP-glucuronosyltransferases (UGT) are involved in aspirin metabolism and clearance, and variant alleles in UGT1A6 have been shown to alter salicylic acid metabolism and risk of colon neoplasia.

Methods

In a randomized, cross-over, placebo-controlled trial of 44 healthy men and women, homozygous for UGT1A6*1 or UGT1A6*2, we explored differences between global epithelial and stromal expression, using Affymetrix U133 + 2.0 microarrays and tested effects of 60-day aspirin supplementation (325 mg/d) on epithelial and stromal gene expression and colon prostaglandin E2 (PGE2) levels.

Results

No statistically significant differences in gene expression were observed in response to aspirin or UGT1A6 genotype, but tissue PGE2 levels were lower with aspirin compared to placebo (p <0.001). Transcripts differentially expressed between epithelium and stroma (N = 4916, P <0.01, false discovery rate <0.001), included a high proportion of genes involved in cell signaling, cellular movement, and cancer. Genes preferentially expressed in epithelium were involved in drug and xenobiotic metabolism, fatty acid and lipid metabolism, apoptosis signaling, and ion transport. Genes preferentially expressed in stroma included those involved in inflammation, cellular adhesion, and extracellular matrix production. Wnt-Tcf4 pathway genes were expressed in both epithelium and stroma but differed by subcellular location.

Conclusions

These results suggest that, in healthy individuals, subtle effects of aspirin on gene expression in normal colon tissue are likely overwhelmed by inter-individual variability in microarray analyses. Differential expression of critical genes between colonic epithelium and stroma suggest that these tissue types need to be considered separately.

Electronic supplementary material

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

miR-182 promotes cell growth and invasion by targeting forkhead box F2 transcription factor in colorectal cancer

Forkhead box F2 transcription factor (FoxF2) has been described to promote organ development, extracellular matrix (ECM) synthesis and epithelial-mesenchymal interaction. Although recent studies reported decreased FoxF2 expression in several types of cancers, indicating its potential role in carcinogenesis, the mechanistic role of FoxF2 is yet to be explored. MicroRNAs (miRNAs) are strongly implicated in carcinogenesis. The oncogenetic properties of miR-182 have been described in multiple cancers. In the present study, we aimed to investigate the role of miR-182 in colorectal cancer (CRC) and identify the regulation of FoxF2 by miR-182. Bioinformatic analyses on gene expression profiling datasets showed decreased FoxF2 expression in colorectal adenomas, primary tumors compared to normal colon epithelial and a negative association between FoxF2 and β-catenin expression. Restoration of FoxF2 in CRC cells suppressed β-catenin expression and simultaneously inhibited cell growth and invasion. Furthermore, we observed that miR-182 was aberrantly upregulated in CRC. Knockdown of miR-182 in CRC cells impeded cell growth and invasion. The direct binding of miR-182 to the 3' untranslated region (3'UTR) of FoxF2 mRNA was confirmed using a luciferase reporter gene assay. Importantly, elevated FoxF2 expression was observed in miR-182-knockdown cells with a simultaneous reduction in β-catenin. In conclusion, the present study describes a potential mechanism underlying an miR-182/FoxF2 link contributing to CRC development. miR-182-induced downregulation of FoxF2 partly accounts for increased activity of β-catenin signaling. Inhibition of miR-182 represents a potential strategy against CRC.

NRMT2 is an N-terminal monomethylase that primes for its homologue NRMT1

NRMT (N-terminal regulator of chromatin condensation 1 methyltransferase) was the first eukaryotic methyltransferase identified to specifically methylate the free α-amino group of proteins. Since the discovery of this N-terminal methyltransferase, many new substrates have been identified and the modification itself has been shown to regulate DNA-protein interactions. Sequence analysis predicts one close human homologue of NRMT, METTL11B (methyltransferase-like protein 11B, now renamed NRMT2). We show in the present paper for the first time that NRMT2 also has N-terminal methylation activity and recognizes the same N-terminal consensus sequences as NRMT (now NRMT1). Both enzymes have similar tissue expression and cellular localization patterns. However, enzyme assays and MS experiments indicate that they differ in their specific catalytic functions. Although NRMT1 is a distributive methyltransferase that can mono-, di- and tri-methylate its substrates, NRMT2 is primarily a monomethylase. Concurrent expression of NRMT1 and NRMT2 accelerates the production of trimethylation, and we propose that NRMT2 activates NRMT1 by priming its substrates for trimethylation.

Case-control study of the PERIOD3 clock gene length polymorphism and colorectal adenoma formation

Clock genes are expressed in a self-perpetuating, circadian pattern in virtually every tissue including the human gastrointestinal tract. They coordinate cellular processes critical for tumor development, including cell proliferation, DNA damage response and apoptosis. Circadian rhythm disturbances have been associated with an increased risk for colon cancer and other cancers. This mechanism has not been elucidated, yet may involve dysregulation of the 'period' (PER) clock genes, which have tumor suppressor properties. A variable number tandem repeat (VNTR) in the PERIOD3 (PER3) gene has been associated with sleep disorders, differences in diurnal hormone secretion, and increased premenopausal breast cancer risk. Susceptibility related to PER3 has not been examined in conjunction with adenomatous polyps. This exploratory case-control study was the first to test the hypothesis that the 5-repeat PER3 VNTR sequence is associated with increased odds of adenoma formation. Information on demographics, medical history, occupation and lifestyle was collected prior to colonoscopy. Cases (n=49) were individuals with at least one histopathologically confirmed adenoma. Controls (n=97) included patients with normal findings or hyperplastic polyps not requiring enhanced surveillance. Unconditional multiple logistic regression was used to calculate odds ratios (ORs) with 95% confidence intervals (CIs), after adjusting for potential confounding. Adenomas were detected in 34% of participants. Cases were more likely to possess the 5-repeat PER3 genotype relative to controls (4/5 OR, 2.1; 95% CI, 0.9-4.8; 5/5 OR, 5.1; 95% CI, 1.4-18.1; 4/5+5/5 OR, 2.5; 95% CI, 1.7-5.4). Examination of the Oncomine microarray database indicated lower PERIOD gene expression in adenomas relative to adjacent normal tissue. Results suggest a need for follow-up in a larger sample.

Personalized Disease Phenotypes from Massive OMICs Data

Application of new high-throughput technologies in molecular medicine collects massive data for hundreds to thousands of persons in large cohort studies by characterizing the phenotype of each individual on a personalized basis. The chapter aims at increasing our understanding of disease genesis and progression and to improve diagnosis and treatment. New methods are needed to handle such “big data.” Machine learning enables one to recognize and to visualize complex data patterns and to make decisions potentially relevant for diagnosis and treatment. The authors address these tasks by applying the method of self-organizing maps and present worked examples from different disease entities of the colon ranging from inflammation to cancer.

HOXD-AS1 is a novel lncRNA encoded in HOXD cluster and a marker of neuroblastoma progression revealed via integrative analysis of noncoding transcriptome

Background

Long noncoding RNAs (lncRNAs) constitute a major, but poorly characterized part of human transcriptome. Recent evidence indicates that many lncRNAs are involved in cancer and can be used as predictive and prognostic biomarkers. Significant fraction of lncRNAs is represented on widely used microarray platforms, however they have usually been ignored in cancer studies.

Results

We developed a computational pipeline to annotate lncRNAs on popular Affymetrix U133 microarrays, creating a resource allowing measurement of expression of 1581 lncRNAs. This resource can be utilized to interrogate existing microarray datasets for various lncRNA studies. We found that these lncRNAs fall into three distinct classes according to their statistical distribution by length. Remarkably, these three classes of lncRNAs were co-localized with protein coding genes exhibiting distinct gene ontology groups. This annotation was applied to microarray analysis which identified a 159 lncRNA signature that discriminates between localized and metastatic stages of neuroblastoma. Analysis of an independent patient cohort revealed that this signature differentiates also relapsing from non-relapsing primary tumors. This is the first example of the signature developed via the analysis of expression of lncRNAs solely. One of these lncRNAs, termed HOXD-AS1, is encoded in HOXD cluster. HOXD-AS1 is evolutionary conserved among hominids and has all bona fide features of a gene. Studying retinoid acid (RA) response of SH-SY5Y cell line, a model of human metastatic neuroblastoma, we found that HOXD-AS1 is a subject to morphogenic regulation, is activated by PI3K/Akt pathway and itself is involved in control of RA-induced cell differentiation. Knock-down experiments revealed that HOXD-AS1 controls expression levels of clinically significant protein-coding genes involved in angiogenesis and inflammation, the hallmarks of metastatic cancer.

Conclusions

Our findings greatly extend the number of noncoding RNAs functionally implicated in tumor development and patient treatment and highlight their role as potential prognostic biomarkers of neuroblastomas.

The molecular signature of the stroma response in prostate cancer-induced osteoblastic bone metastasis highlights expansion of hematopoietic and prostate epithelial stem cell niches

The reciprocal interaction between cancer cells and the tissue-specific stroma is critical for primary and metastatic tumor growth progression. Prostate cancer cells colonize preferentially bone (osteotropism), where they alter the physiological balance between osteoblast-mediated bone formation and osteoclast-mediated bone resorption, and elicit prevalently an osteoblastic response (osteoinduction). The molecular cues provided by osteoblasts for the survival and growth of bone metastatic prostate cancer cells are largely unknown. We exploited the sufficient divergence between human and mouse RNA sequences together with redefinition of highly species-specific gene arrays by computer-aided and experimental exclusion of cross-hybridizing oligonucleotide probes. This strategy allowed the dissection of the stroma (mouse) from the cancer cell (human) transcriptome in bone metastasis xenograft models of human osteoinductive prostate cancer cells (VCaP and C4-2B). As a result, we generated the osteoblastic bone metastasis-associated stroma transcriptome (OB-BMST). Subtraction of genes shared by inflammation, wound healing and desmoplastic responses, and by the tissue type-independent stroma responses to a variety of non-osteotropic and osteotropic primary cancers generated a curated gene signature ("Core" OB-BMST) putatively representing the bone marrow/bone-specific stroma response to prostate cancer-induced, osteoblastic bone metastasis. The expression pattern of three representative Core OB-BMST genes (PTN, EPHA3 and FSCN1) seems to confirm the bone specificity of this response. A robust induction of genes involved in osteogenesis and angiogenesis dominates both the OB-BMST and Core OB-BMST. This translates in an amplification of hematopoietic and, remarkably, prostate epithelial stem cell niche components that may function as a self-reinforcing bone metastatic niche providing a growth support specific for osteoinductive prostate cancer cells. The induction of this combinatorial stem cell niche is a novel mechanism that may also explain cancer cell osteotropism and local interference with hematopoiesis (myelophthisis). Accordingly, these stem cell niche components may represent innovative therapeutic targets and/or serum biomarkers in osteoblastic bone metastasis.

B-Raf inhibitors induce epithelial differentiation in BRAF-mutant colorectal cancer cells

BRAF mutations are associated with aggressive, less-differentiated and therapy-resistant colorectal carcinoma. However, the underlying mechanisms for these correlations remain unknown. To understand how oncogenic B-Raf contributes to carcinogenesis, in particular to aspects other than cellular proliferation and survival, we generated three isogenic human colorectal carcinoma cell line models in which we can dynamically modulate the expression of the B-Raf(V600E) oncoprotein. Doxycyclin-inducible knockdown of endogenous B-Raf(V600E) decreases cellular motility and invasion in conventional and three-dimensional (3D) culture, whereas it promotes cell-cell contacts and induces various hallmarks of differentiated epithelia. Importantly, all these effects are recapitulated by B-Raf (PLX4720, vemurafenib, and dabrafenib) or MEK inhibitors (trametinib). Surprisingly, loss of B-Raf(V600E) in HT29 xenografts does not only stall tumor growth, but also induces glandular structures with marked expression of CDX2, a tumor-suppressor and master transcription factor of intestinal differentiation. By performing the first transcriptome profiles of PLX4720-treated 3D cultures of HT29 and Colo-205 cells, we identify several upregulated genes linked to epithelial differentiation and effector functions, such as claudin-1, a Cdx-2 target gene encoding a critical tight junction component. Thereby, we provide a mechanism for the clinically observed correlation between mutant BRAF and the loss of Cdx-2 and claudin-1. PLX4720 also suppressed several metastasis-associated transcripts that have not been implicated as targets, effectors or potential biomarkers of oncogenic B-Raf signaling so far. Together, we identify a novel facet of clinically applied B-Raf or MEK inhibitors by showing that they promote cellular adhesion and differentiation of colorectal carcinoma cells.

NF-κB-induced KIAA1199 promotes survival through EGFR signalling

Constitutive activation of EGFR- and NF-κB-dependent pathways is a hallmark of cancer, yet signalling proteins that connect both oncogenic cascades are poorly characterized. Here we define KIAA1199 as a BCL-3- and p65-dependent gene in transformed keratinocytes. KIAA1199 expression is enhanced on human papillomavirus (HPV) infection and is aberrantly expressed in clinical cases of cervical (pre)neoplastic lesions. Mechanistically, KIAA1199 binds Plexin A2 and protects from Semaphorin 3A-mediated cell death by promoting EGFR stability and signalling. Moreover, KIAA1199 is an EGFR-binding protein and KIAA1199 deficiency impairs EGF-dependent Src, MEK1 and ERK1/2 phosphorylations. Therefore, EGFR stability and signalling to downstream kinases requires KIAA1199. As such, KIAA1199 promotes EGF-mediated epithelial-mesenchymal transition (EMT). Taken together, our data define KIAA1199 as an oncogenic protein induced by HPV infection and constitutive NF-κB activity that transmits pro-survival and invasive signals through EGFR signalling.

Iro/IRX transcription factors negatively regulate Dpp/TGF-β pathway activity during intestinal tumorigenesis

Activating mutations in Wnt and EGFR/Ras signaling pathways are common in colorectal cancer (CRC). Remarkably, clonal co-activation of these pathways in the adult Drosophila midgut induces "tumor-like" overgrowths. Here, we show that, in these clones and in CRC cell lines, Dpp/TGF-β acts as a tumor suppressor. Moreover, we discover that the Iroquois/IRX-family-protein Mirror downregulates the transcription of core components of the Dpp pathway, reducing its tumor suppressor activity. We also show that this genetic interaction is conserved in human CRC cells, where the Iro/IRX proteins IRX3 and IRX5 diminish the response to TGF-β. IRX3 and IRX5 are upregulated in human adenomas, and their levels correlate inversely with the gene expression signature of response to TGF-β. In addition, Irx5 expression confers a growth advantage in the presence of TGF-β, but is selected against in its absence. Together, our results identify a set of Iro/IRX proteins as conserved negative regulators of Dpp/TGF-β activity. We propose that during the characteristic adenoma-to-carcinoma transition of human CRC, the activity of IRX proteins could reduce the sensitivity to the cytostatic effect of TGF-β, conferring a growth advantage to tumor cells prior to the acquisition of mutations in TGF-β pathway components.

Extracellular superoxide dismutase induces mouse embryonic fibroblast proliferative burst, growth arrest, immortalization, and consequent in vivo tumorigenesis

AIMS

Rat sarcoma virus (RAS)-induced tumorigenesis has been suggested to follow a three-stage model consisting of an initial RAS activation, senescence induction, and evasion of p53-dependent senescence checkpoints. While reactive oxygen species act as second messengers in RAS-induced senescence, they are also involved in oncogenic transformation by inducing proliferation and promoting mutations. In the current work, we investigated the role of extracellular superoxide dismutase (SOD3) in RAS-induced senescence and immortalization in vitro and in vivo. We used a mouse embryonic fibroblast (MEF) primary cell model along with immortalized and transformed human cell lines derived from papillary and anaplastic thyroid cancer.

RESULTS

Based on our data, sod3 RNA interference in H-RasV12-transduced cells markedly inhibited cell growth, while sod3 over-expression in MEFs initially caused a proliferative burst followed by the activation of DNA damage checkpoints, induction of p53-p21 signal transduction, and senescence. Subsequently, sod3-transduced MEF cells developed co-operative p21-p16 down-regulation and acquired transformed cell characteristics such as increased telomerase activity, loss of contact inhibition, growth in low-nutrient conditions, and in vivo tumorigenesis. Interestingly, as previously reported with RAS, we showed a dose-dependent response to SOD3 in vitro and in vivo involving transcriptional and non-transcriptional regulatory mechanisms.

INNOVATION

SOD3 may mediate H-RasV12-induced initiation of primary cell immortalization.

CONCLUSIONS

Our results indicate that SOD3 influences growth signaling in primary and cancer cells downstream of the ras oncogene and could serve as a therapy target at an early tumorigenesis phase.