Comprehensive characterization of the rRNA metabolism-related genes in human cancer

The impact of ribosome biogenesis in cancer: from proliferation to metastasis

The dysregulation of ribosome biogenesis is a hallmark of cancer, facilitating the adaptation to altered translational demands essential for various aspects of tumor progression. This review explores the intricate interplay between ribosome biogenesis and cancer development, highlighting dynamic regulation orchestrated by key oncogenic signaling pathways. Recent studies reveal the multifaceted roles of ribosomes, extending beyond protein factories to include regulatory functions in mRNA translation. Dysregulated ribosome biogenesis not only hampers precise control of global protein production and proliferation but also influences processes such as the maintenance of stem cell-like properties and epithelial-mesenchymal transition, contributing to cancer progression. Interference with ribosome biogenesis, notably through RNA Pol I inhibition, elicits a stress response marked by nucleolar integrity loss, and subsequent G1-cell cycle arrest or cell death. These findings suggest that cancer cells may rely on heightened RNA Pol I transcription, rendering ribosomal RNA synthesis a potential therapeutic vulnerability. The review further explores targeting ribosome biogenesis vulnerabilities as a promising strategy to disrupt global ribosome production, presenting therapeutic opportunities for cancer treatment.

Genomic hallmarks and therapeutic targets of ribosome biogenesis in cancer

Hyperactive ribosome biogenesis (RiboSis) fuels unrestricted cell proliferation, whereas genomic hallmarks and therapeutic targets of RiboSis in cancers remain elusive, and efficient approaches to quantify RiboSis activity are still limited. Here, we have established an in silico approach to conveniently score RiboSis activity based on individual transcriptome data. By employing this novel approach and RNA-seq data of 14 645 samples from TCGA/GTEx dataset and 917 294 single-cell expression profiles across 13 cancer types, we observed the elevated activity of RiboSis in malignant cells of various human cancers, and high risk of severe outcomes in patients with high RiboSis activity. Our mining of pan-cancer multi-omics data characterized numerous molecular alterations of RiboSis, and unveiled the predominant somatic alteration in RiboSis genes was copy number variation. A total of 128 RiboSis genes, including EXOSC4, BOP1, RPLP0P6 and UTP23, were identified as potential therapeutic targets. Interestingly, we observed that the activity of RiboSis was associated with TP53 mutations, and hyperactive RiboSis was associated with poor outcomes in lung cancer patients without TP53 mutations, highlighting the importance of considering TP53 mutations during therapy by impairing RiboSis. Moreover, we predicted 23 compounds, including methotrexate and CX-5461, associated with the expression signature of RiboSis genes. The current study generates a comprehensive blueprint of molecular alterations in RiboSis genes across cancers, which provides a valuable resource for RiboSis-based anti-tumor therapy.

A pan-cancer analysis of Dyskeratosis congenita 1 (DKC1) as a prognostic biomarker

BACKGROUND

Dyskeratosis congenita 1 (DKC1), a critical component of telomerase complex, is highly expressed in a variety of human cancers. However, the association of DKC1 with cancer occurrence and development stages is not clear, making a pan-cancer analysis crucial.

METHODS

We conducted a study using various bioinformatic databases such as TIMER, GEPIA, UALCAN, and KM plotter Analysis to examine the different expressions of DKC1 in multiple tissues and its correlation with pathological stages. Through KEGG analysis, GO enrichment analysis and Venn analysis, we were able to reveal DKC1-associated genes and signaling pathways. In addition, we performed several tests including the CCK, wound healing assay, cell cycle arrest assay, transwell assay and Sa-β-gal staining on DKC1-deleted MDA-231 cells.

RESULTS

Our study demonstrates that DKC1 has relatively low expression specificity in different tissues. Furthermore, we found that in ACC, KICH, KIRP and LIHC, the expression level of DKC1 is positively correlated with pathological stages. Conversely, in NHSC, KIRP, LGG, LIHC, MESO and SARC, we observed a negative influence of DKC1 expression level on the overall survival rate. We also found a significant positive correlation between DKC1 expression and Tumor Mutational Burden in 14 tumors. Additionally, we observed a significantly negative impact of DKC1 DNA methylation on gene expression at the promoter region in BRCA. We also identified numerous phosphorylation sites concentrated at the C-terminus of the DKC1 protein. Our GO analysis revealed a correlation between DKC1 and ribosomal biosynthesis pathways, and the common element UTP14A was identified. We also observed decreased rates of cell proliferation, migration and invasion abilities in DKC1-knockout MDA-MB-231 cell lines. Furthermore, DKC1-knockout induced cell cycle arrest and caused cell senescence.

CONCLUSIONS

Our findings suggest that the precise expression of DKC1 is closely associated with the occurrence and developmental stages of cancer in multiple tissues. Depletion of DKC1 can inhibit the abilities of cancer cells to proliferate, migrate, and invade by arresting the cell cycle and inducing cell senescence. Therefore, DKC1 may be a valuable prognostic biomarker for the diagnosis and treatment of cancer in various tissues.

RTEL1 is upregulated in colorectal cancer and promotes tumor progression

As one of the most common cancers worldwide, the incidence of colorectal cancer (CRC) continues to increase. Metastasis is the leading cause of death for this malignant disease. Regulator of telomere elongation helicase 1 (RTEL1) is a key factor that maintains telomere stability and contributes to DNA repair. Recent advances have suggested that RTEL1 exerts other functions through various mechanisms. However, little is known about its role in human cancers, including CRC. In this study, we revealed that the copy number of RTEL1 was significantly higher in CRC tissues, especially in metastatic CRC tissues, than in paired normal tissues, which was associated with increased expression. Increased RTEL1 expression was significantly correlated with CRC progression and poor survival. Functionally, RTEL1 knockdown suppressed the proliferation and metastasis of CRC both in vitro and in vivo. In addition, multiple signalling pathways, including the mitotic cell cycle, DNA repair, and EMT, were potentially regulated by RTEL1. Notably, GPR17 appeared to be a candidate downstream target gene that partially mediated the tumor-promoting effects of RTEL1 in CRC. Altogether, our results indicate that RTEL1 plays a crucial role in CRC progression and appears to be a promising therapeutic target and prognosis for CRC.

Palmitoylation landscapes across human cancers reveal a role of palmitoylation in tumorigenesis

BACKGROUND

Protein palmitoylation, which is catalyzed by palmitoyl-transferase and de-palmitoyl-transferase, plays a crucial role in various biological processes. However, the landscape and dynamics of protein palmitoylation in human cancers are not well understood.

METHODS

We utilized 23 palmitoyl-acyltransferases and seven de-palmitoyl-acyltransferases as palmitoylation-related genes for protein palmitoylation analysis. Multiple publicly available datasets were employed to conduct pan-cancer analysis, examining the transcriptome, genomic alterations, clinical outcomes, and correlation with c-Myc (Myc) for palmitoylation-related genes. Real-time quantitative PCR and immunoblotting were performed to assess the expression of palmitoylation-related genes and global protein palmitoylation levels in cancer cells treated with Myc depletion or small molecule inhibitors. Protein docking and drug sensitivity analyses were employed to predict small molecules that target palmitoylation-related genes.

RESULTS

We identified associations between palmitoylation and cancer subtype, stage, and patient survival. We discovered that abnormal DNA methylation and oncogenic Myc-driven transcriptional regulation synergistically contribute to the dysregulation of palmitoylation-related genes. This dysregulation of palmitoylation was closely correlated with immune infiltration in the tumor microenvironment and the response to immunotherapy. Importantly, dysregulated palmitoylation was found to modulate canonical cancer-related pathways, thus influencing tumorigenesis. To support our findings, we performed a proof-of-concept experiment showing that depletion of Myc led to reduced expression of most palmitoylation-related genes, resulting in decreased global protein palmitoylation levels. Through mass spectrometry and enrichment analyses, we also identified palmitoyl-acyltransferases ZDHHC7 and ZDHHC23 as significant contributors to mTOR signaling, DNA repair, and immune pathways, highlighting their potential roles in tumorigenesis. Additionally, our study explored the potential of three small molecular (BI-2531, etoposide, and piperlongumine) to modulate palmitoylation by targeting the expression or activity of palmitoylation-related genes or enzymes.

