Next Article in Journal
Disparity of Hepatocellular Carcinoma in Tumor Microenvironment-Related Genes and Infiltrating Immune Cells between Asian and Non-Asian Populations
Previous Article in Journal
Uniparental Lineages from the Oldest Indigenous Population of Ecuador: The Tsachilas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Clinical Relevance of VEGFA (rs3025039) +936 C>T Polymorphism in Primary Myelofibrosis: Susceptibility, Clinical Co-Variates, and Outcomes

1
Center for the Study of Myelofibrosis, Laboratory of Biochemistry, Biotechnology and Advanced Diagnostics, Istituto di Ricovero e Cura a Carattere Scientifico Policlinico S. Matteo Foundation, 27100 Pavia, Italy
2
Laboratory of Biochemistry, Biotechnology and Advanced Diagnostics, Istituto di Ricovero e Cura a Carattere Scientifico Policlinico San Matteo Foundation, 27100 Pavia, Italy
3
Centre for Haematology Research, Department of Immunology and Inflammation, Imperial College London, London SW7 2BU, UK
*
Author to whom correspondence should be addressed.
Submission received: 25 June 2021 / Revised: 17 August 2021 / Accepted: 18 August 2021 / Published: 20 August 2021
(This article belongs to the Section Molecular Genetics and Genomics)

Abstract

:
We evaluated the association of VEGFA rs3025039 polymorphism with clinical co-variates and outcomes in 849 subjects with primary myelofibrosis (PMF) and 250 healthy controls. Minor T-allele frequency was higher in subjects with JAK2V617F compared with those without JAK2V617F (18% vs. 13%; p = 0.014). In subjects with JAK2V617F, the TT genotype was associated at diagnosis with lower platelet concentrations (p = 0.033), higher plasma LDH concentration (p = 0.005), higher blood CD34-positive cells (p = 0.027), lower plasma cholesterol concentration (p = 0.046), and higher concentration of high-sensitivity C-reactive protein (p = 0.018). These associations were not found in subjects with PMF without JAK2V617F. In subjects with the TT genotype, risk of death was higher compared with subjects with CC/CT genotypes (HR = 2.12 [1.03, 4.35], p = 0.041). Finally, the TT genotype was associated with higher frequency of deep vein thrombosis in typical sites (12.5% vs. 2.5%; OR = 5.46 [1.51, 19.7], p = 0.009). In conclusion, in subjects with PMF, the VEGFA rs3025039 CT or TT genotypes are more common in those with JAK2V617F than in those without JAK2V67F mutation and are associated with disease severity, poor prognosis, and risk of deep vein thrombosis.

1. Introduction

Primary myelofibrosis (PMF) is a hematological cancer characterized by abnormal proliferation and differentiation of hematopoietic progenitors, variable degrees of bone marrow fibrosis, cytopenias, elevated blood CD34-positive stem and progenitor hematopoietic cells, splenomegaly, and risk of transformation to a blast phase. It is driven by gain-of-function mutations in Janus kinase-2 (JAK2), calreticulin (CALR), or myeloproliferative leukemia virus (MPL). These persons sometimes have mutations in genes associated with other hematological cancers including ASXL1, EZH2, DMNT3A, IDH1 and IDH2 [1]. However, mutation topography does not completely explain the different phenotypes of PMF, and PMF is associated with different risks of disease progression or transformation.
Angiogenesis is important in the development of cancers and critical in the pathogenesis of PMF [2]. Angiogenesis reflects integrated actions of several vascular growth factors, the most potent being vascular endothelial growth factor A (VEGFA). Considerable data indicate that serum VEGFA concentration is increased in persons with PMF [2,3,4,5,6,7,8,9,10,11]. In humans, VEGFA is located on chromosome 6p21.3. Single-nucleotide variants (SNVs) of VEGFA affect gene expression by altering key regulatory sequences or by altering mRNA. rs3025039 +936 C>T polymorphism in the 3′-UTR has been studied extensively and is reportedly associated with diverse diseases, including cancers [12,13,14,15,16,17].
It is not known whether there is an association between VEGFA rs3025039 genotypes and PMF. We performed a population-based study in 849 subjects with PMF, interrogating correlations between VEGFA rs3025039 genotypes and mutation topography, clinical co-variates and outcomes.

2. Materials and Methods

2.1. Study Population

Stored DNA from blood granulocytes of 849 consecutive subjects with PMF seen at the Center for the Study of Myelofibrosis of the IRCCS Policlinico S. Matteo Foundation in Pavia and included in the institutional database was the primary source material of this study. Clinical data were collected retrospectively on the first visit and prospectively thereafter. In all subjects, the diagnosis was confirmed by a review of the initial bone marrow biopsy. The subjects were classified as PMF based on the WHO diagnostic criteria at the time of their first visit and re-classified according to 2017 revised WHO criteria [18]. After the first examination, visits were scheduled every six months. All subjects gave written informed consent approved by the IRCCS Policlinico S. Matteo Foundation Institutional Ethics Committee for their data to be collected and stored in the database of the Centre for the Study of Myelofibrosis. Controls were healthy Italian subjects belonging to the bone marrow donor transplant registry whose samples were made anonymous for the purpose of the study (N = 250).

2.2. SNV Analyses

DNA was isolated from blood granulocytes, obtained by density gradient centrifugation, using the QIAamp DNA Blood Mini Kit (QIAGENSciences Inc., Germantown, MD, USA). SNP genotyping was conducted using a predesigned, two-labeled (VIC-FAM) TaqMan Assay C_16198794_10 (Applied Biosystems, Foster City, CA, USA) and the reactions were performed on a CFX96 Real-time PCR Detection System (Biorad Company, Hercules, CA, USA) according to the manufacturer’s instructions.

2.3. Data Analyses

