Molecular and Clinical Characterization of a Novel Prognostic and Immunologic Biomarker FAM111A in Diffuse Lower-Grade Glioma

Unravelling the Intricate Roles of FAM111A and FAM111B: From Protease-Mediated Cellular Processes to Disease Implications

Proteases are critical enzymes in cellular processes which regulate intricate events like cellular proliferation, differentiation and apoptosis. This review highlights the multifaceted roles of the serine proteases FAM111A and FAM111B, exploring their impact on cellular functions and diseases. FAM111A is implicated in DNA replication and replication fork protection, thereby maintaining genome integrity. Additionally, FAM111A functions as an antiviral factor against DNA and RNA viruses. Apart from being involved in DNA repair, FAM111B, a paralog of FAM111A, participates in cell cycle regulation and apoptosis. It influences the apoptotic pathway by upregulating anti-apoptotic proteins and modulating cell cycle-related proteins. Furthermore, FAM111B's association with nucleoporins suggests its involvement in nucleo-cytoplasmic trafficking and plays a role in maintaining normal telomere length. FAM111A and FAM111B also exhibit some interconnectedness and functional similarity despite their distinct roles in cellular processes and associated diseases resulting from their dysfunction. FAM111A and FAM111B dysregulation are linked to genetic disorders: Kenny-Caffey Syndrome type 2 and Gracile Bone Dysplasia for FAM111A and POIKTMP, respectively, and cancers. Therefore, the dysregulation of these proteases in diseases emphasizes their potential as diagnostic markers and therapeutic targets. Future research is essential to unravel the intricate mechanisms governing FAM111A and FAM111B and explore their therapeutic implications comprehensively.

N6 -methyladenosine-modified FAM111A-DT promotes hepatocellular carcinoma growth via epigenetically activating FAM111A

As an epitranscriptomic modulation manner, N6 -methyladenosine (m6 A) modification plays important roles in various diseases, including hepatocellular carcinoma (HCC). m6 A modification affects the fate of RNAs. The potential contributions of m6 A to the functions of RNA still need further investigation. In this study, we identified long noncoding RNA FAM111A-DT as an m6 A-modified RNA and confirmed three m6 A sites on FAM111A-DT. The m6 A modification level of FAM111A-DT was increased in HCC tissues and cell lines, and increased m6 A level was correlated with poor survival of HCC patients. m6 A modification increased the stability of FAM111A-DT transcript, whose expression level showed similar clinical relevance to that of the m6 A level of FAM111A-DT. Functional assays found that only m6 A-modified FAM111A-DT promoted HCC cellular proliferation, DNA replication, and HCC tumor growth. Mutation of m6 A sites on FAM111A-DT abolished the roles of FAM111A-DT. Mechanistic investigations found that m6 A-modified FAM111A-DT bound to FAM111A promoter and also interacted with m6 A reader YTHDC1, which further bound and recruited histone demethylase KDM3B to FAM111A promoter, leading to the reduction of the repressive histone mark H3K9me2 and transcriptional activation of FAM111A. The expression of FAM111A was positively correlated with the m6 A level of FAM111A-DT, and the expression of methyltransferase complex, YTHDC1, and KDM3B in HCC tissues. Depletion of FAM111A largely attenuated the roles of m6 A-modified FAM111A-DT in HCC. In summary, the m6 A-modified FAM111A-DT/YTHDC1/KDM3B/FAM111A regulatory axis promoted HCC growth and represented a candidate therapeutic target for HCC.

Quantitative Evaluation of Stem-like Markers of Human Glioblastoma Using Single-Cell RNA Sequencing Datasets

Targeting glioblastoma (GBM) stem-like cells (GSCs) is a common interest in both the laboratory investigation and clinical treatment of GBM. Most of the currently applied GBM stem-like markers lack validation and comparison with common standards regarding their efficiency and feasibility in various targeting methods. Using single-cell RNA sequencing datasets from 37 GBM patients, we obtained a large pool of 2173 GBM stem-like marker candidates. To evaluate and select these candidates quantitatively, we characterized the efficiency of the candidate markers in targeting the GBM stem-like cells by their frequencies and significance of being the stem-like cluster markers. This was followed by further selection based on either their differential expression in GBM stem-like cells compared with normal brain cells or their relative expression level compared with other expressed genes. The cellular location of the translated protein was also considered. Different combinations of selection criteria highlight different markers for different application scenarios. By comparing the commonly used GSCs marker CD133 (PROM1) with markers selected by our method regarding their universality, significance, and abundance, we revealed the limitations of CD133 as a GBM stem-like marker. Overall, we propose BCAN, PTPRZ1, SOX4, etc. for laboratory-based assays with samples free of normal cells. For in vivo targeting applications that require high efficiency in targeting the stem-like subtype, the ability to distinguish GSCs from normal brain cells, and a high expression level, we recommend the intracellular marker TUBB3 and the surface markers PTPRS and GPR56.

