Bioinformatics Analysis of RNA-seq Data Reveals Genes Related to Cancer Stem Cells in Colorectal Cancerogenesis

Cancer stem cells (CSC) play one of the crucial roles in the pathogenesis of various cancers, including colorectal cancer (CRC). Although great efforts have been made regarding our understanding of the cancerogenesis of CRC, CSC involvement in CRC development is still poorly understood. Using bioinformatics and RNA-seq data of normal mucosa, colorectal adenoma, and carcinoma (n = 106) from GEO and TCGA, we identified candidate CSC genes and analyzed pathway enrichment analysis (PEI) and protein–protein interaction analysis (PPI). Identified CSC-related genes were validated using qPCR and tissue samples from 47 patients with adenoma, adenoma with early carcinoma, and carcinoma without and with lymph node metastasis and were compared to normal mucosa. Six CSC-related genes were identified: ANLN, CDK1, ECT2, PDGFD, TNC, and TNXB. ANLN, CDK1, ECT2, and TNC were differentially expressed between adenoma and adenoma with early carcinoma. TNC was differentially expressed in CRC without lymph node metastases whereas ANLN, CDK1, and PDGFD were differentially expressed in CRC with lymph node metastases compared to normal mucosa. ANLN and PDGFD were differentially expressed between carcinoma without and with lymph node metastasis. Our study identified and validated CSC-related genes that might be involved in early stages of CRC development (ANLN, CDK1, ECT2, TNC) and in development of metastasis (ANLN, PDGFD).     

RNA-Sequencing Based microRNA Expression Signature of Colorectal Cancer: The Impact of Oncogenic Targets Regulated by miR-490-3p

To elucidate novel aspects of the molecular pathogenesis of colorectal cancer (CRC), we have created a new microRNA (miRNA) expression signature based on RNA-sequencing. Analysis of the signature showed that 84 miRNAs were upregulated, and 70 were downregulated in CRC tissues. Interestingly, our signature indicated that both guide and passenger strands of some miRNAs were significantly dysregulated in CRC tissues. These findings support our earlier data demonstrating the involvement of miRNA passenger strands in cancer pathogenesis. Our study focused on downregulated miR-490-3p and investigated its tumor-suppressive function in CRC cells. We successfully identified a total of 38 putative oncogenic targets regulated by miR-490-3p in CRC cells. Among these targets, the expression of three genes (IRAK1: p = 0.0427, FUT1: p = 0.0468, and GPRIN2: p = 0.0080) significantly predicted 5-year overall survival of CRC patients. Moreover, we analyzed the direct regulation of IRAK1 by miR-490-3p, and its resultant oncogenic function in CRC cells. Thus, we have clarified a part of the molecular pathway of CRC based on the action of tumor-suppressive miR-490-3p. This new miRNA expression signature of CRC will be a useful tool for elucidating new molecular pathogenesis in this disease.     

Pancancer Analysis of the Oncogenic and Prognostic Role of NOL7: A Potential Target for Carcinogenesis and Survival

Despite growing evidence suggesting the critical function of NOL7 in cancer initiation and development, a systematic pancancer analysis of NOL7 is lacking. Herein, we present a comprehensive study of NOL7 which aimed to explore its potential role and detailed mechanisms across 33 human tumors based on The Cancer Genome Atlas (TCGA) and Clinical Proteomic Tumor Analysis Consortium (CATPAC) databases. As a result, both gene and protein levels of NOL7 were found to be increased in various tumor tissues, including breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), and head and neck squamous cell carcinoma (HNSC) as compared with corresponding normal tissues. Meanwhile, dysregulated NOL7 expression was found to be closely related to pathological stage and prognosis in several cancers, including LIHC, ovarian serous cystadenocarcinoma (OV), and bladder urothelial carcinoma (BLCA). The DNA methylation level of NOL7 was found to be decreased in most cancers and to be negatively associated with NOL7 expression. Furthermore, NOL7 expression was determined to be significantly associated with levels of infiltrating cells and immune checkpoint genes, including HMGB1. Analysis of NOL7-related genes revealed that RNA metabolism pathways, including “ribosome biogenesis”, “spliceosome”, and “RNA transport”, were mainly involved in the functional mechanism of NOL7 in human cancers. In summary, this pancancer study characterized the relationship between NOL7 expression and clinicopathologic features in multiple cancer types and further showed its potential regulatory network in human cancers. It represents a systemic analysis for further functional and therapeutic studies of NOL7 and highlights its predictive value with respect to the carcinogenesis and prognosis of various cancers, especially LIHC.     

