Prolyl 4-hydroxylase subunit alpha 1 (P4HA1) is a biomarker of poor prognosis in primary melanomas, and its depletion inhibits melanoma cell invasion and disrupts tumor blood vessel walls

Hydrogel-based platforms for site-specific doxorubicin release in cancer therapy

Hydrogels are promising candidates for the delivery of therapeutics in the treatment of human cancers. Regarding to the biocomaptiiblity, high drug and encapsulation efficacy and adjustable physico-chemical features, the hydrogels have been widely utilized for the delivery of chemotherapy drugs. Doxorubicin (DOX) is one of the most common chemotherapy drugs used in cancer therapy through impairing topoisomerase II function and increasing oxidative damage. However, the tumor cells have developed resistance into DOX-mediated cytotoxic impacts, requiring the delivery systems to increase internalization and anti-cancer activity of this drug. The hydrogels can deliver DOX in a sustained manner to maximize its anti-cancer activity, improving cancer elimination and reduction in side effects and drug resistance. The natural-based hydrogels such as chitosan, alginate and gelatin hydrogels have shown favourable biocompatibility and degradability in DOX delivery for tumor suppression. The hydrogels are able to co-deliver DOX with other drugs or genes to enhance drug sensitivity and mediate polychemotherapy, synergistically suppressing cancer progression. The incorporation of nanoparticles in the structure of hydrogels can improve the sustained release of DOX and enhancing intracellular internalization, accelerating DOX's cytotoxicity. Furthermore, the stimuli-responsive hydrogels including pH-, redox- and thermo-sensitive platforms are able to improve the specific release of DOX at the tumor site. The DOX-loaded hydrogels can be further employed in the clinic for the treatment of cancer patients and improving efficacy of chemotherapy.

Effects of Continuous Prenatal Low Dose Rate Irradiation on Neurobehavior, Hippocampal Cellularity, Messenger RNA and MicroRNA Expression on B6C3F1 Mice

Epidemiological, experimental, and ecological data have indicated the controversial effect of in utero chronic low dose rate (<6 mGy/h) with accumulative low (≤100 mGy) or high (>100 mGy) dose radiation exposure. Our main goal of this study was to examine if different low dose rates of chronic pre- and/or post-natal radiation exposure with accumulative high doses could induce hippocampal cellular, mRNA, and miRNA changes leading to neuropsychiatric disorders. The comprehensive mouse phenotypic traits, organ weight, pathological, and blood mRNA and miRNA changes were also studied. Using different approaches including SmithKline, Harwell, Imperial College, Royal Hospital, Phenotype Assessment (SHIRPA), neurobehavioral tests, pathological examination, immunohistochemistry, mRNA and miRNA sequencing, and real-time quantitative polymerase chain reaction (qRT-PCR) validation, we found that in prenatally irradiated (100 mGy/d for 18 days with an accumulative dose of 1.8 Gy) 1-year-old mice, no cellular changes, including immature neurons in the subgranular zone, mature neurons and glial cells in the hilus of the dentate gyrus and development of cognitive impairment, neuropsychiatric disorders, occurred. However, a significant reduction in body weight and mass index (BMI) was indicated by the SHIRPA test. A reduced exploratory behavior was shown by an open field test. Organ weights showed significant reductions in the testes, kidneys, heart, liver and epididymides with no abnormal pathology. mRNA and miRNA sequencing and qRT-PCR validation revealed the upregulation of Rubcnl and Abhd14b, and downregulation of Hspa1b, P4ha1, and Banp genes in both the hippocampus and blood of mice prenatally irradiated with 100 mGy/d. Meanwhile, downregulation of miR-448-3p and miR1298-5p in the hippocampus, miR-320-3p, miR-423-5p, miR-486b-5p, miR-486b-3p, miR-423-3p, miR-652-3p, miR-324-3p, miR-181b-5p, miR-let-7b, and miR-6904-5p in the blood was induced. The target scan revealed that Rubcnl is one of the miR-181b-5p targets in the blood. We, therefore, concluded that prenatal chronic irradiation with a low dose rate of 100 mGy/d and accumulative dose of 1.8 Gy or below might not induce significant adverse health effects on the offspring. Further study of different low dose rate radiation exposures with accumulative high doses may provide threshold doses for authorities or regulators to set new radiation safety guidelines to replace those extrapolated from acute high dose/dose rate irradiation to reduce unnecessary emergency evacuation or spending once a nuclear accident or leakage occurs.

