Impact of Genetic Variations on Thromboembolic Risk in Saudis with Sickle Cell Disease
Authors: Alshabeeb Mohammad A. , Alwadaani Deemah , Al Qahtani Farjah H. , Abohelaika Salah , Alzahrani Mohsen , Al Zayed Abdullah , Al Saeed Hussain H. , Al Ajmi Hala , Alsomaie Barrak , Rashid Mamoon , Daly Ann K. 
Published: October 2023

Background: Sickle cell disease (SCD) is a Mendelian disease characterized by multigenic phenotypes. Previous reports indicated a higher rate of thromboembolic events (TEEs) in SCD patients. A number of candidate polymorphisms in certain genes (e.g., FVL, PRT, and MTHFR) were previously reported as risk factors for TEEs in different clinical conditions. This study aimed to genotype these genes and other loci predicted to underlie TEEs in SCD patients. Methodology: A multi-center genome-wide association study (GWAS) involving Saudi SCD adult patients with a history of TEEs (n = 65) and control patients without TEE history (n = 285) was performed. Genotyping used the 10× Affymetrix Axiom array, which includes 683,030 markers. Fisher’s exact test was used to generate p-values of TEE associations with each single-nucleotide polymorphism (SNP). The haplotype analysis software tool version 1.05, designed by the University of Göttingen, Germany, was used to identify the common inherited haplotypes. Results: No association was identified between the targeted single-nucleotide polymorphism rs1801133 in MTHFR and TEEs in SCD (p = 0.79). The allele frequency of rs6025 in FVL and rs1799963 in PRT in our cohort was extremely low (<0.01); thus, both variants were excluded from the analysis as no meaningful comparison was possible. In contrast, the GWAS analysis showed novel genome-wide associations (p < 5 × 10−8) with seven signals; five of them were located on Chr 11 (rs35390334, rs331532, rs317777, rs147062602, and rs372091), one SNP on Chr 20 (rs139341092), and another on Chr 9 (rs76076035). The other 34 SNPs located on known genes were also detected at a signal threshold of p < 5 × 10−6. Seven of the identified variants are located in olfactory receptor family 51 genes (OR51B5, OR51V1, OR51A1P, and OR51E2), and five variants were related to family 52 genes (OR52A5, OR52K1, OR52K2, and OR52T1P). The previously reported association between rs5006884-A in OR51B5 and fetal hemoglobin (HbF) levels was confirmed in our study, which showed significantly lower levels of HbF (p = 0.002) and less allele frequency (p = 0.003) in the TEE cases than in the controls. The assessment of the haplotype inheritance pattern involved the top ten significant markers with no LD (rs353988334, rs317777, rs14788626882, rs49188823, rs139349992, rs76076035, rs73395847, rs1368823, rs8888834548, and rs1455957). A haplotype analysis revealed significant associations between two haplotypes (a risk, TT-AA-del-AA-ins-CT-TT-CC-CC-AA, and a reverse protective, CC-GG-ins-GG-del-TT-CC-TT-GG-GG) and TEEs in SCD (p = 0.024, OR = 6.16, CI = 1.34–28.24, and p = 0.019, OR = 0.33, CI = 0.13–0.85, respectively). Conclusions: Seven markers showed novel genome-wide associations; two of them were exonic variants (rs317777 in OLFM5P and rs147062602 in OR51B5), and less significant associations (p < 5 × 10−6) were identified for 34 other variants in known genes with TEEs in SCD. Moreover, two 10-SNP common haplotypes were determined with contradictory effects. Further replication of these findings is needed.