MicroRNAs as a novel class of diagnostic biomarkers for the detection of osteosarcoma: a meta-analysis

Abstract

MicroRNAs (miRNAs) have been considered as promising diagnostic biomarkers for many diseases, especially for cancers. Numerous studies have reported the value of miRNAs in the diagnosis of osteosarcoma (OS), but the results vary greatly across different studies. Therefore, we conducted this meta-analysis to assess the prospective diagnostic value of miRNAs in diagnosing OS. All relevant articles from prior to July 28, 2017 were selected from PubMed, EMBASE, Web of Science, Cochrane Library, Chinese National Knowledge Infrastructure, and Wan-fang databases. The Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2) was performed to assess the quality of each article. A random-effects model was used to pool the sensitivity and specificity of the positive likelihood ratio (PLR), negative likelihood ratio (NLR) and, diagnostic odds ratio (DOR) together with the area under the curve (AUC) to evaluate diagnostic values. Seventeen studies comprising 2,214 OS patients and 1,534 healthy humans were included in our meta-analysis. The pooled estimations indicated that the miRNAs had a high accuracy for diagnosing OS, with a sensitivity of 0.82, specificity of 0.88, PLR of 10.96, NLR of 0.20, DOR of 54.55, and AUC of 0.93. Twenty-five miRNAs were differentially expressed in OS, including 17 upregulated and 8 downregulated. These miRNAs were correlated with survival time, tumor size, cell differentiation, tumor node metastasis staging, metastasis, tumor/cell invasion, pathological type, and response to radiotherapy and chemotherapy. Several different miRNAs are expressed in OS, and some of them might be potential biomarkers for the early diagnosis of OS.

Keywords:

• osteosarcoma
• miRNAs
• diagnosis
• biomarker
• meta-analysis

Introduction

Osteosarcomas (OS) originate from the mesenchyme tissue and are the most common malicious bone tumor. They occur primarily in the metaphysis of long bones and always affect adolescents and children.^1–3 The incidence of OS iŝ5 per million among Asian populations.⁴ Although some effective treatments, such as advanced surgical methods, adjuvant chemotherapy and radiotherapy, have significantly increased the survival rate of patients with OS, its 5-year survival rate is not high.^1,3,5 The main reasons are highly malignancy, strong invasiveness, and early pulmonary metastasis.^6,7 Early diagnosis and treatment of cancer are still two big challenges for clinicians and, therefore, it has become necessary to find an effective method for the initial diagnosis of OS.

It is well established that the standard OS therapy consists of early operation combined with neoadjuvant and adjuvant chemotherapy, but once metastasis occurs, this type of treatment does not work well. Currently, the diagnosis of osteosarcoma depends primarily on X-ray, localized symptoms, and pathology results, and it is difficult to detect early metastasis. Therefore, it is urgent to find innovative approaches for detecting OS in patients at the initial stage and for developing more efficient therapeutic methods to achieve a better prognosis.

miRNAs are a class of small non-coding RNA molecules (containing ~22 nucleotides) that play a critical role in gene silencing and post-transcriptional regulation of gene expression and occur broadly in plants, animals, and some viruses.^8,9 Researchers have found that ~30% of genes in the human genome are regulated by miRNAs¹⁰ and that miRNAs participate in tumor development, metastasis, angiogenesis, and immune reactions.^4,11,12 Thus, miRNAs play very important roles in homeostasis and tumorigenesis.^13,14 Recently, numerous reports have shown that miRNAs are promising diagnostic and prognostic tools for cancers.^15–17 The miRNAs have been confirmed to be promising diagnostic and prognostic markers in nasopharyngeal carcinoma¹⁸ and glioma¹⁹ and in chronic lymphocytic leukemia,²⁰ acute lymphoblastic leukemia,²¹ breast cancer,^22–24 gastric cancer,^25,26 prostate cancer,^27,28 pancreatic cancer,²⁹ and colorectal cancer.³⁰

