Meta-analysis of gene expression profiles indicates genes in spliceosome pathway are up-regulated in hepatocellular carcinoma (HCC)

Oncofetal SNRPE promotes HCC tumorigenesis by regulating the FGFR4 expression through alternative splicing

BACKGROUND

Due to insufficient knowledge about key molecular events, Hepatocellular carcinoma (HCC) lacks effective treatment targets. Spliceosome-related genes were significantly altered in HCC. Oncofetal proteins are ideal tumor therapeutic targets. Screening of differentially expressed Spliceosome-related oncofetal protein in embryonic liver development and HCC helps discover effective therapeutic targets for HCC.

METHODS

Differentially expressed spliceosome genes were analysis in fetal liver and HCC through bioinformatics analysis. Small nuclear ribonucleoprotein polypeptide E (SNRPE) expression was detected in fetal liver, adult liver and HCC tissues. The role of SNRPE in HCC was performed multiple assays in vitro and in vivo. SNRPE-regulated alternative splicing was recognized by RNA-Seq and confirmed by multiple assays.

RESULTS

We herein identified SNRPE as a crucial oncofetal splicing factor, significantly associated with the adverse prognosis of HCC. SOX2 was identified as the activator for SNRPE reactivation. Efficient knockdown of SNRPE resulted in the complete cessation of HCC tumorigenesis and progression. Mechanistically, SNRPE knockdown reduced FGFR4 mRNA expression by triggering nonsense-mediated RNA decay. A partial inhibition of SNRPE-induced malignant progression of HCC cells was observed upon FGFR4 knockdown.

CONCLUSIONS

Our findings highlight SNRPE as a novel oncofetal splicing factor and shed light on the intricate relationship between oncofetal splicing factors, splicing events, and carcinogenesis. Consequently, SNRPE emerges as a potential therapeutic target for HCC treatment. Model of oncofetal SNRPE promotes HCC tumorigenesis by regulating the AS of FGFR4 pre-mRNA.

The role of alternative pre-mRNA splicing in cancer progression

Alternative pre-mRNA splicing is a critical mechanism that generates multiple mRNA from a single gene, thereby increasing the diversity of the proteome. Recent research has highlighted the significance of specific splicing isoforms in cellular processes, particularly in regulating cell numbers. In this review, we examine the current understanding of the role of alternative splicing in controlling cancer cell growth and discuss specific splicing factors and isoforms and their molecular mechanisms in cancer progression. These isoforms have been found to intricately control signaling pathways crucial for cell cycle progression, proliferation, and apoptosis. Furthermore, studies have elucidated the characteristics and functional importance of splicing factors that influence cell numbers. Abnormal expression of oncogenic splicing isoforms and splicing factors, as well as disruptions in splicing caused by genetic mutations, have been implicated in the development and progression of tumors. Collectively, these findings provide valuable insights into the complex interplay between alternative splicing and cell proliferation, thereby suggesting the potential of alternative splicing as a therapeutic target for cancer.

PPIH gene regulation system and its prognostic significance in hepatocellular carcinoma: a comprehensive analysis

BACKGROUND

Peptidyl-prolyl isomerase H (PPIH) is a member of the cyclophilin protein family, which functions as a molecular chaperone and is involved in the splicing of pre-mRNA. According to reports, the malignant progression of HCC related to hepatitis B virus (HBV) is tightly associated with RNA-binding proteins. Nevertheless, there is no research on PPIH expression or its function in the occurrence and progression of HCC.

RESULTS

We are the first to reveal that the mRNA and protein levels of Ppih are substantially overexpressed in HCC, as the outcomes show. A significant correlation existed between enriched expression of Ppih within HCC and more advanced, poorly differentiated, and TP53-mutated tumors.

CONCLUSION

These findings, which suggest that Ppih may serve as a predictive biomarker for people with HCC, serve as a starting point for further investigation into the function of Ppih in the progression of carcinogenesis.

METHODS

Accordingly, we utilized clinical samples and bioinformatics analysis to assess Ppih's mRNA, protein expression, and gene regulatory system in HCC. Additionally, Wilcoxon signed-rank testing and logistic regression were utilized to inspect the association between clinicopathological factors and Ppih. Clinical pathological traits linked to overall survival (OS) among HCC patients were examined via TCGA data via Cox regression and the Kaplan-Meier approach. Additionally, via TCGA data collection, gene set enrichment assessment was also conducted.

USP39 promotes hepatocellular carcinogenesis through regulating alternative splicing in cooperation with SRSF6/HNRNPC

Abnormal alternative splicing (AS) caused by alterations in spliceosomal factors is implicated in cancers. Standard models posit that splice site selection is mainly determined by early spliceosomal U1 and U2 snRNPs. Whether and how other mid/late-acting spliceosome components such as USP39 modulate tumorigenic splice site choice remains largely elusive. We observed that hepatocyte-specific overexpression of USP39 promoted hepatocarcinogenesis and potently regulated splice site selection in transgenic mice. In human liver cancer cells, USP39 promoted tumor proliferation in a spliceosome-dependent manner. USP39 depletion deregulated hundreds of AS events, including the oncogenic splice-switching of KANK2. Mechanistically, we developed a novel RBP-motif enrichment analysis and found that USP39 modulated exon inclusion/exclusion by interacting with SRSF6/HNRNPC in both humans and mice. Our data represented a paradigm for the control of splice site selection by mid/late-acting spliceosome proteins and their interacting RBPs. USP39 and possibly other mid/late-acting spliceosome proteins may represent potential prognostic biomarkers and targets for cancer therapy.

MKRN1 promotes colorectal cancer metastasis by activating the TGF-β signalling pathway through SNIP1 protein degradation

BACKGROUND

The Makorin ring finger protein 1 (MKRN1) gene, also called RNF61, is located on the long arm of chromosome 7 and is a member of the RING finger protein family. The E3 ubiquitin ligase MKRN1 is closely linked to tumour development, but the exact mechanism needs to be elucidated. In this study, we aimed to investigate the specific mechanism and role of MKRN1 in colorectal cancer (CRC) development.

METHODS

MKRN1 expression in CRC was analysed using the Cancer Cell Line Encyclopaedia and the Cancer Genome Atlas (TCGA) databases. Rectal tumour tissues were frozen to explore the MKRN1 expression in CRC and its clinical significance. The impact of MKRN1 on CRC cell proliferation and migration was observed using CCK8, colony formation, wound healing, and transwell assays. A combination of MKRN1 quantitative proteomics, ubiquitination modification omics analysis, and a string of in vitro and in vivo experiments revealed the potential mechanisms by which MKRN1 regulates CRC metastasis.

RESULTS

MKRN1 expression was significantly elevated in CRC tissues compared to paracancerous tissues and was positively linked with prognosis (P < 0.01). MKRN1 downregulation inhibits CRC cell proliferation, migration, and invasion. Conversely, MKRN1 overexpression promotes the proliferation, migration, and invasion of CRC cells. Mechanistically, MKRN1 induces epithelial-mesenchymal transition (EMT) in CRC cells via ubiquitination and degradation of Smad nuclear-interacting protein 1 (SNIP1). Furthermore, SNIP1 inhibits transforming growth factor-β (TGF-β) signalling, and MKRN1 promotes TGF-β signalling by degrading SNIP1 to induce EMT in CRC cells. Finally, using conditional knockout mice, intestinal lesions and metastatic liver microlesions were greatly reduced in the intestinal knockout MKRN1 group compared to that in the control group.

CONCLUSIONS

High MKRN1 levels promote TGF-β signalling through ubiquitination and degradation of SNIP1, thereby facilitating CRC metastasis, and supporting MKRN1 as a CRC pro-cancer factor. The MKRN1/SNIP1/TGF-β axis may be a potential therapeutic target in CRC.

Identification of MYB gene family and functional analysis of GhMYB4 in cotton (Gossypium spp.)

Myeloblastosis (MYB) transcription factors (TFs) form a large gene family involved in a variety of biological processes in plants. Little is known about their roles in the development of cotton pigment glands. In this study, 646 MYB members were identified in Gossypium hirsutum genome and phylogenetic classification was analyzed. Evolution analysis revealed assymetric evolution of GhMYBs during polyploidization and sequence divergence of MYBs in G. hirustum was preferentially happend in D sub-genome. WGCNA (weighted gene co-expression network analysis) showed that four modules had potential relationship with gland development or gossypol biosynthesis in cotton. Eight differentially expressed GhMYB genes were identified by screening transcriptome data of three pairs of glanded and glandless cotton lines. Of these, four were selected as candidate genes for cotton pigment gland formation or gossypol biosynthesis by qRT-PCR assay. Silencing of GH_A11G1361 (GhMYB4) downregulated expression of multiple genes in gossypol biosynthesis pathway, indicating it could be involved in gossypol biosynthesis. The potential protein interaction network suggests that several MYBs may have indirect interaction with GhMYC2-like, a key regulator of pigment gland formation. Our study was the systematic analysis of MYB genes in cotton pigment gland development, providing candidate genes for further study on the roles of cotton MYB genes in pigment gland formation, gossypol biosynthesis and future crop plant improvement.

