DNA primase polypeptide 1 (PRIM1) involves in estrogen-induced breast cancer formation through activation of the G2/M cell cycle checkpoint

RNA Sequencing and Bioinformatics Analysis Reveals the Downregulation of DNA Replication Genes by Morindone in Colorectal Cancer Cells

Morindone, a natural anthraquinone compound, has been reported to have significant pharmacological properties in different cancers. However, its anticancer effects in colorectal cancer (CRC) and the underlying molecular mechanisms remain obscure. In this study, RNA sequencing was used to assess the differentially expressed genes (DEGs) following morindone treatment in two CRC cell lines, HCT116 and HT29 cells. Functional enrichment analysis of overlapping DEGs revealed that negative regulation of cell development from biological processes and the MAPK signalling pathway were the most significant Gene Ontology terms and Kyoto Encyclopaedia of Genes and Genome pathway, respectively. Seven hub genes were identified among the overlapping genes, including MCM5, MCM6, MCM10, GINS2, POLE2, PRIM1, and WDHD1. All hub genes were found downregulated and involved in DNA replication fork. Among these, GINS2 was identified as the most cancer-dependent gene in both cells with better survival outcomes. Validation was performed on seven hub genes with rt-qPCR, and the results were consistent with the RNA sequencing findings. Collectively, this study provides corroboration of the potential therapeutic benefits and suitable pharmacological targets of morindone in the treatment of CRC.

The selenoenzyme type I iodothyronine deiodinase: a new tumor suppressor in ovarian cancer

The selenoenzyme type I iodothyronine deiodinase (DIO1) catalyzes removal of iodine atoms from thyroid hormones. Although DIO1 action is reported to be disturbed in several malignancies, no work has been conducted in high-grade serous ovarian carcinoma (HGSOC), the most lethal gynecologic cancer. We studied DIO1 expression in HGSOC patients [The Cancer Genome Atlas (TCGA) data and tumor tissues], human cell lines (ES-2 and Kuramochi), normal Chinese hamster ovarian cells (CHO-K1), and normal human fallopian tube cells (FT282 and FT109). To study its functional role, DIO1 was overexpressed, inhibited [by propylthiouracil (PTU)], or knocked down (KD), and cell count, proliferation, apoptosis, cell viability, and proteomics analysis were performed. Lower DIO1 levels were observed in HGSOC compared to normal cells and tissues. TCGA analyses confirmed that low DIO1 mRNA expression correlated with worse survival and therapy resistance in patients. Silencing or inhibiting the enzyme led to enhanced ovarian cancer proliferation, while an opposite effect was shown following DIO1 ectopic expression. Proteomics analysis in DIO1-KD cells revealed global changes in proteins that facilitate tumor metabolism and progression. In conclusion, DIO1 expression and ovarian cancer progression are inversely correlated, highlighting a tumor suppressive role for this enzyme and its potential use as a biomarker in this disease.

YAP/TEAD-induced PRIM1 contributes to the progression and poor prognosis of gastric carcinoma

Gastric carcinoma has a poor prognosis and low survival rate. PRIM1 is closely associated with the origin of DNA replication and serves as a carcinogenic factor in multiple tumors. This study aimed to explore the functions of PRIM1 in the progression of gastric carcinoma. The luciferase reporter assay examined the regulatory effect of YAP1/TEAD4 on PRIM1. A xenograft tumor mouse model was constructed to observe cancer cell proliferation in vivo. The upregulation of PRIM1 was found in gastric carcinoma cells and tissues, and it was associated with poor prognosis. Silencing PRIM1 inhibited cell proliferation, arrested the cell cycle, and upregulated Cdc25, Cyclin B, and Cdc2 expression. In addition, apoptosis was increased upon PRIM1 knockdown, accompanied by increased protein levels of cleaved caspase-3 and caspase-8. In vivo, knockdown of PRIM1 suppressed the growth of xenograft tumors formed by gastric carcinoma cells. Moreover, PRIM1 silencing elevated the chemosensitivity of gastric carcinoma cells. By investigating molecular events downstream of the Hippo signaling pathway, we found that PRIM1 was a target gene of the YAP1/TEAD4 transcriptional regulatory complex. PRIM1 represents a novel target for gastric carcinoma therapeutic approaches.

