Construction of a circRNA-miRNA-mRNA Network Related to Macrophage Infiltration in Hepatocellular Carcinoma

HGATLink: single-cell gene regulatory network inference via the fusion of heterogeneous graph attention networks and transformer

BACKGROUND

Gene regulatory networks (GRNs) involve complex regulatory relationships between genes and play important roles in the study of various biological systems and diseases. The introduction of single-cell sequencing (scRNA-seq) technology has allowed gene regulation studies to be carried out on specific cell types, providing the opportunity to accurately infer gene regulatory networks. However, the sparsity and noise problems of single-cell sequencing data pose challenges for gene regulatory network inference, and although many gene regulatory network inference methods have been proposed, they often fail to eliminate transitive interactions or do not address multilevel relationships and nonlinear features in the graph data well.

RESULTS

On the basis of the above limitations, we propose a gene regulatory network inference framework named HGATLink. HGATLink combines the heterogeneous graph attention network and simplified transformer to capture complex interactions effectively between genes in low-dimensional space via matrix decomposition techniques, which not only enhances the ability to model complex heterogeneous graph structures and alleviate transitive interactions, but also effectively captures the long-range dependencies between genes to ensure more accurate prediction.

CONCLUSIONS

Compared with 10 state-of-the-art GRN inference methods on 14 scRNA-seq datasets under two metrics, AUROC and AUPRC, HGATLink shows good stability and accuracy in gene regulatory network inference tasks.

Comprehensive analysis of circRNA-miRNA-mRNA network related to angiogenesis in recurrent implantation failure

BACKGROUND

Abnormal endometrial blood flow causes a decrease in endometrial receptivity and is considered a relatively independent risk factor for recurrent implantation failure (RIF). This study aimed to explore the potentially functional circRNA-miRNA-mRNA network in RIF, and further explore its mechanism.

METHODS

Datasets were downloaded from the GEO database to identify differentially expressed circRNAs, miRNAs and mRNAs. The circRNA-miRNA-mRNA and PPI networks were constructed using Cytoscape 3.6.0 and the STRING database, the hub genes were identified with the cytoHubba plug-in, and a circRNA-miRNA-hub mRNA regulatory sub-network was constructed. Then, GO and KEGG pathway enrichment analyses of the hub genes were performed to comprehensively analyze the mechanism of hub mRNAs in RIF. Due to the results of circRNAs-miRNAs-hub mRNAs regulatory network, we verified the expression of circRNA_0001721, circRNA_0000714, miR-17-5p, miR-29b-3p, HIF1A and VEGFA in the RIF mouse model by qRT‒PCR and western blotting.

RESULTS

We initially identified 175 DEmRNAs, 48 DEmiRNAs and 56 DEcircRNAs in RIF associated with angiogenesis and constructed a circRNA-miRNA‒mRNA network and PPI network. We further identified six hub genes in the acquired network. Based on these genes, functional enrichment analysis revealed that the HIF-1 signaling pathway plays a vital role in endometrial angiogenesis in RIF. In addition, the interaction networks of circRNA_0001721/miR-17-5p/HIF1A and the circRNA_0000714/miR-29b-3p/VEGFA axis were predicted. In the RIF mouse model, circRNA_0001721, circRNA_0000714, HIF1A and VEGFA were down-regulated, whereas miR-17-5p and miR-29b-3p were up-regulated according to qRT‒PCR and western blotting.

CONCLUSION

This study revealed that the HIF-1 signaling pathway plays a vital role in endometrial angiogenesis in RIF. The circRNA_0001721/miR-17-5p/HIF1A and circRNA_0000714/miR-29b-3p/VEGFA axes might play a role in the pathogenesis of endometrial angiogenesis in RIF.

Construction of the underlying circRNA-miRNA-mRNA regulatory network and a new diagnostic model in ulcerative colitis by bioinformatics analysis

BACKGROUND

Circular RNAs (circRNAs) are involved in the pathogenesis of many diseases through competing endogenous RNA (ceRNA) regulatory mechanisms.

AIM

To investigate a circRNA-related ceRNA regulatory network and a new predictive model by circRNA to understand the diagnostic mechanism of circRNAs in ulcerative colitis (UC).

