A Novel Predictive Model Associated with Osteosarcoma Metastasis

Figures & data

Figure 1 Differential expression profiles. (A) Differential lncRNA between metastatic and non-metastatic OS; (B) Differential mRNA between metastatic and non-metastatic OS.

Figure 2 Identification of the prognosis-related lncRNA. (A and B) LASSO coefficient profiles; (C) multivariate cox analysis of seven model lncRNAs; (D) the risk scores, ROC curve, and Kaplan–Meier curve of patients.

Figure 3 Kaplan–Meier curve of seven prognosis-related lncRNAs. (A) AL512422.1. (B) AL008718.3. (C) AC006033.2. (D) AL357507.1. (E) AL360182.2. (F) C5orf66-AS1. (G) CEBPA-DT.

Table 1 Transcription Factors Associated with 24 Intersecting mRNAs

Category	Term	Count	%	P value	Genes	List Total	Pop Hits	Pop Total	Fold Enrichment	Bonferroni	Benjamini	FDR
UCSC_TFBS	AHRARNT	17	70.83333	0.006029	GADD45GIP1, CGREF1, MESP1, HMGA1, IKZF1, HOXC13, ITGAL, POU2F2, HOXC12, SORL1, VNN1, HRH2, SFRP5, FAM198B, HRC, NRCAM, SLC27A5	23	9110	20,726	1.681583	0.652961	0.601557	0.601557
UCSC_TFBS	AP4	19	79.16667	0.006875	GADD45GIP1, CGREF1, MESP1, GPR65, IKZF1, HOXC13, ITGAL, POU2F2, HOXC12, SORL1, VNN1, HRH2, MYC, SFRP5, FAM198B, HRC, NRCAM, TLR7, SLC27A5	23	11,244	20,726	1.522721	0.700987	0.601557	0.601557
UCSC_TFBS	ARP1	15	62.5	0.018577	GADD45GIP1, CGREF1, MESP1, HMGA1, IKZF1, HOXC13, ITGAL, POU2F2, HOXC12, SORL1, VNN1, HRH2, HRC, TLR7, SLC27A5	23	8167	20,726	1.65507	0.96243	0.87694	0.87694
UCSC_TFBS	NRF2	11	45.83333	0.024629	GADD45GIP1, HRH2, FAM198B, HRC, HMGA1, NRCAM, IKZF1, ITGAL, HOXC12, SORL1, SLC27A5	23	5045	20,726	1.964804	0.987273	0.87694	0.87694
UCSC_TFBS	IK3	14	58.33333	0.025055	GADD45GIP1, MESP1, MMP3, IKZF1, HOXC13, ITGAL, POU2F2, SORL1, MYC, SFRP5, PDPN, FAM198B, NRCAM, TLR7	23	7539	20,726	1.673408	0.988211	0.87694	0.87694
UCSC_TFBS	NF1	12	50	0.043898	GADD45GIP1, CGREF1, HRH2, GPR65, HMGA1, NRCAM, HOXC13, ITGAL, POU2F2, HOXC12, SORL1, SLC27A5	23	6332	20,726	1.707765	0.999612	0.935221	0.935221
UCSC_TFBS	TCF11MAFG	15	62.5	0.047363	CGREF1, MESP1, MMP3, HMGA1, IKZF1, ITGAL, POU2F2, SORL1, VNN1, HRH2, SFRP5, FAM198B, HRC, NRCAM, SLC27A5	23	9047	20,726	1.494082	0.999795	0.935221	0.935221
UCSC_TFBS	HEN1	16	66.66667	0.055655	GADD45GIP1, MESP1, HMGA1, IKZF1, HOXC13, ITGAL, POU2F2, HOXC12, SORL1, VNN1, MYC, SFRP5, HRC, NRCAM, TLR7, SLC27A5	23	10,160	20,726	1.419103	0.999956	0.935221	0.935221
UCSC_TFBS	STAT	11	45.83333	0.069174	CGREF1, MYC, SFRP5, PDPN, MMP3, NRCAM, TLR7, IKZF1, ITGAL, HOXC12, SORL1	23	5930	20,726	1.671574	0.999996	0.935221	0.935221
UCSC_TFBS	CEBPB	15	62.5	0.080449	GADD45GIP1, CGREF1, MESP1, IKZF1, HOXC13, ITGAL, POU2F2, HOXC12, SORL1, VNN1, MYC, FAM198B, HRC, NRCAM, TLR7	23	9630	20,726	1.40363	1	0.935221	0.935221
UCSC_TFBS	E47	16	66.66667	0.081955	GADD45GIP1, CGREF1, MESP1, HMGA1, IKZF1, HOXC13, POU2F2, HOXC12, SORL1, MYC, SFRP5, FAM198B, HRC, NRCAM, TLR7, SLC27A5	23	10,601	20,726	1.360069	1	0.935221	0.935221
UCSC_TFBS	HSF2	12	50	0.094283	MESP1, CGREF1, VNN1, MYC, PDPN, HRC, NRCAM, IKZF1, ITGAL, HOXC12, SORL1, SLC27A5	23	7109	20,726	1.521109	1	0.935221	0.935221
UCSC_TFBS	RP58	14	58.33333	0.098286	GADD45GIP1, CGREF1, HMGA1, IKZF1, HOXC13, ITGAL, POU2F2, HOXC12, SORL1, PDPN, FAM198B, NRCAM, TLR7, SLC27A5	23	8944	20,726	1.410535	1	0.935221	0.935221

Figure 4 Co-expression network and enrichment analysis. (A) Co-expression network of lncRNA and mRNA; (B) Co-expression network of mRNA and TF; (C) GO analysis of the genes in co-expression network.

Figure 5 GSVA enrichment analysis of seven prognosis-related lncRNAs. (A) AL512422.1. (B) AL008718.3. (C) AC006033.2. (D) AL357507.1. (E) AL360182.2. (F) C5orf66-AS1. (G) CEBPA-DT.

Figure 6 LncRNA and immune analysis. (A) Relationship between high/low group and immune cells; (B) Relationship between seven prognosis-related lncRNAs and immune cells. *<0.05, **<0.01.

Figure 7 The qPCR result of seven lncRNAs in five OS highly metastasis and weakly metastasis cell lines. (A) AL512422.1. (B) AL008718.3. (C) AC006033.2. (D) AL357507.1. (E) AL360182.2. (F) C5orf66-AS1. (G) CEBPA-DT. *<0.05, **<0.01, ***<0.001, ****<0.0001.