Identification of RNF213 as a potential suppressor of local invasion in intrahepatic cholangiocarcinoma.

Intrahepatic cholangiocarcinoma (ICC) is a lethal cancer with poor survival especially when it spreads. The histopathology of its rare intraductal papillary neoplasm of the bile duct type (IPNB) characteristically shows cancer cells originating within the confined bile duct space. These cells eventually invade and infiltrate the nearby liver tissues, making it a good model to study the mechanism of local invasion, which is the earliest step of metastasis. To discover potential suppressor genes of local invasion in ICC, we analyzed the somatic mutation profiles and performed clonal evolution analyses of the 11 pairs of macrodissected IPNB tissues with local invasion (LI-IPNB) and IPNB tissues without local invasion from the same patients. We identified a protein-truncating variant (PTV) in an E3 ubiquitin ligase, RNF213 (c.6967C>T; p.Gln2323X; chr17: 78,319,102 [hg19], exon 29), as the most common PTV event in LI-IPNB samples (4/11 patients). Knockdown of RNF213 in HuCCT1 and YSCCC cells showed increased migration and invasion, and reduced vasculogenic mimicry, but maintained normal proliferation. Transcriptomic analysis of the RNF213-knockdown versus control cells was then performed in the HuCCT1, YSCCC, and KKU-100 cells. Gene Ontology (GO) enrichment analysis of the common differentially expressed genes revealed significantly altered cytokine and oxidoreductase-oxidizing metal ions activities, as confirmed by western blotting. Gene Set Enrichment Analysis (GSEA) identified the most enriched pathways being oxidative phosphorylation, fatty acid metabolism, reactive oxygen species, adipogenesis, and angiogenesis. In sum, loss of function of RNF213 is a common genetic alteration in LI-IPNB tissues. RNF213 knockdown leads to increased migration and invasion of ICC cells, potentially through malfunctions of the pathways related inflammation, and energy metabolisms.

PGE2-EP2/EP4 signaling elicits immunosuppression by driving the mregDC-Treg axis in inflammatory tumor microenvironment.

Active inflammation generally promotes immune activation. However, in the tumor microenvironment (TME), active inflammation occurs in parallel with immunosuppression, and both contribute to tumor growth. Why inflammation does not lead to immune activation in TME remains unclear. In this study, using the immune checkpoint inhibitor-insensitive mouse cancer model and single-cell RNA sequencing, we show that PGE2-EP2/EP4 signaling simultaneously promotes active inflammation by inducing expression of the NF-κB genes in myeloid cells and elicits immunosuppression by driving the mregDC (mature DC enriched in immunoregulatory molecules)-Treg (regulatory T cell) axis for Treg recruitment and activation in the tumor. Importantly, the EP2/EP4 expression level is strongly correlated with the gene signatures of both active inflammation and the mregDC-Treg axis and has significant prognosis value in various human cancers. Thus, PGE2-EP2/EP4 signaling functions as the key regulatory node linking active inflammation and immunosuppression in TME, which can be targeted by EP2 and EP4 antagonists for cancer therapeutics.

Capturing tumour heterogeneity in pre- and post-chemotherapy colorectal cancer ascites-derived cells using single-cell RNA-sequencing.

Malignant ascites is an abnormal accumulation of fluid within the peritoneal cavity, caused by metastasis of several types of cancers, including colorectal cancer (CRC). Cancer cells in ascites reflect poor prognosis and serve as a good specimen to study tumour heterogeneity, as they represent a collection of multiple metastatic sites in the peritoneum. In the present study, we have employed single-cell RNA-sequencing (scRNA-seq) to explore and characterise ascites-derived cells from a CRC patient. The samples were prepared using mechanical and enzymatic dissociations, and obtained before and after a chemotherapy treatment. Unbiased clustering of 19,653 cells from four samples reveals 14 subclusters with unique transcriptomic patterns in four major cell types: epithelial cells, myeloid cells, fibroblasts, and lymphocytes. Interestingly, the percentages of cells recovered from different cell types appeared to be influenced by the preparation protocols, with more than 90% reduction in the number of myeloid cells recovered by enzymatic preparation. Analysis of epithelial cell subpopulations unveiled only three out of eleven subpopulations with clear contraction after the treatment, suggesting that the majority of the heterogeneous ascites-derived cells were resistant to the treatment, potentially reflecting the poor treatment outcome observed in the patient. Overall, our study showcases highly heterogeneous cancer subpopulations at single-cell resolution, which respond differently to a particular chemotherapy treatment. All in all, this work highlights the potential benefit of single-cell analyses in planning appropriate treatments and real-time monitoring of therapeutic response in cancer patients through routinely discarded ascites samples.