METHODS

We obtained gene expression profiles of circRNAs, miRNAs, and mRNAs in UC from the Gene Expression Omnibus dataset. The circRNA-miRNA-mRNA network was constructed based on circRNA-miRNA and miRNA-mRNA interactions. Functional enrichment analysis was performed to identify the biological mechanisms involved in circRNAs. We identified the most relevant differential circRNAs for diagnosing UC and constructed a new predictive nomogram, whose efficacy was tested with the C-index, receiver operating characteristic curve (ROC), and decision curve analysis (DCA).

RESULTS

A circRNA-miRNA-mRNA regulatory network was obtained, containing 12 circRNAs, three miRNAs, and 38 mRNAs. Two optimal prognostic-related differentially expressed circRNAs, hsa_circ_0085323 and hsa_circ_0036906, were included to construct a predictive nomogram. The model showed good discrimination, with a C-index of 1(> 0.9, high accuracy). ROC and DCA suggested that the nomogram had a beneficial diagnostic ability.

CONCLUSION

This novel predictive nomogram incorporating hsa_circ_0085323 and hsa_circ_0036906 can be conveniently used to predict the risk of UC. The circRNa-miRNA-mRNA network in UC could be more clinically significant.

Research and experimental verification on the mechanisms of cellular senescence in triple-negative breast cancer

Background

Triple-negative breast cancer (TNBC) is an aggressive breast cancer subtype with high heterogeneity, poor prognosis, and a low 10-year survival rate of less than 50%. Although cellular senescence displays extensive effects on cancer, the comprehensions of cellular senescence-related characteristics in TNBC patients remains obscure.

Method

Single-cell RNA sequencing (scRNA-seq) data were analyzed by Seurat package. Scores for cellular senescence-related pathways were computed by single-sample gene set enrichment analysis (ssGSEA). Subsequently, unsupervised consensus clustering was performed for molecular cluster identification. Immune scores of patients in The Cancer Genome Atlas (TCGA) dataset and associated immune cell scores were calculated using Estimation of STromal and Immune cells in MAlignantTumours using Expression data (ESTIMATE) and Microenvironment Cell Populations-counter (MCP-counter), Tumor Immune Estimation Resource (TIMER) and Estimating the Proportion of Immune and Cancer cells (EPIC) methods, respectively. Immunotherapy scores were assessed using TIDE. Furthermore, feature genes were identified by univariate Cox and Least Absolute Shrinkage and Selection Operator (LASSO) regression analyses; these were used to construct a risk model. Additionally, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and transwell assay were conducted for in vitro validation of hub genes.

Result

TNBC was classified into three subtypes based on cellular senescence-related pathways as clusters 1, 2, and 3. Specifically, cluster 1 showed the best prognosis, followed by cluster 2 and cluster 3. The levels of gene expression in cluster 2 were the lowest, whereas these were the highest in cluster 3. Moreover, clusters 1 and 3 showed a high degree of immune infiltration. TIDE scores were higher for cluster 3, suggesting that immune escape was more likely in patients with the cluster 3 subtype who were less likely to benefit from immunotherapy. Next, the TNBC risk model was constructed and validated. RT-qPCR revealed that prognostic risk genes (MMP28, ACP5 and KRT6A) were up-regulated while protective genes (CT83) were down-regulated in TNBC cell lines, validating the results of the bioinformatics analysis. Meanwhile, cellular experiments revealed that ACP5 could promote the migration and invasion abilities in two TNBC cell lines. Finally, we evaluated the validity of prognostic models for assessing TME characteristics and TNBC chemotherapy response.

Conclusion

In conclusion, these findings help to assess the efficacy of targeted therapies in patients with different molecular subtypes, have practical applications for subtype-specific treatment of TNBC patients, and provide information on prognostic factors, as well as guidance for the revelation of the molecular mechanisms by which senescence-associated genes influence TNBC progression.

LINC02561 promotes metastasis in HCC via HIF1-α/NDRG1/UPF1 axis

