Whole transcriptome analysis of human erythropoietic cells during ontogenesis suggests a role of VEGFA gene as modulator of fetal hemoglobin and pharmacogenomic biomarker of treatment response to hydroxyurea in β-type hemoglobinopathy patients

Milk transcriptome biomarker identification to enhance feed efficiency and reduce nutritional costs in dairy ewes

In recent years, rising prices for high-quality protein-based feeds have significantly increased nutrition costs. Consequently, investigating strategies to reduce these expenses and improve feed efficiency (FE) have become increasingly important for the dairy sheep industry. This research investigates the impact of nutritional protein restriction (NPR) during prepuberty and FE on the milk transcriptome of dairy Assaf ewes (sampled during the first lactation). To this end, we first compared transcriptomic differences between NPR and control ewes. Subsequently, we evaluated gene expression differences between ewes with divergent FE, using feed conversion ratio (FCR), residual feed intake (RFI), and consensus classifications of high- and low-FE animals for both indices. Lastly, we assess milk gene expression as a predictor of FE phenotype using random forest. No effect was found for the prepubertal NPR on milk performance or FE. Moreover, at the milk transcriptome level, only one gene, HBB, was differentially expressed between the NPR (n = 14) and the control group (n = 14). Further, the transcriptomic analysis between divergent FE sheep revealed 114 differentially expressed genes (DEGs) for RFI index (high-FERFI = 10 vs low-FERFI = 10), 244 for FCR (high-FEFCR = 10 vs low-FEFCR = 10), and 1 016 DEGs between divergent consensus ewes for both indices (high-FEconsensus = 8 vs low-FEconsensus = 8). These results underscore the critical role of selected FE indices for RNA-Seq analyses, revealing that consensus divergent animals for both indices maximise differences in transcriptomic responses. Genes overexpressed in high-FEconsensus ewes were associated with milk production and mammary gland development, while low-FEconsensus genes were linked to higher metabolic expenditure for tissue organisation and repair. The best prediction accuracy for FE phenotype using random forest was obtained for a set of 44 genes consistently differentially expressed across lactations, with Spearman correlations of 0.37 and 0.22 for FCR and RFI, respectively. These findings provide insights into potential sustainability strategies for dairy sheep, highlighting the utility of transcriptomic markers as FE proxies.

Do Genetic Polymorphisms Affect Fetal Hemoglobin (HbF) Levels in Patients With Sickle Cell Anemia Treated With Hydroxyurea? A Systematic Review and Pathway Analysis

Hydroxyurea has long been used for the treatment of sickle cell anemia (SCA), and its clinical effectiveness is related to the induction of fetal hemoglobin (HbF), a major modifier of SCA phenotypes. However, there is substantial variability in response to hydroxyurea among patients with SCA. While some patients show an increase in HbF levels and an ameliorated clinical condition under low doses of hydroxyurea, other patients present a poor effect or even develop toxicity. However, the effects of genetic polymorphisms on increasing HbF levels in response to hydroxyurea in patients with SCA (Hb SS) have been less explored. Therefore, we performed a systematic review to assess whether single-nucleotide polymorphisms (SNPs) affect HbF levels in patients with SCA treated with hydroxyurea. Moreover, we performed pathway analysis using the set of genes with SNPs found to be associated with changes in HbF levels in response to hydroxyurea among the included studies. The systematic literature search was conducted on Medline/PubMed, EMBASE, Cochrane Central Register of Controlled Trials, Cumulative Index to Nursing and Allied Health Literature (CINAHL), Scopus, and Web of Science. Seven cohort studies were included following our inclusion and exclusion criteria. From the 728 genetic polymorphisms examined in the included studies, 50 different SNPs of 17 genes were found to be associated with HbF changes in patients with SCA treated with hydroxyurea, which are known to affect baseline HbF but are not restricted to them. Enrichment analysis of this gene set revealed reactome pathways with the lowest adjusted p-values and highest combined scores related to VEGF ligand–receptor interactions (R-HSA-194313; R-HSA-195399) and the urea cycle (R-HSA-70635). Pharmacogenetic studies of response to hydroxyurea therapy in patients with SCA are still scarce and markedly heterogeneous regarding candidate genes and SNPs examined for association with HbF changes and outcomes, suggesting that further studies are needed. The reviewed findings highlighted that similar to baseline HbF, changes in HbF levels upon hydroxyurea therapy are likely to be regulated by multiple loci. There is evidence that SNPs in intron 2 of BCL11A affect HbF changes in response to hydroxyurea therapy, a potential application that might improve the clinical management of SCA.

