Single-cell and spatial transcriptomics reveal metastasis mechanism and microenvironment remodeling of lymph node in osteosarcoma.

BACKGROUND: Osteosarcoma (OS) is the most common primary malignant bone tumor and is highly prone to metastasis. OS can metastasize to the lymph node (LN) through the lymphatics, and the metastasis of tumor cells reestablishes the immune landscape of the LN, which is conducive to the growth of tumor cells. However, the mechanism of LN metastasis of osteosarcoma and remodeling of the metastatic lymph node (MLN) microenvironment is not clear. METHODS: Single-cell RNA sequencing of 18 samples from paracancerous, primary tumor, and lymph nodes was performed. Then, new signaling axes closely related to metastasis were identified using bioinformatics, in vitro experiments, and immunohistochemistry. The mechanism of remodeling of the LN microenvironment in tumor cells was investigated by integrating single-cell and spatial transcriptomics. RESULTS: From 18 single-cell sequencing samples, we obtained 117,964 cells. The pseudotime analysis revealed that osteoblast(OB) cells may follow a differentiation path from paracancerous tissue (PC) → primary tumor (PT) → MLN or from PC → PT, during the process of LN metastasis. Next, in combination of bioinformatics, in vitro and in vivo experiments, and immunohistochemistry, we determined that ETS2/IBSP, a new signal axis, might promote LN metastasis. Finally, single-cell and spatial dissection uncovered that OS cells could reshape the microenvironment of LN by interacting with various cell components, such as myeloid, cancer-associated fibroblasts (CAFs), and NK/T cells. CONCLUSIONS: Collectively, our research revealed a new molecular mechanism of LN metastasis and clarified how OS cells influenced the LN microenvironment, which might provide new insight for blocking LN metastasis.

Single-cell RNA sequencing reveals the immune microenvironment landscape of osteosarcoma before and after chemotherapy.

Chemotherapy, a primary treatment for osteosarcoma (OS), has limited knowledge regarding its impact on tumor immune microenvironment (TIME). Here, tissues from 6 chemotherapy-naive OS patients underwent single-cell RNA sequencing (scRNA-seq) and were analyzed alongside public dataset (GSE152048) containing 7 post-chemotherapy OS tissues. CD45+ (PTPRC+) cells were used for cell clustering and annotation. Changes in immune cell composition pre- and post-chemotherapy were characterized. Totally, 28,636 high-quality CD45+ (PTPRC+) cells were extracted. Following chemotherapy, the proportions of regulatory T cells (Tregs) and activated CD8 T cells decreased, while CD8 effector T cells increased. GO analysis indicated that differentially expressed genes (DEGs) in T cells were associated with cell activation, adaptive immune response, and immune response to tumor cells. Furthermore, the proportions of plasma cells increased, while naive B cells decreased. B cell surface receptors expression was upregulated, and GO analysis revealed DEGs of B cells were mainly enriched in B cell-mediated immunity and B cell activation. Moreover, M2 polarization of macrophages was suppressed post-chemotherapy. Overall, this study elucidates chemotherapy remodels the OS TIME landscape, triggering immune heterogeneity and enhancing anti-tumor properties.

Surgical treatment of an old avulsion fracture of the ischial tuberosity and ischial ramus: A case report.

BACKGROUND: Avulsion fracture of the ischial tuberosity is a relatively clinically rare type of trauma that is mainly incurred by adolescents during competitive sports activities. According to previous literature, the most commonly involved sports are soccer, sprinting, and gymnastics, in descending order. Dance-induced avulsion fracture of the ischial tuberosity and ischial ramus is extremely clinically rare. CASE SUMMARY: A case of a neglected avulsion fracture of the ischial tuberosity and ischial ramus was diagnosed in a young female dancer who complained of pain and restricted movement of her right hip. She stated that she had suffered the injury while performing a split leap during a dance performance 9 mo prior. Eventually, she underwent surgery and obtained satisfactory treatment results. CONCLUSION: Early diagnosis of these fractures is important to ensuring early proper treatment towards a quicker recovery. For old fractures with nonunion and chronic buttock pain, surgery is a preferred therapeutic choice with good treatment outcomes.