Genome-Wide Association Study of Breast Density among Women of African Ancestry

Breast density, the amount of fibroglandular versus fatty tissue in the breast, is a strong breast cancer risk factor. Understanding genetic factors associated with breast density may help in clarifying mechanisms by which breast density increases cancer risk. To date, 50 genetic loci have been associated with breast density, however, these studies were performed among predominantly European ancestry populations. We utilized a cohort of women aged 40-85 years who underwent screening mammography and had genetic information available from the Penn Medicine BioBank to conduct a Genome-Wide Association Study (GWAS) of breast density among 1323 women of African ancestry. For each mammogram, the publicly available "LIBRA" software was used to quantify dense area and area percent density. We identified 34 significant loci associated with dense area and area percent density, with the strongest signals in GACAT3, CTNNA3, HSD17B6, UGDH, TAAR8, ARHGAP10, BOD1L2, and NR3C2. There was significant overlap between previously identified breast cancer SNPs and SNPs identified as associated with breast density. Our results highlight the importance of breast density GWAS among diverse populations, including African ancestry populations. They may provide novel insights into genetic factors associated with breast density and help in elucidating mechanisms by which density increases breast cancer risk.

Cell-cycle and apoptosis related and proteomics-based signaling pathways of human hepatoma Huh-7 cells treated by three currently used multi-RTK inhibitors

Sorafenib, lenvatinib and regorafenib, the multi-RTK inhibitors with potent anti-angiogenesis effects, are currently therapeutic drugs generally recommended for the patients with advanced hepatocellular carcinoma (HCC). To date, however, there have been no published studies on the mechanism underling differential effects of the three drugs on HCC cell proliferation, and the proteomic analysis in HCC cell lines treated by regorafenib or lenvatinib. The present study for the first time performed a direct comparison of the cell cycle arrest and apoptosis induction in the Huh-7 cells caused by sorafenib, regorafenib and lenvatinib at respective IC₅₀ using flow cytometry technique, as well as their pharmacological interventions for influencing whole cell proteomics using tandem mass tag-based peptide-labeling coupled with the nLC-HRMS technique. Sorafenib, regorafenib and lenvatinib at respective IC₅₀ drove the remaining surviving Huh-7 cells into a G₀/G₁ arrest, but lenvatinib and regorafenib were much more effective than sorafenib. Lenvatinib produced a much stronger induction of Huh-7 cells into early apoptosis than sorafenib and regorafenib, while necrotic cell proportion induced by regorafenib was 2.4 times as large as that by lenvatinib. The proteomic study revealed 419 proteins downregulated commonly by the three drugs at respective IC₅₀. KEGG pathway analysis of the downregulated proteins indicated the ranking of top six signaling pathways including the spliceosome, DNA replication, cell cycle, mRNA surveillance, P53 and nucleotide excision repair involved in 33 proteins, all of which were directly related to their pharmacological effects on cell cycle and cell apoptosis. Notably, lenvatinib and regorafenib downregulated the proteins of PCNA, Cyclin B1, BCL-xL, TSP1, BUD31, SF3A1 and Mad2 much more strongly than sorafenib. Moreover, most of the proteins in the P53 signaling pathway were downregulated with lenvatinib and regorafenib by more than 36% at least. In conclusion, lenvatinib and regorafenib have much stronger potency against Huh-7 cell proliferation than sorafenib because of their more potent effects on cell cycle arrest and apoptosis induction. The underling mechanism may be at least due to the 33 downregulated proteins centralizing the signal pathways of cell cycle, p53 and DNA synthesis based on the present proteomics study.

Identification of target and pathway of aspirin combined with Lipitor treatment in prostate cancer through integrated bioinformatics analysis

PURPOSE

Our previous studies have confirmed that aspirin combined with Lipitor inhibited the development of prostate cancer (PCa), but the mechanisms need to be comprehensively expounded. The study aims to screen out the hub genes of combination therapy and to explore their association with the pathogenesis and prognosis of PCa.

METHODS