In the entire world, hepatocellular carcinoma (HCC) is one of the most frequent cancers that lead to death. Experiments on the function of long non-coding RNAs in the emergence of malignancies, including HCC, are ongoing. As a crucial RNA monitoring mechanism in eucaryotic cells, nonsense-mediated mRNA decay (NMD) can recognize and destroy mRNAs, which has an premature termination codons (PTC) in the open reading frame to prevent harmful buildup of truncated protein products in the cells. Nonsense transcript regulator 1 (Up-frameshift suppressor 1, UPF1), as a highly conserved RNA helicase and ATPase, plays a key role in NMD. Our laboratory screened out the highly expressed lncRNA LINC02561 in HCC from the TCGA database. Further research found that LINC02561 enhanced the invasion and transition abilities of liver cancer cells by regulating the protein N-Myc downstream regulated 1 (NDRG1). Hypoxia inducible factor-1 (HIF-1α) can bonded to LINC02561 promoters under hypoxic conditions, thereby promoting the upregulation of LINC02561 expression in liver cancer cells. LINC02561 competes with NDRG1 mRNA to bind UPF1, thereby preventing the degradation of NDRG1 mRNA to facilitate NDRG1 protein level. Taken together, the HIF1α-LINC02561-UPF1-NDRG1 regulatory axis could be an entirely novel target of liver cancer-related treatment.

Prognostic related gene index for predicting survival and immunotherapeutic effect of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common malignant liver tumor. It is an aggressive disease with high mortality rate. In this study, we investigated a new prognosis-related gene index (PRGI) that can predict the survival and efficacy of immunotherapy in patients with HCC. RNA-seq data and clinical data of HCC samples were obtained from the cancer genome atlas and ICGC databases. Prognosis-related genes were obtained using log-rank tests and univariate Cox proportional hazards regression. Univariate and multivariate analyses were performed on the overall survival rate of patients with prognosis-related genes and multiple clinicopathological factors, and a nomogram was constructed. A PRGI was then constructed based on least absolute shrinkage and selection operator or multivariate Cox Iterative Regression. The possible correlation between PRGI and immune cell infiltration or immunotherapy efficacy was discussed. Eight genes were identified to construct the PRGI. PRGI can predict the infiltration of immune cells into the tumor microenvironment of HCC and the response to immunotherapy. PRGI can accurately predict the survival rate of patients with HCC, reflect the immune microenvironment, and predict the efficacy of immunotherapy.

E2F3 renders an immunosuppressive tumor microenvironment in nasopharyngeal carcinoma: Involvements of the transcription activation of PRC1 and BIRC5

BACKGROUND

E2F transcription factors are well-recognized oncogenic molecules, and their correlation with immune cell infiltration has recently been reported. This work studies the impacts and mechanism of E2F transcription factor 3 (E2F3) in the growth and tumor microenvironment (TME) of nasopharyngeal carcinoma (NPC).

METHODS

Aberrantly expressed transcription factors in NPC were screened by abundant bioinformatics analyses. Gene expression in NPC cells was analyzed by reverse transcription-quantitative polymerase chain reaction and Western blot analyses. Malignant behaviors of NPC cells were analyzed by cell counting kit-8, 5-ethynyl-2'-deoxyuridine labeling, Transwell assays, and xenograft tumor models. TPA-induced THP-1 cells (macrophages) were cultured in the conditioned medium of NPC cells to mimic tumor-associated macrophages (TAMs) in vivo, and these TAMs were cocultured with CD8+ T cells. Regulation of E2F3 on protein regulator of cytokinesis 1 (PRC1) and baculoviral IAP repeat containing 5 (BIRC5) was validated by chromatin immunoprecipitation and luciferase reporter assays.

RESULTS

E2F3 was highly expressed in NPC cells, and its knockdown suppressed malignant behavior and tumorigenic ability of the cells. The E2F3 knockdown condition downregulated M2 cytokines CD163 and interleukin-10 in TAMs, which further enhanced proliferation and activation of the cocultured CD8+ T cells. E2F3 promoted transcription of PRC1 and BRIC5. Furthermore, PRC1 or BRIC5 upregulation in NPC cells restored the malignant properties of NPC cells, reprogrammed the TAMs to M2 phenotype, and suppressed the CD8+ T cell proliferation and activation.

CONCLUSION