Prediction of the active compounds and mechanism of Biochanin A in the treatment of Legg-Calvé-Perthes disease based on network pharmacology and molecular docking.

BACKGROUND: Legg-Calvé-Perthes disease is a special self-limited disease in pediatric orthopedics with a high disability rate and a long-term course, and there is still no clear and effective therapeutic drug in clinic. This study aimed to investigate the potential efficacy of biochanin A, a kind of oxygen-methylated isoflavone compound, in treating Perthes disease based on network pharmacology, molecular docking and in vitro experiments. METHODS: IL-6 was used to stimulate human umbilical vein endothelial cells to construct endothelial cell dysfunction model. We demonstrated whether biochanin A could alleviate endothelial dysfunction through CCK8 assay, immunofluorescence. Targets of biochanin A from pharmMappeer, SWISS, and TargetNet databases were screened. Targets of endothelial dysfunction were obtained from Genecards and OMIM databases. Protein-protein interaction, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomics analyses were used to analyze the potential target and the key pathway of the anti-endothelial dysfunction activity of biochanin A. To validate the potential target-drug interactions, molecular docking and molecular dynamics simulations were performed and the result was proved by western blot. RESULTS: It was found that biochanin A can promote the expression of ZO-1, reduce the expression of ICAM-1, which means improving endothelial dysfunction. A total of 585 targets of biochanin A from pharmMappeer, SWISS, and TargetNet databases were screened. A total of 10,832 targets of endothelial dysfunction were obtained from Genecards and OMIM databases. A total of 527 overlapping targets of endothelial dysfunction and biochanin A were obtained. AKT1, TNF-α, VCAM1, ICAM1, and NOS3 might be the key targets of the anti-endothelial dysfunction activity of biochanin A, and the key pathways might be PI3K-Akt and TNF signaling pathways. Molecular docking results indicated that the AKT1 and TNF-α had the highest affinity binding with biochanin A. CONCLUSION: This study indicates that biochanin A can target AKT1 and TNF-α to alleviate endothelial dysfunction induced by IL-6 in Perthes disease, which provides a theoretical basis for the treatment of Perthes disease by using biochanin A.

[Regulation of non-coding RNA in type H vessels angiogenesis of bone].

Objective: To summarize the regulatory effect of non-coding RNA (ncRNA) on type H vessels angiogenesis of bone. Methods: Recent domestic and foreign related literature about the regulation of ncRNA in type H vessels angiogenesis was widely reviewed and summarized. Results: Type H vessels is a special subtype of bone vessels with the ability to couple bone formation. At present, the research on ncRNA regulating type H vessels angiogenesis in bone diseases mainly focuses on microRNA, long ncRNA, and small interfering RNA, which can affect the expressions of hypoxia inducible factor 1α, platelet derived growth factor BB, slit guidance ligand 3, and other factors through their own unique ways of action, thus regulating type H vessels angiogenesis and participating in the occurrence and development of bone diseases. Conclusion: At present, the mechanism of ncRNA regulating bone type H vessels angiogenesis has been preliminarily explored. With the deepening of research, ncRNA is expected to be a new target for the diagnosis and treatment of vascular related bone diseases.

Association of survivin positive circulating tumor cell levels with immune escape and prognosis of osteosarcoma.

PURPOSE: Function of survivin protein (encoded by BIRC5) in circulating tumor cells (CTCs) of osteosarcoma (OS) has not been investigated. The goal of this study is to determine whether the expression of survivin protein of CTCs is associated with circulating immune cell infiltration and disease prognosis of OS. METHODS: Blood samples of 20 patients with OS were collected. CanPatrol™ CTC enrichment technology combined with in situ hybridization (ISH) was applied to enrich and test CTCs and survivin protein. Bioinformation analysis combined with data of routine blood test was used to verify the association between survivin and immune cell infiltration in circulatory system. To screen independent prognostic factors, Kaplan-Meier survival curve, univariate and multivariable Cox regression analyses were performed. RESULTS: Bioinformatics analysis showed that BIRC5 was strongly negatively related to lymphocyte, including T cell, NK cell and B cell, which released that BIRC5 played a key role in immune escape via reducing immune cell infiltration in circulatory system. Meanwhile, the number of survivin+ CTCs was significantly negatively connection with lymphocyte count (R = -0.56, p = 0.011), which was consistent with bioinformatics analysis. Kaplan-Meier curve showed that the overall survival rate in high survivin+ CTCs group was significantly lower than low group (88.9% vs 36.4%, p = 0.04). Multivariable Cox regression analyses showed that survivin+ CTCs were an independent prognostic factor (p = 0.019). CONCLUSION: These findings suggested that survivin protein played a key role in immune escape of CTCs and the presence of survivin+ CTCs might be a promising prognostic factor in OS patients.

