microRNA-Dependent Modulation of Genes Contributes to ESR1's Effect on ERα Positive Breast Cancer

Identification of a hypoxia-suppressed lncRNA RAMP2-AS1 in breast cancer

Hypoxia is a critical feature of solid tumors and exerts crucial roles in cancers, including breast cancer. However, the detailed relationship between lncRNA-miRNA-mRNA triple network and hypoxia in breast cancer is still indistinct. In this study, a series of in silico analyses and online databases or tools were employed to establish a hypoxia-related lncRNA-miRNA-mRNA network in breast cancer based on competing endogenous RNA mechanism at single-cell resolution. RAMP2-AS1 was, eventually, identified as the most potential lncRNA, which was significantly negatively associated with hypoxia in breast cancer. Compared with normal controls, RAMP2-AS1 was markedly downregulated in breast cancer. Moreover, survival analysis revealed favorable prognostic values of RAMP2-AS1 in total or in specific clinicopathological breast cancer patients. Next, miR-660-5p, miR-2277-5p and miR-1301-3p, upregulated and possessed poor prognostic values in breast cancer, were identified as three potential downstream miRNAs of RAMP2-AS1. Then, the most potential downstream hypoxia-related genes (ATM and MYH11) of RAMP2-AS1/miRNA axis in breast cancer were screened out. Intriguingly, in vitro experiments confirmed that RAMP2-AS1 was a hypoxia-suppressed lncRNA and miR-660-5p/ATM was a potential downstream axis of RAMP2-AS1 in breast cancer. Collectively, our current data elucidated a key hypoxia-suppressed lncRNA RAMP2-AS1 and its possible miRNA-mRNA regulatory mechanism in breast cancer.

APOB is a potential prognostic biomarker in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is significantly associated with adverse prognostic outcomes. The development and progression of different types of human tumors are significantly influenced by APOB. Nevertheless, the significance and pathomechanisms of APOB in HCC have not been conclusively determined. We assessed APOB expression levels in HCC using three publicly available databases of TIMER2.0, UALCAN and Human Protein Atlas. To identify the biological function of APOB, we conducted enrichment analysis via LinkedOmics. Moreover, UALCAN was employed to assess the relationship between APOB expression and clinicopathological features among HCC patients. Additionally, the Kaplan-Meier plotter was utilized to investigate the prognostic relevance of APOB in HCC. To explore potential regulatory ncRNAs that could bind to APOB, we utilized StarBase and GEPIA. Furthermore, the correlation between APOB expression and immune cell infiltration, as well as immune checkpoint genes, was investigated using Spearman's correlation analysis in TISIDB, GEPIA, and TIMER2.0. The findings of our investigation showed a notable decrease in the expression levels of APOB among individuals diagnosed with HCC. Moreover, a noteworthy correlation was observed between the expression of APOB and immune checkpoint genes, alongside the occurrence of immune cell infiltration. The levels of APOB expression in HCC tissues also showed correlations with various clinicopathological features. According to Cox regression analysis, decreased APOB expression emerged as a potential autonomous predictor for OS, RFS, DSS, and PFS among HCC patients. Furthermore, we identified six potential pathways associated with non-coding RNA (ncRNA) as the most promising pathway for APOB in HCC. Our results illuminate the possible involvement of APOB in HCC and offer understanding into its governing mechanisms and medical importance.

LncRNAs-Regulated High Expression of LAMC2 Reveals a Prognostic and Immunological Value in Pancreatic Adenocarcinoma

Pancreatic adenocarcinoma (PAAD) is one of the most hazardous cancers in digestive system, and the prognosis is notoriously bad. Increasing evidences indicate that Laminin Subunit Gamma 2 (LAMC2) is critical for the initiation and the growth of various sorts of human cancers. However, the involved molecular pathways of LAMC2 in PAAD are still poorly understood. In this study, prediction programs and databases were employed to conduct pan-cancer analysis. Multiple variations of human malignancies showed increased LAMC2 expression, which was positively correlated to a poor prognosis in PAAD. Moreover, LAMC2 was positively correlated with the biomarkers of immune cells including CD19, CD163, and NOS2 in PAAD. The lncRNA C5orf66 /PTPRG-AS1- miR-128-3p -LAMC2 axis was identified to be a potential upstream regulatory pathway of LAMC2 in PAAD. Furthermore, LAMC2 upregulation in PAAD was associated with PD-L1 expression, indicating promoting carcinoma immune cell infiltration. Our study elucidated prognostic and immunological values of LAMC2 in PAAD, providing a promise target for PAAD treatment.

