Breast cancer patients’ outcomes have improved dramatically in recent years, but relapses and poor prognosis remain common due to its aggressiveness and heterogeneity. The development of reliable biomarkers is still needed for predicting prognosis and treatment effectiveness. Recently, a growing body of research suggests that pseudouridine synthases contribute to the development of many cancers, but their contribution to breast cancer remains largely unknown. Using an integrative analysis, we selected pseudouridine synthase1(PUS1) as the candidate biomarker. A tissue microarray of 131 breast cancer patients was then utilized to determine the clinical significance and prognostic value of PUS1. RNA sequencing analysis was conducted to identify downstream genes that differ between control and PUS1 knockdown cells. The effect of PUS1 on phenotypes of cells was assessed using cell proliferation, colony formation, and transwell invasion assays. We found that breast tumors overexpressed PUS1 compared with paired normal tissues. PUS1 expression was positively correlated with triple-negative breast cancer (TNBC) status (P= 0.020) and tumor grade (P <0.0001), but not with age (P= 0.736), tumor size (P= 0.608), lymph node (P= 0.742), oestrogen receptor (ER) (P= 0.162), progesterone receptor (PR) (P= 0.901), human epidermal growth factor receptor 2 (HER2) (P= 0.608) or tumor stage (P= 0.411). Comparatively, patients with high PUS1 levels had shorter overall survival time (P=0.0001) and relapse-free survival time (P = 0.0093). A univariate and multivariate survival analysis suggested that the overall survival of patients was independently influenced by the PUS1 score (Univariate Cox P <0.0001, HR=5.176, 95% CI =2.420-11.07; Multivariate Cox P = 0.001, HR = 5.291, 95% CI =1.893-14.78). RNA sequencing data revealed the PUS1 knockdown significantly affects a series of cancer related biological process such as regulation of cell proliferation and cell migration, as well as KEGG pathways including Mitophagy and PI3K-Akt signaling. In vitro, knockdown of PUS1 significantly suppressed the proliferation and colony formation abilities of MDA-MB-231 cells and BT-549 cells. Additionally, the ability of tumor cells to invade was remarkably attenuated in low PUS1 expression groups compared with the corresponding control groups. Our results suggested that PUS1 is a novel biomarker that predicts poor outcomes in patients with breast cancer and may prove to be a promising treatment target.

Background: Transmembrane p24 trafficking protein (TMED) family members are implicated in several solid tumors, but their clinical relevance for breast cancer (BC) remains unclear. This study aimed to probe their prognostic values and relations with tumor immunity in BC.Methods: TMED family mRNA expression was assessed in five microarray datasets (GSE65212, GSE42568, GSE5364, GSE22820 and GSE45827) from Gene Expression Omnibus (GEO) database and invasive breast cancer (BRCA) cohort from The Cancer Genome Atlas (TCGA). Receiver operating characteristic (ROC) curve was performed to determine the predictive values of filtered members of the TMED family. The protein expressions of screen genes were validated by Clinical Proteomic Tumor Analysis Consortium (CPTAC) data from University of ALabama at Birmingham CANcer data analysis portal (UALCAN) and detected in the clinical specimens by western blot assay. Clinicopathologic variables were analyzed with bc-GenExMiner, and patient prognostic data were obtained with Kaplan-Meier Plotter.In vitro wound healing and invasion assays were performed on siRNA-transfected BC cell lines. TIMER 2.0, SangerBox, and ImmPort were used to evaluate tumor immune infiltration, immune checkpoints, and other immune-related genes. CbioPortal, Metascape, Expression2kinases, and LinkedOmics were used to explore gene regulatory network.Results: BC tissues expressed TMED2/3/4/9 at a higher level than normal tissues, providing diagnostic potential. All the areas under the ROC curve for TMED2/3/4/9 were more than 0.7. TMED2/3/4/9 correlated with numerous clinical variables, including lymph node status, Scarff-Bloom-Richardson score (SBR), Nottingham Prognostic Index (NPI), estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER-2), and triple-negative breast cancer (TNBC) status, and their high expression predicted the poor prognosis of BC patients. TMED2/3/4/9 knockdown drastically inhibited the migratory and invasive capacities of MDA-MB-231 and HCC1937 cells. TMED2/3/4/9 expressions correlated negatively with the infiltration of tumor-suppressive immune cells such as CD8 + T cells, dendritic cells, and natural killer cells, and was inversely related to a variety of immune checkpoint genes, including programmed cell death 1 (PD-1) and cytotoxic T-lymphocyte associated protein 4 (CTLA4). A set of kinases, transcription factors, and microRNAs (miRNAs) may regulate TMED2/3/4/9 abnormalities at the genome level.