A new insight of immunosuppressive microenvironment in osteosarcoma lung metastasis.

The lung is the primary organ for the metastasis of osteosarcoma. Although the application of neoadjuvant chemotherapy and surgery has remarkably improved the survival rate of patients with osteosarcoma, prognosis is still poor for those patients with metastasis. In this study, we performed further bioinformatics analysis on single-cell RNA sequencing (scRNA-seq) data published before, containing 75,317 cells from two osteosarcoma lung metastasis and five normal lung tissues. First, we classified 17 clusters, including macrophages, T cells, endothelial cells, and so on, indicating highly intratumoral heterogeneity in osteosarcoma lung metastasis. Next, we found macrophages in osteosarcoma lung metastasis did not have significant M1 or M2 polarizations. Then, we identified that T cells occupied the most abundant among all cell clusters, and found CD8+ T cells exhibited a low expression level of immune checkpoints in osteosarcoma lung metastasis. What is more, we identified C2_Malignant cells, and found CD63 might play vital roles in determining the infiltration of T cells and malignant cells in conventional-type osteosarcoma lung metastasis. Finally, we unveiled C1_Therapeutic cluster, a subcluster of malignant cells, was sensitive to oxfendazole and mevastatin, and the potential hydrogen-bond position and binding energy of oxfendazole-KIAA0907 and mevastatin-KIAA0907 were unveiled, respectively. Our results highlighted the power of scRNA-seq technique in identifying the complex tumor microenvironment of osteosarcoma lung metastasis, making it possible to devise precision therapeutic approaches.

Long non-coding RNA and circular RNA and coding RNA profiling of plasma exosomes of osteosarcoma by RNA seq.

Osteosarcoma (OS) is a primary bone tumor with high malignancy and the mechanism of hematogenous metastasis in OS is still not clear. The plasma exosomes derived from osteosarcoma play a key role in the process of tumor metastasis. Here, we established RNA-seq dataset for lncRNAs, circRNAs and mRNAs in plasma exosomes from 10 OS patients and 5 healthy donors. A total of 329.52 Gb of clean data was obtained. Besides, 1754 lincRNAs, 7096 known and 1935 new circRNA was identified. Finally, gene expression profiles and differentially expressed genes (DEGs) were analyzed among these 15 samples. There were 331 DEGs of mRNA, 132 of lincRNA and 489 of circRNA was obtained, respectively. This data set provides a significant resource for relevant researchers to excavate potential dysregulated lncRNAs, circRNAs and mRNAs of plasma exosomes in OS versus normal conditions.

Interaction gene set between osteoclasts and regulatory CD4+ T cells can accurately predict the prognosis of patients with osteosarcoma.

Osteoclasts (OCs) and regulatory CD4+ T cells (CD4+ Tregs) are important components in the tumor microenvironment (TME) of osteosarcoma. In this study, we collected six osteosarcoma samples from our previous study (GSE162454). We also integrated a public database (GSE152048), which included single cell sequencing data of 11 osteosarcoma patients. We obtained 138,192 cells and then successfully identified OCs and CD4+ Tregs. Based on the interaction gene set between OCs and CD4+ Tregs, patients from GSE21257 were distinguished into two clusters by consensus clustering analysis. Both the tumor immune microenvironment and survival prognosis between the two clusters were significantly different. Subsequently, five model genes were identified by protein-protein interaction network based on differentially upregulated genes of cluster 2. Quantitative RT-PCR was used to detect their expression in human osteoblast and osteosarcoma cells. A prognostic model was successfully established using these five genes. Kaplan-Meier survival analysis found that patients in the high-risk group had worse survival (p = 0.029). Therefore, our study first found that cell-cell communication between OCs and CD4+ Tregs significantly alters TME and is connected to poor prognosis of OS. The model we constructed can accurately predict prognosis for osteosarcoma patients.

