"What We Know and What We Do Not Know about Evolutionary Genetic Adaptation to High Altitude Hypoxia in Andean Aymaras".

Three well-studied populations living at high altitudes are Tibetans, Andeans (Aymaras and Quechuas), and Ethiopians. Unlike Tibetans and Ethiopians who have similar hemoglobin (Hb) levels as individuals living at sea level, Aymara Hb levels increase when living at higher altitudes. Our previous whole genome study of Aymara people revealed several selected genes that are involved in cardiovascular functions, but their relationship with Hb levels was not elucidated. Here, we studied the frequencies of known evolutionary-selected variants in Tibetan and Aymara populations and their correlation with high Hb levels in Aymara. We genotyped 177 Aymaras at three different altitudes: 400 m (Santa Cruz), 4000 m (La Paz), and 5000 m (Chorolque), and correlated the results with the elevation of residence. Some of the Tibetan-selected variants also exist in Aymaras, but at a lower prevalence. Two of 10 Tibetan selected variants of EPAS1 were found (rs13005507 and rs142764723) and these variants did not correlate with Hb levels. Allele frequencies of 5 Aymara selected SNPs (heterozygous and homozygous) at 4000 m (rs11578671_BRINP3, rs34913965_NOS2, rs12448902_SH2B1, rs10744822_TBX5, and rs487105_PYGM) were higher compared to Europeans. The allelic frequencies of rs11578671_BRINP3, rs34913965_NOS2, and rs10744822_SH2B1 were significantly higher for Aymaras living at 5000 m than those at 400 m elevation. Variant rs11578671, close to the BRINP3 coding region, correlated with Hb levels in females. Variant rs34913965 (NOS2) correlated with leukocyte counts. Variants rs12448902 (SH2B1) and rs34913965 (NOS2) associated with higher platelet levels. The correlation of these SNPs with blood cell counts demonstrates that the selected genetic variants in Aymara influence hematopoiesis and cardiovascular effects.

Clinical utility of targeted next-generation sequencing panel in routine diagnosis of hereditary hemolytic anemia: A national reference laboratory experience.

INTRODUCTION: Hereditary hemolytic anemias (HHA) comprise a heterogeneous group of disorders resulting from defective red blood cell (RBC) cytoskeleton, RBC enzyme deficiencies, and hemoglobin (Hb) synthesis disorders such as thalassemia or sideroblastic anemia. MATERIALS AND METHODS: Our hemolytic anemia diagnostic next-generation sequencing (NGS) panel includes 28 genes encoding RBC cytoskeletal proteins, membrane transporter, RBC enzymes, and certain bilirubin metabolism genes. The panel covers the complete coding region of these genes, splice junctions, and, wherever appropriate, deep intronic or regulatory regions are also included. Four hundred fifty-six patients with unexplained hemolytic anemia were evaluated using our NGS panel between 2015 and 2019. RESULTS: We identified pathogenic/likely pathogenic variants in 111/456 (24%) patients that were responsible for the disease phenotype (e.g., moderate to severe hemolytic anemia and hyperbilirubinemia). Approximately 40% of the mutations were novel. As expected, 45/456 (10%) patients were homozygous for the promoter polymorphism in the UGT1A1 gene, A(TA)7 TAA (UGT1A1*28). 8/45 homozygous UGT1A1*28 cases were associated with additional pathogenic mutations causing hemolytic anemia, likely exacerbating hyperbilirubinemia. The most common mutated genes were membrane cytoskeleton genes SPTA1, and SPTB, followed by PKLR. Complex interactions between SPTA1 low expression alleles, alpha-LELY and alpha-LEPRA alleles, and intragenic SPTA1 variants were associated with hereditary pyropoikilocytosis and autosomal recessive hereditary spherocytosis in 23/111 patients. CONCLUSIONS: Our results demonstrate that hemolytic anemia is underscored by complex molecular interactions of previously known and novel mutations in RBC cytoskeleton/enzyme genes, and therefore, NGS should be considered in all patients with clinically unexplained hemolytic anemia and in neonates with hyperbilirubinemia. Moreover, low expression alleles alpha-LELY and alpha-LEPRA should be included in all targeted HHA panels.

A Comprehensive Triple-Repeat Primed PCR and a Long-Range PCR Agarose-Based Assay for Improved Genotyping of Guanine-Adenine-Adenine Repeats in Friedreich Ataxia.