The miRNAs in OS have also been widely tested. Cellular miRNAs can be identified rapidly and accurately from tissue and blood, especially serum and plasma,³¹ but the results are controversial. For example, Cai et al showed that miR-195 may be an indicative and predictive marker for OS,³² but other researchers stressed that miR-222 can be used to predict OS.³³ Until now, no confirmed agreement has been reached. In this study, we identified and reviewed the most recent literature on miRNA and performed a meta-analysis to evaluate the diagnostic and prognostic value of miRNAs for OS.

Patients and methods

Search strategy and study selection

We conducted a literature search in the PubMed, EMBASE, Web of Science, Cochrane Library, Chinese National Knowledge Infrastructure, and Wan-fang databases to identify original articles analyzing the diagnostic value of miRNAs in OS. The last retrieval was performed before July 28, 2017. The medical subject headings (Mesh) and the entry terms employed in the literature search included 1) “osteosarcoma” or “osteosarcoma tumor”, 2) “microRNAs” or “miRNAs” or “primary microRNA” or “circulating microRNAs” or “circulating miRNA”, and 3) “diagnosis” or “sensitivity and specificity” or “receiver operating characteristics” or “ROC curve” or “ROC analyses” or “predictive value”. There were no restrictions on the language or publication date. The reference lists of review articles and selected papers were also skimmed to identify any additional eligible studies. Retrieved articles were independently screened and carefully evaluated by two investigators (LH and LZH). Two initial investigators (LP and CL) and another investigator (JC) participated in the data extraction. A fourth investigator (YAX) resolved any disagreements. We extracted data from the selected studies and input them into tables based on dissimilar characteristics of interest.

Selection criteria

This meta-analysis was carried out according to the guidelines of the preferred reporting items for systematic reviews and meta-analyses (PRISMA) statement. Studies that met the following criteria were included: 1) studies that evaluated the diagnostic value of miRNAs in OS patients; 2) studies from which meaningful data, such as sensitivity, specificity, and AUC values, were obtained; 3) studies that included specimens such as serum, plasma, or peripheral blood, and studies that used quantitative real-time polymerase chain reaction; 4) every individual study that involved both patients with OS and healthy control groups; and 5) articles that were written in English or Chinese. Accordingly, studies that met any of the following criteria were excluded: 1) duplicate of data already included in the analysis; 2) non-original paper, such as conference abstracts, letters, and reviews; 3) studies unrelated to OS; 4) animal experiments; or 5) studies that lacked complete data. For studies with overlapping populations or data, only the most complete report was reflected in the present meta-analysis. The selection process and inclusion and exclusion outcomes are shown in Figure 1.

Figure 1 Flow diagram of the selected eligible studies according to the inclusion and exclusion criteria.

Data extraction

Two investigators (LP and JC) examined the full texts of eligible publications and extracted relevant data using a standardized table. From every included study, the following information was extracted: 1) study features, including the first author’s surname, publication year, and country; 2) participants’ general features, including ethnicity, type of case and control groups, and number of specimens; 3) data needed for the diagnostic meta-analysis, including studied miRNAs, specimen, sensibility, and specificity; and 4) information needed for quality assessment.

Quality assessment

In this meta-analysis, we comprehensively evaluated the quality of each included study based on the revised Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2) checklist. The QUADAS-2 tool is composed of 4 key domains: patient selection, index test, reference standard, and flow and timing. These are supported by signaling questions to aid in judging.³⁴ Table 1 presents a summary of all included studies.

Table 1 Characteristics of the 17 reports included in our meta-analysis

Study ID	Publication year	Country	Ethnicity	OS patients	Healthy controls	SEN (%)	SPE (%)	miRNAs profiled	Specimen	QUADAS-2 score
Wang et al^Citation48	2017	China	Asian	72	40	72.00	86.00	miR-491-5p	Serum	7
Fujiwara et al^Citation46	2017	Japan	Asian	14	22	71.40	92.30	miR-25-3p	Serum	6
						64.30	84.60	miR-17-5p	Serum	6