Data collected at diagnosis included sex, age, spleen size (cm from the costal margin and spleen index), complete blood count with differential, and serum lactate dehydrogenase and cholesterol concentrations. Blood concentrations of CD34-positive cells, blood CXCR4 expression on CD34-positive cells, and serum concentration of high-sensitivity C-reactive protein (hs-CRP) were also quantified in most subjects [19,20,21].
JAK2V617F and MPLW515 were detected by real-time polymerase chain reaction (RT-PCR) or high-resolution melting analyses. CALR mutations were identified by capillary electrophoresis and bi-directional sequencing. In selected subjects, we performed cytogenetic testing and tested for mutations in ASXL1, EZH2, DMNT3, IDH1, and IDH2. Cytogenetic analyses and reporting were performed according to the International System for Human Cytogenetic Nomenclature criteria using standard techniques [22]. Next-generation sequencing (NGS) was performed using the diagnostic panel commercially available Oncomine Myeloid Research Assay® (Thermo Fisher, Waltham, MA, USA). Genomic and transcript analyses were performed with IonReporter® software by applying the latest release of the Myeloid workflow (Thermo Fisher, Waltham, MA, USA). Variations causing missense, frameshift, an altered stop/initiation codon, in-frame insertion/deletion or variants affecting splice sites were regarded as mutations. In all subjects, bone marrow biopsies were analyzed at diagnosis and before any therapy started. Grading of bone marrow fibrosis was determined according to the EUMNET consensus [23].

2.4. Statistical Analyses

Continuous variables were presented as median with interquartile range (IQR). Pearson’s χ test with one degree of freedom was employed to compare allele and genotype frequencies between cases and controls. Deviations from the Hardy–Weinberg equilibrium within cases and controls were tested by Fisher’s exact test. Dominant models (presence vs. absence of the minor allele), and recessive models (presence vs. absence of two copies of the minor allele) were used. The independent contribution of rs3025039 SNV to odds of PMF phenotypes was assessed by logistic regression analysis. Because this is an exploratory study, we did not introduce a correction for multiple comparisons [24].
To test whether the VEGFA rs3025039 genotypes correlated with hematological co-variates, we analyzed the risk of developing a hemoglobin concentration <100 g/L, a platelet concentration <150 × 109/L, a WBC <4 × 109/L or >12 × 109/L, blood CD34-positive cells >100 × 106/L or splenomegaly >10 cm below left costal margin. The subjects were censored at the time of therapy or last follow-up visit.
We also analyzed the occurrence of major thromboses including arterial thromboses (myocardial infarction and stroke), venous thromboses (deep vein thromboses and pulmonary embolism), and thromboses in atypical sites (portal vein thrombosis, Budd–Chiari syndrome, and cerebral sinus thrombosis). Survival was calculated as the time from diagnosis until death. A log-rank test was used to test for survival differences for subjects with different VEGFA rs3025039 genotypes. The hazard ratio with 95% confidence interval (CI) for different genotypes was estimated by Cox regression analysis. Computations were performed with STATISTICA© software (Dell Technologies Inc., Round Rock, TX, USA).

3. Results

Table 1 displays demographic and clinical co-variates of the study subjects. In all, 503 subjects were male (59%) with a median age of 52 years (IQR, 46–61 years), and 484 subjects (57%) had prefibrotic myelofibrosis. The IPSS risk distribution was 61% low, 17% intermediate-1, 14% intermediate-2, and 8% high [25].

3.1. Correlation between VEGFA rs3025039 Genotypes and PMF

VEGFA rs3025039 genotypes in subjects with PMF and normals were consistent with the Hardy–Weinberg equilibrium (p = 0.86 and p = 0.54). Neither the frequencies of VEGFA rs3025039 genotypes nor those of the minor T-allele of the SNV differed significantly between the PMF subjects and local healthy controls (Table 2).

3.2. VEGFA rs3025039 Genotypes and Somatic Driver Mutations

PMF cohorts defined by the somatic driver mutations including JAK2V617F, CALR, MPL and the so-called triple negative genotypes were in Hardy–Weinberg equilibrium (p > 0.05). In subjects with JAK2V617F the frequency of the VEGFA rs3025039 minor T-allele was like that of normals and of the entire PMF cohort (Supplementary Materials Table S1). However, T-allele frequency was higher in subjects with JAK2V617F compared with those without JAK2V617F (18% vs. 13%; OR = 1.43 [1.07, 1.91], p = 0.014). In addition, in subjects with JAK2V617F, the frequency of CT and TT genotypes was higher compared with that of the CC genotype (73% vs. 63%; OR = 1.58, 95% [1.41, 2.20], p = 0.006). Adjusting for confounding by age and sex did not alter the findings. In subjects without JAK2V617F, there was a lower frequency of the T-allele genotype compared with the normals (13% vs. 18%; OR = 0.68 [0.48, 0.94], p = 0.02).

3.3. VEGFA rs3025039 Genotypes and Clinical Co-Variates at Diagnosis

Sex, age, hemoglobin concentration, and WBC distributions were not significantly correlated with VEGFA rs3025039 genotypes (Table S2). Consequently, no association was detected between the polymorphism genotypes and IPSS prognostic score. However, the subjects with a TT genotype had a lower platelet concentration (median, 281 vs. 474 × 109/L; p = 0.02) and a higher plasma LDH concentration (1.58 vs. 1.26 times ULN; p = 0.04). In total, 33% of the subjects with the TT genotype had platelet concentration <150 × 109/L, whereas in the subjects with the CC/CT genotypes, this proportion was 14% (OR = 3.10 [1.29, 7.42], p = 0.01). The subjects with a TT genotype had an increased plasma LDH concentration compared with 66% in the subjects with the CC/CT genotypes (p = 0.076). The subjects with the CT/TT genotypes had more frequent bone marrow fibrosis grade 2–3 compared with the subjects with the CC genotype (48% vs. 41; OR, 1.34 [1, 1.81], p = 0.05).
We tested whether these associations were influenced by driver mutation. In the subjects with JAK2V617F these associations persisted (Table S3). In addition, the subjects with JAK2V617F and the TT genotype had higher blood concentration of CD34-positive cells compared with other genotypes (median, 28 vs. 7 × 106; p = 0.03). Moreover, the subjects with CT/TT genotypes had lower plasma cholesterol concentration (median, 146 vs. 160 mg/dL; p = 0.05) and higher plasma hs-CRP concentration compared with the subjects with the CC genotype (median, 0.31 vs. 0.12; p = 0.02). None of these associations were detected in the subjects with driver mutations other than JAK2V617F even when the sample size of the comparisons provides sufficient statistical power (Table S4).