Systems genetics and bioinformatics analyses using ESR1-correlated genes identify potential candidates underlying female bone development

Estrogen receptor α (ESR1) is involved in E2 signaling and plays a major role in postmenopausal bone loss. However, the molecular network underlying ESR1 has not been explored. We used systems genetics and bioinformatics to identify important genes associated with Esr1 in postmenopausal bone loss. We identified ~2300 Esr1-coexpressed genes in female BXD bone femur, functional analysis of which revealed 'osteoblast signaling' as the most enriched pathway. PPI network led to the identification of 25 'female bone candidates'. The gene-regulatory analysis revealed RUNX2 as a key TF. ANKRD1 and RUNX2 were significantly different between osteoporosis patients and healthy controls. Sp7, Col1a1 and Pth1r correlated with multiple femur bone phenotypes in BXD mice. miR-3121-3p targeted Csf1, Ankrd1, Sp7 and Runx2. β-estradiol treatment markedly increased the expression of these candidates in mouse osteoblast. Our study revealed that Esr1-correlated genes Ankrd1, Runx2, Csf1 and Sp7 may play important roles in female bone development.

Parkinson's disease patients combined with constipation tend to have higher serum expression of microRNA 29c, prominent neuropsychiatric disorders, possible RBD conversion, and a substandard quality of life

INTRODUCTION

The symptom of constipation has been confirmed as an early diagnose criteria for Parkinson's disease (PD). Furthermore, evidences suggest that pathogenesis of PD initiates in gut, rather than brain. If so, identifying biomarkers for constipation in PD might have potentials to assist early diagnosis and initial treatment.

METHOD

We first identified that microRNA 29c (miR-29c) was dysregulated both in PD and constipation patients through bioinformatics analysis. Then, serological analysis of the expression of miR-29c in 67 PD patients with constipation (PD-C), 51 PD patients without constipation (PD-NC), and 50 healthy controls (HC) was carried out by qPCR. Demographic and clinical features were also compared. Patients in PD-C group were further classified into two groups: those with prodromal stage constipation (PD-C-Pro) (n = 36) and those with clinical stage constipation (PD-C-Clinic) (n = 31), to explore their different characteristics.

RESULTS

The levels of miR-29c in PD-C group were higher than that in PD-NC group, both higher than HC group. PD-C-Pro group's miR-29c levels were statistically higher compared with PD-C-Clinic group's. What is more, PD-C group had higher scores of MDS-UPDRS-I, NMSS, NMSS3, NMSS4, NMSS6, NMSS9, SCOPA-AUT, HAMD, HAMA, RBDSQ, CSS, and PACQOL compared with PD-NC party. Relative to the PD-C-Clinic, patients in PD-C-Pro group had higher MDS-UPDRS-I, NMSS, NMSS3, HAMD, and HAMA scores, and were more likely to have RBD.

CONCLUSION

Our results indicated that miR-29c seems to be an underlying cause for developing constipation in patients with PD and PD-C identifies a group of patients with more severe non-motor impairment, prominent neuropsychiatric disorders, and possible RBD conversion as well as a substandard quality of life. We further confirmed that there is a close relationship between symptoms representing the same pathological origin, especially constipation and RBD.

Bioinformatics analysis of hedgehog interacting protein in colorectal cancer: a study based on GEO data and TCGA data

Colorectal cancer (CRC) is one of the most prevalent malignancies worldwide. Hedgehog Interacting Protein (HHIP) is evolutionarily conserved protein, which are important morphogens for a wide range of developmental processes,However, the specific role and mechanism of HHIP in CRC remains not fully understood. In this study, we first performed pan-cancer analysis for HHIP's expression via The Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx) data and found that HHIP might be a potential anti-oncogene for CRC. Subsequently, non-coding RNAs (ncRNAs) contributing to down-regulated HHIP expression were identified through a combination of a series of in silico analyses, including expression and correlation analysis. Finally, the LINC02381/miR-577 complex was identified as the top potential upstream regulator of HHIP in CRC. In addition, HHIP expression level was significantly correlated with tumor immune cell infiltration, biomarkers of immune cells, and immune checkpoint expression. Overall, our findings clarified ncRNAs-mediated down-regulation of HHIP which was associated with poor prognosis and tumor immune infiltration in CRC.

ncRNA-mediated low expression of P2RY14 correlates with poor prognosis and tumor immune infiltration in ovarian carcinoma

Background