Bioinformatics analysis and experimental validation of a novel autophagy-related signature relevant to immune infiltration for recurrence prediction after curative hepatectomy

Hepatocellular carcinoma (HCC) remains imposing an enormous economic and healthcare burden worldwide. In this present study, we constructed and validated a novel autophagy-related gene signature to predict the recurrence of HCC patients. A total of 29 autophagy-related differentially expressed genes were identified. A five-gene signature (CLN3, HGF, TRIM22, SNRPD1, and SNRPE) was constructed for HCC recurrence prediction. Patients in high-risk groups exhibited a significantly poor prognosis compared with low-risk patients both in the training set (GSE14520 dataset) and the validation set (TCGA and GSE76427 dataset). Multivariate cox regression analysis demonstrated that the 5-gene signature was an independent risk factor for recurrence-free survival (RFS) in HCC patients. The nomograms incorporating 5-gene signature and clinical prognostic risk factors were able to effectively predict RFS. KEGG and GSEA analysis revealed that the high-risk group was enriched with multiple oncology characteristics and invasive-related pathways. Besides, the high-risk group had a higher level of immune cells and higher levels of immune checkpoint-related gene expression in the tumor microenvironment, suggesting that they might be more likely to benefit from immunotherapy. Finally, the immunohistochemistry and cell experiments confirmed the role of SNRPE, the most significant gene in the gene signature. SNRPE was significantly overexpressed in HCC. After SNRPE knockdown, the proliferation, migration and invasion ability of the HepG2 cell line were significantly inhibited. Our study established a novel five-gene signature and nomogram to predict RFS of HCC, which may help in clinical decision-making for individual treatment.

TRA2: The dominant power of alternative splicing in tumors

The dysregulation of alternative splicing (AS) is frequently found in cancer and considered as key markers for cancer progression and therapy. Transformer 2 (TRA2), a nuclear RNA binding protein, consists of transformer 2 alpha homolog (TRA2A) and transformer 2 beta homolog (TRA2B), and plays a role in the regulation of pre-mRNA splicing. Growing evidence has been provided that TRA2A and TRA2B are dysregulated in several types of tumors, and participate in the regulation of proliferation, migration, invasion, and chemotherapy resistance in cancer cells through alteration of AS of cancer-related genes. In this review, we highlight the role of TRA2 in tumorigenesis and metastasis, and discuss potential molecular mechanisms how TRA2 influences tumorigenesis and metastasis via controlling AS of pre-mRNA. We propose that TRA2Ais a novel biomarker and therapeutic target for cancer progression and therapy.

A Five-LLPS Gene Risk Score Prognostic Signature Predicts Survival in Hepatocellular Carcinoma

Background

Primary liver cancer, dominated by hepatocellular carcinoma (HCC), is one of the most common cancer types and the third leading cause of cancer death in 2020. Previous studies have shown that liquid-liquid phase separation (LLPS) plays an important role in the occurrence and development of cancer including HCC, but its influence on the patient prognosis is still unknown. It is necessary to explore the effect of LLPS genes on prognosis to accurately forecast the prognosis of HCC patients and identify relevant targeted therapeutic sites.

Methods

Using The Cancer Genome Atlas dataset and PhaSepDB dataset, we identified LLPS genes linked to the overall survival (OS) of HCC patients. We applied Least Absolute Shrinkage and Selection Operator (LASSO) Cox penalized regression analysis to choose the best genes for the risk score prognostic signature. We then analysed the validation dataset and evaluated the effectiveness of the risk score prognostic signature. Finally, we performed quantitative real-time PCR experiments to validate the genes in the prognostic signature.

Results

We identified 43 differentially expressed LLPS genes that were associated with the OS of HCC patients. Five of these genes (BMX, FYN, KPNA2, PFKFB4, and SPP1) were selected to generate a prognostic risk score signature. Patients in the low-risk group were associated with better OS than those in the high-risk group in both the training dataset and the validation dataset. We found that BMX and FYN had lower expression levels in HCC tumour tissues, whereas KPNA2, PFKFB4, and SPP1 had higher expression levels in HCC tumour tissues. The validation demonstrated that the five-LLPS gene risk score signature has the capability of predicting the OS of HCC patients.

Conclusion

Our study constructed a five-LLPS gene risk score signature that can be applied as an effective and convenient prognostic tool. These five genes might serve as potential targets for therapy and the treatment of HCC.

Establishment and validation of a ferroptosis-related prognostic signature for hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is the most common type of primary liver cancer with high heterogeneity. The prognosis of HCC is quite poor and the prognostic prediction also has challenges. Ferroptosis is recently recognized as a kind of iron-dependent cell death, which is involved in tumor progression. However, further study is needed to validate the influence of drivers of ferroptosis (DOFs) on the prognosis of HCC.

Methods

The FerrDb database and the Cancer Genome Atlas (TCGA) database were applied to retrieve DOFs and information of HCC patients respectively. HCC patients were randomly divided into training and testing cohorts with a 7:3 ratio. Univariate Cox regression, LASSO and multivariate Cox regression analyses were carried out to identify the optimal prognosis model and calculate the risk score. Then, univariate and multivariate Cox regression analyses were performed to assess the independence of the signature. At last, gene functional, tumor mutation and immune-related analyses were conducted to explore the underlying mechanism. Internal and external databases were used to confirm the results. Finally, the tumor tissue and normal tissue from HCC patients were applied to validate the gene expression in the model.

Results

Five genes were identified to develop as a prognostic signature in the training cohort relying on the comprehensive analysis. Univariate and multivariate Cox regression analyses confirmed that the risk score was able to be an independent factor for the prognosis of HCC patients. Low-risk patients showed better overall survival than high-risk patients. Receiver operating characteristic (ROC) curve analysis confirmed the signature's predictive capacity. Furthermore, internal and external cohorts were consistent with our results. There was a higher proportion of nTreg cell, Th1 cell, macrophage, exhausted cell and CD8⁺T cell in the high-risk group. The Tumor Immune Dysfunction and Exclusion (TIDE) score suggested that high-risk patients could respond better to immunotherapy. Besides, the experimental results showed that some genes were differentially expressed between tumor and normal tissues.

Conclusion

In summary, the five ferroptosis gene signature showed potential in prognosis of patients with HCC and could also be regarded as a value biomarker for immunotherapy response in these patients.

scBPGRN: Integrating single-cell multi-omics data to construct gene regulatory networks based on BP neural network

The deterioration and metastasis of cancer involve various aspects of genomic changes, including genomic DNA changes, epigenetic modifications, gene expression, and other complex interactions. Therefore, integrating single-cell multi-omics data to construct gene regulatory networks containing more omics information is of great significance for understanding the pathogenesis of cancer. In this article, an algorithm integrating single-cell RNA sequencing data and DNA methylation data to construct a gene regulatory network based on the back-propagation (BP) neural network (scBPGRN) is proposed. This algorithm uses biweight extreme correlation coefficients to measure the correlation between factors and uses neural networks to calculate generalized weights to construct gene regulation networks. Finally, the node strength is calculated to identify the genes associated with cancer. We apply the scBPGRN algorithm to hepatocellular carcinoma (HCC) data. We construct a regulatory network and identify top-ranked genes, such as MYCBP, KLHL35, PRKCZ, and SERPINA6, as the key HCC-related genes. We analyze the top 100 genes, and the HCC-related genes are concentrated in the top 20. In addition, the single cell data is found to consist of two subpopulations. We also apply scBPGRN to two subpopulations. We analyze the top 50 genes in them, and the HCC-related genes are concentrated in the top 20. The consequences of functional enrichment analysis indicate that the gene regulatory network we have constructed is valid. Our results have been verified in several pieces of literature. This study provides a reference for the integration of single-cell multi-omics data to construct gene regulatory networks.

Dysregulation of global circular RNA abundance regulated by spliceosomes predicts prognosis in hepatocellular carcinoma

CircRNAs have been reported to play crucial roles in tumor progression and recurrence, showing potential as biomarkers in cancer. However, the global abundance of circRNA and their involvement in hepatocellular carcinoma (HCC) development have not been fully explored. Whole transcriptome sequencing was performed on tumor and peritumor from 60 patients with HCC to quantify the expression of circRNAs, and the global circRNA abundance was calculated by circRNA index (CRI). Gene-set enrichment analysis and weighted gene co-expression network analysis were used to reveal the biological signaling pathways associated with the global circRNA abundance. The correlation between the global circRNA abundance and the infiltration level of CD8⁺ T cells was explored by immunohistochemical assays. Small interfering RNA was used to knock down the pre-messenger RNA spliceosome in HCC cell lines to verify the regulation of spliceosome in global circRNA abundance. We found that dysregulation of global circRNA abundance in both tumor and peritumor could lead to worse prognosis. The immunohistochemical assay further revealed that the dysregulation of global circRNA abundance in both tumor and peritumor would obstruct the CD8⁺ T cells from invading into the tumor, which might explain its correlation with HCC prognosis. We also demonstrated that the spliceosome genes were the main factors to regulate the global circRNA abundance in HCC, and these results were also confirmed by knockdown experiments. Conclusion: This study revealed the association between the global circRNA abundance and patients' prognosis and its underlying mechanism.