CONCLUSIONS

Overall, our findings underscore the critical role of dysregulated palmitoylation in tumorigenesis and the response to immunotherapy, mediated through classical cancer-related pathways and immune cell infiltration. Additionally, we propose that the aforementioned three small molecule hold promise as potential therapeutics for modulating palmitoylation, thereby offering novel avenues for cancer therapy.

ELOA promotes tumor growth and metastasis by activating RBP1 in gastric cancer

BACKGROUND

Elongin A (ELOA), our previous work revealed, serves as a novel tumor suppressor in colorectal cancer. However, the function and mechanism of ELOA in other cancer types, including gastric cancer (GC), remain to be elucidated.

METHODS

The expression of ELOA was measured by quantitative reverse transcription-polymerase chain reaction and western blot. The effects of ELOA on GC growth and metastasis were assessed through a series of in-vitro and in-vivo assays. Furthermore, the potential mechanism of ELOA was revealed by RNA sequencing, dual luciferase reporter assay, chromatin immunoprecipitation, and rescue experiments in GC.

RESULTS

We uncovered increased expression of ELOA in GC tissues compared with paired normal tissues via bioinformatic analyses and our sample detection. Enhanced ELOA expression in GC tissues was obviously correlated with poor tumor differentiation, lymph node metastasis, advanced tumor stage, and a poor prognosis. A series of functional experiments showed that ELOA promoted the proliferation and metastasis of GC. Mechanistically, we revealed that the decreased levels of miR-490-3p caused the upregulation of ELOA in GC. Both RNA-seq and ChIP assays revealed that ELOA transcriptionally activated retinol-binding protein 1 (RBP1) by binding to its promotor. Furthermore, specific knockdown of RBP1 reduced the tumor-promoting ability of ELOA in GC cells.

CONCLUSIONS

In summary, our findings demonstrate that ELOA exerts oncogenic properties by activating RBP1 expression, providing the basis for a promising therapeutic target in GC.

Infiltration of a Unique CD8+CD274+ Cell Subgroup in Hepatocellular Carcinoma is Associated with Poor Clinical Outcomes

Introduction

Immune checkpoint (IC) inhibitor-related immunotherapies have attracted considerable attention in hepatocellular carcinoma (HCC). High IC expression and high tumor infiltrating lymphocyte levels are the current indicators of sensitivity to IC inhibitors. Thus, it is imperative to apply precision medicine strategies for patient selection.

Methods

Six independent HCC cohorts were used for analysis at the single-cell and tissue levels. Multiplex immunofluorescence and immunochemistry staining assays were used to validate our results. A series of methodologies were used for immune-related evaluations.

Results

Herein, we uncovered a unique CD8⁺CD274⁺ cell subpopulation that is associated with tumor progression and poor survival in HCC at the single-cell level. We assessed this subset at the tissue level and found that the prognostic significance of CD274 is dependent on CD8A expression in HCC. Subsequently, we identified a unique high-risk subpopulation that showed high CD8A expression coupled with intense CD274 expression in multiple HCC cohorts. CD8A^HighCD274^High* subgroup was correlated with malignant indexes and remained an independent prognostic factor when considering the influence of these indexes. Molecular characteristic analyses showed that the CD8A^HighCD274^High* subgroup harbored more mutations, had higher immune response activity and presented enrichment of cancer-related biological processes. Moreover, this high-risk subpopulation in HCC was characterized by high immune cell infiltration, low tumor purity, and enrichment of cancer-related signatures. Finally, cases with this phenotype demonstrated higher immunomodulator and IC levels and greater sensitivity to IC inhibitors.

Conclusion

Our findings illustrate that some HCC patients may have a poor prognosis despite high CD8⁺ T-cell infiltration. These patients would probably benefit from IC inhibitor-based combination treatment.

Development of an endoplasmic reticulum stress-related signature with potential implications in prognosis and immunotherapy in head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is a multisite malignancy that responds well to immunotherapy. Despite the initial enthusiasm, the clinical benefits of immunotherapy in HNSCC patients are overall limited. Endoplasmic reticulum stress (ERS) has been indicated to play a key role in the process of anti-tumor immune response mediation. However, ERS-related biomarkers which can accurately predict prognosis and immunotherapy response in HNSCC are still lacking.

METHODS AND RESULTS

In this study, we identify and validate an ERS-related signature comprises of six genes (ASNS, EXOSC6, BAK1, TPP1, EXOSC8, and TATDN2) that can predict the prognosis of HNSCC patients. GSEA analysis indicates that the ERS-related signature is significantly correlated with tumor immunity in HNSCC. Moreover, the infiltration of naive B cells and CD8 + T cells are significantly diminished in patients with high-risk scores compared to those with low-risk scores, while macrophages and activated mast cells are remarkably enhanced. Furthermore, the ERS-related signature also displays a tremendous potential for predicting immunotherapy response in HNSCC.

CONCLUSIONS

Our study identifies an ERS-related signature that can predict the prognosis of HNSCC patients and highlights its potential value as a predictive biomarker of immunotherapy response, potentially enabling more precise and personalized immunotherapy response and paving the way for further investigation of the prognostic and therapeutic potentials of ERS.

Integration of single sample and population analysis for understanding immune evasion mechanisms of lung cancer

A deep understanding of the complex interaction mechanism between the various cellular components in tumor microenvironment (TME) of lung adenocarcinoma (LUAD) is a prerequisite for understanding its drug resistance, recurrence, and metastasis. In this study, we proposed two complementary computational frameworks for integrating multi-source and multi-omics data, namely ImmuCycReg framework (single sample level) and L0Reg framework (population or subtype level), to carry out difference analysis between the normal population and different LUAD subtypes. Then, we aimed to identify the possible immune escape pathways adopted by patients with different LUAD subtypes, resulting in immune deficiency which may occur at different stages of the immune cycle. More importantly, combining the research results of the single sample level and population level can improve the credibility of the regulatory network analysis results. In addition, we also established a prognostic scoring model based on the risk factors identified by Lasso-Cox method to predict survival of LUAD patients. The experimental results showed that our frameworks could reliably identify transcription factor (TF) regulating immune-related genes and could analyze the dominant immune escape pathways adopted by each LUAD subtype or even a single sample. Note that the proposed computational framework may be also applicable to the immune escape mechanism analysis of pan-cancer.

Integrative analysis of WDR12 as a potential prognostic and immunological biomarker in multiple human tumors

Background: Mammalian WD-repeat protein 12 (WDR12), a family member of proteins containing repeats of tryptophan-aspartic acid (WD), is a potential homolog of yeast Ytm1p and consists of seven repeats of WD. Aim of the study: This study aims to investigate the potential oncogenic effects of WDR12 in various human malignancies throughout a pan-cancer analysis that has been carried out to examine the various patterns in which this gene is expressed and behaves in tumor tissues. Methods: Herein, we used The Cancer Genome Atlas (TCGA) and various computational tools to explore expression profiles, prognostic relevance, genetic mutations, immune cell infiltration, as well as the functional characteristics of WDR12 in multiple human cancers. Results: We found that WDR12 was inconsistently expressed in various cancers and that variations in WDR12 expression predicted survival consequences for cancer patients. Furthermore, we observed a significant correlation between WDR12 gene mutation levels and the prognosis of some tumors. Furthermore, significant correlations were found between WDR12 expression patterns and cancer-associated fibroblast (CAF) infiltration, myeloid-derived suppressor cells (MDSCs), tumor mutation burden, microsatellite instability and immunoregulators. Ultimately, pathway enrichment analysis revealed that WDR12-related pathways are involved in carcinogenesis. Conclusions: The findings of our study are stisfactory, demonstrating that WDR12 could serve as a promising reliable prognostic biomarker, as well as a therapeutic target for novel cancer therapeutic approaches.