Xiao et al^Citation47	2016	China	Asian	20	20	90.00	95.00	miR-125b	Serum	6
Liu et al^Citation49	2016	China	Asian	114	114	84.20	88.60	miR-300	Serum	6
Cao et al^Citation35	2016	China	Asian	60	20	83.70	94.50	miR-326	Serum	6
Zhang et al^Citation33	2016	China	Asian	57	57	66.70	84.20	miR-222	Serum	7
Dong et al^Citation36	2016	China	Asian	112	50	89.50	97.20	miR-223	Serum	7
Tang et al^Citation42	2015	China	Asian	166	60	70.01	98.30	miR-27a	Serum	5
Yang et al^Citation44	2015	China	Asian	108	50	65.70	100.00	miR-221	Serum	7
Zhou et al^Citation45	2015	China	Asian	88	60	88.33	76.67	miR-199a-5p	Serum	7
Wang et al^Citation39	2015	China	Asian	80	20	96.20	92.50	miR-152	Serum	7
Wang et al^Citation43	2015	China	Asian	100	20	74.00	100.00	miR-191	Serum	4
Lian et al^Citation37	2015	China	Asian	90	90	91.10	94.40	miR-195-5p	Plasma	7
								miR-199a-3p
								miR-320a
								miR-374a-5p
Cai et al^Citation32	2015	China	Asian	166	60	88.00	83.30	miR-195	Serum	4
Ma et al^Citation38	2014	China	Asian	89	89	69.70	83.10	miR-148a	Plasma	6
Ren et al^Citation41	2014	China	Asian	22	30	96.00	97.00	miR-199b-5p/miR-124	Plasma	5
						27.00	90.00	miR-181b/miR-124	Plasma	5
						59.00	97.00	miR-451/miR-124	Plasma	5
Ouyang et al^Citation40	2013	China	Asian	40	40	90.50	93.80	miR-21	Plasma	6
								miR-199a-3p
								miR-143

Abbreviations: miRNA, microRNA; OS, osteosarcoma; SEN, sensitivity; SPE, specificity; QUADAS-2, Quality Assessment of Diagnostic Accuracy Studies-2.

Statistical analysis

The analytical software Meta-Disc version 1.4 and the statistical software Review Manager (RevMan 5.3) were utilized to analyze the diagnostic value and to calculate the 95% CI. The data were considered statistically significant when two-sided p<0.05. To indicate substantial heterogeneity inconsistency, I² was computed. An I² value >50% indicated substantial heterogeneity. Based on the heterogeneity values, we used the random-effect model throughout the statistical tests to compute the data. To investigate the possible sources of heterogeneity, meta-regression analyses were carried out. A subgroup analysis of the specimen type was performed for sensitivity, specificity, positive likelihood ratio (PLR), negative likelihood ratio (NLR) and, diagnostic odds ratio (DOR). The summary receiver operating characteristic (SROC) curve and the area under the curve (AUC) were also used to predict sensitivity and specificity. Additional data from validation studies that were not appropriate for pooling are listed in Table 2.

Table 2 Summary estimates of diagnostic criteria and their 95% CIs

Subgroup	Sensitivity (95% CI)	Specificity (95% CI)	PLR (95% CI)	NLR (95% CI)	DOR (95% CI)	AUC (95% CI)
miRNA profiling single	0.79 (0.77–0.82)	0.82 (0.78–0.85)	9.27 (3.02–28.46)	0.25 (0.19–0.32)	36.13 (14.71–88.71)	0.90 (0.88–0.92)
Sample types
Serum base	0.80 (0.77–0.82)	0.83 (0.79–0.88)	9.46 (2.96–30.20)	0.25 (0.19–0.33)	36.78 (14.79–91.5)	0.89 (0.87–0.91)
Plasma base	0.85 (0.82–0.87)	0.93 (0.91–0.95)	12.32 (7.28–20.83)	0.15 (0.07–0.32)	86.05 (32.76–225.99)	0.98 (0.96–0.99)
Overall studies	0.82 (0.80–0.85)	0.88 (0.86–0.90)	10.96 (4.67–25.76)	0.20 (0.15–0.28)	54.55 (27.38–108.67)	0.93 (0.92–0.94)

Abbreviations: miRNA, microRNA; PLR, positive likelihood ratio; NLR, negative likelihood ratio; DOR, diagnostic odds ratio; AUC, area under the curve; CI, confidence interval.

