A genome-wide association study of total bilirubin and cholelithiasis risk in sickle cell anemia

Genetic Variation and Sickle Cell Disease Severity: A Systematic Review and Meta-Analysis

Importance

Sickle cell disease (SCD) is a monogenic disorder, yet clinical outcomes are influenced by additional genetic factors. Despite decades of research, the genetics of SCD remain poorly understood.

Objective

To assess all reported genetic modifiers of SCD, evaluate the design of associated studies, and provide guidelines for future analyses according to modern genetic study recommendations.

Data Sources

PubMed, Web of Science, and Scopus were searched through May 16, 2023, identifying 5290 publications.

Study Selection

At least 2 reviewers identified 571 original, peer-reviewed English-language publications reporting genetic modifiers of human SCD phenotypes, wherein the outcome was not treatment response, and the comparison was not between SCD subtypes or including healthy controls.

Data Extraction and Synthesis

Data relevant to all genetic modifiers of SCD were extracted, evaluated, and presented following STREGA and PRISMA guidelines. Weighted z score meta-analyses and pathway analyses were conducted.

Main Outcomes and Measures

Outcomes were aggregated into 25 categories, grouped as acute complications, chronic conditions, hematologic parameters or biomarkers, and general or mixed measures of SCD severity.

Results

The 571 included studies reported on 29 670 unique individuals (50% ≤ 18 years of age) from 43 countries. Of the 17 757 extracted results (4890 significant) in 1552 genes, 3675 results met the study criteria for meta-analysis: reported phenotype and genotype, association size and direction, variability measure, sample size, and statistical test. Only 173 results for 62 associations could be cross-study combined. The remaining associations could not be aggregated because they were only reported once or methods (eg, study design, reporting practice) and genotype or phenotype definitions were insufficiently harmonized. Gene variants regulating fetal hemoglobin and α-thalassemia (important markers for SCD severity) were frequently identified: 19 single-nucleotide variants in BCL11A, HBS1L-MYB, and HBG2 were significantly associated with fetal hemoglobin (absolute value of Z = 4.00 to 20.66; P = 8.63 × 10-95 to 6.19 × 10-5), and α-thalassemia deletions were significantly associated with increased hemoglobin level and reduced risk of albuminuria, abnormal transcranial Doppler velocity, and stroke (absolute value of Z = 3.43 to 5.16; P = 2.42 × 10-7 to 6.00 × 10-4). However, other associations remain unconfirmed. Pathway analyses of significant genes highlighted the importance of cellular adhesion, inflammation, oxidative and toxic stress, and blood vessel regulation in SCD (23 of the top 25 Gene Ontology pathways involve these processes) and suggested future research areas.

Conclusions and Relevance

The findings of this comprehensive systematic review and meta-analysis of all published genetic modifiers of SCD indicated that implementation of standardized phenotypes, statistical methods, and reporting practices should accelerate discovery and validation of genetic modifiers and development of clinically actionable genetic profiles.

Genome-wide meta-analysis identifies new candidate genes for sickle cell disease nephropathy

Sickle cell disease nephropathy (SCDN), a common SCD complication, is strongly associated with mortality. Polygenic risk scores calculated from recent transethnic meta-analyses of urinary albumin-to-creatinine ratio and estimated glomerular filtration rate (eGFR) trended toward association with proteinuria and eGFR in SCD but the model fit was poor (R2 < 0.01), suggesting that there are likely unique genetic risk factors for SCDN. Therefore, we performed genome-wide association studies (GWAS) for 2 critical manifestations of SCDN, proteinuria and decreased eGFR, in 2 well-characterized adult SCD cohorts, representing, to the best of our knowledge, the largest SCDN sample to date. Meta-analysis identified 6 genome-wide significant associations (false discovery rate, q ≤ 0.05): 3 for proteinuria (CRYL1, VWF, and ADAMTS7) and 3 for eGFR (LRP1B, linc02288, and FPGT-TNNI3K/TNNI3K). These associations are independent of APOL1 risk and represent novel SCDN loci, many with evidence for regulatory function. Moreover, GWAS SNPs in CRYL1, VWF, ADAMTS7, and linc02288 are associated with gene expression in kidney and pathways important to both renal function and SCD biology, supporting the hypothesis that SCDN pathophysiology is distinct from other forms of kidney disease. Together, these findings provide new targets for functional follow-up that could be tested prospectively and potentially used to identify patients with SCD who are at risk, before onset of kidney dysfunction.

Determinants of severity in sickle cell disease

Sickle cell disease is a very variable condition, with outcomes ranging from death in childhood to living relatively symptom free into the 8^th decade. Much of this variability is unexplained. The co-inheritance of α thalassaemia and factors determining HbF levels significantly modify the phenotype, but few other significant genetic variants have been identified, despite extensive studies. Environmental factors are undoubtedly important, with socio-economics and access to basic medical care explaining the huge differences in outcomes between many low- and high-income countries. Exposure to cold and windy weather seems to precipitate acute complications in many people, although these effects are unpredictable and vary with geography. Many studies have tried to identify prognostic factors which can be used to predict outcomes, particularly when applied in infancy. Overall, low haemoglobin, low haemoglobin F percentage and high reticulocytes in childhood are associated with worse outcomes, although again these effects are fairly weak and inconsistent.

A polygenic score for acute vaso-occlusive pain in pediatric sickle cell disease

Individuals with monogenic disorders can experience variable phenotypes that are influenced by genetic variation. To investigate this in sickle cell disease (SCD), we performed whole-genome sequencing (WGS) of 722 individuals with hemoglobin HbSS or HbSβ0-thalassemia from Baylor College of Medicine and from the St. Jude Children's Research Hospital Sickle Cell Clinical Research and Intervention Program (SCCRIP) longitudinal cohort study. We developed pipelines to identify genetic variants that modulate sickle hemoglobin polymerization in red blood cells and combined these with pain-associated variants to build a polygenic score (PGS) for acute vaso-occlusive pain (VOP). Overall, we interrogated the α-thalassemia deletion -α3.7 and 133 candidate single-nucleotide polymorphisms (SNPs) across 66 genes for associations with VOP in 327 SCCRIP participants followed longitudinally over 6 years. Twenty-one SNPs in 9 loci were associated with VOP, including 3 (BCL11A, MYB, and the β-like globin gene cluster) that regulate erythrocyte fetal hemoglobin (HbF) levels and 6 (COMT, TBC1D1, KCNJ6, FAAH, NR3C1, and IL1A) that were associated previously with various pain syndromes. An unweighted PGS integrating all 21 SNPs was associated with the VOP event rate (estimate, 0.35; standard error, 0.04; P = 5.9 × 10-14) and VOP event occurrence (estimate, 0.42; standard error, 0.06; P = 4.1 × 10-13). These associations were stronger than those of any single locus. Our findings provide insights into the genetic modulation of VOP in children with SCD. More generally, we demonstrate the utility of WGS for investigating genetic contributions to the variable expression of SCD-associated morbidities.