Glycoproteomic and Phenotypic Elucidation of B4GALNT2 Expression Variants in the SID Histo-Blood Group System

The Sd^a histo-blood group antigen (GalNAcβ1-4(NeuAcα2-3)Galβ-R) is implicated in various infections and constitutes a potential biomarker for colon cancer. Sd(a−) individuals (2–4% of Europeans) may produce anti-Sd^a, which can lead to incompatible blood transfusions, especially if donors with the high-expressing Sd(a++)/Cad phenotype are involved. We previously reported the association of B4GALNT2 mutations with Sd(a−), which established the SID blood-group system. The present study provides causal proof underpinning this correlation. Sd(a−) HEK293 cells were transfected with different B4GALNT2 constructs and evaluated by immunostaining and glycoproteomics. The predominant SID^null candidate allele with rs7224888:T>C (p.Cys406Arg) abolished Sd^a synthesis, while this antigen was detectable as N- or O-glycans on glycoproteins following transfection of wildtype B4GALNT2. Surprisingly, two rare missense variants, rs148441237:A>G and rs61743617:C>T, found in a Sd(a−) compound heterozygote, gave results similar to wildtype. To elucidate on whether Sd(a++)/Cad also depends on B4GALNT2 alterations, this gene was sequenced in five individuals. No Cad-specific changes were identified, but a detailed erythroid Cad glycoprotein profile was obtained, especially for glycophorin-A (GLPA) O-glycosylation, equilibrative nucleoside transporter 1 (S29A1) O-glycosylation, and band 3 anion transport protein (B3AT) N-glycosylation. In conclusion, the p.Cys406Arg β4GalNAc-T2 variant causes Sd^a-deficiency in humans, while the enigmatic Cad phenotype remains unresolved, albeit further characterized.     

HOXB9 Overexpression Promotes Colorectal Cancer Progression and Is Associated with Worse Survival in Liver Resection Patients for Colorectal Liver Metastases

As is known, HOXB9 is an important factor affecting disease progression and overall survival (OS) in cancer. However, its role in colorectal cancer (CRC) remains unclear. We aimed to explore the role of HOXB9 in CRC progression and its association with OS in colorectal liver metastases (CRLM). We analysed differential HOXB9 expression in CRC using the Tissue Cancer Genome Atlas database (TCGA). We modulated HOXB9 expression in vitro to assess its impact on cell proliferation and epithelial-mesenchymal transition (EMT). Lastly, we explored the association of HOXB9 protein expression with OS, using an institutional patient cohort (n = 110) who underwent liver resection for CRLM. Furthermore, HOXB9 was upregulated in TCGA-CRC (n = 644) vs. normal tissue (n = 51) and its expression levels were elevated in KRAS mutations (p < 0.0001). In vitro, HOXB9 overexpression increased cell proliferation (p < 0.001) and upregulated the mRNA expression of EMT markers (VIM, CDH2, ZEB1, ZEB2, SNAI1 and SNAI2) while downregulated CDH1, (p < 0.05 for all comparisons). Conversely, HOXB9 silencing disrupted cell growth (p < 0.0001). High HOXB9 expression (HR = 3.82, 95% CI: 1.59–9.2, p = 0.003) was independently associated with worse OS in CRLM-HOXB9-expressing patients after liver resection. In conclusion, HOXB9 may be associated with worse OS in CRLM and may promote CRC progression, whereas HOXB9 silencing may inhibit CRC growth.     

Aberrant Methylation of 20 miRNA Genes Specifically Involved in Various Steps of Ovarian Carcinoma Spread: From Primary Tumors to Peritoneal Macroscopic Metastases