Results

Basic characteristics of included studies

As shown in Figure 1, a total of 1,167 published records related to the diagnostic value of miRNAs in OS were retrieved by a primary literature search (n=1,167) and hand search (n=0). Duplicate records (n=288) were removed. Then, 831 of the remaining 879 articles that were about animal experiments, that included other diseases, that did not include blood samples, and that were case reports or reviews were excluded, and full texts of the remaining 48 articles were evaluated for eligibility. The data were not complete in 31 of the 48 remaining articles. Therefore, this meta-analysis was conducted based on the remaining 17 articles.^{32,33,35–49} The main characteristics of the included articles are presented in Table 1.

The publication years of the 17 included articles were 2017 (n=2), 2016 (n=5), 2015 (n=6), 2014 (n=3), and 2012 (n=1). The dominant ethnicity of the study subjects was Asian. The expression levels of miRNAs of the plasma or serum samples were analyzed by quantitative real-time polymerase chain reaction.

Diagnostic accuracy of miRNAs in distinguishing OS from healthy controls

Table 2 presents the accuracy of the miRNAs in differentiating OS from healthy controls. A total of 17 studies covering 25 types of miRNAs and involving 3,748 participants (2,214 OS and 1,534 healthy controls) were included in the pooled analysis. The pooled sensitivity (Figure 2A) and specificity (Figure 2B) of the miRNAs for the diagnosis of OS were 0.82 (95% CI: 0.80–0.85) and 0.88 (95% CI: 0.86–0.90), respectively. The PLR, NLR, DOR, and AUC for the total miRNA levels in the 17 studies (Figures 3 and 4), with the corresponding 95% CIs, were 10.96 (95% CI: 4.67–25.76), 0.20 (95% CI: 0.15–0.28), 54.55 (95% CI: 27.38–108.67), and 0.93 (95% CI: 0.92–0.94), respectively.

Figure 2 Forest plot of diagnostic value indices of miRNA expression in OS patients. (A) Sensitivity and (B) specificity of miRNA expression in diagnosis of OS in all studies.

Abbreviations: CI, confidence interval; miRNA, microRNA; OS, osteosarcoma.

Figure 3 Forest plot of diagnostic value indices of miRNA expression in OS patients. (A) Positive and (B) negative LR of miRNA expression in diagnosis of OS in all studies.

Abbreviations: CI, confidence interval; miRNA, microRNA; OS, osteosarcoma; LR, likelihood ratio.

Figure 4 Forest plot of diagnostic value indices of miRNA expression in OS patients. (A) Diagnostic odds ratio and (B) SROC of miRNA expression in diagnosis of OS in all studies.

Abbreviations: CI, confidence interval; miRNA, microRNA; OR, odds ratio; OS, osteosarcoma; SROC, summary receiver operating characteristic curve; AUC, area under the curve.

In this meta-analysis, we conducted subgroup analyses of blood and plasma samples. The analyses results showed between-study heterogeneity. For a single miRNA assay, the pooled sensitivity, specificity, PLR, NLR, DOR, and AUC were 0.79 (95% CI: 0.77–0.82), 0.82 (95% CI: 0.78–0.85), 9.27 (95% CI: 3.02–28.46), 0.25 (95% CI: 0.19–0.32), 36.13 (95% CI: 14.71–88.71), and 0.90 (95% CI: 0.88–0.92), respectively. To detect whether plasma-based assays had a better diagnostic accuracy than serum-based assays, we determined the number of studies that used plasma/serum sample: of the 25 studies that assessed miRNAs, 11 used plasma samples and 14 used serum samples. The sensitivity of the studies that used plasma was 0.85 (95% CI: 0.82–0.87) and the specificity was 0.93 (95% CI: 0.91–0.95), with a pooled DOR of 86.05 (95% CI: 32.76–225.99) and AUC of 0.98 (95% CI: 0.96–0.99). Nevertheless, for serum-based studies, the sensitivity was 0.80 (95% CI: 0.77–0.82) and the specificity was 0.83 (95% CI: 0.79–0.88), with a pooled DOR of 36.78 (95% CI: 14.79–91.5) and AUC of 0.89 (95% CI: 0.87–0.91).