Diagnosis of Sickle Cell Disease and HBB Haplotyping in the Era of Personalized Medicine: Role of Next Generation Sequencing

Hemoglobin genotype and HBB haplotype are established genetic factors that modify the clinical phenotype in sickle cell disease (SCD). Current methods of establishing these two factors are cumbersome and/or prone to errors. The throughput capability of next generation sequencing (NGS) makes it ideal for simultaneous interrogation of the many genes of interest in SCD. This study was designed to confirm the diagnosis in patients with HbSS and Sβ-thalassemia, identify any ß-thal mutations and simultaneously determine the ß^S HBB haplotype. Illumina Ampliseq custom DNA panel was used to genotype the DNA samples. Haplotyping was based on the alleles on five haplotype-specific SNPs. The patients studied included 159 HbSS patients and 68 Sβ-thal patients, previously diagnosed using high performance liquid chromatography (HPLC). There was considerable discordance between HPLC and NGS results, giving a false +ve rate of 20.5% with a sensitivity of 79% for the identification of Sβthal. Arab/India haplotype was found in 81.5% of β^S chromosomes, while the two most common, of the 13 β-thal mutations detected, were IVS-1 del25 and IVS-II-1 (G>A). NGS is very versatile and can be deployed to simultaneously screen multiple gene loci for modifying polymorphisms, to afford personalized, evidence-based counselling and early intervention.

Plasma Metabolomic Signatures of Chronic Obstructive Pulmonary Disease and the Impact of Genetic Variants on Phenotype-Driven Modules

Background: Small studies have recently suggested that there are specific plasma metabolic signatures in chronic obstructive pulmonary disease (COPD), but there have been no large comprehensive study of metabolomic signatures in COPD that also integrate genetic variants. Materials and Methods: Fresh frozen plasma from 957 non-Hispanic white subjects in COPDGene was used to quantify 995 metabolites with Metabolon's global metabolomics platform. Metabolite associations with five COPD phenotypes (chronic bronchitis, exacerbation frequency, percent emphysema, post-bronchodilator forced expiratory volume at one second [FEV1]/forced vital capacity [FVC], and FEV1 percent predicted) were assessed. A metabolome-wide association study was performed to find genetic associations with metabolite levels. Significantly associated single-nucleotide polymorphisms were tested for replication with independent metabolomic platforms and independent cohorts. COPD phenotype-driven modules were identified in network analysis integrated with genetic associations to assess gene-metabolite-phenotype interactions. Results: Of metabolites tested, 147 (14.8%) were significantly associated with at least 1 COPD phenotype. Associations with airflow obstruction were enriched for diacylglycerols and branched chain amino acids. Genetic associations were observed with 109 (11%) metabolites, 72 (66%) of which replicated in an independent cohort. For 20 metabolites, more than 20% of variance was explained by genetics. A sparse network of COPD phenotype-driven modules was identified, often containing metabolites missed in previous testing. Of the 26 COPD phenotype-driven modules, 6 contained metabolites with significant met-QTLs, although little module variance was explained by genetics. Conclusion: A dysregulation of systemic metabolism was predominantly found in COPD phenotypes characterized by airflow obstruction, where we identified robust heritable effects on individual metabolite abundances. However, network analysis, which increased the statistical power to detect associations missed previously in classic regression analyses, revealed that the genetic influence on COPD phenotype-driven metabolomic modules was modest when compared with clinical and environmental factors.

Data science and precision health care

Precision health care plays a crucial role in an elderly society by providing personalized health care plans for improving an individual's health conditions and preventing disease. To realize precision health care, data science is key; it allows for analyses of health-related big data. In this article, an actual analysis of time-series health check-up data is presented and as is a discussion of how personalized simulation models of health conditions are constructed and used to modify individual behavior. Future directions for precision health care based on the integration of genetic variations and the microbiome are also discussed.

A Genome-wide Association Study of Circulating Levels of Atorvastatin and Its Major Metabolites

Atorvastatin (ATV) is frequently prescribed and generally well  tolerated, but can lead to myotoxicity, especially at higher doses. A genome-wide association study of circulating levels of ATV, 2-hydroxy (2-OH) ATV, ATV lactone (ATV L), and 2-OH ATV L was performed in 590 patients who had been hospitalized with a non-ST elevation acute coronary syndrome 1 month earlier and were on high-dose ATV (80 mg or 40 mg daily). The UGT1A locus (lead single nucleotide polymorphism, rs887829) was strongly associated with both increased 2-OH ATV/ATV (P = 7.25 × 10^-16 ) and 2-OH ATV L/ATV L (P = 3.95 × 10^-15 ) metabolic ratios. Moreover, rs45446698, which tags CYP3A7*1C, was nominally associated with increased 2-OH ATV/ATV (P = 6.18 × 10^-7 ), and SLCO1B1 rs4149056 with increased ATV (P = 2.21 × 10^-6 ) and 2-OH ATV (P = 1.09 × 10^-6 ) levels. In a subset of these patients whose levels of ATV and metabolites had also been measured at 12 months after hospitalization (n = 149), all of these associations remained, except for 2-OH ATV and rs4149056 (P = 0.057). Clinically, rs4149056 was associated with increased muscular symptoms (odds ratio (OR) 3.97; 95% confidence interval (CI) 1.29-12.27; P = 0.016) and ATV intolerance (OR 1.55; 95% CI 1.09-2.19; P = 0.014) in patients (n = 870) primarily discharged on high-dose ATV. In summary, both novel and recognized genetic associations have been identified with circulating levels of ATV and its major metabolites. Further study is warranted to determine the clinical utility of genotyping rs4149056 in patients on high-dose ATV.

Prediction of blood test values under different lifestyle scenarios using time-series electronic health record

Owing to increasing medical expenses, researchers have attempted to detect clinical signs and preventive measures of diseases using electronic health record (EHR). In particular, time-series EHRs collected by periodic medical check-up enable us to clarify the relevance among check-up results and individual environmental factors such as lifestyle. However, usually such time-series data have many missing observations and some results are strongly correlated to each other. These problems make the analysis difficult and there exists strong demand to detect clinical findings beyond them. We focus on blood test values in medical check-up results and apply a time-series analysis methodology using a state space model. It can infer the internal medical states emerged in blood test values and handle missing observations. The estimated models enable us to predict one's blood test values under specified condition and predict the effect of intervention, such as changes of body composition and lifestyle. We use time-series data of EHRs periodically collected in the Hirosaki cohort study in Japan and elucidate the effect of 17 environmental factors to 38 blood test values in elderly people. Using the estimated model, we then simulate and compare time-transitions of participant's blood test values under several lifestyle scenarios. It visualizes the impact of lifestyle changes for the prevention of diseases. Finally, we exemplify that prediction errors under participant's actual lifestyle can be partially explained by genetic variations, and some of their effects have not been investigated by traditional association studies.