Recent Advances and Therapeutic Implications of 2-Oxoglutarate-Dependent Dioxygenases in Ischemic Stroke

Ischemic stroke is a common disease with a high disability rate and mortality, which brings heavy pressure on families and medical insurance. Nowadays, the golden treatments for ischemic stroke in the acute phase mainly include endovascular therapy and intravenous thrombolysis. Some drugs are used to alleviate brain injury in patients with ischemic stroke, such as edaravone and 3-n-butylphthalide. However, no effective neuroprotective drug for ischemic stroke has been acknowledged. 2-Oxoglutarate-dependent dioxygenases (2OGDDs) are conserved and common dioxygenases whose activities depend on O2, Fe2+, and 2OG. Most 2OGDDs are expressed in the brain and are essential for the development and functions of the brain. Therefore, 2OGDDs likely play essential roles in ischemic brain injury. In this review, we briefly elucidate the functions of most 2OGDDs, particularly the effects of regulations of 2OGDDs on various cells in different phases after ischemic stroke. It would also provide promising potential therapeutic targets and directions of drug development for protecting the brain against ischemic injury and improving outcomes of ischemic stroke.

The dynamic transcriptomic response of the goldfish brain under chronic hypoxia

Oxygen is essential to fuel aerobic metabolism. Some species evolved mechanisms to tolerate periods of severe hypoxia and even anoxia in their environment. Among them, goldfish (Carassius auratus) are unique, in that they do not enter a comatose state under severely hypoxic conditions. There is thus significant interest in the field of comparative physiology to uncover the mechanistic basis underlying hypoxia tolerance in goldfish, with a particular focus on the brain. Taking advantage of the recently published and annotated goldfish genome, we profile the transcriptomic response of the goldfish brain under normoxic (21 kPa oxygen saturation) and, following gradual reduction, constant hypoxic conditions after 1 and 4 weeks (2.1 kPa oxygen saturation). In addition to analyzing differentially expressed protein-coding genes and enriched pathways, we also profile differentially expressed microRNAs (miRs). Using in silico approaches, we identify possible miR-mRNA relationships. Differentially expressed transcripts compared to normoxia were either common to both timepoints of hypoxia exposure (n = 174 mRNAs; n = 6 miRs), or exclusive to 1-week (n = 441 mRNAs; n = 23 miRs) or 4-week hypoxia exposure (n = 491 mRNAs; n = 34 miRs). Under chronic hypoxia, an increasing number of transcripts, including those of paralogous genes, was downregulated over time, suggesting a decrease in transcription. GO-terms related to the vascular system, oxidative stress, stress signalling, oxidoreductase activity, nucleotide- and intermediary metabolism, and mRNA posttranscriptional regulation were found to be enriched under chronic hypoxia. Known 'hypoxamiRs', such as miR-210-3p/5p, and miRs such as miR-29b-3p likely contribute to posttranscriptional regulation of these pathways under chronic hypoxia in the goldfish brain.

The role of collagen triple helix repeat containing 1 (CTHRC1) in cancer development and progression

INTRODUCTION

Collagen triple helix repeat containing 1 (CTHRC1) is a protein that has been implicated in pro-migratory pathways, arterial tissue-repair processes, and inhibition of collagen deposition via the regulation of multiple signaling cascades. Studies have also demonstrated an upregulation of CTHRC1 in multiple cancers where it has been linked to enhanced proliferation, invasion, and metastasis. However, the understanding of the exact role and mechanisms of CTHRC1 in cancer is far from complete.

AREAS COVERED