Significance of Spliceosome-Related Genes in the Prediction of Prognosis and Treatment Strategies for Lung Adenocarcinoma

BACKGROUND

The leading cause of cancer-related fatalities globally is lung cancer; lung adenocarcinoma (LUAD) is the most common histological type in it. The spliceosome plays an important role in a majority of malignancies. However, it is yet unclear how spliceosome-related genes affect patients with LUAD in terms of treatment course and prognosis.

METHODS

Spliceosome-related genes were assessed from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database to obtain clinical information and gene expression in patients with LUAD. A spliceosome-related gene signature and prognostic model were constructed by using the least absolute shrinkage and selection operator (LASSO), time-dependent receiver operating characteristic (ROC), and nomogram. Immune infiltrate levels, mutation analysis, and pathway enrichment were predicted potential mechanisms of the signature by using single-sample gene set enrichment analysis (ssGSEA), Gene Set Cancer Analysis (GSCA) database, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Ontology (GO) database. Then, a protein-protein interaction (PPI) network and transcription factor- (TF-) hub gene and drug mining network were also established by Cytoscape software.

RESULTS

Firstly, we constructed a prognostic model for 11 spliceosome signature genes. Based on the prognostic risk score, we stratified patients with LUAD into high- and low-risk groups. The high- and low-risk groups were closely related to the OS, tumor immune infiltration level, immune checkpoint molecules, and tumor mutation burden (TMB) of LUAD patients. Based on PPI networks, we also predict relevant TF genes that may regulate signature prognostic genes. Finally, drugs including oxaliplatin, arsenic trioxide, cisplatin, and sunitinib were excavated for the treatment of the 11 spliceosome signature genes in LUAD patients.

CONCLUSION

In conclusion, this study is the first to explore the importance of spliceosome-related genes in the prognosis and treatment of LUAD. Through our study, we have innovatively provided potential prognosis genes and new therapeutic drug targets for the treatment of LUAD patients.

SF3B4 Depletion Retards the Growth of A549 Non-Small Cell Lung Cancer Cells via UBE4B-Mediated Regulation of p53/p21 and p27 Expression

Splicing factor B subunit 4 (SF3B4), a component of the U2-pre-mRNA spliceosomal complex, contributes to tumorigenesis in several types of tumors. However, the oncogenic potential of SF3B4 in lung cancer has not yet been determined. The in vivo expression profiles of SF3B4 in non-small cell lung cancer (NSCLC) from publicly available data revealed a significant increase in SF3B4 expression in tumor tissues compared to that in normal tissues. The impact of SF3B4 deletion on the growth of NSCLC cells was determined using a siRNA strategy in A549 lung adenocarcinoma cells. SF3B4 silencing resulted in marked retardation of the A549 cell proliferation, accompanied by the accumulation of cells at the G0/G1 phase and increased expression of p27, p21, and p53. Double knockdown of SF3B4 and p53 resulted in the restoration of p21 expression and partial recovery of cell proliferation, indicating that the p53/p21 axis is involved, at least in part, in the SF3B4-mediated regulation of A549 cell proliferation. We also provided ubiquitination factor E4B (UBE4B) is essential for p53 accumulation after SF3B4 depletion based on followings. First, co-immunoprecipitation showed that SF3B4 interacts with UBE4B. Furthermore, UBE4B levels were decreased by SF3B4 depletion. UBE4B depletion, in turn, reproduced the outcome of SF3B4 depletion, including reduction of polyubiquitinated p53 levels, subsequent induction of p53/p21 and p27, and proliferation retardation. Collectively, our findings indicate the important role of SF3B4 in the regulation of A549 cell proliferation through the UBE4B/p53/p21 axis and p27, implicating the therapeutic strategies for NSCLC targeting SF3B4 and UBE4B.

RNA splicing: a dual-edged sword for hepatocellular carcinoma

RNA splicing is the fundamental process that brings diversity at the transcriptome and proteome levels. The spliceosome complex regulates minor and major processes of RNA splicing. Aberrant regulation is often associated with different diseases, including diabetes, stroke, hypertension, and cancer. In the majority of cancers, dysregulated alternative RNA splicing (ARS) events directly affect tumor progression, invasiveness, and often lead to poor survival of the patients. Alike the rest of the gastrointestinal malignancies, in hepatocellular carcinoma (HCC), which alone contributes to ~ 75% of the liver cancers, a large number of ARS events have been observed, including intron retention, exon skipping, presence of alternative 3'-splice site (3'SS), and alternative 5'-splice site (5'SS). These events are reported in spliceosome and non-spliceosome complexes genes. Molecules such as MCL1, Bcl-X, and BCL2 in different isoforms can behave as anti-apoptotic or pro-apoptotic, making the spliceosome complex a dual-edged sword. The anti-apoptotic isoforms of such molecules bring in resistance to chemotherapy or cornerstone drugs. However, in contrast, multiple malignant tumors, including HCC that target the pro-apoptotic favoring isoforms/variants favor apoptotic induction and make chemotherapy effective. Herein, we discuss different splicing events, aberrations, and antisense oligonucleotides (ASOs) in modulating RNA splicing in HCC tumorigenesis with a possible therapeutic outcome.

Transformer 2 alpha homolog is a downstream gene of hypoxia-inducible factor 1 subunit alpha and is involved in the progression of pancreatic cancer

Intratumoral hypoxia is a common feature of pancreatic cancer (PC) and also plays a role in its progression. However, hypoxia-regulated signatures in PC are still not completely understood. This study aimed to identify core hypoxia-associated genes and determine their underlying molecular mechanisms in PC cells. Transformer 2 alpha homolog (TRA2A) was found to be an important hypoxia-associated gene, which was upregulated in PC tissues and in PC cells cultured under hypoxia. High TRA2A expression was associated with advanced stage, poor differentiation, and lymph node metastasis. Under normoxic and hypoxic conditions, knockdown of TRA2A both markedly suppressed PC cell proliferation and motility in vitro and in vivo, as well as activation of the AKT pathway. Hypoxia-inducible factor 1 subunit alpha (HIF1α) upregulated the transcription of TRA2A by directly binding to its promoter. TRA2A showed a co-expression relationship with HIF1α in PC tissues. Overexpression of TRA2A alleviated the pro-inhibitive functions of HIF1α-inhibition on PC cell proliferation and motility under hypoxia. In conclusion, TRA2A is a crucial downstream gene of HIF1α that accelerates the proliferation and motility of PC cells. TRA2A may be a novel and practical molecular target for investigating the hypoxic response of PC cells.Abbreviations: TRA2A, transformer 2A protein; PC, pancreatic cancer; HIF1α, hypoxia-inducible factor 1-alpha; GEO, Gene Expression Omnibus; IHC, immunohistochemical staining.

Association of Glyphosate Exposure with Blood DNA Methylation in a Cross-Sectional Study of Postmenopausal Women

BACKGROUND

Glyphosate is the most commonly used herbicide in the world and is purported to have a variety of health effects, including endocrine disruption and an elevated risk of several types of cancer. Blood DNA methylation has been shown to be associated with many other environmental exposures, but to our knowledge, no studies to date have examined the association between blood DNA methylation and glyphosate exposure.

OBJECTIVE

We conducted an epigenome-wide association study to identify DNA methylation loci associated with urinary glyphosate and its metabolite aminomethylphosphonic acid (AMPA) levels. Secondary goals were to determine the association of epigenetic age acceleration with glyphosate and AMPA and develop blood DNA methylation indices to predict urinary glyphosate and AMPA levels.

METHODS

For 392 postmenopausal women, white blood cell DNA methylation was measured using the Illumina Infinium MethylationEPIC BeadChip array. Glyphosate and AMPA were measured in two urine samples per participant using liquid chromatography-tandem mass spectrometry. Methylation differences at the probe and regional level associated with glyphosate and AMPA levels were assessed using a resampling-based approach. Probes and regions that had an false discovery rate q < 0.1 in ≥ 90 % of 1,000 subsamples of the study population were considered differentially methylated. Differentially methylated sites from the probe-specific analysis were combined into a methylation index. Epigenetic age acceleration from three epigenetic clocks and an epigenetic measure of pace of aging were examined for associations with glyphosate and AMPA.

RESULTS

We identified 24 CpG sites whose methylation level was associated with urinary glyphosate concentration and two associated with AMPA. Four regions, within the promoters of the MSH4, KCNA6, ABAT, and NDUFAF2/ERCC8 genes, were associated with glyphosate levels, along with an association between ESR1 promoter hypomethylation and AMPA. The methylation index accurately predicted glyphosate levels in an internal validation cohort. AMPA, but not glyphosate, was associated with greater epigenetic age acceleration.