Gene expressions were quantified by RNA sequencing (RNA-seq). Altered biological function, pathways of differentially expressed genes (DEGs), protein-protein interaction network, the filtering of hub genes, gene co-expression and the pathogenesis and prognosis were revealed by bioinformatics analysis. The correlation between hub gene expression and patient survival was validated by Kaplan-Meier. The effects of silent DNA replication and sister chromatid cohesion 1 (siDSCC1) combined with Lipitor and aspirin on DSCC1 expression, viability, invasion and migration of PCa cells were detected by qRT-PCR, Wound healing and transwell assays.

RESULTS

157 overlapped DEGs involved in FoxO, PI3K-Akt and p53 signaling pathways were identified. Ten hub genes (NEIL3, CDC7, DSCC1, CDC25C, PRIM1, MCM10, FBXO5, DTL, SERPINE1, EXO1) were verified to be correlated with the pathology and prognosis of PCa. DSCC1 silencing not only inhibited the viability, migration and invasion of PCa cells, but also strengthened the suppressing effects of Lipitor and aspirin alone or in combination on PCa cells.

CONCLUSION

The enrichment pathways and targets of Lipitor combined with aspirin in PCa are discovered, and DSCC1 silencing can potentiate the effect of Lipitor combined with aspirin in the treatment of PCa.

Statistical Bioinformatics to Uncover the Underlying Biological Mechanisms That Linked Smoking with Type 2 Diabetes Patients Using Transcritpomic and GWAS Analysis

Type 2 diabetes (T2D) is a chronic metabolic disease defined by insulin insensitivity corresponding to impaired insulin sensitivity, decreased insulin production, and eventually failure of beta cells in the pancreas. There is a 30–40 percent higher risk of developing T2D in active smokers. Moreover, T2D patients with active smoking may gradually develop many complications. However, there is still no significant research conducted to solve the issue. Hence, we have proposed a highthroughput network-based quantitative pipeline employing statistical methods. Transcriptomic and GWAS data were analysed and obtained from type 2 diabetes patients and active smokers. Differentially Expressed Genes (DEGs) resulted by comparing T2D patients’ and smokers’ tissue samples to those of healthy controls of gene expression transcriptomic datasets. We have found 55 dysregulated genes shared in people with type 2 diabetes and those who smoked, 27 of which were upregulated and 28 of which were downregulated. These identified DEGs were functionally annotated to reveal the involvement of cell-associated molecular pathways and GO terms. Moreover, protein–protein interaction analysis was conducted to discover hub proteins in the pathways. We have also identified transcriptional and post-transcriptional regulators associated with T2D and smoking. Moreover, we have analysed GWAS data and found 57 common biomarker genes between T2D and smokers. Then, Transcriptomic and GWAS analyses are compared for more robust outcomes and identified 1 significant common gene, 19 shared significant pathways and 12 shared significant GOs. Finally, we have discovered protein–drug interactions for our identified biomarkers.

Analysis of core genes for colorectal cancer prognosis based on immune and stromal scores

Background

Colorectal cancer (CRC) is one of the most common malignancies.An early diagnosis and an accurate prognosis are major focuses of CRC research. Tumor microenvironment cells and the extent of infiltrating immune and stromal cells contribute significantly to the tumor prognosis.

Methods

Immune and stromal scores were calculated based on the ESTIMATE algorithm using the sample expression profile of the The Cancer Genome Atlas (TCGA) database. GSE102479 was used as the validation database. Differentially expressed genes whose expression was significantly associated with the prognosis of CRC patients were identified based on the immune matrix score. Survival analysis was conducted on the union of the differentially expressed genes. A protein–protein interaction (PPI) network was constructed using the STRING database to identify the closely connected modules. To conduct functional enrichment analysis of the relevant genes, GO and KEGG pathway analyses were performed with Cluster Profiler. Pivot analysis of the ncRNAs and TFs was performed by using the RAID2.0 database and TRRUST v2 database. TF-mRNA regulatory relationships were analyzed in the TRRUST V2 database. Hubgene targeting relationships were screened in the TargetScan, miRTarBase and miRDB databases. The SNV data of the hub genes were analyzed by using the R maftools package. A ROC curve was drawn based on the TCGA database. The proportion of immune cells was estimated using CIBERSORT and the LM22 feature matrix.

Results