This review focuses on analyzing the role of CTHRC1 in cancer as well as its associations with clinicopathologies and cancer-related processes and signaling. We have also summarized the available literature information regarding the role of CTHRC1 in tumor microenvironment and immune signaling. Finally, we have discussed the mechanisms associated with CTHRC1 regulations, and opportunities and challenges regarding the development of CTHRC1 as a potential target for cancer management.

EXPERT OPINION

CTHRC1 is a multifaceted protein with critical roles in cancer progression and other pathological conditions. Its association with lower overall survival in various cancers, and impact on the tumor immune microenvironment make it an intriguing target for further research and potential therapeutic interventions in cancer.

Multi-omics analysis reveals cuproptosis and mitochondria-based signature for assessing prognosis and immune landscape in osteosarcoma

Background

Osteosarcoma (OSA), the most common primary mesenchymal bone tumor, is a health threat to children and adolescents with a dismal prognosis. While cuproptosis and mitochondria dysfunction have been demonstrated to exert a crucial role in tumor progression and development, the mechanisms by which they are regulated in OSA still await clarification.

Methods

Two independent OSA cohorts containing transcriptome data and clinical information were collected from public databases. The heterogeneity of OSA were evaluated by single cell RNA (scRNA) analysis. To identify a newly molecular subtype, unsupervised consensus clustering was conducted. Cox relevant regression methods were utilized to establish a prognostic gene signature. Wet lab experiments were performed to confirm the effect of model gene in OSA cells.

Results

We determined 30 distinct cell clusters and assessed OSA heterogeneity and stemness scRNA analysis. Then, univariate Cox analysis identified 24 candidate genes which were greatly associated with the prognosis of OSA. Based on these prognostic genes, we obtained two molecular subgroups. After conducting step Cox regression, three model genes were selected to construct a signature showing a favorable performance to forecast clinical outcome. Our proposed signature could also evaluate the response to chemotherapy and immunotherapy of OSA cases.

Conclusion

We generated a novel risk model based on cuproptosis and mitochondria-related genes in OSA with powerful predictive ability in prognosis and immune landscape.

Identification and validation of INHBE and P4HA1 as hub genes in non-alcoholic fatty liver disease

PURPOSE

Non-alcoholic fatty liver disease (NAFLD) is currently the most prevalent type of liver disease and a worldwide disease threatening human health. This study aims to identify the novel diagnostic biomarkers of NAFLD by comprehensive bioinformatics and machine learning, and to validate our results in hepatocyte and animal models.

METHODS

We used Gene Expression Omnibus (GEO) databases on NAFLD patients for differential gene expression analyses. Intersections were taken with genes from the key modules of WGCNA and differentially expressed genes (DEGs). Machine learning algorithms like LASSO regression analysis, SVM-RFE, and RandomForest were used to screen hub genes. In addition, a nomogram model and calibration curves were built in order to forecast the probability of NAFLD occurrence. Then, the relationship between hub genes and immune cells was verified using Spearman analysis. Finally, we further verified the expression of key genes by constructing a steatosis hepatocyte model and animal model.

RESULTS

Key genes (INHBE and P4HA1) were identified by comprehensive bioinformatics analysis and machine learning. INHBE and P4HA1 were up-regulated and down-regulated in the steatosis hepatocyte model, respectively. Animal experiments also showed that INHBE was up-regulated in the liver of mice fed with high fat diet (HFD).

CONCLUSION

INHBE and P4HA1 are the hub genes of NAFLD. Our findings may contribute to a greater understanding of the occurrence and development of NAFLD and provide potential biomarkers and possible therapeutic targets for future clinical diagnosis and treatment.

Mitochondrial-related genes markers that predict survival in patients with head and neck squamous cell carcinoma affect immunomodulation through hypoxia, glycolysis, and angiogenesis pathways