This work suggests that E2F3 renders an immunosuppressive TME in NPC by activating PRC1 and BIRC5. Suppression of any member involved might favor tumor elimination.

Circular RNAs in hepatocellular carcinoma: biogenesis, function, and pathology

Hepatocellular carcinoma (HCC) is one of the most common causes of cancer-related death worldwide. Both genetic and environmental factors through a multitude of underlying molecular mechanisms participate in the pathogenesis of HCC. Recently, numerous studies have shown that circular RNAs (circRNAs), an emerging class of non-coding RNAs characterized by the presence of covalent bonds linking 3' and 5' ends, play an important role in the initiation and progression of cancers, including HCC. In this review, we outline the current status of the field of circRNAs, with an emphasis on the functions and mechanisms of circRNAs in HCC and its microenvironment. We also summarize and discuss recent advances of circRNAs as biomarkers and therapeutic targets. These efforts are anticipated to throw new insights into future perspectives about circRNAs in basic, translational and clinical research, eventually advancing the diagnosis, prevention and treatment of HCC.

Noncoding RNA-mediated macrophage and cancer cell crosstalk in hepatocellular carcinoma

The tumor microenvironment (TME) is a well-recognized system that plays an essential role in tumor initiation, development, and progression. Intense intercellular communication between tumor cells and other cells (especially macrophages) occurs in the TME and is mediated by cell-to-cell contact and/or soluble messengers. Emerging evidence indicates that noncoding RNAs (ncRNAs) are critical regulators of the relationship between cells within the TME. In this review, we provide an update on the regulation of ncRNAs (primarily micro RNAs [miRNAs], long ncRNAs [lncRNAs], and circular RNAs [circRNAs]) in the crosstalk between macrophages and tumor cells in hepatocellular carcinoma (HCC). These ncRNAs are derived from macrophages or tumor cells and act as oncogenes or tumor suppressors, contributing to tumor progression not only by regulating the physiological and pathological processes of tumor cells but also by controlling macrophage infiltration, activation, polarization, and function. Herein, we also explore the options available for clinical therapeutic strategies targeting crosstalk-related ncRNAs to treat HCC. A better understanding of the relationship between macrophages and tumor cells mediated by ncRNAs will uncover new diagnostic biomarkers and pharmacological targets in cancer.

Role of Circular RNAs in the Regulation of Immune Cells in Response to Cancer Therapies

Circular RNAs (CircRNAs) are a class of small endogenous noncoding RNA that are formed by means of either the spliceosome or lariat-type splicing. CircRNAs have multiple regulatory functions and have been detected in different cell types, like normal, tumor and immune cells. CircRNAs have been suggested to regulate T cell functions in response to cancer. CircRNAs can enter into T cells and promote the expression of molecules that either trigger antitumoral responses or promote suppression and the consequent evasion to the immune response. Additionally, circRNAs may promote tumor progression and resistance to anticancer treatment in different types of neoplasias. In this minireview we discuss the impact of circRNAs and its function in the regulation of the T-cells in immune response caused by cancer therapies.

Non-coding RNAs as biomarkers for hepatocellular carcinoma-A systematic review

Hepatocellular carcinoma (HCC) is the sixth most common malignancy in the world and the fourth leading cause of cancer-related death, and its incidence is increasing globally. Despite significant advances in treatment strategies for HCC, the prognosis is still poor due to its high recurrence rate. Therefore, there is an urgent need to understand the pathogenesis of HCC and further develop new therapies to improve the prognosis and quality of life of HCC patients. MicroRNAs (miRNAs, miRs) are small non-coding RNAs involved in post-transcriptional regulation of gene expression that is abnormally expressed in cancer-associated genomic regions or vulnerable sites. More and more findings have shown that miRNAs are important regulatory factors of mRNA expression in HCC, and they are receiving more and more attention as a possible key biomarker of HCC. This review mainly summarizes the potential applied value on miRNAs as diagnostic, drug resistant, prognostic, and therapeutic biomarkers in the diagnosis, therapy, and prognosis of HCC. Also, we summarize the research value of long non-coding RNA (lncRNAs), circular RNAs (circRNAs), and miRNAs network in HCC as novel biomarkers, aiming at providing some references for the therapy of HCC.