DISCUSSION

Glyphosate and AMPA exposure were associated with DNA methylation differences that could promote the development of cancer and other diseases. Further studies are warranted to replicate our results, determine the functional impact of glyphosate- and AMPA-associated differential DNA methylation, and further explore whether DNA methylation could serve as a biomarker of glyphosate exposure. https://doi.org/10.1289/EHP10174.

Meta-Analysis of Randomized Clinical Trials Evaluating Effectiveness of a Multivitamin Supplementation against Oxidative Stress in Healthy Subjects

A meta-analysis has been widely applied to draw general conclusions using a set of studies with similar purposes and designs. This study aimed to perform a meta-analysis of six randomized placebo-controlled trials, independently conducted for the relationship between a plant-based multivitamin/mineral supplementation (PMS) and oxidative stress for 6 to 8 weeks, to provide overall estimates of those effects. In detail, linear mixed model analysis was first conducted on each study to obtain individual estimates; then, two types of meta-analysis were applied to combine the individual estimates from all available studies (overall meta-analysis) and region-specific studies (subgroup meta-analysis). In the meta-analysis, we selected 19 biomarker variables that overlapped in at least two studies and found 6 variables significant in at least one meta-analysis. The overall estimates of beta coefficients were 0.17 for vitamin C, 0.80 for vitamin B₆, 0.46 for vitamin B₁₂, 0.81 for folate, 0.36 for β-carotene, and -0.17 for oxidized LDL (ox-LDL). Subsequent association analysis revealed significant negative correlations between plasma free radical scavenging nutrients and plasma ox-LDL levels, indicating a general benefit of PMS in alleviating oxidative stress by providing exogenous oxidant scavengers.

Identification of SSBP1 as a prognostic marker in human lung adenocarcinoma using bioinformatics approaches

OBJECTIVE

Single-stranded DNA-binding protein 1 (SSBP1) plays an important role in DNA repair processes and the maintenance of genomic stability. The aim of this study was to evaluate the expression of SSBP1 and its prognostic value in lung adenocarcinoma (LUAD) using bioinformatics approaches.

METHODS

We applied databases including UALCAN, Kaplan-Meier plotter, LinkedOmics, Webgestalt, cBioPortal and TIMER2.0 in this study.

RESULTS

We found that SSBP1 expression was up-regulated in LUAD samples and was correlated with clinicopathological features including age, cancer stage, and nodal metastasis status by the UALCAN analysis. Multivariate Cox regression analysis by the Kaplan-Meier plotter showed that high SSBP1 expression was independently correlated with poor overall survival (hazard ratio = 1.63, 95% confidence interval: 1.08-2.46, logrank P = 0.02). The LinkedOmics analysis showed that 5078 genes were positively correlated with SSBP1 expression, whereas 7905 genes were negatively correlated with SSBP1 in LUAD. Functional enrichment analysis using the Webgestalt tool showed that for SSBP1 and the genes positively correlating with it, the significantly enriched biological process was ribosomal large subunit biogenesis, and the significantly enriched pathway was proteasome. According to the cBioPortal database, the frequency of SSBP1 alterations was 1.7% in LUAD patients, and patients with SSBP1 alterations had worse prognosis (logrank P = 4.26e-05) compared with those unaltered for SSBP1. Finally, SSBP1 expression was negatively correlated with B cell infiltration level (Rho = -0.193, P = 1.54e-05) and the expression of B cell biomarkers including CD79A and CD19.

CONCLUSION

Our results suggest that SSBP1 may be a prognostic marker for human LUAD.

Impact of Alternative Splicing Variants on Liver Cancer Biology

Simple Summary

Among the top ten deadly solid tumors are the two most frequent liver cancers, hepatocellular carcinoma, and intrahepatic cholangiocarcinoma, whose development and malignancy are favored by multifactorial conditions, which include aberrant maturation of pre-mRNA due to abnormalities in either the machinery involved in the splicing, i.e., the spliceosome and associated factors, or the nucleotide sequences of essential sites for the exon recognition process. As a consequence of cancer-associated aberrant splicing in hepatocytes- and cholangiocytes-derived cancer cells, abnormal proteins are synthesized. They contribute to the dysregulated proliferation and eventually transformation of these cells to phenotypes with enhanced invasiveness, migration, and multidrug resistance, which contributes to the poor prognosis that characterizes these liver cancers.

Abstract

The two most frequent primary cancers affecting the liver, whose incidence is growing worldwide, are hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA), which are among the five most lethal solid tumors with meager 5-year survival rates. The common difficulty in most cases to reach an early diagnosis, the aggressive invasiveness of both tumors, and the lack of favorable response to pharmacotherapy, either classical chemotherapy or modern targeted therapy, account for the poor outcome of these patients. Alternative splicing (AS) during pre-mRNA maturation results in changes that might affect proteins involved in different aspects of cancer biology, such as cell cycle dysregulation, cytoskeleton disorganization, migration, and adhesion, which favors carcinogenesis, tumor promotion, and progression, allowing cancer cells to escape from pharmacological treatments. Reasons accounting for cancer-associated aberrant splicing include mutations that create or disrupt splicing sites or splicing enhancers or silencers, abnormal expression of splicing factors, and impaired signaling pathways affecting the activity of the splicing machinery. Here we have reviewed the available information regarding the impact of AS on liver carcinogenesis and the development of malignant characteristics of HCC and iCCA, whose understanding is required to develop novel therapeutical approaches aimed at manipulating the phenotype of cancer cells.

Comprehensive analysis of N6 -methyladenosine-related long non-coding RNAs for prognosis prediction in liver hepatocellular carcinoma

Background

Liver hepatocellular carcinoma (LIHC) is a lethal cancer. This study aimed to identify the N⁶‐methyladenosine (m⁶A)‐targeted long non‐coding RNA (lncRNA) related to LIHC prognosis and to develop an m⁶A‐targeted lncRNA model for prognosis prediction in LIHC.

Methods

The expression matrix of mRNA and lncRNA was obtained, and differentially expressed (DE) mRNAs and lncRNAs between tumor and normal samples were identified. Univariate Cox and pathway enrichment analyses were performed on the m⁶A‐targeted lncRNAs and the LIHC prognosis‐related m⁶A‐targeted lncRNAs. Prognostic analysis, immune infiltration, and gene DE analyses were performed on LIHC subgroups, which were obtained from unsupervised clustering analysis. Additionally, a multi‐factor Cox analysis was used to construct a prognostic risk model based on the lncRNAs from the LASSO Cox model. Univariate and multivariate Cox analyses were used to assess prognostic independence.

Results

A total of 5031 significant DEmRNAs and 292 significant DElncRNAs were screened, and 72 LIHC‐specific m⁶A‐targeted binding lncRNAs were screened. Moreover, a total of 29 LIHC prognosis‐related m⁶A‐targeted lncRNAs were obtained and enriched in cytoskeletal, spliceosome, and cell cycle pathways. An 11‐m⁶A‐lncRNA prognostic model was constructed and verified; the top 10 lncRNAs included LINC00152, RP6‐65G23.3, RP11‐620J15.3, RP11‐290F5.1, RP11‐147L13.13, RP11‐923I11.6, AC092171.4, KB‐1460A1.5, LINC00339, and RP11‐119D9.1. Additionally, the two LIHC subgroups, Cluster 1 and Cluster 2, showed significant differences in the immune microenvironment, m⁶A enzyme genes, and prognosis of LIHC.

Conclusion

The m⁶A‐lncRNA prognostic model accurately and effectively predicted the prognostic survival of LIHC. Immune cells, immune checkpoints (ICs), and m⁶A enzyme genes could act as novel therapeutic targets for LIHC.

Potential Role and Clinical Value of PPP2CA in Hepatocellular Carcinoma

BACKGROUND AND AIMS

Protein phosphatase 2A (PP2A) is associated with many cancers. This study aimed to clarify whether PPP2CA, which encodes the alpha isoform of the catalytic subunit of PP2A, plays a role in hepatocellular carcinoma (HCC) and to identify the potential underlying molecular pathways.

METHODS

Based on bioinformatics, public databases and our in-house RNA-Seq database, we analyzed the clinical value and molecular mechanism of PPP2CA in HCC.

RESULTS

Data were analyzed from 2,545 patients with HCC and 1,993 controls without HCC indexed in The Cancer Genome Atlas database, the Gene Expression Omnibus database and our in-house RNA-Seq database. PPP2CA expression was significantly higher in HCC tissue than in non-cancerous tissues (standardized mean difference: 0.69, 95% confidence interval [CI]: 0.50-0.89). PPP2CA expression was able to differentiate HCC from non-HCC, with an area under the summary receiver operator characteristic curve of 0.79 (95% CI: 0.75-0.83). Immunohistochemistry of tissue sections confirmed that PPP2CA protein was up-regulated in HCC tissues. High PPP2CA expression in HCC patients was associated with shorter overall, progression-free and disease-free survival. Potential molecular pathways through which PPP2CA may be involved in HCC were determined using miRWalk 2.0 as well as analysis of Gene Ontology categories, Kyoto Encyclopedia of Genes and Genomes pathways, and protein-protein interaction networks.