Ovarian cancer (OV) has been puzzling clinicians because of its poor prognosis. More and more evidence show that the G protein coupled receptor P2RY14 plays a key role in the initiation and progression of various types of human cancer. The purpose of our study is to explore the correlation between P2RY14 and the prognosis of ovarian cancer patients and the relevant mechanism.

Methods

First, the differentially expressed gene P2RY14 was screened from The Cancer Genome Atlas (TCGA) database. Explored possible P2RY14 related miRNAs and lncRNAs through multiple public databases, predicted and analyzed the expression level of candidate miRNAs and candidate lncRNAs that can bind to candidate miRNAs in OV through StarBase database. The TIMER database was used to comprehensively analyze the expression of tumor infiltrating immune cells, and to analyze the correlation between the expression level of P2RY14 and the level of immune cell infiltration in OV or the expression level of immune checkpoints.

Results

Patients with P2RY14 overexpression had better overall survival (OS) and progression-free interval (PFI). In the Targetscan database, 22 upstream miRNAs that may bind to P2RY14 were predicted. According to the regulatory network constructed by the Cytoscape software, correlation analysis and the role of miRNAs in the prognosis of OV, we first determined that the candidate miRNAs were miR-34c-5p. Then, we predicted the upstream lncRNAs of miR-34c-5p in the StarBase database, the expression level of these lncRNAs in OV in the Gene Expression Profiling Interactive Analysis (GEPIA) database, and the role in prognosis. We determined that LINC00665 is the most potential lncRNA upstream of ovarian cancer miRNA (hsa-miR-34c-5p)-P2RY14. Then, we analyzed the results in the Timer database, suggesting that P2RY14 expression was positively correlated with CD8⁺T Cell, CD4⁺T Cell, Macrophage, Neutral and Dendritic cells, and negatively correlated with B cells. Meanwhile, P2RY14 was positively correlated with CD274 and PDCD1.

Conclusions

P2RY14 can be used as a new predictive biomarker of ovarian cancer. Intervention of P2RY14 can affect the prognosis of ovarian cancer by affecting LINC00665-miR-34c-5p-P2RY14 axis. These findings provide a potential target for the development of anti-cancer strategies for ovarian cancer.

MACC1 Correlates with Tumor Progression and Immune Cell Infiltration of Colon Adenocarcinoma and is Regulated by the lncRNA ZFAS1/miR-642a-5p Axis

Colon adenocarcinoma (COAD) is the most common pathologic type of colon cancer. Metastasis is responsible for the high mortality rate of patients with COAD. The gene, metastasis-associated in colon cancer 1 (MACC1), is a biomarker predictive of both metastatic and metastasis-free survival in patients with colon cancer and other solid tumors. However, the underlying mechanism by which MACC1 affect COAD progression and metastasis remains unknown. In this study, we analyzed the expression level and prognostic value of MACC1, as well as their correlation, in patients with various types of cancer included in The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. High MACC1 expression was found to be significantly associated with poor prognosis in patients with COAD. Analysis of the potential upstream miRNA of MACC1 showed that miR-642a-5p was downregulated in COAD and was negatively correlated with MACC1 expression. Analysis of the upstream regulators of miR-642a-5p showed that the long non-coding RNA (lncRNA) ZFAS1was the most likely upstream regulator of miR-642a-5p. In addition, the expression of MACC1 correlated positively with tumor immune cell infiltration, as well as with the levels of biomarkers of five kinds of immune cells. In summary, these findings suggest that MACC1 contributes to COAD progression and immune cell infiltration via the ZFAS1/miR-642a-5p/MACC1 axis.

Two novel lncRNAs AF111167.2 and AL162377.1 targeting miR-21-5p mediated down expression of SYDE2 correlates with poor prognosis and tumor immune infiltration of ccRCC

Advanced clear cell Renal Cell Carcinoma (ccRCC) is notoriously known for its poor prognosis. Synapse defective protein 1 homolog 2 encoded by the SYDE2 gene is a Rho GTPase-activating protein whose functional tumorigenic significance is still unclear. Recent pan-cancer analysis using the Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) data showed the potential tumor-suppressing effects of SYDE2 in ccRCC. Subsequently, the TCGA, GTEx data, and human protein atlas were employed to assess the correlation between the SYDE2 expression, clinical data, and overall survival (OS) in ccRCC patients. Furthermore, microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) contributing to SYDE2 down expression were identified by expression, relationship, and survival analysis. Eventually, two novel lncRNAs, AL162377.1 and AF111167.2, targeting the miR-21-5p axis, were identified in the SYDE2 upstream non-coding RNAs (ncRNAs)-related pathway in ccRCC. The expression level of SYDE2 highly depends on the tumor immune cell infiltration and immune checkpoint expression. In summary, these data demonstrated that lncRNAs/miRNAs-mediated down-regulation of SYDE2 is related to the tumor immune infiltration. Hence, giving an insight into the prognosis of ccRCC.