Uridine diphosphate glucuronosyl transferase 1A (UGT1A1) promoter polymorphism in young patients with sickle cell anaemia: report of the first cohort study from Nigeria

Background

(TA) n repeat sequence (rs8175347) of UGT1A1 gene promoter polymorphism is associated with serum bilirubin levels and gallstones among different sickle cell anaemia (SCA) populations. There are no data on UGT1A1 polymorphisms and their impact on Nigerian SCA patients. In this study, we determined the distribution of the UGT1A1 (TA) n genotypes among a group of young Nigerian SCA patients and healthy controls. In addition, the influence of UGT1A1 (TA) n genotypes on the laboratory and clinical events among the patients was determined.

Methods

The distribution of the UGT1A1 (TA) n genotypes among 101 young Nigerian SCA patients and 64 normal appropriate controls were determined and studied. The UGT1A1 (TA) n genotypes were further classified into subgroups and used to differentiate the clinical events and laboratory parameters of the patients.

Results

Four (TA) n alleles:(TA)5, 6, 7, and 8 were found. These were associated with 10 genotypes: TA5/5, 5/6, 5/7, 5/8, 6/6, 6/7, 6/8, 7/7, 7/8, 8/8. The normal (wild-type)-(TA) 6/6), low- (TA) 7/7, 7/8, 8/8), intermediate- (TA) 5/7, 5/8, 6/7, 6/8), and high-activity (TA) 5/5, 5/6,) genotypes were found in 24.8, 24.8, 41.5, and 8.9% patients and 20.3, 15.6, 61, and 3.1% controls respectively. The general genotype distribution of the patients and control group were not significantly different. There were significant differences in serum bilirubin and lactate dehydrogenase (LDH) of the patients when differentiated by the UGT1A1 (TA) n genotypes (p<0.05). Asymptomatic gallstones were found in 5.9% of patients and were significantly of the low-activity genotypes sub-group 5 (20%) vs 1(1.3%) p = 0.0033. Although, bilirubin and fetal hemoglobin (HbF) of patients with gallstones were significantly different from those without gallstone, only the serum bilirubin was associated with UGT1A1 (TA) n genotypes on multivariate analysis (p < 0.0001).

Conclusion

This study highlights the contribution of UGT1A1 polymorphisms, a non-globin genetic factor, to the laboratory and clinical manifestations of young Nigerian SCA patients for the first time. It also shows that children with co-inheritance of low UGT1A1 (TA) n affinity genotypes may be at risk of gallstone, hence the need to follow them up.

Incidence and Risk of Gallstone Disease in Gilbert's Syndrome Patients in Indian Population

BACKGROUND/OBJECTIVES

Individuals with Gilbert's syndrome (GS) harbor mutations in the UGT1A1 gene and are known to have elevated levels of bilirubin, which enhances the risk for gall stone formation. The aim of this study is to screen Indian patients with GS for the incidence of gall stone disease.

METHODS

Individuals with persistently elevated serum bilirubin levels were genotyped for two polymorphisms (rs8175347; rs4148323) in UGT1A1 gene to confirm GS in them. Flanking regions of the above polymorphisms were amplified followed by direct sequencing. Ultrasonography was done to detect gallstone disease. Clinical data, including assessment of liver function, circulating levels of total and direct bilirubin, as well as routine hematological parameters were obtained as per standard procedures (Autoanalyzer).

RESULTS

Of the total 1621 individuals subjected to genotyping, 1191 (1149 males of 29.6 ± 11.3 years with mean BMI of 22.1 ± 3.7 kg/m² and 42 females of 30.8 ± 14.8 years with mean BMI of 20.9 ± 3.7 kg/m²) were confirmed to have GS. Gall bladder abnormalities including cholelithiasis (n = 106/1191; 8.9%), polyps (n = 18/1191; 1.5%) and gallbladder wall thickening (n = 17/1191; 1.4%) were noted. Incidence of gall stone disease was observed in 103 males (out of 1149) and 3 females (out of 42) indicating the risk of the disease to be 9.0% and 7.1% respectively in males and females with GS.

CONCLUSION

Early recognition of GS by genetic analysis is required before these patients with intermittent episodes of jaundice run the risk of unnecessary operations on their bile ducts from the mistaken assumption ascribing the jaundice to a stone which has been left behind.

Sickle Cell Anemia and Its Phenotypes

In the 100 years since sickle cell anemia (SCA) was first described in the medical literature, studies of its molecular and pathophysiological basis have been at the vanguard of scientific discovery. By contrast, the translation of such knowledge into treatments that improve the lives of those affected has been much too slow. Recent years, however, have seen major advances on several fronts. A more detailed understanding of the switch from fetal to adult hemoglobin and the identification of regulators such as BCL11A provide hope that these findings will be translated into genomic-based approaches to the therapeutic reactivation of hemoglobin F production in patients with SCA. Meanwhile, an unprecedented number of new drugs aimed at both the treatment and prevention of end-organ damage are now in the pipeline, outcomes from potentially curative treatments such as allogeneic hematopoietic stem cell transplantation are improving, and great strides are being made in gene therapy, where methods employing both antisickling β-globin lentiviral vectors and gene editing are now entering clinical trials. Encouragingly, after a century of neglect, the profile of the vast majority of those with SCA in Africa and India is also finally improving.

Genetic modifiers of severity in sickle cell disease

Sickle cell disease (SCD) is one of the most common single disease disorders world-wide. It is remarkable for its clinical heterogeneity, even among individuals with identical genotypes. Some individuals experience morbidity and mortality in early childhood, while others have a relatively mild course, and normal or near normal life expectancy. Many clinical complications are associated with SCD; most notably frequent pain episodes, stroke, acute chest syndrome, avascular necrosis, nephropathy, retinopathy and pulmonary hypertension. While the effects of higher fetal hemoglobin (HbF) levels, UGTA1A polymorphisms, alpha-thalassemia and G6PD deficiency on SCD has been extensively studied, these variables do not explain all of the clinical heterogeneity of SCD. It is not known why some patients develop certain complications, and it is difficult to predict which complications a particular patient will experience. Much work has been done to identify genetic variants associated with these disease complications; many associations remain unvalidated. As the field continues to move beyond small sample collections and candidate gene approaches into whole genome sequencing and merging of samples from all over the world, we will identify more genetic variants associated with development of specific SCD related complications, and hopefully leverage this knowledge into targeted therapies.