Friedreich ataxia is a rare autosomal recessive, neuromuscular degenerative disease caused by an expansion of a trinucleotide [guanine-adenine-adenine (GAA)] repeat in intron 1 of the FXN gene. It is common in the White population, characterized by progressive gait and limb ataxia, lack of tendon reflexes in the legs, loss of position sense, and hypertrophic cardiomyopathy. Detection and genotyping of the trinucleotide repeat length is important for the diagnosis and prognosis of the disease. A two-tier genotyping assay with an improved triple-repeat primed PCR (TR-PCR) for alleles <200 GAA repeats (±1 to 5 repeats) and an agarose gel-based, long-range PCR (LR-PCR) assay to genotype expanded alleles >200 GAA repeats (±50 repeats) is described. Of the 1236 DNA samples tested using TR-PCR, 31 were identified to have expanded alleles >200 repeats and were reflexed to the LR-PCR procedure for confirmation and quantification. The TR-PCR assay described herein is a diagnostic genotyping assay that reduces the need for further testing. The LR-PCR component is a confirmatory test for true homozygous and heterozygous samples with normal and expanded alleles, as indicated by the TR-PCR assay. The use of this two-tier method offers a comprehensive evaluation to detect and genotype the smallest and largest number of GAA repeats, improving the classification of FXN alleles as normal, mutable normal, borderline, and expanded alleles.

Oral Methylene Blue Treatment in A Dog with Cytochrome B5 Reductase Deficiency And 78, XX Testicular Disorder of Sex Development.