c-Myc-IMPDH1/2 axis promotes tumourigenesis by regulating GTP metabolic reprogramming

BACKGROUND

Metabolic reprogramming is a hallmark of cancer. Metabolic rate-limiting enzymes and oncogenic c-Myc (Myc) play critical roles in metabolic reprogramming to affect tumourigenesis. However, a systematic assessment of metabolic rate-limiting enzymes and their relationship with Myc in human cancers is lacking.

METHODS

Multiple Pan-cancer datasets were used to develop the transcriptome, genomic alterations, clinical outcomes and Myc correlation landscapes of 168 metabolic rate-limiting enzymes across 20 cancers. Real-time quantitative PCR and immunoblotting were, respectively, used to examine the mRNA and protein of inosine monophosphate dehydrogenase 1 (IMPDH1) in human colorectal cancer (CRC), azoxymethane/dextran sulphate sodium-induced mouse CRC and spontaneous intestinal tumours from APCMin/+ mice. Clone formation, CCK-8 and subcutaneous xenograft model were applied to investigate the possible mechanisms connecting IMPDH1 to CRC growth. Co-immunoprecipitation and protein half-life assay were used to explore the mechanisms underlying the regulation of IMPDH1.

RESULTS

We explored the global expression patterns, dysregulation profiles, genomic alterations and clinical relevance of 168 metabolic rate-limiting enzymes across human cancers. Importantly, a series of enzymes were associated with Myc, especially top three upregulated enzymes (TK1, RRM2 and IMPDH1) were positively correlated with Myc in multiple cancers. As a proof-of-concept exemplification, we demonstrated that IMPDH1, a rate-limiting enzyme in GTP biosynthesis, is highly upregulated in CRC and promotes CRC growth in vitro and in vivo. Mechanistically, IMPDH2 stabilizes IMPDH1 by decreasing the polyubiquitination levels of IMPDH1, and Myc promotes the de novo GTP biosynthesis by the transcriptional activation of IMPDH1/2. Finally, we confirmed that the Myc-IMPDH1/2 axis is dysregulated across human cancers.

CONCLUSIONS

Our study highlights the essential roles of metabolic rate-limiting enzymes in tumourigenesis and their crosstalk with Myc, and the Myc-IMPDH1/2 axis promotes tumourigenesis by altering GTP metabolic reprogramming. Our results propose the inhibition of IMPDH1 as a viable option for cancer treatment.

EXOSC8 promotes colorectal cancer tumorigenesis via regulating ribosome biogenesis-related processes

Extensive protein synthesis is necessary for uncontrolled cancer cell proliferation, requiring hyperactive ribosome biogenesis. Our previous Pan-cancer study has identified EXOSC8 as a potential copy number variation (CNV)-driven rRNA metabolism-related oncogene in colorectal cancer (CRC). Herein, we further investigated proliferation-prompting functions and mechanisms of EXOSC8 in CRC by performing in silico analyses and wet-lab experiments. We uncovered that increased EXOSC8 expression and CNV levels are strongly associated with ribosome biogenesis-related factor levels in CRC, including ribosome proteins (RPs), eukaryotic translation initiation factors and RNA polymerase I/III. EXOSC8 silence decreases nucleolar protein and proliferation marker levels, as well as rRNA/DNA and global protein syntheses. Clinically, EXOSC8 is upregulated across human cancers, particularly CNV-driven upregulation in CRC was markedly associated with poor clinical outcomes. Mechanistically, EXOSC8 knockdown increased p53 levels in CRC, and the oncogenic proliferation phenotypes of EXOSC8 depended on p53 in vitro and in vivo. We discovered that EXOSC8 knockdown in CRC cells triggers ribosomal stress, nucleolar RPL5/11 being released into the nucleoplasm and "hijacking" Mdm2 to block its E3 ubiquitin ligase function, thus releasing and activating p53. Furthermore, our therapeutic experiments provided initial evidence that EXOSC8 might serve as a potential therapeutic target in CRC. Our findings revealed, for the first time, that the RNA exosome gene (EXOSC8) promotes CRC tumorigenesis by regulating cancer-related ribosome biogenesis in CRC. This study further extends our previous Pan-cancer study of the rRNA metabolism-related genes. The inhibition of EXOSC8 is a novel therapeutic strategy for the RPs-Mdm2-p53 ribosome biogenesis surveillance pathway in CRC.

A novel high-risk subpopulation identified by CTSL and ZBTB7B in gastric cancer

BACKGROUND

Gastric cancer (GC) is characterised by a heterogeneous tumour microenvironment (TME) that is closely associated with the response to treatment, especially immunotherapies. However, most previous GC molecular subtyping systems need complex gene signatures and examination methods, restricting their clinical applications. Thus, we developed a new TME-based molecular subtype using only two genes.

METHODS

Nine independent GC cohorts at the tissue- or single-cell level with more than 2000 patients were used in this study, including data we examined by single-cell sequencing, quantitative RT-PCR and immunochemistry/immunofluorescence staining. Nine different methods, five existing molecular subtypes and a series of signatures were used to evaluate the TME and molecular characteristics of GC.

RESULTS

We established a CTSL/ZBTB7B subtyping system and uncovered the novel CTSLHighZBTB7BLow high-risk subgroup, but characterised by relative higher immune cell infiltration and lower tumour purity. This subgroup demonstrate higher levels of immune checkpoints and more enrichment of cancer-related pathways compared with other cases.

CONCLUSIONS

We identified a high-risk subpopulation with unique TME features based on expressions of CTSL and ZBTB7B, suggesting a counterbalancing phenotype between immunostimulatory and immunosuppressive mechanisms. This subtyping system could be used to select treatment and management strategies for GC.

Identification and validation of a novel tumor driver gene signature for diagnosis and prognosis of head and neck squamous cell carcinoma

Objective: Head and neck squamous cell carcinoma (HNSCC) is a common heterogeneous cancer with complex carcinogenic factors. However, the current TNM staging criteria to judge its severity to formulate treatment plans and evaluate the prognosis are particularly weak. Therefore, a robust diagnostic model capable of accurately diagnosing and predicting HNSCC should be established. Methods: Gene expression and clinical data were retrieved from The Cancer Genome Atlas and Gene Expression Omnibus databases. Key prognostic genes associated with HNSCC were screened with the weighted gene co-expression network analysis and least absolute shrinkage and selection operator (LASSO) Cox regression model analysis. We used the timeROC and survival R packages to conduct time-dependent receiver operating characteristic curve analyses and calculated the area under the curve at different time points of model prediction. Patients in the training and validation groups were divided into high- and low-risk subgroups, and Kaplan-Meier (K-M) survival curves were plotted for all subgroups. Subsequently, LASSO and support vector machine algorithms were used to screen genes to construct diagnostic model. Furthermore, we used the Wilcoxon signed-rank test to compare the half-maximal inhibitory concentrations of common chemotherapy drugs among patients in different risk groups. Finally, the expression levels of eight genes were measured using quantitative real-time polymerase chain reaction and immunohistochemistry. Results: Ten genes (SSB, PFKP, NAT10, PCDH9, SHANK2, PAX8, CELSR3, DCLRE1C, MAP2K7, and ODF4) with prognostic potential were identified, and a risk score was derived accordingly. Patients were divided into high- and low-risk groups based on the median risk score. The K-M survival curves confirmed that patients with high scores had significantly worse overall survival. Receiver operating characteristic curves proved that the prognostic signature had good sensitivity and specificity for predicting the prognosis of patients with HNSCC. Univariate and multivariate Cox regression analyses confirmed that the gene signature was an independent prognostic risk factor for HNSCC. Diagnostic model was built by identifying eight genes (SSB, PFKP, NAT10, PCDH9, CELSR3, DCLRE1C, MAP2K7, and ODF4). The high-risk group showed higher sensitivity to various common chemotherapeutic drugs. DCLRE1C expression was higher in normal tissues than in HNSCC tissues. Conclusion: Our study identified the important role of tumor-driver genes in HNSCC and their potential clinical diagnostic and prognostic values to facilitate individualized management of patients with HNSCC.

