P.O. Box 3660, Riyadh 11481, Mail Code 1515 (KAIMRC)
+966 (11) 429-4444
+966 (11) 429-4440

Publication Details

Title :

Utilizing Whole-Exome Sequencing to Characterize the Phenotypic Variability of Sickle Cell Disease


Genet Test Mol Biomarkers.

Impact Factor:



Alsultan A1, Al-Suliman AM2, Aleem A3, AlGahtani FH3, Alfadhel M4,5.


1 Department of Pediatrics, College of Medicine, King Saud University , Riyadh, Saudi Arabia .

2 Department of Medicine, King Fahad Hospital , Hofuf, Saudi Arabia .

3 Department of Internal Medicine, College of Medicine, King Saud University , Riyadh, Saudi Arabia .

4 Department of Pediatrics, King Abdullah Specialist Children’s Hospital , King Abdullah International Medical Research Centre, King Abdulaziz Medical City, Ministry of National Guard Health Affairs, Riyadh, Saudi Arabia .

5 King Saud bin Abdulaziz University for Health Sciences , Riyadh, Saudi Arabia .

Year of Publication:





Sickle cell disease (SCD) is a monogenic disease that has wide variety of phenotypes with both and environmental factors contributing to its severity.

To assess the incidence of neonatal thrombocytopenia and identify factors associated with its outcomes, namely time to disease onset, recovery duration, and platelet count.

We performed whole-exome sequencing (WES) in 22 Saudi SCD patients to identify variants that could explain differences in disease phenotypes. All variants, except those that were benign and likely benign, described in the ClinVar database, were considered in our analysis. Gene-based association testing using sequence kernel association optimal unified test (SKAT-O) with small sample adjustment was performed to evaluate the effect of multiple variants in genes on SCD phenotypes.

The mean age of participants was 28 (range, 10-48 years). All patients were homozygous for the sickle cell mutation. The Benin haplotype was present in 15 patients and the Arab-Indian haplotype in 7 patients. One patient who had both SCD and CHARGE association was heterozygous for pathogenic mutation p.Arg987Ter in the CHD7 gene. One SCD individual who had a stroke was a carrier of the pathogenic variant p.Asp36Tyr in the VKORC1 gene which is, associated with warfarin resistance. Two patients with steady hemoglobin levels of 7.5 and 7.1 g/dL were carriers of the pathogenic mutation p.Gly140Ser in the RPL5 gene that is associated with Diamond-Blackfan anemia. None of the patients were transfusion dependent. A heterozygous pathogenic mutation in the LDLR gene associated with autosomal dominant familial hypercholesterolemia was present in one patient with deep venous thrombosis, although their cholesterol level was normal. One individual with stroke was a carrier for the p.Arg284Ter variant in the NLRP12 gene, which is associated with familial cold autoinflammatory syndrome 2. Another patient with stroke and a pulmonary embolism was heterozygous for the p.Pro106Leu variant of the MPL gene, which has been associated with thrombocytosis. Coding variants in the GOLGB1, ENPP1, and PON1 genes showed no association with stroke in our study. SKAT-O analysis did not explain SCD heterogeneity.

WES provided limited information to explain the severity of SCD. Whole genome sequencing, epigenetic studies, and assessment of environmental factors might expand our knowledge of SCD heterogeneity.