A 6-month-old mixed breed dog was referred for evaluation of a potential disorder of sex development (DSD) and lower than expected energy level. Genitourinary examination revealed ambiguous external genitalia, hypospadias, and a subtle pouch of skin that resembled an empty scrotum. Corrective surgery was planned and subsequently aborted after cyanosis was identified preoperatively and an arterial blood gas analysis by co-oximetry identified increased methemoglobin (MetHb) concentration (35%, normal <2%) with normal arterial oxygen tension. Ensuing investigations confirmed hereditary methemoglobinemia caused by cytochrome b5 reductase (CYB5R) deficiency via molecular genetic (Arg219Pro homozygous variant in CYB5R3 gene) and biochemical (cytochrome b5 reductase enzyme activity of 8% [normal, 100% activity] testing. Karyotyping and molecular analysis of sex chromosomes revealed the dog was genetically female with a normal female karyotype (78,XX), and was negative for the Y-linked SRY gene and positive for the X-linked androgen receptor gene. Methylene blue (MB, 3.3 mg/kg per os [PO] q24 h) was administered and the MetHb concentration decreased to 9% within 14 days. Urogenital revision surgery proceeded without complication and the dog was maintained on MB (3-4 mg/kg PO q24 h) long-term without adverse effects. This is the first report to describe the use of PO MB to decrease MetHb concentrations in a dog with CYB5R deficiency in preparation for anesthesia and highlights its potential as a viable alternative to the intravenous formulation for elective procedures. In addition, this report describes the clinical, molecular, imaging, surgical, and macroscopic and microscopic pathological features of a dog with SRY-negative, 78,XX testicular DSD.

A computational study of structural differences of binding of NADP+ and G6P substrates to G6PD Mediterraneanc.563T, G6PD A-c.202A/c.376G, G6PD Cairoc.404C and G6PD Gazac.536A mutations.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a common X-linked inherited enzymopathic disorder that may lead to transfusion-requiring acute hemolytic anemia (AHA) triggered by fava beans ingestion, infection or some drugs. The gene encoding for G6PD carries a large number of genetic variants that have varying pathogenicity. We reported on three G6PD variants in the Gaza Strip Palestinian population with differing clinical impacts and frequencies: G6PD Mediterraneanc.563T, African G6PD A-c.202A/c.376G, and G6PD Cairoc.404C. We also identified a novel G6PD missense (Ser179Asn) mutation c.536G > A "G6PD Gaza". In this work we explore the effect of these four genetic variants on the structural and substrate (NADP+ and G6P) binding characteristics of the G6PD enzyme using the Monte Carlo (MC) flexible docking and molecular dynamics (MD) simulation approaches. We report that G6PD A-c.202A/c.376G, G6PD Mediterraneanc.563T, G6PD Cairoc.404C and G6PD Gazac.536A mutations cause significant structural changes in G6PD enzyme to induce conformational instability leading to the loss of binding of one or both substrates and are causative of G6PD deficiency.

Clinical, metabolic, and molecular genetic characterization of hereditary methemoglobinemia caused by cytochrome b5 reductase deficiency in 30 dogs.

Genotype-phenotype correlations of humans and dogs with hereditary methemoglobinemia are not yet well characterized. We determined total hemoglobin and methemoglobin (MetHb) concentrations, cytochrome b5 reductase (CYB5R) enzyme activities, genotypes, and clinical signs in 30 dogs with persistent cyanosis without cardiopulmonary disease. Erythrocytic CYB5R enzyme activities were low in all dogs assayed. Owner-reported quality of life ranged from subclinical to occasional exertional syncope. Two previously reported and two novel CYB5R3 missense variants were identified among the methemoglobinemic cohort and were predicted to impair enzyme function. Two variants were recurrent: a homozygous Ile194Leu substitution was found in Pomeranians and other small dogs, and a homozygous Arg219Pro change occurred predominately in pit bull terriers. The other two variants were Thr202Ala and Gly76Ser substitutions in single dogs. Of the two common CYB5R3 genotypes, Arg219Pro was associated with a more severe metabolic phenotype. We conclude that CYB5R3 deficiency is the predominate cause of canine hereditary methemoglobinemia. Although this finding is unlikely to alter the clinical approach to hereditary methemoglobinemia in dogs, it demonstrates the possibility of how genotype-phenotype cohort analysis might facilitate precision medicine in the future in veterinary medicine.

A Novel Variant in G6PD (c.1375C>G) Identified from a Hispanic Neonate with Extreme Hyperbilirubinemia and Low G6PD Enzymatic Activity.

We report a novel glucose-6-phosphate dehydrogenase (G6PD) variant (c.1375C>G) discovered in a 3-day-old Hispanic male child from Salt Lake City, UT, USA. This newborn presented with severe hyperbilirubinemia (29.8 mg/dL or 510 µmol/L) and marked hemolysis evidenced by elevated end-tidal carbon monoxide concentration (5.9 ppm, normal <1.7 ppm). Despite a very low prevalence of G6PD deficiency in Hispanic populations, we pursued testing for this condition and found he had low erythrocyte G6PD enzyme activity (2.8 U/g Hb, normal 9.9-16.6 U/g Hb) and a novel G6PD variant. His mother was heterozygous for this same variant, and she had a moderate decrease in G6PD enzyme activity (7.1 U/g Hb). On the basis of these findings, we propose this variant as a novel pathogenic mutation.

The CYB5R3c .350C>G and G6PD A alleles modify severity of anemia in malaria and sickle cell disease.

Genetic modifiers of anemia in Plasmodium falciparum infection and sickle cell disease (SCD) are not fully known. Both conditions are associated with oxidative stress, hemolysis and anemia. The CYB5R3 gene encodes cytochrome b5 reductase 3, which converts methemoglobin to hemoglobin through oxidation of NADH. CYB5R3c.350C > G encoding CYB5R3T117S , the most frequent recognized African-specific polymorphism, does not have known functional significance, but its high allele frequency (23% in African Americans) suggests a selection advantage. Glucose-6-phosphate dehydrogenase (G6PD) is essential for protection from oxidants; its African-polymorphic X-linked A+ and A- alleles, and other variants with reduced activity, coincide with endemic malaria distribution, suggesting protection from lethal infection. We examined the association of CYB5R3c.350C > G with severe anemia (hemoglobin <5 g/dL) in the context of G6PD A+ and A- status among 165 Zambian children with malaria. CYB5R3c.350C > G offered protection against severe malarial anemia in children without G6PD deficiency (G6PD wild type or A+/A- heterozygotes) (odds ratio 0.29, P = .022) but not in G6PD A+ or A- hemizygotes/homozygotes. We also examined the relationship of CYB5R3c.350C > G with hemoglobin concentration among 267 children and 321 adults and adolescents with SCD in the US and UK and found higher hemoglobin in SCD patients without G6PD deficiency (ß = 0.29, P = .022 children; ß = 0.33, P = .004 adults). Functional studies in SCD erythrocytes revealed mildly lower activity of native CYB5R3T117S compared to wildtype CYB5R3 and higher NADH/NAD+ ratios. In conclusion, CYB5R3c.350C > G appears to ameliorate anemia severity in malaria and SCD patients without G6PD deficiency, possibly accounting for CYB5R3c.350C > G selection and its high prevalence.

Molecular Understanding of Severe Cases of ß-Thalassemia in the Nablus Region, West Bank, Palestine.

Hemoglobinopathies are common inherited monogenic diseases that are likely to remain a serious regional health problem where thalassemias and sickle cell disease are prevalent. In regions where recessive alleles for hemoglobinopathy disorders are present with high consanguinity rates, such as in Palestine, coinheritance of two different genetic defects becomes anticipated and prevalent. In this report, we characterize the molecular variants of the HBB gene for 16 patients with transfusion-dependent anemia registered at the Thalassemia Patient Friends Society in Nablus governorate, West Bank, Palestine. Analysis revealed that 63.0% (10/16) of the patients were homozygous for ß-thalassemia (ß-thal), IVS-I-6 (T>C) (HBB: c.92+6T>C) or IVS-I-110 (G>A) (HBB: c.93-21G>A); 19.0% (3/16) homozygous for sickle cell disease or Hb S (HBB: c.20A>T, p.Glu6Val); 13.0% (2/16) were double heterozygotes for Hb S/ß-thal, (HBB: c.20A>T/HBB: c.92G>C) and HBB: c.20A>T/HBB: c.321_322insG; and one case was a compound heterozygote for ß-thal, codon 39 (C>T) (HBB: c.118C>T) and IVS-I-110. The most common mutation reported in the 16 patients was IVS-I-6 (0.38), followed by IVS-I-110 (0.28) Hb S (0.25) and 0.03 each for codon 39, codons 106/107 (HBB: c.321_322insG) and Hb Monroe (HBB: c.92G>C). In conclusion, in Palestine, a variety of intricate inheritance patterns are encountered in clinical practice.

Genetic predisposition in type 2 diabetes: A promising approach toward a personalized management of diabetes.

Diabetes mellitus, also known simply as diabetes, has been described as a chronic and complex endocrine metabolic disorder that is a leading cause of death across the globe. It is considered a key public health problem worldwide and one of four important non-communicable diseases prioritized for intervention through world health campaigns by various international foundations. Among its four categories, Type 2 diabetes (T2D) is the commonest form of diabetes accounting for over 90% of worldwide cases. Unlike monogenic inherited disorders that are passed on in a simple pattern, T2D is a multifactorial disease with a complex etiology, where a mixture of genetic and environmental factors are strong candidates for the development of the clinical condition and pathology. The genetic factors are believed to be key predisposing determinants in individual susceptibility to T2D. Therefore, identifying the predisposing genetic variants could be a crucial step in T2D management as it may ameliorate the clinical condition and preclude complications. Through an understanding the unique genetic and environmental factors that influence the development of this chronic disease individuals can benefit from personalized approaches to treatment. We searched the literature published in three electronic databases: PubMed, Scopus and ISI Web of Science for the current status of T2D and its associated genetic risk variants and discus promising approaches toward a personalized management of this chronic, non-communicable disorder.

Clinical, metabolic, and genetic characterization of hereditary methemoglobinemia caused by cytochrome b5 reductase deficiency in cats.

Two non-pedigreed male castrated cats had persistent cyanosis over a 3-year observation period. Clinical cardiopulmonary evaluations did not reveal abnormalities, but the blood remained dark after exposure to air. Erythrocytic methemoglobin concentrations were high (~40% of hemoglobin) and cytochrome b5 reductase (CYB5R) activities in erythrocytes were low (≤15% of control). One cat remained intolerant of exertion, and the other cat developed anemia and died due to an unidentified comorbidity. Whole-genome sequencing revealed a homozygous c.625G>A missense variant (B4:137967506) and a c.232-1G>C splice acceptor variant (B4:137970815) in CYB5R3, respectively, which were absent in 193 unaffected additional cats. The p.Gly209Ser missense variant likely disrupts a nicotinamide adenine dinucleotide (NADH)-binding domain, while the splicing error occurs at the acceptor site for exon 4, which likely affects downstream translation of the protein. The 2 novel CYB5R3 variants were associated with methemoglobinemia using clinical, biochemical, genomics, and in silico protein studies. The variant prevalence is unknown in the cat population.

Three Novel Spectrin Variants in Jaundiced Neonates.

Various mutations in the genes encoding alpha spectrin (SPTA1) or beta spectrin (SPTB) are known to cause erythrocyte membrane disorders, sometimes associated with severe neonatal jaundice and anemia. We used a next-generation sequencing panel to evaluate 3 unrelated neonates who had puzzling cases of nonimmune hemolytic jaundice. In each case, we identified novel mutations in either SPTA1 or SPTB. Correlating erythrocyte morphology, clinical course, and computational analysis, we submit that each of the 3 variants is a probable pathogenic cause of the hereditary hemolytic conditions in these patients. We hope other pediatric practitioners caring for neonates with what appears to be idiopathic severe neonatal hyperbilirubinemia will look for spectrin variants as a possible cause, because additional cases with these specific variants along with this clinical phenotype are needed to confirm our postulate that these 3 cases are indeed pathogenic mutations.

Gain-of-function EGLN1 prolyl hydroxylase (PHD2 D4E:C127S) in combination with EPAS1 (HIF-2α) polymorphism lowers hemoglobin concentration in Tibetan highlanders.

Tibetans have lived at high altitude for generations and are thought to be genetically adapted to hypoxic environments. Most are protected from hypoxia-induced polycythemia, and a haplotype of EPAS1, encoding hypoxia-inducible factor (HIF-2α), has been associated with lower hemoglobin levels. We earlier reported a Tibetan-specific EGLN1 haplotype encoding PHD2 which abrogates HIF augmentation in hypoxia. We genotyped 347 Tibetan individuals from varying altitudes for both the Tibetan-specific EGLN1 haplotype and 10 candidate SNPs in the EPAS1 haplotype and correlated their association with hemoglobin levels. The effect of the EGLN1 haplotype on hemoglobin exhibited age dependency at low altitude, while at higher altitudes, it showed a trend to lower hemoglobin levels in the presence of the Tibetan-selected EPAS1 rs142764723 C/C allele. The observed gene-environment and gene-gene interactions and the moderate effect of the EGLN1 and EPAS1 haplotypes on hemoglobin indicate that other modifiers exist. It remains to be determined whether a blunting of erythropoiesis or other physiological consequences of HIF downregulation are the primary drivers of these genetic adaptations among Tibetans. KEY MESSAGE: Most Tibetans are protected from polycythemia while living in high altitude. An EGLN1 co-adapted haplotype, EGLN1 c.12C>G, c.380G>C is uniquely Tibetan. The Tibetan EPAS1 haplotype has introgressed from the Denisovan genome. While EGLN1 and EPAS1 genotypes lower Hb, this study indicates additional Hb modifiers.

Possible association of 3' UTR +357 A>G, IVS11-nt 93 T>C, c.1311 C>T polymorphism with G6PD deficiency.

BACKGROUND: Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a common X-linked inherited enzymopathic disorder affecting more than 500 million people worldwide. It has so far been linked to 217 distinct genetic variants in the exons and exon-intron boundaries of the G6PD gene, giving rise to a wide range of biochemical heterogeneity and clinical manifestations. OBJECTIVES: Reports from different settings suggested the association of intronic and other mutations outside the reading frame of the G6PD gene with reduced enzyme activity and presenting clinical symptoms. The present study aimed to investigate any association of other variations apart of the exonic or exonic intronic boundaries in the development of G6PD deficiency. METHODS: Sixty-seven unrelated Palestinian children admitted to the pediatric hospital with hemolytic crises due to G6PD deficiency were studied. RESULTS: In our Palestinian cohort of 67 [59 males (M) and 8 females (F)] G6PD-deficient children, previously hospitalized for acute hemolytic anemia due to favism, molecular sequencing of the G6PD gene revealed four cases (3M and 1F) that did not have any of the variants known to cause G6PD deficiency, but the 3' UTR c.*+357A>G (rs1050757) polymorphism in association with IVS 11 (c.1365-13T>C; rs2071429), and c.1311C>T (rs2230037). CONCLUSION: We now provide an additional evidence form Palestinian G6PD-deficient subjects for a possible role of 3' UTR c.*+357 A>G, c.1365-13T>C, and/or c.1311C>T polymorphism for G6PD deficiency, suggesting that not only a single variation in the exonic or exonic intronic boundaries, but also a haplotype of G6PD should considered as a cause for G6PD deficiency.

Evolutionary selected Tibetan variants of HIF pathway and risk of lung cancer.

Tibetans existed in high altitude for ~25 thousand years and have evolutionary selected unique haplotypes assumed to be beneficial to hypoxic adaptation. EGLN1/PHD2 and EPAS1/HIF-2α, both crucial components of hypoxia sensing, are the two best-established loci contributing to high altitude adaptation. The co-adapted Tibetan-specific haplotype encoding for PHD2:p.[D4E/C127S] promotes increased HIF degradation under hypoxic conditions. The Tibetan-specific 200 kb EPAS1 haplotype introgressed from an archaic human population related to Denisovans which underwent evolutionary decay; however, the functional variant(s) responsible for high-altitude adaptation at EPAS1/HIF-2α have not yet been identified. Since HIF modulates the behavior of cancer cells, we hypothesized that these Tibetan selected genomic variants may modify cancer risk predisposition. Here, we ascertained the frequencies of EGLN1D4E/C127S and EGLN1C127S variants and ten EPAS1/HIF-2α variants in lung cancer patients and controls in Nepal, whose population consists of people with Indo-Aryan origin and Tibetan-related Mongoloid origin. We observed a significant association between the selected Tibetan EGLN1/PHD2 haplotype and lung cancer (p=0.0012 for D4E, p=0.0002 for C127S), corresponding to a two-fold increase in lung cancer risk. We also observed a two-fold or greater increased risk for two of the ten EPAS1/HIF-2α variants, although the association was not significant after correcting for multiple comparisons (p=0.12). Although these data cannot address the role of these genetic variants on lung cancer initiation or progression, we conclude that some selected Tibetan variants are strongly associated with a modified risk of lung cancer.

Loss of Major DNase I Hypersensitive Sites in Duplicated ß-globin Gene Cluster Incompletely Silences HBB Gene Expression.

We report an infant with sickle cell disease phenotype by biochemical analysis whose ß-globin gene (HBB) sequencing showed sickle cell mutation (HBBS ) heterozygosity. The proband has a unique head-to-tail duplication of the ß-globin gene cluster having wild-type (HBBA ) and HBBS alleles inherited from her father; constituting her HBBS /HBBS -HBBA genotype. Further analyses revealed that proband's duplicated ß-globin gene cluster (∼650 kb) encompassing HBBA does not include the immediate upstream locus control region (LCR) or 3' DNase I hypersensitivity (HS) element. The LCR interacts with ß-globin gene cluster involving long range DNA interactions mediated by various transcription factors to drive the regulation of globin genes expression. However, a low level of HBBA transcript was clearly detected by digital PCR. In this patient, the observed transcription from the duplicated, distally displaced HBBA cluster demonstrates that the loss of LCR and flanking 3'HS sites do not lead to complete silencing of HBB transcription.

Favism, the commonest form of severe hemolytic anemia in Palestinian children, varies in severity with three different variants of G6PD deficiency within the same community.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a common genetic abnormality known to predispose to acute hemolytic anemia (AHA), which can be triggered by certain drugs or infection. However, the commonest trigger is fava beans (Vicia faba) ingestion, causing AHA (favism), which may be life-threatening especially in children. G6PD deficiency is genetically highly heterogeneous, as nearly 200 different mutations have been observed. We have investigated the hematological features of acute favism in the Palestinian Gaza community that is characterized by the polymorphic coexistence of three different G6PD deficiency genes (G6PD A-, G6PD Cairo, G6PD Med). We have found by comparison to the general population (485 adults and 466 newborns) that children with favism, in terms of relative frequency, G6PD A- was under-represented, whereas G6PD Med was over-represented. We also found that the severity of anemia was significantly greater with G6PD Med and G6PD Cairo than with G6PD A-; and with G6PD Cairo, compared to the other two variants, there was greater hyperbilirubinemia, as well as persistence of mild anemia and reticulocytosis for as long as 4months after recovery from favism. This is the first report determining a differential impact of different G6PD mutations on the clinical features of favism in the same population and the same environment.