Clinical studies that supplied the ratio of miRNAs in OS patients

The results showed that different miRNAs were differentially expressed in patients with OS and healthy controls, with 17 upregulated miRNAs (including miR-222, miR-27a, and miR-221), 8 downregulated miRNAs (including miR-326, miR-223, and miR-15), and the same miRNA expressed in different tissues in distinctly different patterns. Our research shows that the expression of miR-195 and miR-199a-3p is much higher in patients with OS and that their expression in plasma and serum exhibits distinctive difference between OS and control groups. The detailed results are shown in Table 3.

Table 3 Comparison of miRNAs in blood samples from patients with OS and healthy humans

Upregulated miRNAs (n=17)	Downregulated miRNAs (n=8)
miR-222; miR-27a; miR-221; miR-199a-5p;	miR-326; miR-223;
miR-191; miR-195-5p; miR-199a-3p; miR-320a;	miR-152; miR-195;
miR-374a-5p; miR-148a; miR-199b-5p/miR-124;	miR-199a-3p; miR-143;
miR-21; miR-25-3p; miR-17-5p; miR-300;	miR-491-5p; miR-125b
miR-451/miR-124; miR-181b/miR-124

Abbreviation: miRNA, microRNA; OS, osteosarcoma.

Association between clinicopathological factors and differentially expressed miRNAs

The association between differentially expressed miRNAs and clinicopathological factors of OS patients is presented in Table 4. The individually expressed miRNAs were associated with survival time, tumor size, cell differentiation, tumor node metastasis staging, metastasis, tumor/cell invasion, pathological type, and response to radiotherapy and chemotherapy.

Table 4 Association between the expression of miRNAs and clinicopathological features of OS patients

Clinicopathological features	Expression of miRNAs
Survival time	miR-195, miR-27a, miR-221, miR-326, miR-148a, miR-152, miR-222, miR-223, miR-191, miR-491-5p, miR-300
Tumor size	miR-148a, miR-222, miR-191
Cell differentiation	miR-326
TNM staging	miR-195, miR-27a, miR-221, miR-326, miR-222, miR-223, miR-191, miR-491-5p, miR-300
Metastasis	miR-195, miR-27a, miR-221, miR-326, miR-21, miR-143, miR-148a, miR-222, miR-223, miR-191, miR-195-5p, miR-199a-3p, miR-491-5p, miR-300
Tumor/cell invasion	miR-326, miR-223
Pathological type	miR-21, miR-143, miR-199a-3p, miR-199a-3p, miR-320a
Radiotherapy/chemotherapy	High miR-27a expression poorly responded to chemotherapy

Abbreviations: miRNA, microRNA; OS, osteosarcoma; TNM, tumor node metastasis.

Discussion

Over the past several decades, we have achieved a lot in terms of the diagnostic methods for OS and significantly improved the treatment and prognosis of patients suffering from OS. However, the diagnosis of early stage OS and early lung metastasis is still challenging for clinicians. Finding other more sensitive and noninvasive methods for diagnosis is particularly important for the future treatment and diagnosis of OS.

The functional contributions of miRNAs in the development and progression of malignancies have resulted in the development of new diagnostic approaches. Many publications reported that miRNAs showed high diagnostic accuracy in distinguishing between patients with OS and healthy humans. Cao et al³⁵ noted that miR-326 showed 83.7% sensitivity and 94.5% specificity. Dong et al³⁶ reported the sensitivity and specificity for diagnostic accuracy of miR-223 to be 89.50 and 97.20%, respectively. Lian et al³⁷ revealed that the miRNA panel (miR-195-5p, miR-199a-3p, miR-320a, and miR-374a-5p) could provide a high accuracy of OS in patients, with 91.10% sensitivity, 94.40% specificity, and an AUC of 96.08%. A reasonable diagnostic accuracy in discerning patients with OS and healthy humans is 83.7% sensitivity and 94.5% specificity. The AUC value is an effective indicator for the assessment system. An ideal test with perfect discrimination is at an AUC of 1.0. As the AUC value of a test gets closer to 1.0, the overall efficacy of the test will increase. Through the collection of relevant literature, we found that miRNAs could be used to screen OS patients compared to healthy controls with an AUC value of 0.90, which is very close to 1.0. This suggests that miRNAs have a relatively high ability to determine OS patients from healthy humans.