Pan-cancer analyses of classical protein tyrosine phosphatases and phosphatase-targeted therapy in cancer

Protein tyrosine phosphatases function in dephosphorylating target proteins to regulate signaling pathways that control a broad spectrum of fundamental physiological and pathological processes. Detailed knowledge concerning the roles of classical PTPs in human cancer merits in-depth investigation. We comprehensively analyzed the regulatory mechanisms and clinical relevance of classical PTPs in more than 9000 tumor patients across 33 types of cancer. The independent datasets and functional experiments were employed to validate our findings. We exhibited the extensive dysregulation of classical PTPs and constructed the gene regulatory network in human cancer. Moreover, we characterized the correlation of classical PTPs with both drug-resistant and drug-sensitive responses to anti-cancer drugs. To evaluate the PTP activity in cancer prognosis, we generated a PTPscore based on the expression and hazard ratio of classical PTPs. Our study highlights the notable role of classical PTPs in cancer biology and provides novel intelligence to improve potential therapeutic strategies based on pTyr regulation.

Integrated Bioinformatic Investigation of EXOSCs in Hepatocellular Carcinoma Followed by the Preliminary Validation of EXOSC5 in Cell Proliferation

The Exosome complex (EXOSC) is a multiprotein complex that was originally discovered as the machinery of RNA degradation. Interestingly, recent studies have reported that EXOSC family members (EXOSCs) are associated with various human diseases, including cancers. It will be interesting to investigate whether EXOSCs are related to the processes of hepatocellular carcinoma (HCC). In this study, multiple public databases and experimental validation were utilized to systemically investigate the role of EXOSCs, especially EXOSC5, in HCC. It is worth considering that the mRNA and protein levels of many EXOSCs were elevated in HCC, although there were some differences in the results from different database analyses. The over-expression of EXOSCs could predict HCC to some extent, as evidenced by the positive correlation between the elevated EXOSCs and alpha fetoprotein (AFP) levels, as well as with a high accuracy, as shown by the receiver operating characteristic curve analysis. Additionally, higher mRNA expressions of specific EXOSCs were significantly related to clinical cancer stage, shorter overall survival and disease-free survival in HCC patients. A moderate mutation rate of EXOSCs was also observed in HCC. Furthermore, a gene functional enrichment analysis indicated that EXOSCs were mainly involved in the metabolism of RNA. Moreover, we revealed that the expression of EXOSCs is remarkably related to immune cell infiltration. Finally, EXOSC5 was upregulated in HCC tissues and cell lines, promoting cell growth and proliferation via activated signal transducer and activator of transcription 3 (STAT3). The bioinformatic analyses, following verification in situ and in vitro, provided a direction for further functions and underlying mechanism of EXOSCs in HCC.

A pan-cancer analysis confirms PTPN11's potential as a prognostic and immunological biomarker

Protein tyrosine phosphatase, non-receptor type 11 (PTPN11) is a multifunctional tyrosine phosphatase and has a significant part in many types of tumors. As of yet, neither the expression profile of PTPN11 nor its significance in pan-cancer diagnosis has been clarified. With the assistance of The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), we have comprehensively mapped the expression profiles, prognostic significance, genetic alteration, phosphorylation status, infiltration of immune cells, and functional properties of PTPN11 in 33 human tumors. There was an inconsistent expression of PTPN11 in different tumors, and the alteration of PTPN11 expression predicted the survival outcomes of cancer patients. A significant association was found between the genetic alteration levels of PTPN11 and some tumor predictions. Besides, the reduced PTPN11 phosphorylation levels were observed in breast cancer, clear cell RCC, head and neck carcinoma, and lung adenocarcinoma (LUAD). Furthermore, there was a significant association between PTPN11 expression and infiltration of cancer-associated fibroblasts and endothelial cells, along with tumor mutation burden, microsatellite instability, mismatch repair genes, and immunoregulators. Finally, pathway enrichment analysis demonstrated that PTPN11-associated terms and pathways were involved in malignancy. Taken together, PTPN11 may become a new biomarker and target for cancer therapy.

Transcriptomic and Functional Evidence for Differential Effects of MCF-7 Breast Cancer Cell-Secretome on Vascular and Lymphatic Endothelial Cell Growth

Vascular and lymphatic vessels drive breast cancer (BC) growth and metastasis. We assessed the cell growth (proliferation, migration, and capillary formation), gene-, and protein-expression profiles of Vascular Endothelial Cells (VECs) and Lymphatic Endothelial Cells (LECs) exposed to a conditioned medium (CM) from estrogen receptor-positive BC cells (MCF-7) in the presence or absence of Estradiol. We demonstrated that MCF-7-CM stimulated growth and capillary formation in VECs but inhibited LEC growth. Consistently, MCF-7-CM induced ERK1/2 and Akt phosphorylation in VECs and inhibited them in LECs. Gene expression analysis revealed that the LECs were overall (≈10-fold) more sensitive to MCF-7-CM exposure than VECs. Growth/angiogenesis and cell cycle pathways were upregulated in VECs but downregulated in LECs. An angiogenesis proteome array confirmed the upregulation of 23 pro-angiogenesis proteins in VECs. In LECs, the expression of genes related to ATP synthesis and the ATP content were reduced by MCF-7-CM, whereas MTHFD2 gene, involved in folate metabolism and immune evasion, was upregulated. The contrasting effect of MCF-7-CM on the growth of VECs and LECs was reversed by inhibiting the TGF-β signaling pathway. The effect of MCF-7-CM on VEC growth was also reversed by inhibiting the VEGF signaling pathway. In conclusion, BC secretome may facilitate cancer cell survival and tumor growth by simultaneously promoting vascular angiogenesis and inhibiting lymphatic growth. The differential effects of BC secretome on LECs and VECs may be of pathophysiological relevance in BC.

Ferroptosis-Associated Molecular Features to Aid Patient Clinical Prognosis and Therapy Across Human Cancers

Ferroptosis is a new non-apoptotic form that regulates cell death and is mainly dependent on iron-mediated oxidative damage and subsequent cell membrane damage. Ferroptosis may be a potential therapeutic strategy for immunotherapy, chemotherapy, and radiotherapy in human cancers. Numerous studies have analyzed ferroptosis-correlated signatures or genes, but a systematic landscape of associations among tumor ferroptosis, clinical outcomes, tumor microenvironment, and therapies in human cancers is lacking. Here, we developed a relative ferroptosis level (RFL) combined with drive/suppress regulators and validated it in the Gene Expression Omnibus datasets of ferroptotic drug treatment. Based on this effective evaluation method, we classified about 7,000 tumor samples into high and low RFL groups in each cancer type and observed that high RFL cases demonstrate favorable survival outcomes in nine cancer types from The Cancer Genome Atlas. Then, several RFL-correlated candidate genes that have not been reported to be ferroptosis-related were selected and experimentally validated in five cancer cell lines using Erastin treatment. We further showed that both immunostimulatory and immunosuppressive phenotypes were observed in high RFL tumors, suggesting that the consideration of ferroptosis could be a potential strategy in cancer immunotherapy. Moreover, we found that high RFL cases/cells showed responder or sensitivity to chemotherapy and radiotherapy. Our study provides a comprehensive molecular-level understanding of ferroptosis and may have practical implications for clinical cancer therapies, including immunotherapy, chemotherapy, and radiotherapy.