The results showed that the matrix score was significantly correlated with colorectal cancer stage T. A total of 789 differentially expressed genes and 121 survival-related prognostic genes were identified. The PPI network showed that 22 core genes were related to the CRC prognosis. Furthermore, four ncRNAs that regulated the core prognosis genes, 11 TFs with regulatory effects on the core prognosis genes, and two drugs, quercetin and pseudoephedrine, that have regulatory effects on colorectal cancer were also identified.

Conclusions

We obtained a list of tumor microenvironment-related genes for CRC patients. These genes could be useful for determining the prognosis of CRC patients. To confirm the function of these genes, additional experiments are necessary.

Bioinformatics Analysis of Potential Biomarkers and Pathway Identification for Major Depressive Disorder

Aiming at a more comprehensive understanding of the molecular biomarkers and potential mechanisms of major depressive disorder (MDD), from the Gene Expression Omnibus (GEO) database, we first obtained mRNA expression profiles and identified 585 differentially expressed genes (DEGs) through the R software, including 263 upregulated genes and 322 downregulated genes. Then, through the Kyoto Encyclopedia of Genome and Genome (KEGG) pathway and biological process (BP) analysis, we found that the upregulated and downregulated DEGs were abundant in different pathways, respectively. It was noteworthy that upregulated DEGs were the most significantly enriched in the mTOR signaling pathway. Subsequently, through the protein-protein interaction (PPI) network, we identified seven hub genes, namely, EXOSC2, CAMK2A, PRIM1, SMC4, TYMS, CDK6, and RPA2. Finally, through gene set enrichment analysis (GSEA), we obtained that hypoxia, epithelial-mesenchymal transition, hedgehog signaling, and reactive oxygen species pathway were the enriched pathways for MDD patients. The above data results would provide a new direction for the treatment of MDD patients.

Screening of genes related to survival prognosis of cervical squamous cell carcinoma and construction of prognosis prediction model

AIM

We aimed to screen for the genes related to survival prognosis of cervical squamous cell carcinoma (CSCC) and then constructed a prognosis prediction model.

METHODS

The GSE63514 dataset was obtained from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO). The CSCC gene dataset and the GSE44001 dataset were obtained from The Cancer Genome Atlas and NCBI GEO, respectively. The Kaplan-Meier (KM) curve was used to evaluate the association between high and low prognosis that was with the actual survival prognosis information. The Cox proportional hazards model was used to screen out the optimized prognostic-related signature differentially expressed gene (DEG) combinations. Gene set enrichment analysis was used to perform pathway enrichment annotation analysis for DEGs that were related to risk grouping.

RESULTS

In total, 16 399 DEGs were obtained and 23 gene ontology biological processes and 8 Kyoto Encyclopedia of Genes and Genomes pathways were screened. Nine optimized DEG groups related to independent prognosis were selected. The KM curves of pathologic N0 and N1 showed that low-risk group were associated with a better overall survival (p = 1.518e; p = 1.704e-01). The pathways related to risk grouping were cytokine-cytokine receptor interaction, JAK stat signaling pathway, and glycolysis-gluconeogenesis.

CONCLUSION

On the basis of this study, we established a prognostic risk model, which provided a reliable prognostic tool and was of great significance for locating the biomarkers related to survival prognosis in CSCC.

Six Metabolism Related mRNAs Predict the Prognosis of Patients With Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is one of the most common aggressive solid malignant tumors and current research regards HCC as a type of metabolic disease. This study aims to establish a metabolism-related mRNA signature model for risk assessment and prognosis prediction in HCC patients. Methods: HCC data were obtained from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC) and Gene Enrichment Analysis (GSEA) website. Least absolute shrinkage and selection operator (LASSO) was used to screen out the candidate mRNAs and calculate the risk coefficient to establish the prognosis model. A high-risk group and low-risk group were separated for further study depending on their median risk score. The reliability of the prediction was evaluated in the validation cohort and the whole cohort. Results: A total of 548 differential mRNAs were identified from HCC samples (n = 374) and normal controls (n = 50), 45 of which were correlated with prognosis. A total of 373 samples met the screening criteria and there were randomly divided into the training cohort (n = 186) and the validation cohort (n = 187). In the training cohort, six metabolism-related mRNAs were used to construct a prognostic model with a LASSO regression model. Based on the risk model, the overall survival rate of the high-risk cohort was significantly lower than that of the low-risk cohort. The results of a time-ROC curve proved that the risk score (AUC = 0.849) had a higher prognostic value than the pathological grade, clinical stage, age or gender. Conclusion: The model constructed by the six metabolism-related mRNAs has a significant value for survival prediction and can be applied to guide the evaluation of HCC and the designation of clinical therapy.

