The androgen-regulated Calcium-Activated Nucleotidase 1 (CANT1) is commonly overexpressed in prostate cancer and is tumor-biologically relevant in vitro

WTAP-mediated N6-methyladenine Modification of circEEF2 Promotes Lung Adenocarcinoma Tumorigenesis by Stabilizing CANT1 in an IGF2BP2-dependent Manner

N6-methyladenosine (m6A) is a common posttranscriptional RNA modification and plays an important role in cancer biology. Circular RNAs (circRNAs) are also reported to participate in lung adenocarcinoma (LUAD) progression. Here, we aimed to investigate the functions of Wilms tumor 1-associating protein (WTAP) methyltransferase and circEEF2 in LUAD cell tumorigenesis, and probe whether circEEF2 functioned through WTAP-induced m6A modification and its potential mechanisms. Functional analyses were conducted by tube formation, sphere formation, 5-ethynyl-2'-deoxyuridine (EdU), flow cytometry, and transwell assays in vitro as well as tumor formation experiments in mice, respectively. The N6-methyladenine (m6A) modification in circEEF2 mRNA was determined by RNA immunoprecipitation (Me-RIP) assay. The interaction between IGF2BP2 (Insulin Like Growth Factor 2 MRNA-Binding Protein 2) and circEEF2 or Calcium-activated nucleotidase 1 (CANT1) mRNA was confirmed using RIP assay. LUAD tissues and cells showed high circEEF2 expression, and the deficiency of circEEF2 suppressed LUAD cell angiogenesis, stemness, proliferation, migration, and invasion. WTAP induced circEEF2 m6A modification. WTAP silencing repressed the oncogenic phenotypes of LUAD cells via stabilizing circEEF2 in an m6A-dependent manner. IGF2BP2 interacted with circEEF2 and CANT1, and WTAP and circEEF2 could regulate CANT1 expression through IGF2BP2. The inhibition of LUAD cell oncogenic phenotypes caused by circEEF2 deficiency was abolished by CANT1 overexpression. In addition, WTAP silencing impeded LUAD growth via modulating circEEF2 and CANT1 in vivo. WTAP-mediated m6A modification of circEEF2 promotes lung adenocarcinoma growth and tumorigenesis by stabilizing CANT1 through IGF2BP2.

Cant1 Affects Cartilage Proteoglycan Properties: Aggrecan and Decorin Characterization in a Mouse Model of Desbuquois Dysplasia Type 1

Desbuquois dysplasia type 1 (DBQD1) is a recessive chondrodysplasia caused by mutations in the CANT1 gene, encoding for the Golgi Calcium-Activated Nucleotidase 1 (CANT1). The enzyme hydrolyzes UDP, the by-product of glycosyltransferase reactions, but it might play other roles in different cell types. Using a Cant1 knock-out mouse, we demonstrated that CANT1 is crucial for glycosaminoglycan (GAG) synthesis; however, its impact on the biochemical properties of cartilage proteoglycans remains unknown. Thus, in this work, we characterized decorin and aggrecan from primary chondrocyte cultures and cartilage biopsies of mutant mice at post-natal day 4 by Western blots and further investigated their distribution in the cartilage extracellular matrix (ECM) by immunohistochemistry. We demonstrated that the GAG synthesis defect caused by CANT1 impairment led to the synthesis and secretion of proteoglycans with shorter GAG chains compared with wild-type animals. However, this alteration did not result in the synthesis and secretion of decorin and aggrecan in the unglycanated form. Interestingly, the defect was not cartilage-specific since also skin decorin showed a reduced hydrodynamic size. Finally, immunohistochemical studies in epiphyseal sections of mutant mice demonstrated that the proteoglycan structural defect moderately affected decorin distribution in the ECM.

Co-expression in tissue-specific gene networks links genes in cancer-susceptibility loci to known somatic driver genes

BACKGROUND

The genetic background of cancer remains complex and challenging to integrate. Many somatic mutations within genes are known to cause and drive cancer, while genome-wide association studies (GWAS) of cancer have revealed many germline risk factors associated with cancer. However, the overlap between known somatic driver genes and positional candidate genes from GWAS loci is surprisingly small. We hypothesised that genes from multiple independent cancer GWAS loci should show tissue-specific co-regulation patterns that converge on cancer-specific driver genes.

RESULTS

We studied recent well-powered GWAS of breast, prostate, colorectal and skin cancer by estimating co-expression between genes and subsequently prioritising genes that show significant co-expression with genes mapping within susceptibility loci from cancer GWAS. We observed that the prioritised genes were strongly enriched for cancer drivers defined by COSMIC, IntOGen and Dietlein et al. The enrichment of known cancer driver genes was most significant when using co-expression networks derived from non-cancer samples of the relevant tissue of origin.

CONCLUSION

We show how genes within risk loci identified by cancer GWAS can be linked to known cancer driver genes through tissue-specific co-expression networks. This provides an important explanation for why seemingly unrelated sets of genes that harbour either germline risk factors or somatic mutations can eventually cause the same type of disease.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

Pan-cancer analysis predicts CANT1 as a potential prognostic, immunologic biomarker

BACKGROUND

CANT1, calcium-activated nucleotidase 1, was reported to be upregulated in certain tumors. However, the function mechanism of CANT1 in pan-cancer is still unclear.