Nuclear RNA Exosome and Pervasive Transcription: Dual Sculptors of Genome Function

The genome is pervasively transcribed across various species, yielding numerous non-coding RNAs. As a counterbalance for pervasive transcription, various organisms have a nuclear RNA exosome complex, whose structure is well conserved between yeast and mammalian cells. The RNA exosome not only regulates the processing of stable RNA species, such as rRNAs, tRNAs, small nucleolar RNAs, and small nuclear RNAs, but also plays a central role in RNA surveillance by degrading many unstable RNAs and misprocessed pre-mRNAs. In addition, associated cofactors of RNA exosome direct the exosome to distinct classes of RNA substrates, suggesting divergent and/or multi-layer control of RNA quality in the cell. While the RNA exosome is essential for cell viability and influences various cellular processes, mutations and alterations in the RNA exosome components are linked to the collection of rare diseases and various diseases including cancer, respectively. The present review summarizes the relationships between pervasive transcription and RNA exosome, including evolutionary crosstalk, mechanisms of RNA exosome-mediated RNA surveillance, and physiopathological effects of perturbation of RNA exosome.

FAM65A as a novel prognostic biomarker in human tumors reveal by a pan-cancer analysis

BACKGROUND

Family with sequence similarity 65 member A (FAM65A), also known as RIPOR1, is differentially expressed between human tumor and non-tumor tissues in kinds of cancers. In addition, it was reported that the product of FAM65A may be a biomarker for cholangiocarcinoma patients. However, there is still no evidence on the relationship between the FAM65A and different types of tumors. Our study is mainly for exploring the prognostic values of FAM65A in pan-cancer and for further discovering a potential therapeutics target.

METHODS

We analyzed FAM65A expression, prognostic values, genetic alteration, protein phosphorylation, immune infiltration and enrichment analysis across different types of human malignant tumors based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets. Additionally, Real-Time PCR (RT-qPCR) was performed to further confirm the roles of FAM65A in the pathogenesis of colorectal cancer.

RESULTS

We found that FAM65A expression was associated with the prognosis of multiple human tumors, especially colorectal cancer. Moreover, we also observed that FAM65A was highly expressed in colorectal cancer through RT-qPCR. We observed that decreasing phosphorylation level of the S351 locus in colon adenocarcinoma, uterine corpus endometrial carcinoma and lung adenocarcinoma. And the expression of FAM65A was positively related to cancer-associated fibroblasts (CAFs) infiltration in many tumors, such as colon adenocarcinoma. Therefore, FAM65A may be a potential prognostic biomarker of human tumors.

The membrane-associated E3 ubiquitin ligase MARCH3 downregulates the IL-6 receptor and suppresses colitis-associated carcinogenesis

The IL-6-STAT3 axis is critically involved in inflammation-associated carcinogenesis (IAC). How this axis is regulated to modulate IAC remains unknown. Here, we show that the plasma membrane-associated E3 ubiquitin ligase MARCH3 negatively regulates STAT3 activation triggered by IL-6, as well as another IL-6 subfamily member, Oncostatin M (OSM). MARCH3 is associated with the IL-6 receptor α-chain (IL-6Rα) and its coreceptor gp130. Biochemical experiments indicated that MARCH3 mediates the polyubiquitination of IL-6Rα at K401 and gp130 at K849 following IL-6 stimulation, leading to their translocation to and degradation in lysosomes. MARCH3 deficiency increases IL-6- and OSM-triggered activation of STAT3 and induction of downstream effector genes in various cell types. MARCH3 deficiency enhances dextran sulfate sodium (DSS)-induced STAT3 activation, increases the expression of inflammatory cytokines, and exacerbates colitis, as well as azoxymethane (AOM)/DSS-induced colitis-associated cancer in mice. In addition, MARCH3 is downregulated in human colorectal cancer tissues and associated with poor survival across different cancer types. Our findings suggest that MARCH3 is a pivotal negative regulator of IL-6-induced STAT3 activation, inflammation, and inflammation-associated carcinogenesis.

Role of FRG1 in predicting the overall survivability in cancers using multivariate based optimal model

FRG1 has a role in tumorigenesis and angiogenesis. Our preliminary analysis showed that FRG1 mRNA expression is associated with overall survival (OS) in certain cancers, but the effect varies. In cervix and gastric cancers, we found a clear difference in the OS between the low and high FRG1 mRNA expression groups, but the difference was not prominent in breast, lung, and liver cancers. We hypothesized that FRG1 expression level could affect the functionality of the correlated genes or vice versa, which might mask the effect of a single gene on the OS analysis in cancer patients. We used the multivariate Cox regression, risk score, and Kaplan Meier analyses to determine OS in a multigene model. STRING, Cytoscape, HIPPIE, Gene Ontology, and DAVID (KEGG) were used to deduce FRG1 associated pathways. In breast, lung, and liver cancers, we found a distinct difference in the OS between the low and high FRG1 mRNA expression groups in the multigene model, suggesting an independent role of FRG1 in survival. Risk scores were calculated based upon regression coefficients in the multigene model. Low and high-risk score groups showed a significant difference in the FRG1 mRNA expression level and OS. HPF1, RPL34, and EXOSC9 were the most common genes present in FRG1 associated pathways across the cancer types. Validation of the effect of FRG1 mRNA expression level on these genes by qRT-PCR supports that FRG1 might be an upstream regulator of their expression. These genes may have multiple regulators, which also affect their expression, leading to the masking effect in the survival analysis. In conclusion, our study highlights the role of FRG1 in the survivability of cancer patients in tissue-specific manner and the use of multigene models in prognosis.

Nanobody-Engineered Natural Killer Cell Conjugates for Solid Tumor Adoptive Immunotherapy

Cancer immunotherapy based on natural killer (NK) cells is demonstrated to be a promising strategy. However, NK cells are deficient in ligands that target specific tumors, resulting in limited antitumor efficacy. Here, a glycoengineering approach to imitate the chimeric antigen receptor strategy and decorate NK cells with nanobodies to promote NK-based immunotherapy in solid tumors is proposed. Nanobody 7D12, which specifically recognizes the human epidermal growth factor receptor (EGFR) that is overexpressed on many solid tumors, is coupled to the chemically synthesized DBCO-PEG₄ -GGG-NH₂ by sortase A-mediated ligation to generate DBCO-7D12. The NK92MI cells bearing azide groups are then equipped with DBCO-7D12 via bioorthogonal click chemistry. The resultant 7D12-NK92MI cells exhibit high specificity and affinity for EGFR-overexpressing tumor cells in vitro and in vivo by the 7D12-EGFR interaction, causing increased cytokine secretion to more effectively kill EGFR-positive tumor cells, but not EGFR-negative cancer cells. Importantly, the 7D12-NK92MI cells also show a wide anticancer spectrum and extensive tumor penetration. Furthermore, mouse experiments reveal that 7D12-NK92MI treatment achieves excellent therapeutic efficacy and outstanding safety. The authors' works provide a cell modification strategy using specific protein ligands without genetic manipulation and present a potential novel method for cancer-targeted immunotherapy by NK cells.

LncRNAs specifically overexpressed in endocervical adenocarcinoma are associated with an unfavorable recurrence prognosis and the immune response