Identification of Gene Signature as Diagnostic and Prognostic Blood Biomarker for Early Hepatocellular Carcinoma Using Integrated Cross-Species Transcriptomic and Network Analyses

Background: Hepatocellular carcinoma (HCC) is considered the most common type of liver cancer and the fourth leading cause of cancer-related deaths in the world. Since the disease is usually diagnosed at advanced stages, it has poor prognosis. Therefore, reliable biomarkers are urgently needed for early diagnosis and prognostic assessment. Methods: We used genome-wide gene expression profiling datasets from human and rat early HCC (eHCC) samples to perform integrated genomic and network-based analyses, and discovered gene markers that are expressed in blood and conserved in both species. We then used independent gene expression profiling datasets for peripheral blood mononuclear cells (PBMCs) for eHCC patients and from The Cancer Genome Atlas (TCGA) database to estimate the diagnostic and prognostic performance of the identified gene signature. Furthermore, we performed functional enrichment, interaction networks and pathway analyses. Results: We identified 41 significant genes that are expressed in blood and conserved across species in eHCC. We used comprehensive clinical data from over 600 patients with HCC to verify the diagnostic and prognostic value of 41-gene-signature. We developed a prognostic model and a risk score using the 41-geneset that showed that a high prognostic index is linked to a worse disease outcome. Furthermore, our 41-gene signature predicted disease outcome independently of other clinical factors in multivariate regression analysis. Our data reveals a number of cancer-related pathways and hub genes, including EIF4E, H2AFX, CREB1, GSK3B, TGFBR1, and CCNA2, that may be essential for eHCC progression and confirm our gene signature's ability to detect the disease in its early stages in patients' biological fluids instead of invasive procedures and its prognostic potential. Conclusion: Our findings indicate that integrated cross-species genomic and network analysis may provide reliable markers that are associated with eHCC that may lead to better diagnosis, prognosis, and treatment options.

The Role of Non-Coding RNAs in the Sorafenib Resistance of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the second most common cause of cancer-related death. Sorafenib is approved by the U.S. Food and Drug Administration to be a first-line chemotherapy agent for patients with advanced HCC. A portion of advanced HCC patients can benefit from the treatment with sorafenib, but many patients ultimately develop sorafenib resistance, leading to a poor prognosis. The molecular mechanisms of sorafenib resistance are sophisticated and indefinite. Notably, non-coding RNAs (ncRNAs), which include long ncRNAs (lncRNAs), microRNAs (miRNAs) and circular RNAs (circRNAs), are critically participated in the occurrence and progression of tumors. Moreover, growing evidence has suggested that ncRNAs are crucial regulators in the development of resistance to sorafenib. Herein, we integrally and systematically summarized the molecular mechanisms and vital role of ncRNAs impact sorafenib resistance of HCC, and ultimately explored the potential clinical administrations of ncRNAs as new prognostic biomarkers and therapeutic targets for HCC.

Identification of RNA binding protein interacting with circular RNA and hub candidate network for hepatocellular carcinoma

The interaction between RNA binding protein (RBP) and circular RNA (circRNA) is important for the regulation of tumor progression. This study aimed to identify the RBP-circRNA network in hepatocellular carcinoma (HCC). 22 differentially expressed (DE) circRNAs in HCC were screened out from Gene Expression Omnibus (GEO) database and their binding RBPs were predicted by Circular RNA Interactome. Among them, 17 DERBPs, which were commonly dysregulated in HCC from The Clinical Proteomic Tumor Analysis Consortium (CPTAC), The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) projects, were utilized to construct the RBP-circRNA network. Through survival analysis, we found TARDBP was the only prognostic RBP for HCC in CPTAC, TCGA and ICGC projects. High expression of TARDBP was correlated with high grade, advanced stage and low macrophage infiltration of HCC. Additionally, gene set enrichment analysis showed that dysregulated TARDBP might be involved in some pathways related to the HCC pathogenesis. Therefore, a hub RBP-circRNA network was generated based on TARDBP. RNA immunoprecipitation and RNA pull-down confirmed that hsa_circ_0004913 binds to TARDBP. These findings indicated certain RBP-circRNA regulatory network potentially involved in the pathogenesis of HCC, which provides novel insights into the mechanism study and biomarker identification for HCC.