METHODS

We utilized the Cancer Genome Atlas Program (TCGA) and UALCAN databases to analyze CANT1 expression at the level of mRNA, protein, and promoter methylation in pan-cancer, and the cBioportal database to study the frequency of gene changes for CANT1. Wilcoxon test was applied to discuss the correlation between CANT1 and clinicopathological features in different tumor types. The prognosis of CANT1 in pan-cancer was discussed by Cox regression. Spearman analysis was applied to discuss the relationship of CANT1 with tumor mutation burden(TMB), microsatellite instability(MSI), immune cell infiltration, and immune checkpoints. The association between CANT1 expression and drug sensitivity for pan-cancer was investigated by the GSCALite database. In addition, we collected 40 cases of lung adenocarcinoma (LUAD) and adjacent normal tissues for immunohistochemical staining and investigated the relationship between CANT1 and clinicopathology and prognosis in LUAD. Finally, the molecular pathways involved in CANT1-related genes in LUAD were analyzed by gene set enrichment analysis(GSEA).

RESULTS

The CANT1 mRNA level was significant higher in 14 tumors, and CANT1 protein level was significant higher in 7 tumors compared with normal tissues. CANT1 expression was linked with the T stage, N stage, and pathological stage in some tumors, and overexpression CANT1 was associated with adverse overall survival(OS) and disease-specific survival(DSS) in kidney renal papillary cell carcinoma(KIRP), brain lower grade glioma(LGG), and LUAD. By Spearman correlation analysis, the results showed that CANT1 had a positive correlation with TMB and MSI in bladder urothelial carcinoma(BLCA), breast invasive carcinoma(BRCA), esophageal carcinoma(ESCA), LGG, and sarcoma(SARC). Furthermore, CANT1 was related to immune cell infiltration and immune checkpoints in several cancers. Drug sensitivity analysis suggested that CANT1 was inversely linked with three drugs. Immunohistochemical staining analysis showed that CANT1 expression was higher in LUAD than in normal tissues, and was related to N stage and pathological stage. Survival curves showed that CANT1 overexpression had poor OS and DSS. Time-dependent ROC curves revealed that the 1-year, 5-year, and 10-year OS and DSS in LUAD were above 0.5. CANT1-related genes in LUAD mainly participated in the pathway of dorso ventral axis formation, small cell lung cancer, DNA replication, O-glycan biosynthesis, and cell cycle.

CONCLUSION

CANT1 is considered a potential marker for prognosis in several tumors, and a promising target for tumor immunological treatment.

Uncovering hidden cancer self-dependencies through analysis of shRNA-level dependency scores

Large-scale short hairpin RNA (shRNA) screens on well-characterized human cancer cell lines have been widely used to identify novel cancer dependencies. However, the off-target effects of shRNA reagents pose a significant challenge in the analysis of these screens. To mitigate these off-target effects, various approaches have been proposed that aggregate different shRNA viability scores targeting a gene into a single gene-level viability score. Most computational methods for discovering cancer dependencies rely on these gene-level scores. In this paper, we propose a computational method, named NBDep, to find cancer self-dependencies by directly analyzing shRNA-level dependency scores instead of gene-level scores. The NBDep algorithm begins by removing known batch effects of the shRNAs and selecting a subset of concordant shRNAs for each gene. It then uses negative binomial random effects models to statistically assess the dependency between genetic alterations and the viabilities of cell lines by incorporating all shRNA dependency scores of each gene into the model. We applied NBDep to the shRNA dependency scores available at Project DRIVE, which covers 26 different types of cancer. The proposed method identified more well-known and putative cancer genes compared to alternative gene-level approaches in pan-cancer and cancer-specific analyses. Additionally, we demonstrated that NBDep controls type-I error and outperforms statistical tests based on gene-level scores in simulation studies.

Upregulated CANT1 is correlated with poor prognosis in hepatocellular carcinoma

BACKGROUND

CANT1, as calcium-activated protein nucleotidase 1, is a kind of phosphatase. It is overexpressed in some tumors and related to poor prognosis, but few studies explore its function and carcinogenic mechanism in hepatocellular carcinoma (HCC).

METHODS

The expression of CANT1 mRNA and protein was analyzed by the Cancer Genome Atlas (TCGA) database and immunohistochemistry(IHC) staining. The relationship between CANT1 expression and clinicopathology was evaluated by various public databases. The receiver operating characteristic (ROC) curve was used to assess the diagnostic accuracy of CANT1 by the area under curve (AUC). Univariate, multivariate Cox regression and Kaplan-Meier curves were applied to evaluate the predictive value of CANT1 on the prognosis of HCC. Methsurv was used to analyze gene changes and DNA methylation, and its impact on prognosis. The enrichment analysis of DEGs associated with CANT1 revealed the biological process of CANT1 based on Gene Set Enrichment Analysis (GSEA). The relationship between immune cell infiltration level and CANT1 expression in HCC was investigated using the single-sample GSEA (ssGSEA) method and the Tumor Immune Estimation Resource (TIMER) database. Finally, the association between CANT1 and immune checkpoints and drug sensitivity was also analyzed.

RESULTS