Biomarker signatures of sickle cell disease severity

Identifying sickle cell disease patients at high risk of complications could lead to personalized treatment and better prognosis but despite many advances prediction of the clinical course of these patients remains elusive. We propose a system-type approach to discover profiles of multiple, common biomarkers that correlate with morbidity and mortality in sickle cell disease. We used cluster analysis to discover 17 signatures of 17 common circulating biomarkers in 2320 participants of the Cooperative Study of Sickle Cell Disease, and evaluated the association of these signatures with risk for stroke, pain, leg ulceration, acute chest syndrome, avascular necrosis, seizure, death, and trend of fetal hemoglobin and hemolysis using longitudinally collected data. The analysis shows that some of the signatures are associated with reduced risk for complications, while others are associated with increased risk for complications. We also show that these signatures repeat in two more contemporary studies of sickle cell disease and correlate with recently discovered biomarkers of pulmonary vascular disease. With replication and further study, these biomarker signatures could become an important and affordable precision medicine tool to aid treatment and management of the disease.

Multiple Testing in the Context of Gene Discovery in Sickle Cell Disease Using Genome-Wide Association Studies

The issue of multiple testing, also termed multiplicity, is ubiquitous in studies where multiple hypotheses are tested simultaneously. Genome-wide association study (GWAS), a type of genetic association study that has gained popularity in the past decade, is most susceptible to the issue of multiple testing. Different methodologies have been employed to address the issue of multiple testing in GWAS. The purpose of the review is to examine the methodologies employed in dealing with multiple testing in the context of gene discovery using GWAS in sickle cell disease complications.

UGT1A1 gene linkage analysis: application of polymorphic markers rs4148326/rs4124874 in the Iranian population

OBJECTIVES

Mutations in the UGT1A1 gene are responsible for hyperbilirubinemia syndromes including Crigler-Najjar type 1 and 2 and Gilbert syndrome. In view of the genetic heterogeneity and involvement of large numbers of the disease causing mutations, the application of polymorphic markers in the UGTA1 gene could be useful in molecular diagnosis of the disease.

MATERIALS AND METHODS

In the present study, two polymorphic markers including rs4148326 and rs4124874 in the UGT1A1 gene region were characterized. The markers were selected using bioinformatics analysis of the UGT1A1 gene region and genotyped in 212 unrelated healthy individuals and 13 family trios in the Iranian population using Tetra-Primer ARMS PCR technique. The allele frequency and population status of the alleles were estimated using GENEPOP, FBAT, PowerMarker and Arlequin software.

RESULTS

The results indicated that in the case of rs4148326 marker, allele frequency for T and C allele was 66.04% and 33.96%, respectively. For rs4124874 marker, allele frequency for G and T alleles was 39.4% and 60.6%, respectively. The values of heterozygosity index for the markers examined were 64.1 for rs4148326 and 72.1 for rs4124874, respectively. The haplotype estimation analysis of the markers resulted in three informative haplotypes with frequencies ≥0.05. Moreover, the results suggested the presence of linkage disequilibrium between two markers.

CONCLUSION

Altogether, the data suggested that rs4148326 and rs4124874 could be introduced as informative markers for molecular diagnosis of Crigler-Najjar type 1 and 2 and Gilbert syndrome in the Iranian population.

Do Alpha Thalassemia, Fetal Hemoglobin, and the UGT1A1 Polymorphism have an Influence on Serum Bilirubin Levels and Cholelithiasis in Patients with Sickle Cell Disease?

BACKGROUND

Increased destruction of erythrocytes in patients with sickle cell disease results in chronic hyperbilirubinemia and leads to the formation of gallstones.

OBJECTIVES

The objective of this study was to determine the combined influence of alpha thalassemia, fetal hemoglobin, and the UGT1A1 polymorphism on serum bilirubin levels and cholelithiasis in patients with sickle cell disease.

METHODS

We analyzed 72 patients treated in the outpatient hematology unit of the Clinical Hospital of Porto Alegre. The alpha thalassemia trait was determined by multiplex polymerase chain reaction and the polymorphisms of UGT1A1 by capillary electrophoresis with tagged primers.

RESULTS

Total and indirect bilirubin levels differed significantly between genotypes TA7/TA7 and TA6/TA6 (p < 0.05). Bilirubin levels were influenced by the UGT1A1 polymorphism but not by alpha thalassemia and fetal hemoglobin. There was no association between cholelithiasis and any of the variables studied.

CONCLUSION

These preliminary findings suggest that the UGT1A1 gene can influence serum bilirubin levels in sickle cell anemia and serve as a tool to differentiate an acute hemolytic condition from a pre-existing condition of hyperbilirubinemia.

UGT1A1 (TA)n genotype is not the major risk factor of cholelithiasis in sickle cell disease children

OBJECTIVES

Because of the increased hemolytic rate, a significant proportion of patients with sickle cell disease (SCD) are prone to develop cholelithiasis. The present study investigated the role of several genetic factors (UGT1A1 promoter (TA)_n repeat polymorphism, alpha-globin status), hematological parameters, clinical severity, and hydroxyurea (HU) therapy on the occurrence of cholelithiasis in SCD.

METHODS

One hundred and fifty-eight children (2-18 yr old) regularly followed at the University Hospital of Lyon (France) were included. A multivariate Cox model was used to test the associations between cholelithiasis and the different parameters analyzed.

RESULTS

We confirmed that alpha-thalassemia and low basal reticulocyte (RET) count were independent protective factors for cholelithiasis while 7/7, 8/8 and 7/8 UGT1A1 (TA)_n genotypes were independent predisposing factors for this complication. We also showed for the first time that HU treatment decreased the risk for cholelithiasis while frequent vaso-occlusive crises and acute chest syndrome events increased that risk.

CONCLUSIONS

Our findings demonstrate that UGT1A1 (TA)_n polymorphism is not the only factor triggering gallstone formation in SCD. Cholelithiasis is also modulated by RET count, the number of deleted alpha-genes, HU therapy and the frequency of vaso-occlusive events.

Genetic Basis and Genetic Modifiers of β-Thalassemia and Sickle Cell Disease

β-thalassemia and sickle cell disease (SCD) are prototypical Mendelian single gene disorders, both caused by mutations affecting the adult β-globin gene. Despite the apparent genetic simplicity, both disorders display a remarkable spectrum of phenotypic severity and share two major genetic modifiers-α-globin genotype and innate ability to produce fetal hemoglobin (HbF, α₂γ₂).This article provides an overview of the genetic basis for SCD and β-thalassemia, and genetic modifiers identified through phenotype correlation studies. Identification of the genetic variants modifying HbF production in combination with α-globin genotype provide some prediction of disease severity for β-thalassemia and SCD but generation of a personalized genetic risk score to inform prognosis and guide management requires a larger panel of genetic modifiers yet to be discovered.Nonetheless, genetic studies have been successful in characterizing some of the key variants and pathways involved in HbF regulation, providing new therapeutic targets for HbF reactivation.

