The extent of linkage disequilibrium caused by selection on G6PD in humans

The gene coding for glucose-6-phosphate dehydrogenase (G6PD) is subject to positive selection by malaria in some human populations. The G6PD A- allele, which is common in sub-Saharan Africa, is associated with deficient enzyme activity and protection from severe malaria. To delimit the impact of selection on patterns of linkage disequilibrium (LD) and nucleotide diversity, we resequenced 5.1 kb at G6PD and approximately 2-3 kb at each of eight loci in a 2.5-Mb region roughly centered on G6PD in a diverse sub-Saharan African panel of 51 unrelated men (including 20 G6PD A-, 11 G6PD A+, and 20 G6PD B chromosomes). The signature of selection is evident in the absence of genetic variation at G6PD and at three neighboring loci within 0.9 Mb from G6PD among all individuals bearing G6PD A- alleles. A genomic region of approximately 1.6 Mb around G6PD was characterized by long-range LD associated with the A- alleles. These patterns of nucleotide variability and LD suggest that G6PD A- is younger than previous age estimates and has increased in frequency in sub-Saharan Africa due to strong selection (0.1 < s < 0.2). These results also show that selection can lead to nonrandom associations among SNPs over great physical and genetic distances, even in African populations.     

Contrasting histories of G6PD molecular evolution and malarial resistance in humans and chimpanzees

Although mutations in the glucose-6-phosphate dehydrogenase (G6PD) gene result in several blood-related diseases in humans, they also confer resistance to malarial infection. This association between G6PD and malaria was supported by population genetic analyses of the G6PD locus, which indicated that these mutations may have recently risen in frequency in certain geographic regions as a result of positive selection. Here we characterize nucleotide sequence variation in a 5.2-kb region of the G6PD locus in a population sample of 56 chimpanzees, as well as among 7 other nonhuman primates, to compare with that in humans in determining whether other primates that are impacted by malaria also exhibit patterns of G6PD polymorphism or divergence consistent with positive selection. We find that chimpanzees have several amino acid variants but that the overall pattern at G6PD in chimpanzees, as well as in Old and New World primates in general, can be explained by recent purifying selection as well as strong functional constraint dating back to at least 30-40 MYA. These comparative analyses suggest that the recent signature of positive selection at G6PD in humans is unique.     

Nucleotide variability at G6pd and the signature of malarial selection in humans

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common enzymopathy in humans. Deficiency alleles for this X-linked disorder are geographically correlated with historical patterns of malaria, and the most common deficiency allele in Africa (G6PD A-) has been shown to confer some resistance to malaria in both hemizygous males and heterozygous females. We studied DNA sequence variation in 5.1 kb of G6pd from 47 individuals representing a worldwide sample to examine the impact of selection on patterns of human nucleotide diversity and to infer the evolutionary history of the G6PD A- allele. We also sequenced 3.7 kb of a neighboring locus, L1cam, from the same set of individuals to study the effect of selection on patterns of linkage disequilibrium. Despite strong clinical evidence for malarial selection maintaining G6PD deficiency alleles in human populations, the overall level of nucleotide heterozygosity at G6pd is typical of other genes on the X chromosome. However, the signature of selection is evident in the absence of genetic variation among A- alleles from different parts of Africa and in the unusually high levels of linkage disequilibrium over a considerable distance of the X chromosome. In spite of a long-term association between Plasmodium falciparum and the ancestors of modern humans, patterns of nucleotide variability and linkage disequilibrium suggest that the A- allele arose in Africa only within the last 10,000 years and spread due to selection.     

The functional impact of Pgm amino acid polymorphism on glycogen content in Drosophila melanogaster

Earlier studies of the common PGM allozymes in Drosophila melanogaster reported no in vitro activity differences. However, our study of nucleotide variation observed that PGM allozymes are a heterogeneous mixture of amino acid polymorphisms. In this study, we analyze 10 PGM protein haplotypes with respect to PGM activity, thermostability, and adult glycogen content. We find a twofold difference in activity among PGM protein haplotypes that is associated with a threefold difference in glycogen content. The latitudinal clines for several Pgm amino acid polymorphisms show that high PGM activity, and apparently higher flux to glycogen synthesis, parallel the low activity clines at G6PD for reduced pentose shunt flux in northern latitudes. This suggests that amino acid polymorphism is under selection at this branch point and may be favored for increased metabolic storage associated with stress resistance and adaptation to temperate regions.     