CANT1 was highly expressed in 22 cancers, including HCC, and CANT1 overexpression in HCC was confirmed by IHC. The expression of CANT1 was correlated with clinical features, such as histologic grade. Highly expressed CANT1 caused poor overall survival (OS) of HCC patients. Univariate and multivariate regression analysis suggested that CANT1 was an independent prognostic marker. Of the 31 DNA methylation at CpG sites, three CpG sites were associated with the prognosis of HCC. GSEA indicated that CANT1 was mainly involved in the cell cycle, DNA replication, and etc. Moreover, CANT1 expression was correlated with immune cell infiltration and independently associated with the prognosis of HCC patients. Finally, CANT1 expression was correlated with most immune checkpoints and drug sensitivity.

CONCLUSION

CANT1 may be a latent oncogene of HCC, and associated with immune cells and immune checkpoints, which may assist in HCC treatment.

Circ_0001715 accelerated lung adenocarcinoma process by the miR-1322/CANT1 axis

Lung adenocarcinoma (LUAD) is a type of lung cancer, which belongs to non-small cell lung cancer and has seriously endangered the physical and mental health of people. The study of circRNAs (circRNAs) has been increasingly hot in recent years, in which circRNAs also play an important regulatory role in cancer. The aim of this study was to investigate the biological molecular mechanisms of circ_0001715 in the progression of LUAD. The expression of circ_0001715, miR-1322 and calcium-activated nucleotidase 1 (CANT1) in LUAD tissues and cell lines was assessed by quantitative reverse transcription PCR (RT-qPCR) and western bot assay. Clone formation assay, 5-Ethynyl-2'-Deoxyuridine (EDU) assay and wound healing assay were used to verify the proliferation ability of cells. Dual-luciferase reporter assay and RNA pull-down assay were performed to characterize the interactions between the three factors. Finally, a mouse tumor model was constructed to assess the tumorigenicity of circ_0001715. RT-qPCR assay results showed that circ_0001715 expression was significantly increased in LUAD tissues and cell lines. Finally, knockdown of circ_0001715 could inhibit tumor growth in vivo. Circ_0001715 regulated the progression of LUAD through the miR-1322/CANT1 axis. The results of this study provided ideas for understanding the molecular mechanisms of circ_0001715 in LUAD.

Identification of DNA repair gene signature and potential molecular subtypes in hepatocellular carcinoma

DNA repair is a critical factor in tumor progression as it impacts tumor mutational burden, genome stability, PD-L1 expression, immunotherapy response, and tumor-infiltrating lymphocytes (TILs). In this study, we present a prognostic model for hepatocellular carcinoma (HCC) that utilizes genes related to the DNA damage response (DDR). Patients were stratified based on their risk score, and groups with lower risk scores demonstrated better survival rates compared to those with higher risk scores. The prognostic model's accuracy in predicting 1-, 3-, and 5-year survival rates for HCC patients was analyzed using receiver operator curve analysis (ROC). Results showed good accuracy in predicting survival rates. Additionally, we evaluated the prognostic model's potential as an independent factor for HCC prognosis, along with tumor stage. Furthermore, nomogram was employed to determine the overall survival year of patients with HCC based on this independent factor. Gene set enrichment analysis (GSEA) revealed that in the high-risk group, apoptosis, cell cycle, MAPK, mTOR, and WNT cascades were highly enriched. We used training and validation datasets to identify potential molecular subtypes of HCC based on the expression of DDR genes. The two subtypes differed in terms of checkpoint receptors for immunity and immune cell filtration capacity.Collectively, our study identified potential biomarkers of HCC prognosis, providing novel insights into the molecular mechanisms underlying HCC.

CANT1 serves as a potential prognostic factor for lung adenocarcinoma and promotes cell proliferation and invasion in vitro

Background

Calcium-activated nucleotidase 1 (CANT1), functions as a calcium-dependent nucleotidase with a preference for UDP. However, the potential clinical value of CANT1 in lung adenocarcinoma (LA) has not been fully clarified. Thus, we sought to identify its potential biological function and mechanism through bioinformatics analysis and in vitro experiments in LA.

Methods

In the present study, we comprehensively investigated the prognostic role of CANT1 in LA patients through bioinformatics analysis and in vitro experiments. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were utilized to analyze the expression of CANT1 in LA patients and their clinical-prognostic value. The immunohistochemistry staining was obtained from the Human Protein Atlas (HPA). A Cox regression model was used to evaluate prognostic factors. Gene ontology (GO) and Gene set enrichment analysis (GSEA) was performed to explore the potential regulatory mechanism of CANT1 in the development of LA. Moreover, we also examined the relationship between CANT1 expression and DNA methylation. Finally, we did in vitro experiments to evaluate the biological behavior and role of CANT1 in LA cells (LACs).

Results

Our study showed that the CANT1 expression was significantly elevated in the LA tissues compared with the normal lung tissues. Increased CANT1 expression was significantly associated with the TN stage. A univariate Cox analysis indicated that high CANT1 expression levels were correlated with poor overall survival (OS) in LA. Besides, CANT1 expression was independently associated with OS in multivariate analysis. GO and GSEA analysis showed the enrichment of mitotic nuclear division, DNA methylation, and DNA damage. Then we found that the high expression of CANT1 is positively correlated with hypomethylation. The methylation level was associated with prognosis in LA patients. Finally, in vitro experiments indicated that knockdown of CANT1 resulted in decreased cell proliferation, invasion, and G1 phase cell-cycle arrest in LACs.