NcRNA-mediated upregulation of CAMK2N1 is associated with poor prognosis and tumor immune infiltration of gastric cancer

Gastric cancer (GC) is still notorious for its poor prognosis and aggressive characteristics. Though great developments have been made in diagnosis and therapy for GC, the prognosis of patient is still perishing. In this study, differentially expressed genes (DEGs) in GC were first screened using three Gene Expression Omnibus (GEO) datasets (GSE13911, GSE29998, and GSE26899). Second, The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) data were used to validate expression of these DEGs and perform survival analysis. We selected seven candidate genes (CAMK2N1, OLFML2B, AKR7A3, CYP4X1, FMO5, MT1H, and MT1X) to carry out the next analysis. To construct the ceRNA network, we screened the most potential upstream ncRNAs of the candidate genes. A series of bioinformatics analyses, including expression analysis, correlation analysis, and survival analysis, revealed that the SNHG10–hsa-miR-378a-3p might be the most potential regulatory axis in GC. Then, the expression of CAMK2N1, miR-378a-3p, and SNHG10 was verified in GC cell lines (GES-1, MGC-803, BGC-823, HGC-27, MKN-45, and AGS) by qRT-PCR and Western blotting. We found that SNHG10 and CAMK2N1 were highly expressed in gastric cancer lines, and the miR-378a-3p was lowly expressed in BGC-823, HGC-27, and MKN-45. Furthermore, CAMK2N1 levels were significantly negatively associated with tumor immune cell infiltration, biomarkers of immune cells, and immune checkpoint expression. In summary, our results suggest that the ncRNA-mediated high expression of CAMK2N1 is associated with poor prognosis and tumor immune infiltration of GC.

A Comprehensive Analysis for Expression, Diagnosis, and Prognosis of m5C Regulator in Breast Cancer and Its ncRNA–mRNA Regulatory Mechanism

Recent studies have well demonstrated that 5-methylcytosine (m⁵C) regulators play pivotal roles in pathological conditions, including cancer. This study first tried to identify potential 5-methylcytosine (m⁵C) regulators in breast cancer by combination of expression, diagnosis, and survival analyses, and then established an ncRNA–mRNA network accounting for m⁵C regulators’ roles in breast cancer. Among 13 m⁵C regulators, DNMT3B and ALYREF were significantly upregulated in breast cancer and their high expression indicated unfavorable prognosis. Both DNMT3B and ALYREF possessed the statistical abilities to distinguish breast cancer from normal breast samples. Moreover, five potential upstream miRNAs (let-7b-5p, miR-195-5p, miR-29a-3p, miR-26a-5p, and miR-26b-5p) of m⁵C regulators could not only serve as independent prognostic predictors but also together made up a promising miRNA prognostic signature in breast cancer. Next, upstream potential lncRNAs of the five miRNAs were predicted and analyzed. Pathway enrichment analysis revealed that the target genes of these miRNAs were markedly enriched in some cancer-related pathways, and further investigation indicated VEGFA and EZH2 were found to be the most potential target genes in the m⁵C regulators-related ncRNA–mRNA network in breast cancer. These findings comprehensively provided key clues for developing m⁵C regulators-related effective therapeutic targets and promising diagnostic biomarkers in breast cancer.

SQLE, A Key Enzyme in Cholesterol Metabolism, Correlates With Tumor Immune Infiltration and Immunotherapy Outcome of Pancreatic Adenocarcinoma

Background

Pancreatic adenocarcinoma (PAAD) is a treatment-refractory cancer with poor prognosis. Accumulating evidence suggests that squalene epoxidase (SQLE) plays a pivotal role in the development and progression of several cancer types in humans. However, the function and underlying mechanism of SQLE in PAAD remain unclear.

Methods

SQLE expression data were downloaded from The Cancer Genome Atlas and the Genotype-Tissue Expression database. SQLE alterations were demonstrated based on the cBioPortal database. The upstream miRNAs regulating SQLE expression were predicted using starBase. The function of miRNA was validated by Western blotting and cell proliferation assay. The relationship between SQLE expression and biomarkers of the tumor immune microenvironment (TME) was analyzed using the TIMER and TISIDB databases. The correlation between SQLE and immunotherapy outcomes was assessed using Tumor Immune Dysfunction and Exclusion. The log-rank test was performed to compare prognosis between the high and low SQLE groups.