Natural and synthetic alleles provide complementary insights into the nature of selection acting on the Men polymorphism of Drosophila melanogaster

Two malic enzyme alleles, Men^113A and Men^113G, occur at approximately equal frequency in North American populations of Drosophila melanogaster, while only Men^113A occurs in African populations. We investigated the population genetics, biochemical characteristics, and selective potential of these alleles. Comparable levels of nucleotide polymorphism in both alleles suggest that the Men^113G allele is not recently derived, but we find no evidence in the DNA sequence data for selection maintaining the polymorphism. Interestingly, the alleles differ in both V_max and K_m for the substrate malate. Triglyceride concentration and isocitrate dehydrogenase (IDH) and glucose-6-phosphate dehydrogenase (G6PD) activities are negatively correlated with the in vivo activities of the Men alleles. We examined the causality of the observed correlations using P-element excision-derived knockout alleles of the Men gene and found significant changes in the maximum activities of both IDH and G6PD, but not in triglyceride concentration, suggesting compensatory interactions between MEN, IDH, and G6PD. Additionally, we found significantly higher than expected levels of MEN activity in knockout heterozygotes, which we attribute to transvection effects. The distinct differences in biochemistry and physiology between the naturally occurring alleles and between the engineered alleles suggest the potential for selection on the Men locus.     

Extensive genome-wide variability of human cytomegalovirus in congenitally infected infants

Research has shown that RNA virus populations are highly variable, most likely due to low fidelity replication of RNA genomes. It is generally assumed that populations of DNA viruses will be less complex and show reduced variability when compared to RNA viruses. Here, we describe the use of high throughput sequencing for a genome wide study of viral populations from urine samples of neonates with congenital human cytomegalovirus (HCMV) infections. We show that HCMV intrahost genomic variability, both at the nucleotide and amino acid level, is comparable to many RNA viruses, including HIV. Within intrahost populations, we find evidence of selective sweeps that may have resulted from immune-mediated mechanisms. Similarly, genome wide, population genetic analyses suggest that positive selection has contributed to the divergence of the HCMV species from its most recent ancestor. These data provide evidence that HCMV, a virus with a large dsDNA genome, exists as a complex mixture of genome types in humans and offer insights into the evolution of the virus.     

Extensive amino acid polymorphism at the pgm locus is consistent with adaptive protein evolution in Drosophila melanogaster

PGM plays a central role in the glycolytic pathway at the branch point leading to glycogen metabolism and is highly polymorphic in allozyme studies of many species. We have characterized the nucleotide diversity across the Pgm gene in Drosophila melanogaster and D. simulans to investigate the role that protein polymorphism plays at this crucial metabolic branch point shared with several other enzymes. Although D. melanogaster and D. simulans share common allozyme mobility alleles, we find these allozymes are the result of many different amino acid changes at the nucleotide level. In addition, specific allozyme classes within species contain several amino acid changes, which may explain the absence of latitudinal clines for PGM allozyme alleles, the lack of association of PGM allozymes with the cosmopolitan In(3L)P inversion, and the failure to detect differences between PGM allozymes in functional studies. We find a significant excess of amino acid polymorphisms within D. melanogaster when compared to the complete absence of fixed replacements with D. simulans. There is also strong linkage disequilibrium across the 2354 bp of the Pgm locus, which may be explained by a specific amino acid haplotype that is high in frequency yet contains an excess of singleton polymorphisms. Like G6pd, Pgm shows strong evidence for a branch point enzyme that exhibits adaptive protein evolution.     

Selection on amino acid substitutions in Arabidopsis