3.4. VEGFA rs3025039 Genotypes and Outcomes

After a median follow-up of 77 months (IQR, 36–150 months), 364 subjects (43%) developed a hemoglobin concentration <100 g/L, 42% splenomegaly >10 cm below the left costal margin, 41%, a WBC > 12 × 109/L, 30% platelet concentration <150 × 109/L, 15%, a WBC < 4 × 109/L, and 34%, blood CD34-positive cells >100 × 106/L. A total of 102 subjects (12%) had a splenectomy, 93 (11%) received a hematopoietic cell transplant, 155 (18%) had blast transformation, and 226 (27%) died. Median survival was 21 years (IQR, 10, 32 years).
Hazards of developing a hemoglobin concentration <100 g/L, splenomegaly >10 cm from the costal margin, WBC > 12 × 109/L, or WBC < 4 × 109/L, did not differ between VEGFA rs3025039 genotypes. However, T-allele frequency stratified the subjects into two well-delineated risk groups for the development of thrombocytopenia (recessive model, HR = 2.17 [1.25, 3.85], p = 0.0064; Table 3). In total, 11 of 24 subjects with the TT genotype had a platelet concentration <150 × 109/L at diagnosis or developed it after diagnosis compared with 238 of 821 subjects with CC/CT genotypes. Association of the TT genotype with risk of thrombocytopenia was significant only in the subjects with JAK2V617F (HR = 2.5 [1.35, 5.0], p = 0.004; Table S5). In addition, in these subjects, the CT/TT genotype was associated with the risk of developing a spleen >10 cm below the costal margin (HR = 1.31 [1.01, 1.72], p = 0.05).
Risk of death was higher in the subjects with the TT genotype compared to the subjects with the CC or CT genotype (OR = 1.92; [1.06, 3.45], p = 0.03; Table 3; Figure 1). However, this increased risk was evident only in the subjects with JAK2V617F (Table S5). The few events in the subjects without JAK2V617F precluded critical testing.

3.5. VEGFA rs3025039 Genotypes and Thromboses

In total, 170 subjects (20%) had a major thrombotic event (Table 4), 97 (57%) of which were vein thrombosis in atypical sites, 49 (29%) arterial thrombosis, and 24 (14%) vein thrombosis in typical sites. Incidences did not correlate with VEGFA rs3025039 genotype. However, a higher frequency of deep vein thrombosis in typical sites was detected in the subjects with the TT genotype compared with other genotypes (12.5% vs. 2.5%; OR = 5.46 [1.51, 19.7], p = 0.0096).

4. Discussion

We report correlations between VEGFA rs3025039 genotypes with disease frequency, clinical co-variates at diagnosis, and outcomes in subjects with PMF. According to the WHO classification, the PMF cohort includes prefibrotic myelofibrosis and overt myelofibrosis. We show that PMF subjects with the rs3025039 minor T-allele have an increased susceptibility to the JAK2V617F driver mutation, in keeping with data indicating that other genetic polymorphisms in JAK2, MECOM, TERT, TET2, HBS1L-MYB, and the corticosteroid receptor predispose to acquiring JAK2V617F [26,27,28,29,30,31,32,33].
In the subjects with JAK2V617F, the TT genotype was associated with more severe disease at diagnosis, i.e., lower platelet concentration, higher plasma LDH concentration, splenomegaly, increased blood concentration of CD34-positive cells, lower serum cholesterol concentration, more inflammation evidenced by higher plasma hs-CRP concentration, and an increased risk of developing thrombocytopenia. These factors translated to worse survival in subjects with JAK2V617F.
These data suggest a functional role for VEGFA rs3025039 SNV in PMF severity and progression. They also suggest that this influence operates predominantly in persons with JAK2V617F but not in other driver mutations. However, this observation needs confirmation because of the relatively few subjects with MPL mutations and those who were triple negative.
How VEGFA expression affects PMF disease activity and progression is unknown. The minor T-allele of the VEGFA rs3025039 in the 3′-UTR is associated with lower plasma VEGFA concentrations in normals and disease [34,35,36]. Other data suggest a discordance between VEGFA expression and VEGFA plasma and that the other genes involved in VEGFA signaling are also important in cancer susceptibility and phenotype [37,38,39].
Pathways active by JAK2V617F and CALR mutation in PMF are similar but not identical. For example, CALR mutations interact with the thrombopoietin receptor (MPL), a protein involved in calcium regulation and processing of IL-6 receptors, interactions not found in JAK2V617F [40,41]. In addition, different mechanisms of evasion of immune surveillance are reported [42,43,44,45]. We recently reported that myeloid expansion, blast transformation, and survival correlate with sIL-2Rα-mediated immune escape in subjects with JAK2V617F but not those with a CALR mutation [46].
The subjects with the VEGFA rs3025039 TT genotype had a higher risk of deep vein thrombosis in typical sites. Current knowledge about the potential impact of the rs3025039 SNV in VEGFA on angiogenesis and endothelial function is still very limited; thus, the explanation for the observed relationship with thrombosis remains highly speculative.
Our study has limitations. It is retrospective and predominantly descriptive, limiting the certainty of some of our conclusions. Moreover, the numbers of subjects with MPL mutations or those who were triple negative were small, limiting our power to detect some effects. Further studies are needed to explore the effects of VEGFA genotypes and gene–gene interactions in PMF and to validate our findings [47].

5. Conclusions