However, at this time, no agreement has been reached among several miRNA expression–profiling studies. Different protocols, measurement platforms, and small sample sizes could result in incomparable expression levels. Thus, we conducted this meta-analysis to systematically evaluate the diagnostic value of miRNAs for OS diagnosis. In our study, the results showed that there were 17 unregulated and 8 downregulated miRNAs in OS patients. Additionally, we found that miRNA assays using serum specimens, with a specificity of 0.83 and DOR of 36.78, greatly improved the diagnostic accuracy compared to plasma-based assays, with a specificity of 0.91 and DOR of 86.05. Although the results are inconsistent with previous research, it was indicated that sample types might play a critical role in investigating the utility of miRNAs as biomarkers in disease diagnosis. There were some specific miRNAs that were expressed differently in different studies. Liu et al⁴⁰ reported that expression of miR-199a-3p was decreased in osteosarcoma patients compared to controls, whereas Lian et al³⁷ reported that expression of the same miRNA was significantly higher in osteosarcoma patients than controls. The reason may be that the blood samples were taken at different times. Liu et al reported that fasting blood samples were obtained from patients prior to surgery and/or chemotherapy, whereas Lian et al³⁷ reported that plasma samples were collected before surgery and 1 month after surgery. To confirm the best time for collecting the samples, future studies are needed. How the pathology, and response to radiotherapy and chemotherapy affect the expression of miRNAs and the specificity of these miRNA biomarkers should be further confirmed by comparing samples from patients with OS and patients with other types of cancers.

It was reported that miRNA expression profiles were considerably different between plasma and serum samples. Some researchers found that miRNA concentrations were higher in plasma.⁵⁰ We also found that miRNAs assays using serum specimen, with a specificity of 0.83 and DOR of 36.78, greatly improved the diagnostic accuracy over plasma-based assays, with a specificity of 0.93 and DOR of 86.05. Based on these results, it is concluded that sample types may play a critical role in investigating the utility of circulating miRNAs as biomarkers in disease diagnosis.

Detecting the altered expressions of miRNAs could have several advantages: 1) we can get stable miRNAs from fresh tissues (cells and peripheral blood);^51–53 2) each specific tissue is specifically expressed in miRNA; 3) miRNAs have been shown to be differentially expressed in patients with OS and normal humans; and 4) the expression of miRNAs combined with tumor stages has been determined.^54,55 Thus, continued research on miRNAs could enable us to find innovative biomarkers and to better understand the molecular mechanisms of OS.

Our meta-analysis has multiple merits compared with previous studies. First, we included numerous miRNAs examined in preceding studies, which has increased the reliability of our assessment of their overall diagnostic value, compared with previous meta-analyses.⁵⁶ Second, we conducted subgroup analyses to find potential sources of heterogeneity. Third, additional sensitivity analyses were utilized to confirm the strength of our results. However, in spite of these benefits, some limitations should be acknowledged. First, the literature sample size is small, which may decrease the level of evidence, so studies using larger sample sizes are needed to address the issue. Second, the biological functions of blood miRNAs are not very clear, which makes some data interpretation difficult. Third, although presented here as a potential diagnostic biomarker for OS, miRNAs possess numerous benefits and it is questionable whether these results can be implemented in clinical practice. Fourth, all of the data in this study are related to Asian countries, which is expected to result in heterogeneity. Finally, we could not obtain primary information regarding circulating miRNA expression data from these studies.

Conclusion

In conclusion, we found that many miRNAs play an important role in the early diagnosis of OS. Meanwhile, the miRNAs taken from plasma and serum samples could be promising noninvasive biomarkers for the early diagnosis of OS. Although the scheme seemed workable on paper, more experimental and clinical data need to be studied and collected step by step.

Disclosure

The authors report no conflicts of interest in this work.