Fetal hemoglobin in sickle cell anemia: genome-wide association studies suggest a regulatory region in the 5' olfactory receptor gene cluster

Identification of genetic variants associated with clinical features of sickle cell disease

Sickle cell disease (SCD) is an inherited blood disorder marked by homozygosity of hemoglobin S, which is a defective hemoglobin caused by a missense mutation in the β-globin gene. However, clinical phenotypes of SCD vary among patients. To investigate genetic variants associated with various clinical phenotypes of SCD, we genotyped DNA samples from 520 SCD subjects and used a genome-wide association study (GWAS) approach to identify genetic variants associated with phenotypic features of SCD. For HbF levels, the previously reported 2p16.1 locus (BCL11A) reached genome significance (rs1427407, P = 8.58 × 10-10) in our GWAS as expected. In addition, we found a new genome-wide significance locus at 15q14 (rs8182015, P = 2.07 × 10-8) near gene EMC7. GWAS of acute chest syndrome (ACS) detected a locus (rs79915189, P = 3.70 × 10-8) near gene IDH2 at 15q26.1. The SNP, rs79915189, is also an expression quantitative trait locus (eQTL) of IDH2 in multiple tissues. For vasoocclusive episode (VOE), GWAS detected multiple significant signals at 2p25.1 (rs62118798, P = 4.27 × 10-8), 15q26.1 (rs62020555, P = 2.04 × 10-9) and 15q26.3 (rs117797325, P = 4.63 × 10-8). Our findings provide novel insights into the genetic mechanisms of SCD suggesting that common genetic variants play an important role in the presentation of the clinical phenotypes of patients with SCD.

Pharmacogenomics of Drugs Used in β-Thalassemia and Sickle-Cell Disease: From Basic Research to Clinical Applications

In this short review we have presented and discussed studies on pharmacogenomics (also termed pharmacogenetics) of the drugs employed in the treatment of β-thalassemia or Sickle-cell disease (SCD). This field of investigation is relevant, since it is expected to help clinicians select the appropriate drug and the correct dosage for each patient. We first discussed the search for DNA polymorphisms associated with a high expression of γ-globin genes and identified this using GWAS studies and CRISPR-based gene editing approaches. We then presented validated DNA polymorphisms associated with a high HbF production (including, but not limited to the HBG2 XmnI polymorphism and those related to the BCL11A, MYB, KLF-1, and LYAR genes). The expression of microRNAs involved in the regulation of γ-globin genes was also presented in the context of pharmacomiRNomics. Then, the pharmacogenomics of validated fetal hemoglobin inducers (hydroxyurea, butyrate and butyrate analogues, thalidomide, and sirolimus), of iron chelators, and of analgesics in the pain management of SCD patients were considered. Finally, we discuss current clinical trials, as well as international research networks focusing on clinical issues related to pharmacogenomics in hematological diseases.

The Impact of Patterns in Linkage Disequilibrium and Sequencing Quality on the Imprint of Balancing Selection