Overexpression of Pyruvate Carboxylase Is Correlated With Colorectal Cancer Progression and Supports Growth of Invasive Colon Cancer HT-29 Cell Line.

BACKGROUND/AIM: Pyruvate carboxylase (PC) is a major anaplerotic enzyme for generating oxaloacetate for the TCA cycle and also a key enzyme in gluconeogenesis, de novo fatty acid and amino acid synthesis in normal cells. Recent studies have identified PC overexpression in different cancers, such as breast and lung. However, the involvement of PC in colorectal cancer (CRC) is unclear. Our purpose was to investigate the PC expression levels and its correlations with potentially relevant clinical-pathological parameters in CRC. MATERIALS AND METHODS: PC expression levels in tissues from 60 Thai CRC patients were investigated by immunohistochemistry while a clonogenic assay was performed for determining cell growth of HT-29 cells with PC knockdown. RESULTS: Our results showed for the first time that high PC expression levels were significantly correlated with late stage of the cancer, perineural invasion and lymph node metastasis. The overexpression of PC was also significantly associated with poor overall and disease-free survival times of CRC patients. In addition, suppression of cancer cell growth was found in PC-deficient cell lines using CRISPR-Cas9. CONCLUSION: The overexpression levels of PC were correlated with CRC progression and survival times. Therefore, PC might serve as a potential clinical prognostic marker for colorectal cancer.

Overexpression of Holocarboxylase Synthetase Predicts Lymph Node Metastasis and Unfavorable Prognosis in Breast Cancer.

BACKGROUND/AIM: Holocarboxylase synthetase (HLCS) catalyzes the specific attachment of biotin onto biotin-dependent carboxylases (BDCs) which play important roles in intermediary metabolism. Previous studies show that BDCs are overexpressed in many cancer types. However, expression of HLCS in cancerous tissues has not been reported. MATERIALS AND METHODS: Immunohistochemistry was used to investigate HLCS expression in breast tissue obtained from 65 Thai patients, and the correlation between its expression and key clinical-pathological parameters was assessed. The role of HLCS in supporting invasion was investigated in HLCS-knockdown MCF-7 cells. RESULTS: Overexpression of HLCS was significantly associated with metastasis of breast cancer cells to other lymph nodes but not the sentinel and axillary lymph nodes - a finding supported in cellular invasion assays using HLCS knockdown cells. Furthermore, overexpression of HLCS reduced survival time of patients with breast cancer. CONCLUSION: HLCS appears to be a prognostic marker for patients with breast cancer.

IGSF3 mutation identified in patient with severe COPD alters cell function and motility.

Cigarette smoking (CS) and genetic susceptibility determine the risk for development, progression, and severity of chronic obstructive pulmonary diseases (COPD). We posited that an incidental balanced reciprocal chromosomal translocation was linked to a patient's risk of severe COPD. We determined that 46,XX,t(1;4)(p13.1;q34.3) caused a breakpoint in the immunoglobulin superfamily member 3 (IGSF3) gene, with markedly decreased expression. Examination of COPDGene cohort identified 14 IGSF3 SNPs, of which rs1414272 and rs12066192 were directly and rs6703791 inversely associated with COPD severity, including COPD exacerbations. We confirmed that IGSF3 is a tetraspanin-interacting protein that colocalized with CD9 and integrin B1 in tetraspanin-enriched domains. IGSF3-deficient patient-derived lymphoblastoids exhibited multiple alterations in gene expression, especially in the unfolded protein response and ceramide pathways. IGSF3-deficient lymphoblastoids had high ceramide and sphingosine-1 phosphate but low glycosphingolipids and ganglioside levels, and they were less apoptotic and more adherent, with marked changes in multiple TNFRSF molecules. Similarly, IGSF3 knockdown increased ceramide in lung structural cells, rendering them more adherent, with impaired wound repair and weakened barrier function. These findings suggest that, by maintaining sphingolipid and membrane receptor homeostasis, IGSF3 is required for cell mobility-mediated lung injury repair. IGSF3 deficiency may increase susceptibility to CS-induced lung injury in COPD.

c-Myc directly targets an over-expression of pyruvate carboxylase in highly invasive breast cancer.