Our work aimed to differentiate 20 aberrantly methylated miRNA genes that participate at different stages of development and metastasis of ovarian carcinoma (OvCa) using methylation-specific qPCR in a representative set of clinical samples: 102 primary tumors without and with metastases (to lymph nodes, peritoneum, or distant organs) and 30 peritoneal macroscopic metastases (PMM). Thirteen miRNA genes (MIR107, MIR124-2, MIR124-3, MIR125B-1, MIR127, MIR129-2, MIR130B, MIR132, MIR193A, MIR339, MIR34B/C, MIR9-1, and MIR9-3) were hypermethylated already at the early stages of OvCa, while hypermethylation of MIR1258, MIR137, MIR203A, and MIR375 was pronounced in metastatic tumors, and MIR148A showed high methylation levels specifically in PMM. We confirmed the significant relationship between methylation and expression levels for 11 out of 12 miRNAs analyzed by qRT-PCR. Moreover, expression levels of six miRNAs were significantly decreased in metastatic tumors in comparison with nonmetastatic ones, and downregulation of miR-203a-3p was the most significant. We revealed an inverse relationship between expression levels of miR-203a-3p and those of ZEB1 and ZEB2 genes, which are EMT drivers. We also identified three miRNA genes (MIR148A, MIR9-1, and MIR193A) that likely regulate EMT–MET reversion in the colonization of PMM. According to the Kaplan–Meier analysis, hypermethylation of several examined miRNA genes was associated with poorer overall survival of OvCa patients, and high methylation levels of MIR130B and MIR9-1 were related to the greatest relative risk of death.     

Upregulation of TLRs and IL-6 as a Marker in Human Colorectal Cancer

Toll-like receptors (TLRs) not only form an important part of the innate immune system but also serve to activate the adaptive immune system in response to cancer. Real-time PCR; immunohistochemical stain and Western blotting analyses were performed to clarify molecular alterations in colorectal cancer (CRC) patients. We identified Toll-like receptor 1 (TLR1), TLR2, TLR4 and TLR8 gene expression levels and downstream gene, i.e., interleukin-6 (IL-6), IL-8, interferon-α (IFN-α) and myeloid differentiation primary-response protein-88 (MyD88), expression levels in CRC patients and in cancer cell lines. CRC tissues have higher TLR1, TLR2, TLR4, TLR8, IL-6 and IL-8 gene expression levels than do the normal colon mucosa (p < 0.05). TLR2 expression varied in different cell types (mucosa and lymphocytes). There was no difference in the MyD88 and IFN-α gene expression levels between cancerous and normal colon mucosa. CRC patients had higher levels of IL-6 (p = 0.002) and IL-8 (p = 0.038) expression than healthy volunteers did; and higher IL-6 and IL-8 expression was also found to signify a higher risk of recurrence. CL075 (3M002) treatments can reduce the production of IL-8 in different cancer cell lines. The signaling pathway of TLRs in cancer tissue is different from that in normal cells; and is MyD88-independent. Higher expression levels of TLR1, TLR2, TLR 4 and TLR 8 mRNA were related to upregulation inflammatory cytokines IL-6 and IL-8 gene expression in tissue and to the upregulation of IL-6 in blood. The concentration of IL-6 in serum can be used as an indicator of the possibility of CRC recurrence. Treatment with 3M002 can reduce IL-6 production in vitro and may prevent CRC recurrence. Our findings provide evidence that TLR1, TLR2, TLR4 and TLR8 gene expression induce downstream IL-6 and IL-8 gene expression; detection of these expression levels can serve as a CRC marker.     

B4GALNT2 Controls Sda and SLex Antigen Biosynthesis in Healthy and Cancer Human Colon

The Sd^a carbohydrate antigen and the corresponding biosynthetic enzyme B4GALNT2 are primarily expressed in human normal colonic mucosa and are down-regulated to variable degrees in colon cancer. On the other hand, the tumor associated antigen SLe^x is not detected in the healthy colon and is upregulated in colon cancer. High level of B4GALNT2 gene expression appears to be a good marker of prognosis in colon cancer; however, the molecular mechanisms regulating these carbohydrate antigens’ expression are still poorly understood. We review here the most recent progress made towards understanding this balanced expression of blood group carbohydrate epitopes Sd^a and SLe^x. In particular in recent years, we have attained a better understanding of genetic and epigenetic regulation of the B4GALNT2 gene and of the subcellular fate of B4GALNT2 isoforms.     