A novel molecular signature for predicting prognosis and immunotherapy response in osteosarcoma based on tumor-infiltrating cell marker genes.

Background: Tumor infiltrating lymphocytes (TILs), the main component in the tumor microenvironment, play a critical role in the antitumor immune response. Few studies have developed a prognostic model based on TILs in osteosarcoma. Methods: ScRNA-seq data was obtained from our previous research and bulk RNA transcriptome data was from TARGET database. WGCNA was used to obtain the immune-related gene modules. Subsequently, we applied LASSO regression analysis and SVM algorithm to construct a prognostic model based on TILs marker genes. What's more, the prognostic model was verified by external datasets and experiment in vitro. Results: Eleven cell clusters and 2044 TILs marker genes were identified. WGCNA results showed that 545 TILs marker genes were the most strongly related with immune. Subsequently, a risk model including 5 genes was developed. We found that the survival rate was higher in the low-risk group and the risk model could be used as an independent prognostic factor. Meanwhile, high-risk patients had a lower abundance of immune cell infiltration and many immune checkpoint genes were highly expressed in the low-risk group. The prognostic model was also demonstrated to be a good predictive capacity in external datasets. The result of RT-qPCR indicated that these 5 genes have differential expression which accorded with the predicting outcomes. Conclusions: This study developed a new molecular signature based on TILs marker genes, which is very effective in predicting OS prognosis and immunotherapy response.

Isosinensetin alleviates estrogen deficiency-induced osteoporosis via suppressing ROS-mediated NF-κB/MAPK signaling pathways.

The formation of osteoclasts and their hyperactive bone resorption are related to the aggregation of intracellular reactive oxygen species (ROS). Flavonoids, derived from plant active ingredients, can alleviate the symptoms of osteoporosis (OP). Isosinensetin (Iss) is a flavonoid with antioxidant effects obtained mainly from citrus fruits, and its effect on osteoclastogenesis has not been reported. In this study, we investigated the antioxidant activity of Iss on osteoclast differentiation and function, as well as the therapeutic impact of Iss on OP. We found that Iss inhibited osteoclastogenesis and suppressed the bone resorption function of osteoclasts. Additionally, Iss reduced receptor activator of nuclear factor-κB ligand (RANKL)-induced intracellular ROS. Using quantitative real-time polymerase chain reaction and western blot, we further found that Iss inhibited osteoclast-specific genes and related proteins, while promoting the expression of antioxidant enzyme-related genes and proteins. Mechanistically, Iss reduces intracellular ROS by activating nuclear factor-erythroid 2-related factor 2 (Nrf2) and its related antioxidant enzymes and inhibits the downstream nuclear factor-κB (NF-κB) and mitogen-activated protein kinase (MAPK) signaling pathways of ROS, which in turn inhibits nuclear factor of activated T cells 1 (NFATc1), and ultimately inhibits osteoclastogenesis. In vivo, by micro-computed tomography (Micro-CT) assay and histological analyses, we found that Iss could reduce bone loss in ovariectomized (OVX) mice. Therefore, Iss has the potential as an OP preventative and therapeutic drug option.

Analysis of scRNA-seq and bulk RNA-seq demonstrates the effects of EVI2B or CD361 on CD8+ T cells in osteosarcoma.

Osteosarcoma (OS) is a common primary malignant tumor of the bone in children and adolescents. The five-year survival rate is estimated to be ~70% based on the currently available treatment modalities. It is well known that tumor-infiltrating immune cells (TIICs) that are the most important components in the tumor microenvironment can exert a killing effect on tumor cells. Therefore, in the present study, 85 RNA-sequencing OS samples were categorized into high- and low-immune score groups with ESTIAMATE. Based on the immune score groups, 474 differentially expressed genes (DEGs) were acquired using the LIMMA package of R language. Subsequently, 86 DEGs were taken through univariate COX regression analysis, of which 14 were screened out by least absolute shrinkage and selection operator regression analysis. Furthermore, multivariate COX regression analysis was performed to obtain 4 DEGs. Finally, ecotropic virus integration site 2B (EVI2B) or CD361 gene was screened out via Kaplan-Meier analysis. In addition, CIBERSORT algorithm was used to evaluate the proportion of 22 kinds of TIICs in OS. Correlation analysis revealed that the high expression level of EVI2B can elevate the infiltrated proportion of CD8+ T cells. Moreover, analysis of single cell RNA-sequencing transcriptome datasets and immunohistochemical staining uncovered that EVI2B was mainly expressed on CD8+ T cells and that EVI2B could promote the expression of granzyme A and K of CD8+ T cells to exhibit a potent killing effect on tumor cells. Therefore, EVI2B was identified as a protective immune-related gene and contributed to good prognosis in OS patients.