N6-Methyladenosine-modified lncRNA LINREP promotes Glioblastoma progression by recruiting the PTBP1/HuR complex

Glioblastoma multiforme (GBM) is acknowledged as the most aggressive primary brain tumor in adults. It is typically characterized by the high heterogeneity which corresponds to extensive genetic mutations and complex alternative splicing (AS) profiles. Known as a major repressive splicing factor in AS, polypyrimidine tract-binding protein 1 (PTBP1) is involved in the exon skipping events of multiple precursor mRNAs (pre-mRNAs) in GBM. However, precise mechanisms that modulate the expression and activity of PTBP1 remain to be elucidated. In present study, we provided evidences for the role of a long intergenic noncoding RNA (LINREP) implicated in the regulation of PTBP1-induced AS. LINREP interacted with PTBP1 and human antigen R (HuR, ELAVL1) protein complex and protected PTBP1 from the ubiquitin-proteasome degradation. Consequently, a broad spectrum of PTBP1-induced spliced variants was generated by exon skipping, especially for the skipping of reticulon 4 (RTN4) exon 3. Interestingly, LINREP also promoted the dissociation of nuclear UPF1 from PTBP1, which increased the binding of PTBP1 to RTN4 transcripts, thus enhancing the skipping of RTN4 exon 3 to some extent. Besides, HuR recruitment was essential for the stabilization of LINREP via a manner dependent on N6-methyladenosine (m6A) formation and identification. Taken together, our results demonstrated the functional significance of LINREP in human GBM for its dual regulation of PTBP1-induced AS and its m6A modification modality, implicating that HuR/LINREP/PTBP1 axis might serve as a potential therapeutic target for GBM.

Imputation of single-cell transcriptome data enables the reconstruction of networks predictive of breast cancer metastasis

Single-cell transcriptome data provide a unique opportunity to explore the gene networks of a particular cell type. However, insufficient capture rate and high dimensionality of single-cell RNA sequencing (scRNA-seq) data challenge cell-type-specific gene network (CGN) reconstruction. Here, we demonstrated that the imputation of scRNA-seq data enables reconstruction of CGNs by effective retrieval of gene functional associations. We reconstructed CGNs for seven primary and nine metastatic breast cancer cell lines using scRNA-seq data with imputation. Key genes for primary or metastatic cell lines were prioritized based on network centrality measures and CGN hub genes that were presumed to be the major determinant of cell type characteristics. To identify novel genes in breast cancer metastasis, we used the average rank difference of centrality between the primary and metastatic cell lines. Genes predicted using CGN centrality analysis were more enriched for known breast cancer metastatic genes than those predicted using differential expression. The molecular chaperone CCT2 was identified as a novel gene for breast metastasis during knockdown assays of several candidate genes. Overall, our study demonstrated an effective CGN reconstruction technique with imputation of scRNA-seq data and the feasibility of identifying key genes for particular cell subsets using single-cell network analysis.

Identification of circRNA-lncRNA-miRNA-mRNA competitive endogenous RNA networks as prognostic bio-markers in head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is one of the most prevalent neoplasms of malignancy. The pathogenesis of HNSCC is involved in the change of gene expression and the abnormal interaction between genes.

OBJECTIVE

The purpose of our research is to find genes that affect the survival time of patients and construct circRNA-lncRNA-miRNA-mRNA network. This paper aims to reveal the development mechanism of HNSCC and provide new ideas for clinical prognosis prediction.

METHODS

Transcriptome data were obtained from TCGA and GEO databases. A total of 556 samples were included. We screened the differentially expressed lncRNAs, circRNAs and mRNAs to predict interactions between lncRNA-miRNA, circRNA-miRNA and miRNA-mRNA. Then, we selected the hub mRNAs to predict the prognosis of patients.

RESULTS

We identified 8 hub prognostic mRNAs (FAM111A, IFIT2, CAV1, KLF9, OCIAD2, SLCO3A1, S100A16 and APOL3) by a series of bioinformatics analyses and established the prognostic ceRNA network of 8 mRNAs, 3 miRNAs, 2 circRNAs and 1 lncRNAs according to the targeting relationship by using databases.

CONCLUSION

We established the circRNA-lncRNA-miRNA-mRNA gene interaction network in HNSCC. We illuminated the molecular mechanism underlying the gene regulation associated with the pathogenesis of HNSCC and predicted the biomarkers related to prognosis.

Comprehensive Analysis of the Prognostic Value and Immune Infiltration of Butyrophilin Subfamily 2/3 (BTN2/3) Members in Pan-Glioma