Emerging point-of-care technologies for sickle cell disease screening and monitoring

INTRODUCTION

Sickle Cell Disease (SCD) affects 100,000 Americans and more than 14 million people globally, mostly in economically disadvantaged populations, and requires early diagnosis after birth and constant monitoring throughout the life-span of the patient. Areas covered: Early diagnosis of SCD still remains a challenge in preventing childhood mortality in the developing world due to requirements of skilled personnel and high-cost of currently available modalities. On the other hand, SCD monitoring presents insurmountable challenges due to heterogeneities among patient populations, as well as in the same individual longitudinally. Here, we describe emerging point-of-care micro/nano platform technologies for SCD screening and monitoring, and critically discuss current state of the art, potential challenges associated with these technologies, and future directions. Expert commentary: Recently developed microtechnologies offer simple, rapid, and affordable screening of SCD and have the potential to facilitate universal screening in resource-limited settings and developing countries. On the other hand, monitoring of SCD is more complicated compared to diagnosis and requires comprehensive validation of efficacy. Early use of novel microdevices for patient monitoring might come in especially handy in new clinical trial designs of emerging therapies.

Minireview: Clinical severity in sickle cell disease: the challenges of definition and prognostication

Sickle cell disease (SCD) is a monogenic, yet highly phenotypically variable disease with multisystem pathology. This manuscript provides an overview of many of the known determinants, modifiers, and correlates of disease severity in SCD. Despite this wealth of data, modeling the variable and multisystem pathology of SCD continues to be difficult. The current status of prediction of specific adverse outcomes and global disease severity in SCD is also reviewed, highlighting recent successes and ongoing challenges.

Minireview: Genetic basis of heterogeneity and severity in sickle cell disease

Sickle cell disease, a common single gene disorder, has a complex pathophysiology that at its root is initiated by the polymerization of deoxy sickle hemoglobin. Sickle vasoocclusion and hemolytic anemia drive the development of disease complications. In this review, we focus on the genetic modifiers of disease heterogeneity. The phenotypic heterogeneity of disease is only partially explained by genetic variability of fetal hemoglobin gene expression and co-inheritance of α thalassemia. Given the complexity of pathophysiology, many different definitions of severity are possible complicating a full understanding of its genetic foundation. The pathophysiological complexity and the interlocking nature of the biological processes underpinning disease severity are becoming better understood. Nevertheless, useful genetic signatures of severity, regardless of how this is defined, are insufficiently developed to be used for treatment decisions and for counseling.

Genomic approaches to identifying targets for treating β hemoglobinopathies

Sickle cell disease and β thalassemia are common severe diseases with little effective pathophysiologically-based treatment. Their phenotypic heterogeneity prompted genomic approaches to identify modifiers that ultimately might be exploited therapeutically. Fetal hemoglobin (HbF) is the major modulator of the phenotype of the β hemoglobinopathies. HbF inhibits deoxyHbS polymerization and in β thalassemia compensates for the reduction of HbA. The major success of genomics has been a better understanding the genetic regulation of HbF by identifying the major quantitative trait loci for this trait. If the targets identified can lead to means of increasing HbF to therapeutic levels in sufficient numbers of sickle or β-thalassemia erythrocytes, the pathophysiology of these diseases would be reversed. The availability of new target loci, high-throughput drug screening, and recent advances in genome editing provide the opportunity for new approaches to therapeutically increasing HbF production.

Genetics of human metabolism: an update

Genome-wide association studies with metabolomics (mGWAS) identify genetically influenced metabotypes (GIMs), their ensemble defining the heritable part of every human's metabolic individuality. Knowledge of genetic variation in metabolism has many applications of biomedical and pharmaceutical interests, including the functional understanding of genetic associations with clinical end points, design of strategies to correct dysregulations in metabolic disorders and the identification of genetic effect modifiers of metabolic disease biomarkers. Furthermore, it has been shown that GIMs provide testable hypotheses for functional genomics and metabolomics and for the identification of novel gene functions and metabolite identities. mGWAS with growing sample sizes and increasingly complex metabolic trait panels are being conducted, allowing for more comprehensive and systems-based downstream analyses. The generated large datasets of genetic associations can now be mined by the biomedical research community and provide valuable resources for hypothesis-driven studies. In this review, we provide a brief summary of the key aspects of mGWAS, followed by an update of recently published mGWAS. We then discuss new approaches of integrating and exploring mGWAS results and finish by presenting selected applications of GIMs in recent studies.

Exome-Wide Association Study Identifies New Low-Frequency and Rare UGT1A1 Coding Variants and UGT1A6 Coding Variants Influencing Serum Bilirubin in Elderly Subjects: A Strobe Compliant Article

Genome-wide association studies (GWASs) have identified loci contributing to total serum bilirubin level. However, no exome-wide approaches have been performed to address this question. Using exome-wide approach, we assessed the influence of protein-coding variants on unconjugated, conjugated, and total serum bilirubin levels in a well-characterized cohort of 773 ambulatory elderly subjects from Italy. Coding variants were replicated in 227 elderly subjects from the same area. We identified 4 missense rare (minor allele frequency, MAF < 0.5%) and low-frequency (MAF, 0.5%-5%) coding variants located in the first exon of the UGT1A1 gene, which encodes for the substrate-binding domain (rs4148323 [MAF = 0.06%; p.Gly71Arg], rs144398951 [MAF = 0.06%; p.Ile215Val], rs35003977 [MAF = 0.78%; p.Val225Gly], and rs57307513 [MAF = 0.06%; p.Ser250Pro]). These variants were in strong linkage disequilibrium with 3 intronic UGT1A1 variants (rs887829, rs4148325, rs6742078), which were significantly associated with total bilirubin level (P = 2.34 × 10(-34), P = 7.02 × 10(-34), and P = 8.27 × 10(-34)), as well as unconjugated, and conjugated bilirubin levels. We also identified UGT1A6 variants in association with total (rs6759892, p.Ser7Ala, P = 1.98 × 10(-26); rs2070959, p.Thr181Ala, P = 2.87 × 10(-27); and rs1105879, p.Arg184Ser, P = 3.27 × 10(-29)), unconjugated, and conjugated bilirubin levels. All UGT1A1 intronic variants (rs887829, rs6742078, and rs4148325) and UGT1A6 coding variants (rs6759892, rs2070959, and rs1105879) were significantly associated with gallstone-related cholecystectomy risk. The UGT1A6 variant rs2070959 (p.Thr181Ala) was associated with the highest risk of gallstone-related cholecystectomy (OR, 4.58; 95% CI, 1.58-13.28; P = 3.21 × 10(-3)). Using an exome-wide approach we identified coding variants on UGT1A1 and UGT1A6 genes in association with serum bilirubin level and hyperbilirubinemia risk in elderly subjects. UGT1A1 intronic single-nucleotide polymorphisms (SNPs) (rs6742078, rs887829, rs4148324) serve as proxy markers for the low-frequency and rare UGT1A1 variants, thereby providing mechanistic explanation to the relationship between UGT1A1 intronic SNPs and the UGT1A1 enzyme activity. UGT1A1 and UGT1A6 variants might be potentially associated with gallstone-related cholecystectomy risk.