Corrigendum: Single-cell transcriptomics reveals the complexity of the tumor microenvironment of treatment-naive osteosarcoma.

[This corrects the article DOI: 10.3389/fonc.2021.709210.].

Construction of a transcription factor-miRNA-mRNA interactive network elucidates underlying pathogenesis for osteosarcoma and validation by qRT-PCR.

PURPOSE: Osteosarcoma is characterized by features of rapid growth and early metastasis with a poor prognosis. The aim of our research is to investigate the potential transcription factor (TF)-miRNA-mRNA regulatory mechanism in osteosarcoma utilizing bioinformatics methods and validate by qRT-PCR. METHODS: The microRNA (miRNA) expression profiling datasets (GSE28423 and GSE65071) and mRNA expression profiling dataset GSE33382 were collected from the Gene Expression Omnibus (GEO) database. Differentially expressed miRNAs (DEMs) and differentially expressed genes (DEGs) were screened using the limma package. Then, the TransmiR v2.0, miRDB, and Targetscan 7.2 database were applied for the acquisition of TF-miRNA and miRNA-mRNA interaction relationships, respectively. Finally, we built a TF-miRNA-mRNA interactive network. Furthermore, survival analysis was performed to identify sub-network with prognostic value and validate through qRT-PCR. RESULTS: Eight overlapping DEMs and 682 DEGs were identified. Based on bioinformatics methods, 30 TF-miRNA interaction pairs and 25 miRNA-mRNA interaction pairs were screened. Finally, we constructed a TF-miRNA-mRNA regulatory network. Furthermore, laminin subunit gamma 1 (LAMC1) and thrombospondin-1 (THBS1), which involved in the network, were determined to have prognostic value and the corresponding subnetwork was identified. qRT-PCR results showed that LAMC1 mRNA expression was higher in osteosarcoma cells. CONCLUSION: Based on the survival analysis, a TF-miRNA-mRNA sub-network, that is TFs (SPI1, HEY1, and CEBPB)-hsa-miR-338-3p-target genes (LAMC1 and THBS1) was established. In conclusion, the construction of a potential TF-related regulatory network will help elucidate the underlying pathological mechanisms of osteosarcoma, and may provide novel insights for the diagnosis and treatment of osteosarcoma.

Isoliensinine Suppresses Osteoclast Formation Through NF-κB Signaling Pathways and Relieves Ovariectomy-Induced Bone Loss.

Osteoporosis is among the major contributors of pathologic fracture in postmenopausal women, which is caused by the bone metabolic disorder owing to the over-activation of osteoclasts. Inhibition of osteoclast differentiation and maturation has become a mainstream research interest in the prevention of osteoporosis. Isoliensinine (Iso) is a dibenzyl isoquinoline alkaloid with antioxidant, anti-inflammatory, and anti-cancer activities. However, whether it can be used as a potential treatment for osteoporosis remains undiscovered. Here, we investigated whether Iso might suppress the differentiation of osteoclasts in vitro and in vivo to play an anti-osteoporosis role. Our results showed that Iso inhibits the formation of mature multinuclear osteoclasts induced by RANKL, the bone resorption, and the osteoclast-specific genes expression by blocking the nuclear translocation of NF-κB p65, and the effect was in a dosage-dependent way. Furthermore, we investigated the therapeutic effect of Iso on osteoporosis in ovariectomized (OVX) mice. We found that Iso attenuated bone loss in the OVX mice and significantly promoted BS, Conn. DN, Tb.Th, TB.N, and BV/TV Index. All in all, Iso showed a prominent effect of osteoclast inhibition, with great promise for treating osteoporosis.