The BTN2/3 subfamilies are overexpressed in many cancers, including pan-glioma (low- and high-grade gliomas). However, the expression and prognosis of BTN2/3 subfamilies and tumor-infiltrating lymphocytes in pan-glioma remain unknown. In the present study, we systematically explored and validated the expression and prognostic value of BTN2/3 subfamily members in pan-glioma [The Cancer Genome Atlas–glioblastoma and low-grade glioma (TCGA-GBMLGG) merge cohort] using multiple public databases. We used clinical specimens for high-throughput verification and cell lines for qRT-PCR verification, which confirmed the expression profiles of BTN2/3 subfamilies. In addition, the function of the BTN2/3 subfamily members and the correlations between BTN2/3 subfamily expression and pan-glioma immune infiltration levels were investigated. We found that BTN2/3 subfamily members were rarely mutated. BTN2/3 subfamilies were overexpressed in pan-glioma; high expression of BTN2/3 subfamily members was correlated with poor prognosis. In addition, BTN2/3 subfamilies might positively regulate proliferation, and the overexpression of BTN2/3 subfamilies influenced cell cycle, differentiation, and glioma stemness. In terms of immune infiltrating levels, BTN2/3 subfamily expression was positively associated with CD4+ T-cell, B-cell, neutrophil, macrophage, and dendritic cell infiltrating levels. These findings suggest that BTN2/3 subfamily expression is correlated with prognosis and immune infiltration levels in glioma. Therefore, the BTN2/3 subfamilies can be used as biomarkers for pan-glioma and prognostic biomarkers for determining the prognosis and immune infiltration levels in pan-glioma.

Development and Validation of an Inflammatory Response-Related Gene Signature for Predicting the Prognosis of Pancreatic Adenocarcinoma

Pancreatic adenocarcinoma (PAAD) is a highly dangerous malignant tumor of the digestive tract, and difficult to diagnose, treat, and predict the prognosis. As we all know, tumor and inflammation can affect each other, and thus the inflammatory response in the microenvironment can be used to affect the prognosis. So far, the prognostic value of inflammatory response-related genes in PAAD is still unclear. Therefore, this study aimed to explore the inflammatory response-related genes for predicting the prognosis of PAAD. In this study, the mRNA expression profiles of PAAD patients and the corresponding clinical characteristics data of PAAD patients were downloaded from the public database. The least absolute shrinkage and selection operator (LASSO) Cox analysis model was used to identify and construct the prognostic gene signature in The Cancer Genome Atlas (TCGA) cohort. The PAAD patients used for verification are from the International Cancer Genome Consortium (ICGC) cohort. The Kaplan-Meier method was used to compare the overall survival (OS) between the high- and low-risk groups. Univariate and multivariate Cox analyses were performed to identify the independent predictors of OS. Gene set enrichment analysis (GSEA) was performed to obtain gene ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and the correlation between gene expression and immune infiltrates was investigated via single sample gene set enrichment analysis (ssGSEA). The GEPIA database was performed to examine prognostic genes in PAAD. LASSO Cox regression analysis was used to construct a model of inflammatory response-related gene signature. Compared with the low-risk group, patients in the high-risk group had significantly lower OS. The receiver operating characteristic curve (ROC) analysis confirmed the signature's predictive capacity. Multivariate Cox analysis showed that risk score is an independent predictor of OS. Functional analysis shows that the immune status between the two risk groups is significantly different, and the cancer-related pathways were abundant in the high-risk group. Moreover, the risk score is significantly related to tumor grade, stage, and immune infiltration types. It was also obtained that the expression level of prognostic genes was significantly correlated with the sensitivity of cancer cells to anti-tumor drugs. In addition, there are significant differences in the expression of PAAD tissues and adjacent non-tumor tissues. The novel signature constructed from five inflammatory response-related genes can be used to predict prognosis and affect the immune status of PAAD. In addition, suppressing these genes may be a treatment option.

The Circ_0001367/miR-545-3p/LUZP1 Axis Regulates Cell Proliferation, Migration and Invasion in Glioma Cells

Glioma is the most common primary intracranial malignant tumour in adults. It has a high incidence and poses a serious threat to human health. Circular RNA is a hotspot of cancer research. In this study, we aimed to explore the role of circ_0001367 in gliomagenesis and the underlying mechanism. First, qRT-PCR was conducted, which showed that circ_0001367 level was downregulated in glioma tissues and cells. Next, gain-of-function and loss-of-function assays were performed, which indicated that circ_0001367 inhibited the proliferation, migration and invasion of glioma cells. Subsequent bioinformatics analysis, dual-luciferase reporter assays, RNA immunoprecipitation assays and cell function assays demonstrated that circ_0001367 inhibited the proliferation, migration and invasion of glioma cells by absorbing miR-545-3p and thereby regulating the expression of leucine zipper protein (LUZP1). Finally, an in vivo experiment was conducted, which demonstrated that circ_0001367 inhibited glioma growth in vivo by modulating miR-545-3p and LUZP1. Taken together, the results of this study demonstrate that the circ_0001367/miR-545-3p/LUZP1 axis may be a novel target for glioma therapy.