Conclusion

The present study suggested that CANT1 may serve as a potential prognosis biomarker in patients with LA. High CANT1 expression and promoter demethylation was associated with worse outcome. Finally, in vitro experiments verified the biological functions and behaviors of CANT1 in LA.

Calcium-activated nucleotides 1 (CANT1)-driven nuclear factor-k-gene binding (NF-ĸB) signaling pathway facilitates the lung cancer progression

Dysregulation of calcium-activated nucleotides 1 (CANT1) has been observed in different organs. Thus, its biological function in cancer has increasingly attracted researchers. The current work aims to study the CANT1 role in lung cancer and understand the underlying pathological mechanisms. High amplification of CANT1 was observed in lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) tissues compared to normal tissues. The high-CANT1 patients showed a dismal prognosis in comparison with the low-CANT1 patients. Highly expressed CANT1 was significantly associated with the N stage of LUSC patients. Ectopic expression of CANT1 conspicuously increased the proliferation and viability of A549 cells. Conversely, CANT1 depletion resulted in adverse effects in H1299 cells. CANT1 depletion also resulted in the retardation of tumor growth in vivo. Mechanically, we found that CANT1 could elevate NF-ĸB (nuclear factor-k-gene binding) transcriptional activity in a concentration-dependent manner. This regulatory relationship was also established by the Western blot technique. Inhibiting NF-ĸB can significantly blunt the increased NF-κ-B Inhibitor-α (IκBα) expression caused by CANT1 overexpression in A549 cells. In conclusion, highly amplified CANT1 promotes the proliferation and viability of lung cancer cells. We also elucidate a new signaling axis of CANT1-NF-ĸB in lung cancer. This approach might be a promising strategy for lung cancer treatment.

The effect and mechanism of miR-607/CANT1 axis in lung squamous carcinoma

Lung squamous carcinoma (LUSC) is the second most frequent subtype of non-small cell lung cancer. Rarely gene alterations are identified in LUSC. Therefore, identifying LUSC-related genes to explain the relevant molecular mechanism is urgently needed. A potential biomarker, calcium-activated nucleotidase 1 (CANT1), was elevated in tissues of LUSC patients relative to normal cases based on the TCGA and/or GTEx database. CCK-8 and transwell tests were then implemented to measure the proliferative, invasive and migratory capacities, and showed that knockdown of CANT1 blocked LUSC cells proliferation. miR-607, predicted as an upstream factor for CANT1, was declined in LUSC using TargetScan analysis and luciferase activity test. Low miR-607 expression was related with unfavorable outcomes of LUSC patients. Moreover, miR-607 downregulation elevated cell viability, invasion and migration in LUSC cells, which was antagonized by si-CANT1. GEPIA website was accessed to estimate the relevance between CANT1 and epithelial-mesenchymal transition (EMT)-related positive factors. The protein levels of Fibronectin, Vimentin, Snail and β-catenin were altered due to the abnormal CANT1 and miR-607 expression. Together, these data unveiled that miR-607/CANT1 pair may exert a vital role in the progression of LUSC through mediating EMT process, which would furnish an available therapeutic therapy for LUSC.

A pharmacophore-based classification better predicts the outcomes of HER2-negative breast cancer patients receiving the anthracycline- and/or taxane-based neoadjuvant chemotherapy

AIMS

Prognosis of patients for human epidermal growth factor receptor 2 (HER2)-negative breast cancer post neoadjuvant chemotherapy is not well understood. The aim of this study was to develop a novel pharmacophore-based signature to better classify and predict the risk of HER2-negative patients after anthracycline-and/or taxane-based neoadjuvant chemotherapy (NACT).

MAIN METHODS

Anthracycline and taxane pharmacophore-based genes were obtained from PharmMapper. Drug-targeted genes (DTG) related clinical and bioinformatic analyses were undertaken in four GEO datasets.

KEY FINDINGS

We used 12 genes from the pharmacophore to develop a DTG score (DTG-S). The DTG-S classification exhibited significant prognostic ability with respect to disease free survival (DFS) for HER2-negative patients who receive at least one type of neoadjuvant chemotherapy that included anthracycline and/or taxane. DTG-S associated with a high predictive ability for pathological complete response (pCR) as well as for prognosis of breast cancer. Using the DTG-S classification in other prediction models may improve the reclassification accuracy for DFS. Combining the DTG-S with other clinicopathological factors may further improve its predictive ability of patients' outcomes. Gene ontology and KEGG pathway analysis showed that the biological processes of DTG-S high group were associated with the cell cycle, cell migration, and cell signal transduction pathways. Targeted drug analysis shows that some CDK inhibitors and PI3K-AKT pathway inhibitors may be useful for high DTG-S patients.

SIGNIFICANCE

The DTG-S classification adds prognostic and predictive information to classical parameters for HER2-negative patients who receive anthracycline-and/or taxane-based NACT, which could improve the patients' risk stratification and may help guide adjuvant treatment.

Identification and Characterization of Osmoregulation Related MicroRNAs in Gills of Hybrid Tilapia Under Three Types of Osmotic Stress