Results

We demonstrated a potential oncogenic role of SQLE. SQLE expression was upregulated in PAAD, and it predicted poor disease-free survival (DFS) and overall survival (OS) in patients with PAAD. "Amplification" was the dominant type of SQLE alteration. In addition, this alteration was closely associated with the OS, disease-specific survival, DFS, and progression-free survival of patients with PAAD. Subsequently, hsa-miR-363-3p was recognized as a critical microRNA regulating SQLE expression and thereby influencing PAAD patient outcome. In vitro experiments suggested that miR-363-3p could knock down the expression of SQLE and inhibit the proliferation of PANC-1. SQLE was significantly associated with tumor immune cell infiltration, immune checkpoints (including PD-1 and CTLA-4), and biomarkers of the TME. KEGG and GO analyses indicated that cholesterol metabolism-associated RNA functions are implicated in the mechanisms of SQLE. SQLE was inversely associated with cytotoxic lymphocytes and predicted immunotherapy outcomes.

Conclusions

Collectively, our results indicate that cholesterol metabolism-related overexpression of SQLE is strongly correlated with tumor immune infiltration and immunotherapy outcomes in patients with PAAD.

NcRNAs-mediated P2RX1 expression correlates with clinical outcomes and immune infiltration in patients with breast invasive carcinoma

The development of novel treatments for breast invasive carcinoma (BC) has been stagnant. P2RX1, a member of the purinergic receptor family, has been found to have a prognostic impact in several tumors. Therefore, we analyzed the expression pattern of P2RX1 in pan-cancers including BC and its impact on survival and found that the expression level of P2RX1 was lower in BC compared with para-cancerous tissues, and higher P2RX1 expression indicated better prognoses. But real-time quantitative reverse transcription PCR (RT-qPCR) and Western blot detected that the P2RX1 expression in normal mammary epithelial cells was lower than that in tumor cells. Then we comprehensively analyzed the regulatory mechanism and protein-protein interaction network, and found that P2RX1 was significantly positively linked with immune cell infiltration and immune checkpoints.

Multiple types of noncoding RNA are involved in potential modulation of PTTG1's expression and function in breast cancer

Breast cancer is a malignant type with morbidity ranking the first of women globally. As widely acknowledged, there exist close links between ncRNA-mRNA axis and breast cancer. In this study, we first overviewed expression and prognostic values of pituitary tumor transforming gene (PTTGs) in breast cancer. Next, two binding miRNAs (miR-186-5p and miR-655-3p) of PTTG1 in breast cancer were identified. Subsequently, several potential upstream ncRNAs of PTTG1-miR-186-5p/miR-655-3p axis in breast cancer were successively screened out, consisting of 11 lncRNAs, 17 circRNAs and 12 pseudogene-derived RNAs. Enrichment analysis for downstream target genes of PTTG1-miR-186-5p/miR-655-3p axis revealed that this axis is associated with TGF-beta signaling and MAPK signaling pathways. Further investigation demonstrated AURKA was one of the most key hub genes. Collectively, we established a potential PTTG1-related ncRNA-mRNA regulatory network in breast cancer.

COMMD2 Upregulation Mediated by an ncRNA Axis Correlates With an Unfavorable Prognosis and Tumor Immune Infiltration in Liver Hepatocellular Carcinoma

Liver hepatocellular carcinoma (LIHC) seriously endangers the health and quality of life of individuals worldwide. Increasing evidence has underscored that the copper metabolism MURR1 domain (COMMD) family plays important roles in tumorigenesis. However, the specific role, biological function, mechanism and prognostic value of COMMD2 and its correlation with immune cell infiltration in LIHC remain unknown. In this study, we first determined the expression and prognostic potential of COMMD2 in human tumors using The Cancer Genome Atlas (TCGA) data and identified COMMD2 as a potential oncogene in LIHC. High COMMD2 expression was associated with pathological tumor stage and metastasis. Subsequently, noncoding RNAs (ncRNAs) upregulating COMMD2 expression were identified by performing expression, correlation, and survival analyses in combination. The CRNDE/LINC00511/SNHG17/HCG18-miR-29c-3p axis was identified as the most likely ncRNA-associated pathway upstream of COMMD2 in LIHC. Next, the expression profiles of COMMD2 and ncRNAs were validated in LIHC tissues and adjacent normal tissues. Furthermore, COMMD2 was significantly positively correlated with tumor immune cell infiltration, immune cell biomarkers, and immune checkpoint molecule expression. Importantly, COMMD2 potentially influenced prognosis by regulating immune cell infiltration in LIHC. Finally, COMMD2 was knocked down in LIHC cell lines using siRNAs for functional assays in vitro, resulting in suppressed cell proliferation and migration. In summary, our findings showed that the ncRNA-mediated upregulation of COMMD2 was associated with an unfavorable prognosis correlated with immune cell infiltration in LIHC.