In conclusion, we found that persons with PMF and VEGFA rs3025039 minor T-allele genotypes are more likely to have a JAK2V617F compared to other driver mutations, and different clinical co-variates at diagnosis and clinical outcomes. These findings help to enrich our knowledge about the genetic basis of PMF and, if validated, they have relevant clinical and prognostic implications.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/genes12081271/s1, Table S1: Frequency of the genotypes and of the T-allele of the VEGFA rs3025039 polymorphism in subjects with primary myelofibrosis (PMF) stratified according to the somatic driver mutations. Table S2: Clinical and laboratory co-variates by VEGFA rs3025039 genotype in subjects with primary myelofibrosis (PMF). Table S3: Clinical and laboratory co-variates of PMF subjects with JAK2V617F and VEGFA rs3025039 genotype Table S4: Clinical and laboratory co-variates of PMF subjects without JAK2V617F and VEGFA rs3025039 genotype. Table S5: Hazard ratio (HR) of the outcomes of PMF subjects with JAK2V617F and VEGFA rs3025039 genotypes.

Author Contributions

L.V. and G.B. designed the project; V.R., M.M., R.C., A.C., and C.A. led the database sample collection and clinical characterization efforts; P.C. and L.V. did genotyping for the dataset; G.B. and R.P.G. prepared the typescript. All authors reviewed and approved the typescript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board statement

The research was conducted in accordance with the World Medical Association Declaration of Helsinki. All subjects gave written informed consent that their data could be collected and stored in the database of the Center for the Study of Myelofibrosis, and the informed consent form was approved by the IRCCS Policlinico S. Matteo Foundation Institutional Ethics Committee. The Ethics Committee of the Hospital also approved a written informed consent for patients to donate samples for molecular research (reference number 20110004143 of 26 September 2011).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Supported by AIRC 5 × 1000 called “Metastatic disease: the key unmet need in oncology” to the MYNERVA project, #21267 (MYeloid Research Venture AIRC), and by Ricerca Corrente IRCCS Policlinico San Matteo Foundation, Pavia, Italy, project number 874, code number 08054517, received by Vittorio Rosti (www.sanmatteo.org).

Conflicts of Interest

R.P.G. is a consultant to BeiGene Ltd., Fusion Pharma LLC, LaJolla NanoMedical Inc., Mingsight Pharmaceuticals Inc., Kite Pharma, and CStone Pharmaceuticals; advisor to Antegene Biotech LLC; medical director at FFF Enterprises Inc.; partner in AZACA Inc.; on the Board of Directors, RakFond Foundation for Cancer Research Support; and on the Scientific Advisory Board, StemRad Ltd.