CONCLUSIONS

PPP2CA is up-regulated in HCC and higher expression correlates with worse prognosis. PPP2CA shows potential as a diagnostic marker for HCC. Future studies should examine whether PPP2CA contributes to HCC through the candidate microRNAs, pathways and hub genes identified in this study.

Highly expressed centromere protein L indicates adverse survival and associates with immune infiltration in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is characterized by rapid progression, high recurrence rate and poor prognosis. Early prediction for the prognosis and immunotherapy efficacy is of great significance to improve the survival of HCC patients. However, there is still no reliable predictor at present. This study is aimed to explore the role of centromere protein L (CENPL) in predicting prognosis and its association with immune infiltration in HCC.

METHODS

The expression of CENPL was identified through analyzing the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) data. The association between CENPL expression and clinicopathological features was investigated by the Wilcoxon signed-rank test or Kruskal Wallis test and logistic regression. The role of CENPL in prognosis was examined via Kaplan-Meier method and Log-rank test as well as univariate and multivariate Cox regression analysis. Besides, in TIMER and GEPIA database, we investigated the correlation between CENPL level and immunocyte and immunocyte markers, and the prognostic-related methylation sites in CENPL were identified by MethSurv.

RESULTS

CENPL had a high expression in HCC samples. Increased CENPL was prominently associated with unfavorable survival, and maybe an independent prognostic factor of worse overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), progression-free interval (PFI). Additionally, CENPL expression was significantly correlated with immune cell infiltration and some markers. CENPL also contained a methylation site that was notably related to poor prognosis.

CONCLUSIONS

Elevated CENPL may be a promising prognostic marker and associate with immune infiltration in HCC.

circRPS16 Promotes Proliferation and Invasion of Hepatocellular Carcinoma by Sponging miR-876-5p to Upregulate SPINK1

The roles of serine protease inhibitor Kazal type 1 (SPINK1) in multiple types of cancers have been significantly documented. However, its specific roles in hepatocellular carcinoma (HCC) remain to be investigated. This study found that SPINK1 is upregulated in HCC and its upregulation correlates with poor prognosis. Besides, functional assays revealed that SPINK1 promotes cell proliferation, cell cycle, and invasion in vitro. Through bioinformatics analysis, we speculate that circRPS16 regulates SPINK1 expression by sponging miR-876-5p. This was further verified by the dual-luciferase reporter and fluorescent in situ hybridization (FISH) assays. Subsequently, rescue assays verified that circRPS16 promotes cell proliferation, cell cycle, and invasion through miR-876-5p. Importantly, silencing circRPS16 inhibited tumor growth by downregulating SPINK1 expression in vivo. Collectively, our results confirm that SPINK1 is a downstream target of circRPS16. Besides, circRPS16 and SPINK1 are oncogenic factors in HCC progression; they provide novel diagnostic and therapeutic targets for HCC patients.

ZEB1 serves an oncogenic role in the tumourigenesis of HCC by promoting cell proliferation, migration, and inhibiting apoptosis via Wnt/β-catenin signaling pathway

Zinc finger E-box-binding homeobox 1 (ZEB1), a functional protein of zinc finger family, was aberrant expressed in many kinds of liver disease including hepatic fibrosis and Hepatitis C virus. Bioinformatics results showed that ZEB1 was abnormally expressed in HCC tissues. However, to date, the potential regulatory role and molecular mechanisms of ZEB1 are still unclear in the occurrence and development of HCC. This study demonstrated that the expression level of ZEB1 was significantly elevated both in liver tissues of HCC patients and cell lines (HepG2 and SMMC-7721 cells). Moreover, ZEB1 could promote the proliferation, migration, and invasion of HCC cells. On the downstream regulation mechanism, ZEB1 could activate the Wnt/β-catenin signaling pathway by upregulating the protein expression levels of β-catenin, c-Myc, and cyclin D1. Novel studies showed that miR-708 particularly targeted ZEB1 3'-UTR regions and inhibited the HCC cell proliferation, migration, and invasion. Furthermore, results of nude mice experiments of HCC model indicated that miR-708 could inhibit tumor growth and xenograft metastasis model was established to validate that miR-708 could inhibit HCC cell metastasis through tail-vein injection in vivo. Together, the study suggested that ZEB1 modulated by miR-708 might be a potential therapeutic target for HCC therapy.

Potential Role of NEU1 in Hepatocellular Carcinoma: A Study Based on Comprehensive Bioinformatical Analysis

Background: Aberrant expression of NEU1 has been identified in many malignancies. Nevertheless, the clinical significance of NEU1 in hepatocellular carcinoma (HCC) has not been fully elucidated.

Methods: In our study, multiple databases, including ONCOMINE, The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), International Cancer Genome Consortium (ICGC), Cancer Cell Line Encyclopedia (CCLE), Human Protein Atlas (HPA), Kaplan-Meier (KM) plotter, MethSurv, Gene Expression Profiling Interactive Analysis (GEPIA), and Metascape, etc., were utilized to investigate the expression, prognostic value, and function of NEU1 in HCC.

Results: ONCOMINE, GEO, and TCGA analyses revealed that NEU1 was more highly expressed in HCC compared to normal tissues. Additionally, the mRNA and protein expression levels of NEU1 were increased in liver cancer cell lines and HCC tissues, respectively. Moreover, a trend toward increased NEU1 expression with advanced stage or grade was found. Furthermore, higher mRNA expression of NEU1 was found to be remarkably correlated with worse survival in HCC patients, and multivariate Cox analysis indicated that high mRNA expression of NEU1 was an independent prognostic factor for poor prognosis of HCC patients. Also, 21 methylated CpGs were found to be significantly related to HCC prognosis. Besides, functional enrichment analyses indicated that high NEU1 expression group had lower levels of B cells, CD8⁺ T cells, NK cells, and T helper cells, etc. than the low NEU1 expression group, and NEU1 may regulate a variety of tumor-related proteins and pathways, including lysosome, spliceosome, mTOR signaling pathway and so on.

Conclusion: High expression level of NEU1 was positively correlated with unfavorable prognosis of HCC patients, which may be related to the regulation of cancer-associated pathways and the inhibition of immune function by NEU1. Thus, NEU1 could be used as a potential prognostic biomarker and target for HCC.

TRA2A Binds With LncRNA MALAT1 To Promote Esophageal Cancer Progression By Regulating EZH2/β-catenin Pathway

The RNA binding protein TRA2A, a member of the transformer 2 homolog family, plays a crucial role in the alternative splicing of pre-mRNA. However, it remains unclear whether TRA2A is involved in non-coding RNA regulation and, if so, what are the functional consequences. By analyzing expression profiling data, we found that TRA2A is highly expressed in esophageal cancer and is associated with disease-free survival and overall survival time. Subsequent gain- and loss-of-function studies demonstrated that TRA2A promotes proliferation and migration of esophageal squamous cell carcinoma and adenocarcinoma cells. RNA immunoprecipitation and RNA pull-down assay indicated that TRA2A can directly bind specific sites on MALAT1 in cells. In addition, ectopic expression or depletion of TRA2A leads to MALAT expression changes accordingly, thus modulates EZH2/β-catenin pathway. Together, these findings elucidated that TRA2A triggers carcinogenesis via MALAT1 mediated EZH2/β-catenin axis in esophageal cancer cells.

Quantitative proteomics identifies a plasma multi-protein model for detection of hepatocellular carcinoma

More efficient biomarkers are needed to facilitate the early detection of hepatocellular carcinoma (HCC). We aimed to identify candidate biomarkers for HCC detection by proteomic analysis. First, we performed a global proteomic analysis of 10 paired HCC and non-tumor tissues. Then, we validated the top-ranked proteins by targeted proteomic analyses in another tissue cohort. At last, we used enzyme-linked immunosorbent assays to validate the candidate biomarkers in multiple serum cohorts including HCC cases (HCCs), cirrhosis cases (LCs), and normal controls (NCs). We identified and validated 33 up-regulated proteins in HCC tissues. Among them, eight secretory or membrane proteins were further evaluated in serum, revealing that aldo-keto reductase family 1 member B10 (AKR1B10) and cathepsin A (CTSA) can distinguish HCCs from LCs and NCs. The area under the curves (AUCs) were 0.891 and 0.894 for AKR1B10 and CTSA, respectively, greater than that of alpha-fetoprotein (AFP; 0.831). Notably, combining the three proteins reached an AUC of 0.969, which outperformed AFP alone (P < 0.05). Furthermore, the serum AKR1B10 levels dramatically decreased after surgery. AKR1B10 and CTSA are potential serum biomarkers for HCC detection. The combination of AKR1B10, CTSA, and AFP may improve the HCC diagnostic efficacy.

The Prediction of Survival in Hepatocellular Carcinoma Based on A Four Long Non-coding RNAs Expression Signature