Researchers have increasingly suggested that microRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression and protein translation in organs and respond to abiotic and biotic stressors. To understand the function of miRNAs in osmotic stress regulation of the gills of hybrid tilapia (Oreochromis mossambicus ♀ × Oreochromis urolepis hornorum ♂), high-throughput Illumina deep sequencing technology was used to investigate the expression profiles of miRNAs under salinity stress (S, 25‰), alkalinity stress (A, 4‰) and salinity-alkalinity stress (SA, S: 15‰, A: 4‰) challenges. The results showed that 31, 41, and 27 upregulated and 33, 42, and 40 downregulated miRNAs (P < 0.05) were identified in the salt stress, alkali stress, and saline-alkali stress group, respectively, which were compared with those in the control group (C). Fourteen significantly differently expressed miRNAs were selected randomly and then validated by a quantitative polymerase chain reaction. On the basis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, genes related to osmoregulation and biosynthesis were enriched in the three types of osmotic stress. In addition, three miRNAs and three predicted target genes were chosen to conduct a quantitative polymerase chain reaction in the hybrid tilapia and its parents during 96-h osmotic stress. Differential expression patterns of miRNAs provided the basis for research data to further investigate the miRNA-modulating networks in osmoregulation of teleost.

Prognostic role of TSPAN1, KIAA1324 and ESRP1 in prostate cancer

The aim of this study was to validate prostate cancer‐associated genes on transcript level and to assess the prognostic value of the most promising markers by immunohistochemistry. Based on differentially expressed genes found in a previous study, 84 genes were further validated using mRNA expression data and follow‐up information from the Cancer Genome Atlas (TCGA) prostate cancer cohort (n = 497). Immunohistochemistry was used for validation of three genes in an independent, clinically annotated prostatectomy patient cohort (n = 175) with biochemical relapse as endpoint. Also, associations with clinicopathological variables were evaluated. Eleven protein‐coding genes from the list of 84 genes were associated with biochemical recurrence‐free survival on mRNA expression level in multivariate Cox‐analyses. Three of these genes (TSPAN1, ESRP1 and KIAA1324) were immunohistochemically validated using an independent cohort of prostatectomy patients. Both ESRP1 and KIAA1324 were independently associated with biochemical recurrence‐free survival. TSPAN1 was univariately prognostic but failed significance on multivariate analysis, probably due to its strong correlation with high Gleason scores. Multistep filtering using the publicly available TCGA cohort, data of an earlier expression profiling study which profiled 3023 cancer‐associated transcripts in 42 primary prostate cancer cases, identified two novel candidate prognostic markers (ESRP1 and KIAA1324) of primary prostate cancer for further study.

Human genes differ by their UV sensitivity estimated through analysis of UV-induced silent mutations in melanoma

We hypothesized that human genes differ by their sensitivity to ultraviolet (UV) exposure. We used somatic mutations detected by genome-wide screens in melanoma and reported in the Catalog Of Somatic Mutations In Cancer. As a measure of UV sensitivity, we used the number of silent mutations generated by C>T transitions in pyrimidine dimers of a given transcript divided by the number of potential sites for this type of mutations in the transcript. We found that human genes varied by UV sensitivity by two orders of magnitude. We noted that the melanoma-associated tumor suppressor gene CDKN2A was among the top five most UV-sensitive genes in the human genome. Melanoma driver genes have a higher UV-sensitivity compared with other genes in the human genome. The difference was more prominent for tumor suppressors compared with oncogene. The results of this study suggest that differential sensitivity of human transcripts to UV light may explain melanoma specificity of some driver genes. Practical significance of the study relates to the fact that differences in UV sensitivity among human genes need to be taken into consideration whereas predicting melanoma-associated genes by the number of somatic mutations detected in a given gene.

Prognosis and personalized treatment prediction in TP53-mutant hepatocellular carcinoma: an in silico strategy towards precision oncology

TP53 mutation is one of the most common genetic changes in hepatocellular carcinoma (HCC). It is of great clinical significance to tailor specialized prognostication approach and to explore more therapeutic options for TP53-mutant HCCs. In this study, a total of 1135 HCC patients were retrospectively analyzed. We developed a random forest-based prediction model to estimate TP53 mutational status, tackling the problem of limited sample size in TP53-mutant HCCs. A multi-step process was performed to develop robust poor prognosis-associated signature (PPS). Compared with previous established population-based signatures, PPS manifested superior ability to predict survival in TP53-mutant patients. After in silico screening of 2249 drug targets and 1770 compounds, we found that three targets (CANT1, CBFB and PKM) and two agents (irinotecan and YM-155) might have potential therapeutic implications in high-PPS patients. The results of drug targets prediction and compounds prediction complemented each other, presenting a comprehensive view of potential treatment strategy. Overall, our study has not only provided new insights into personalized prognostication approaches, but also thrown light on integrating tailored risk stratification with precision therapy.

An atlas of transposable element-derived alternative splicing in cancer