Mitochondria play a crucial role in the occurrence and development of tumors. We used mitochondria-related genes for consistent clustering to identify three stable molecular subtypes of head and neck squamous cell carcinoma (HNSCC) with different prognoses, mutations, and immune characteristics. Significant differences were observed in clinical characteristics, immune microenvironment, immune cell infiltration, and immune cell scores. TP53 was the most significantly mutated; cell cycle-related pathways and tumorigenesis-related pathways were activated in different subtypes. Risk modeling was conducted using a multifactor stepwise regression method, and nine genes were identified as mitochondria-related genes affecting prognosis (DKK1, EFNB2, ITGA5, AREG, EPHX3, CHGB, P4HA1, CCND1, and JCHAIN). Risk score calculations revealed significant differences in prognosis, immune cell scores, immune cell infiltration, and responses to conventional chemotherapy drugs. Glycolysis, angiogenesis, hypoxia, and tumor-related pathways were positively correlated with the RiskScore. Clinical samples were subjected to qPCR to validate the results. In this work, we constructed a prognostic model based on the mitochondrial correlation score, which well reflects the risk and positive factors for the prognosis of patients with HNSCC. This model can be used to guide individualized adjuvant and immunotherapy in patients with HNSCC.

The Role of P4HA1 in Multiple Cancer Types and its Potential as a Target in Renal Cell Carcinoma

Background: Prolyl 4-hydroxylase subunit alpha 1 (P4HA1) provides the majority of the catalytic site of the active P4H enzyme. Emerging evidence has revealed that P4HA1 participates in the initiation and development of several malignant tumors. However, a pan-cancer analysis of P4HA1 has not been performed.

Methods: In this study, we carried out an in-depth analysis of the expression patterns and prognostic value of P4HA1 using the datasets of The Cancer Genome Atlas (TCGA) and Kaplan–Meier Plotter. Genomic and epigenetic alterations of P4HA1 and the correlation of P4HA1 with DNA methylation in different cancers were also analyzed across multiple databases. In addition, the purity-adjusted partial Spearman’s correlation test was utilized to evaluate the correlation between P4HA1 expression and immune cell infiltration. We also further explored the biological function and mechanism of P4HA1 in renal cell carcinoma (RCC).

Results: We characterized the expression profiles and prognostic values of P4HA1 in multiple cancer types. P4HA1 expression was increased in clear cell renal cell carcinoma (RCC) compared to adjacent normal tissues, and P4HA1 positively correlated with the overall survival (OS) and disease-free survival (DFS) in papillary RCC. In addition, a positive correlation between P4HA1 expression and immune cell infiltration was observed in clear cell RCC. We also identified a strong correlation between P4HA1 expression and immune checkpoint gene expression, microsatellite instability, and tumor mutation burden in chromophobe RCC. Finally, the results of in vitro experiments verified that overexpression of P4HA1 promoted the proliferation, migration, invasion, and epithelial–mesenchymal transition of RCC cells.

Conclusion: Overall, our study has suggested that P4HA1 might play a significant role in tumorigenesis in RCC and may be a prognostic biomarker and therapeutic target for several malignant tumors, including RCC.

P4HA1 Regulates CD31 via COL6A1 in the Transition of Glioblastoma Stem-Like Cells to Tumor Endothelioid Cells

Glioblastoma multiforme (GBM) is a common intracranial malignancy characterized by abundant and aberrant vasculature. The efficiency of existing antivascular treatments remains unsatisfactory. The transition of glioblastoma stem-like cells (GSCs) into tumor endothelioid cells (ECs) has been thought to cause glioma neovascularization and anti-angiogenesis tolerance, but the mechanisms regulating glioma transdifferentiation remains unclear. Our previous study found that P4HA1 regulates GSCs vascular mimicry in a hypoxic microenvironment, but the detailed molecular mechanism has not been determined. In this study, candidate protein COL6A1 was screened by mass spectrometry. In vitro experiments show that P4HA1 regulates the expression of CD31 via COL6A1, with the levels of expression of P4HA1, COL6A1 and the vascular endothelial molecular markers CD31 showing positive correlations in vivo assay. Altering the expression of P4HA1 in GSCs altered the expression of COL6A1 and CD31, thereby inducing glioma angiogenesis. In conclusion, this study revealed that the P4HA1/COL6A1 axis modulates the transdifferentiation process of GSCs into ECs. Interrupting this signaling axis can inhibit glioma angiogenesis, suggesting that this axis may be a novel target for antivascular therapy in patients with glioma.

