Abstract
Haemolytic anaemia is variable among patients with sickle cell anaemia and can be estimated by reticulocyte count, lactate dehydrogenase, aspartate aminotransferase and bilirubin levels. Using principal component analysis of these measurements we computed a haemolytic score that we used as a subphenotype in a genome-wide association study. We identified in one cohort and replicated in two additional cohorts the association of a single nucleotide polymorphism in NPRL3 (rs7203560; chr16p13·3) (P = 6·04 × 10(-07) ). This association was validated by targeted genotyping in a fourth independent cohort. The HBA1/HBA2 regulatory elements, hypersensitive sites (HS)-33, HS-40 and HS-48 are located in introns of NPRL3. Rs7203560 was in perfect linkage disequilibrium (LD) with rs9926112 (r(2) = 1) and in strong LD with rs7197554 (r(2) = 0·75) and rs13336641 (r(2) = 0·77); the latter is located between HS-33 and HS-40 sites and next to a CTCF binding site. The minor allele for rs7203560 was associated with the -∝(3·7) thalassaemia gene deletion. When adjusting for HbF and ∝ thalassaemia, the association of NPRL3 with the haemolytic score was significant (P = 0·00375) and remained significant when examining only cases without gene deletion∝ thalassaemia (P = 0·02463). Perhaps by independently down-regulating expression of the HBA1/HBA2 genes, variants of the HBA1/HBA2 gene regulatory loci, tagged by rs7203560, reduce haemolysis in sickle cell anaemia.