Hsa_circ_0091581 promotes glioma progression by regulating RMI1 via sponging miR-1243-5p

Glioma is a pervasive malignancy and the main cause of cancer-related deaths worldwide. Circular RNA is an important subject of cancer research, and its role and function in glioma are poorly understood. This study demonstrated that hsa_circ_0091581 is upregulated in glioma tissues and cells. The results of the CCK-8, EdU, and transwell assays indicated that hsa_circ_0091581 promotes proliferation, migration, and invasion of glioma cells. The results of the luciferase reporter and RNA immunoprecipitation assays indicated that the mechanism of the effects of hsa_circ_0091581 on glioma cells involves sponging miR-1243-5p to regulate RMI1. The results of the rescue experiments indicated that hsa_circ_0091581 regulates proliferation, migration, and invasion of glioma cells by targeting RMI1 in a miR-1243-5p dependent manner. The results of the nude mice xenograft assays showed that knockdown of hsa_circ_0091581 inhibits glioma growth in vivo. Thus, our study determined the role of hsa_circ_0091581/miR-1243-5p/RMI1 in glioma and suggests that this axis may be a novel therapeutic target in glioma.

Circular RNAs Sparkle in the Diagnosis and Theranostics of Hepatocellular Carcinoma

Exonic circular RNAs (circRNAs) are a novel subgroup of non-coding RNAs, which are generated by a back-splicing mechanism of the exons or introns. Unlike the linear RNA, circRNA forms a covalently closed loop, and it normally appears more abundant than the linear products of its host gene. Due to the relatively high specificity and stability of circular RNAs in tissues and body fluid, circular RNAs have attracted widely scientific interest for its potential application in cancer diagnosis and as a guide for preclinical therapy, especially for hard-to-treat cancers with high heterogeneity, such as hepatocellular carcinoma (HCC). Thus, we summarize the updated knowledge of circular RNAs, including the mechanism of the generation of endogenous circular RNAs and their regulatory, diagnostic, and therapeutic roles in HCC.

Construction and analysis of macrophage infiltration related circRNA-miRNA-mRNA regulatory networks in hepatocellular carcinoma

BACKGROUND

Macrophage play a crucial role in regulating tumor progression. This study intended to investigate the circular RNA (circRNA) regulatory network associated with macrophage infiltration in hepatocellular carcinoma (HCC).

METHODS

The immune cell fractions of HCC from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium were calculated by Estimation of the Proportion of Immune and Cancer cells algorithm. The differentially expressed mRNAs (DEmRNAs), microRNAs (DEmiRNAs) and circRNAs (DEcircRNAs) were identified from HCC and adjacent non-tumor cases of TCGA or Gene Expression Omnibus database. The DEmRNAs related to macrophage were selected by weighted gene co-expression network analysis and then utilized to generate the circRNA-miRNA-mRNA network. A hub circRNA regulatory network was established based on the co-expressed DEmiRNAs and DEmRNAs owning contrary correlation with the clinical characteristics, survival and macrophage infiltration level. A gene signature based on the DEmRNAs in hub network was also generated for further evaluation. The circRNA binding bite for miRNA was detected by luciferase assay.

RESULTS

High macrophage fraction predicted good survival for HCC. A circRNA-miRNA-mRNA network was constructed by 27 macrophage related DEmRNAs, 21 DEmiRNAs, and 15 DEcircRNAs. Among this network, the expression of hsa-miR-139-5p was negatively correlated with CDCA8, KPNA2, PRC1 or TOP2A. Hsa-miR-139-5p low or targeted DEmRNA high expression was associated with low macrophage infiltration, high grade, advanced stage and poor prognosis of HCC. Additionally, the risk score generated by 4-DEmRNA signature could reflect the macrophage infiltration status and function as an independent prognostic factor for HCC. Finally, hsa_circ_0007456 acting on hsa-miR-139-5p related network was viewed as the hub circRNA regulatory network. Taken together, some circRNA regulatory networks may be associated with macrophage infiltration, which provides clues for mechanism study and therapeutic strategies of HCC.