References

  1. Grinfeld, J.; Nangalia, J.; Baxter, E.J.; Wedge, D.C.; Angelopoulos, N.; Cantrill, R.; Godfrey, A.L.; Papaemmanuil, E.; Gundem, G.; MacLean, C.; et al. Classification and Personalized Prognosis in Myeloproliferative Neoplasms. N. Engl. J. Med. 2018, 379, 1416–1430. [Google Scholar] [CrossRef] [PubMed]
  2. Gadomska, G.; Stankowska, K.; Boinska, J.; Ślusarz, R.; Tylicka, M.; Michalska, M.; Jachalska, A.; Rość, D. VEGF-A, sVEGFR-1, and sVEGFR-2 in BCR-ABL negative myeloproliferative neoplasms. Medicina 2017, 53, 34–39. [Google Scholar] [CrossRef] [PubMed]
  3. Medinger, M.; Passweg, J. Angiogenesis in myeloproliferative neoplasms, new markers and future directions. Memo-Mag. Eur. Med. Oncol. 2014, 7, 206–210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Boiocchi, L.; Vener, C.; Savi, F.; Bonoldi, E.; Moro, A.; Fracchiolla, N.S.; Iurlo, A.; Deliliers, G.L.; Coggi, G.; Bosari, S.; et al. Increased expression of vascular endothelial growth factor receptor 1 correlates with VEGF and microvessel density in Philadelphia chromosome-negative myeloproliferative neoplasms. J. Clin. Pathol. 2011, 64, 226–231. [Google Scholar] [CrossRef]
  5. Medinger, M.; Skoda, R.; Gratwohl, A.; Theocharides, A.; Buser, A.; Heim, D.; Dirnhofer, S.; Tichelli, A.; Tzankov, A. Angiogenesis and vascular endothelial growth factor-/receptor expression in myeloproliferative neoplasms: Correlation with clinical parameters and JAK2-V617F mutational status. Br. J. Haematol. 2009, 146, 150–157. [Google Scholar] [CrossRef]
  6. Alonci, A.; Allegra, A.; Bellomo, G.; Penna, G.; D’Angelo, A.; Quartarone, E.; Musolino, C. Evaluation of circulating endothelial cells, VEGF and VEGFR2 serum levels in patients with chronic myeloproliferative diseases. Hematol. Oncol. 2008, 26, 235–239. [Google Scholar] [CrossRef]
  7. Gianelli, U.; Vener, C.; Raviele, P.R.; Savi, F.; Somalvico, F.; Calori, R.; Iurlo, A.; Radaelli, F.; Fermo, E.; Bucciarelli, P.; et al. VEGF expression correlates with microvessel density in Philadelphia chromosome-negative chronic myeloproliferative disorders. Am. J. Clin. Pathol. 2007, 128, 966–973. [Google Scholar] [CrossRef] [Green Version]
  8. Panteli, K.; Bai, M.; Hatzimichael, E.; Zagorianakou, N.; Agnantis, N.J.; Bourantas, K. Serum levels, and bone marrow immunohistochemical expression of, vascular endothelial growth factor in patients with chronic myeloproliferative diseases. Hematology 2007, 12, 481–486. [Google Scholar] [CrossRef]
  9. Steurer, M.; Zoller, H.; Augustin, F.; Fong, D.; Heiss, S.; Strasser-Weippl, K.; Gastl, G.; Tzankov, A. Increased angiogenesis in chronic idiopathic myelofibrosis: Vascular endothelial growth factor as a prominent angiogenic factor. Hum. Pathol. 2007, 38, 1057–1064. [Google Scholar] [CrossRef]
  10. Ho, C.L.; Arora, B.; Hoyer, J.D.; Wellik, L.E.; Mesa, R.A.; Tefferi, A. Bone marrow expression of vascular endothelial growth factor in myelofibrosis with myeloid metaplasia. Eur. J. Haematol. 2005, 74, 35–39. [Google Scholar] [CrossRef]
  11. Di Raimondo, F.; Azzaro, M.P.; Palumbo, G.A.; Bagnato, S.; Stagno, F.; Giustolisi, G.M.; Cacciola, E.; Sortino, G.; Guglielmo, P.; Giustolisi, R. Elevated vascular endothelial growth factor (VEGF) serum levels in idiopathic myelofibrosis. Leukemia 2001, 15, 976–980. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Grimm, C.; Watrowski, R.; Polterauer, S.; Baumühlner, K.; Natter, C.; Rahhal, J.; Heinze, G.; Schuster, E.; Hefler, L.; Reinthaller, A. Vascular endothelial growth factor gene polymorphisms and risk of cervical intraepithelial neoplasia. Int. J. Gynecol. Cancer 2011, 21, 597–601. [Google Scholar] [CrossRef] [PubMed]
  13. Hou, Q.; Li, M.-Y.; Huang, W.-T.; Wei, F.-F.; Peng, J.-P.; Lou, M.; Qiu, J.-G. Association between three VEGF polymorphisms and renal cell carcinoma susceptibility: A meta-analysis. Oncotarget 2017, 8, 50061–50070. [Google Scholar] [CrossRef]
  14. Kim, D.H.; Lee, N.Y.; Lee, M.-H.; Sohn, S.K.; Do, Y.R.; Park, J.Y. Vascular endothelial growth factor (VEGF) gene (VEGFA) polymorphism can predict the prognosis in acute myeloid leukaemia patients. Br. J. Haematol. 2007, 140, 71–79. [Google Scholar] [CrossRef] [Green Version]
  15. Mandal, R.K.; Yadav, S.S.; Panda, A.K.; Khattri, S. Vascular endothelial growth factor 936 c>T polymorphism increased oral cancer risk: Evidence from a meta-analysis. Genet. Test. Mol. Biomark. 2013, 17, 543–547. [Google Scholar] [CrossRef]
  16. Rodrigues, P.; Furriol, J.; Tormo, E.; Ballester, S.; Lluch, A.; Eroles, P. The single-nucleotide polymorphisms +936 C/T VEGF and 710 C/T VEGFR1 are associated with breast cancer protection in a Spanish population. Breast Cancer Res. Treat. 2012, 133, 769–778. [Google Scholar] [CrossRef]
  17. Zhang, J.; Yang, J.; Chen, Y.; Mao, Q.; Li, S.; Xiong, W.; Lin, Y.; Chen, J.; Ge, J. Genetic Variants of VEGF (rs201963 and rs3025039) and KDR (rs7667298, rs2305948, and rs1870377) Are Associated with Glioma Risk in a Han Chinese Population: A Case-Control Study. Mol. Neurobiol. 2015, 53, 2610–2618. [Google Scholar] [CrossRef]
  18. Swerdow, S.H.; Campo, E.; Harris, N.L.; Jaffe, E.S.; Pileri, S.A.; Stein, H.; Thiele, J. WHO Classification of Tumours of Haematopoietic and Lymphoid Tissues (Revised 4th Edition); IARC: Lyon, France, 2017. [Google Scholar]
  19. Barosi, G.; Viarengo, G.; Pecci, A.; Rosti, V.; Piaggio, G.; Marchetti, M.; Frassoni, F. Diagnostic and clinical relevance of the number of circulating CD34(+) cells in myelofibrosis with myeloid metaplasia. Blood J. Am. Soc. Hematol. 2001, 98, 3249–3255. [Google Scholar] [CrossRef]
  20. Barosi, G.; Rosti, V.; Catarsi, P.; Villani, L.; Abbà, C.; Carolei, A.; Magrini, U.; Gale, R.P.; Massa, M.; Campanelli, R. Reduced CXCR4-expression on CD34-positive blood cells predicts outcomes of persons with primary myelofibrosis. Leukemia 2021, 35, 468–475. [Google Scholar] [CrossRef] [PubMed]
  21. Barosi, G.; Massa, M.; Campanelli, R.; Fois, G.; Catarsi, P.; Viarengo, G.; Villani, L.; Poletto, V.; Bosoni, T.; Magrini, U.; et al. Primary myelofibrosis: Older age and high JAK2V617F allele burden are associated with elevated plasma high-sensitivity C-reactive protein levels and a phenotype of progressive disease. Leuk. Res. 2017, 60, 8–23. [Google Scholar] [CrossRef]
  22. Larson, D.P.; Akkari, Y.M.; Van Dyke, D.L.; Raca, G.; Gardner, J.A.; Rehder, C.W.; Kaiser-Rogers, K.A.; Eagle, P.; Yuhas, J.A.; Gu, J.; et al. Conventional Cytogenetic Analysis of Hematologic Neoplasms: A 20-Year Review of Proficiency Test Results From the College of American Pathologists/American College of Medical Genetics and Genomics Cytogenetics Committee. Arch. Pathol. Lab. Med. 2021, 145, 176–190. [Google Scholar] [CrossRef]
  23. Thiele, J.; Kvasnicka, H.M. Myelofibrosis—What’s in a name? Consensus on definition and EUMNET grading. Pathobiology 2007, 74, 89–96. [Google Scholar] [CrossRef] [PubMed]
  24. Streiner, D.L.; Norman, G.R. Correction for multiple testing: Is there a resolution? Chest 2011, 140, 16–18. [Google Scholar] [CrossRef] [PubMed]
  25. Cervantes, F.; Dupriez, B.; Pereira, A.; Passamonti, F.; Reilly, J.T.; Morra, E.; Vannucchi, A.M.; Mesa, R.A.; Demory, J.L.; Barosi, G.; et al. New prognostic scoring system for primary myelofibrosis based on a study of the International Working Group for Myelofibrosis Research and Treatment. Blood J. Am. Soc. Hematol. 2009, 113, 2895–2901. [Google Scholar] [CrossRef] [PubMed]
  26. Jones, A.V.; Chase, A.; Silver, R.T.; Oscier, D.; Zoi, K.; Wang, Y.L.; Cario, H.; Pahl, H.L.; Collins, A.; Reiter, A.; et al. JAK2 haplotype is a major risk factor for the development of myeloproliferative neoplasms. Nat. Genet. 2009, 41, 446–449. [Google Scholar] [CrossRef] [Green Version]
  27. Olcaydu, D.; Harutyunyan, A.; Jäger, R.; Berg, T.; Gisslinger, B.; Pabinger, I.; Gisslinger, H.; Kralovics, R. A common JAK2 haplotype confers susceptibility to myeloproliferative neoplasms. Nat. Genet. 2009, 41, 450–454. [Google Scholar] [CrossRef] [PubMed]
  28. Kilpivaara, O.; Mukherjee, S.; Schram, A.M.; Wadleigh, M.; Mullally, A.; Ebert, B.L.; Bass, A.; Marubayashi, S.; Heguy, A.; Garcia-Manero, G.; et al. A germline JAK2 SNP is associated with predisposition to the development of JAK2(V617F)-positive myeloproliferative neoplasms. Nat. Genet. 2009, 41, 455–459. [Google Scholar] [CrossRef] [Green Version]
  29. Tapper, W.; Jones, A.V.; Kralovics, R.; Harutyunyan, A.S.; Zoi, K.; Leung, W.; Godfrey, A.L.; Guglielmelli, P.; Callaway, A.; Ward, D.; et al. Genetic variation at MECOM, TERT, JAK2 and HBS1L-MYB predisposes to myeloproliferative neoplasms. Nat. Commun. 2015, 6, 6691. [Google Scholar] [CrossRef] [Green Version]
  30. Trifa, A.P.; Bănescu, C.; Bojan, A.S.; Voina, C.M.; Popa, Ș.; Vișan, S.; Ciubean, A.D.; Tripon, F.; Dima, D.; Popov, V.M.; et al. MECOM, HBS1L-MYB, THRB-RARB, JAK2, and TERT polymorphisms defining the genetic predisposition to myeloproliferative neoplasms: A study on 939 patients. Am. J. Hematol. 2018, 93, 100–106. [Google Scholar] [CrossRef] [Green Version]
  31. Lighezan, D.L.; Bojan, A.S.; Iancu, M.; Pop, R.M.; Gligor-Popa, Ș.; Tripon, F.; Cosma, A.S.; Tomuleasa, C.; Dima, D.; Zdrenghea, M.; et al. TET2 rs1548483 SNP Associating with Susceptibility to Molecularly Annotated Polycythemia Vera and Primary Myelofibrosis. J. Pers. Med. 2020, 10, 259. [Google Scholar] [CrossRef]
  32. Tefferi, A.; Lasho, T.L.; Mudireddy, M.; Finke, C.M.; Hanson, C.A.; Ketterling, R.P.; Gangat, N.; Pardanani, A. The germline JAK2 GGCC (46/1) haplotype and survival among 414 molecularly-annotated patients with primary myelofibrosis. Am. J. Hematol. 2019, 94, 299–305. [Google Scholar] [CrossRef] [PubMed]
  33. Poletto, V.; Rosti, V.; Villani, L.; Catarsi, P.; Carolei, A.; Campanelli, R.; Massa, M.; Martinetti, M.; Viarengo, G.; Malovini, A.; et al. A3669G polymorphism of glucocorticoid receptor is a susceptibility allele for primary myelofibrosis and contributes to phenotypic diversity and blast transformation. Blood J. Am. Soc. Hematol. 2012, 120, 3112–3117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Dong, P.P. Association of vascular endothelial growth factor expression and polymorphisms with the risk of gestational diabetes mellitus. J. Clin. Lab. Anal. 2019, 33, e22686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Bhaskari, J.; Premalata, C.S.; Shilpa, V.; Rahul, B.; Pallavi, V.R.; Ramesh, G.; Krishnamoorthy, L. Vascular endothelial growth factor polymorphisms and a synchronized examination of plasma and tissue expression in epithelial ovarian cancers. Tumour. Biol. 2016, 37, 1017–1023. [Google Scholar] [CrossRef]
  36. Breunis, W.B.; Biezeveld, M.H.; Geissler, J.; Ottenkamp, J.; Kuipers, I.M.; Lam, J.; Hutchinson, A.; Welch, R.; Chanock, S.J.; Kuijpers, T.W. Vascular endothelial growth factor gene haplotypes in Kawasaki disease. Arthritis Rheum. Off. J. Am. Coll. Rheumatol. 2006, 54, 1588–1594. [Google Scholar] [CrossRef] [PubMed]
  37. Pasqualetti, G.; Danesi, R.; Del Tacca, M.; Bocci, G. Vascular endothelial growth factor pharmacogenetics: A new perspective for anti-angiogenic therapy. Pharmacogenomics 2007, 8, 49–66. [Google Scholar] [CrossRef]
  38. Sellami, N.; Lamine, L.B.; Turki, A.; Sarray, S.; Jailani, M.; Al-Ansari, A.K.; Ghorbel, M.; Mahjoub, T.; Almawi, W.Y. Association of VEGFA variants with altered VEGF secretion and type 2 diabetes: A case-control study. Cytokine 2018, 106, 29–34. [Google Scholar] [CrossRef]
  39. Al-Habboubi, H.H.; Sater, M.S.; Almawi, A.W.; Al-Khateeb, G.M.; Almawi, W.Y. Contribution of VEGF polymorphisms to variation in VEGF serum levels in a healthy population. Eur. Cytokine Netw. 2011, 22, 154–158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Mullaly, A. Both sides now: Losses and gains of mutant CALR. Blood J. Am. Soc. Hematol. 2020, 135, 82–83. [Google Scholar] [CrossRef]
  41. Balliu, M.; Calabresi, L.; Bartalucci, N.; Romagnoli, S.; Maggi, L.; Manfredini, R.; Lulli, M.; Guglielmelli, P.; Vannucchi, A.M. Activated IL-6 signaling contributes to the pathogenesis of, and is a novel therapeutic target for, CALR-mutated MPNs. Blood Adv. 2021, 5, 2184–2195. [Google Scholar] [CrossRef] [PubMed]
  42. Leone, P.; Shin, E.C.; Perosa, F.; Vacca, A.; Dammacco, F.; Racanelli, V. MHC class I antigen processing and presenting machinery: Organization, function, and defects in tumor cells. J. Natl. Cancer Inst. 2013, 105, 1172–1187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Beatty, G.L.; Gladney, W.L. Immune escape mechanisms as a guide for cancer immunotherapy. Clin. Cancer Res. 2015, 21, 687–692. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Barkal, A.A.; Weiskopf, K.; Kao, K.S.; Gordon, S.R.; Rosental, B.; Yiu, Y.Y.; George, B.M.; Markovic, M.; Ring, N.G.; Tsai, J.M.; et al. Engagement of MHC class I by the inhibitory receptor LILRB1 suppresses macrophages and is a target of cancer immunotherapy. Nat. Immunol. 2018, 19, 76–84. [Google Scholar] [CrossRef] [PubMed]
  45. Hofmann, S.; Babiak, A.; Greiner, J. Immunotherapy for myeloproliferative neoplasms (MPN). Curr. Cancer Drug Targets 2011, 11, 72–84. [Google Scholar] [CrossRef] [PubMed]
  46. Barosi, G.; Campanelli, R.; Catarsi, P.; De Amici, M.; Abbà, C.; Viarengo, G.; Villani, L.; Gale, R.P.; Rosti, V.; Massa, M. Plasma sIL-2Rα levels are associated with disease progression in myelofibrosis with JAK2V617F but not CALR mutation. Leuk. Res. 2020, 90, 106319. [Google Scholar] [CrossRef] [PubMed]
  47. Nishi, A.; Milner, D.A., Jr.; Giovannucci, E.L.; Nishihara, R.; Tan, A.S.; Kawachi, I.; Ogino, S. Integration of molecular pathology, epidemiology and social science for global precision medicine. Expert Rev. Mol. Diagn. 2016, 16, 11–23. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (A) Probability of survival in subjects with PMF stratified for the VEGFA rs3025039 polymorphism genotypes. (B) Probability of survival in subjects with PMF bearing the JAK2V617F mutation, stratified for the VEGFA rs3025039 polymorphism genotypes. The subjects with the TT genotype had a shorter survival than those with the CC or CT genotype (p = 0.029 and p = 0.036, respectively).