NRMT1 knockout mice exhibit phenotypes associated with impaired DNA repair and premature aging

Though defective genome maintenance and DNA repair have long been known to promote phenotypes of premature aging, the role protein methylation plays in these processes is only now emerging. We have recently identified the first N-terminal methyltransferase, NRMT1, which regulates protein-DNA interactions and is necessary for both accurate mitotic division and nucleotide excision repair. To demonstrate if complete loss of NRMT1 subsequently resulted in developmental or aging phenotypes, we constructed the first NRMT1 knockout (Nrmt1(-/-)) mouse. The majority of these mice die shortly after birth. However, the ones that survive, exhibit decreased body size, female-specific infertility, kyphosis, decreased mitochondrial function, and early-onset liver degeneration; phenotypes characteristic of other mouse models deficient in DNA repair. The livers from Nrmt1(-/-) mice produce less reactive oxygen species (ROS) than wild type controls, and Nrmt1(-/-) mouse embryonic fibroblasts show a decreased capacity for handling oxidative damage. This indicates that decreased mitochondrial function may benefit Nrmt1(-/-) mice and protect them from excess internal ROS and subsequent DNA damage. These studies position the NRMT1 knockout mouse as a useful new system for studying the effects of genomic instability and defective DNA damage repair on organismal and tissue-specific aging.

Genome-wide association studies in Africans and African Americans: expanding the framework of the genomics of human traits and disease

Genomic research is one of the tools for elucidating the pathogenesis of diseases of global health relevance and paving the research dimension to clinical and public health translation. Recent advances in genomic research and technologies have increased our understanding of human diseases, genes associated with these disorders, and the relevant mechanisms. Genome-wide association studies (GWAS) have proliferated since the first studies were published several years ago and have become an important tool in helping researchers comprehend human variation and the role genetic variants play in disease. However, the need to expand the diversity of populations in GWAS has become increasingly apparent as new knowledge is gained about genetic variation. Inclusion of diverse populations in genomic studies is critical to a more complete understanding of human variation and elucidation of the underpinnings of complex diseases. In this review, we summarize the available data on GWAS in recent African ancestry populations within the western hemisphere (i.e. African Americans and peoples of the Caribbean) and continental African populations. Furthermore, we highlight ways in which genomic studies in populations of recent African ancestry have led to advances in the areas of malaria, HIV, prostate cancer, and other diseases. Finally, we discuss the advantages of conducting GWAS in recent African ancestry populations in the context of addressing existing and emerging global health conditions.

Combined effects of smoking and bilirubin levels on the risk of lung cancer in Korea: the severance cohort study

BACKGROUND

Smoking is a major risk factor for lung cancer. Bilirubin, an antioxidant, is inversely associated with the risk of diseases related to oxidative stress. This study was conducted to determine the influence of smoking and bilirubin levels on the risk of lung cancer in the Severance cohort study.

METHODS

This study included 68,676 Korean who received a health examination at Severance Health Promotion Center from 1994 to 2004. Serum bilirubin measurements within normal range were divided into tertiles whereas smoking states were divided as never-smokers, former smokers and current smokers. A diagnosis of lung cancer was coded as occurring based on the report from the National Cancer Registry. Hazard ratios (HRs) and 95% confidence intervals (95% CIs) were calculated using Cox proportional hazards model.

RESULTS

At the end of the study period, 240 patients (men: 181, women: 59) developed lung cancer. Compared to those with bilirubin levels ≥ 1.0 mg/dL, HRs (95% CI) for lung cancer were 2.8 (1.8-4.2) for subjects having bilirubin levels from 0.2 to 0.7 mg/dL in men. When we stratified our analysis by smoking status, bilirubin consistently showed a protective effect on the risk of lung cancer on both never- and current smokers. Current smokers having bilirubin levels from 0.2 to 0.7 mg/dL had a risk of lung cancer by 6.0-fold higher than never-smokers with bilirubin levels ≥ 1.0 mg/dL in men.

CONCLUSION

In this large prospective study, higher baseline bilirubin level in the normal range was associated with low risk of lung cancer. Smoking and low bilirubin levels were cumulatively associated with a higher risk of lung cancer.

Genetic factors affecting gene transcription and catalytic activity of UDP-glucuronosyltransferases in human liver

The aim of this study was to discover cis- and trans-acting factors significantly affecting mRNA expression and catalytic activity of human hepatic UDP-glucuronosyltransferases (UGTs). Transcription levels of five major hepatic UGT1A (UGT1A1, UGT1A3, UGT1A4, UGT1A6 and UGT1A9) and five UGT2B (UGT2B4, UGT2B7, UGT2B10, UGT2B15 and UGT2B17) genes were quantified in human liver tissue samples (n = 125) using real-time PCR. Glucuronidation activities of 14 substrates were measured in 47 livers. We genotyped 167 tagSNPs (single-nucleotide polymorphisms) in UGT1A (n = 43) and UGT2B (n = 124), as well as the known functional UGT1A1*28 and UGT2B17 CNV (copy number variation) polymorphisms. Transcription levels of 15 transcription factors (TFs) known to regulate these UGTs were quantified. We found that UGT expression and activity were highly variable among the livers (median and range of coefficient of variations: 135%, 74-217% and 52%, 39-105%, respectively). CAR, PXR and ESR1 were found to be the most important trans-regulators of UGT transcription (median and range of correlation coefficients: 46%, 6-58%; 47%, 9-58%; and 52%, 24-75%, respectively). Hepatic UGT activities were mainly determined by UGT gene transcription levels. Twenty-one polymorphisms were significantly (FDR-adjusted P < 0.05) associated with mRNA expression and/or activities of UGT1A1, UGT1A3 and UGT2B17. We found novel SNPs in the UGT2B17 CNV region accounting for variability in UGT2B17 gene transcription and testosterone glucuronidation rate, in addition to that attributable to the UGT2B17 CNV. Our study discovered novel pharmacogenetic markers and provided detailed insight into the genetic network regulating hepatic UGTs.