Studies of nucleotide diversity have found an excess of low-frequency amino acid polymorphisms segregating in Arabidopsis thaliana, suggesting a predominance of weak purifying selection acting on amino acid polymorphism in this inbreeding species. Here, we investigate levels of diversity and divergence at synonymous and nonsynonymous sites in 6 circumpolar populations of the outbreeding Arabidopsis lyrata and compare these results with A. thaliana, to test for differences in mutation and selection parameters across genes, populations, and species. We find that A. lyrata shows an excess of low-frequency nonsynonymous polymorphisms both within populations and species wide, consistent with weak purifying selection similar to the patterns observed in A. thaliana. Furthermore, nonsynonymous polymorphisms tend to be more restricted in their population distribution in A. lyrata, consistent with purifying selection preventing their geographic spread. Highly expressed genes show a reduced ratio of amino acid to synonymous change for both polymorphism and fixed differences, suggesting a general pattern of stronger purifying selection on high-expression proteins.     

Evidence of amino acid diversity-enhancing selection within humans and among primates at the candidate sperm-receptor gene PKDREJ

Sperm-egg interaction is a crucial step in fertilization, yet the identity of most interacting sperm-egg proteins that mediate this process remains elusive. Rapid evolution of some fertilization proteins has been observed in a number of species, including evidence of positive selection in the evolution of components of the mammalian egg coat. The rapid evolution of the egg-coat proteins could strongly select for changes on the sperm receptor, to maintain the interaction. Here, we present evidence that positive selection has driven the evolution of PKDREJ, a candidate sperm receptor of mammalian egg-coat proteins. We sequenced PKDREJ from a panel of 14 primates, including humans, and conducted a comparative maximum-likelihood analysis of nucleotide changes and found evidence of positive selection. An additional panel of 48 humans was surveyed for nucleotide polymorphisms at the PKDREJ locus. The regions predicted to have been subject to adaptive evolution among primates show several amino acid polymorphisms within humans. The distribution of polymorphisms suggests that balancing selection may maintain diverse PKDREJ alleles in some populations. It remains unknown whether there are functional differences associated with these diverse alleles, but their existence could have consequences for human fertility.     

Clinal variation for amino acid polymorphisms at the Pgm locus in Drosophila melanogaster

Clinal variation is common for enzymes in the glycolytic pathway for Drosophila melanogaster and is generally accepted as an adaptive response to different climates. Although the enzyme phosphoglucomutase (PGM) possesses several allozyme polymorphisms, it is unique in that it had been reported to show no clinal variation. Our recent DNA sequence investigation of Pgm found extensive cryptic amino acid polymorphism segregating with the allozyme alleles. In this study, we characterize the geographic variation of Pgm amino acid polymorphisms at the nucleotide level along a latitudinal cline in the eastern United States. A survey of 15 SNPs across the Pgm gene finds significant clinal differentiation for the allozyme polymorphisms as well as for many of the cryptic amino acid polymorphisms. A test of independence shows that pervasive linkage disequilibrium across this gene region can explain many of the amino acid clines. A single Pgm haplotype defined by two amino acid polymorphisms shows the strongest correlation with latitude and the steepest change in allele frequency across the cline. We propose that clinal selection at Pgm may in part explain the extensive amino acid polymorphism at this locus and is consistent with a multilocus response to selection in the glycolytic pathway.     

Complex selection on a regulator of social cognition: Evidence of balancing selection,regulatory interactions and population differentiation in the prairie vole Avpr1a locus

Adaptive variation in social behaviour depends upon standing genetic variation, but we know little about how evolutionary forces shape genetic diversity relevant to brain and behaviour. In prairie voles (Microtus ochrogaster), variants at the Avpr1a locus predict expression of the vasopressin 1a receptor in the retrosplenial cortex (RSC), a brain region that mediates spatial and contextual memory; cortical V1aR abundance in turn predicts diversity in space use and sexual fidelity in the field. To examine the potential contributions of adaptive and neutral forces to variation at the Avpr1a locus, we explore sequence diversity at the Avpr1a locus and throughout the genome in two populations of wild prairie voles. First, we refine results demonstrating balancing selection at the locus by comparing the frequency spectrum of variants at the locus to a random sample of the genome. Next, we find that the four single nucleotide polymorphisms that predict high V1aR expression in the RSC are in stronger linkage disequilibrium than expected by chance despite high recombination among intervening variants, suggesting that epistatic selection maintains their association. Analysis of population structure and a haplotype network for two populations revealed that this excessive LD was unlikely to be due to admixture alone. Furthermore, the two populations differed considerably in the region shown to be a regulator of V1aR expression despite the extremely low levels of genomewide genetic differentiation. Together, our data suggest that complex selection on Avpr1a locus favours specific combinations of regulatory polymorphisms, maintains the resulting alleles at population‐specific frequencies, and may contribute to unique patterns of spatial cognition and sexual fidelity among populations.     

African Glucose-6-Phosphate Dehydrogenase Alleles Associated with Protection from Severe Malaria in Heterozygous Females in Tanzania

X-linked Glucose-6-phosphate dehydrogenase (G6PD) A- deficiency is prevalent in sub-Saharan Africa populations, and has been associated with protection from severe malaria. Whether females and/or males are protected by G6PD deficiency is uncertain, due in part to G6PD and malaria phenotypic complexity and misclassification. Almost all large association studies have genotyped a limited number of G6PD SNPs (e.g. G6PD202 / G6PD376), and this approach has been too blunt to capture the complete epidemiological picture. Here we have identified 68 G6PD polymorphisms and analysed 29 of these (i.e. those with a minor allele frequency greater than 1%) in 983 severe malaria cases and controls in Tanzania. We establish, across a number of SNPs including G6PD376, that only female heterozygotes are protected from severe malaria. Haplotype analysis reveals the G6PD locus to be under balancing selection, suggesting a mechanism of protection relying on alleles at modest frequency and avoiding fixation, where protection provided by G6PD deficiency against severe malaria is offset by increased risk of life-threatening complications. Our study also demonstrates that the much-needed large-scale studies of severe malaria and G6PD enzymatic function across African populations require the identification and analysis of the full repertoire of G6PD genetic markers.     

Parallel introgression and selection on introduced alleles in a native species

As humans cause the redistribution of species ranges, hybridization between previously allopatric species is on the rise. Such hybridization can have complex effects on overall fitness of native species as new allelic combinations are tested. Widespread species introductions provide a unique opportunity to study selection on introgressed alleles in independent, replicated populations. We examined selection on alleles that repeatedly introgressed from introduced rainbow trout (Oncorhynchus mykiss) into native westslope cutthroat trout (Oncorhynchus clarkii lewisi) populations in western Canada. We found that the degree of introgression of individual single nucleotide polymorphisms from the invasive species into the native is correlated between independent watersheds. A number of rainbow trout alleles have repeatedly swept to high frequency in native populations, suggesting parallel adaptive advantages. Using simulations, we estimated large selection coefficients up to 0.05 favoring several rainbow trout alleles in the native background. Although previous studies have found reduced hybrid fitness and genome‐wide resistance to introgression in westslope cutthroat trout, our results suggest that some introduced genomic regions are strongly favored by selection. Our study demonstrates the utility of replicated introductions as case studies for understanding parallel adaptation and the interactions between selection and introgression across the genome. We suggest that understanding this variation, including consideration of beneficial alleles, can inform management strategies for hybridizing species.     

A New Glucose-6-Phosphate Dehydrogenase Variant, G6PD Orissa (44 Ala→Gly), is the Major Polymorphic Variant in Tribal Populations in India

Deficiency of glucose-6-phosphate dehydrogenase (G6PD) is usually found at high frequencies in areas of the world where malaria has been endemic. The frequency and genetic basis of G6PD deficiency have been studied in Africa, around the Mediterranean, and in the Far East, but little such information is available about the situation in India. To determine the extent of heterogeneity of G6PD, we have studied several different Indian populations by screening for G6PD deficiency, followed by molecular analysis of deficient alleles. The frequency of G6PD deficiency varies between 3% and 15% in different tribal and urban groups. Remarkably, a previously unreported deficient variant, G6PD Orissa (44 Ala→Gly), is responsible for most of the G6PD deficiency in tribal Indian populations but is not found in urban populations, where most of the G6PD deficiency is due to the G6PD Mediterranean (188 Ser→Phe) variant. The K of G6PD Orissa is fivefold higher than that of the normal enzyme. This may be due to the fact that the alanine residue that is replaced by glycine is part of a putative coenzyme-binding site.     