Background

Cervical cancer is the fourth most common gynecological tumor in terms of both the incidence and mortality of females worldwide. Cervical squamous cell carcinoma (CSCC) accounts for 70–80% of cervical cancers, and endocervical adenocarcinoma (EAC) accounts for 20–25%. Unlike CSCC, EAC has worse clinical outcomes and prognosis. In this study, we explored the relationship between various types of long noncoding RNAs (lncRNAs) and pathological types of cervical cancer.

Methods

RNA sequencing (RNA-Seq) and clinical data from The Cancer Genome Atlas (TCGA) were used in this study. A single-sample gene set enrichment analysis (ssGSEA) and the ESTIMATE package were used to assess lncRNA activity and immune responses, respectively. RT-qPCR was performed to verify our findings.

Results

We explored the relationship between various types of lncRNAs and pathological types of cervical cancer. A series of long intergenic noncoding RNAs (lincRNAs) and antisense RNAs, which are the major types of lncRNAs, were identified to be specifically expressed in EAC and associated with a poor recurrence prognosis in patients with cervical cancer, suggesting that they might serve as independent prognostic markers of recurrence in patients with cervical cancer. RT-qPCR was performed to verify the 10 EAC-specific lncRNAs in cervical cancer samples we collected. Furthermore, the overexpression of these lncRNAs was positively correlated with EAC pathology levels but negatively correlated with immune responses in the microenvironment of cervical cancer.

Conclusions

These lncRNAs potentially represent new biomarkers for the prediction of the recurrence prognosis and help obtain deeper insights into potential immunotherapeutic approaches for treating cervical cancer.

Identification of prognostic long intergenic non-coding RNAs as competing endogenous RNAs with KRAS mutations in colorectal cancer

Colorectal cancer (CRC) is recognized as a common type of human cancer, and KRAS mutations are correlated with poor CRC survival outcomes. The evaluation and prediction of CRC results remain challenging. In the present study, RNA sequencing and clinical data from The Cancer Genome Atlas were used to identify KRAS mutation-related prognostic long intergenic non-coding RNAs (lincRNAs) in CRC. Significantly dysregulated lincRNAs and independent prognostic lincRNAs with KRAS mutations in CRC were identified. Two lincRNAs with KRAS mutations, LINC00265 and AL390719.2, were selected as key prognostic lincRNAs for both 10- and 5-year survival rates. In addition, competing endogenous (ce)RNA models were constructed to comprehensively assess the oncogenic performance of the two key lincRNAs. The ceRNA models suggested that LINC00265 and AL390719.2 are critical for the cell cycle and cancer pathways. Finally, reverse transcription-quantitative PCR was used to validate the ceRNA models in 12 pairs of CRC tissue samples. These prognostic lincRNAs may provide novel biomarkers for the prognostic prediction of CRC. The ceRNA model may also demonstrate the underlying mechanism of these lincRNAs in CRC.

Effect of ISM1 on the Immune Microenvironment and Epithelial-Mesenchymal Transition in Colorectal Cancer

Background: An increasing number of studies have shown that Isthmin 1 (ISM1), a secreted protein, is important in tumorigenesis and invasion, including in colorectal cancer (CRC). However, the mechanisms are still unclear. This study aims to explore the function and prognosis capacity of ISM1 in CRC.

Methods: We investigated the expression of ISM1 in 18 CRC tissues vs. adjacent normal tissues from GSE50760, 473 CRC tissues vs. 41 normal tissues from The Cancer Genome Atlas (TCGA), and across gastrointestinal cancer types. Differences were further confirmed in CRC tissues via quantitative real-time polymerase chain reaction (qRT-PCR). Then, we analyzed correlations between clinicopathologic features and ISM1 expression, including prognostic prediction value, using the Kaplan–Meier method and multivariate Cox regression. Gene set enrichment analysis (GSEA) was performed to identify ISM1-related pathways. In vitro experiments were performed to verify the role of ISM1 in epithelial-mesenchymal transition (EMT) and CRC progression.

Results: Multiple datasets showed that ISM1 is upregulated in CRC tissues, which was validated. Patients with higher ISM1 expression had shorter overall survival (OS), and ISM1 expression served as an independent prognostic factor. Enrichment analysis showed that ISM1 upregulation was positively correlated with cancer-related pathways, such as EMT, hypoxia, and the Notch and KRAS signaling pathways. We were exclusively interested in the connection between ISM1 and EMT because 71% of genes in this pathway were significantly positively co-expressed with ISM1, which may account for why patients with higher ISM1 expression are prone to regional lymph node involvement and progression to advanced stages. In addition, we found that ISM1 was positively correlated with multiple immunosuppressive pathways such as IL2/STAT5, TNF-α/NF-κB, and TGF-β, and immune checkpoints, including PD-L1, PD-1, CTLA-4, and LAG3, which may account for upregulation of ISM1 in immunotherapy-resistant patients. Notably, through in vitro experiments, we found that ISM1 promoted EMT and colon cancer cell migration and proliferation.

Conclusion: ISM1 is critical for CRC development and progression, which enhances our understanding of the low response rate of CRC to immunotherapy via immunosuppressive signaling pathways.

FBXL6 degrades phosphorylated p53 to promote tumor growth

The ubiquitin-proteasome system regulates many distinct biological processes. Its dysregulation causes various diseases, including but not limited to cancer. In this study, based on the analysis of gene expression in several colorectal cancer (CRC) datasets, we show that FBXL6, a poorly-characterized F-box protein, is amplified, over-expressed, and highly correlated with poor prognosis in human CRC patients. Mechanistically, FBXL6 targets phospho-p53 (S315) to mediate its polyubiquitination and proteasomal degradation, thereby inhibiting p53 signaling. FBXL6 depletion inhibits proliferation of p53 wild-type (WT) CRC cells by inducing cell cycle arrest and apoptosis. Furthermore, p53 transcriptionally suppresses FBXL6 expression by binding its core promoter region. Taken together, these results identify the feed-forward loop of FBXL6-p53 as a potential therapeutic target for CRC treatments.

Pan-cancer analysis reveals an immunological role and prognostic potential of PXN in human cancer

Paxillin (PXN) is a protein involved in numerous physiological processes, and its presence is closely related to the occurrence and development of many types of tumors. However, no studies have analyzed PXN from a pan-cancer perspective. We analyzed PXN expression, immune cell infiltration, prognosis, and biological function across different types of tumors included in The Cancer Genome Atlas and Gene Expression Omnibus datasets. The results showed that expression of PXN varies in different tumors. Expression of PXN strongly correlated with prognosis in patients with tumors; higher PXN expression usually was linked to poor overall and disease-free survival. Expression of PXN in breast invasive carcinoma and lymphoid neoplasm diffuse large B-cell lymphoma was related to the degree of CD8+ T-cell infiltration, and infiltration of cancer-associated fibroblasts, such as kidney renal papillary cell carcinoma and brain lower-grade glioma, was also observed in other tumors. The results of pan-cancer analysis showed that abnormal PXN expression was related to poor prognosis, immune infiltration, and protein phosphorylation in different tumor types. Therefore, the PXN gene may become a potential biomarker of clinical tumor prognosis.

Comprehensive Characterization of Integrin Subunit Genes in Human Cancers