Prognostic stratification in hepatocellular carcinoma (HCC) patients is still challenging. Long non-coding RNAs (lncRNAs) have been proven to play a crucial role in tumorigenesis and progression of cancers. The aim of this study is to develop a useful prognostic index based on lncRNA signature to identify patients at high risk of disease progression. We obtained lncRNA expression profiles from three publicly available datasets from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). By the risk scoring method, we built an individualized four-lncRNA signature (HCCLnc-4) to predict survival of HCC patients in the discovery set (ROC curve, AUC: 0.83, 95% CI: 0.65-1.00, P < 0.05, Kaplan-Meier analysis and log-rank test, P < 0.01). Similar prognostic value of HCCLnc-4 has been further verified in two other independent sets. Stratified analysis and multivariate Cox regression analysis suggested the independence of HCCLnc-4 for prediction of HCC patient survival from traditional clinicopathological factors. Area under curve (AUC) analysis suggested that HCCLnc-4 could compete sufficiently with, or might be even better than classical pathological staging systems to predict HCC patient prognosis in the same data sets. Functional analysis and network analysis suggested the potential implication of lncRNA biomarkers. Our study developed and validated the lncRNA prognostic index of HCC patients, warranting further clinical evaluation and preventive interventions.

Development of an immune-related prognostic index associated with hepatocellular carcinoma

Liver hepatocellular carcinoma (LIHC), an inflammation-associated cancer induced by a variety of etiological factors, is still one of the most prevalent and lethal cancers in human population. In this study, the expression profiles of immune-related genes (IRGs) were integrated with the overall survival (OS) of 378 LIHC patients based on the Cancer Genome Atlas (TCGA) dataset. Moreover, the differentially expressed and survival related IRGs among LIHC patients were predicted through the computational difference algorithm and COX regression analysis. As a result, 7 genes, including HSPA4, S100A10, FABP6, CACYBP, HDAC1, FCGR2B and SHC1, were retrieved to construct a predictive model associated with the overall survival (OS) of LIHC patients. Typically, the as-constructed model performed moderately in predicting prognosis, which was also correlated with tumor grade. Functional enrichment analysis revealed that the genes of high-risk group were actively involved in mRNA binding and the spliceosome pathway. Intriguingly, the prognostic index established based on IRGs reflected infiltration by multiple types of immunocytes. Our findings screen several IRGs with clinical significance, reveal the drivers of immune repertoire, and illustrate the importance of a personalized, IRG-based immune signature in LIHC recognition, surveillance, and prognosis prediction.

The SF3b complex: splicing and beyond

The SF3b complex is an intrinsic component of the functional U2 small nuclear ribonucleoprotein (snRNP). As U2 snRNP enters nuclear pre-mRNA splicing, SF3b plays key roles in recognizing the branch point sequence (BPS) and facilitating spliceosome assembly and activation. Since the discovery of SF3b, substantial progress has been made in elucidating its molecular mechanism during splicing. In addition, numerous recent studies indicate that SF3b and its components are engaged in various molecular and cellular events that are beyond the canonical role in splicing. This review summarizes the current knowledge on the SF3b complex and highlights its multiple roles in splicing and beyond.

SF3b4: A Versatile Player in Eukaryotic Cells

Spliceosomes are large protein-RNA complexes regulating pre-mRNA processing in eukaryotes. SF3b4 encodes a core subunit of the U2-type spliceosome, loss- or gain-of-function of which often associates with abnormal cell growth, leading to tumorigenesis. Homologs of SF3b4 in other phyla are also essential. In this review, we summarize recent findings on the function of SF3b4. Importantly, we highlight the versatile roles of SF3b4, not only as a component for pre-mRNA splicing, but also as a regulator for transcription, translation, and cell signaling. Recent studies of SF3b4 homologs in different species across evolution will facilitate a better understanding of human diseases caused by the malfunction of SF3b4, such as Nager syndrome (NS) and cancer, in the future.

Development and internal validation of a novel model and markers to identify the candidates for lymph node metastasis in patients with prostate cancer

BACKGROUND

High-grade prostate cancer (PCa) has a poor prognosis, and up to 15% of patients worldwide experience lymph node invasion (LNI). To further improve the prediction lymph node invasion in prostate cancer, we adopted risk scores of the genes expression based on the nomogram in guidelines.

METHODS

We analyzed clinical data from 320 PCa patients from the Cancer Genome Atlas database. Weighted gene coexpression network analysis was used to identify the genes that were significantly associated with LNI in PCa (n = 390). Analyses using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases were performed to identify the activated signaling pathways. Univariate and multivariate logistic regression analyses were performed to identify the independent risk factors for the presence of LNI.

RESULTS

We found that patients with actual LNI and predicted LNI had the worst survival outcomes. The 7 most significant genes (CTNNAL1, ENSA, MAP6D1, MBD4, PRCC, SF3B2, TREML1) were selected for further analysis. Pathways in the cell cycle, DNA replication, oocyte meiosis, and 9 other pathways were dramatically activated during LNI in PCa. Multivariate analyses identified that the risk score (odds ratio [OR] = 1.05 for 1% increase, 95% confidence interval [CI]: 1.04-1.07, P < .001), serum PSA level, clinical stage, primary biopsy Gleason grade (OR = 2.52 for a grade increase, 95% CI: 1.27-5.22, P = .096), and secondary biopsy Gleason grade were independent predictors of LNI. A nomogram built using these predictive variables showed good calibration and a net clinical benefit, with an area under the curve (AUC) value of 90.2%.

CONCLUSIONS

In clinical practice, the application of our nomogram might contribute significantly to the selection of patients who are good candidates for surgery with extended pelvic lymph node dissection.

Gene Expression Pathways across Multiple Tissues in Antineutrophil Cytoplasmic Antibody-associated Vasculitis Reveal Core Pathways of Disease Pathology

OBJECTIVE

To identify commonalities in gene expression data across all antineutrophil cytoplasmic antibody (ANCA)-associated vasculitis (AAV) tissues thus far characterized.

METHODS

Gene expression data were collected from the 3 AAV tissues thus far characterized (orbit, peripheral leukocytes, and sinus brushings). These data were analyzed to identify commonly expressed genes and disease pathways. The pathways data were adjusted for multiple comparisons using a combined local false discovery rate, which estimates the probability of a false discovery of a given pathway in all 3 tissues analyzed.

RESULTS

Only 4 genes were upregulated in all 3 tissues - IL1RN, TLR2, SLC11A1, and MMP9. After multiple comparison adjustments, the network pathway analysis revealed 28 pathways associated with all 3 tissues. The most strongly associated pathway for all 3 tissues was the neutrophil degranulation pathway [multidimensional local false discovery (md-locfdr) = 1.05 × 10-12], followed by the osteoclast differentiation (md-locfdr = 3.8 × 10-05), cell surface interactions at the vascular wall (md-locfdr = 4.2 × 10-04), signaling by interleukins (md-locfdr = 6.1 × 10-04), and phagosome (md-locfdr = 0.003) pathways. There were no downregulated genes or pathways common to all 3 tissues.

CONCLUSION

This analysis identified individual genes and pathways of disease common to all AAV tissues thus far characterized. The use of a network pathway analysis allowed us to identify pathologic mechanisms that were not readily apparent in the commonly expressed genes alone. Many of these pathways are consistent with current theories about infectious drivers and the crossroads of innate and adaptive immune mechanisms. In addition, this analysis highlights novel pathways, such as vessel wall interactions and platelet activation, which require further investigation.

Meta-analysis of mammary RNA seq datasets reveals the molecular understanding of bovine lactation biology

A better understanding of the biology of lactation, both in terms of gene expression and the identification of candidate genes for the production of milk and its components, is made possible by recent advances in RNA seq technology. The purpose of this study was to understand the synthesis of milk components and the molecular pathways involved, as well as to identify candidate genes for milk production traits within whole mammary transcriptomic datasets. We performed a meta-analysis of publically available RNA seq transcriptome datasets of mammary tissue/milk somatic cells. In total, 11 562 genes were commonly identified from all RNA seq based mammary gland transcriptomes. Functional annotation of commonly expressed genes revealed the molecular processes that contribute to the synthesis of fats, proteins, and lactose in mammary secretory cells and the molecular pathways responsible for milk synthesis. In addition, we identified several candidate genes responsible for milk production traits and constructed a gene regulatory network for RNA seq data. In conclusion, this study provides a basic understanding of the lactation biology of cows at the gene expression level.

SF3B4 is regulated by microRNA-133b and promotes cell proliferation and metastasis in hepatocellular carcinoma

BACKGROUND

Splicing factor 3b subunit 4 (SF3B4) is a splicing factor and potential oncogene in hepatocellular carcinoma (HCC); however, its regulatory mechanism is yet unclear. We aimed to determine the role of SF3B4 in HCC and the underlying mechanism.

METHODS