Expression and Roles of Individual HIF Prolyl 4-Hydroxylase Isoenzymes in the Regulation of the Hypoxia Response Pathway along the Murine Gastrointestinal Epithelium

The HIF prolyl 4-hydroxylases (HIF-P4H) control hypoxia-inducible factor (HIF), a powerful mechanism regulating cellular adaptation to decreased oxygenation. The gastrointestinal epithelium subsists in “physiological hypoxia” and should therefore have an especially well-designed control over this adaptation. Thus, we assessed the absolute mRNA expression levels of the HIF pathway components, Hif1a, HIF2a, Hif-p4h-1, 2 and 3 and factor inhibiting HIF (Fih1) in murine jejunum, caecum and colon epithelium using droplet digital PCR. We found a higher expression of all these genes towards the distal end of the gastrointestinal tract. We detected mRNA for Hif-p4h-1, 2 and 3 in all parts of the gastrointestinal tract. Hif-p4h-2 had significantly higher expression levels compared to Hif-p4h-1 and 3 in colon and caecum epithelium. To test the roles each HIF-P4H isoform plays in the gut epithelium, we measured the gene expression of classical HIF target genes in Hif-p4h-1^−/−, Hif-p4h-2 hypomorph and Hif-p4h-3^−/− mice. Only Hif-p4h-2 hypomorphism led to an upregulation of HIF target genes, confirming a predominant role of HIF-P4H-2. However, the abundance of Hif-p4h-1 and 3 expression in the gastrointestinal epithelium implies that these isoforms may have specific functions as well. Thus, the development of selective inhibitors might be useful for diverging therapeutic needs.     

Stress Modifies the Expression of Glucocorticoid-Responsive Genes by Acting at Epigenetic Levels in the Rat Prefrontal Cortex: Modulatory Activity of Lurasidone

Epigenetics is one of the mechanisms by which environmental factors can alter brain function and may contribute to central nervous system disorders. Alterations of DNA methylation and miRNA expression can induce long-lasting changes in neurobiological processes. Hence, we investigated the effect of chronic stress, by employing the chronic mild stress (CMS) and the chronic restraint stress protocol, in adult male rats, on the glucocorticoid receptor (GR) function. We focused on DNA methylation specifically in the proximity of the glucocorticoid responsive element (GRE) of the GR responsive genes Gadd45β, Sgk1, and Gilz and on selected miRNA targeting these genes. Moreover, we assessed the role of the antipsychotic lurasidone in modulating these alterations. Chronic stress downregulated Gadd45β and Gilz gene expression and lurasidone normalized the Gadd45β modification. At the epigenetic level, CMS induced hypermethylation of the GRE of Gadd45β gene, an effect prevented by lurasidone treatment. These stress-induced alterations were still present even after a period of rest from stress, indicating the enduring nature of such changes. However, the contribution of miRNA to the alterations in gene expression was moderate in our experimental conditions. Our results demonstrated that chronic stress mainly affects Gadd45β expression and methylation, effects that are prolonged over time, suggesting that stress leads to changes in DNA methylation that last also after the cessation of stress procedure, and that lurasidone is a modifier of such mechanisms.     

Expression Pattern of Purinergic Signaling Components in Colorectal Cancer Cells and Differential Cellular Outcomes Induced by Extracellular ATP and Adenosine

The purine nucleotide adenosine triphosphate (ATP) is known for its fundamental role in cellular bioenergetics. However, in the last decades, different works have described emerging functions for ATP, such as that of a danger signaling molecule acting in the extracellular space on both tumor and stromal compartments. Beside its role in immune cell signaling, several studies have shown that high concentrations of extracellular ATP can directly or indirectly act on cancer cells. Accordingly, it has been reported that purinergic receptors are widely expressed in tumor cells. However, their expression pattern is often associated with contradictory cellular outcomes. In this work, we first investigated gene expression profiles through “RNA-Sequencing” (RNA Seq) technology in four colorectal cancer (CRC) cell lines (HT29, LS513, LS174T, HCT116). Our results demonstrate that CRC cells mostly express the A2B, P2X4, P2Y1, P2Y2 and P2Y11 purinergic receptors. Among these, the P2Y1 and P2Y2 coding genes are markedly overexpressed in all CRC cells compared to the HCEC-1CT normal-like colonic cells. We then explored the cellular outcomes induced by extracellular ATP and adenosine. Our results show that in terms of cell death induction extracellular ATP is consistently more active than adenosine against CRC, while neither compound affected normal-like colonic cell survival. Intriguingly, while for the P2Y2 receptor pharmacological inhibition completely abolished the rise in cytoplasmic Ca²⁺ observed after ATP exposure in all CRC cell lines, Ca²⁺ mobilization only impacted the cellular outcome for HT29. In contrast, non-selective phosphodiesterase inhibition completely abolished the effects of extracellular ATP on CRC cells, suggesting that cAMP and/or cGMP levels might determine cellular outcome. Altogether, our study provides novel insights into the characterization of purinergic signaling in CRC.     