Figure 1. (A) Probability of survival in subjects with PMF stratified for the VEGFA rs3025039 polymorphism genotypes. (B) Probability of survival in subjects with PMF bearing the JAK2V617F mutation, stratified for the VEGFA rs3025039 polymorphism genotypes. The subjects with the TT genotype had a shorter survival than those with the CC or CT genotype (p = 0.029 and p = 0.036, respectively).
Genes 12 01271 g001
Table 1. Baseline co-variates (N = 849).
Table 1. Baseline co-variates (N = 849).
Demographic Co-Variates
Age, y, median (IQR)52 (46–61)
Male, N (%)503 (59)
Laboratory Co-Variates
Hemoglobin, g/L, median (IQR)131 (109–148)
WBC × 109/L, median (IQR)8.5 (6.4–11.6)
Platelets × 109/L, median (IQR)467 (246–711)
Monocytes × 109/L, median (IQR)496 (332–688)
Spleen size, cm2, median (IQR) a120 (90–160)
Plasma LDH, × upper limit of normal (ULN), median (IQR) b1.29 (0.92–1.95)
Serum cholesterol, mg/dL, median (IQR) c158 (130–183)
Plasma high-sensitivity C-reactive protein, mg/dL, median (IQR) d10 (4–44)
Blood CD34-positive cells, × 106/L, median (IQR) e10 (4–43)
CXCR4/CD34, %, median (IQR) f41(21–63)
Molecular Co-Variates
JAK2V617F-positive, N (%)544 (66)
JAK2V617F allele frequency, median (IQR)41 (22–68)
CALR mutation, N (%)171 (21)
MPL mutation, N (%)44 (5)
Triple negative, N (%)68 (8)
NGS detected mutations, N (%) g46 (20)
Cytogenetic abnormalities, N (%) h87 (30)
Bone Marrow Fibrosis Grade
0, N (%)258 (30)
1, N (%)226 (27)
2, N (%)245 (29)
3, N (%)117 (14)
a Spleen size was measured using the spleen index calculated by multiplying the length of the longitudinal axis by the transverse axis. b Plasma lactic dehydrogenase activity (LDH) was available in 474 subjects. c Serum cholesterol concentration was available in 422 subjects. d Plasma high-sensitivity C-reactive protein was available in 237 subjects. e Blood CD34-positive cell concentration was available in 402 subjects. f CXCR4/CD34 ratio was available in 298 subjects. g Next-generation sequencing (NGS) data were available in 242 subjects. h Cytogenetics was available in 291 subjects.
Table 2. VEGFA rs3025039 genotypes.
Table 2. VEGFA rs3025039 genotypes.
CCCTTTCC/CTCT/TTT-Allele Frequency
PMF, N (%)849601
(70.8)
224
(26.4)
24
(2.8)
825
(97.2)
248
(29.2)
272/1698
(16)
Controls, N (%)250165
(66)
79
(31.6)
6
(2.4)
244
(97.6)
85
(34)
91/500
(18.2)
Table 3. Hazard ratio of the outcomes of subjects with PMF stratified according to the VEGFA 3025039 genotype.
Table 3. Hazard ratio of the outcomes of subjects with PMF stratified according to the VEGFA 3025039 genotype.
CT/TT (N = 248) vs. CC (N = 601)TT (N = 24) vs. CC/CT (N = 825)
Outcome HR (95% CI) p-Value HR (95% CI) p-Value
Hemoglobin < 100 g/L1.16 (0.92, 1.45)0.201.41 (0.80, 2.44)0.23
Spleen > 10 cm below left costal margin1.22 (0.97, 1.51)0.101.01 (0.51, 1.96)0.99
WBC > 12 × 109/L1.10 (0.87, 1.37)0.451.02 (0.52, 1.98)0.94
WBC < 4 × 109/L1.07 (0.73, 1.58)0.721.26 (0.46, 3.45)0.65
Platelets < 150 × 109/L1.12 (0.93, 1.51)0.392,17 (1.25, 3.85)0.006
Blood CD34-positive cells > 100 × 106/L1.25 (0.93, 1.51)0.111.11 (0.46, 2.70)0.81
Transplant1.49 (0.98, 2.27)0.071.39 (0.44, 4.35)0.58
Blast transformation1.15 (0.80, 1.65)0.431.22 (0.49, 2.94)0.67
Death1.31 (1.00, 1.72)0.121.92 (1.06, 3.45)0.03
Table 4. Major thrombotic events stratified by the VEGFA rs3025039 genotype.
Table 4. Major thrombotic events stratified by the VEGFA rs3025039 genotype.
VEGFA rs3025039 GenotypeCT/TT vs. CC
OR (95% CI)
TT vs. CC/CT
OR (95% CI)
All Subjects (N = 848)CC
(N = 600)
CT
(N = 224)
TT
(N = 24)
CC/CT
(N = 824)
CT/TT (N = 248)
Thrombotic events, N (%)170
(20)
124
(20.7)
41
(18.3)
5
(20.8)
165
(20)
46
(18.5)
OR = 0.87
(0.60, 1.27)
p = 0.48
OR = 1.05
(0.38, 2.85)
p = 0.92
Arterial thrombosis, N (% of PMF cases)49
(5.8)
36
(6)
12
(5.3)
1
(4.2)
48
(5.8)
13
(5.2)
OR = 0.87
(0.45, 1.66)
p = 0.67
OR = 0.70
(0.09, 5.31)
p = 0.73
- In the year before diagnosis, N (% of thromboses)14
(28.6)
10
(27.7)
4
(33.3)
0
(0)
14
(29.2)
4
(30.7)
- At diagnosis, N (% of thromboses)12
(24.5)
9
(25)
2
(16.7)
1
(100)
11
(22-9)
3
(23.1)
- After diagnosis, N (% of thromboses)23
(46.9)
17
(47.2)
6
(50)
0
(0)
23
(48)
6
(46.1)
Deep vein thrombosis in typical sites, N (% of PMF cases)24
(2.8)
17
(2.8)
4
(1.8)
3
(12.5)
21
(2.5)
7
(2.8)
OR = 0.99
(0.41, 2.43)
p = 0.99
OR = 5.46
(1.51, 19.7)
p = 0.0096
- In the year before diagnosis, N (% of thromboses)3
(12.5)
1
(5.9)
1
(25)
1
(33.3)
2
(9.5)
2
(28.6)
- At diagnosis, N (% of thromboses)5
(28.8)
5
(29.4)
0
(0)
0
(0)
5
(23.8)
0
(0)
- After diagnosis, N (% of thromboses)16
(66.6)
11
(64.7)
3
(75)
2
(66.6)
14
(66.6)
5
(71.4)
Vein thrombosis in atypical sites, N (% of PMF cases)97
(11.4)
71
(11.8)
25
(11.1)
1
(4.2)
96
(11.6)
26
(9.7)
OR = 0.87
(0.54, 1.40)
p = 0.57
OR = 0.33
(0.04, 2.47)
p = 0.28
- In the year before diagnosis, N (% of thromboses)8
(8.2)
6
(8.5)
2
(8)
0
(0)
8
(8.3)
2
(7.7)
- At diagnosis, N (% of thromboses)73
(75.2)
53
(74.7)
19
(76)
1
(100)
72
(75)
20
(76.9)
- After diagnosis, N (% of thromboses)16
(16.5)
12
(16.9)
4
(16)
0
(0)
16
(16.7)
4
(15.4)
PMF = primary myelofibrosis; OR = odds ratio.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Villani, L.; Carolei, A.; Rosti, V.; Massa, M.; Campanelli, R.; Catarsi, P.; Abbà, C.; Gale, R.P.; Barosi, G. Clinical Relevance of VEGFA (rs3025039) +936 C>T Polymorphism in Primary Myelofibrosis: Susceptibility, Clinical Co-Variates, and Outcomes. Genes 2021, 12, 1271. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12081271

AMA Style

Villani L, Carolei A, Rosti V, Massa M, Campanelli R, Catarsi P, Abbà C, Gale RP, Barosi G. Clinical Relevance of VEGFA (rs3025039) +936 C>T Polymorphism in Primary Myelofibrosis: Susceptibility, Clinical Co-Variates, and Outcomes. Genes. 2021; 12(8):1271. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12081271

Chicago/Turabian Style

Villani, Laura, Adriana Carolei, Vittorio Rosti, Margherita Massa, Rita Campanelli, Paolo Catarsi, Carlotta Abbà, Robert Peter Gale, and Giovanni Barosi. 2021. "Clinical Relevance of VEGFA (rs3025039) +936 C>T Polymorphism in Primary Myelofibrosis: Susceptibility, Clinical Co-Variates, and Outcomes" Genes 12, no. 8: 1271. https://0-doi-org.brum.beds.ac.uk/10.3390/genes12081271

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop