PK-LR gene mutations in pyruvate kinase deficient Portuguese patients

PYK-SubstitutionOME: an integrated database containing allosteric coupling, ligand affinity and mutational, structural, pathological, bioinformatic and computational information about pyruvate kinase isozymes

Interpreting changes in patient genomes, understanding how viruses evolve and engineering novel protein function all depend on accurately predicting the functional outcomes that arise from amino acid substitutions. To that end, the development of first-generation prediction algorithms was guided by historic experimental datasets. However, these datasets were heavily biased toward substitutions at positions that have not changed much throughout evolution (i.e. conserved). Although newer datasets include substitutions at positions that span a range of evolutionary conservation scores, these data are largely derived from assays that agglomerate multiple aspects of function. To facilitate predictions from the foundational chemical properties of proteins, large substitution databases with biochemical characterizations of function are needed. We report here a database derived from mutational, biochemical, bioinformatic, structural, pathological and computational studies of a highly studied protein family-pyruvate kinase (PYK). A centerpiece of this database is the biochemical characterization-including quantitative evaluation of allosteric regulation-of the changes that accompany substitutions at positions that sample the full conservation range observed in the PYK family. We have used these data to facilitate critical advances in the foundational studies of allosteric regulation and protein evolution and as rigorous benchmarks for testing protein predictions. We trust that the collected dataset will be useful for the broader scientific community in the further development of prediction algorithms. Database URL https://github.com/djparente/PYK-DB.

Laboratory Approach to Hemolytic Anemia

Hemolytic anemias are a group of disorders with varied clinical and molecular heterogeneity. They are characterized by decreased levels of circulating erythrocytes in blood. The pathognomic finding is a reduced red cell life span with severe anemia or, compensated hemolysis accompanied by reticulocytosis. The diagnostic workup or laboratory approach for hemolytic anemias is based on methodical step-wise testing which includes red blood cell morphology, hematological indices with increased reticulocyte count along with clinical features of hemolytic anemias. If conventional laboratory tests are unable to detect the underlying cause of hemolysis, genetic testing is recommended. Sanger sequencing along with conventional testing is the most efficient way to diagnose the underlying genetic causes, especially in thalassemias/hemoglobinopathies, if required. However, hemolytic anemias being highly heterogeneous disorders, next-generation sequencing-based screening is rapidly becoming an efficient way to decipher the etiologies where common causes have been excluded.

Red blood cell PK deficiency: An update of PK-LR gene mutation database

Pyruvate kinase (PK) deficiency is known as being the most common cause of chronic nonspherocytic hemolytic anemia (CNSHA). Clinical PK deficiency is transmitted as an autosomal recessive trait, that can segregate neither in homozygous or in a compound heterozygous modality, respectively. Two PK genes are present in mammals: the pyruvate kinase liver and red blood cells (PK-LR) and the pyruvate kinase muscle (PK-M), of which only the first encodes for the isoenzymes normally expressed in the red blood cells (R-type) and in the liver (L-type). Several reports have been published describing a large variety of genetic defects in PK-LR gene associated to CNSHA. Herein, we present a review of about 250 published mutations and six polymorphisms in PK-LR gene with the corresponding clinical and molecular data. We consulted the PubMed website for searching mutations and papers, along with two main databases: the Leiden Open Variation Database (LOVD, https://grenada.lumc.nl/LOVD2/mendelian_genes/home.php?select_db=PKLR) and Human Gene Mutation Database (HGMD, http://www.hgmd.cf.ac.uk/ac/gene.php?gene=PKLR) for selecting, reviewing and listing the annotated PK-LR gene mutations present in literature. This paper is aimed to provide useful information to clinicians and laboratory professionals regarding overall reported PK-LR gene mutations, also giving the opportunity to harmonize data regarding PK-deficient individuals.

Molecular and clinical heterogeneity in pyruvate kinase deficiency in India

We studied the PK-LR gene in 10 unrelated Indian patients with congenital haemolytic anemia associated with erythrocyte pyruvate kinase deficiency. The patients had a variable presentation ranging from a very mild compensated hemolysis to severe anemia. Nine different mutations were detected among the 20 mutated alleles identified: one deletion (c.1042-1044del) p.Lys348del and eight single-nucleotide (nt) substitutions resulting in amino acid exchanges c.397A>G (p.Asn133Asp), c.992A>G (p.Asp331Gly), c.1072G>A (p.Gly358Arg), c.1076G>A (p.Arg359His), c.1219G>A (p.Glu407Lys), c.1241C>T (p.Pro414Leu), c.1436G>A (p.Arg479His) and c.1529G>A (p.Arg510Gln) were identified. Although all the exons, the flanking regions and the promoter region were sequenced in all cases, we failed to detect the second expected mutation in two subjects. Two mutations [c.397A>G; c.1241C>T] were novel. These novel missense mutations involved highly conserved amino acids. Two mutations were identified for the first time in the homozygous state globally (c1042-1044del; c.1072G>A) and two other mutations were identified for the first time in our population (c.1076G>A; c.1529G>A). This study along with our earlier report suggests that the most frequent mutations in India would appear to be c.1436G>A (18.33%), followed by c.992A>G (11.66%) and c.1456C>T (11.66%). Structural implications of amino acid substitutions were correlated with the clinical phenotypes seen.

Pyruvate kinase deficiency in sub-Saharan Africa: identification of a highly frequent missense mutation (G829A;Glu277Lys) and association with malaria

Background

Pyruvate kinase (PK) deficiency, causing hemolytic anemia, has been associated to malaria protection and its prevalence in sub-Saharan Africa is not known so far. This work shows the results of a study undertaken to determine PK deficiency occurrence in some sub-Saharan African countries, as well as finding a prevalent PK variant underlying this deficiency.

Materials and Methods

Blood samples of individuals from four malaria endemic countries (Mozambique, Angola, Equatorial Guinea and Sao Tome and Principe) were analyzed in order to determine PK deficiency occurrence and detect any possible high frequent PK variant mutation. The association between this mutation and malaria was ascertained through association studies involving sample groups from individuals showing different malaria infection and outcome status.

Results

The percentage of individuals showing a reduced PK activity in Maputo was 4.1% and the missense mutation G829A (Glu277Lys) in the PKLR gene (only identified in three individuals worldwide to date) was identified in a high frequency. Heterozygous carrier frequency was between 6.7% and 2.6%. A significant association was not detected between either PK reduced activity or allele 829A frequency and malaria infection and outcome, although the variant was more frequent among individuals with uncomplicated malaria.

Conclusions

This was the first study on the occurrence of PK deficiency in several areas of Africa. A common PKLR mutation G829A (Glu277Lys) was identified. A global geographical co-distribution between malaria and high frequency of PK deficiency seems to occur suggesting that malaria may be a selective force raising the frequency of this 277Lys variant.

Malaria: looking for selection signatures in the human PKLR gene region

The genetic component of susceptibility to malaria is both complex and multigenic and the better-known protective polymorphisms are those involving erythrocyte-specific structural proteins and enzymes. In vivo and in vitro data have suggested that pyruvate kinase deficiency, which causes a nonspherocytic haemolytic anaemia, could be protective against malaria severity in humans, but this hypothesis remains to be tested. In the present study, we conducted a combined analysis of Short Tandem Repeats (STRs) and Single Nucleotide Polymorphisms (SNPs) in the pyruvate kinase-encoding gene (PKLR) and adjacent regions (chromosome 1q21) to look for malaria selective signatures in two sub-Saharan African populations from Angola and Mozambique, in several groups with different malaria infection outcome. A European population from Portugal, including a control and a pyruvate kinase-deficient group, was used for comparison. Data from STR and SNP loci spread along the PKLR gene region showed a considerably higher differentiation between African and Portuguese populations than that usually found for neutral markers. In addition, a wider region showing strong linkage disequilibrium was found in an uncomplicated malaria group, and a haplotype was found to be associated with this clinical group. Altogether, this data suggests that malaria selective pressure is acting in this genomic region.

Analysis of malaria associated genetic traits in Cabo Verde, a melting pot of European and sub Saharan settlers

Malaria has occurred in the Cabo Verde archipelago with epidemic characteristics since its colonization. Nowadays, it occurs in Santiago Island alone and though prophylaxis is not recommended by the World Health Organization, studies have highlight the prospect of malaria becoming a serious public health problem as a result of the presence of antimalarial drug resistance associated with mutations in the parasite populations and underscore the need for tighter surveillance. Despite the presumptive weak immune status of the population, severe symptoms of malaria are not observed and many people present a subclinical course of the disease. No data on the prevalence of sickle-cell trait and red cell glucose-6-phosphate dehydrogenase deficiency (two classical genetic factors associated with resistance to severe malaria) were available for the Cabo Verde archipelago and, therefore, we studied the low morbidity from malaria in relation to the particular genetic characteristics of the human host population. We also included the analysis of the pyruvate kinase deficiency associated gene, reported as putatively associated with resistance to the disease. Allelic frequencies of the polymorphisms examined are closer to European than to African populations and no malaria selection signatures were found. No association was found between the analyzed human factors and infection but one result is of high interest: a linkage disequilibrium test revealed an association of distant loci in the PKLR gene and adjacent regions, only in non-infected individuals. This could mean a more conserved gene region selected in association to protection against the infection and/or the disease.

Spectrum of novel mutations in the human PKLR gene in pyruvate kinase-deficient Indian patients with heterogeneous clinical phenotypes

Eighteen unrelated pyruvate kinase (PK)-deficient Indian patients were identified in the past 4 years with varied clinical phenotypes ranging from a mild chronic haemolytic anaemia to a severe transfusion-dependent disorder. We identified 17 different mutations in the PKLR gene among the 36 mutated alleles. Ten novel mutations were identified: 427G>A, 499C>A, 1072G>A, 1180G>T, 1216G>A, 1220A>G, 644delG, IVS5 (+20) C>A, IVS9 (+44) C>T, and IVS9 (+93) A>C. A severe syndrome was commonly associated with some mutations, 992A>G, 1436G>A, 1220A>G, 644delG and IVS9 (+93) A>C, in the PKLR gene. Molecular graphics analysis of human red blood cell PK (RPK), based on the crystal structure of human PK, shows that mutations located near the substrate or fructose 1,6-diphosphate binding site may change the conformation of the active site, resulting in very low PK activity and severe clinical symptoms. The mutations target distinct regions of RPK structure, including domain interfaces and catalytic and allosteric sites. In particular, the 1216G>A and 1219G>A mutations significantly affect the interdomain interaction because they are located near the catalytic site in the A/B interface domains. The most frequent mutations in the Indian population appear to be 1436G>A (19.44%), followed by 1456C>T (16.66%) and 992A>G (16.66%). This is the first study to correlate the clinical profile with the molecular defects causing PK deficiency from India where 10 novel mutations that produce non-spherocytic haemolytic anaemia were identified.

Pyruvate kinase deficiency associated with severe liver dysfunction in the newborn

Significant hyperbilirubinaemia, anemia, and splenomegaly are common features in patients with severe haemolysis due to pyruvate kinase (PK) deficiency. Until now, severe neonatal PK deficiency has not been associated with fatal liver disease at this age. We present two neonatal cases of severe PK deficiency complicated with progressive fatal liver disease. The patients presented with severe haemolysis, progressive cholestasis, and hepatosplenomegaly, and both patients ultimately developed liver failure at a very young age. Despite extensive investigations, no specific explanation for liver disease and failure was found. We suggest that the PK deficiency itself directly led to liver dysfunction.

Pyruvate kinase deficiency: the genotype-phenotype association

Red cell pyruvate kinase (PK) deficiency is the most frequent enzyme abnormality of glycolysis causing chronic non-spherocytic haemolytic anaemia. The disease is transmitted as an autosomal recessive trait, clinical symptoms usually occurring in compound heterozygotes for two mutant alleles and in homozygotes. The severity of haemolysis is highly variable, ranging from very mild or fully compensated forms to life-threatening neonatal anaemia necessitating exchange transfusions. Erythrocyte PK is synthesised under the control of the PK-LR gene located on chromosome 1. One hundred eighty different mutations in PK-LR gene, mostly missense, have been so far reported associated to PK deficiency. First attempts to delineate the genotype-phenotype association were mainly based on the analysis of the enzyme's three-dimensional structure and the observation of the few homozygous patients. More recently, the comparison of the recombinant mutants of human red cell PK with the wild-type enzyme has enabled the effects of amino acid replacements on the enzyme molecular properties to be determined. However, the clinical manifestations of red cell enzyme defects are not merely dependent on the molecular properties of the mutant protein but rather reflect the complex interactions of additional factors, including genetic background, concomitant functional polymorphisms of other enzymes, posttranslational or epigenetic modifications, ineffective erythropoiesis and differences in splenic function.

Prenatal diagnosis for a novel homozygous mutation inPKLR gene in an Indian family

OBJECTIVE

To provide prenatal diagnosis of pyruvate kinase deficiency by direct DNA analysis in an Indian family.

MATERIALS AND METHOD

This case report describes diagnosis of a novel homozygous mutation in PKLR gene that subsequently helped the family in the next pregnancy.

RESULTS

Advancement in molecular genetics has resulted in the prenatal diagnosis of relatively uncommon genetic disorders like pyruvate kinase deficiency.

CONCLUSION

This case reiterates the importance of application of molecular genetics in clinical practice and prenatal diagnosis especially for rare, incurable genetic disorders.

Red cell pyruvate kinase deficiency: molecular and clinical aspects

Red cell pyruvate kinase (PK) deficiency is the most frequent enzyme abnormality of the glycolytic pathway causing hereditary non-spherocytic haemolytic anaemia. The degree of haemolysis varies widely, ranging from very mild or fully compensated forms, to life-threatening neonatal anaemia and jaundice necessitating exchange transfusions. Erythrocyte PK is synthesized under the control of the PK-LR gene located on chromosome 1. To date, more than 150 different mutations in the PK-LR gene have been associated with PK deficiency. First attempts to delineate the biochemical and clinical consequences of the molecular defect were mainly based on the observation of the few homozygous patients and on the analysis of the three-dimensional structure of the enzyme. More recently, the comparison of the recombinant mutants of human red cell PK with the wild-type enzyme has enabled the effects of amino acid replacements on the enzyme molecular properties to be determined and help to correlate genotype to clinical phenotype.

Molecular characterization of the PK-LR gene in sixteen pyruvate kinase-deficient patients

We studied the PK-LR gene in 16 unrelated patients with congenital haemolytic anaemia associated with erythrocyte pyruvate kinase deficiency. Fifteen different mutations were detected among the 28 mutated alleles identified: two deletions (del 1010G, del 1042--1044); one four nucleotide duplication (nt 1515--1518, GGTC); one splice site [IVS6(-2)t]; nine missense (991A, 1003A, 1151T, 1160G, 1181T, 1181A, 1456T, 1483A, 1529A); and two nonsense (721T, 1675T) mutations. Eight of them [del 1010G, del 1042--1044, dupl 1515--1518, IVS6(-2)t, 1003A, 1160G, 1181T, 1181A] were novel. The deletion 1042-1044 causes the loss of Lys 348. Deletion 1010G and duplication 1515-1518 determine a frameshift and the creation of a stop codon at nucleotides 1019 and 1554 respectively. Mutation IVS6(-2)t leads to an alteration of the 5' and 3' splice site consensus sequence; the cDNA analysis shows a 67-bp deletion in the first part of exon 11 (del 1437--1503). All the four new missense mutations involve highly conserved amino acids. The most frequent mutation in Italy would appear to be 1456T. Correlation was made between mutations, biochemical characteristics of the enzyme and clinical course of the disease.

A new PKLR gene mutation in the R-type promoter region affects the gene transcription causing pyruvate kinase deficiency

Mutations in the PKLR gene responsible for pyruvate kinase (PK)-deficient anaemia are mainly located in the coding regions: 11 are in the splicing sites and, recently, three mutations have been described in the promoter region. We now report a novel point mutation A-->G on nucleotide 72, upstream from the initiation codon of the PKLR gene, in four Portuguese PK-deficient patients. This new regulatory mutation occurs within the most proximal of the four GATA motifs (GATA-A element) in the R-type promoter region. In two patients who were homozygous for this mutation, a semiquantitative reverse transcription polymerase chain reaction (PCR) procedure was used to evaluate the amount of R-PK mRNA transcript in the reticulocytes. The mRNA level was about five times lower than in normal controls, demonstrating that the PKLR gene transcription is severely affected, most probably because the -72A-->G point mutation disables the binding of the erythroid transcription factor GATA-1 to the GATA-A element. Supporting these data, the two patients homozygous for the -72A-->G mutation had severe haemolytic anaemia and were transfusion dependent until splenectomy. Two other patients who were compound heterozygous for this mutation and the previously described missense mutation 1456C-->T had a mild condition.

Estimating the prevalence of pyruvate kinase deficiency from the gene frequency in the general white population

Pyruvate kinase (PK) deficiency is the most common cause of hereditary nonspherocytic hemolytic anemia. The prevalence of this deficiency is unknown, though some estimates have been made based on the frequency of low red cell PK activity in the population. An additional 20 patients with hereditary nonspherocytic hemolytic anemia caused by PK deficiency have been genotyped. One previously unreported mutation 1153C→T (R385W) was encountered. The relative frequency of PK mutations in patients with hemolytic anemia caused by PK deficiency was calculated from the 18 white patients reported here and from 102 patients previously reported in the literature. DNA samples from 3785 subjects from different ethnic groups have been screened for the 4 more frequently encountered mutations—c.1456 C→T(1456T), c.1468 C→T(1468T), c.1484 C→T(1484T), and c.1529 G6A (1529A)—by allele-specific oligonucleotide hybridization. Among white patients the frequency of the 1456T mutation was 3.50 × 10−3; that of the 1529A mutation was 2.03 × 10−3. Among African Americans the frequency of the 1456T mutation was 3.90 × 10−3 The only mutation found in the limited number of Asians tested was 1468T at a frequency of 7.94 × 10−3. Based on the gene frequency of the 1529A mutation in the white population and on its relative abundance in patients with hemolytic anemia caused by PK deficiency, the prevalence of PK deficiency is estimated at 51 cases per million white population. This number would be increased by inbreeding and decreased by failure of patients with PK deficiency to survive.

Estimating the prevalence of pyruvate kinase deficiency from the gene frequency in the general white population

Pyruvate kinase (PK) deficiency is the most common cause of hereditary nonspherocytic hemolytic anemia. The prevalence of this deficiency is unknown, though some estimates have been made based on the frequency of low red cell PK activity in the population. An additional 20 patients with hereditary nonspherocytic hemolytic anemia caused by PK deficiency have been genotyped. One previously unreported mutation 1153C→T (R385W) was encountered. The relative frequency of PK mutations in patients with hemolytic anemia caused by PK deficiency was calculated from the 18 white patients reported here and from 102 patients previously reported in the literature. DNA samples from 3785 subjects from different ethnic groups have been screened for the 4 more frequently encountered mutations—c.1456 C→T(1456T), c.1468 C→T(1468T), c.1484 C→T(1484T), and c.1529 G6A (1529A)—by allele-specific oligonucleotide hybridization. Among white patients the frequency of the 1456T mutation was 3.50 × 10−3; that of the 1529A mutation was 2.03 × 10−3. Among African Americans the frequency of the 1456T mutation was 3.90 × 10−3 The only mutation found in the limited number of Asians tested was 1468T at a frequency of 7.94 × 10−3. Based on the gene frequency of the 1529A mutation in the white population and on its relative abundance in patients with hemolytic anemia caused by PK deficiency, the prevalence of PK deficiency is estimated at 51 cases per million white population. This number would be increased by inbreeding and decreased by failure of patients with PK deficiency to survive.

Eight novel mutations and consequences on mRNA and protein level in pyruvate kinase-deficient patients with nonspherocytic hemolytic anemia

Pyruvate kinase (PK) deficiency (PKD) is an autosomal recessive disorder with the typical manifestation of nonspherocytic hemolytic anemia. We analyzed the mutant enzymes of 10 unrelated patients with PKD, whose symptoms ranged from a mild, chronic hemolytic anemia to a severe anemia, by sequence analysis for the presence of alterations in the PKLR gene. In all cases the patients were shown to be compound heterozygous. Eight novel mutations were identified: 458T-->C (Ile153Thr), 656T-->C (Ile219Thr), 877G-->A (Asp293Asn), 991G-->A (Asp331Asn), 1055C-->A (Ala352Asp), 1483G-->A (Ala495Thr), 1649A-->T (Asp550Val), and 183-184ins16bp. This 16 bp duplication produces a frameshift and subsequent stop codon resulting in a drastically reduced mRNA level, and probably in an unstable gene product. Surprisingly, the existence of M2-type PK could be demonstrated in the patient's red blood cells. The study of different polymorphic sites revealed, with one exception, a strict linkage of the 1705C, 1738T, IVS5(+51)T, T(10) polymorphisms and the presence of 14 ATT repeats in intron 11. Our analyses show the consequences of a distorted structure on enzyme function and we discuss the correlations between the mutations identified and the parameters indicative for enzyme function.

Biochemical and genetic basis of red cell enzyme deficiencies

Enzyme deficiencies have been identified in all erythrocyte pathways. Their frequencies differ with respect to the affected enzyme, the severity of the clinical manifestations and the geographical distribution. Most mutations are found within the coding sequences of genes, missense mutations occurring more often than deletions, insertions, splice site defects or premature stop codons. Promoter mutations are rare. The clinical manifestations are chronic or non-chronic haemolytic anaemias. The first of these are characterized by an impairment of cell function at normal values of the external load parameters kATPase and kGSHox. Haemolysis with a non-chronic course is induced only at enhanced values of the load parameters, caused by free radical generation by oxidative drugs, fava beans, infections, fever and physical exercise. The development of secondary haemochromatosis is the most common cause of mortality in patients suffering from severe chronic non-spherocytic haemolytic anaemia. Intracellular iron deposits must be prevented by timely treatment with effective chelating agents.

Red cell pyruvate kinase deficiency: from genetics to clinical manifestations

Pyruvate kinase deficiency is the most frequent enzyme abnormality of the Embden-Meyerhof pathway causing hereditary non-spherocytic haemolytic anaemia. The degree of haemolysis varies widely, ranging from very mild or fully compensated forms, to life-threatening neonatal anaemia and jaundice necessitating exchange transfusions. Splenectomy should be reserved for young patients who require regular blood transfusions. The gene encoding for pyruvate kinase (PK-LR) has been localized to the long arm of chromosome I; the cDNA of R-type is 2060 bp long and codes for 574 amino acids. More than 130 different mutations, mostly missense, have so far been described in association with PK deficiency, 1529A and 1456T being considered to be the most common mutations in Caucasians. Analysis of the three-dimensional structure of the enzyme may help in predicting the severity of the molecular defect. Further data on clinical features of homozygous patients are needed, at least for some mutations, to allow a more precise genotype/phenotype correlation.