Applications of high-throughput DNA sequencing to benign hematology

The development of novel technologies for high-throughput DNA sequencing is having a major impact on our ability to measure and define normal and pathologic variation in humans. This review discusses advances in DNA sequencing that have been applied to benign hematologic disorders, including those affecting the red blood cell, the neutrophil, and other white blood cell lineages. Relevant examples of how these approaches have been used for disease diagnosis, gene discovery, and studying complex traits are provided. High-throughput DNA sequencing technology holds significant promise for impacting clinical care. This includes development of improved disease detection and diagnosis, better understanding of disease progression and stratification of risk of disease-specific complications, and development of improved therapeutic strategies, particularly patient-specific pharmacogenomics-based therapy, with monitoring of therapy by genomic biomarkers.

Early complication in sickle cell anemia children due to A(TA)nTAA polymorphism at the promoter of UGT1A1 gene

Aim. To determine the implication of the polymorphism, namely, A(TA)_nTAA of UGT1A1 in lithogenesis for the first time in Tunisia among sickle cell anemia (SCA) children patients. Material and Methods. Our study was performed in 2010 and it involved 76 subjects chosen as control group characterized with normal hemoglobin status and presence of cholelithiasis and 102 SCA pediatric patients among whom 52 have cholelithiasis. We analyzed the polymorphism A(TA)_nTAA at the UGT1A1 promoter and the relationships between the various A(TA)_nTAA genotypes and alleles and bilirubin levels and occurrence of cholelithiasis. Results and Discussion. The repartition of genotypes found according to serum bilirubin level shows a significant association between genotypes carrying variant (TA)₇ and hyperbilirubinemia (P < 0.05). We demonstrated the association of two genotypes with gallstones formation among SCA children patients: (TA)₇/(TA)₇ and (TA)₇/(TA)₈ with P = 8.1 × 10⁻⁸ and P = 0.01, respectively. (TA)₇ and (TA)₈ allele variants act as a risk factor for early gallstones formation in SCA patients with P = 5.8 × 10⁻⁹ and P = 0.01, respectively. As for the control group only the genotype (TA)₇/(TA)₇ presented a risk factor for gallstones formation. Conclusion. The novelty of this report is that it is the first time that a similar study was made on the Tunisian children sickle cell population and that the results show a clear association of (TA)₇ variant in early gallstones formation in Tunisian SCA children. Interestingly our findings highlighted the association of (TA)₈ variant as well, which was not found in previous studies.

Gallstone disease in Swedish twins is associated with the Gilbert variant of UGT1A1

BACKGROUND & AIMS

The Gilbert syndrome-associated functional TATA box variant UGT1A1*28 (A(TA)7TAA) was found to increase susceptibility to pigment gallstone formation in patients with haemolytic anaemia. Further studies in extensive cohorts demonstrated an increased risk of this variant for cholesterol gallstone disease (GD). We now investigated this polymorphism as a determinant of symptomatic GD in Swedish twins.

METHODS

The Swedish Twin Registry was merged with the Hospital Discharge and Causes of Death Registries and searched for GD-related diagnoses among monozygotic (MZ) twins living in the Stockholm area. In addition, we screened the TwinGene database for GD. In total, we found 44 MZ twin pairs with and eight MZ twins without GD to be evaluable. GD-free twins from TwinGene (109 concordantly MZ and 126 independent DZ) served as controls. UGT1A1*28 genotyping was performed using TaqMan assays.

RESULTS

Overall, 58 and 8 of 106 twins with GD were hetero- and homozygous UGT1A1 risk allele carriers respectively. The case-control association tests showed a significantly (P < 0.05) increased risk of developing GD (OR = 1.62, 95% CI 1.00-2.63) in heterozygotes carriers and in addition, a trend (P = 0.075) for an increased risk among carriers (OR = 1.52, 95% CI 0.97-2.44) of the risk allele.

CONCLUSION

These data from Swedish twins confirm the Gilbert variant as risk factor for GD. Our observation is in line with nucleation in bilirubin supersaturated bile representing an initial step in cholelithogenesis.

UGT1A1 promoter polymorphism associated with serum bilirubin level in Saudi patients with sickle cell disease

BACKGROUND AND OBJECTIVES

Polymorphism in (TA)n of the UGT1A1 promoter influences bilirubin level and risk of gallstones in patients with sickle cell disease (SCD) of African descent. Modifiers of bilirubin level and gallstones in Saudi patients with SCD are not known.

DESIGN AND SETTINGS

Patients with SCD presenting to participating institutions between July 2009 and July 2012 were enrolled in our study.

METHODS

A total of 223 SCD patients were enrolled. Laboratory workup at steady state included complete blood count, reticulocytes, serum bilirubin, lactate dehydrogenase (LDH), G6PD level, and hemoglobin (Hb) electrophoresis. The (TA)n UGT1A1 promoter polymorphism and presence of a-thalassemia were also deter.mined.

RESULTS

TA6/6 in the UGT1A1 promoter was identified in 189 patients (84.7%), TA7/7 in 26 (11.7%), TA5/5 in 6 (2.7%), and TA5/6 in 2 (0.9%). Increased (TA)n of the UGT1A1 promoter (P < .0001), male gender (P=.02), higher LDH (P=.001), and lower Hb level (P=.009) were associated with higher bilirubin level, while the co-inheritance of a-thalassemia (P=.003) was linked with lower bilirubin level. UGT1A1 (TA)n (P < .0001) and Hb level (P=.005) remained significant on multivariate analysis. Gallstones were more frequent in patients with TA7/7 (72%) compared to patients with TA6/6 (57%) and TA5/5 or 5/6 (37%); however, this difference was not statistically significance (P=.18). Older age (P=.0001) and absence of a-thalassemia (P=.03) were associated with higher risk of gallstones.

CONCLUSION

(TA)n in the UGT1A1 promoter and intensity of hemolysis modify steady-state serum bilirubin level in SCD. Co-inheritance of a-thalassemia reduces the risk of gallstones in Saudi patients with SCD.

The genetics of complex cholestatic disorders