To investigate the association between alternative splicing events and miRNAs, putative miRNAs were screened using TargetScan. Expression levels of and prognostic information for SF3B4 and miRNAs were determined based on public genomic data and clinical samples. Then, we examined the possible roles of SF3B4 and miRNA-133b in HCC cells and a xenograft mouse model. Pearson correlation analysis and in vitro experiments verified SF3B4 as a miRNA-133b target. Protein levels of key targets from the SF3B4 signaling pathway were estimated using western blotting.

FINDINGS

The expression of SF3B4 was upregulated in HCC tissues and cell lines whereas, the expression of miRNA-133b was downregulated. MiRNA-133b negatively regulated the expression of SF3B4. Effects of SF3B4 overexpression were partially abolished by miRNA-133b mimics, confirming that SF3B4 is a target of miRNA-133b. Moreover, molecules associated with SF3B4, including KLF4, KIP1, and SNAI2, were also modulated by miRNA-133b.

INTERPRETATION

SF3B4 plays a crucial role in HCC and is negatively regulated by miRNA-133b. The miRNA-133b/ SF3B4 axis may serve as a new therapeutic target for HCC treatment. FUND: China National Funds for Distinguished Young Scientists (No.81425019), the State Key Program of National Natural Science Foundation of China (No.81730076), Shanghai Science and Technology Committee Program (No.18XD1405300) and Specially-Appointed Professor Fund of Shanghai (GZ2015009). China National Funds for National Natural Science Fund (No.81672899).

A new metabolic gene signature in prostate cancer regulated by JMJD3 and EZH2

Histone methylation is essential for gene expression control. Trimethylated lysine 27 of histone 3 (H3K27me3) is controlled by the balance between the activities of JMJD3 demethylase and EZH2 methyltransferase. This epigenetic mark has been shown to be deregulated in prostate cancer, and evidence shows H3K27me3 enrichment on gene promoters in prostate cancer. To study the impact of this enrichment, a transcriptomic analysis with TaqMan Low Density Array (TLDA) of several genes was studied on prostate biopsies divided into three clinical grades: normal (n = 23) and two tumor groups that differed in their aggressiveness (Gleason score ≤ 7 (n = 20) and >7 (n = 19)). ANOVA demonstrated that expression of the gene set was upregulated in tumors and correlated with Gleason score, thus discriminating between the three clinical groups. Six genes involved in key cellular processes stood out: JMJD3, EZH2, MGMT, TRA2A, U2AF1 and RPS6KA2. Chromatin immunoprecipitation demonstrated collocation of EZH2 and JMJD3 on gene promoters that was dependent on disease stage. Gene set expression was also evaluated on prostate cancer cell lines (DU 145, PC-3 and LNCaP) treated with an inhibitor of JMJD3 (GSK-J4) or EZH2 (DZNeP) to study their involvement in gene regulation. Results showed a difference in GSK-J4 sensitivity under PTEN status of cell lines and an opposite gene expression profile according to androgen status of cells. In summary, our data describe the impacts of JMJD3 and EZH2 on a new gene signature involved in prostate cancer that may help identify diagnostic and therapeutic targets in prostate cancer.

Barrier to autointegration factor 1, procollagen-lysine, 2-oxoglutarate 5-dioxygenase 3, and splicing factor 3b subunit 4 as early-stage cancer decision markers and drivers of hepatocellular carcinoma

An accurate tool enabling early diagnosis of hepatocellular carcinoma (HCC) is clinically important, given that early detection of HCC markedly improves survival. We aimed to investigate the molecular markers underlying early progression of HCC that can be detected in precancerous lesions. We designed a gene selection strategy to identify potential driver genes by integrative analysis of transcriptome and clinicopathological data of human multistage HCC tissues, including precancerous lesions, low- and high-grade dysplastic nodules. The gene selection process was guided by detecting the selected molecules in both HCC and precancerous lesion. Using various computational approaches, we selected 10 gene elements as a candidate and, through immunohistochemical staining, showed that barrier to autointegration factor 1 (BANF1), procollagen-lysine, 2-oxoglutarate 5-dioxygenase 3 (PLOD3), and splicing factor 3b subunit 4 (SF3B4) are HCC decision markers with superior capability to diagnose early-stage HCC in a large cohort of HCC patients, as compared to the currently popular trio of HCC diagnostic markers: glypican 3, glutamine synthetase, and heat-shock protein 70. Targeted inactivation of BANF1, PLOD3, and SF3B4 inhibits in vitro and in vivo liver tumorigenesis by selectively modulating epithelial-mesenchymal transition and cell-cycle proteins. Treatment of nanoparticles containing small-interfering RNAs of the three genes suppressed liver tumor incidence as well as tumor growth rates in a spontaneous mouse HCC model. We also demonstrated that SF3B4 overexpression triggers SF3b complex to splice tumor suppressor KLF4 transcript to nonfunctional skipped exon transcripts. This contributes to malignant transformation and growth of hepatocyte through transcriptional inactivation of p27^Kip1 and simultaneously activation of Slug genes.

CONCLUSION

The findings suggest molecular markers of BANF1, PLOD3, and SF3B4 indicating early-stage HCC in precancerous lesion, and also suggest drivers for understanding the development of hepatocarcinogenesis. (Hepatology 2018;67:1360-1377).

Aberrantly activated Gli2-KIF20A axis is crucial for growth of hepatocellular carcinoma and predicts poor prognosis

Glioma-associated oncogene 2 (Gli2), a primary transcriptional regulator of Hedgehog (Hh) signaling, is essential for hepatocellular carcinoma (HCC) growth and survival. However, the underlying molecular mechanism and crucial downstream targets of Gli2 in human HCC are not fully understood. Here, we report the identification of kinesin family member 20A (KIF20A) as a novel downstream target of Gli2, which is important for HCC proliferation and tumor growth. Inhibition of Hh signaling leads to a remarkable decrease of KIF20A expression in HCC cells, whereas overexpression of Gli2 elevates KIF20A expression by activating Forkhead Box M1 (FoxM1)-MMB complex-mediated transcription of this kinesin gene. Gli2-induced HCC cell growth requires enhanced expression of KIF20A, and knockdown of Gli2 or KIF20A represses the proliferation of HCC cells in vitro and in vivo. Correlated with these results, analyses of clinical HCC samples show that Gli2, FoxM1 and KIF20A are highly elevated in primary HCC samples and represent significant risk factors for HCC recurrence and survival. Conclusion: KIF20A is an important downstream target gene of Hh signaling. And, the Gli2-KIF20A axis is essential for the proliferation and growth of human HCC cells. Our study also suggests Gli2-KIF20A axis as a potential target for future therapeutic intervention and as an independent prognostic biomarker for HCC.

DEIsoM: a hierarchical Bayesian model for identifying differentially expressed isoforms using biological replicates

Motivation

High-throughput mRNA sequencing (RNA-Seq) is a powerful tool for quantifying gene expression. Identification of transcript isoforms that are differentially expressed in different conditions, such as in patients and healthy subjects, can provide insights into the molecular basis of diseases. Current transcript quantification approaches, however, do not take advantage of the shared information in the biological replicates, potentially decreasing sensitivity and accuracy.

Results

We present a novel hierarchical Bayesian model called Differentially Expressed Isoform detection from Multiple biological replicates (DEIsoM) for identifying differentially expressed (DE) isoforms from multiple biological replicates representing two conditions, e.g. multiple samples from healthy and diseased subjects. DEIsoM first estimates isoform expression within each condition by (1) capturing common patterns from sample replicates while allowing individual differences, and (2) modeling the uncertainty introduced by ambiguous read mapping in each replicate. Specifically, we introduce a Dirichlet prior distribution to capture the common expression pattern of replicates from the same condition, and treat the isoform expression of individual replicates as samples from this distribution. Ambiguous read mapping is modeled as a multinomial distribution, and ambiguous reads are assigned to the most probable isoform in each replicate. Additionally, DEIsoM couples an efficient variational inference and a post-analysis method to improve the accuracy and speed of identification of DE isoforms over alternative methods. Application of DEIsoM to an hepatocellular carcinoma (HCC) dataset identifies biologically relevant DE isoforms. The relevance of these genes/isoforms to HCC are supported by principal component analysis (PCA), read coverage visualization, and the biological literature.

Availability and implementation

The software is available at https://github.com/hao-peng/DEIsoM

Supplementary information

Supplementary data are available at Bioinformatics online.

Serine peptidase inhibitor Kazal type 1 (SPINK1) as novel downstream effector of the cadherin-17/β-catenin axis in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is the most common type of liver cancer worldwide. Previously, we reported that cadherin-17 (CDH17) and its related CDH17/β-catenin axis may be responsible for inducing HCC in a subset of patients exhibiting CDH17 over-expression. Here we aimed at obtaining a better understanding of the CDH17-related HCC biology and to obtain further indications for the design of targeted therapies in CDH17 over-expressing HCC patients.

RESULTS