Transposable element (TE)-derived sequences comprise more than half of the human genome, and their presence has been documented to alter gene expression in a number of different ways, including the generation of alternatively spliced transcript isoforms. Alternative splicing has been associated with tumorigenesis for a number of different cancers. The objective of this study was to broadly characterize the role of human TEs in generating alternatively spliced transcript isoforms in cancer. To do so, we screened for the presence of TE-derived sequences co-located with alternative splice sites that are differentially used in normal versus cancer tissues. We analysed a comprehensive set of alternative splice variants characterized for 614 matched normal-tumour tissue pairs across 13 cancer types, resulting in the discovery of 4820 TE-generated alternative splice events distributed among 723 cancer-associated genes. Short interspersed nuclear elements (Alu) and long interspersed nuclear elements (L1) were found to contribute the majority of TE-generated alternative splice sites in cancer genes. A number of cancer-associated genes, including MYH11, WHSC1 and CANT1, were shown to have overexpressed TE-derived isoforms across a range of cancer types. TE-derived isoforms were also linked to cancer-specific fusion transcripts, suggesting a novel mechanism for the generation of transcriptome diversity via trans-splicing mediated by dispersed TE repeats. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Classification of early and late stage liver hepatocellular carcinoma patients from their genomics and epigenomics profiles

Background

Liver Hepatocellular Carcinoma (LIHC) is one of the major cancers worldwide, responsible for millions of premature deaths every year. Prediction of clinical staging is vital to implement optimal therapeutic strategy and prognostic prediction in cancer patients. However, to date, no method has been developed for predicting the stage of LIHC from the genomic profile of samples.

Methods

The Cancer Genome Atlas (TCGA) dataset of 173 early stage (stage-I), 177 late stage (stage-II, Stage-III and stage-IV) and 50 adjacent normal tissue samples for 60,483 RNA transcripts and 485,577 methylation CpG sites, was extensively analyzed to identify the key transcriptomic expression and methylation-based features using different feature selection techniques. Further, different classification models were developed based on selected key features to categorize different classes of samples implementing different machine learning algorithms.

Results

In the current study, in silico models have been developed for classifying LIHC patients in the early vs. late stage and cancerous vs. normal samples using RNA expression and DNA methylation data. TCGA datasets were extensively analyzed to identify differentially expressed RNA transcripts and methylated CpG sites that can discriminate early vs. late stages and cancer vs. normal samples of LIHC with high precision. Naive Bayes model developed using 51 features that combine 21 CpG methylation sites and 30 RNA transcripts achieved maximum MCC (Matthew’s correlation coefficient) 0.58 with an accuracy of 78.87% on the validation dataset in discrimination of early and late stage. Additionally, the prediction models developed based on 5 RNA transcripts and 5 CpG sites classify LIHC and normal samples with an accuracy of 96–98% and AUC (Area Under the Receiver Operating Characteristic curve) 0.99. Besides, multiclass models also developed for classifying samples in the normal, early and late stage of cancer and achieved an accuracy of 76.54% and AUC of 0.86.

Conclusion