Diagnosis and treatment of primary hyperparathyroidism with pathological fracture of the limbs: A retrospective observational study.

Primary hyperparathyroidism (PHPT) with pathological fracture is rare, and the early symptoms of PHPT lack specificity, leading to misdiagnosis. Therefore, this study aimed to summarize the clinical characteristics and treatment of PHPT patients with pathological fractures and to improve the attention of orthopedic clinicians to PHPT. It is a retrospective study, 2226 patients with hyperparathyroidism in our hospital from 2009 to 2019 were screened, excluding secondary hyperparathyroidism and patients without limb fracture, and the remaining 20 patients with PHPT accompanied by pathological fractures were finally analyzed. Parathyroid hormone (PTH) and calcium levels were compared on the first postoperative day, and the prognosis of the patients was assessed by bone mineral density and Visual Analogue Scale scores at 3 and 12 months postoperatively. The early symptoms of PHPT patients in this study included urinary calculi (80%), bone pain (30%), and digestive tract symptoms (25%). Fourteen (70%) cases were misdiagnosed at the initial diagnosis. After parathyroidectomy, the blood calcium and PTH levels decreased significantly in all patients (P < .05). For the treatment of fracture, 9 of the patients underwent surgical treatment of the fracture, while the remaining patients received splint external fixation. The follow-up time was 4.60 ± 0.62 years (1-10 years). All patients recovered well from the fracture, the symptoms of systemic bone pain were markedly improved, and bone mineral density was significantly improved after surgery. Orthopedic surgeons need to avoid misdiagnosis and pay attention to the early symptoms in PHPT patients with pathological fracture, and better therapeutic effects can be obtained by combining parathyroidectomy with fractures fixation.

Influence of ulnar bow sign on surgical treatment of missed Bado type I Monteggia fracture in children.

The present study aimed to explore the influence of ulnar bow on the surgical treatment of Bado type I missed Monteggia fracture in children. A retrospective review of 28 patients was conducted between November 2010 and June 2020. All patients were treated with open reduction of the radial head and ulnar opening wedge osteotomy without annular ligament reconstruction. Four months (range 1-12 months) was the mean interval between injury onset and surgery. The average age of patients at the time of surgery was 6.1 years old (range 2-10 years old). The maximum ulnar bow (MUB) and MUB position (P-MUB) via radiography were evaluated. The patients were divided into two groups according to P-MUB, as follows: middle group (A) included 17 cases, and the MUB was located at 40-60% of the distal ulna; and distal group (B) included 11 cases, and the MUB was located at 20-40% from the distal end of the ulna. The mean follow-up period was 33 months (range 6-102 months). At the last follow-up, all the children showed stable reduction of the radial head, and the flexion function of elbow joint improved after the operation (P < 0.05). Group A presented a larger ratio of maximum ulnar bow (R-MUB) and angle of ulnar osteotomy (OA) than group B (P < 0.05). The osteotomy angle was positively correlated with the R-MUB (R2 = 0.394, P = 0.038). The osteotomy angle was positively correlated with the P-MUB (R2 = 0.683, P = 0.000). The R-MUB was proportional to the P-MUB (R2 = 0.459, P < 0.0001). The regression equation of P-MUB and osteotomy angle was as follows: OA = 32.64* P-MUB + 7.206. If the ulnar bow was positioned at the middle ulna, then a stable reduction of radial head needed to be achieved through a large angle in the ulnar osteotomy. If the position of maximum ulnar bow (P-MUB) was closer to the middle of the ulna, or the ratio of maximum ulnar bow (R-MUB) was larger, then the osteotomy angle was larger.

The regulatory mechanism of LncRNA-mediated ceRNA network in osteosarcoma.