Bioinformatics Analysis Combined With Experiments Predicts PUDP as a Potential Prognostic Biomarker for Hepatocellular Carcinoma Through Its Interaction With Tumor Microenvironment

Hepatocellular carcinoma (HCC) is one of the deadliest tumors in the world and is notorious for poor prognosis. There is mounting evidence that pseudouridine performs key functions in the initiation and progression of several cancers. A previous study demonstrated that Pseudouridine 5’-phosphatase (PUDP) may be a novel prognostic biomarker in colorectal cancer. However, in the past, we have paid little attention to PUDP and we are still not clear about its function and role in cancer. In this study, a pan-cancer analysis of PUDP expression and prognosis was performed firstly using The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) data and we found that PUDP may be a potential oncogene for HCC. Then the most potential upstream microRNA contributing to PUDP was identified as let-7c-5p through expression analysis, correlation analysis, and survival analysis. Subsequently, the result of single cell RNA sequencing (scRNA-seq) demonstrated that PUDP was significantly highly expressed on malignant cells. In addition, there are significantly positive correlations between PUDP and tumor immune cell infiltration, biomarkers of immune cells, and immune checkpoint expression, especially with tumor-promoting immune cells such as T cell regulatory (Treg), Myeloid-derived suppressor cell (MDSC), cancer-associated fibroblast (CAF). Moreover, we found the methylation level of three loci was positively correlated with PUDP expression and four loci were negatively correlated. 15 pairs of HCC and normal adjacent tissues from HCC patients who were treated at our center were used to verify the results of the bioinformatics analysis and the results of experiments are similar to the bioinformatics analysis. Our study demonstrated that HCC patients with high PUDP expression are less likely to benefit from immunotherapy, and in addition, we explored the relationship between PUDP and anticancer drugs. Finally, we explored the clinical relevance of PUDP, identified PUDP as an independent risk factor for HCC patients and constructed a prognostic model, used International Cancer Genome Consortium (ICGC) data to do external validation. Collectively, our study demonstrated that high expression of PUDP suggested a poor prognosis and low response to immunotherapy, providing new insight into the treatment and prognosis of HCC.

FAT10 is a Prognostic Biomarker and Correlated With Immune Infiltrates in Skin Cutaneous Melanoma

Background: Skin Cutaneous Melanoma (SKCM) is the deadliest cutaneous neoplasm. Previous studies have proposed ubiquitin-like protein FAT10 plays key roles in the initiation and progression of several types of human cancer, but little is known about the interrelation between FAT10 gene expression, tumor immunity, and prognosis of patients with SKCM.

Methods: Here, we first performed pan-cancer analysis for FAT10’s expression and prognosis using the Cancer Genome Atlas and the Genotype-Tissue Expression data. Subsequently, we investigated the mRNA expression level, prognostic value, and gene-gene interaction network of FAT10 in SKCM using the Oncomine databases, GEPIA, TIMER, UALCAN, and starBase. The relationship between FAT10 expression and tumor immune invasion was studied by using the TIMER database. Additionally, the expression and functional status of FAT10 in SKCM were evaluated by the single-cell RNA sequencing and CancerSEA databases.

Results: In this study, we found that FAT10 expression was increased in SKCM and was correlated with a better survival rate in patients with SKCM. Moreover, we identified FAT10 level was significantly positively associated with immune infiltrates, biomarkers of immune cells, and immune checkpoint expression, and negatively correlated with tumor cell invasion and DNA damage, indicating that increased FAT10 expression in SKCM was a favorable response to immune checkpoint inhibitors.

Conclusion: Our findings suggest that upregulation of FAT10 correlated with better prognosis and tumor immune infiltration in SKCM.

Prediction of Mechanosensitive Genes in Vascular Endothelial Cells Under High Wall Shear Stress

