Myeloproliferative neoplasms have characteristic alterations in laboratory exams, as well as genetic findings that permit their identification and differentiation. Findings involving genetic alterations in introns are not yet fully understood, but this scenario is becoming of increasing interest for understanding the etiopathogenic aspects and the role of these DNA regions in these diseases.
Essential thrombocythemia proved to be the most frequent myeloproliferative neoplasm, which are findings that align with the premises established by Torres15, who studied a population with BCR::ABL1-negative myeloproliferative neoplasms in the state of Amazonas (Brazil). Similar data were described by Macedo16, who reported a similar scenario in patients from the states of Paraná and São Paulo who had the same hematologic malignancy, and these data converge with descriptions found in other countries17, 18
The age range of individuals was between the fifth and seventh decades of life, which is consistent with what is stated in other studies19, 20. The progressive accumulation of genetic variations in hematopoietic stem cells and the biological machinery of the DNA repair system21, 22, an increase or decrease in telomeres23, 24 and cumulative exposure to risk factors throughout life, such as smoking and obesity25, 26, may explain the prevalence of this age group in the context of myeloproliferative neoplasms.
Regarding clinical characteristics, polycythemia vera (PV) showed an equal proportion of men and women, while essential thrombocythemia (ET) revealed a majority of cases involving women, and these data are in line with the literature27, 28. Some studies have demonstrated that women have an increased risk of developing myeloproliferative neoplasms29 and a higher likelihood of developing cardiovascular complications and splenomegaly26. The reason for this risk is uncertain, but changes in sex chromosomes, hormonal factors and gene expression may be possible contributors to this process28. Laboratory data, and thrombotic and hemorrhagic events presented as expected for each neoplasm: PV demonstrated a higher prevalence of increased erythrogram values and ET showed changes in the megakaryocytic series, with a higher risk of hemorrhagic events, as described by the World Health Organization3, and in other studies on the subject27, 30.
Regarding the genetic findings, PV demonstrates a higher prevalence of positive cases for the JAK2 V617F variant, since it is directly associated with the specific pathogenesis of this hematologic malignancy36 and plays a role in the constitutive activation of the JAK-STAT pathway5. It is interesting to note that 58% of our PV population was positive for the variant, which may initially differ from findings commonly described in the literature that point to JAK2 V617F frequencies of over 70% in Brazilian, Korean, Chinese, Japanese, and European patients31–35.
The limited number of PV_JAK2 V617F+ patients identified in this study is related to cytoreductive therapy. We emphasize that 66.6% of patients with PV were on cytoreductive therapy, which acts to suppress and/or decrease the variant burden of JAK2 V617F through the inhibition of the myeloproliferative process of mutated hematopoietic cells, as noted in a recent study36. This directly affects the sensitivity of the molecular detection methods used to identify the variant, underscoring the increased importance of incorporating molecular analysis for JAK2 V617F in the initial suspicions for MPNs.
In the literature, the germline haplotype 46/1, identified by the rs10974944 (C > G) variant, has a well-documented association with JAK2 V617F14, 37–39 as also observed in our study. The high frequency of the G allele of rs10974944 in individuals positive for JAK2 V617F contributes to discussions about the non-random correlation between these two genetic alterations13, 40 This relationship is in line with another finding from our study, haplotype 2 (rs10974944G/rs10815151C/rs1011004A/rs77375493T), which strengthens concepts based on the interaction between rs10974944 (C > G) and JAK2 V617F (rs77375493 - G > T). These propositions are in agreement with findings involving haplotype 46/1 in other Brazilian, Taiwanese, European, Chinese, and Japanese populations16, 32–34, 41, indicating that the possible mechanisms preceding the acquisition of JAK2 V617F are not limited to a specific ethnic group; therefore, its evolutionary basis can be considered as a genetic predisposition factor for the disease8.
Studies report a higher risk of individuals with the GG genotype of rs10974944 being positive for JAK2 V617F14, 40, 42. Consistent with the results of the aforementioned studies, our population exhibited a four-fold increase in the risk of positive JAK2 V617F in individuals with the GG genotype of rs10974944. (OR: 4.1; 95% CI: 8-13.9). These findings support the hypothesis of hypermutability, which establishes haplotype 46/1 as a dysregulating agent of the JAK2 gene, which increases the risk of DNA replication errors and conditions a mutagenic scenario for the acquisition of variants with selective advantages, such as JAK2 V617F43–45
The association of rs10974944 (G) and the JAK2 V617F VAF suggests a possible involvement of haplotype 46/1 in clonal expansion. We identified a six-fold higher risk of individuals carrying the G allele of rs10974944 and JAK2 V617F VAF of ≥ 50%. Our data indicate that the marker of haplotype 46/1 may play a role not only in the acquisition of JAK2 V617F but is also attributed to clonal expansion, maintenance, and survival. Tefferi46 suggests that JAK2 V617F is not the initial clonogenic event in MPNs but rather one of several subclones derived from an ancestral clone. This is in accordance with the notes of Pardanani et al.47, which support the hypothesis that this haplotype is located in a favorable cis regulatory environment, which facilitates the acquisition of JAK2 V617F, and which, in turn, is responsible for clonal expansion and the development of MPNs.
Furthermore, the possible role of acquired uniparental disomy, a genetic event that leads to mitotic recombination associated with neutral loss of heterozygosity of chromosome 9p in MPN patients, reducing both the haplotype and JAK2 V617F to a homozygous state14, 48, 49, cannot be ruled out. In this context, cells with both variants theoretically have a selective advantage, which conditions greater myeloproliferative potential and favors the establishment of variant cells over healthy cells, thus explaining the increased VAF in individuals with the combination rs10974944 (G) + rs77375493 (T) (JAK2 V617F) in homozygosity.
Association between the elevation of hematological indices and the presence of 46/1 is observed in the literature16, 33, 50; however, this is not a consensus among the scientific community8,53,61. Our data show significant differences in MCV, MCH values in the PV group, and RBC, Hb, and Ht in TE carriers of the G allele of rs10974944, which has been observed in previous studies7, 42, 51.
The present research is the first to analyze the 46/1 haplotype using the rs10974944 variant, present in intron 12 of JAK2, in a population from the Brazilian Amazon. The results of this study show that the rs10974944 (G) variant is associated with BCR::ABL1-negative myeloproliferative neoplasms, in patients positive for JAK2 V617F, especially those with PV, and a high allelic variant burden in these patients, and hematological alterations. Furthermore, the haplotype rs10974944G/rs10815151C/rs1011004A/rs77375493T was identified as a factor related to PV.