Our study reveals stage prediction of LIHC samples with high accuracy based on the genomics and epigenomics profiling is a challenging task in comparison to the classification of cancerous and normal samples. Comprehensive analysis, differentially expressed RNA transcripts, methylated CpG sites in LIHC samples and prediction models are available from CancerLSP (http://webs.iiitd.edu.in/raghava/cancerlsp/).

Associations of Genetic Variations in MicroRNA Seed Regions With Acute Adverse Events and Survival in Patients With Rectal Cancer Receiving Postoperative Chemoradiation Therapy

PURPOSE

The aim of this study was to investigate the associations between single nucleotide polymorphisms (SNPs) in the seed regions of microRNAs and acute adverse events (AEs) and survival in patients with rectal cancer receiving postoperative chemoradiation therapy.

METHODS AND MATERIALS

Eighteen SNPs were genotyped in 365 patients with rectal cancer receiving postoperative chemoradiation therapy. The associations between genotypes and AEs were estimated by odds ratios and 95% confidence intervals (CIs), which were computed by using multivariate logistic regression models. The hazard ratios and 95% CIs to assess the death of patients for different genotypes were calculated by Cox proportional regression models. Overall survival and disease-free survival of patients with different genotypes were estimated by Kaplan-Meier plots, and the statistical significance was determined by using the log-rank test.

RESULTS

In these patients, the most common grade ≥2 AEs were diarrhea (44.1%), leukopenia (29.6%), and dermatitis (18.9%). With false discovery rate correction, SNP rs2273626 was significantly associated with a decreased risk of grade ≥2 leukopenia (odds ratio, 0.48; 95% CI, 0.31-0.74; P = .0009). In addition, SNP rs202195689 was associated with overall survival and disease-free survival in patients receiving postoperative chemoradiation therapy, with the hazard ratios for death being 2.02 (95% CI, 1.36-3.01; P = .0006) and 1.91 (95% CI, 1.36-2.70; P = .0002), respectively. However, no significant association between these SNPs and diarrhea and dermatitis was observed.

CONCLUSIONS

These results suggest that rs2273626 and rs202195689 in microRNA seed regions might serve as independent biomarkers for predicting AEs and prognosis in patients with rectal cancer receiving postoperative chemoradiation therapy. Independent replication of these findings is required to confirm these results.

The role of tumor metabolism as a driver of prostate cancer progression and lethal disease: results from a nested case-control study

Background

Understanding the biologic mechanisms underlying the development of lethal prostate cancer is critical for improved therapeutic and prevention strategies. In this study we explored the role of tumor metabolism in prostate cancer progression using mRNA expression profiling of seven metabolic pathways; fatty acid metabolism, glycolysis/gluconeogenesis, oxidative phosphorylation, pentose phosphate, purine metabolism, pyrimidine metabolism and the tricarboxylic acid cycle.

Methods

The study included 404 men with archival formalin-fixed, paraffin-embedded prostate tumor tissue from the prospective Health Professionals Follow-up Study and Physicians’ Health Study. Lethal cases (n = 113) were men who experienced a distant metastatic event or died of prostate cancer during follow-up. Non-lethal controls (n = 291) survived at least 8 years post-diagnosis without metastases. Of 404 men, 202 additionally had matched normal tissue (140 non-lethal, 62 lethal). Analyses compared expression levels between tumor and normal tissue, by Gleason grade and by lethal status. Secondary analyses considered the association with biomarkers of cell proliferation, apoptosis and angiogenesis.

Results

Oxidative phosphorylation and pyrimidine metabolism were identified as the most dysregulated pathways in lethal tumors (p < 0.007), and within these pathways, a number of novel differentially expressed genes were identified including POLR2K and APT6V1A. The associations were tumor specific as there was no evidence any pathways were altered in the normal tissue of lethal compared to non-lethal cases.

Conclusions

The results suggest prostate cancer progression and lethal disease are associated with alterations in key metabolic signaling pathways. Pathways supporting proliferation appeared to be of particular importance in prostate tumor aggressiveness.

Electronic supplementary material

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

Functional annotation of HOT regions in the human genome: implications for human disease and cancer

Advances in genome-wide association studies (GWAS) and large-scale sequencing studies have resulted in an impressive and growing list of disease- and trait-associated genetic variants. Most studies have emphasised the discovery of genetic variation in coding sequences, however, the noncoding regulatory effects responsible for human disease and cancer biology have been substantially understudied. To better characterise the cis-regulatory effects of noncoding variation, we performed a comprehensive analysis of the genetic variants in HOT (high-occupancy target) regions, which are considered to be one of the most intriguing findings of recent large-scale sequencing studies. We observed that GWAS variants that map to HOT regions undergo a substantial net decrease and illustrate development-specific localisation during haematopoiesis. Additionally, genetic risk variants are disproportionally enriched in HOT regions compared with LOT (low-occupancy target) regions in both disease-relevant and cancer cells. Importantly, this enrichment is biased toward disease- or cancer-specific cell types. Furthermore, we observed that cancer cells generally acquire cancer-specific HOT regions at oncogenes through diverse mechanisms of cancer pathogenesis. Collectively, our findings demonstrate the key roles of HOT regions in human disease and cancer and represent a critical step toward further understanding disease biology, diagnosis, and therapy.

Comprehensive validation of published immunohistochemical prognostic biomarkers of prostate cancer -what has gone wrong? A blueprint for the way forward in biomarker studies

Background:

Treatment planning of localised prostate cancer remains challenging. Besides conventional parameters, a wealth of prognostic biomarkers has been proposed so far. None of which, however, have successfully been implemented in a routine setting so far. The aim of our study was to systematically verify a set of published prognostic markers for prostate cancer.

Methods:

Following an in-depth PubMed search, 28 markers were selected that have been proposed as multivariate prognostic markers for primary prostate cancer. Their prognostic validity was examined in a radical prostatectomy cohort of 238 patients with a median follow-up of 60 months and biochemical progression as endpoint of the analysis. Immunohistochemical evaluation was performed using previously published cut-off values, but allowing for optimisation if necessary. Univariate and multivariate Cox regression were used to determine the prognostic value of biomarkers included in this study.

Results:

Despite the application of various cut-offs in the analysis, only four (14%) markers were verified as independently prognostic (AKT1, stromal AR, EZH2, and PSMA) for PSA relapse following radical prostatectomy.

Conclusions:

Apparently, many immunohistochemistry-based studies on prognostic markers seem to be over-optimistic. Codes of best practice, such as the REMARK guidelines, may facilitate the performance of conclusive and transparent future studies.

Restoration of compact Golgi morphology in advanced prostate cancer enhances susceptibility to galectin-1-induced apoptosis by modifying mucin O-glycan synthesis

Prostate cancer progression is associated with upregulation of sialyl-T antigen produced by β-galactoside α-2,3-sialyltransferase-1 (ST3Gal1) but not with core 2-associated polylactosamine despite expression of core 2 N-acetylglucosaminyltransferase-L (C2GnT-L/GCNT1). This property allows androgen-refractory prostate cancer cells to evade galectin-1 (LGALS1)-induced apoptosis, but the mechanism is not known. We have recently reported that Golgi targeting of glycosyltransferases is mediated by golgins: giantin (GOLGB1) for C2GnT-M (GCNT3) and GM130 (GOLGA2)-GRASP65 (GORASP1) or GM130-giantin for core 1 synthase. Here, we show that for Golgi targeting, C2GnT-L also uses giantin exclusively whereas ST3Gal1 uses either giantin or GM130-GRASP65. In addition, the compact Golgi morphology is detected in both androgen-sensitive prostate cancer and normal prostate cells, but fragmented Golgi and mislocalization of C2GnT-L are found in androgen-refractory cells as well as primary prostate tumors (Gleason grade 2-4). Furthermore, failure of giantin monomers to be phosphorylated and dimerized prevents Golgi from forming compact morphology and C2GnT-L from targeting the Golgi. On the other hand, ST3Gal1 reaches the Golgi by an alternate site, GM130-GRASP65. Interestingly, inhibition or knockdown of non-muscle myosin IIA (MYH9) motor protein frees up Rab6a GTPase to promote phosphorylation of giantin by polo-like kinase 3 (PLK3), which is followed by dimerization of giantin assisted by protein disulfide isomerase A3 (PDIA3), and restoration of compact Golgi morphology and targeting of C2GnT-L. Finally, the Golgi relocation of C2GnT-L in androgen-refractory cells results in their increased susceptibility to galectin-1-induced apoptosis by replacing sialyl-T antigen with polylactosamine.

IMPLICATIONS

This study demonstrates the importance of Golgi morphology and regulation of glycosylation and provides insight into how the Golgi influences cancer progression and metastasis.

Objective assessment of cancer genes for drug discovery

Selecting the best targets is a key challenge for drug discovery, and achieving this effectively, efficiently and systematically is particularly important for prioritizing candidates from the sizeable lists of potential therapeutic targets that are now emerging from large-scale multi-omics initiatives, such as those in oncology. Here, we describe an objective, systematic, multifaceted computational assessment of biological and chemical space that can be applied to any human gene set to prioritize targets for therapeutic exploration. We use this approach to evaluate an exemplar set of 479 cancer-associated genes, reveal the tension between biological relevance and chemical tractability, and describe major gaps in available knowledge that could be addressed to aid objective decision-making. We also propose drug repurposing opportunities and identify potentially druggable cancer-associated proteins that have been poorly explored with regard to the discovery of small-molecule modulators, despite their biological relevance.

Comparison of p40 (ΔNp63) and p63 expression in prostate tissues--which one is the superior diagnostic marker for basal cells?

AIMS

p63 is one of the standard markers for basal cells of the prostate gland. Recently, it has been suggested that the p63 isoform p40 might be more specific as a basal cell marker. In this study we compare the staining characteristics of p63 and p40 in normal and malignant prostate tissues.

METHODS AND RESULTS

A prostatectomy cohort (n = 640) in tissue microarray format was evaluated for p63 (clone 4A4) and for p40 (rabbit polyclonal) immunoreactivity in malignant and normal tissues. Immunoreactivity of basal and secretory cells was evaluated in a semiquantitative manner and compared case-wise. In normal tissues, p40 showed highly similar immunoreactivity compared to p63. The staining patterns were absolutely identical in 88% of cases. Additional cytoplasmic p40 staining in tumour cells occurred in 59.6% of cancer cases. Differences were seen in nuclear staining of carcinomas: 1.4% of carcinomas were p63-positive, whereas 0.6% were p40-positive.

CONCLUSIONS

In most cases, p40 stains prostatic basal cells as reliably as p63, with only minor differences. Aberrant staining of tumour cells is seen more rarely than with p63 (clone 4A4), establishing its higher specificity. However, additional cytoplasmic immunoreactivity narrows its eligibility for antibody cocktails (e.g. with alpha-methylacyl-CoA racemase).

Revelation of a catalytic calcium-binding site elucidates unusual metal dependence of a human apyrase

Human soluble calcium-activated nucleotidase 1 (hSCAN-1) represents a new family of apyrase enzymes that catalyze the hydrolysis of nucleotide di- and triphosphates, thereby modulating extracellular purinergic and pyrimidinergic signaling. Among well-characterized phosphoryl transfer enzymes, hSCAN-1 is unique not only in its unusual calcium-dependent activation, but also in its novel phosphate-binding motif. Its catalytic site does not utilize backbone amide groups to bind phosphate, as in the common P-loop, but contains a large cluster of acidic ionizable side chains. By employing a state-of-the-art computational approach, we have revealed a previously uncharacterized catalytic calcium-binding site in hSCAN-1, which elucidates the unusual calcium-dependence of its apyrase activity. In a high-order coordination shell, the newly identified calcium ion organizes the active site residues to mediate nucleotide binding, to orient the nucleophilic water, and to facilitate the phosphoryl transfer reaction. From ab initio QM/MM molecular dynamics simulations with umbrella sampling, we have characterized a reverse protonation catalytic mechanism for hSCAN-1 and determined its free energy reaction profile. Our results are consistent with available experimental studies and provide new detailed insight into the structure-function relationship of this novel calcium-activated phosphoryl transfer enzyme.

Implications of the Golgi apparatus in prostate cancer

The classical view of the Golgi apparatus is of a small membranous organelle involved in protein transport and secretion. Recent descriptions of the molecular network connecting the Golgi to other organelles demonstrate the essential roles of the Golgi in cellular activities as a stress sensor, apoptosis trigger, lipid/protein modifier, mitotic checkpoint, and a mediator of malignant transformation. Thus, the Golgi function should have a fundamental impact on cancer cell survival. Prostate cancer is initially responsive to androgenic hormones; however, it almost invariably progresses to a castration-refractory or hormone-insensitive state. Nevertheless, androgen signaling remains active at this stage and is important as a therapeutic target. Certain Golgi-associated molecules have recently been demonstrated to be regulated by androgen action, and the Golgi is emerging as a new therapeutic target in prostate cancer. The key Golgi-associated molecules essential for prostate cancer development and the potential therapeutic options targeting the Golgi apparatus are discussed.