Objective: The vulnerability of atherosclerotic plaques is among the leading cause of ischemic stroke. High wall shear stress (WSS) promotes the instability of atherosclerotic plaques by directly imparting mechanical stimuli, but the specific mechanisms remain unclear. We speculate that modulation of mechanosensitive genes may play a vital role in accelerating the development of plaques. The purpose of this study was to find mechanosensitive genes in vascular endothelial cells (ECs) through combining microarray data with bioinformatics technology and further explore the underlying dynamics–related mechanisms that cause the progression and destabilization of atherosclerotic plaques.

Methods: Microarray data sets for human vascular ECs under high and normal WSS were retrieved from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified through the R language. The performance of enrichment analysis and protein–protein interaction (PPI) network presented the biological function and signaling pathways of the DEGs. Hub genes were identified based on the PPI network and validated by GEO data sets. Predicted transcription factor (TF) genes and miRNAs interaction with potential mechanosensitive genes were identified by NetworkAnalyst.

Results: A total of 260 DEGs, 121 upregulated and 139 downregulated genes, were screened between high and normal WSS from GSE23289. A total of 10 hub genes and four cluster modules were filtered out based on the PPI network. The enrichment analysis showed that the biological functions of the hub genes were mainly involved in responses to unfolded protein and topologically incorrect protein, and t to endoplasmic reticulum stress. The significant pathways associated with the hub genes were those of protein processing in the endoplasmic reticulum, antigen processing, and presentation. Three out of the 10 hub genes, namely, activated transcription factor 3 (ATF3), heat shock protein family A (Hsp70) member 6 (HSPA6), and dual specificity phosphatase 1 (DUSP1, also known as CL100, HVH1, MKP-1, PTPN10), were verified in GSE13712. The expression of DUSP1 was higher in the senescent cell under high WSS than that of the young cell. The TF–miRNA–mechanosensitive gene coregulatory network was constructed.

Conclusion: In this work, we identified three hub genes, ATF3, HSPA6, and DUSP1, as the potential mechanosensitive genes in the human blood vessels. DUSP1 was confirmed to be associated with the senescence of vascular ECs. Therefore, these three mechanosensitive genes may have emerged as potential novel targets for the prediction and prevention of ischemic stroke. Furthermore, the TF–miRNA–mechanosensitive genes coregulatory network reveals an underlying regulatory mechanism and the pathways to control disease progression.

Multi-omics-based identification of atopic dermatitis target genes and their potential associations with metabolites and miRNAs

Atopic dermatitis (AD), or atopic eczema, is one of the most common inflammatory skin diseases with up to 10% prevalence in adults, and approximately 15-20% in children in industrialized countries. As a result, there is an unmet need for faster, safer, and effective treatments for AD. AD pathogenesis represents a complex interplay between multiple factors, such as environmental factors or stimuli, genetic factors, immune dysfunctions. However, although multi-omics label studies have been very useful in understanding the pathophysiological mechanisms of AD and its clinical manifestations, there have been very few studies that integrate different labels of omics data. Here, we attempted to integrate gene expression and metabolomics datasets from multiple different publicly available AD cohort datasets and conduct an integrated systems-level AD analysis. We used four different GEO transcriptome data sets and, by applying an elastic net machine learning algorithm, identified robust hub genes that can be used as signatures, for example, H2AFX, MCM7, ESR1 and SF3A2. Moreover, we investigated potential associations of those genes by applying a pathway-based approach over metabolomics and miRNA datasets. Our results revealed potential novel associations between fatty acids and peroxisomal lipid metabolism pathways, as well as with several microRNAs.

MUC14-Related ncRNA-mRNA Network in Breast Cancer

Abstract: Background Growing evidences have showed that mucins (MUCs) are linked to occurrence and progression of human cancers. However, a comprehensive study regarding the expression, diagnosis, prognosis and mechanism of MUCs in breast cancer remains absent. Methods: A series of in silico analyses were employed in this study. Results: After performing comprehensive analysis for MUCs, MUC14 was identified as the most potential regulator in breast cancer, with downregulated expression in both mRNA and protein levels and significant diagnostic and prognostic values in breast cancer. Mechanistic exploration revealed that a potential ncRNA-mRNA axis, involving LINC01128/LINC01140/SGMS1-AS1/LINC00667-miR-137/miR-429-BCL2, might be partially responsible for MUC14′s functions in breast cancer. Conclusions: Collectively, our study elucidated a key role of MUC14 in breast cancer and also provided some clues for explanation of the molecular action mechanism of MUC14 in breast cancer.