DNA primase subunit 1 deteriorated progression of hepatocellular carcinoma by activating AKT/mTOR signaling and UBE2C-mediated P53 ubiquitination

Background

DNA primase subunit 1 (PRIM1) has been reported as a novel oncogene in several cancer types. However, its roles in hepatocellular carcinoma (HCC) remain unclear. This study aimed to investigate underlying mechanisms of PRIM1 and identify it as a potential molecular target for HCC.

Methods

Hub genes were screened between HCC tissues and normal liver tissues in 3 gene expression omnibus (GEO) datasets and the cancer genome atlas (TCGA). The expression features and prognostic value of one of the hub genes PRIM1 were analyzed by bioinformatic analyses and immunohistochemistry. Loss-of-function and gain-of-function studies were used to investigate the regulatory role of PRIM1 in HCC cells. Real-time (RT)-qPCR, western blotting, and ubiquitin immunoprecipitation assays were performed to explore the underlying mechanisms. The xenograft model was employed to detect the roles of PRIM1 in tumor growth in vivo. Finally, the 3D spheroid model was conducted to validate the role of PRIM1 in tumor growth and sorafenib resistance.

Results

The hub genes of HCC were screened in multiple bioinformatic datasets. PRIM1, as one of the hub genes, was significantly overexpressed in HCC tissues in mRNA and protein levels. In addition, high expression of PRIM1 indicated poor prognosis of HCC patients in TCGA, ICGC, and Nantong cohorts. Overexpression of PRIM1 promoted the proliferation, migration/invasion, and sorafenib resistance of HCC cells, with the decrease in apoptosis and cell cycle arrest. Mechanically, PRIM1 facilitated epithelial-mesenchymal transition (EMT) process and the activity of PI3K/AKT/mTOR signaling of HCC cells. Additionally, PRIM1 could cause the ubiquitination and degradation of P53 by upregulating Ubiquitin Conjugating Enzyme E2 C (UBE2C). Furthermore, knockdown of PRIM1 significantly inhibited the growth of xenograft tumors and HCC cells-derived spheroids with enhanced sorafenib resistance.

Conclusion

This study implies that PRIM1 may play a key role in the progression of HCC and may serve as a potential target for HCC treatment.

Identification of estrogen receptor target genes involved in gonadal feminization caused by estrogen in Xenopus laevis

Estrogens and estrogenic endocrine disrupting chemicals can cause gonadal feminization in some vertebrates mainly through estrogen receptor (ER), but the underlying molecular mechanisms are unclear. The present study aimed to identify ER target genes involved in estrogen-caused gonadal feminization in Xenopus laevis. Based on our recent transcriptomic data that 10 nM 17β-estradiol (E2) altered gene transcription in feminizing gonads of male X. laevis at NF stages 48, 50, and 52, we searched estrogen response element (ERE) using the Dragon ERE Finder software in the promoter region of all the E2-regulated genes. As a result, 163 genes containing ERE sequence were identified as predicted ER target genes at NF stage 50 (on the 14th day postfertilization), a crucial stage for gonadal feminization. Then, some of these predicted ER target genes were further investigated, mainly including the genes that were suggested to be involved in E2-caused gonadal feminization and genes being dramatically up or down-regulated by E2. Fifteen genes were demonstrated to be responsive to E2, in turn ER antagonist blocked the E2-regulated transcription. Finally, we identified 10 genes that can bind to ERα by a chromatin immunoprecipitation-qPCR. Taken together, we identified the 10 genes that contain predicted ERE sequences, are responsive to estrogen and ER antagonist, and have ability to bind to ER as ER target genes, including pglyrp2, apoa1, fgb, tdo2, ca6, nags, cpb2, tmprss6, nudc, zwilch. Our results could help to improve the understanding of the molecular mechanisms for gonadal feminization caused by estrogenic endocrine disrupting chemicals in X. laevis, and even in other species.