Protein Expression of AEBP1, MCM4, and FABP4 Differentiate Osteogenic, Adipogenic, and Mesenchymal Stromal Stem Cells

Mesenchymal stem cells (MSCs) gain an increasing focus in the field of regenerative medicine due to their differentiation abilities into chondrocytes, adipocytes, and osteoblastic cells. However, it is apparent that the transformation processes are extremely complex and cause cellular heterogeneity. The study aimed to characterize differences between MSCs and cells after adipogenic (AD) or osteoblastic (OB) differentiation at the proteome level. Comparative proteomic profiling was performed using tandem mass spectrometry in data-independent acquisition mode. Proteins were quantified by deep neural networks in library-free mode and correlated to the Molecular Signature Database (MSigDB) hallmark gene set collections for functional annotation. We analyzed 4108 proteins across all samples, which revealed a distinct clustering between MSCs and cell differentiation states. Protein expression profiling identified activation of the Peroxisome proliferator-activated receptors (PPARs) signaling pathway after AD. In addition, two distinct protein marker panels could be defined for osteoblastic and adipocytic cell lineages. Hereby, overexpression of AEBP1 and MCM4 for OB as well as of FABP4 for AD was detected as the most promising molecular markers. Combination of deep neural network and machine-learning algorithms with data-independent mass spectrometry distinguish MSCs and cell lineages after adipogenic or osteoblastic differentiation. We identified specific proteins as the molecular basis for bone formation, which could be used for regenerative medicine in the future.

Identification of a novel glycolysis-related gene signature for predicting the survival of patients with colon adenocarcinoma

BACKGROUND

The most frequent histologic subtype of colon cancer is colon adenocarcinoma (COAD). A major problem in the diagnosis and treatment of COAD is that there is lack of new biomarkers to indicate the early stage of COAD. Compared with normally differentiated cells, the glycolytic pathways of tumor cells are more active, thus making them more adaptable to the hypoxic environment of solid tumors, which is known as the Warburg effect. Therefore, establishing a diagnostic and prognostic model based on glycolysis-related genes may provide guidance for the precise treatment of colon cancer.

METHODS

The Cancer Genome Atlas (TCGA) mRNA data were used to identify differentially expressed genes (DEGs). The glycolysis-related DEGs were identified using Gene Set Enrichment Analysis (GSEA) with HALLMARK gene sets. Combined with clinical data, we identified prognostic genes in glycolysis-related DEGs based on Cox regression analysis. Four glycolysis-related genes were identified and a predictive model was developed using univariate and multivariate Cox regression analysis. cBioPortal investigated the chromosomal variations of these genes. Following that, survival analysis and receiver operating characteristic (ROC) curve validation were carried out. The correlations between glycolysis-related gene signatures and molecular features and cancer subtypes were analyzed.

RESULTS

We discovered five genes (SPAG4, P4HA1, STC2, ENO3, and GPC1) that are associated with COAD patients' prognosis. The risk score was more accurate in predicting prognosis when based on this gene signature in COAD patients. Furthermore, multivariate Cox regression analysis demonstrated that the glycolysis-related gene signature's predictive value was independent of clinical variables.

CONCLUSION

We identified a glycolysis-related five-gene signature and developed a risk staging model potentially valuable for the clinical management of COAD patients. Our results suggest that prognostic markers based on glycolysis-related genes may be a reliable predictive tool for the prognosis of COAD patients.

Polysaccharide hydrogels regulate macrophage polarization and enhance the anti-tumor efficacy of melanoma