MiR-192-5p inhibits proliferation, migration, and invasion in papillary thyroid carcinoma cells by regulation of SH3RF3

Background: The decreased level of miR-192-5p has been reported in several kinds of cancers, including bladder, colon, ovarian, and non-small cell lung cancer. However, the expression and function of miR-192-5p in papillary thyroid carcinoma/cancer (PTC) remains unknown.

Objective: The present study aimed to explore the function and underlying mechanism of miR-192-5p in PTC development.

Methods: PTC tissues and relative normal controls from PTC patients were collected. qRT-PCR analysis was performed to measure miR-192-5p and SH3RF3 mRNA level in PTC tissues and cell lines. CCK-8 method and FCM assay were used to test cell proliferation and apoptosis in TPC-1 cells, respectively. The abilities of cell migration and invasion were detected by wound healing and transwell assays, respectively. The protein expression was evaluated by Western blot. The interaction between miR-192-5p and Src homology 3 (SH3) domain containing ring finger 3 (SH3RF3) were confirmed by dual-luciferase reporter assay.

Results: MiR-192-5p level was obviously decreased in PTC tissues and cell lines. Overexpression of miR-192-5p suppressed proliferation, migration, invasion, and EMT process, while induced apoptosis in TPC-1 cells. In addition, miR-192-5p negatively modulated SH3RF3 expression by binding to its 3′-untranslated region (3′UTR). Silencing SH3RF3 inhibited the migration, invasion, and EMT of TPC-1 cells. In the meantime, matrine, an alkaloid extracted from herb, exerted its anti-cancer effects in PTC cells dependent on increase in miR-192-5p expression and decrease in SH3RF3 expression.

Conclusion: We firstly declared that miR-192-5p played a tumor suppressive role in PTC via targeting SH3RF3. Moreover, matrine exerted its anti-cancer effects in PTC via regulating miR-192-5p/SH3RF3 pathway.

ncRNAs-mediated high expression of SEMA3F correlates with poor prognosis and tumor immune infiltration of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is notorious for its poor prognosis. Increasing evidence has demonstrated that semaphorin 3F (SEMA3F) plays key roles in initiation and progression of several types of human cancer. However, the specific role and mechanism of SEMA3F in HCC remains not fully determined. In this study, we first performed pan-cancer analysis for SEMA3F's expression and prognosis using The Cancer Genome Atlas (TCGA) and The Genotype-Tissue Expression (GTEx) data and found that SEMA3F might be a potential oncogene in HCC. Subsequently, noncoding RNAs (ncRNAs) contributing to SEMA3F overexpression were identified by a combination of a series of in silico analyses, including expression analysis, correlation analysis, and survival analysis. Finally, the TMPO-AS1/SNHG16-let-7c-5p axis was identified as the most potential upstream ncRNA-related pathway of SEMA3F in HCC. Moreover, SEMA3F level was significantly positively associated with tumor immune cell infiltration, biomarkers of immune cells, and immune checkpoint expression. Collectively, our findings elucidated that ncRNAs-mediated upregulation of SEMA3F correlated with poor prognosis and tumor immune infiltration in HCC.

Construction of a lncRNA/pseudogene-hsa-miR-30d-5p-GJA1 regulatory network related to metastasis of pancreatic cancer

Pancreatic cancer, the most lethal malignant tumor, is notorious for its poor prognosis and metastatic potential. Non-coding RNAs (ncRNAs) are reported to play key roles in cancer metastasis. In this study, miRNA and gene expression profiles between metastatic pancreatic cancer cell M8 and its parental cell BxPC.3 were determined. Using differential expression analysis, survival analysis, target gene prediction, pathway enrichment analysis, intersection analysis and correlation analysis, hsa-miR-30d-5p/GJA1 axis was identified as the most potential pathway involved in metastasis of pancreatic cancer. Subsequently, two upstream lncRNAs (HELLPAR and OIP-AS1) and four upstream pseudogenes (AC093616.1, AC009951.1, TMEM183B and PABPC1P4) of hsa-miR-30d-5p/GJA1 axis were predicted and were then identified via assessment of RNA-RNA expression relationship. Furthermore, CTNNA1, CTNNB1 and CTNND1 were regarded as three crucial molecules to be participated in hsa-miR-30d-5p/GJA1-mediated metastatic potential in pancreatic cancer. In conclusion, we established a novel lncRNA/pseudogene-hsa-miR-30d-5p-GJA1 regulatory network linked to metastasis of pancreatic cancer.