miR-452 Reverses Abnormal Glycosylation Modification of ERα and Estrogen Resistance in TNBC (Triple-Negative Breast Cancer) Through Targeting UGT1A1

Background: The breast epithelial cells in patients with triple-negative breast cancer (TNBC) actually have specific estrogen receptor (ER) expression, and the abnormal glycosylation of UGT1A1 in TNBC cells resulted in abnormal expression and function of ERα through regulating the modification of ERα. Therefore, our study targets the role of UGT1A1 expression, then glycosylation modification of ERα (estrogen receptor α) and estrogen resistance in development of TNBC.

Methods: The differential expression of mRNA and miRNA in TNBC tissues was tested. Luciferase activity was analyzed in TNBC cells treated with miR-452. Moreover, the human mammary gland and TNBC cell lines were dealt with estrogen and miR-452 or its inhibitors, then proliferation ability was further determined. Moreover, the role of interaction between UGT1A1 and ERα in the glycosylation modification of ERα and UGT activity, and metabolism of estrogen were assessed. The effects of miR-452 on TNBC by improving abnormal glycosylation modification of ERα by targeting UGT1A1 and estrogen resistance were studied in vitro and in vivo.

Results: The expression level of UGT1A1 in TNBC tumor tissues was higher than its matched para-tumorous tissues, but the miR-452 expression was opposite. The glycosylation modification site of ERα expressed in TNBC cells was different from that of normal mammary epithelial cells. The estrogen 17β-estradiol (E2) significantly promoted mitotic entry of TNBC cells. The interaction between UGT1A1 and ERα affected the expression level of each other, as well as the UGT enzyme activity and proliferation of TNBC cells. UGT1A1 induced production of intracellular estrogens and TNBC proliferation, but it could be reversed by overexpression of ERα. Upregulation of ERα caused the downregulation of UGT1A1 and marked decrease of intracellular estrogen products, and then suppressed TNBC proliferation. Moreover, UGT1A1 was the target gene of miR-452; miR-452 antagomir restrained TNBC xenograft.

Conclusion: Our results demonstrated that estrogen was a positive factor in the proliferation of TNBC cells at onset of mitosis through accentuating the expression and enzyme activity of UGT1A1. However, miR-452 targeted to UGT1A1, then regulated glycosylation modification of ERα, estrogen metabolism, and TNBC development associated with estrogen resistance.

PRIM1 promotes the proliferation of hepatocellular carcinoma cells in vitro and in vivo

PRIM1 plays an important role during oncogenesis, however it has never been reported in liver cancer, and thus our objective is to explore the role of PRIM1 in liver cancer. We selected RNAseq data of 50 paired liver cancer samples from the Cancer Gene Atlas (TCGA), and then bioinformatics methods and Mann-Whitney U test were used to analyze the correlation between PRIM1 and the clinical pathological stage of liver cancer. Quantitative polymerase chain reaction (QPCR) was used to detect mRNA expression of PRIM1 in BEL-7404, BCL-7402, HepG2 and SMMC-7721 cell lines. LV-PRIM1-RNAi was transfected into BEL-7404 and SMMC-7721 cells by lentivirus, and then Celigo imaging cytometer, Caspase3/7 Assay, flow cytometry and MTT assay were used to detect the proliferation and apoptosis of BEL-7404 and SMMC-7721 cells with ≥50% gene reduction rate after lentivirus transfection detected by QPCR. BEL-7404 and SMMC-7721 carrying PRIM1 gene were used for oncogenesis in vitro to observe the weight and fluorescence intensity of the tumor. Bioinformatics method was used to obtain the information about PRIM1 gene, and the correlation between PRIM1 and clinical pathological stage of liver cancer was analyzed by Mann-Whitney U test. QPCR results showed that PRIM1 was expressed in BEL-7404, BCL-7402, HepG2 and SMMC-7721 cell lines, which was highest in BCL-7404 cell line. Celigo imaging cytometer, Caspase3/7 Assay, flow cytometry and MTT assay showed that the proliferative ability of BEL-7404 and SMMC-7721 were decreased after LV-PRIM1-RNAi transfection. Furthermore, the weight and the fluorescence intensity of the tumors in vitro formed by LV-PRIM1-RNAi cells on SCID mice were decreased. So, interference of PRIM1 expression can inhibit the proliferation of BEL-7404 and SMMC-7721 cells, as well as induce the apoptosis of liver cancer cells.