Regions under balancing selection are characterized by dense polymorphisms and multiple persistent haplotypes, along with other sequence complexities. Successful identification of these patterns depends on both the statistical approach and the quality of sequencing. To address this challenge, at first, a new statistical method called LD-ABF was developed, employing efficient Bayesian techniques to effectively test for balancing selection. LD-ABF demonstrated the most robust detection of selection in a variety of simulation scenarios, compared against a range of existing tests/tools (Tajima's D, HKA, Dng, BetaScan, and BalLerMix). Furthermore, the impact of the quality of sequencing on detection of balancing selection was explored, as well, using: (i) SNP genotyping and exome data, (ii) targeted high-resolution HLA genotyping (IHIW), and (iii) whole-genome long-read sequencing data (Pangenome). In the analysis of SNP genotyping and exome data, we identified known targets and 38 new selection signatures in genes not previously linked to balancing selection. To further investigate the impact of sequencing quality on detection of balancing selection, a detailed investigation of the MHC was performed with high-resolution HLA typing data. Higher quality sequencing revealed the HLA-DQ genes consistently demonstrated strong selection signatures otherwise not observed from the sparser SNP array and exome data. The HLA-DQ selection signature was also replicated in the Pangenome samples using considerably less samples but, with high-quality long-read sequence data. The improved statistical method, coupled with higher quality sequencing, leads to more consistent identification of selection and enhanced localization of variants under selection, particularly in complex regions.

CRISPR-based precision medicine for hematologic disorders: Advancements, challenges, and prospects

The development of programmable nucleases to introduce defined alterations in genomic sequences has been a powerful tool for precision medicine. While several nucleases such as zinc-finger nucleases (ZFN), transcriptor activator-like effector nucleases (TALEN), and meganucleases have been explored, the advent of CRISPR/Cas9 technology has revolutionized the field of genome engineering. In addition to disease modeling, the CRISPR/Cas9 technology has contributed to safer and more effective treatment strategies for hematologic diseases and personalized T-cell-based therapies. Here we discuss the applications of the CRISPR technology in the treatment of hematologic diseases, their efficacy, and ongoing clinical trials. We examine the obstacles to their successful use and the approaches investigated to overcome these challenges. Finally, we provide our perspectives to improve this genome editing tool for targeted therapies.

Impact of Genetic Variations on Thromboembolic Risk in Saudis with Sickle Cell Disease

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.

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.

Precision Editing as a Therapeutic Approach for β-Hemoglobinopathies

Beta-hemoglobinopathies are the most common genetic disorders worldwide, caused by a wide spectrum of mutations in the β-globin locus, and associated with morbidity and early mortality in case of patient non-adherence to supportive treatment. Allogeneic transplantation of hematopoietic stem cells (allo-HSCT) used to be the only curative option, although the indispensable need for an HLA-matched donor markedly restricted its universal application. The evolution of gene therapy approaches made possible the ex vivo delivery of a therapeutic β- or γ- globin gene into patient-derived hematopoietic stem cells followed by the transplantation of corrected cells into myeloablated patients, having led to high rates of transfusion independence (thalassemia) or complete resolution of painful crises (sickle cell disease-SCD). Hereditary persistence of fetal hemoglobin (HPFH), a syndrome characterized by increased γ-globin levels, when co-inherited with β-thalassemia or SCD, converts hemoglobinopathies to a benign condition with mild clinical phenotype. The rapid development of precise genome editing tools (ZFN, TALENs, CRISPR/Cas9) over the last decade has allowed the targeted introduction of mutations, resulting in disease-modifying outcomes. In this context, genome editing tools have successfully been used for the introduction of HPFH-like mutations both in HBG1/HBG2 promoters or/and in the erythroid enhancer of BCL11A to increase HbF expression as an alternative curative approach for β-hemoglobinopathies. The current investigation of new HbF modulators, such as ZBTB7A, KLF-1, SOX6, and ZNF410, further expands the range of possible genome editing targets. Importantly, genome editing approaches have recently reached clinical translation in trials investigating HbF reactivation in both SCD and thalassemic patients. Showing promising outcomes, these approaches are yet to be confirmed in long-term follow-up studies.

Are Genetic Modifiers the Answer to Different Responses to Hydroxyurea Treatment?-A Pharmacogenetic Study in Sickle Cell Anemia Angolan Children

Sickle cell anemia (SCA) is an inherited disease affecting the hemoglobin that is particularly common in sub-Saharan Africa. Although monogenic, phenotypes are markedly heterogeneous in terms of severity and life span. Hydroxyurea is still the most common treatment for these patients, and the response to treatment is highly variable and seems to be an inherited trait. Therefore, identifying the variants that might predict hydroxyurea response is important for identifying patients who will have a poorer or non-response to treatment, and the ones that are more prone to suffer from severe side effects. In the present pharmacogenetic study, we analyzed the exons of 77 genes described in the literature as potentially associated with hydroxyurea metabolism in Angolan children treated with hydroxyurea and evaluated the drug response considering fetal hemoglobin levels, other hematological and biochemical parameters, hemolysis, number of vaso-occlusive crises and hospitalizations. Thirty variants were identified in 18 of those genes as possibly associated with drug response, five of them in gene DCHS2. Other polymorphisms in this gene were also associated with hematological, biochemical and clinical parameters. Further research examining the maximum tolerated dose and fixed dose with a larger sample size is necessary to corroborate these findings.

l-glutamine, crizanlizumab, voxelotor, and cell-based therapy for adult sickle cell disease: Hype or hope?

For more than two decades, hydroxyurea was the only therapeutic agent approved by the Food and Drug Administration (FDA) for sickle cell disease (SCD). Although curative allogeneic hematopoietic stem cell transplants (allo-HSCT) were also available, only very few patients underwent the procedure due to lack of matched-related donors. However, therapeutic options for SCD patients increased dramatically in the last few years. Three new agents, l-glutamine, crizanlizumab, and voxelotor, were approved by the FDA for use in SCD patients. The number of SCD patients who underwent allo-HSCT also increased as a result of advances in the prevention of graft failure and graft-versus-host disease from using mismatched donor HSC. More recently gene therapy was made available on clinical trials. The increased treatment options for SCD have led to a sense of optimism and excitement among many physicians that these new approaches would alter the clinical course and disease burden. Although these newer agents do provide hope to SCD patients, the hyped-up responses need to be evaluated in the context of reality. In this review, we will discuss and compare these new agents and cell-based therapy, evaluate their clinical and economic impacts, and examine their roles in reducing the disease burden.

Unique Polymorphisms at BCL11A, HBS1L-MYB and HBB Loci Associated with HbF in Kuwaiti Patients with Sickle Cell Disease

Patients with sickle cell disease (SCD) in Kuwait have elevated HbF levels ranging from ~10-44%; however, the modulating factors are unclear. We investigated the association of single nucleotide polymorphisms (SNPs) at BCL11A, HBS1L-MYB and HBB with HbF levels in 237 Kuwaiti SCD patients, divided into 3 subgroups according to their HbF levels. Illumina Ampliseq custom DNA panel was used for genotyping and confirmed by arrayed primer extension or Sanger sequencing. In the BCL11A locus, the CC genotype of rs7606173 [χ² = 16.5] and (GG) of rs10195871 [χ² = 15.0] were associated with Hb-F1 and HbF-2 subgroups, unlike rs1427404-T [χ² = 17.3], which showed the highest association across the three subgroups. HBS1L-MYB locus revealed 2 previously-described SNPs (rs66650371 [χ² = 9.5] and rs35795442 [χ² = 9.2]) and 2 previously-unreported SNPs, (rs13220662 [χ² = 6.2] and rs1406811 [χ² = 6.7]) that were associated with the HbF-3 subgroup, making this the key locus elevating HbF to the highest levels. HBB cluster variants were associated with lower levels of HbF (β = -1.1). We report four previously-unpublished variants showing significant association with HbF. Each of the three quantitative trait loci affects HbF levels differently; unique SNPs, especially in HBS1L-MYB, elevate HbF to the highest levels.

Haemoglobin switching modulator SNPs rs5006884 is associated with increased HbA2 in β-thalassaemia carriers

INTRODUCTION

Haemoglobin A₂ (HbA₂), the tetramer of α- and δ-globin chains, is used as a diagnostic biomarker for β-thalassaemia carriers. The HbA₂ levels are regulated by the presence of HPFH, δ-thalassaemia, HbA_1/2 gene triplication, and variants of KLF1, β-globin gene, and HbF regulating QTLs. Saudi Arabia has a high incidence of borderline HbA₂ levels, thereby making it difficult to classify the haemoglobinopathies. This study aims to investigate the association of known HbF enhancer QTL gene SNPs with HbA₂ levels.

MATERIAL AND METHODS

14 Specific SNPs in BCL11A, HMIP, OR51B6, HBBP1, and HBG2 loci were genotyped in 164 Saudi β-thalassaemia carriers by TaqMan assay to validate their role as regulators of HbA₂ levels. HbA2 levels were determined using the Variant II β-Thalassemia Short Program Recorder kit. The non-random association of these SNPs was tested using HaploView software. Protein interaction was assessed using 3D structure modelling for OR51B6 (rs5006884), comparative energy minimisation, and root-mean-square deviation (RMSD) prediction.

RESULTS

Elevated HbA₂ levels were associated with SNPs in HBBP1, OR51B6, and TCT haplotype from HBG2 promoter region. The bioinformatics modelling and prediction revealed that the exonic rs5006884 had RMSD value deviations and significantly varied binding energy minimisation. α-globin variations were found in 57.92% of individuals but were not associated with elevated HbA₂.

CONCLUSIONS

The haemoglobin switching modulators rs2071348, rs7482144, and rs5006884 are involved in regulation of HbA₂ level with rs5006884 influencing the tetramer formation. Screening for haemoglobinopathies should take these SNPs into consideration, specifically in borderline HbA₂ cases. Assiduous analysis of rs5006884 as HbA₂ modulator for amelioration of disease severity is recommended.

F cell numbers are associated with an X-linked genetic polymorphism and correlate with haematological parameters in patients with sickle cell disease

Patients with sickle cell disease (SCD) with high fetal haemoglobin (HbF) tend to have a lower incidence of complications and longer survival due to inhibition of deoxyhaemoglobin S (HbS) polymerisation by HbF. HbF-containing cells, namely F cells, are strongly influenced by genetic factors. We measured the percentage of F cells (Fcells%) in 222 patients with SCD to evaluate the association of (i) Fcells% with genetic HbF-modifier variants and (ii) Fcells% with haematological parameters. There was a different distribution of Fcells% in females compared to males. The association of the B-cell lymphoma/leukaemia 11A (BCL11A) locus with Fcells% (β = 8·238; P < 0·001) and with HbF% (β = 2·490; P < 0·001) was significant. All red cell parameters except for Hb and mean corpuscular Hb concentration levels in males and females were significantly different. The Fcells% was positively associated with mean cell Hb, mean cell volume and reticulocytes. To explain the significant gender difference in Fcells%, we tested for associations with single nucleotide polymorphisms on the X chromosomal region Xp22.2, where a genetic determinant of HbF had been previously hypothesised. We found in males a significant association with a SNP in FERM and PDZ domain-containing protein 4 (FRMPD4) and adjacent to male-specific lethal complex subunit 3 (MSL3). Thus, we have identified an X-linked locus that could account for a significant fraction of the Fcells% variation in patients with SCD.

Fetal Hemoglobin in Sickle Hemoglobinopathies: High HbF Genotypes and Phenotypes

Fetal hemoglobin (HbF) usually consists of 4 to 10% of total hemoglobin in adults of African descent with sickle cell anemia. Rarely, their HbF levels reach more than 30%. High HbF levels are sometimes a result of β-globin gene deletions or point mutations in the promoters of the HbF genes. Collectively, the phenotype caused by these mutations is called hereditary persistence of fetal hemoglobin, or HPFH. The pancellularity of HbF associated with these mutations inhibits sickle hemoglobin polymerization in most sickle erythrocytes so that these patients usually have inconsequential hemolysis and few, if any, vasoocclusive complications. Unusually high HbF can also be associated with variants of the major repressors of the HbF genes, BCL11A and MYB. Perhaps most often, we lack an explanation for very high HbF levels in sickle cell anemia.

Association of Exon 14 of the SOX6 Gene Sequence Variations with Response to Hydroxyurea Therapy in Patients Carrying Non Transfusion-Dependent Thalassemia

Hydroxyurea (HU) activates the γ-globin gene, resulting in increased Hb F synthesis. The SOX6 gene is a member of the Sox (Sry-type HMG box) family of transcription factors, characterized by minor groove binding domain. The DNA binding domain of this gene is encoded by exon 14. We assessed the relationship between response to HU and exon 14 of the SOX6 gene sequence variations in patients with non transfusion-dependent thalassemia (NTDT). One hundred NTDT patients from southern Iran underwent HU therapy randomly participated in this cross-sectional study between February 2013 and October 2014. Based on response to HU therapy, the patients were divided into two groups: good and poor responder. Sequence variations of exon 14 of the SOX6 gene was assayed by the Sanger sequencing technique. From all evaluated single nucleotide polymorphisms (SNPs) as above, we found no significant association between sequence variations of exon 14 of the SOX6 gene and response to HU therapy (p > 0.05). It seems that no SNPs in exon 14 of the SOX6 gene is associated with response to HU in NTDT patients, but more studies are needed for further evaluation.

Potential causal role of l-glutamine in sickle cell disease painful crises: A Mendelian randomization analysis

In a recent clinical trial, the metabolite l-glutamine was shown to reduce painful crises in sickle cell disease (SCD) patients. To support this observation and identify other metabolites implicated in SCD clinical heterogeneity, we profiled 129 metabolites in the plasma of 705 SCD patients. We tested correlations between metabolite levels and six SCD-related complications (painful crises, cholecystectomy, retinopathy, leg ulcer, priapism, aseptic necrosis) or estimated glomerular filtration rate (eGFR), and used Mendelian randomization (MR) to assess causality. We found a potential causal relationship between l-glutamine levels and painful crises (N = 1278, odds ratio (OR) [95% confidence interval] = 0.68 [0.52-0.89], P = 0.0048). In two smaller SCD cohorts (N = 299 and 406), the protective effect of l-glutamine was observed (OR = 0.82 [0.50-1.34]), although the MR result was not significant (P = 0.44). We identified 66 significant correlations between the levels of other metabolites and SCD-related complications or eGFR. We tested these correlations for causality using MR analyses and found no significant causal relationship. The baseline levels of quinolinic acid were associated with prospectively ascertained survival in SCD patients, and this effect was dependent on eGFR. Metabolomics provide a promising approach to prioritize small molecules that may serve as biomarkers or drug targets in SCD.

Association between BCL11A, HSB1L-MYB, and XmnI γG-158 (C/T) gene polymorphism and hemoglobin F level in Egyptian sickle cell disease patients

Sickle cell disease (SCD) is a monogenic disease characterized by multisystem morbidity and highly variable clinical course. Inter-individual variability in hemoglobin F (HbF) levels is one of the main modifiers that account for the clinical heterogeneity in SCD. HbF levels are affected by, among other factors, single nucleotide polymorphisms (SNPs) at the BCL11A gene and the HBS1L-MYB intergenic region and Xmn1 gene. Our aim was to investigate HbF-enhancer haplotypes at these loci to obtain a first overview of the genetic situation of SCD patients in Egypt and its impact on the severity of the disease. The study included 100 SCD patients and 100 matched controls. Genotyping of BCL11A (rs1886868 C/T), HBS1L-MYB (rs9389268 A/G) and Xmn1 γ^G158 (rs7842144 C/T) SNPs showed no statistically significant difference between SCD patients and controls except for the hetero-mutant genotypes of BCL11A which was significantly higher in SCD patients compared with controls. Baseline HbF levels were significantly higher in those with co-inheritance of polymorphic genotypes of BCL11A + HSB1L-MYB and BCL11A + Xmn1. Steady-state HbF levels, used as an indicator of disease severity, were significantly higher in SCD-Sβ patients having the polymorphic genotypes of HSB1L-MYB. Fold change of HbF in both patient groups did not differ between those harboring the wild and the polymorphic genotypes of the studied SNPs. In conclusion, BCL11A, HSB1L, and Xmn1 genetic polymorphisms had no positive impact on baseline HbF levels solely but had if coexisted. Discovery of the molecular mechanisms controlling HbF production could provide a more effective strategy for HbF induction.

The Effect of Five Single Nucleotide Polymorphisms on Hb F Variation of β-Thalassemia Traits and Hematologically Normal Individuals in Southeast Turkey

β-Thalassemia (β-thal) is caused by deficiency of β-globin chain synthesis and leads to the accumulation of unstable globin chain production. This results in a higher Hb F level in order to neutralize the excess α chains. In addition, γ-globin gene expression, due to genetic factors after birth, leads to increased Hb F levels in adulthood [hereditary persistence of fetal hemoglobin (Hb) (HPFH)]. In this study, the relationship between β-thal trait and individuals with suspected HPFH and a control group was investigated in Adıyaman, Turkey. Single nucleotide polymorphism (SNP) analyses were performed in five different polymorphic regions using real-time polymerase chain reaction (qPCR) methods [rs4671393 (G>A), rs766432 (A>C), rs9402686 (G>A), rs28384513 (T>G), rs1609812 (A>G)]. No significant difference was found between the control and β-thal group in the codominant inheritance model in the rs1609812 (A>G) polymorphism region only, while all the other polymorphic regions were found to be statistically significant. It was found that different genotype models increased Hb F levels between 1.6- and 3.06-fold in four studied polymorphic regions [rs4671393 (G>A), rs766432 (A>C), rs9402686 (G>A), rs28384513 (T>G)]. All of the polymorphic regions increased the Hb F levels from 1.86- to 24.76-fold, except rs9402686 (G>A) and rs28384513 (T>G) over dominant and rs1609812 (A>G) codominant inheritance models. The AC and AA genotypes increased Hb F levels in the B-cell CLL/lymphoma 11 A haplotype studies. It was determined that both haplotypes 2 and 4 increased Hb F levels. As a result, SNPs strongly affect the Hb F levels in both healthy individuals and β-thal trait.

Functional polymorphisms of BCL11A and HBS1L-MYB genes affect both fetal hemoglobin level and clinical outcomes in a cohort of children with sickle cell anemia

Fetal hemoglobin (HbF) ameliorates clinical severity of sickle cell anemia (SCA). The major loci regulating HbF levels are HBB cluster, BCL11A, and HMIP-2 (HBS1L-MYB). However, the impact of noncoding single-nucleotide polymorphisms (SNPs) in these loci on clinical outcomes and their functional role on regulating HbF levels should be better elucidated. Therefore, we performed comprehensive association analyses of 14 noncoding SNPs in five loci with HbF levels and with clinical outcomes in a cohort of 250 children with SCA from Southeastern Brazil, and further performed functional annotation of these SNPs. We found SNPs independently associated with HbF levels: rs4671393 in BCL11A (β-coefficient = 0.28), rs9399137 in HMIP-2A (β-coefficient = 0.16), and rs4895441 in HMIP-2B (β-coefficient = 0.15). Patients carrying minor (HbF-boosting) alleles for rs1427407, rs93979137, rs4895441, rs9402686, and rs9494145 showed reduced count of reticulocytes (p < 0.01), while those carrying the T allele of rs9494145 showed lower white blood cell count (p = 0.002). Carriers of the minor allele for rs9402686 showed higher peripheral saturation of oxygen (p = 0.002). Patients carrying minor alleles in BCL11A showed lower risk of transfusion incidence rate ratio (IRR ≥ 1.3; p < 0.0001). This effect was independent of HbF effect (p = 0.005). Carriers of minor alleles for rs9399137 and rs9402686 showed lower risk of acute chest syndrome (IRR > 1.3; p ≤ 0.01). Carriers of the reference allele for rs4671393 showed lower risk of infections (IRR = 1.16; p = 0.01). In conclusion, patients carrying HbF-boosting alleles of BCL11A and HMIP-2 were associated with milder clinical phenotypes. Higher HbF concentration may underlie this effect.

Genetic variants in the G gamma-globin promoter modulate fetal hemoglobin expression in the Colombian population

Fetal hemoglobin (HbF) is a determining factor for the development of sickle cell anemia. High HbF levels lower the intensity of symptoms of this disease. HbF levels can vary in patients with sickle cell anemia and individuals without the disease. The purpose of this study was to identify the genetic variants in the G gamma-globin gene promoter that can modulate HbF expression in patients with sickle cell anemia and healthy individuals from Colombia. In total, 413 bp of the G gamma-globin gene promoter were sequenced in 60 patients with sickle cell anemia and 113 healthy individuals. The allelic and genotype frequencies of the identified variants were compared between individuals with low and high HbF for both patients and healthy individuals. In total, we identified 15 variants in both groups, only three of which were shared between patients and healthy individuals. In healthy individuals, sites -16 and -309 (rs112479156) exhibited differences in allele frequencies. The mutant allele of -16 lowered the production of HbF, whereas the mutant allele of -309 increased its production. These results reveal the presence of different mechanisms of HbF regulation between patients with sickle cell and healthy individuals.

STEPS: an efficient prospective likelihood approach to genetic association analyses of secondary traits in extreme phenotype sequencing

It has been well acknowledged that methods for secondary trait (ST) association analyses under a case-control design (ST$_{\text{CC}}$) should carefully consider the sampling process to avoid biased risk estimates. A similar situation also exists in the extreme phenotype sequencing (EPS) designs, which is to select subjects with extreme values of continuous primary phenotype for sequencing. EPS designs are commonly used in modern epidemiological and clinical studies such as the well-known National Heart, Lung, and Blood Institute Exome Sequencing Project. Although naïve generalized regression or ST$_{\text{CC}}$ method could be applied, their validity is questionable due to difference in statistical designs. Herein, we propose a general prospective likelihood framework to perform association testing for binary and continuous STs under EPS designs (STEPS), which can also incorporate covariates and interaction terms. We provide a computationally efficient and robust algorithm to obtain the maximum likelihood estimates. We also present two empirical mathematical formulas for power/sample size calculations to facilitate planning of binary/continuous STs association analyses under EPS designs. Extensive simulations and application to a genome-wide association study of benign ethnic neutropenia under an EPS design demonstrate the superiority of STEPS over all its alternatives above.

Sickle cell disease nephropathy: an update on risk factors and potential biomarkers in pediatric patients

One of the major chronic complications of sickle cell disease (SCD) is sickle cell nephropathy. The aim of this review is to discuss the pathophysiology, natural history, clinical manifestations, risk factors, biomarkers and therapeutic approaches for sickle cell nephropathy, focusing on studies with pediatric patients. The earliest manifestation of renal disease is an increase in the glomerular filtration rate. A finding that may also be observed in early childhood is microalbuminuria. Nephrin, KIM-1, VGFs, chemokines and renin-angiotensin system molecules have emerged as potential early markers of renal dysfunction in SCD. In regards to a therapeutic approach, renin-angiotensin system inhibitors and angiotensin receptor blockers seem to be effective for the control of albuminuria in adults with SCD, although new studies in children are needed. The precise moment to begin renoprotection in SCD patients who should be treated remains to be determined.

Genetic comparison of sickle cell anaemia cohorts from Brazil and the United States reveals high levels of divergence

Genetic analysis of admixed populations raises special concerns with regard to study design and data processing, particularly to avoid population stratification biases. The point mutation responsible for sickle cell anaemia codes for a variant hemoglobin, sickle hemoglobin or HbS, whose presence drives the pathophysiology of disease. Here we propose to explore ancestry and population structure in a genome-wide study with particular emphasis on chromosome 11 in two SCA admixed cohorts obtained from urban populations of Brazil (Pernambuco and São Paulo) and the United States (Pennsylvania). Ancestry inference showed different proportions of European, African and American backgrounds in the composition of our samples. Brazilians were more admixed, had a lower African background (43% vs. 78% on the genomic level and 44% vs. 76% on chromosome 11) and presented a signature of positive selection and Iberian introgression in the HbS region, driving a high differentiation of this locus between the two cohorts. The genetic structures of the SCA cohorts from Brazil and US differ considerably on the genome-wide, chromosome 11 and HbS mutation locus levels.

Metastatic group 3 medulloblastoma is driven by PRUNE1 targeting NME1-TGF-β-OTX2-SNAIL via PTEN inhibition

Genetic modifications during development of paediatric groups 3 and 4 medulloblastoma are responsible for their highly metastatic properties and poor patient survival rates. PRUNE1 is highly expressed in metastatic medulloblastoma group 3, which is characterized by TGF-β signalling activation, c-MYC amplification, and OTX2 expression. We describe the process of activation of the PRUNE1 signalling pathway that includes its binding to NME1, TGF-β activation, OTX2 upregulation, SNAIL (SNAI1) upregulation, and PTEN inhibition. The newly identified small molecule pyrimido-pyrimidine derivative AA7.1 enhances PRUNE1 degradation, inhibits this activation network, and augments PTEN expression. Both AA7.1 and a competitive permeable peptide that impairs PRUNE1/NME1 complex formation, impair tumour growth and metastatic dissemination in orthotopic xenograft models with a metastatic medulloblastoma group 3 cell line (D425-Med cells). Using whole exome sequencing technology in metastatic medulloblastoma primary tumour cells, we also define 23 common 'non-synonymous homozygous' deleterious gene variants as part of the protein molecular network of relevance for metastatic processes. This PRUNE1/TGF-β/OTX2/PTEN axis, together with the medulloblastoma-driver mutations, is of relevance for future rational and targeted therapies for metastatic medulloblastoma group 3.10.1093/brain/awy039_video1awy039media15742053534001.

Perspective: A Novel Prognostic for Sickle Cell Disease

Sickle hemoglobin (α₂β^S ₂) polymerization drives disease pathophysiology in sickle cell anemia. Fetal hemoglobin (α₂γ₂) restricts disease severity by inhibiting the polymerization of sickle hemoglobin in a concentration-dependent manner. Clinical decision-making relies on diagnostic technologies evaluating fetal hemoglobin as mean percent or mean quantity in blood. Limitation of this approach is exemplified by patients with significant high fetal hemoglobin levels and severe disease, suggesting that fetal hemoglobin is unevenly distributed across F-cells. Therefore, determination of fetal hemoglobin/F-cell would provide a new paradigm for ascertaining prognosis and response to fetal hemoglobin-inducing agents. Measurement of fetal hemoglobin/F-cell, ultimately adapted to widespread standardized analytical use, is a promising fetal hemoglobin-related prognostic approach to monitor the severity of sickle cell disease and the best "phenotype" to follow when developing new candidate fetal hemoglobin inducers or titrating hydroxyurea in treated sickle cell patients.

Metformin induces FOXO3-dependent fetal hemoglobin production in human primary erythroid cells

Induction of red blood cell (RBC) fetal hemoglobin (HbF; α2γ2) ameliorates the pathophysiology of sickle cell disease (SCD) by reducing the concentration of sickle hemoglobin (HbS; α2β^S2) to inhibit its polymerization. Hydroxyurea (HU), the only US Food and Drug Administration (FDA)-approved drug for SCD, acts in part by inducing HbF; however, it is not fully effective, reflecting the need for new therapies. Whole-exome sequence analysis of rare genetic variants in SCD patients identified FOXO3 as a candidate regulator of RBC HbF. We validated these genomic findings through loss- and gain-of-function studies in normal human CD34⁺ hematopoietic stem and progenitor cells induced to undergo erythroid differentiation. FOXO3 gene silencing reduced γ-globin RNA levels and HbF levels in erythroblasts, whereas overexpression of FOXO3 produced the opposite effect. Moreover, treatment of primary CD34⁺ cell-derived erythroid cultures with metformin, an FDA-approved drug known to enhance FOXO3 activity in nonerythroid cells, caused dose-related FOXO3-dependent increases in the percentage of HbF protein and the fraction of HbF-immunostaining cells (F cells). Combined HU and metformin treatment induced HbF additively and reversed the arrest in erythroid maturation caused by HU treatment alone. HbF induction by metformin in erythroid precursors was dependent on FOXO3 expression and did not alter expression of BCL11A, MYB, or KLF1. Collectively, our data implicate FOXO3 as a positive regulator of γ-globin expression and identify metformin as a potential therapeutic agent for SCD.

Statistical selection of biological models for genome-wide association analyses

Genome-wide association studies have discovered many biologically important associations of genes with phenotypes. Typically, genome-wide association analyses formally test the association of each genetic feature (SNP, CNV, etc) with the phenotype of interest and summarize the results with multiplicity-adjusted p-values. However, very small p-values only provide evidence against the null hypothesis of no association without indicating which biological model best explains the observed data. Correctly identifying a specific biological model may improve the scientific interpretation and can be used to more effectively select and design a follow-up validation study. Thus, statistical methodology to identify the correct biological model for a particular genotype-phenotype association can be very useful to investigators. Here, we propose a general statistical method to summarize how accurately each of five biological models (null, additive, dominant, recessive, co-dominant) represents the data observed for each variant in a GWAS study. We show that the new method stringently controls the false discovery rate and asymptotically selects the correct biological model. Simulations of two-stage discovery-validation studies show that the new method has these properties and that its validation power is similar to or exceeds that of simple methods that use the same statistical model for all SNPs. Example analyses of three data sets also highlight these advantages of the new method. An R package is freely available at www.stjuderesearch.org/site/depts/biostats/maew.

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.

Common α-globin variants modify hematologic and other clinical phenotypes in sickle cell trait and disease

Co-inheritance of α-thalassemia has a significant protective effect on the severity of complications of sickle cell disease (SCD), including stroke. However, little information exists on the association and interactions for the common African ancestral α-thalassemia mutation (-α3.7 deletion) and β-globin traits (HbS trait [SCT] and HbC trait) on important clinical phenotypes such as red blood cell parameters, anemia, and chronic kidney disease (CKD). In a community-based cohort of 2,916 African Americans from the Jackson Heart Study, we confirmed the expected associations between SCT, HbC trait, and the -α3.7 deletion with lower mean corpuscular volume/mean corpuscular hemoglobin and higher red blood cell count and red cell distribution width. In addition to the recently recognized association of SCT with lower estimated glomerular filtration rate and glycated hemoglobin (HbA1c), we observed a novel association of the -α3.7 deletion with higher HbA1c levels. Co-inheritance of each additional copy of the -α3.7 deletion significantly lowered the risk of anemia and chronic kidney disease among individuals with SCT (P-interaction = 0.031 and 0.019, respectively). Furthermore, co-inheritance of a novel α-globin regulatory variant was associated with normalization of red cell parameters in individuals with the -α3.7 deletion and significantly negated the protective effect of α-thalassemia on stroke in 1,139 patients with sickle cell anemia from the Cooperative Study of Sickle Cell Disease (CSSCD) (P-interaction = 0.0049). Functional assays determined that rs11865131, located in the major alpha-globin enhancer MCS-R2, was the most likely causal variant. These findings suggest that common α- and β-globin variants interact to influence hematologic and clinical phenotypes in African Americans, with potential implications for risk-stratification and counseling of individuals with SCD and SCT.

Whole exome sequencing reveals intertumor heterogeneity and distinct genetic origins of sporadic synchronous colorectal cancer

Sporadic synchronous colorectal cancer (CRC) refers to more than one primary tumor detected in a single patient at the time of the first diagnosis without predisposition of cancer development. Given the same genetic and microenvironment they raise, sporadic synchronous CRC is a unique model to study CRC tumorigenesis. We performed whole exome sequencing in 32 fresh frozen tumor lesions from 15 patients with sporadic synchronous CRC to compare their genetic alterations. This approach identified ubiquitously mutated genes in the range from 0.34% to 4.22% and from 0.8% to 7.0% in non-hypermutated tumors and hypermutated tumors, respectively, in a single patient. We show that both ubiquitously mutated genes and candidate cancer genes from different tumors in the same patient mutated at different sites. Consistently, obvious differences in somatic copy number variations (SCNV) were found in most patients with non-hypermutated tumor lesions, which had ubiquitous copy number amplification rates ranging from 0% to 8.8% and ubiquitous copy number deletion rates ranging from 0% to 8.2%. Hypermutated lesions were nearly diploid with 0% to 18.8% common copy number aberrations. Accordingly, clonal structures, altered signaling pathways and druggable genes in a single patient with synchronous CRC varied significantly. Taken together, the disparate SCNVs and mutations in synchronous CRC supported the field effect theory of tumorigenesis. Moreover, the intertumor heterogeneity of synchronous CRCs implies that analysis of all tumor lesions from the same patient is necessary for appropriate clinical treatment decisions.

Fetal hemoglobin in sickle cell anemia: The Arab-Indian haplotype and new therapeutic agents

Fetal hemoglobin (HbF) has well-known tempering effects on the symptoms of sickle cell disease and its levels vary among patients with different haplotypes of the sickle hemoglobin gene. Compared with sickle cell anemia haplotypes found in patients of African descent, HbF levels in Saudi and Indian patients with the Arab-Indian (AI) haplotype exceed that in any other haplotype by nearly twofold. Genetic association studies have identified some loci associated with high HbF in the AI haplotype but these observations require functional confirmation. Saudi patients with the Benin haplotype have HbF levels almost twice as high as African patients with this haplotype but this difference is unexplained. Hydroxyurea is still the only FDA approved drug for HbF induction in sickle cell disease. While most patients treated with hydroxyurea have an increase in HbF and some clinical improvement, 10 to 20% of adults show little response to this agent. We review the genetic basis of HbF regulation focusing on sickle cell anemia in Saudi Arabia and discuss new drugs that can induce increased levels of HbF.

A phased SNP-based classification of sickle cell anemia HBB haplotypes

Background

Sickle cell anemia causes severe complications and premature death. Five common β-globin gene cluster haplotypes are each associated with characteristic fetal hemoglobin (HbF) levels. As HbF is the major modulator of disease severity, classifying patients according to haplotype is useful. The first method of haplotype classification used restriction fragment length polymorphisms (RFLPs) to detect single nucleotide polymorphisms (SNPs) in the β-globin gene cluster. This is labor intensive, and error prone.

Methods

We used genome-wide SNP data imputed to the 1000 Genomes reference panel to obtain phased data distinguishing parental alleles.

Results

We successfully haplotyped 813 sickle cell anemia patients previously classified by RFLPs with a concordance >98%. Four SNPs (rs3834466, rs28440105, rs10128556, and rs968857) marking four different restriction enzyme sites unequivocally defined most haplotypes. We were able to assign a haplotype to 86% of samples that were either partially or misclassified using RFLPs.

Conclusion

Phased data using only four SNPs allowed unequivocal assignment of a haplotype that was not always possible using a larger number of RFLPs. Given the availability of genome-wide SNP data, our method is rapid and does not require high computational resources.

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.

Variability of homozygous sickle cell disease: The role of alpha and beta globin chain variation and other factors

The single base molecular substitution characterizing sickle cell haemoglobin, β⁶glu→val, might be expected to result in predictable haematological and clinical features. However, the disease manifests remarkable diversity believed to reflect the interaction with other genetic and environmental factors. Some of the genetic modifiers include the beta globin haplotypes, alpha thalassaemia, factors influencing the persistence of fetal haemoglobin and the effects of the environment are addressed in this review. It is concluded that much of the genetic data present conflicting results. Environmental factors such as climate and infections, and psychological, educational and social support mechanisms also influence expression of the disease. These interactions illustrate how the expression of a 'single gene' disorder may be influenced by a variety of other genetic and environmental factors.

Genome-wide association study of erythrocyte density in sickle cell disease patients

Deoxy-hemoglobin S polymerization into rigid fibers is the direct cause of the clinical sequelae observed in sickle cell disease (SCD). The rate of polymerization of sickle hemoglobin is determined primarily by intracellular hemoglobin concentration, itself dependent on the amount of sickle hemoglobin and on red blood cell (RBC) volume. Dense, dehydrated RBC (DRBC) are observed in SCD patients, and their number correlates with hemolytic parameters and complications such as renal dysfunction, leg ulcers and priapism. To identify new genes involved in RBC hydration in SCD, we performed the first genome-wide association study for DRBC in 374 sickle cell anemia (HbSS) patients. We did not find genome-wide significant results, indicating that variants that modulate DRBC have modest-to-weak effects. A secondary analysis demonstrated a nominal association (P=0.003) between DRBC in SCD patients and a variant associated with mean corpuscular hemoglobin concentration (MCHC) in non-anemic individuals. This intronic variant controls the expression of ATP2B4, the main calcium pump in erythrocytes. Our study highlights ATP2B4 as a promising target for modulation of RBC hydration in SCD patients.

The diversified function and potential therapy of ectopic olfactory receptors in non-olfactory tissues

Olfactory receptors (ORs) are mainly distributed in olfactory neurons and play a key role in detecting volatile odorants, eventually resulting in the production of smell perception. Recently, it is also reported that ORs are expressed in non-olfactory tissues including heart, lung, sperm, skin, and cancerous tissues. Interestingly, ectopic ORs are associated with the development of diseases in non-olfactory tissues. For instance, ectopic ORs initiate the hypoxic ventilatory responses and maintain the oxygen homeostasis of breathing in the carotid body when oxygen levels decline. Ectopic ORs induce glucose homeostasis in diabetes. Ectopic ORs regulate systemic blood pressure by increasing renin secretion and vasodilation. Ectopic ORs participate in the process of tumor cell proliferation, apoptosis, metastasis, and invasiveness. Ectopic ORs accelerate the occurrence of obesity, angiogenesis and wound-healing processes. Ectopic ORs affect fetal hemoglobin levels in sickle cell anemia and thalassemia. Finally, we also elaborate some ligands targeting for ORs. Obviously, the diversified function and related signal pathway of ectopic ORs may play a potential therapeutic target in non-olfactory tissues. Thus, this review focuses on the latest research results about the diversified function and therapeutic potential of ectopic ORs in non-olfactory tissues.

Relationship Between Some Single-nucleotide Polymorphism and Response to Hydroxyurea Therapy in Iranian Patients With β-Thalassemia Intermedia

OBJECTIVE

To evaluate the possible relationship between hydroxyurea (HU) response and some single-nucleotide polymorphism (SNP) in patients affected by β-thalassemia intermedia.

MATERIALS AND METHODS

In this cross-sectional study, 100 β-thalassemia intermedia patients who were taking HU with a dose of 8 to 15 mg/kg body weight per day for a period of at least 6 months were randomly selected between February 2013 and October 2014 in southern Iran. HU response was defined based on decrease or cessation of the blood transfusion need and evaluation of Hb level.

RESULTS

In univariate analysis, from all evaluated SNPs, only rs10837814 SNP of olfactory receptors (ORs) OR51B2 showed a significant association with HU response (P=0.038) and from laboratory characteristics, only nucleated red blood cells showed significant associations (116%±183%) in good responders versus (264%±286%) in poor responders (P=0.045). In multiple logistic regression, neither laboratory variables nor different SNPs, showed significant association with HU response. Three novel nucleotide variations (-665 [A→C], -1301 [T→G],-1199 delA) in OR51B2 gene were found in good responders.

CONCLUSIONS

None of the evaluated SNPs in our study showed significant association with HU response. Further larger studies and evaluation of other genes are suggested.

Genetic Variants at BCL11A and HBS1L-MYB loci Influence Hb F Levels in Chinese Zhuang β-Thalassemia Intermedia Patients

Increased Hb F levels can ameliorate the symptoms of β-thalassemia (β-thal). Due to the genetic heterogenicity of β-thal, the relationship between genetic variants in modifier genes and Hb F level has been studied in different populations. The Chinese Zhuang has the second largest population in China and has 6.78% prevalence of β-thal. However, the effects of these single nucleotide polymorphism (SNP) variants on the Hb F levels of β-thal intermedia (β-TI) patients in this population have not been reported. To explore the association between modifier loci (β-globin gene cluster, HBS1L-MYB intergenic region and BCL11A) and Hb F levels in Chinese Zhuang β-TI patients, 96 unrelated β-TI patients (50 males and 46 females) with different Hb F levels were recruited and genotyped by mass spectrometry. A total of 13 SNPs were confirmed to be in a significant relationship with Hb F levels in this population. Of these, high-risk genotypes of six Hb F-associated SNPs, rs9376090, rs7776054, rs9399137, rs9389268, rs9402685 in the HBS1L-MYB intergenic region and rs189984760 in the BCL11A locus, showed association with high Hb F levels, especially for SNPs in linkage disequilibrium. One novel Hb F-associated SNP, rs189984760, was identified in our study. Our findings will be of valuable reference for correlation between modifier genes and Hb F in Chinese Zhuang populations and may lead to better understand the modifying mechanisms for β-thal.

Protective BCL11A and HBS1L-MYB polymorphisms in a cohort of 102 Congolese patients suffering from sickle cell anemia

BACKGROUND

We aimed to investigate the distribution of selected BCL11A and HMIP polymorphisms (SNP's), and to assess the correlation with HPFH in a cohort of sickle cell patients.

METHODS

A preliminary cross-sectional study was conducted in 102 patients. Group 1 was composed of patients with HPFH and Group 2 consisted of patients without HbF. We assessed 8 SNPs previously associated with HPFH in cohorts genetically close to the Congolese population. Observed frequencies were compared to expected frequencies.

RESULTS

In the group 1, at rs7606173, the observed frequency for the genotype GG was significantly higher and the genotype GC was significantly lower than their respective expected frequencies. At rs9399137, the observed frequency of the genotype TT was significantly lower than expected. Conversely, the observed frequency of the genotype TC was significantly higher than expected. The observed frequency of the genotype TT at rs11886868 was significantly lower than the expected whereas the frequency of the genotype TC was significantly higher than observed. The lowest HbF level was recorded in patients with genotype CC at rs11886868.

CONCLUSION

In this preliminary study, the results demonstrate that alleles of some of the 8 studied SNPs are not randomly distributed among patients with or without HPFH in this cohort.

Variant-aware saturating mutagenesis using multiple Cas9 nucleases identifies regulatory elements at trait-associated loci

Cas9-mediated, high-throughput, saturating in situ mutagenesis permits fine-mapping of function across genomic segments. Disease- and trait-associated variants from genome-wide association studies largely cluster in regulatory DNA. Here we demonstrate the use of multiple designer nucleases and variant-aware library design to interrogate trait-associated regulatory DNA at high resolution. We developed a computational tool for the creation of saturating mutagenesis libraries with single or combinatorial nucleases with incorporation of variants. We applied this methodology to the HBS1L-MYB intergenic region, a locus associated with red blood cell traits, including fetal hemoglobin levels. This approach identified putative regulatory elements that control MYB expression. Analysis of genomic copy number highlighted potential false positive regions, which emphasizes the importance of off-target analysis in design of saturating mutagenesis experiments. Taken together, these data establish a widely applicable high-throughput and high-resolution methodology to reliably identify minimal functional sequences within large regions of disease- and trait-associated DNA.

Existence of HbF Enhancer Haplotypes at HBS1L-MYB Intergenic Region in Transfusion-Dependent Saudi β-Thalassemia Patients

Background and Objectives. β-Thalassemia and sickle cell disease are genetic disorders characterized by reduced and abnormal β-globin chain production, respectively. The elevation of fetal hemoglobin (HbF) can ameliorate the severity of these disorders. In sickle cell disease patients, the HbF level elevation is associated with three quantitative trait loci (QTLs), BCL11A, HBG2 promoter, and HBS1L-MYB intergenic region. This study elucidates the existence of the variants in these three QTLs to determine their association with HbF levels of transfusion-dependent Saudi β-thalassemia patients. Materials and Methods. A total of 174 transfusion-dependent β-thalassemia patients and 164 healthy controls from Eastern Province of Saudi Arabia were genotyped for fourteen single nucleotide polymorphisms (SNPs) from the three QTL regions using TaqMan assay on real-time PCR. Results. Genotype analysis revealed that six alleles of HBS1L-MYB QTL (rs9376090C p = 0.0009, rs9399137C p = 0.008, rs4895441G p = 0.004, rs9389269C p = 0.008, rs9402686A p = 0.008, and rs9494142C p = 0.002) were predominantly associated with β-thalassemia. In addition, haplotype analysis revealed that haplotypes of HBS1L-MYB (GCCGCAC p = 0.022) and HBG2 (GTT p = 0.009) were also predominantly associated with β-thalassemia. Furthermore, the HBS1L-MYB region also exhibited association with the high HbF cohort. Conclusion. The stimulation of HbF gene expression may provide alternative therapies for the amelioration of the disease severity of β-thalassemia.

GWAS study using DNA pooling strategy identifies association of variant rs4910623 in OR52B4 gene with anti-VEGF treatment response in age-related macular degeneration

Pooled DNA based GWAS to determine genetic association of SNPs with visual acuity (VA) outcome in anti-vascular endothelial growth factor (anti-VEGF) treated neovascular age-related macular degeneration (nAMD) patients. We performed pooled DNA based GWAS on 285 anti-VEGF treated nAMD patients using high density Illumina 4.3 M array. Primary outcome was change in VA in Early Treatment Diabetic Retinopathy Study (ETDRS) letters after 6 months of anti-VEGF treatment (patients who lost ≥5 ETDRS letters classified as non-responders and all remaining classified as responders). GWAS analysis identified 44 SNPs of interest: 37 with strong evidence of association (p < 9 × 10⁻⁸), 2 in drug resistance genes (p < 5 × 10⁻⁶) and 5 nonsynonymous changes (p < 1 × 10⁻⁴). In the validation phase, individual genotyping of 44 variants showed three SNPs (rs4910623 p = 5.6 × 10⁻⁵, rs323085 p = 6.5 × 10⁻⁴ and rs10198937 p = 1.30 × 10⁻³) remained associated with VA response at 6 months. SNP rs4910623 also associated with treatment response at 3 months (p = 1.5 × 10⁻³). Replication of these three SNPs in 376 patients revealed association of rs4910623 with poor VA response after 3 and 6 months of treatment (p = 2.4 × 10⁻³ and p = 3.5 × 10⁻², respectively). Meta-analysis of both cohorts (673 samples) confirmed association of rs4910623 with poor VA response after 3 months (p = 1.2 × 10⁻⁵) and 6 months (p = 9.3 × 10⁻⁶) of treatment in nAMD patients.

Studies of novel variants associated with Hb F in Sardinians and Tanzanians in sickle cell disease patients from Cameroon

High level of Hb F has been shown to improve survival in sickle cell disease. Among 453 Cameroonians with sickle cell disease, we have investigated 18 selected single-nucleotide polymorphisms (SNPs) in novel and suggestive loci associated with Hb F level identified through a genomewide association study in sickle cell disease patients in Tanzania, and whole-genome sequencing of a population from Sardinia. Seven of 10 variants reported in Sardinians were either monomorphic or very rare in the Cameroonians. No associations were observed with any SNPs and Hb F levels in Cameroonians affected by sickle cell disease. The present study illustrates the complexity of replicating Hb F-promoting variants association results across populations.

A candidate transacting modulator of fetal hemoglobin gene expression in the Arab-Indian haplotype of sickle cell anemia

Fetal hemoglobin (HbF) levels are higher in the Arab-Indian (AI) β-globin gene haplotype of sickle cell anemia compared with African-origin haplotypes. To study genetic elements that effect HbF expression in the AI haplotype we completed whole genome sequencing in 14 Saudi AI haplotype sickle hemoglobin homozygotes-seven selected for low HbF (8.2% ± 1.3%) and seven selected for high HbF (23.5% ± 2.6%). An intronic single nucleotide polymorphism (SNP) in ANTXR1, an anthrax toxin receptor (chromosome 2p13), was associated with HbF. These results were replicated in two independent Saudi AI haplotype cohorts of 120 and 139 patients, but not in 76 Saudi Benin haplotype, 894 African origin haplotype and 44 AI haplotype patients of Indian origin, suggesting that this association is effective only in the Saudi AI haplotype background. ANTXR1 variants explained 10% of the HbF variability compared with 8% for BCL11A. These two genes had independent, additive effects on HbF and together explained about 15% of HbF variability in Saudi AI sickle cell anemia patients. ANTXR1 was expressed at mRNA and protein levels in erythroid progenitors derived from induced pluripotent stem cells (iPSCs) and CD34⁺ cells. As CD34⁺ cells matured and their HbF decreased ANTXR1 expression increased; as iPSCs differentiated and their HbF increased, ANTXR1 expression decreased. Along with elements in cis to the HbF genes, ANTXR1 contributes to the variation in HbF in Saudi AI haplotype sickle cell anemia and is the first gene in trans to HBB that is associated with HbF only in carriers of the Saudi AI haplotype. Am. J. Hematol. 91:1118-1122, 2016. © 2016 Wiley Periodicals, Inc.

An Expert Review of Pharmacogenomics of Sickle Cell Disease Therapeutics: Not Yet Ready for Global Precision Medicine

Sickle cell disease (SCD) is a blood disease caused by a single nucleotide substitution (T > A) in the beta globin gene on chromosome 11. The single point mutation (Glu6Val) promotes polymerization of hemoglobin S (HbS) and causes sickling of erythrocytes. Vaso-occlusive painful crises are associated with recurrent and long-term use of analgesics/opioids and hydroxyurea (HU) by people living with SCD. The present analysis offers a state-of-the-art expert review of the effectiveness of pharmacogenomics/genetics of pain management in SCD, with specific focus on HU and opioids. The literature search used the following keywords: SCD, pharmacogenomics, pharmacogenetics, pain, antalgics, opioids, morphine, and HU. The literature was scanned until March 2016, with specific inclusion of targeted landmark and background articles on SCD. Surprisingly, our review identified only a limited number of studies that addressed the genetic/genomic basis of variable responses to pain (e.g., variants in OPRM1, HMOX-1, GCH1, VEGFA COMT genes), and pharmacogenomics of antalgics and opioids (e.g., variants in OPRM1, STAT6, ABCB1, and COMT genes) in SCD. There has been greater progress made toward identifying the key genomic variants, mainly in BCL11A, HBS1L-MYB, or SAR1, which contribute to response to HU treatment. However, the complete picture on pharmacogenomic determinants of the above therapeutic phenotypes remains elusive. Strikingly, no study has been conducted in sub-Saharan Africa where majority of the patients with SCD live. This alerts the broader global life sciences community toward the existing disparities in optimal and ethical targeting of research and innovation investments for SCD specifically and precision medicine and pharmacology research broadly.

A genome-wide association study of kynurenic acid in cerebrospinal fluid: implications for psychosis and cognitive impairment in bipolar disorder

Elevated cerebrospinal fluid (CSF) levels of the glia-derived N-methyl-D-aspartic acid receptor antagonist kynurenic acid (KYNA) have consistently been implicated in schizophrenia and bipolar disorder. Here, we conducted a genome-wide association study based on CSF KYNA in bipolar disorder and found support for an association with a common variant within 1p21.3. After replication in an independent cohort, we linked this genetic variant-associated with reduced SNX7 expression-to positive psychotic symptoms and executive function deficits in bipolar disorder. A series of post-mortem brain tissue and in vitro experiments suggested SNX7 downregulation to result in a caspase-8-driven activation of interleukin-1β and a subsequent induction of the brain kynurenine pathway. The current study demonstrates the potential of using biomarkers in genetic studies of psychiatric disorders, and may help to identify novel drug targets in bipolar disorder.

Genome-wide association study suggests common variants within RP11-634B7.4 gene influencing severe pre-treatment pain in head and neck cancer patients

Pain is often one of the first signs of squamous cell carcinoma of the head and neck (HNSCC). Pain at diagnosis is an important prognostic marker for the development of chronic pain, and importantly, for the overall survival time. To identify variants influencing severe pre-treatment pain in 1,368 patients newly diagnosed with HNSCC, we conducted a genome-wide association study based on 730,525 tagging SNPs. The patients were all previously untreated for cancer. About 15% of the patients had severe pre-treatment pain, defined as pain score ≥7 (0 = “no pain” and 10 = “worst pain”). We identified 3 common genetic variants in high linkage disequilibrium for severe pre-treatment pain, representing one genomic region at 1q44 (rs3862188, P = 3.45 × 10⁻⁸; rs880143, P = 3.45 × 10⁻⁸; and rs7526880, P = 4.92 × 10⁻⁸), which maps to the RP11-634B7.4 gene, a novel antisense gene to three olfactory receptor genes. Olfactory receptor genes, upstream effectors of the MAPK signaling cascade, might be novel target genes for pain in HNSCC patients. Future experimental validation to explore biological mechanisms will be key to defining the role of the intronic variants and non-coding RNA for pain in patients with HNSCC.