Here we showed that the c-Myc oncogene is responsible for overexpression of pyruvate carboxylase (PC) in highly invasive MDA-MB-231 cells. Pharmacological inhibition of c-Myc activity with 10074-G5 compound, resulted in a marked reduction of PC mRNA and protein, concomitant with reduced cell growth, migration and invasion. This growth inhibition but not migration and invasion can be partly restored by overexpression of PC, indicating that PC is a c-Myc-regulated pro-proliferating enzyme. Analysis of chromatin immunoprecipitation sequencing of c-Myc bound promoters revealed that c-Myc binds to two canonical c-Myc binding sites, locating at nucleotides -417 to -407 and -301 to -291 in the P2 promoter of human PC gene. Mutation of either c-Myc binding site in the P2 promoter-luciferase construct resulted in 50-60% decrease in luciferase activity while double mutation of c-Myc binding sites further decreased the luciferase activity in MDA-MB-231 cells. Overexpression of c-Myc in HEK293T cells that have no endogenous c-Myc resulted in 250-fold increase in luciferase activity. Mutation of either E-boxes lowered luciferase activity by 50% and 25%, respectively while double mutation of both sites abolished the c-Myc transactivation response. An electrophoretic mobility shift assay using nuclear proteins from MDA-MB-231 confirmed binding of c-Myc to both c-Myc binding sites in the P2 promoter. Bioinformatic analysis of publicly available transcriptomes from the cancer genome atlas (TCGA) dataset revealed an association between expression of c-Myc and PC in primary breast, as well as in lung and colon cancer tissues, suggesting that overexpression of PC is deregulated by c-Myc in these cancers.

Dysregulated microRNA expression profiles in cholangiocarcinoma cell-derived exosomes.

AIM: Cholangiocarcinoma (CCA) is a malignant tumor of bile duct epithelial cells. The prognosis of CCA is poor due to lack of effective therapeutic targets and detection at an advanced stage. Exosomes are secreted nano-sized vesicles and contribute to the malignancy of several cancers via transferring their miRNAs between cells. Thus, exosomal miRNAs may serve as new therapeutic targets and potential biomarkers for CCA. MAIN METHODS: Exosomes were isolated from three different CCA cell lines and normal human cholangiocyte cells, followed by miRNA profiling analysis. Potential role of dysregulated miRNA was investigated by knockdown experiment. KEY FINDINGS: We found that 38 and 460 miRNAs in CCA exosomes were significantly up- and down-regulated, respectively. Of these differentially expressed miRNAs, the hsa-miR-205-5p and miR-200 family members were markedly up-regulated for 600-1500 folds, whereas the miR-199 family members and their clustered miRNA, hsa-miR-214-3p, were down-regulated for 1000-2000 folds. The expression patterns of these representative exosomal miRNAs were similar to those observed in all types of CCA cells. The target genes of the top ten most up- and down-regulated miRNAs are significantly associated with well-characterized cancer-related pathways. Consistently, knockdown of the most up-regulated miRNA, miR-205-5p, reduced KKU-M213 cell invasion and migration. SIGNIFICANCE: We have demonstrated the distinct miRNA signatures in exosomes released from CCA cells, compared to normal human cholangiocyte cells. These exosomal miRNAs may have the potential to be novel therapeutic targets and biomarkers for CCA.

Emergence of Intrahepatic Cholangiocarcinoma: How High-Throughput Technologies Expedite the Solutions for a Rare Cancer Type.

Intrahepatic cholangiocarcinoma (ICC) is the cancer of the intrahepatic bile ducts, and together with hepatocellular carcinoma (HCC), constitute the majority of primary liver cancers. ICC is a rare disorder as its overall incidence is < 1/100,000 in the United States and Europe. However, it shows much higher incidence in particular geographical regions, such as northeastern Thailand, where liver fluke infection is the most common risk factor of ICC. Since the early stages of ICC are often asymptomatic, the patients are usually diagnosed at advanced stages with no effective treatments available, leading to the high mortality rate. In addition, unclear genetic mechanisms, heterogeneous nature, and various etiologies complicate the development of new efficient treatments. Recently, a number of studies have employed high-throughput approaches, including next-generation sequencing and mass spectrometry, in order to understand ICC in different biological aspects. In general, the majority of recurrent genetic alterations identified in ICC are enriched in known tumor suppressor genes and oncogenes, such as mutations in TP53, KRAS, BAP1, ARID1A, IDH1, IDH2, and novel FGFR2 fusion genes. Yet, there are no major driver genes with immediate clinical solutions characterized. Interestingly, recent studies utilized multi-omics data to classify ICC into two main subgroups, one with immune response genes as the main driving factor, while another is enriched with driver mutations in the genes associated with epigenetic regulations, such as IDH1 and IDH2. The two subgroups also show different hypermethylation patterns in the promoter regions. Additionally, the immune response induced by host-pathogen interactions, i.e., liver fluke infection, may further stimulate tumor growth through alterations of the tumor microenvironment. For in-depth functional studies, although many ICC cell lines have been globally established, these homogeneous cell lines may not fully explain the highly heterogeneous genetic contents of this disorder. Therefore, the advent of patient-derived xenograft and 3D patient-derived organoids as new disease models together with the understanding of evolution and genetic alterations of tumor cells at the single-cell resolution will likely become the main focus to fill the current translational research gaps of ICC in the future.