DNA Primase Subunit 1 Expression in Hepatocellular Carcinoma and Its Clinical Implication

DNA Primase Subunit 1 (PRIM1) is crucial for cancer development and progression. However, there remains a lack of comprehension concerning the clinical implication of PRIM1 in HCC. Here, aberrant expression of PRIM1 was identified in HCC according to available databases. The prognostic value of PRIM1 in patients presenting with HCC was further assessed based on TCGA data. Gene set enrichment analysis (GSEA) was subsequently conducted to investigate the potential function of PRIM1. Additionally, the correlations between tumor-infiltrating immune cells (TIICs) and PRIM1 expression were evaluated. The data from TCGA, GEO, ONCOMINE, and HCCDB databases illustrated that PRIM1 was overexpressed in HCC tissues, compared to normal liver tissues (all p < 0.05). Kaplan-Meier analysis revealed that high PRIM1 expression in HCC was closely correlated with worse overall survival (p < 0.05). The univariate and multivariate analyses illustrated that PRIM1 expression was an independent novel prognostic indicator in HCC. Additionally, the area under the receiver operating characteristic (AUROC) curve for PRIM1 reached 0.8651, indicating the diagnostic significance of PRIM1 in patients with HCC. GSEA showed that PRIM1 overexpression was significantly enriched in several tumor-related signaling pathways. Besides, TIIC analysis clarified the association between PRIM1 expression and TIICs in HCC. The findings disclose that PRIM1 profoundly implicated in promoting tumorigenesis might work as a desirable biomarker for HCC.

Comprehensive analysis reveals a metabolic ten-gene signature in hepatocellular carcinoma

Background

Due to the complicated molecular and cellular heterogeneity in hepatocellular carcinoma (HCC), the morbidity and mortality still remains high level in the world. However, the number of novel metabolic biomarkers and prognostic models could be applied to predict the survival of HCC patients is still small. In this study, we constructed a metabolic gene signature by systematically analyzing the data from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC).

Methods

Differentially expressed genes (DEGs) between tumors and paired non-tumor samples of 50 patients from TCGA dataset were calculated for subsequent analysis. Univariate cox proportional hazard regression and LASSO analysis were performed to construct a gene signature. The Kaplan–Meier analysis, time-dependent receiver operating characteristic (ROC), Univariate and Multivariate Cox regression analysis, stratification analysis were used to assess the prognostic value of the gene signature. Furthermore, the reliability and validity were validated in four types of testing cohorts. Moreover, the diagnostic capability of the gene signature was investigated to further explore the clinical significance. Finally, Go enrichment analysis and Gene Set Enrichment Analysis (GSEA) have been performed to reveal the different biological processes and signaling pathways which were active in high risk or low risk group.

Results

Ten prognostic genes were identified and a gene signature were constructed to predict overall survival (OS). The gene signature has demonstrated an excellent ability for predicting survival prognosis. Univariate and Multivariate analysis revealed the gene signature was an independent prognostic factor. Furthermore, stratification analysis indicated the model was a clinically and statistically significant for all subgroups. Moreover, the gene signature demonstrated a high diagnostic capability in differentiating normal tissue and HCC. Finally, several significant biological processes and pathways have been identified to provide new insights into the development of HCC.

Conclusion

The study have identified ten metabolic prognostic genes and developed a prognostic gene signature to provide more powerful prognostic information and improve the survival prediction for HCC.

Structural Basis for Inhibition of Human Primase by Arabinofuranosyl Nucleoside Analogues Fludarabine and Vidarabine