Although integrin subunit genes (ITGs) have been reported to be associated with some human cancer types, a systematic assessment of ITGs across human cancers is lacking. Hence, we performed comprehensive analyses to investigate mRNA expression, copy number variation (CNV), DNA methylation, mutation, and clinical landscapes of ITGs in more than 8000 cancer patients from The Cancer Genome Atlas (TCGA) dataset. Landscapes of ITGs were established across 20 human cancer types. We observed that ITGs are extensively dysregulated with heterogeneity in different system cancer types, part of which are driven by CNV, DNA hypomethylation or mutation. Furthermore, dysregulated prognosis-related ITGs were systematically identified in each cancer type, including ITGA11 in stomach adenocarcinoma (STAD). The models based on dysregulated ITGs with clinical relevance and TNM staging indexes are good indicators in STAD and head and neck squamous cell carcinoma. Finally, ITGA11 is overexpressed and associated with poor survival in STAD cases from the TCGA and additionally Gene Expression Omnibus cohorts. Functionally, ITGA11 knockdown inhibits malignant phenotypes in STAD cell lines AGS and MKN45, demonstrating the oncogenic role of ITGA11 in STAD. Together, this study highlights the important roles of ITGs in tumorigenesis as potential prognostic biomarkers, and provide an effective resource that identifies cancer-related genes of ITGs in human cancers.

Comprehensive molecular and clinical characterization of SLC1A5 in human cancers

Although SLC1A5 has been reported to be closely associated with some cancer types, a comprehensive and systematic assessment of SLC1A5 across human cancers is lacking. Thus, Pan-cancer analysis of SLC1A5 was performed across 30 types of human cancers in this study. We examined mRNA expression, protein expression, copy number variation (CNV), DNA methylation, clinical relevance, cell functions, drug response and total immune infiltrates of SLC1A5 in more than 9000 patients across 30 human cancer types from The Cancer Genome Atlas (TCGA) dataset. Additionally, nine independent Gene Expression Omnibus datasets, more than 800 cancer cell lines from the Cancer Cell Line Encyclopedia dataset and the Project Achilles dataset were used to validate our findings in the TCGA dataset. Landscapes of SLC1A5 were established across multiple cancers. We showed that SLC1A5 is upregulated in multiple cancers, particularly in digestive and respiratory system cancers. SLC1A5 upregulation may be driven by CNV gain and DNA hypomethylation in human cancers. Furthermore, SLC1A5 overexpression is associated with tumor progression and poor survival in multiple cancers. Moreover, we systematically explored the potential effects of SLC1A5 expression on cell functions and drug response in human cancers. SLC1A5 knockdown showed significant proliferation-inhibiting effects in most human cancer types, especially in the digestive system and KRAS-mutant cancers. SLC1A5 expression is associated with proliferation activities of KRAS-mutant cancer cell lines and drug response of many anti-cancer drugs. Finally, we demonstrated that SLC1A5-realted tumor immune microenvironment characteristics showed strong heterogeneity in human cancers. Taken together, our findings highlight the important roles of SLC1A5 in tumorigenesis, progression, prognosis and therapy.

Identification and Validation of a Novel RNA-Binding Protein-Related Gene-Based Prognostic Model for Multiple Myeloma

Background

Multiple myeloma (MM) is a malignant hematopoietic disease that is usually incurable. RNA-binding proteins (RBPs) are involved in the development of many tumors, but their prognostic significance has not been systematically described in MM. Here, we developed a prognostic signature based on eight RBP-related genes to distinguish MM cohorts with different prognoses.

Method

After screening the differentially expressed RBPs, univariate Cox regression was performed to evaluate the prognostic relevance of each gene using The Cancer Genome Atlas (TCGA)-Multiple Myeloma Research Foundation (MMRF) dataset. Lasso and stepwise Cox regressions were used to establish a risk prediction model through the training set, and they were validated in three Gene Expression Omnibus (GEO) datasets. We developed a signature based on eight RBP-related genes, which could classify MM patients into high- and low-score groups. The predictive ability was evaluated using bioinformatics methods. Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and gene set enrichment analyses were performed to identify potentially significant biological processes (BPs) in MM.

Result

The prognostic signature performed well in the TCGA-MMRF dataset. The signature includes eight hub genes: HNRNPC, RPLP2, SNRPB, EXOSC8, RARS2, MRPS31, ZC3H6, and DROSHA. Kaplan-Meier survival curves showed that the prognosis of the risk status showed significant differences. A nomogram was constructed with age; B2M, LDH, and ALB levels; and risk status as prognostic parameters. Receiver operating characteristic (ROC) curve, C-index, calibration analysis, and decision curve analysis (DCA) showed that the risk module and nomogram performed well in 1, 3, 5, and 7-year overall survival (OS). Functional analysis suggested that the spliceosome pathway may be a major pathway by which RBPs are involved in myeloma development. Moreover, our signature can improve on the R-International Staging System (ISS)/ISS scoring system (especially for stage II), which may have guiding significance for the future.

Conclusion

We constructed and verified the 8-RBP signature, which can effectively predict the prognosis of myeloma patients, and suggested that RBPs are promising biomarkers for MM.

Multi-omics characterization and validation of invasiveness-related molecular features across multiple cancer types

Background

Tumor invasiveness reflects many biological changes associated with tumorigenesis, progression, metastasis, and drug resistance. Therefore, we performed a systematic assessment of invasiveness-related molecular features across multiple human cancers.

Materials and methods

Multi-omics data, including gene expression, miRNA, DNA methylation, and somatic mutation, in approximately 10,000 patients across 30 cancer types from The Cancer Genome Atlas, Gene Expression Omnibus, PRECOG, and our institution were enrolled in this study.

Results

Based on a robust gene signature, we established an invasiveness score and found that the score was significantly associated with worse prognosis in almost all cancers. Then, we identified common invasiveness-associated dysregulated molecular features between high- and low-invasiveness score group across multiple cancers, as well as investigated their mutual interfering relationships thus determining whether the dysregulation of invasiveness-related genes was caused by abnormal promoter methylation or miRNA expression. We also analyzed the correlations between the drug sensitivity data from cancer cell lines and the expression level of 685 invasiveness-related genes differentially expressed in at least ten cancer types. An integrated analysis of the correlations among invasiveness-related genetic features and drug response were conducted in esophageal carcinoma patients to outline the complicated regulatory mechanism of tumor invasiveness status in multiple dimensions. Moreover, functional enrichment suggests the invasiveness score might serve as a predictive biomarker for cancer patients receiving immunotherapy.

Conclusion

Our pan-cancer study provides a comprehensive atlas of tumor invasiveness and may guide more precise therapeutic strategies for tumor patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-02773-x

Identification of an immune overdrive high-risk subpopulation with aberrant expression of FOXP3 and CTLA4 in colorectal cancer

Colorectal cancer (CRC) is characterized by a heterogeneous tumor microenvironment (TME) that regulates cancer progression and therapeutic response. Overexpression of FOXP3 and CTLA4 is associated with immunosuppressive TME and poor prognosis in many cancer types. However, opposite results were reported in CRC. Thus, we performed comprehensive analyses to evaluate the exact prognostic value of FOXP3 and CTLA4 in CRC. Here, the expression levels of FOXP3 and CTLA4 were used to construct a subtyping system based on >1200 CRC patients from multiple independent public datasets. We revealed that, in CRC patients with relatively high expression of FOXP3, there exist two different subpopulations with opposite survival patterns according to CLTA4 expression. We further established a method for evaluating all cohorts and identified a novel FOXP3^HighCTLA4^High* CRC risk subpopulation that accounts for 5-10% of CRC patients. Moreover, different methods of functional enrichment and immune evaluation were used to analyze the TME characteristics of different FOXP3/CTLA4 subtypes. The FOXP3^HighCTLA4^High* CRC risk subpopulation was characterized by an immune overdrive TME phenotype, including high immune cell infiltration, low tumor purity, high immune checkpoint levels, and TGF-β activation. Finally, the constructed FOXP3/CTLA4 subtyping system was further validated by quantitative RT-PCR, immunochemistry staining, and multicolor immunofluorescence in an independent CRC cohort we collected. This high-risk subpopulation was also observed in kidney cancers and low-grade glioma patients by a Pan-cancer analysis. Together, our study revealed that the established FOXP3/CTLA4 molecular subtyping system could be used to select treatment and management strategies for CRC and other cancers.