Systematic Review Registration: (https://www.crd.york.ac.uk/prospero/display_record.php?RecordID=208790).

Genome-based therapeutic interventions for β-type hemoglobinopathies

For decades, various strategies have been proposed to solve the enigma of hemoglobinopathies, especially severe cases. However, most of them seem to be lagging in terms of effectiveness and safety. So far, the most prevalent and promising treatment options for patients with β-types hemoglobinopathies, among others, predominantly include drug treatment and gene therapy. Despite the significant improvements of such interventions to the patient's quality of life, a variable response has been demonstrated among different groups of patients and populations. This is essentially due to the complexity of the disease and other genetic factors. In recent years, a more in-depth understanding of the molecular basis of the β-type hemoglobinopathies has led to significant upgrades to the current technologies, as well as the addition of new ones attempting to elucidate these barriers. Therefore, the purpose of this article is to shed light on pharmacogenomics, gene addition, and genome editing technologies, and consequently, their potential use as direct and indirect genome-based interventions, in different strategies, referring to drug and gene therapy. Furthermore, all the latest progress, updates, and scientific achievements for patients with β-type hemoglobinopathies will be described in detail.

Epigenetic inactivation of ERF reactivates γ-globin expression in β-thalassemia

The fetal-to-adult hemoglobin switch is regulated in a developmental stage-specific manner and reactivation of fetal hemoglobin (HbF) has therapeutic implications for treatment of β-thalassemia and sickle cell anemia, two major global health problems. Although significant progress has been made in our understanding of the molecular mechanism of the fetal-to-adult hemoglobin switch, the mechanism of epigenetic regulation of HbF silencing remains to be fully defined. Here, we performed whole-genome bisulfite sequencing and RNA sequencing analysis of the bone marrow-derived GYPA⁺ erythroid cells from β-thalassemia-affected individuals with widely varying levels of HbF groups (HbF ≥ 95th percentile or HbF ≤ 5th percentile) to screen epigenetic modulators of HbF and phenotypic diversity of β-thalassemia. We identified an ETS2 repressor factor encoded by ERF, whose promoter hypermethylation and mRNA downregulation are associated with high HbF levels in β-thalassemia. We further observed that hypermethylation of the ERF promoter mediated by enrichment of DNMT3A leads to demethylation of γ-globin genes and attenuation of binding of ERF on the HBG promoter and eventually re-activation of HbF in β-thalassemia. We demonstrated that ERF depletion markedly increased HbF production in human CD34⁺ erythroid progenitor cells, HUDEP-2 cell lines, and transplanted NCG-Kit-V831M mice. ERF represses γ-globin expression by directly binding to two consensus motifs regulating γ-globin gene expression. Importantly, ERF depletion did not affect maturation of erythroid cells. Identification of alterations in DNA methylation of ERF as a modulator of HbF synthesis opens up therapeutic targets for β-hemoglobinopathies.

Genomic variants in members of the Krüppel-like factor gene family are associated with disease severity and hydroxyurea treatment efficacy in β-hemoglobinopathies patients

Aim: β-Type hemoglobinopathies are characterized by vast phenotypic diversity as far as disease severity is concerned, while differences have also been observed in hydroxyurea (HU) treatment efficacy. These differences are partly attributed to the residual expression of fetal hemoglobin (HbF) in adulthood. The Krüppel-like family of transcription factors (KLFs) are a set of zinc finger DNA-binding proteins which play a major role in HbF regulation. Here, we explored the possible association of variants in KLF gene family members with response to HU treatment efficacy and disease severity in β-hemoglobinopathies patients. Materials & methods: Six tag single nucleotide polymorphisms, located in four KLF genes, namely KLF3, KLF4, KLF9 and KLF10, were analyzed in 110 β-thalassemia major patients (TDT), 18 nontransfusion dependent β-thalassemia patients (NTDT), 82 sickle cell disease/β-thalassemia compound heterozygous patients and 85 healthy individuals as controls. Results: Our findings show that a KLF4 genomic variant (rs2236599) is associated with HU treatment efficacy in sickle cell disease/β-thalassemia compound heterozygous patients and two KLF10 genomic variants (rs980112, rs3191333) are associated with persistent HbF levels in NTDT patients. Conclusion: Our findings provide evidence that genomic variants located in KLF10 gene may be considered as potential prognostic biomarkers of β-thalassemia clinical severity and an additional variant in KLF4 gene as a pharmacogenomic biomarker, predicting response to HU treatment.

Omics Studies in Hemoglobinopathies

Hemoglobinopathies include all genetic diseases of hemoglobin and are grouped into thalassemia syndromes and structural hemoglobin variants. The β-thalassemias constitute a group of severe anemias with monogenic inheritance, caused by β-globin gene mutations. This review is focused on omics studies in hemoglobinopathies and mainly β-thalassemia, and discusses genomic, epigenomic, transcriptomic, proteomic and metabolomic findings. Omics analyses have identified various disease modifiers with an impact on disease severity and efficacy of treatments. These modifiers have contributed to the understanding of globin genes regulation/hemoglobin switching and the development of novel therapies. How omics data and their integration can contribute to efficient patient stratification, therapeutic management, improvements in existing treatments and application of novel personalized therapies is discussed.

Role of Genomic Biomarkers in Increasing Fetal Hemoglobin Levels Upon Hydroxyurea Therapy and in β-Thalassemia Intermedia: A Validation Cohort Study

Hemoglobinopathies exhibit a remarkable phenotypic diversity in terms of disease severity, while individual genetic background plays a key role in differential response to drug treatment. In the last decade, genomic variants in genes located within, as well as outside the human β-globin cluster have been shown to be significantly associated with Hb F increase, in relation to hydroxyurea (HU) therapy in patients with these diseases. Here, we aim to determine the effect of genomic variants located in genes, such as MAP3K5, ASS1, NOS2A, TOX, PDE7B, NOS1, FLT1 and ARG2, previously shown to modulate fetal hemoglobin (Hb F) levels in patients with β type hemoglobinopathies and reflecting disease severity and response to HU therapy in an independent cohort of Greek patients with these diseases. We recruited and genotyped 45 β-thalassemia patients (β-thal), either transfusion-dependent (TDT) or non transfusion-dependent (NTDT), 42 Hb S (HBB: c.20A>T)-β-thal compound heterozygotes, who were treated with HU, as well as 53 healthy individuals, all of Hellenic origin. Our study showed that genomic variants of the MAP3K5, NOS2A and ARG2 gene are associated with HU therapy efficacy in Hb S-β-thal compound heterozygotes. We have also shown that FLT1 and ARG2 genomic variants are associated with the mild phenotype of NTDT patients. Our findings provide evidence that MAP3K5, NOS2A, ARG2 and FLT1 genomic variants could be considered as genomic biomarkers to predict HU therapy efficacy in Hb S-β-thal compound heterozygotes and also to describe disease severity in patients with β type hemoglobinopathies.

Modelling human haemoglobin switching

Genetic lesions of the β-globin gene result in haemoglobinopathies such as β-thalassemia and sickle cell disease. To discover and test new molecular medicines for β-haemoglobinopathies, cell-based and animal models are now being widely utilised. However, multiple in vitro and in vivo models are required due to the complex structure and regulatory mechanisms of the human globin gene locus, subtle species-specific differences in blood cell development, and the influence of epigenetic factors. Advances in genome sequencing, gene editing, and precision medicine have enabled the first generation of molecular therapies aimed at reactivating, repairing, or replacing silenced or damaged globin genes. Here we compare and contrast current animal and cell-based models, highlighting their complementary strengths, reflecting on how they have informed the scope and direction of the field, and describing some of the novel molecular and precision medicines currently under development or in clinical trial.