Cholestatic liver diseases are caused by a range of hepatobiliary insults and involve complex interactions among environmental and genetic factors. Little is known about the pathogenic mechanisms of specific cholestatic diseases, which has limited our ability to manage patients with these disorders. However, recent genome-wide studies have provided insight into the pathogenesis of gallstones, primary biliary cirrhosis, and primary sclerosing cholangitis. A lithogenic variant in the gene that encodes the hepatobiliary transporter ABCG8 has been identified as a risk factor for gallstone disease; this variant has been associated with altered cholesterol excretion and metabolism. Other variants of genes encoding transporters that affect the composition of bile have been associated with cholestasis, namely ABCB11, which encodes the bile salt export pump, and ABCB4, which encodes hepatocanalicular phosphatidylcholine floppase. In contrast, studies have associated primary biliary cirrhosis and primary sclerosing cholangitis with genes encoding major histocompatibility complex proteins and identified loci associated with microbial sensing and immune regulatory pathways outside this region, such as genes encoding IL12, STAT4, IRF5, IL2 and its receptor (IL2R), CD28, and CD80. These discoveries have raised interest in the development of reagents that target these gene products. We review recent findings from genetic studies of patients with cholestatic liver disease. Future characterization of genetic variants in animal models, stratification of risk alleles by clinical course, and identification of interacting environmental factors will increase our understanding of these complex cholestatic diseases.

Gene-centric association study of acute chest syndrome and painful crisis in sickle cell disease patients

Patients with sickle cell disease (SCD) present with a wide range of clinical complications. Understanding this clinical heterogeneity offers the prospects to tailor the right treatments to the right patients and also guide the development of novel therapies. Several environmental (eg, nutrition) and nonenvironmental (eg, fetal hemoglobin levels, α-thalassemia status) factors are known to modify SCD severity. To find new genetic modifiers of SCD severity, we performed a gene-centric association study in 1514 African American participants from the Cooperative Study of Sickle Cell Disease (CSSCD) for acute chest syndrome (ACS) and painful crisis. From the initial results, we selected 36 single nucleotide polymorphism (SNPs) and genotyped them for replication in 387 independent patients from the CSSCD, 318 SCD patients recruited at Georgia Health Sciences University, and 449 patients from the Duke SCD cohort. In the combined analysis, an association between ACS and rs6141803 reached array-wide significance (P = 4.1 × 10(-7)). This SNP is located 8.2 kilobases upstream of COMMD7, a gene highly expressed in the lung that interacts with nuclear factor-κB signaling. Our results provide new leads to gaining a better understanding of clinical variability in SCD, a "simple" monogenic disease.

Genetic determinants of haemolysis in sickle cell anaemia

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.

Genetic association studies in β-hemoglobinopathies

Characterization of the molecular basis of the β-thalassemias and sickle cell disease (SCD) clearly showed that individuals with the same β-globin genotypes can have extremely diverse clinical severity. Two key modifiers, an innate ability to produce fetal hemoglobin and coinheritance of α-thalassemia, both derived from family and population studies, affect the pathophysiology of both disorders at the primary level. In the past 2 decades, scientific research had applied genetic approaches to identify additional genetic modifiers. The review summarizes recent genetic studies and key genetic modifiers identified and traces the story of fetal hemoglobin genetics, which has led to an emerging network of globin gene regulation. The discoveries have provided insights on new targets for therapeutic intervention and raise possibilities of developing fetal hemoglobin predictive diagnostics for predicting disease severity in the newborn and for integration into prenatal diagnosis to better inform genetic counseling.

The search for genetic modifiers of disease severity in the β-hemoglobinopathies

Sickle cell disease (SCD) and β-thalassemia, two monogenic diseases caused by mutations in the β-globin gene, affect millions of individuals worldwide. These hemoglobin disorders are characterized by extreme clinical heterogeneity, complicating patient management and treatment. A better understanding of this patient-to-patient clinical variability would dramatically improve care and might also guide the development of novel therapies. Studies of the natural history of these β-hemoglobinopathies have identified fetal hemoglobin levels and concomitant α-thalassemia as important modifiers of disease severity. Several small-scale studies have attempted to identify additional genetic modifiers of SCD and β-thalassemia, without much success. Fortunately, improved knowledge of the human genome and the development of new genomic tools, such as genome-wide genotyping arrays and next-generation DNA sequencers, offer new opportunities to use genetics to better understand the causes of the many complications observed in β-hemoglobinopathy patients. Here I discuss the most important factors to consider when planning an experiment to find associations between β-hemoglobinopathy-related complications and DNA sequence variants, with a focus on how to successfully perform a genome-wide association study. I also review the literature and explain why most published findings in the field of SCD modifier genetics are likely to be false-positive reports, with the goal to draw lessons allowing investigators to design better genetic experiments.

Bayesian methods for multivariate modeling of pleiotropic SNP associations and genetic risk prediction

Genome-wide association studies (GWAS) have identified numerous associations between genetic loci and individual phenotypes; however, relatively few GWAS have attempted to detect pleiotropic associations, in which loci are simultaneously associated with multiple distinct phenotypes. We show that pleiotropic associations can be directly modeled via the construction of simple Bayesian networks, and that these models can be applied to produce single or ensembles of Bayesian classifiers that leverage pleiotropy to improve genetic risk prediction. The proposed method includes two phases: (1) Bayesian model comparison, to identify Single-Nucleotide Polymorphisms (SNPs) associated with one or more traits; and (2) cross-validation feature selection, in which a final set of SNPs is selected to optimize prediction. To demonstrate the capabilities and limitations of the method, a total of 1600 case-control GWAS datasets with two dichotomous phenotypes were simulated under 16 scenarios, varying the association strengths of causal SNPs, the size of the discovery sets, the balance between cases and controls, and the number of pleiotropic causal SNPs. Across the 16 scenarios, prediction accuracy varied from 90 to 50%. In the 14 scenarios that included pleiotropically associated SNPs, the pleiotropic model search and prediction methods consistently outperformed the naive model search and prediction. In the two scenarios in which there were no true pleiotropic SNPs, the differences between the pleiotropic and naive model searches were minimal. To further evaluate the method on real data, a discovery set of 1071 sickle cell disease (SCD) patients was used to search for pleiotropic associations between cerebral vascular accidents and fetal hemoglobin level. Classification was performed on a smaller validation set of 352 SCD patients, and showed that the inclusion of pleiotropic SNPs may slightly improve prediction, although the difference was not statistically significant. The proposed method is robust, computationally efficient, and provides a powerful new approach for detecting and modeling pleiotropic disease loci.

Genetic modifiers of sickle cell disease

Sickle cell anemia is associated with unusual clinical heterogeneity for a Mendelian disorder. Fetal hemoglobin concentration and coincident α thalassemia, both which directly affect the sickle erythrocyte, are the major modulators of the phenotype of disease. Understanding the genetics underlying the heritable subphenotypes of sickle cell anemia would be prognostically useful, could inform personalized therapeutics, and might help the discovery of new "druggable" pathophysiologic targets. Genotype-phenotype association studies have been used to identify novel genetic modifiers. In the future, whole genome sequencing with its promise of discovering hitherto unsuspected variants could add to our understanding of the genetic modifiers of this disease.