We found that SPINK1 acts as a downstream effector of the CDH17/β-catenin axis in HCC. In addition, we found that SPINK1 expression exhibited a positive correlation with CDH17 expression in human HCCs and was over-expressed in up to 70% of the tumors. We identified SPINK1 as a downstream effector of the CDH17/β-catenin axis using a spectrum of in vitro assays, including gene expression modulation and inhibitor assays, bioinformatics analyses and luciferase reporter assays. These in vitro results were validated in primary human HCCs, including the observation that alteration in β-catenin expression (a core component of the CDH17/β-catenin axis) in tumors affects SPINK1 serum levels in HCC patients. Similar to CDH17, SPINK1 expression in HCC cells was found to be associated with specific tumor-related properties via activating the c-Raf/MEK/ERK pathway.

CONCLUSIONS

Our current data substantiate our knowledge on the role of CDH17 in the biology of HCC and suggest that components of the CDH17/β-catenin axis may serve as therapeutic targets in CDH17 over-expressing HCC patients.

Metformin alters DNA methylation genome-wide via the H19/SAHH axis

The molecular mechanisms underlying the antineoplastic properties of metformin, a first-line drug for type 2 diabetes, remain elusive. Here we report that metformin induces genome-wide alterations in DNA methylation by modulating the activity of S-adenosylhomocysteine hydrolase (SAHH). Exposing cancer cells to metformin leads to hypermethylation of tumor-promoting pathway genes and concomitant inhibition of cell proliferation. Metformin acts by upregulating microRNA let-7 through AMPK activation, leading to degradation of H19 long noncoding RNA, which normally binds to and inactivates SAHH. H19 knockdown activates SAHH, enabling DNA methyltransferase 3B to methylate a subset of genes. This metformin-induced H19 repression and alteration of gene methylation are recapitulated in endometrial cancer tissue samples obtained from patients treated with antidiabetic doses of metformin. Our findings unveil a novel mechanism of action for the drug metformin with implications for the molecular basis of epigenetic dysregulation in cancer. This novel mechanism of action also may be occurring in normal cells.

SF3B4 is decreased in pancreatic cancer and inhibits the growth and migration of cancer cells

Splicing factor 3b subunit 4, a critical component of pre-message RNA splicing complex, has been reported to play an important part in the tumorigenesis. However, the expression pattern and biological role of splicing factor 3b subunit 4 in pancreatic cancer have never been investigated. In this study, we found that both the messenger RNA ( p < 0.001) and protein level of splicing factor 3b subunit 4 were decreased significantly in pancreatic cancer specimens compared with their adjacent normal tissues. Overexpression of splicing factor 3b subunit 4 in pancreatic cancer cells inhibited cell growth and motility in vitro, while suppressing splicing factor 3b subunit 4 expression promoted the proliferation and migration of pancreatic cancer cells. In addition, splicing factor 3b subunit 4 was found to inhibit the activity of signal transducer and activator of transcription 3 signaling via downregulating the phosphorylation of signal transducer and activator of transcription 3 on a tyrosine residue at position 705. Taken together, these findings demonstrated that splicing factor 3b subunit 4 acted as a suppressive role in pancreatic cancer and indicated that restoring the function of splicing factor 3b subunit 4 might be a strategy for cancer therapy.

Pediatric pineal germinomas: Epigenetic and genomic approach

OBJECTIVE

We identify and correlate chromosomal alterations, methylation patterns and gene expression in pediatric pineal germinomas.

METHODS

CGH microarray, methylation and gene expression were performed through the Agilent platform. The results were analyzed with MatLab software, MapViewer, DAVID, GeneCards and Hippie.

RESULTS

Amplifications were found in 1q24.2, 1q31.3, 2p11.2, 3p22.2, 7p13, 7p15.2, 8p22, 12p13.2, 14q24.3 y 22q12; and deletions were found in 1q21.2, 9p24.1, 10q11.22, 11q11, 15q11.2 and 17q21.31. In the methylation analysis, we observed 10,428 CpG Islands with a modified methylation status that may affect 11,726 genes. We identified 1260 overexpressed genes and 470 underexpressed genes. The genes RUNDC3A, CDC247, CDCA7L, ASAH1, TRA2A, LPL and NPC2 were altered among the three levels.

CONCLUSIONS

We identified the 1q24.2 and 1q31.3 amplified regions and the 1q21.3 and 11q11 deleted regions as the most important aims. The genes NPC2 and ASAH1 may play an important role in the development, progression and tumor maintenance. The ASAH1 gene is an ideal candidate to identify drug responses. These genomic and epigenetic studies may help to characterize the formation of pineal germ cell tumors to determine prognostic markers and also to identify shared characteristics in gonadal and extragonadal tumors.

Ethoxyquin provides neuroprotection against cisplatin-induced neurotoxicity

Ethoxyquin was recently identified as a neuroprotective compound against toxic neuropathies and efficacy was demonstrated against paclitaxel-induced neurotoxicity in vivo. In this study we examined the efficacy of ethoxyquin in preventing neurotoxicity of cisplatin in rodent models of chemotherapy-induced peripheral neuropathy and explored its mechanism of action. Ethoxyquin prevented neurotoxicity of cisplatin in vitro in a sensory neuronal cell line and primary rat dorsal root ganglion neurons. In vivo, chronic co-administration of ethoxyquin partially abrogated cisplatin-induced behavioral, electrophysiological and morphological abnormalities. Furthermore, ethoxyquin did not interfere with cisplatin's ability to induce tumor cell death in ovarian cancer cell line in vitro and in vivo. Finally, ethoxyquin reduced the levels of two client proteins (SF3B2 and ataxin-2) of a chaperone protein, heat shock protein 90 (Hsp90) when co-administered with cisplatin in vitro. These results implied that the neuroprotective effect of ethoxyquin is mediated through these two client proteins of Hsp90. In fact, reducing levels of SF3B2 in tissue-cultured neurons was effective against neurotoxicity of cisplatin. These findings suggest that ethoxyquin or other compounds that inhibit chaperone activity of Hsp90 and reduce levels of its client protein, SF3B2 may be developed as an adjuvant therapy to prevent neurotoxicity in cisplatin-based chemotherapy protocols.

Gene expression and pathway analysis of human hepatocellular carcinoma cells treated with cadmium

Cadmium (Cd) is a toxic and carcinogenic metal naturally occurring in the Earth's crust. A common route of human exposure is via diet and cadmium accumulates in the liver. The effects of Cd exposure on gene expression in human hepatocellular carcinoma (HepG2) cells were examined in this study. HepG2 cells were acutely-treated with 0.1, 0.5, or 1.0 μM Cd for 24h; or chronically-treated with 0.01, 0.05, or 0.1 μM Cd for three weeks and gene expression analysis was performed using Affymetrix GeneChip® Human Gene 1.0 ST Arrays. Acute and chronic exposures significantly altered the expression of 333 and 181 genes, respectively. The genes most upregulated by acute exposure included several metallothioneins. Downregulated genes included the monooxygenase CYP3A7, involved in drug and lipid metabolism. In contrast, CYP3A7 was upregulated by chronic Cd exposure, as was DNAJB9, an anti-apoptotic J protein. Genes downregulated following chronic exposure included the transcriptional regulator early growth response protein 1. Ingenuity Pathway Analysis revealed that the top networks altered by acute exposure were lipid metabolism, small molecule biosynthesis, cell morphology, organization, and development; while top networks altered by chronic exposure were organ morphology, cell cycle, cell signaling, and renal and urological diseases/cancer. Many of the dysregulated genes play important roles in cellular growth, proliferation, and apoptosis, and may be involved in carcinogenesis. In addition to gene expression changes, HepG2 cells treated with cadmium for 24h indicated a reduction in global levels of histone methylation and acetylation that persisted 72 h post-treatment.

Identification of Commonly Dysregulated Genes in Non-small-cell Lung Cancer by Integrated Analysis of Microarray Data and qRT-PCR Validation

BACKGROUND

Non-small-cell lung cancer (NSCLC), the most common lung cancer, leads to the largest number of cancer-related deaths worldwide. There are many studies to identify the differentially expressed genes (DEGs) between NSCLC and normal control (NC) tissues by means of microarray technology. Because of the inconsistency of the microarray data sets, we performed an integrated analysis to identify DEGs and analyzed their biological function.

METHODS AND RESULTS

We combined 15 microarray data sets and identified 1063 DEGs between NSCLC and NC tissues; in addition, we found that the DEGs were enriched in regulation of cell proliferation process and focal adhesion signaling pathway. The protein-protein interaction network analysis for the top 20 significantly DEGs revealed that CAV1, COL1A1, and ADRB2 were the significant hub proteins. Finally, we employed qRT-PCR to validate the meta-analysis approach by determining the expression of the top 10 most significantly DEGs and found that the expression of these genes were significantly different between tumor and NC tissues, in accordance with the results of meta-analysis.

CONCLUSION

qRT-PCR results indicated that the meta-analysis approach in our study was acceptable. Our data suggested that some of the DEGs, including MMP12, COL11A1, THBS2, FAP, and CAV1, may participate in the pathology of NSCLC and could be applied as potential markers or therapeutic targets for NSCLC.