Chemotherapy occupies a prominent position in combination treatments of melanoma. However, the severe systemic side effects and the pro-tumorigenic microenvironment limited its therapeutic efficacy. In the present study, polysaccharide hydrogels (SCOD) were constructed by N-succinyl chitosan and oxidized dextran through Schiff-base formation to deliver doxorubicin (Dox) locally. The gelation time and mechanical properties of SCOD hydrogels could be fine-tuned by varying concentration of precursor solutions. Rheological data revealed that SCOD hydrogels possessed injectable shear-shinning property and remarkable self-healing capability. It also could be degraded by lysozyme widely present in body fluids. Moreover, SCOD hydrogels were readily loaded with Dox in precursor solutions and released drug over 1 week in a pH-dependent manner. The ability of Dox-loaded SCOD hydrogels to inhibit the growth of murine B16 and human A375 melanoma was verified by in vitro assays. Strikingly, Dox-loaded SCOD hydrogels were found to efficiently induce polarization of tumor-associated macrophages towards M1 phenotype that favors an anti-tumorigenic tumor microenvironment. Notably, in vivo experiments demonstrated that locally injected Dox-loaded SCOD hydrogels exhibited excellent anti-tumor activity against B16 melanoma, outperforming Dox at equivalent doses administrated intravenously. Therefore, the injectable and self-healing polysaccharide hydrogels are a promising strategy to improve locoregional control in melanoma.

Constructe a novel 5 hypoxia genes signature for cervical cancer

Background

Hypoxia, which affects the development, metastasis and prognosis of cancer, represents a key feature of cancer. This study describe a hypoxia risk factor model, with predicting the prognosis of cervical cancer.

Methods

Based on hypoxia pathway related genes, we divided cervical cancer samples into high and low expression groups. A cox analysis was then performed. Genes from these cervical cancer samples showing a significant impact on OS were selected for cluster analysis to obtain two subtypes. The TPM dataset of TCGA was divided into training and validation sets. For the training set, a lasso analysis was conducted as based on cox analysis of meaningful genes and a risk factor model was constructed. The constructed model was verified in internal and external data sets. Finally, RT-PCR, immunohistochemistry were used to detect the expression of relative genes or proteins and functional assays were used to evaluate the biological function of signature genes.

Results

Two molecular subtypes were obtained, Cluster2 vs Cluster1.These subtypes were obtained by clustering with a total of 149 DEGs (Differential expressed genes) being in line with this standard, of which 27 were up-regulated and 122 were down-regulated. The five genes with lambda = 0.0571 were selected to construct the model, the RiskScore = AK4*0.042 + HK2*0.021 + P4HA1*0.22 + TGFBI*0.1 + VEGFA*0.077. Further, in order to verify the signature, we used TCGA-test and GSE44001 chip datasets to test, and finally got a good risk prediction effect in those datasets. Moreover, the result of RT-PCR and immunohistochemistry demonstrated that AK4, HK2, P4HA1, TGFBI and VEGFA were all highly expressed in these cervical cancer tissue samples. The functional study shown that expression of AK4, HK2, P4HA1, TGFBI and VEGFA can regulate the proliferation, migration, and invasion ability of cervical cancer cells in vitro.

Conclusions

In summary, we developed a 5-gene signature prognostic hierarchical system based on the hypoxic pathway of cervical cancer, which is independent of clinical characteristics. And also conducted experimental verifications on these signature gene. Therefore, we propose that use of this classifier as a molecular diagnostic test can provide an effective means for evaluating the prognostic risk of cervical cancer patients, and provide potential targets for the treatment of cervical cancer patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02050-3.

Prognostic Value of Eight-Gene Signature in Head and Neck Squamous Carcinoma