Systemic transcriptome comparison between early- And late-onset pre-eclampsia shows distinct pathology and novel biomarkers

Objectives

Pre‐eclampsia is a leading cause of morbidity and mortality during pregnancy. Although the two forms of this disorder, early‐ (EOPE) and late‐onset of pre‐eclampsia (LOPE) are different, the underlying pathology remains elusive. We aim to unravel the difference and to identify novel biomarkers for EOPE and LOPE.

Materials and Methods

A complete comparison of both placental and peripheral blood transcriptomes was performed to investigate the pathology of pre‐eclampsia. Single‐cell transcriptomics of the maternal‐fetal interface were integrated to identify novel biomarkers for EOPE and LOPE which were further verified at protein or mRNA level in patients.

Results

We found that the transcriptomes of placentae from EOPE, but not LOPE, were significantly different from their respective controls. Conversely, the transcriptomes of peripheral blood from LOPE were more different from their controls than EOPE. Importantly, we identified that several classical biomarkers of pre‐eclampsia were expressed specifically in extravillous trophoblast and syncytiotrophoblast and only upregulated in EOPE, suggesting they should not be applied to all pre‐eclampsia patients in general. We further identified novel biomarkers for EOPE and LOPE from differentially expressed genes (DEGs) of placental and peripheral blood, respectively. The new biomarkers EBI3, IGF2, ORMDL3, GATA2 and KIR2DL4 were experimentally verified with patient blood samples.

Conclusion

Our data demonstrate distinct pathology of EOPE and LOPE, and uncover new biomarkers that can be applied in diagnosis for pre‐eclampsia.

Development and validation of a novel metabolic signature for predicting prognosis in patients with laryngeal cancer

PURPOSE

Despite advances in the development of treatments for laryngeal cancer (LC), including surgical treatments and radio-chemotherapy, the survival rate of LC remains low. Therefore, novel metabolic signatures are urgently needed to evaluate the prognosis of LC patients.

METHODS

Differentially expressed metabolic genes were extracted via bioinformatics analysis from the raw data of The Cancer Genome Atlas and Gene Expression Omnibus (GEO) databases. Univariate Cox regression and LASSO analyses were performed to identify metabolic genes that were significantly correlated with overall survival (OS). Using the Kaplan-Meier analysis and receiver operating characteristics, the prognostic power of candidate signatures was evaluated in the two databases. Gene Set Enrichment Analysis (GSEA) was performed to explore significant signaling pathways and underlying mechanisms in the high- and low-risk groups.

RESULTS

Thirteen metabolism genes showed superior ability to predict OS for LC when compared to clinical variables, and patients in the high-risk group showed significantly poorer OS than those in the low-risk group. The area under the curve of receiver operating curves for 5- and 3-year OS was 0.929 and 0.899, respectively, which were better than the OS obtained with clinicopathological variables. Similar results obtained in the GEO cohort indicated that this gene signature could differentiate between LC patients with and without recurrence.

CONCLUSION

To our knowledge, this study is the first to report that the 13 metabolic genes could serve as an independent biomarker for LC, which could provide vital prognostic information and prediction for personalized treatment of LC.

Multigenerational maternal obesity increases the incidence of HCC in offspring via miR-27a-3p

BACKGROUND & AIMS

Obesity is an independent risk factor for malignancies, including hepatocellular carcinoma (HCC). However, it remains unknown whether maternal obesity affects the incidence of HCC in offspring. Thus, we aimed to investigate this association and its underlying mechanisms.

METHODS

Diethylnitrosamine (DEN) was used to induce HCC in a high-fat diet (HFD)-induced multigenerational obesity model. RNA-sequencing was performed to identify the genes and microRNAs (miRNAs) that were altered over generations. The role of the miR-27a-3p-Acsl1/Aldh2 axis in HCC was evaluated in cell lines and HCC-bearing nude mice, and its intergenerational impact was studied in pregnant mice and their offspring.

RESULTS

Under HFD stress, maternal obesity caused susceptibility of offspring to DEN-induced HCC, and such susceptibility was cumulative over generations. We identified that Acsl1 and Aldh2, direct targets of miR-27a-3p, were gradually changed over generations. Under hyperlipidemic conditions, downregulation of Acsl1 and Aldh2 increased cell proliferation (in vitro) or tumor growth (in vivo) in synergy. Intratumor injection of an miR-27a-3p agomir exacerbated tumor growth by downregulating Acsl1 and Aldh2; while intratumor injection of an miR-27a-3p antagomir had the opposite effect. Moreover, serum miR-27a-3p levels gradually increased in the HFD-fed maternal lineage over generations. Injecting pregnant mice with an miR-27a-3p agomir not only upregulated hepatic miR-27a-3p and downregulated Acsl1/Aldh2 in offspring (fetus, young and adult stages), but also exacerbated HCC development in DEN-treated offspring. In human HCC, upregulated miR-27a-3p and downregulated Acsl1/Aldh2 were negatively correlated with survival on TCGA analysis; while, hepatic miR-27a-3p was negatively correlated with Acsl1/Aldh2 expression in tumor/non-tumor tissues from fatty/non-fatty livers.

CONCLUSIONS

Maternal obesity plays a role in regulating cumulative susceptibility to HCC development in offspring over multiple generations through the miR-27a-3p-Acsl1/Aldh2 axis.

LAY SUMMARY

It is not currently known how maternal obesity affects the incidence of liver cancer in offspring. In this study, we identified a microRNA (miR-27a-3p) that was upregulated in obese mothers and could be passed on to their offspring. This microRNA enhanced the risk of liver cancer in offspring by regulating 2 genes (Acsl1 and Aldh2). This mechanism could be a future therapeutic target.

Effect of bovine leukemia virus (BLV) infection on bovine mammary epithelial cells RNA-seq transcriptome profile

Bovine leukemia virus (BLV) is a δ-retrovirus responsible for Enzootic Bovine Leukosis (EBL), a lymphoproliferative disease that affects cattle. The virus causes immune system deregulation, favoring the development of secondary infections. In that context, mastitis incidence is believed to be increased in BLV infected cattle. The aim of this study was to analyze the transcriptome profile of a BLV infected mammary epithelial cell line (MAC-T). Our results show that BLV infected MAC-T cells have an altered expression of IFN I signal pathway and genes involved in defense response to virus, as well as a collagen catabolic process and some protooncogenes and tumor suppressor genes. Our results provide evidence to better understand the effect of BLV on bovine mammary epithelial cell's immune response.

GLIS2 promotes colorectal cancer through repressing enhancer activation

Gene transcription is coordinately regulated by multiple transcription factors. However, a systematic approach is still lacking to identify co-regulators for transcription factors. Here, we performed ChIP-Seq analysis and predicted the regulators for p53-mediated transcription process, from which we confirmed the roles of GLIS2, MAZ and MEF2A in regulating p53 target genes. We revealed that GLIS2 selectively regulates the transcription of PUMA but not p21. GLIS2 deficiency caused the elevation of H3K27ac and p53 binding on the PUMA enhancer, and promoted PUMA expression. It increased the rate of apoptosis, but not cell cycle. Moreover, GLIS2 represses H3K27ac level on enhancers, regulates the gene expression related with focal adhesion and promotes cell migration, through inhibiting p300. Big data analysis supports GLIS2 as an oncogene in colon cancer, and perhaps other cancers. Taken together, we have predicted candidates for p53 transcriptional regulators, and provided evidence for GLIS2 as an oncogene through repressing enhancer activation.