Multiple glucose 6-phosphate dehydrogenase-deficient variants correlate with malaria endemicity in the Vanuatu archipelago (southwestern Pacific).

In studying the relationship between genetic abnormalities of red blood cells and malaria endemicity in the Vanuatu archipelago in the southwestern Pacific, we have found that of 1,442 males tested, 98 (6.8%) were G6PD deficient. The prevalence of GdPD deficiency varied widely (0%-39%), both from one island to another and in different parts of the same island, and generally correlated positively with the degree of malaria transmission. The properties of G6PD from GdPD-deficient subjects were analyzed in a subset of 53 samples. In all cases the residual red-blood-cell activity was < 10%. There were three phenotypic patterns. PCR amplification and sequencing of the entire coding region of the G6PD gene showed that the first of these patterns corresponded to G6PD Union (nucleotide 1360C-->T; amino acid 454Arg-->Cys), previously encountered elsewhere. Analysis of samples exhibiting the second pattern revealed two new mutants: G6PD Vanua Lava (nucleotide 383T-->C; amino acid 128Leu-->Pro) and G6PD Namoru (nucleotide 208T-->C; amino acid 70Tyr-->His); in three samples, the underlying mutation has not yet been identified. Analysis of the sample exhibiting the third pattern revealed another new mutant: G6PD Naone (nucleotide 497G-->A; amino acid 166Arg-->His). Of the four mutations, G6PD Union and G6PD Vanua Lava have a polymorphic frequency in more than one island; and G6PD Vanua Lava has also been detected in a sample from Papua New Guinea. G6PD deficiency is of clinical importance in Vanuatu because it is a cause of neonatal jaundice and is responsible for numerous episodes of drug-induced acute hemolytic anemia.     

Modeling genetic imprinting effects of DNA sequences with multilocus polymorphism data

Single nucleotide polymorphisms (SNPs) represent the most widespread type of DNA sequence variation in the human genome and they have recently emerged as valuable genetic markers for revealing the genetic architecture of complex traits in terms of nucleotide combination and sequence. Here, we extend an algorithmic model for the haplotype analysis of SNPs to estimate the effects of genetic imprinting expressed at the DNA sequence level. The model provides a general procedure for identifying the number and types of optimal DNA sequence variants that are expressed differently due to their parental origin. The model is used to analyze a genetic data set collected from a pain genetics project. We find that DNA haplotype GAC from three SNPs, OPRKG36T (with two alleles G and T), OPRKA843G (with alleles A and G), and OPRKC846T (with alleles C and T), at the kappa-opioid receptor, triggers a significant effect on pain sensitivity, but with expression significantly depending on the parent from which it is inherited (p = 0.008). With a tremendous advance in SNP identification and automated screening, the model founded on haplotype discovery and statistical inference may provide a useful tool for genetic analysis of any quantitative trait with complex inheritance.     

Molecular genetics of the glucose-6-phosphate dehydrogenase (G6PD) Mediterranean variant and description of a new G6PD mutant, G6PD Andalus1361A

Glucose-6-phosphate dehydrogenase (G6PD; E.C.1.1.1.49) deficiency is the most common human enzymopathy; more than 300 different biochemical variants of the enzyme have been described. In many parts of the world the Mediterranean type of G6PD deficiency is prevalent. However, G6PD Mediterranean has come to be regarded as a generic term applied to similar G6PD mutations thought, however, to represent a somewhat heterogeneous group. A C----T mutation at nucleotide 563 of G6PD Mediterranean has been identified by Vulliamy et al., and the same mutation has been found by De Vita et al. in G6PD Mediterranean, G6PD Sassari, and G6PD Cagliari. The latter subjects had an additional mutation, at nucleotide 1311, that did not produce a coding change. We have examined genomic DNA of five patients--four of Spanish origin and one of Jewish origin--having enzymatically documented G6PD Mediterranean. All had both the mutation at nucleotide 563 and that at nucleotide 1311. A sixth sample, resembling G6PD Mediterranean kinetically but with a slightly rapid electrophoretic mobility, was designated G6PD Andalus and was found to have a different mutation, a G----A transition at nucleotide 1361, producing an arginine-to-histidine substitution. These studies suggest that G6PD Mediterranean is, after all, relatively homogeneous.     