Aberrantly expressed lncRNAs have been reported to be closely related to the oncogenesis and development of osteosarcoma. However, the role of a dysregulated lncRNA-miRNA-mRNA network in osteosarcoma in the same individual needs to be further investigated. Whole transcriptome sequencing was performed on the tumour tissues and matched paratumour tissues of three patients with confirmed osteosarcoma. Two divergent lncRNA-miRNA-mRNA regulatory networks were constructed in accordance with their biological significance. The GO and KEGG analysis results of the mRNAs in the two networks revealed that the aberrantly expressed lncRNAs were involved in regulating bone growth and development, epithelial cell proliferation, cell cycle arrest and the N-terminal acetylation of proteins. The survival analysis results of the two networks showed that patients with high expression of GALNT3, FAM91A1, STC2 and SLC7A1 end in poorer prognosis. Likewise, patients with low expression of IGF2, BLCAP, ZBTB47, THRB, PKIA and MITF also had poor prognosis. A subnetwork was then constructed to demonstrate the key genes regulated by aberrantly expressed lncRNAs at the posttranscriptional level via the ceRNA network. Aberrantly expressed lncRNAs in osteosarcoma tissues regulate genes involved in cellular proliferation, differentiation, angiogenesis and the cell cycle via the ceRNA network.

ATG16L1 is a Potential Prognostic Biomarker and Immune Signature for Osteosarcoma: A Study Based on Bulk RNA and Single-Cell RNA-Sequencing.

BACKGROUND: Osteosarcoma is a common solid malignancy of the bone in children and adolescents, and its metastasis and recurrence are the principal causes of poor treatment outcomes. METHODS: Autophagy-related genes were used to cluster osteosarcoma patients by consensus clustering analysis using the GSE21257 database. Differentially expressed genes (DEGs) were identified by limma package. Multiple-gene risk signature was constructed using least absolute shrinkage and selection operator (LASSO) analysis and Cox regression analyses. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to determine gene expression levels. Then, single-cell RNA-sequencing dataset GSE152048 were used to identify the correlation between the DEGs and effector molecules expressed in specific tumor-infiltrating immune cells. RESULTS: Two clusters were identified in the consensus clustering analysis, which were confirmed by principal component analysis. Limma analysis revealed that 15 genes were related, and 9 genes were screened using protein-protein interaction network and LASSO regression analysis. Cox regression analyses identified 5 genes. Combined with survival analysis, only the autophagy related 16 like 1 gene (ATG16L1) was significant. The results of qRT-PCR showed low expression levels of ATG16L1 in tumor cells group. Immune infiltration analysis revealed significantly lower expression of CD8+ T cells in the high ATG16L1 gene expression group. ScRNA-seq revealed that in the ATG16L1+ CD8 + T cell group, the expression of GZMB was lower, whereas the expression of ITGA1 was higher. These results showed that ATG16L1 is an immune-related gene, which is associated with poor prognosis in patients with osteosarcoma. CONCLUSION: ATG16L1 is a potential prognostic biomarker and immune signature and may be a therapeutic target for osteosarcoma.

Significance of the neutrophil-to-lymphocyte ratio in predicting the response to neoadjuvant chemotherapy in extremity osteosarcoma: a multicentre retrospective study.

BACKGROUND: At present, no predictive factor has been validated for the early efficacy of neoadjuvant chemotherapy (NACT) in osteosarcoma. The purpose of this study was to investigate the significance of the neutrophil-to-lymphocyte ratio (NLR) in predicting the response to NACT in extremity osteosarcoma. METHODS: Pathological complete response (pCR) was used to assess the efficacy of NACT. Receiver operating characteristic (ROC) curves and the Youden index (sensitivity + specificity-1) were used to determine the optimal cut-off values of the NLR. Univariate and multivariate analyses using logistic regression models were conducted to confirm the independent factors affecting the efficacy of NACT. RESULTS: The optimal NLR cut-off value was 2.36 (sensitivity, 80.0%; specificity, 71.3%). Univariate analysis revealed that patients with a smaller tumour volume, lower stage, lower NLR and lower PLR were more likely to achieve pCR. Multivariate analyses confirmed that the NLR before treatment was an independent risk factor for pCR. Compared to patients with a high NLR, those with a low NLR showed a more than 2-fold higher likelihood of achieving pCR (OR 2.82, 95% CI 1.36-5.17, p = 0.02). CONCLUSION: The NLR is a novel and effective predictive factor for the response to NACT in extremity osteosarcoma patients. Patients with a higher NLR showed a lower percentage of pCR after NACT.