Using VBIM Technique to Discover ARMC4/ODAD2 as a Novel Negative Regulator of NF-κB and a New Tumor Suppressor in Colorectal Cancer

Since nuclear factor (NF) κB plays pivotal roles in inflammation and cancer, understanding its regulation holds great promise for disease therapy. Using the powerful validation-based insertional mutagenesis (VBIM) technique established by us previously, we discovered armadillo repeat-containing protein 4 (ARMC4)/outer dynein arm docking complex subunit 2 (ODAD2), a rarely studied protein known to date, as a novel negative regulator of NF-κB in colorectal cancer (CRC). High expression of ARMC4 downregulated the expression of NF-κB-dependent genes, dramatically reduced NF-κB activity, cellular proliferation, anchorage-independent growth, and migratory ability in vitro, and significantly decreased xenograft tumor growth in vivo. Co-immunoprecipitation experiments demonstrated that ARMC4 forms a complex with NF-κB. Importantly, the lower ARMC4 expression in patient tumors than normal tissues indicates its potential tumor suppressor function in CRC. Collectively, we uncovered a completely new facet of ARMC4 function by identifying it as a novel NF-κB negative regulator, thus uncovering ARMC4 as a potential new therapeutic target in CRC.     

One-Step Generation of Multiple Gene-Edited Pigs by Electroporation of the CRISPR/Cas9 System into Zygotes to Reduce Xenoantigen Biosynthesis

Xenoantigens cause hyperacute rejection and limit the success of interspecific xenografts. Therefore, genes involved in xenoantigen biosynthesis, such as GGTA1, CMAH, and B4GALNT2, are key targets to improve the outcomes of xenotransplantation. In this study, we introduced a CRISPR/Cas9 system simultaneously targeting GGTA1, CMAH, and B4GALNT2 into in vitro-fertilized zygotes using electroporation for the one-step generation of multiple gene-edited pigs without xenoantigens. First, we optimized the combination of guide RNAs (gRNAs) targeting GGTA1 and CMAH with respect to gene editing efficiency in zygotes, and transferred electroporated embryos with the optimized gRNAs and Cas9 into recipient gilts. Next, we optimized the Cas9 protein concentration with respect to the gene editing efficiency when GGTA1, CMAH, and B4GALNT2 were targeted simultaneously, and generated gene-edited pigs using the optimized conditions. We achieved the one-step generation of GGTA1/CMAH double-edited pigs and GGTA1/CMAH/B4GALNT2 triple-edited pigs. Immunohistological analyses demonstrated the downregulation of xenoantigens; however, these multiple gene-edited pigs were genetic mosaics that failed to knock out some xenoantigens. Although mosaicism should be resolved, the electroporation technique could become a primary method for the one-step generation of multiple gene modifications in pigs aimed at improving pig-to-human xenotransplantation.     

A Presurgical Study of Curcumin Combined with Anthocyanin Supplements in Patients with Colorectal Adenomatous Polyps