Bridging the gap between clinicians and systems biologists: from network biology to translational biomedical research.

With the wealth of data accumulated from completely sequenced genomes and other high-throughput experiments, global studies of biological systems, by simultaneously investigating multiple biological entities (e.g. genes, transcripts, proteins), has become a routine. Network representation is frequently used to capture the presence of these molecules as well as their relationship. Network biology has been widely used in molecular biology and genetics, where several network properties have been shown to be functionally important. Here, we discuss how such methodology can be useful to translational biomedical research, where scientists traditionally focus on one or a small set of genes, diseases, and drug candidates at any one time. We first give an overview of network representation frequently used in biology: what nodes and edges represent, and review its application in preclinical research to date. Using cancer as an example, we review how network biology can facilitate system-wide approaches to identify targeted small molecule inhibitors. These types of inhibitors have the potential to be more specific, resulting in high efficacy treatments with less side effects, compared to the conventional treatments such as chemotherapy. Global analysis may provide better insight into the overall picture of human diseases, as well as identify previously overlooked problems, leading to rapid advances in medicine. From the clinicians' point of view, it is necessary to bridge the gap between theoretical network biology and practical biomedical research, in order to improve the diagnosis, prevention, and treatment of the world's major diseases.

MicroRNAs and oncogenic transcriptional regulatory networks controlling metabolic reprogramming in cancers.

Altered cellular metabolism is a fundamental adaptation of cancer during rapid proliferation as a result of growth factor overstimulation. We review different pathways involving metabolic alterations in cancers including aerobic glycolysis, pentose phosphate pathway, de novo fatty acid synthesis, and serine and glycine metabolism. Although oncoproteins, c-MYC, HIF1α and p53 are the major drivers of this metabolic reprogramming, post-transcriptional regulation by microRNAs (miR) also plays an important role in finely adjusting the requirement of the key metabolic enzymes underlying this metabolic reprogramming. We also combine the literature data on the miRNAs that potentially regulate 40 metabolic enzymes responsible for metabolic reprogramming in cancers, with additional miRs from computational prediction. Our analyses show that: (1) a metabolic enzyme is frequently regulated by multiple miRs, (2) confidence scores from prediction algorithms might be useful to help narrow down functional miR-mRNA interaction, which might be worth further experimental validation. By combining known and predicted interactions of oncogenic transcription factors (TFs) (c-MYC, HIF1α and p53), sterol regulatory element binding protein 1 (SREBP1), 40 metabolic enzymes, and regulatory miRs we have established one of the first reference maps for miRs and oncogenic TFs that regulate metabolic reprogramming in cancers. The combined network shows that glycolytic enzymes are linked to miRs via p53, c-MYC, HIF1α, whereas the genes in serine, glycine and one carbon metabolism are regulated via the c-MYC, as well as other regulatory organization that cannot be observed by investigating individual miRs, TFs, and target genes.

Genomic repertoires of DNA-binding transcription factors across the tree of life.

Sequence-specific transcription factors (TFs) are important to genetic regulation in all organisms because they recognize and directly bind to regulatory regions on DNA. Here, we survey and summarize the TF resources available. We outline the organisms for which TF annotation is provided, and discuss the criteria and methods used to annotate TFs by different databases. By using genomic TF repertoires from ∼700 genomes across the tree of life, covering Bacteria, Archaea and Eukaryota, we review TF abundance with respect to the number of genes, as well as their structural complexity in diverse lineages. While typical eukaryotic TFs are longer than the average eukaryotic proteins, the inverse is true for prokaryotes. Only in eukaryotes does the same family of DNA-binding domain (DBD) occur multiple times within one polypeptide chain. This potentially increases the length and diversity of DNA-recognition sequence by reusing DBDs from the same family. We examined the increase in TF abundance with the number of genes in genomes, using the largest set of prokaryotic and eukaryotic genomes to date. As pointed out before, prokaryotic TFs increase faster than linearly. We further observe a similar relationship in eukaryotic genomes with a slower increase in TFs.