Nucleoside analogues are widely used in clinical practice as chemotherapy drugs. Arabinose nucleoside derivatives such as fludarabine are effective in the treatment of patients with acute and chronic leukemias and non-Hodgkin's lymphomas. Although nucleoside analogues are generally known to function by inhibiting DNA synthesis in rapidly proliferating cells, the identity of their in vivo targets and mechanism of action are often not known in molecular detail. Here we provide a structural basis for arabinose nucleotide-mediated inhibition of human primase, the DNA-dependent RNA polymerase responsible for initiation of DNA synthesis in DNA replication. Our data suggest ways in which the chemical structure of fludarabine could be modified to improve its specificity and affinity toward primase, possibly leading to less toxic and more effective therapeutic agents.

The Y-linked proto-oncogene TSPY contributes to poor prognosis of the male hepatocellular carcinoma patients by promoting the pro-oncogenic and suppressing the anti-oncogenic gene expression

Background

Liver cancer is one of the major causes of cancer death worldwide, with significantly higher incidence and mortality among the male patients. Although sex hormones and their receptors could contribute to such sex differences, the story is incomplete. Genes on the male-specific region of the Y chromosome could play a role(s) in this cancer. TSPY is the putative gene for the gonadoblastoma locus on the Y chromosome (GBY) that is ectopically expressed in a subset of male hepatocellular carcinomas (HCCs). Although various studies showed that TSPY expression is associated with poor prognosis in the patients and its overexpression promotes cell proliferation of various cancer cell lines, it remains unclear how TSPY contributes to the clinical outcomes of the HCC patients. Identifying the downstream genes and pathways of TSPY actions would provide novel insights on its contribution(s) to male predominance in this deadly cancer.

Results

To determine the effects of TSPY on HCC, a TSPY transgene was introduced to the HCC cell line, HuH-7, and studied with RNA-Seq transcriptome analysis. The results showed that TSPY upregulates various genes associated with cell-cycle and cell-viability, and suppresses cell-death related genes. To correlate the experimental observations with those of clinical specimens, transcriptomes of male HCCs with high TSPY expression were analyzed with reference to those with silent TSPY expression from the Cancer Genome Atlas (TCGA). The comparative analysis identified 49 genes, which showed parallel expression patterns between HuH-7 cells overexpressing TSPY and clinical specimens with high TSPY expression. Among these 49 genes, 16 likely downstream genes could be associated with survival rates in HCC patients. The major upregulated targets were cell-cycle related genes and growth factor receptor genes, including CDC25B and HMMR, whose expression levels are negatively correlated with the patient survival rates. In contrast, PPARGC1A, SLC25A25 and SOCS2 were downregulated with TSPY expression, and possess favorable prognoses for HCC patients.

Conclusion

We demonstrate that TSPY could exacerbate the oncogenesis of HCC by differentially upregulate the expression of pro-oncogenic genes and downregulate those of anti-oncogenic genes in male HCC patients, thereby contributing to the male predominance in this deadly cancer.

Electronic supplementary material

The online version of this article (10.1186/s13578-019-0287-x) contains supplementary material, which is available to authorized users.

Inactivation of PRIM1 Function Sensitizes Cancer Cells to ATR and CHK1 Inhibitors

The phosphoinositide 3-kinase-related kinase ATR is a central regulator of the DNA damage response. Its chemical inhibition eliminates subsets of cancer cells in various tumor types. This effect is caused at least partly by the synthetically lethal relationship between ATR and certain DNA repair genes. In a previous screen using an siRNA library against DNA repair genes, we identified PRIM1, a part of the polymerase α-primase complex, as acting synthetically lethal with ATR. Applying a genetic ATR knock-in model of colorectal cancer cells, we confirmed that PRIM1 depletion inhibited proliferation of ATR-deficient cells and excluded artifacts due to clonal variation using an ATR reexpressing cell clone. We expanded these data by demonstrating in different cell lines that also chemical inhibition of ATR or its main effector kinase CHK1 reduces proliferation upon depletion of PRIM1. Mechanistically, PRIM1 depletion in ATR-deficient cells caused S-phase stasis in the absence of increased DNA damage followed by Wee1-mediated activation of caspase 8 and apoptosis. As PRIM1 inactivation sensitizes cancer cells to ATR and CHK1 inhibitors, mutations in PRIM1 or other components of the polymerase α-primase complex could represent novel targets for individualized tumor therapeutic approaches using ATR/CHK1 inhibitors, as has been previously demonstrated for POLD1, the catalytic subunit of polymerase δ.