Head and neck cancer (HNC) is the fifth most common cancer worldwide. In this study, we performed an integrative analysis of the discovery set and established an eight-gene signature for the prediction of prognosis in patients with head and neck squamous cell carcinoma (HNSCC). Univariate Cox analysis was used to identify prognosis-related genes (with P < 0.05) in the GSE41613, GSE65858, and TCGA-HNSC RNA-Seq datasets after data collection. We performed LASSO Cox regression analysis and identified eight genes (CBX3, GNA12, P4HA1, PLAU, PPL, RAB25, EPHX3, and HLF) with non-zero regression coefficients in TCGA-HNSC datasets. Survival analysis revealed that the overall survival (OS) of GSE41613 and GSE65858 datasets and the progression-free survival(DFS)of GSE27020 and GSE42743 datasets in the low-risk group exhibited better survival outcomes compared with the high-risk group. To verify that the eight-mRNA prognostic model was independent of other clinical features, KM survival analysis of the specific subtypes with different clinical characteristics was performed. Univariate and multivariate Cox regression analyses were used to identify three independent prognostic factors to construct a prognostic nomogram. Finally, the GSVA algorithm identified six pathways that were activated in the intersection of the TCGA-HNSC, GSE65858, and GSE41613 datasets, including early estrogen response, cholesterol homeostasis, oxidative phosphorylation, fatty acid metabolism, bile acid metabolism, and Kras signaling. However, the epithelial–mesenchymal transition pathway was inhibited at the intersection of the three datasets. In conclusion, the eight-gene prognostic signature proved to be a useful tool in the prognostic evaluation and facilitate personalized treatment of HNSCC patients.

Identification of a novel glycolysis-related gene signature for predicting the prognosis of osteosarcoma patients

Glycolysis ensures energy supply to cancer cells, thereby facilitating tumor progression. Here, we identified glycolysis-related genes that could predict the prognosis of patients with osteosarcoma. We examined 198 glycolysis-related genes that showed differential expression in metastatic and non-metastatic osteosarcoma samples in the TARGET database, and identified three genes (P4HA1, ABCB6, and STC2) for the establishment of a risk signature. Based on the signature, patients in the high-risk group had poor outcomes. An independent Gene Expression Omnibus database GSE21257 was selected as the validation cohort. Receiver operating characteristic curve analysis was performed and the accuracy of predicting the 1- and 3-year survival rates was shown by the areas under the curve. The results were 0.884 and 0.790 in the TARGET database, and 0.740 and 0.759 in the GSE21257, respectively. Furthermore, we applied ESTIMATE algorithm and performed single sample gene set enrichment analysis to compare tumor immunity between high- and low-risk groups. We found that the low-risk group had higher immune scores and immune infiltration levels than the high-risk group. Finally, we chose P4HA1 as a representative gene to verify the function of risk genes in vitro and in vivo and found that P4HA1 could promote the metastasis of osteosarcoma cells. Our study established a novel glycolysis-related risk signature that could predict the prognosis of patients with osteosarcoma.

Gene expression profiles for low-dose exposure to diethyl phthalate in rodents and humans: a translational study with implications for breast carcinogenesis

Phthalates are commonly included as ingredients in personal care products such as cosmetics, shampoos and perfumes. Diethyl phthalate (DEP) has been found to be anti-androgenic and linked with adverse reproductive effects on males, but effects on females are poorly understood. We designed an integrative and translational study to experimentally examine the effects of DEP exposure at a human-equivalent dose on the mammary transcriptome in rats and to subsequently examine the DEP gene signature in breast tissues (both pre-malignant and tumor) from a population study. In Sprague-Dawley rats treated orally with DEP from birth to adulthood, we identified a signature panel of 107 genes predominantly down-regulated by DEP exposure. Univariate analysis of this 107 DEP gene signature in pre-malignant breast tissues revealed that six genes (P4HA1, MPZL3, TMC4, PLEKHA6, CA8, AREG) were inversely associated with monoethyl phthalate (MEP; the urinary metabolite of DEP) concentration (p < 0.05) among postmenopausal women; all six genes loaded on to one of seven factors identified by factor analysis. Transcription factor enrichment analysis revealed that genes in this factor were enriched for androgen receptor binding sites. These six genes were also significantly down-regulated in pre-malignant adjacent tissues compared to the corresponding tumor tissues in pair-wise analyses (p < 0.05). Results from our translational study indicate that low level exposure to diethyl phthalate results in measurable genomic changes in breast tissue with implications in breast carcinogenesis.