Adenomatous polyps are precancerous lesions associated with a higher risk of colorectal cancer (CRC). Curcumin and anthocyanins have shown promising CRC-preventive activity in preclinical and epidemiological studies. The objective of this window-of-opportunity, proof-of principle trial was to evaluate the effect of curcumin combined with anthocyanin supplements on tissue biomarkers of colorectal adenomatous polyps. Eligible patients received either anthocyanin and curcumin supplementation or related matching placebo for 4–6 weeks before polyp removal. Adenomatous polyps and adjacent tissue biopsies were collected at baseline and after supplementation for immunohistochemical assessment of β-catenin, NF-kappa B (NF-κB), Ki-67, P53, and dysplasia. No differences were observed in baseline biomarker expression between normal and dysplastic tissues. The combination of anthocyanins and curcumin resulted in a significant borderline reduction of NF-κB immunohistochemistry (IHC) expression in adenoma tissue (geometric mean ratio (GMR): 0.72; 95% confidence interval (CI): 0.51–1.00; p-value: 0.05) and a trend to a reduction of Ki-67 (GMR: 0.73; 95% CI: 0.50–1.08; p-value: 0.11). No significant modulation of biomarkers in normal adjacent mucosa was observed. We concluded that the combined supplementation of anthocyanins and curcumin seems to lead to a potentially favorable modulation of tissue biomarkers of inflammation and proliferation in colon adenomas.     

The Profile of MicroRNA Expression and Potential Role in the Regulation of Drug-Resistant Genes in Cisplatin- and Paclitaxel-Resistant Ovarian Cancer Cell Lines

Ovarian cancer is the most lethal gynecological malignancy. The high mortality results from late diagnosis and the development of drug resistance. Drug resistance results from changes in the expression of different drug-resistance genes that may be regulated miRNA. The main aim of our study was to detect changes in miRNA expression levels in two cisplatin (CIS) and two paclitaxel (PAC)—resistant variants of the A2780 drug-sensitive ovarian cancer cell line—by miRNA microarray. The next goal was to identify miRNAs responsible for the regulation of drug-resistance genes. We observed changes in the expression of 46 miRNA that may be related to drug resistance. The overexpression of miR-125b-5p, miR-99a-5p, miR-296-3p, and miR-887-3p and downregulation of miR-218-5p, miR-221-3p, and miR-222-3p was observed in both CIS-resistant cell lines. In both PAC-resistant cell lines, we observed the upregulation of miR-221-3p, miR-222-3p, and miR-4485, and decreased expression of miR-551b-3p, miR-551b-5p, and miR-218-5p. Analysis of targets suggest that expression of important drug-resistant genes like protein Tyrosine Phosphatase Receptor Type K (PTPRK), receptor tyrosine kinase—EPHA7, Semaphorin 3A (SEMA3A), or the ATP-binding cassette subfamily B member 1 gene (ABCB1) can be regulated by miRNA.     

Significant Overexpression of DVL1 in Taiwanese Colorectal Cancer Patients with Liver Metastasis

Undetected micrometastasis plays a key role in the metastasis of cancer in colorectal cancer (CRC) patients. The aim of this study is to identify a biomarker of CRC patients with liver metastasis through the detection of circulating tumor cells (CTCs). Microarray and bioinformatics analysis of 10 CRC cancer tissue specimens compared with normal adjacent tissues revealed that 31 genes were up-regulated (gene expression ratio of cancer tissue to paired normal tissue > 2) in the cancer patients. We used a weighted enzymatic chip array (WEnCA) including 31 prognosis-related genes to investigate CTCs in 214 postoperative stage I–III CRC patients and to analyze the correlation between gene expression and clinico-pathological parameters. We employed the immunohistochemistry (IHC) method with polyclonal mouse antibody against DVL1 to detect DVL1 expression in 60 CRC patients. CRC liver metastasis occurred in 19.16% (41/214) of the patients. Using univariate analysis and multivariate proportional hazards regression analysis, we found that DVL1 mRNA overexpression had a significant, independent predictive value for liver metastasis in CRC patients (OR: 5.764; 95% CI: 2.588–12.837; p < 0.0001 on univariate analysis; OR: 3.768; 95% CI: 1.469–9.665; p = 0.006 on multivariate analysis). IHC staining of the immunoreactivity of DVL1 showed that DVL1 was localized in the cytoplasm of CRC cells. High expression of DVL1 was observed in 55% (33/60) of CRC tumor specimens and was associated significantly with tumor depth, perineural invasion and liver metastasis status (all p < 0.05). Our experimental results demonstrated that DVL1 is significantly overexpressed in CRC patients with liver metastasis, leading us to conclude that DVL1 could be a potential prognostic and predictive marker for CRC patients.