ABSTRACT
Osteosarcoma is the most common aggressive primary bone tumor in children and adolescents. Angiogenesis, induced by vascular endothelial growth factor (VEGF) and its receptor (kinase insert domain receptor - KDR) is involved in tumor development. Both genes (VEGF and KDR) are polymorphic, and their association with osteosarcoma remains unclear. A systematic review of observational studies was performed to evaluate the association between polymorphisms in these genes and osteosarcoma development. Pubmed, Medline, Lilacs, and Scielo databases were searched for observational studies published up to April 2024. Eight publications of case-control studies were included, with quality ranging from 82% to 95%. All subjects were from the Chinese population: 1,681 cases and 2,049 controls. A total of six VEGF polymorphisms were analyzed. For osteosarcoma susceptibility, three studies found an increased risk with VEGF rs699947, four with VEGF rs2010963, two with VEGF rs3025039, one with VEGF rs833061, and no studies found an association with the VEGF rs1570360 and VEGF rs10434 SNPs. In addition, no studies evaluated SNPs in the KDR gene and osteosarcoma susceptibility. Further studies in diverse populations, particularly in Brazil, are necessary to clarify the role of VEGF and KDR polymorphisms in osteosarcoma development and prognosis. Level of evidence III; Systematic Review.
Keywords:
Osteosarcoma; Genetic Polymorphism; Angiogenesis; VEGF; Systematic Review
RESUMO
O osteossarcoma é o tumor ósseo primário agressivo mais comum em crianças e adolescentes. A angiogênese, induzida pelo fator de crescimento endotelial vascular (VEGF) e seu receptor domínio quinase (KDR), está envolvida no desenvolvimento tumoral. Ambos os genes (VEGF e KDR) são polimórficos, e sua associação com o osteossarcoma permanece incerta. Uma revisão sistemática de estudos observacionais foi realizada para avaliar a associação de polimorfismos nesses genes com o desenvolvimento de osteossarcoma. As bases de dados Pubmed, Medline, Lilacs e Scielo foram consultadas em busca de estudos observacionais publicados até abril de 2024. Foram incluídas oito publicações de estudos de caso-controle, com qualidade variando de 82% a 95%. Todos os participantes eram da população chinesa: 1.681 casos e 2.049 controles. Foram analisados seis polimorfismos do VEGF. Três estudos encontraram um risco aumentado de desenvolvimento do osteossarcoma com VEGF rs699947, quatro com VEGF rs2010963, dois com VEGF rs3025039, um com VEGF rs833061, e nenhum estudo encontrou associação com os SNPs VEGF rs1570360 e VEGF rs10434. Nenhum estudo avaliou SNPs no gene KDR. Estudos adicionais em populações diversas, especialmente no Brasil, são necessários para esclarecer o papel dos polimorfismos nos genes VEGF e KDR no desenvolvimento e também no prognóstico do osteossarcoma. Nível de Evidência III; Revisão Sistemática.
Descritores:
Osteossarcoma; Polimorfismo Genético; Angiogênese; Fator de Crescimento Endotelial Vascular; Revisão Sistemática
INTRODUCTION
Osteosarcoma is a rare malignant primary bone tumor that usually affects long bones such as the distal femur, tibia, and proximal humerus.1,2 Its age distribution is bimodal, with a higher peak between 10 and 20 years of age and a second peak between 70 and 80 years of age.1 In childhood, the incidence is higher in females, but between 15 and 19 years of age, it becomes higher in males.3 The most common symptoms of osteosarcoma are pain and local swelling. The bone becomes fragile due to the tumor presence,4 and approximately 17% of patients present with pathologic fractures.5 Imaging and laboratory tests aid in the diagnostic investigation; however, the definitive diagnosis is confirmed by histologic evaluation.6 The standard treatment consists of neoadjuvant and adjuvant chemotherapy combined with surgery.7 Survival rates for patients with localized disease range from 60-75%.8-10 On average, 20% of patients have pulmonary metastases at diagnosis.11 The presence of metastases significantly decreases overall survival in 10 years, approximately 12.5% in Brazil11 and 25% in the United States and Europe.12
The mechanism of osteosarcoma development is complex, and its exact etiology is still unknown. However, environmental and genetic factors have been associated with the development of this type of tumor.13 Additionally, angiogenesis is a crucial process in the development, and proliferation of tumor cells, cancer progression, and metastasis development.14 Vascular endothelial growth factor (VEGF) is one of the main pro-angiogenic factors, and high expression of VEGF and its kinase insert domain receptor (KDR) has been observed in osteosarcoma tissue samples compared to adjacent tissues.15,16 VEGF expression was significantly correlated with histological grade, tumor stage, distant metastasis, and lower survival.17-19 Positive VEGF status in tumor samples from osteosarcoma patients negatively impacts overall survival and disease-free survival, suggesting this factor is an effective biological marker for prognosis.20,21 In this regard, anti-angiogenic inhibitor drugs targeting VEGF/VEGFR successfully inhibited osteosarcoma angiogenesis and cell proliferation in vivo and in vitro,22 supporting the relevance of VEGF and its receptor in osteosarcoma.
The KDR and VEGF genes are located on chromosomes 4q11-q12 and 6p21.1, respectively, and have several single nucleotide polymorphisms (SNPs) in the coding and regulatory regions that may affect enzyme activity or expression.23-25 In addition, different populations show increased frequencies of these SNPs.26-28 Additionally, VEGF SNPs have been associated with osteosarcoma risk in the Chinese population29-34 although the literature reports show controversial results. Different isoforms of VEGF have been associated with osteosarcoma progression,35 metastases development,36 and poor response to chemotherapy.37 Among these, the most important isoform involved in angiogenesis and vasculogenesis is VEGF-A,38 which is the subject of the current study. In this context, the aim of this study was to perform a systematic review of observational studies evaluating VEGF and KDR SNPs and osteosarcoma susceptibility to discuss the causes of the previously described conflicting results.
METHODS
Systematic Review Search Strategy
A literature search was conducted in the Pubmed, Medline, Lilacs, and Scielo databases to identify all articles that evaluated SNPs in the VEGF and KDR genes associated with the risk of developing osteosarcoma. These last two databases were consulted with specific interest in studies involving Latin American populations. All studies published up to April 2024 were evaluated, considering the following combined descriptors in both English and Portuguese: osteosarcoma and (“polymorphism” or “SNP” or “genetic polymorphism”) and (“VEGF” or “Vascular endothelial growth factor” or “VEGFR-2” or “Vascular endothelial growth factor-2” or “KDR” or “Kinase Insert Domain Receptor”). Additionally, the reference list of review/meta-analysis articles and of each selected article were reviewed to identify any other articles that might not have been included in the initial data search.
Inclusion and Exclusion Criteria
The articles selected for this review followed the following inclusion criteria: (i) observational studies that evaluated the association of SNPs in the VEGF and KDR genes with the development of osteosarcoma; (ii) studies published in English or Portuguese; (iii) studies published up to April 2024. The exclusion criteria were: (i) publications that did not include osteosarcoma patients, (ii) did not analyze VEGF and KDR SNPs; (iii) analyzed only disease progression; and (iv) publications without completely accessible texts.
Data Extraction
After selecting the studies according to the inclusion and exclusion criteria, the following information was extracted from each article by two reviewers (JVC and JPN): first author and year of publication; study population; the number of osteosarcoma cases and controls; the age of both groups; gender; type of control recruited; patient-reported family history of cancer; location of the osteosarcoma; presence or absence of metastasis; clinical staging of the osteosarcoma cases; identification of SNPs; genotyping technique for SNP identification; SNP allele frequency data; and Pearson’s chi-square test and odds ratio (OR) values with their respective 95% confidence intervals (CI 95%).
Data management
The data extracted from the studies included in this review were organized and analyzed using Excel spreadsheets. Percentages were used as statistical measures to compare the frequency of each variable extracted from the studies, followed by a descriptive analysis of the data. To aid in data synthesis and visualization, Excel tools were employed for basic calculations such as sums and percentages, as well as for arranging the characteristics of each study. All data were reviewed by two researchers (JVC and JPN). When an article did not report the frequency of the alleles for the polymorphisms included, the frequencies were calculated based on genotype information. Since each individual has two alleles, the total number of alleles was determined by summing the set of genotypes.
Quality Assessment of Included Studies
All included studies were independently analyzed by two reviewers (JVC and JPN) according to the STROBE checklist used for observational studies. This quality assessment tool has twenty-two items related to: title and abstract (item 1), introduction (items 2 and 3), methodology (items 4 to 12), results (items 13 to 17), discussion (items 18 to 21), and funding of each study (item 22), assigning a score of 0 or 1 for each item. In the end, studies with a score greater than 50% were approved for evaluation in this review according to the guideline.39
RESULTS
The flowchart of the included articles is illustrated in Figure 1. Out of the 206 publications found, we excluded 189 duplicates and were left with 17 valid publications. One additional article31 was identified after reviewing the reference list of the 17 selected articles, making a total of 18 articles for evaluation of eligibility criteria. Subsequently, 10 articles were excluded: 4 meta-analyses,40-43 3 studies that analyzed only disease progression,16,44,45 2 studies that did not analyze SNPs in the VEGF or KDR genes,46,47 and one study that did not include osteosarcoma patients.48 Consequently, eight publications29-34,49,50 were selected for full-text evaluation and included in this review. All eight studies investigated SNPs in the VEGF gene only. No publication was found that described an association of SNPs in the KDR gene with osteosarcoma development (Figure 1).
Table 1 describes the characteristics of osteosarcoma cases from the studies included in this review. Regarding the STROBE rating of the included articles, the studies had a quality score of over 80%, with Wang et al.50 having the highest score (21 out of 22 items) and the studies by Hu et al.29 Zhang et al.33 and Zhang et al.34 having the lowest scores (18 out of 22 items).
The study population in all articles was Chinese, with a total of 1,681 osteosarcoma cases analyzed, with the majority (75%) of articles presenting more individuals under the age of 20 years,29,31-33,49,50 and males being the most common sex in all articles (100%). Only the studies by Hu et al.29 and Liu et al.31 did not include information on family history of cancer (25%). Among the other studies, four (50%) observed a high frequency of cancer in the family30,32-34 and two (25%) reported a low frequency of this characteristic.49,50 In all articles, osteosarcoma was more frequently located in the long bones, and only in the article by Liu et al.31 about 77% of patients had metastasis, while the other studies (n = 6; 75%) observed a much lower frequency (ranging from 19.4 to 34.6%).29,32-34,40,50 The study by Li-Lian et al.30 did not include this information. In addition, four (50%) studies29,30,32,33 reported stage I-II being the most common. Wang et al.50 Liu et al.31 Zhang et al.34 and Cao et al.49 did not report the tumor staging, although Liu et al.31 reported collecting this information along with histologic subtype (osteoblastic, chondroblastic, fibroblastic, and mixed types). However, the frequency of both characteristics was not found in the text.
The demographic characteristics and inclusion and exclusion criteria of the control participants from the selected studies (a total of 2,049 controls) are also described in this review (Table 2). The study period of the included articles ranged from 2008 to 2015, with seven being hospital-based (87.5%) and one being population-based.50 For controls (healthy individuals), most studies (n = 7; 87.5%) reported a negative personal history of cancer and matched controls to cases by age and sex. In the article by Cao et al.49 the authors only reported a negative personal history of cancer for healthy individuals. In the study by Wang et al.50 controls were also matched by area of residence. Li-Lian et al.30 and Zhang et al.34 also reported a negative personal history of chronic diseases, while only Liu et al.31 reported that the controls had no family history of cancer. In most articles (75%), the controls were under 20 years of age,9,31-34,49,50 with a higher frequency of males observed in all studies. Four studies (50%) reported a high frequency of family history of cancer in the control group.30,32-34 This characteristic was less frequent only in the studies of Wang et al.50 and Cao et al.49 Hu et al.29 did not provide this information.
Figure 2 displays the frequencies of the variant alleles of each VEGF gene SNP found in the eligible studies of this systematic review. For the rs699947 C>A SNP, the variant allele frequency ranged from 37.4% to 42.6% in the case group and from 31.2% to 33.6% in the control group. Of the four published articles on this SNP,30-32,34 three30,32,34 found a positive association with the development of osteosarcoma (Figure 2A; Table 3). For the rs1570360 A>G SNP, the variant allele frequency ranged from 26.8% to 31.7% in cases and from 24.5% to 26.2% in controls, with no significant differences observed between groups observed in any study (Figure 2B; Table 3).30-32 Regarding the rs833061 T>C SNP, the variant allele frequency ranged from 39.3% to 47.4% in cases and from 36.8% to 38.3% in controls, with only Li-Lian et al.30 finding an increased risk of developing osteosarcoma with the CC genotype or the dominant model (TC+CC) (Figure 2C; Table 3). Seven articles evaluated the rs2010963 G>C SNP,29-33,49,50 with four reporting a positive association with disease development in the presence of the G allele or GG genotype.31-33,49 The variant allele frequency varied from 40.2% to 58.8% in the case group and from 39.1% to 61.8% in the control group across published studies (Figure 2D; Table 3).
Frequency of variant alleles of the SNPs VEGF rs699947 A (A), VEGF rs1570360 G (B), VEGF rs833061 C (C), VEGF rs2010963 C (D), VEGF rs3025039 T (E), and VEGF rs10434 A (F) found in the studies included in this review.
Zhang H. et al. 2015;34 Zhang G. et al, 2015.33 * P-value of Person's Chi-square test < 0.05.
All articles in this review examined the rs3025039 C>T SNP, but only Wang et al.50 and Hu et al.29 observed an approximately twofold increased risk of developing osteosarcoma. The frequency of the rs3025039 T variant allele ranged from 20.1% to 37.8% in cases and from 15.4% to 34.6% in controls (Figure 2E; Table 3). Finally, six articles examined the rs10434 G>A SNP,29-33,50 and none of them found a significant association with tumor development. The frequency of the variant allele ranged from 33.5% to 47.4% in cases and from 32.2% to 46.5% in controls (Figure 2F; Table 3). It is worth noting that, of the articles included in this review, only one evaluated both the impact of polymorphisms on the risk of developing osteosarcoma and their prognostic value. The authors31 associated the polymorphism rs2010963 (-634GG) with a shorter overall survival (OS) time (HR, 3.10; 95% CI, 1.17-8.38).
Table 3 details the characteristics of the VEGF gene SNPs and the data from the association analyses of each study. The eight studies included in this review, identified six VEGF gene variants: rs699947 C>A, rs1570360 A>G, rs833061 T>C, rs2010963 G>C, rs3025039 C>T, and rs10434 G>A. The first three variants are located in the promoter region, one in the 5’ UTR, and the remaining two in the 3’ UTR, respectively. All studies used the same genotyping technique for these SNPs, polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP).
Additionally, Li-Lian et al.30 Tie et al.32 and Zhang G. et al.33 performed stratification analyses based on environmental characteristics for the SNPs in their respective studies. Li-Lian et al.30 found no associations for any variables, whereas Tie et al.32 reported that the rs699947 AA and rs2010963 GG genotypes were associated with an increased risk of osteosarcoma in individuals younger than 20 years and those with a family history of cancer. Furthermore, the rs699947 and rs2010963 genotypes were associated with an increased risk of developing osteosarcoma in males and females, respectively. Zhang et al.33 also found that the rs2010963 GG genotype was associated with an increased risk in individuals younger than 20 years, in both sexes and with a family history of cancer. In addition, Liu et al.31 evaluated the survival of osteosarcoma cases and found that the rs2010963 SNP could be considered as a prognostic factor for the disease, with the rs2010963 genotype being associated with a threefold higher risk of decreased survival (OR = 3.1; 95% CI = 1.17-8.38).
DISCUSSION
This study conducted a systematic review of observational studies investigating the potential relationship between SNPs in the VEGF and KDR genes and the risk of developing osteosarcoma. However, at the time of the literature review, none of the published studies had evaluated the association between SNPs in the KDR gene and osteosarcoma susceptibility. Eight studies examined the interaction of VEGF SNPs and osteosarcoma susceptibility,29-34,49,50 and only four SNPs (rs699947, rs833061, rs2010963, and rs3025039) were associated with osteosarcoma risk.29-34,49,50 Furthermore, one VEGF SNP (rs2010963) was associated with shorter survival.31
In most articles, osteosarcoma cases involved men under the age of 20, which is consistent with other studies conducted in different populations.1,51 A study by our group of Brazilian patients with osteosarcoma recruited from a reference orthopedic hospital included samples from 37 histologically confirmed patients (24 males and 13 females) with a mean age of 18 years,52 which is consistent with the current systematic review. A family history of cancer was highly prevalent (87.6-91.8% and 91.8-96.9%) among cases and controls, respectively, in four of the selected studies.30,32-34 Various studies have found that a family history of cancer increases the likelihood of developing endometrial and prostate cancer by 3.1% and 59.7%, respectively, and that 15-20% of breast cancer cases are associated with a family history of cancer.53,54 For osteosarcoma, a 3.6-fold risk was observed, associated with maternal endocrine gland cancer.55
In all reviewed articles, osteosarcoma was more commonly found in long bones, which is consistent with the literature.1,51 The growth plate has a higher rate of cell proliferation, making cells more susceptible to mutations.51 Similarly, our group’s study showed that most tumors were located in the knee region (81%) and the most common histologic variant was high-grade central conventional osteosarcoma (83.8%), with 27% presenting metastases.52 Most patients in this review did not have metastases and had stage I-II tumors, corroborating other studies.56,57 Advanced age, axial skeletal location, lesion size, and residence in disadvantaged areas are risk factors associated with the presence of metastases at diagnosis.56
The angiogenic process involving VEGF and its receptors plays an important role in bone tissue regeneration and in the growth, progression, and metastasis of bone tumors.58 Studies suggest an imbalance between angiogenic inducers and inhibitors, leading to an angiogenic tumor phenotype and tumor development.15,59 In osteosarcoma, VEGF expression is associated with poor prognosis and metastasis development.60 In vitro, VEGF silencing inhibits proliferation and promotes apoptosis in SaOS-2 cells, while in vivo, it suppresses tumor growth and angiogenesis.61 Zheng et al.62 showed higher expression of VEGF receptor (VEGFR2) in metastatic tumor samples (lung) compared to the primary tumor, and receptor inhibition reduces metastatic ability in vitro and significantly attenuates metastatic capacity in vivo.
In 2016, a meta-analysis of seven articles to provide a comprehensive assessment of the associations of SNPs in the VEGF gene with osteosarcoma susceptibility found that the VEGF rs3025039 C>T and rs2010963 G>C SNPs were associated with increased risk and protection, respectively, and that the rs10434 G>A SNP had no association with disease development.40 Then, in 2017, another meta-analysis observed that the VEGF rs699947 C>A, rs1570360 A>G, rs833061 T>C, rs2010963 G>C, and rs3025039 C>T SNPs were associated with osteosarcoma risk in the Han Chinese population, while the VEGF rs10434 G>A SNP showed no significant association with the disease.43 Recently, in 2022, another meta-analysis highlighted that the VEGF rs1570360 A>G, rs2010963 G>C, and rs3025039 C>T SNPs were associated with osteosarcoma risk, while the C allele of the VEGF rs699947 C>A SNP had a protective effect.42 More recently, a meta-analysis that included the analysis of five SNPs, four of which were also evaluated in this review. The data from Hassanain et al.63 support the risk analysis for the VEGF rs3025039 TT and rs2010963 GG SNPs. In contrast to our results, the leave-one-out analysis showed that the VEGF rs699947, rs10434, and rs1570360 SNPs were not significant. The VEGF rs833061 CC SNP was excluded from the meta-analysis because of the low quality of the articles. However, in our review, the association (2-fold risk) was based on the study by Li-Lian et al.30 which was not cited in the meta-analysis by Hassanain et al. Similarly, contrary to our findings, the authors reported a risk associated with the VEGF rs157090360 based on Zhao’s meta-analysis,43 which was not included in this review once we excluded meta-analysis studies.
The VEGF rs699947 C>A, rs1570360 A>G, and rs833061 T>C SNPs may affect transcriptional activity and lead to increased protein levels because they are located in the promoter region of the gene.64 The VEGF rs2010963 G>C SNP, located in the 5’UTR, is significantly positively correlated with VEGF protein synthesis by peripheral blood mononuclear cells (PBMCs).24 Meanwhile, the VEGF rs3025039 C>T SNP, located in the 3’UTR, may influence plasma VEGF levels.23 For VEGF rs10434 G>A SNP, also located in the 3’UTR, no studies have yet confirmed its influence on the gene; however, this type of SNP may affect mRNA stability and alter the ability of microRNAs to interact with their target due to its position in the gene.65
Although most of the studies included in this review reported similar frequencies of VEGF SNPs, they were all conducted in the Chinese population, which could lead to the inconsistencies in the association results. These inconsistencies may also be explained by (i) variability in sample size, which may affect risk estimates, and (ii) the type of control used, which may also be associated with the variants studied. In case-control studies, the selection of the control group is critical, as these individuals should not have the chance of developing the outcome or have risk factors similar to the outcome, such as SNPs in the VEGF gene, which could bias risk estimates. Most articles in this review recruited healthy hospital-based controls matched to the case group by age and sex and with a negative personal history of cancer, which is a positive aspect as the ideal control group is directly determined by the definition and selection of the case group and recruited from the same source.66 A limitation of the present review, common to other literature reviews, is the possibility that we did not include all the literature and missed some relevant studies. We included three databases, English-language articles only, and did not include grey literature.
Comparing the frequencies of VEGF rs699947 A, rs1570360 G, rs833061 C, and rs2010963 C SNPs in populations from different regions worldwide (according to data from the National Center for Biotechnology Information - NCBI), with the Chinese population (from this review) and the Brazilian population,67 which has an extensive admixture of three major ancestral roots (indigenous, European and African) and is predominant in today’s Brazilian population.68 For VEGF rs699947 A SNP, the frequencies in Chinese (31.2-33.6%), Asian (28-34%) (NCBI), and Brazilian (31.9%) populations67 were found to differ from those in African American (23%) and Sub-Saharan African (12-17%) populations (NCBI). For VEGF rs1570360 G SNP, the frequencies in European (67-69%), African (91-97%), American (78-80%), Latin American (75-80%) (NCBI), and Brazilian populations (83.8%)67 differ from those in the Chinese population (24.5-26.2%). For the VEGF rs833061 C SNP, the frequencies in European (47%) (NCBI) and Brazilian (41.7%) populations67 differ from those in Asian (27-34%), Chinese (38.3-36.8%), African American (33%), and Sub-Saharan African (30%) populations (NCBI). Finally, for the VEGF rs2010963 C and rs3025039 SNPs, the frequencies in European (27-42%; 10-23%, respectively), African American (24-36%; 5-12%, respectively) (NCBI), and Brazilian (38.6% and 12.5%, respectively) populations67 differed from those in the Chinese population (39.1-61.8%; 15.4-34.6%, respectively).
The inconsistencies in SNP frequencies reinforce the idea that the variants studied in this review may exist at variable frequencies in both more homogeneous populations, such as the Chinese, and more heterogeneous populations, such as the Brazilians. This phenomenon may occur because of factors such as long-distance migration and the encounter of different populations during demographic expansion, resulting in genetic mixing or complete population replacement.69,70 In addition, genetic interactions with environmental factors may influence these frequencies.71 Therefore, extrapolation of data from Chinese SNP association studies to mixed populations is likely to lead to conclusions that are inconsistent with reality. In addition, the paucity of studies evaluating SNPs in the KDR gene identified in this review highlights the need for further investigation of this gene, given its importance in VEGF signaling and consequently in the angiogenic process.
CONCLUSIONS
For osteosarcoma susceptibility, three studies found an increased risk with VEGF rs699947, four with VEGF rs2010963, two with VEGF rs3025039, one with VEGF rs833061, and no studies found an association with the VEGF rs1570360 and VEGF rs10434 SNPs. In addition, no studies evaluated SNPs in the KDR gene and osteosarcoma susceptibility.
REFERENCES
-
1 Sadykova LR, Ntekim AI, Muyangwa-Semenova M, Rutland CS, Jeyapalan JN, Blatt N, Rizvanov AA. Epidemiology and Risk Factors of Osteosarcoma. Cancer Invest. 2020;38(5):259-269. doi: 10.1080/07357907.2020.1768401.
» https://doi.org/10.1080/07357907.2020.1768401 -
2 Aran V, Devalle S, Meohas W, Heringer M, Cunha Caruso A, Pinheiro Aguiar D, Leite Duarte ME, Moura Neto V. Osteosarcoma, chondrosarcoma and Ewing sarcoma: Clinical aspects, biomarker discovery and liquid biopsy. Crit Rev Oncol Hematol. 2021;162:103340. doi: 10.1016/j.critrevonc.2021.103340.
» https://doi.org/10.1016/j.critrevonc.2021.103340 -
3 Balmant NV, Reis RS, Santos MO, Maschietto M, Camargo B. Incidence and mortality of bone cancer among children, adolescents and young adults of Brazil. Clinics (Sao Paulo). 2019;74:e858. doi: 10.6061/clinics/2019/e858.
» https://doi.org/10.6061/clinics/2019/e858 -
4 Taran SJ, Taran R, Malipatil NB. Pediatric Osteosarcoma: An Updated Review. Indian J Med Paediatr Oncol. 2017;38(1):33-43. doi: 10.4103/0971-5851.203513.
» https://doi.org/10.4103/0971-5851.203513 - 5 Scully SP, Ghert MA, Zurakowski D, Thompson RC, Gebhardt MC. Pathologic fracture in osteosarcoma: prognostic importance and treatment implications. J Bone Joint Surg Am. 2002;84(1):49-57.
-
6 Abarrategi A, Tornin J, Martinez-Cruzado L, Hamilton A, Martinez-Campos E, Rodrigo JP, et al. Osteosarcoma: Cells-of-Origin, Cancer Stem Cells, and Targeted Therapies. Stem Cells Int. 2016;2016:3631764. doi: 10.1155/2016/3631764.
» https://doi.org/10.1155/2016/3631764 -
7 Ando K, Heymann MF, Stresing V, Mori K, Rédini F, Heymann D. Current therapeutic strategies and novel approaches in osteosarcoma. Cancers (Basel). 2013;5(2):591-616. doi: 10.3390/cancers5020591.
» https://doi.org/10.3390/cancers5020591 -
8 Bielack S, Jürgens H, Jundt G, Kevric M, Kühne T, Reichardt P, et al. Osteosarcoma: the COSS experience. Cancer Treat Res. 2009;152:289-308. doi: 10.1007/978-1-4419-0284-9_15.
» https://doi.org/10.1007/978-1-4419-0284-9_15 -
9 Li X, Ashana AO, Moretti VM, Lackman RD. The relation of tumour necrosis and survival in patients with osteosarcoma. Int Orthop. 2011;35(12):1847-53. doi: 10.1007/s00264-011-1209-7.
» https://doi.org/10.1007/s00264-011-1209-7 -
10 Jaffe N. Historical perspective on the introduction and use of chemotherapy for the treatment of osteosarcoma. Adv Exp Med Biol. 2014;804:1-30. doi: 10.1007/978-3-319-04843-7_1.
» https://doi.org/10.1007/978-3-319-04843-7_1 -
11 Petrilli AS, de Camargo B, Filho VO, Bruniera P, Brunetto AL, Jesus-Garcia R, et al. Results of the Brazilian Osteosarcoma Treatment Group Studies III and IV: prognostic factors and impact on survival. J Clin Oncol. 2006;24(7):1161-8. doi: 10.1200/JCO.2005.03.5352.
» https://doi.org/10.1200/JCO.2005.03.5352 -
12 Bielack SS, Kempf-Bielack B, Delling G, Exner GU, Flege S, Helmke K, et al. Prognostic factors in high-grade osteosarcoma of the extremities or trunk: an analysis of 1,702 patients treated on neoadjuvant cooperative osteosarcoma study group protocols. J Clin Oncol. 2002;20(3):776-90. doi: 10.1200/JCO.2002.20.3.776.
» https://doi.org/10.1200/JCO.2002.20.3.776 -
13 Longhi A, Pasini A, Cicognani A, Baronio F, Pellacani A, Baldini N, et al. Height as a risk factor for osteosarcoma. J Pediatr Hematol Oncol. 2005;27(6):314-8. doi: 10.1097/01.mph.0000169251.57611.8e.
» https://doi.org/10.1097/01.mph.0000169251.57611.8e -
14 Kim KJ, Li B, Winer J, Armanini M, Gillett N, Phillips HS, Ferrara N. Inhibition of vascular endothelial growth factor-induced angiogenesis suppresses tumour growth in vivo. Nature. 1993;362(6423):841-4. doi: 10.1038/362841a0.
» https://doi.org/10.1038/362841a0 -
15 Rastogi S, Kumar R, Sankineani SR, Marimuthu K, Rijal L, Prakash S, et al. Role of vascular endothelial growth factor as a tumour marker in osteosarcoma: a prospective study. Int Orthop. 2012;36(11):2315-21. doi: 10.1007/s00264-012-1663-x.
» https://doi.org/10.1007/s00264-012-1663-x -
16 Liu JY, Zhu BR, Wang YD, Sun X. The efficacy and safety of Apatinib mesylate in the treatment of metastatic osteosarcoma patients who progressed after standard therapy and the VEGFR2 gene polymorphism analysis. Int J Clin Oncol. 2020;25(6):1195-1205. doi: 10.1007/s10147-020-01644-7.
» https://doi.org/10.1007/s10147-020-01644-7 -
17 Han G, Wang Y, Bi W, Jia J, Wang W, Xu M. Effects of vascular endothelial growth factor expression on pathological characteristics and prognosis of osteosarcoma. Clin Exp Med. 2016;16(4):577-584. doi: 10.1007/s10238-015-0382-1.
» https://doi.org/10.1007/s10238-015-0382-1 -
18 Wu H, Zhang J, Dai R, Xu J, Feng H. Transferrin receptor-1 and VEGF are prognostic factors for osteosarcoma. J Orthop Surg Res. 2019;14(1):296. doi: 10.1186/s13018-019-1301-z.
» https://doi.org/10.1186/s13018-019-1301-z -
19 Mohamed FEA, Khalil EZI, Toni NDM. Caveolin-1 Expression Together with VEGF can be a Predictor for Lung Metastasis and Poor Prognosis in Osteosarcoma. Pathol Oncol Res. 2020;26(3):1787-1795. doi: 10.1007/s12253-019-00755-5.
» https://doi.org/10.1007/s12253-019-00755-5 -
20 Chen D, Zhang YJ, Zhu KW, Wang WC. A systematic review of vascular endothelial growth factor expression as a biomarker of prognosis in patients with osteosarcoma. Tumour Biol. 2013;34(3):1895-9. doi: 10.1007/s13277-013-0733-z.
» https://doi.org/10.1007/s13277-013-0733-z -
21 Xin S, Wei G. Correlation of vascular endothelial growth factor with survival and pathological characteristics of patients with osteosarcoma: A systematic review and meta-analysis. Eur J Cancer Care (Engl). 2022;31(5):e13629. doi: 10.1111/ecc.13629.
» https://doi.org/10.1111/ecc.13629 -
22 Wang BD, Yu XJ, Hou JC, Fu B, Zheng H, Liu QK, et al. Bevacizumab attenuates osteosarcoma angiogenesis by suppressing MIAT encapsulated by serum-derived extracellular vesicles and facilitating miR-613-mediated GPR158 inhibition. Cell Death Dis. 2022;13(3):272. doi: 10.1038/s41419-022-04620-3.
» https://doi.org/10.1038/s41419-022-04620-3 -
23 Renner W, Kotschan S, Hoffmann C, Obermayer-Pietsch B, Pilger E. A common 936 C/T mutation in the gene for vascular endothelial growth factor is associated with vascular endothelial growth factor plasma levels. J Vasc Res. 2000;37(6):443-8. doi: 10.1159/000054076.
» https://doi.org/10.1159/000054076 -
24 Watson CJ, Webb NJ, Bottomley MJ, Brenchley PE. Identification of polymorphisms within the vascular endothelial growth factor (VEGF) gene: correlation with variation in VEGF protein production. Cytokine. 2000;12(8):1232-5. doi: 10.1006/cyto.2000.0692.
» https://doi.org/10.1006/cyto.2000.0692 -
25 Wang Y, Zheng Y, Zhang W, Yu H, Lou K, Zhang Y, et al. Polymorphisms of KDR gene are associated with coronary heart disease. J Am Coll Cardiol. 2007;50(8):760-7. doi: 10.1016/j.jacc.2007.04.074.
» https://doi.org/10.1016/j.jacc.2007.04.074 - 26 Cardoso JV, Machado DE, Ferrari R, Silva MC, Berardo PT, Perini JA. Polymorphisms in VEGF and KDR Genes in the Development of Endometriosis: A Systematic Review. Rev. Bras. Saude Mater. Infant. 2016;16(3):219-232. doi:10.1590/1806-93042016000300002.
-
27 Wang J, Shen C, Fu Y, Yu T, Song J. The associations between five polymorphisms of vascular endothelial growth factor and renal cell carcinoma risk: an updated meta-analysis. Onco Targets Ther. 2017;10:1725-1734. doi: 10.2147/OTT.S125965.
» https://doi.org/10.2147/OTT.S125965 -
28 Li XY, Wang YL, Yang S, Liao CS, Li SF, Han PF. Correlation between vascular endothelial growth factor A gene polymorphisms and tendon and ligament injury risk: a systematic review and meta-analysis. J Orthop Surg Res. 2024;19(1):122. doi: 10.1186/s13018-024-04589-z.
» https://doi.org/10.1186/s13018-024-04589-z -
29 Hu GL, Ma G, Ming JH. Impact of common SNPs in VEGF gene on the susceptibility of osteosarcoma. Genet Mol Res. 2015;14(4):14561-6. doi: 10.4238/2015.November.18.19.
» https://doi.org/10.4238/2015.November.18.19 -
30 Hu GL, Ma G, Ming JH. Impact of common SNPs in VEGF gene on the susceptibility of osteosarcoma. Genet Mol Res. 2015;14(4):14561-6. doi: 10.4238/2015.November.18.19.
» https://doi.org/10.4238/2015.November.18.19 -
31 Liu JQ, Bai X, Duan DC, Dou AX. Role of five small nucleotide polymorphisms in the VEGF gene on the susceptibility to osteosarcoma and overall survival of patients. Oncol Lett. 2015;10(3):1481-1486. doi: 10.3892/ol.2015.3396.
» https://doi.org/10.3892/ol.2015.3396 - 32 Tie Z, Bai R, Zhai Z, Zhang G, Zhang H, Zhao Z, et al. Single nucleotide polymorphisms in VEGF gene are associated with an increased risk of osteosarcoma. Int J Clin Exp Pathol. 2014;7(11):8143-9.
-
33 Zhang G, Bai R, Zhang T, Zhang H, Wen SZ, Jiang DM. Investigation of the role of VEGF gene polymorphisms in the risk of osteosarcoma. Genet Mol Res. 2015;14(3):8283-9. doi: 10.4238/2015.July.27.16.
» https://doi.org/10.4238/2015.July.27.16 -
34 Zhang HF, Yan JP, Zhuang YS, Han GQ. Association between angiogenic growth factor genetic polymorphisms and the risk of osteosarcoma. Genet Mol Res. 2015;14(3):10524-9. doi: 10.4238/2015.September.8.14.
» https://doi.org/10.4238/2015.September.8.14 -
35 Matsumoto I, Oda M, Yachi T, Tsuchiya H, Zen Y, Watanabe G. Outcome prediction of pulmonary metastasectomy can be evaluated using metastatic lesion in osteosarcoma patients. World J Surg. 2013;37(8):1973-80. doi: 10.1007/s00268-013-2022-9.
» https://doi.org/10.1007/s00268-013-2022-9 -
36 Fan H, Lu S, Wang S, Zhang S. Identification of critical genes associated with human osteosarcoma metastasis based on integrated gene expression profiling. Mol Med Rep. 2019;20(2):915-930. doi: 10.3892/mmr.2019.10323.
» https://doi.org/10.3892/mmr.2019.10323 -
37 Abdeen A, Chou AJ, Healey JH, Khanna C, Osborne TS, Hewitt SM, et al. Correlation between clinical outcome and growth factor pathway expression in osteogenic sarcoma. Cancer. 2009;115(22):5243-50. doi: 10.1002/cncr.24562.
» https://doi.org/10.1002/cncr.24562 -
38 Assi T, Watson S, Samra B, Rassy E, Le Cesne A, Italiano A, Mir O. Targeting the VEGF Pathway in Osteosarcoma. Cells. 2021;10(5):1240. doi: 10.3390/cells10051240.
» https://doi.org/10.3390/cells10051240 -
39 von Elm E, Altman DG, Egger M, Pocock SJ, Gøtzsche PC, Vandenbroucke JP; STROBE Initiative. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. BMJ. 2007;335(7624):806-8. doi: 10.1136/bmj.39335.541782.AD.
» https://doi.org/10.1136/bmj.39335.541782.AD -
40 Hu YY, Du XY, Zhan AL, Zhou L, Jiang Q, Niu YM, et al. Vascular endothelial growth factor polymorphisms are associated with osteosarcoma susceptibility. Oncotarget. 2016;7(30):47711-47719. doi: 10.18632/oncotarget.10278.
» https://doi.org/10.18632/oncotarget.10278 -
41 Wang X, Liu Z. Systematic meta-analysis of genetic variants associated with osteosarcoma susceptibility. Medicine (Baltimore). 2018;97(38):e12525. doi: 10.1097/MD.0000000000012525.
» https://doi.org/10.1097/MD.0000000000012525 -
42 Yuan D, Tian J, Fang X, Xiong Y, Banskota N, Kuang F, et al. Epidemiological Evidence for Associations Between Genetic Variants and Osteosarcoma Susceptibility: A Meta-Analysis. Front Oncol. 2022;12:912208. doi: 10.3389/fonc.2022.912208.
» https://doi.org/10.3389/fonc.2022.912208 -
43 Yuan D, Tian J, Fang X, Xiong Y, Banskota N, Kuang F, et al. Epidemiological Evidence for Associations Between Genetic Variants and Osteosarcoma Susceptibility: A Meta-Analysis. Front Oncol. 2022;12:912208. doi: 10.3389/fonc.2022.912208.
» https://doi.org/10.3389/fonc.2022.912208 -
44 Dong-Ju Z, Ai-Ju X, Yun-Jiao T, Ming-Qiu Z. Polymorphisms of vascular endothelial growth factor on prognosis in osteosarcoma patients. Pak J Med Sci. 2014;30(5):1072-6. doi: 10.12669/pjms.305.5170.
» https://doi.org/10.12669/pjms.305.5170 -
45 Li XD, Yu DD, Lu JR, Wu C, Jin WX. Effect of polymorphisms of vascular endothelial growth factor on prognosis in osteosarcoma patients. Genet Mol Res. 2015;14(2):4354-60. doi: 10.4238/2015.April.30.8.
» https://doi.org/10.4238/2015.April.30.8 -
46 Negri GL, Grande BM, Delaidelli A, El-Naggar A, Cochrane D, Lau CC, et al. Integrative genomic analysis of matched primary and metastatic pediatric osteosarcoma. J Pathol. 2019;249(3):319-331. doi: 10.1002/path.5319.
» https://doi.org/10.1002/path.5319 -
47 Sasa K, Saito T, Kurihara T, Hasegawa N, Sano K, Kubota D, et al Establishment of Rapid and Accurate Screening System for Molecular Target Therapy of Osteosarcoma. Technol Cancer Res Treat. 2022;21:15330338221138217. doi: 10.1177/15330338221138217.
» https://doi.org/10.1177/15330338221138217 -
48 Niini T, Lahti L, Michelacci F, Ninomiya S, Hattinger CM, Guled M, et al. Array comparative genomic hybridization reveals frequent alterations of G1/S checkpoint genes in undifferentiated pleomorphic sarcoma of bone. Genes Chromosomes Cancer. 2011;50(5):291-306. doi: 10.1002/gcc.20851.
» https://doi.org/10.1002/gcc.20851 - 49 Cao L, Zhang S, Ma W. The Vascular Endothelial Growth Factor (VEGF) Gene Rs2010963 and Rs3025039 Polymorphisms and Risk of Osteosarcoma in Chinese Population: Evidence from a Case-Control Study and a Meta-Analysis. Int J Clin Exp Pathol. 2016;9(11):11276-11288.
-
50 Wang Z, Wen P, Luo X, Fang X, Wang Q, Ma F, et al. Association of the vascular endothelial growth factor (VEGF) gene single-nucleotide polymorphisms with osteosarcoma susceptibility in a Chinese population. Tumour Biol. 2014;35(4):3605-10. doi: 10.1007/s13277-013-1475-7.
» https://doi.org/10.1007/s13277-013-1475-7 - 51 Ferguson JL, Turner SP. Bone Cancer: Diagnosis and Treatment Principles. Am Fam Physician. 2018;98(4):205-213.
-
52 Fortuna-Costa A, Granato RA, Meohas W, Lopes ACS, Caruso AC, Castro E Silva Pinheiro R, et al. An association between successful engraftment of osteosarcoma patient-derived xenografts and clinicopathological findings. Histol Histopathol. 2020;35(11):1295-1307. doi: 10.14670/HH-18-256.
» https://doi.org/10.14670/HH-18-256 -
53 Kicinski M, Vangronsveld J, Nawrot TS. An epidemiological reappraisal of the familial aggregation of prostate cancer: a meta-analysis. PLoS One. 2011;6(10):e27130. doi: 10.1371/journal.pone.0027130.
» https://doi.org/10.1371/journal.pone.0027130 -
54 Phi XA, Houssami N, Hooning MJ, Riedl CC, Leach MO, Sardanelli F, et al. Accuracy of screening women at familial risk of breast cancer without a known gene mutation: Individual patient data meta-analysis. Eur J Cancer. 2017;85:31-38. doi: 10.1016/j.ejca.2017.07.055.
» https://doi.org/10.1016/j.ejca.2017.07.055 -
55 Ji J, Hemminki K. Familial risk for histology-specific bone cancers: an updated study in Sweden. Eur J Cancer. 2006;42(14):2343-9. doi: 10.1016/j.ejca.2005.11.043.
» https://doi.org/10.1016/j.ejca.2005.11.043 -
56 Miller BJ, Cram P, Lynch CF, Buckwalter JA. Risk factors for metastatic disease at presentation with osteosarcoma: an analysis of the SEER database. J Bone Joint Surg Am. 2013;95(13):e89. doi: 10.2106/JBJS.L.01189.
» https://doi.org/10.2106/JBJS.L.01189 -
57 Zhang C, Guo X, Xu Y, Han X, Cai J, Wang X, et al. Lung metastases at the initial diagnosis of high-grade osteosarcoma: prevalence, risk factors and prognostic factors. A large population-based cohort study. Sao Paulo Med J. 2019;137(5):423-429. doi: 10.1590/1516-3180.2018.0381120619.
» https://doi.org/10.1590/1516-3180.2018.0381120619 -
58 Tabone MD, Brugières L, Piperno-Neumann S, Selva MA, Marec-Bérard P, Pacquement H, et al. Prognostic impact of blood and urinary angiogenic factor levels at diagnosis and during treatment in patients with osteosarcoma: a prospective study. BMC Cancer. 2017;17(1):419. doi: 10.1186/s12885-017-3409-z.
» https://doi.org/10.1186/s12885-017-3409-z -
59 Shaikh AB, Li F, Li M, He B, He X, Chen G, et al. Present Advances and Future Perspectives of Molecular Targeted Therapy for Osteosarcoma. Int J Mol Sci. 2016;17(4):506. doi: 10.3390/ijms17040506.
» https://doi.org/10.3390/ijms17040506 -
60 Topkas E, Cai N, Cumming A, Hazar-Rethinam M, Gannon OM, Burgess M, et al. Auranofin is a potent suppressor of osteosarcoma metastasis. Oncotarget. 2016;7(1):831-44. doi: 10.18632/oncotarget.5704.
» https://doi.org/10.18632/oncotarget.5704 - 61 Peng N, Gao S, Guo X, Wang G, Cheng C, Li M, et al. Silencing of VEGF inhibits human osteosarcoma angiogenesis and promotes cell apoptosis via VEGF/PI3K/AKT signaling pathway. Am J Transl Res. 2016;8(2):1005-15.
-
62 Zheng B, Zhou C, Qu G, Ren C, Yan P, Guo W, et al. VEGFR2 Promotes Metastasis and PD-L2 Expression of Human Osteosarcoma Cells by Activating the STAT3 and RhoA-ROCK-LIMK2 Pathways. Front Oncol. 2020;10:543562. doi: 10.3389/fonc.2020.543562.
» https://doi.org/10.3389/fonc.2020.543562 -
63 Hassanain O, Alaa M, Khalifa MK, Kamal N, Albagoury A, El Ghoneimy AM. Genetic variants associated with osteosarcoma risk: a systematic review and meta-analysis. Sci Rep. 2024;14(1):3828. doi: 10.1038/s41598-024-53802-w.
» https://doi.org/10.1038/s41598-024-53802-w -
64 Brogan IJ, Khan N, Isaac K, Hutchinson JA, Pravica V, Hutchinson IV. Novel polymorphisms in the promoter and 5' UTR regions of the human vascular endothelial growth factor gene. Hum Immunol. 1999;60(12):1245-9. doi: 10.1016/s0198-8859(99)00132-9.
» https://doi.org/10.1016/s0198-8859(99)00132-9 -
65 Maney DL. Polymorphisms in sex steroid receptors: From gene sequence to behavior. Front Neuroendocrinol. 2017;47:47-65. doi: 10.1016/j.yfrne.2017.07.003.
» https://doi.org/10.1016/j.yfrne.2017.07.003 -
66 Zondervan KT, Cardon LR, Kennedy SH. What makes a good case-control study? Design issues for complex traits such as endometriosis. Hum Reprod. 2002;17(6):1415-23. doi: 10.1093/humrep/17.6.1415.
» https://doi.org/10.1093/humrep/17.6.1415 -
67 Cardoso JV, Abrão MS, Vianna-Jorge R, Ferrari R, Berardo PT, Machado DE, et al. Combined effect of vascular endothelial growth factor and its receptor polymorphisms in endometriosis: a case-control study. Eur J Obstet Gynecol Reprod Biol. 2017;209:25-33. doi: 10.1016/j.ejogrb.2016.10.046.
» https://doi.org/10.1016/j.ejogrb.2016.10.046 -
68 Pena SD, Di Pietro G, Fuchshuber-Moraes M, Genro JP, Hutz MH, Kehdy Fde S, et al. The genomic ancestry of individuals from different geographical regions of Brazil is more uniform than expected. PLoS One. 2011;6(2):e17063. doi: 10.1371/journal.pone.0017063.
» https://doi.org/10.1371/journal.pone.0017063 -
69 Pickrell JK, Reich D. Toward a new history and geography of human genes informed by ancient DNA. Trends Genet. 2014;30(9):377-89. doi: 10.1016/j.tig.2014.07.007.
» https://doi.org/10.1016/j.tig.2014.07.007 -
70 Ma WQ, Wang Y, Han XQ, Zhu Y, Liu NF. Association of genetic polymorphisms in vascular endothelial growth factor with susceptibility to coronary artery disease: a meta-analysis. BMC Med Genet. 2018;19(1):108. doi: 10.1186/s12881-018-0628-3.
» https://doi.org/10.1186/s12881-018-0628-3 -
71 Angioni S, D'Alterio MN, Coiana A, Anni F, Gessa S, Deiana D. Genetic Characterization of Endometriosis Patients: Review of the Literature and a Prospective Cohort Study on a Mediterranean Population. Int J Mol Sci. 2020;21(5):1765. doi: 10.3390/ijms21051765.
» https://doi.org/10.3390/ijms21051765
-
The study was conducted at Instituto Nacional de Traumatologia e Ortopedia (INTO), Av. Brasil, 500, Caju, Rio de Janeiro, RJ, Brazil. 20940-070.
Supplementary Figure 1
Publication Dates
-
Publication in this collection
13 Oct 2025 -
Date of issue
Sept 2025
History
-
Received
05 Nov 2024 -
Accepted
13 Feb 2025