Targeted disruption of the housekeeping gene encoding glucose 6-phosphate dehydrogenase (G6PD): G6PD is dispensable for pentose synthesis but essential for defense against oxidative stress.

Glucose 6-phosphate dehydrogenase (G6PD) is a housekeeping enzyme encoded in mammals by an X-linked gene. It has important functions in intermediary metabolism because it catalyzes the first step in the pentose phosphate pathway and provides reductive potential in the form of NADPH. In human populations, many mutant G6PD alleles (some present at polymorphic frequencies) cause a partial loss of G6PD activity and a variety of hemolytic anemias, which vary from mild to severe. All these mutants have some residual enzyme activity, and no large deletions in the G6PD gene have ever been found. To test which, if any, function of G6PD is essential, we have disrupted the G6PD gene in male mouse embryonic stem cells by targeted homologous recombination. We have isolated numerous clones, shown to be recombinant by Southern blot analysis, in which G6PD activity is undetectable. We have extensively characterized individual clones and found that they are extremely sensitive to H2O2 and to the sulfydryl group oxidizing agent, diamide. Their markedly impaired cloning efficiency is restored by reducing the oxygen tension. We conclude that G6PD activity is dispensable for pentose synthesis, but is essential to protect cells against even mild oxidative stress.     

Isolation of human glucose-6-phosphate dehydrogenase (G6PD) cDNA clones: primary structure of the protein and unusual 5'' non-coding region.

Glucose-6-phosphate dehydrogenase (G6PD) is an ubiquitous enzyme which by determining the NADPH level has a crucial role in NADPH-mediated reductive processes in all cells (1). The structural gene for G6PD, Gd, is X-linked in mammals and on the basis of its expression in many tissues, it can be regarded as a typical "housekeeping" gene (2). Over 300 variants of the protein are known, many of which have deficient enzyme activity. Nearly 100 of these variants are polymorphic in various populations (3). The mammalian enzyme is a homodimer or a homotetramer with a subunit molecular weight of approximately 56000 daltons (4). Here we report the isolation of cDNA clones from HeLa cells, SV40-transformed human fibroblasts, human placenta and human teratocarcinoma cell lines. These clones have enabled us to sequence the entire coding region of Gd. Thus, the entire amino acid sequence of human G6PD is provided for the first time. This work is the first step for structural analysis of G6PD variants and for an understanding of the biological features of this enzyme at the molecular level.     

Human mutations in glucose 6-phosphate dehydrogenase reflect evolutionary history.

Glucose 6-phosphate dehydrogenase (G6PD) is a cytosolic enzyme encoded by a housekeeping X-linked gene whose main function is to produce NADPH, a key electron donor in the defense against oxidizing agents and in reductive biosynthetic reactions. Inherited G6PD deficiency is associated with either episodic hemolytic anemia (triggered by fava beans or other agents) or life-long hemolytic anemia. We show here that an evolutionary analysis is a key to understanding the biology of a housekeeping gene. From the alignment of the amino acid (aa) sequence of 52 glucose 6-phosphate dehydrogenase (G6PD) species from 42 different organisms, we found a striking correlation between the aa replacements that cause G6PD deficiency in humans and the sequence conservation of G6PD: two-thirds of such replacements are in highly and moderately conserved (50-99%) aa; relatively few are in fully conserved aa (where they might be lethal) or in poorly conserved aa, where presumably they simply would not cause G6PD deficiency. This is consistent with the notion that all human mutants have residual enzyme activity and that null mutations are lethal at some stage of development. Comparing the distribution of mutations in a human housekeeping gene with evolutionary conservation is a useful tool for pinpointing amino acid residues important for the stability or the function of the corresponding protein. In view of the current explosive increase in full genome sequencing projects, this tool will become rapidly available for numerous other genes.