Haplotype diversity and linkage disequilibrium at human G6PD: recent origin of alleles that confer malarial resistance

Host preference patterns in domestic and wild settings: Insights into Anopheles feeding behavior

The adaptation of Anopheles malaria vectors to domestic settings is directly linked to their ability to feed on humans. The strength of this species-habitat association is unequal across the species within the genus, with the major vectors being particularly dependent on humans. However, our understanding of how blood-feeding behavior interacts with and adapts to environmental settings, including the presence of humans, remains limited. Using a field-based approach, we first investigated Anopheles community structure and feeding behavior patterns in domestic and sylvatic settings in La Lopé National Park in Gabon, Central Africa. We characterized the preference indices using a dual-host choice sampling approach across mosquito species, habitats, and seasons. We then quantified the plastic biting behavior of mosquito species in each habitat. We collected individuals from 16 Anopheles species that exhibited significant differences in species composition and abundance between sylvatic and domestic settings. The host-seeking behavior also varied among the seven most abundant species. The general attractiveness to each host, human or animal, remained relatively constant for each species, but with significant variations between habitats across species. These variations, to more generalist and to more anthropophilic behavior, were related to seasonal changes and distance from the village, respectively. Finally, we pointed out that the host choice of major malaria vectors changed in the absence of humans, revealing a plastic feeding behavior of these species. This study highlights the effect of humans on Anopheles distribution and feeding evolution. The characterization of feeding behavior in wild and domestic settings provides opportunities to better understand the interplay between genetic determinants of host preference and ecological factors. Our findings suggest that protected areas may offer alternative thriving conditions to major malaria vectors.

scAI-SNP: a method for inferring ancestry from single-cell data

Collaborative efforts, such as the Human Cell Atlas, are rapidly accumulating large amounts of single-cell data. To ensure that single-cell atlases are representative of human genetic diversity, we need to determine the ancestry of the donors from whom single-cell data are generated. Self-reporting of race and ethnicity, although important, can be biased and is not always available for the datasets already collected. Here, we introduce scAI-SNP, a tool to infer ancestry directly from single-cell genomics data. To train scAI-SNP, we identified 4.5 million ancestry-informative single-nucleotide polymorphisms (SNPs) in the 1000 Genomes Project dataset across 3201 individuals from 26 population groups. For a query single-cell data set, scAI-SNP uses these ancestry-informative SNPs to compute the contribution of each of the 26 population groups to the ancestry of the donor from whom the cells were obtained. Using diverse single-cell data sets with matched whole-genome sequencing data, we show that scAI-SNP is robust to the sparsity of single-cell data, can accurately and consistently infer ancestry from samples derived from diverse types of tissues and cancer cells, and can be applied to different modalities of single-cell profiling assays, such as single-cell RNA-seq and single-cell ATAC-seq. Finally, we argue that ensuring that single-cell atlases represent diverse ancestry, ideally alongside race and ethnicity, is ultimately important for improved and equitable health outcomes by accounting for human diversity.

Ancient genomes illuminate Eastern Arabian population history and adaptation against malaria

The harsh climate of Arabia has posed challenges in generating ancient DNA from the region, hindering the direct examination of ancient genomes for understanding the demographic processes that shaped Arabian populations. In this study, we report whole-genome sequence data obtained from four Tylos-period individuals from Bahrain. Their genetic ancestry can be modeled as a mixture of sources from ancient Anatolia, Levant, and Iran/Caucasus, with variation between individuals suggesting population heterogeneity in Bahrain before the onset of Islam. We identify the G6PD Mediterranean mutation associated with malaria resistance in three out of four ancient Bahraini samples and estimate that it rose in frequency in Eastern Arabia from 5 to 6 kya onward, around the time agriculture appeared in the region. Our study characterizes the genetic composition of ancient Arabians, shedding light on the population history of Bahrain and demonstrating the feasibility of studies of ancient DNA in the region.

Genetic basis and spatial distribution of glucose-6-phosphate dehydrogenase deficiency in ecuadorian ethnic groups: a malaria perspective

BACKGROUND

Glucose-6-phosphate dehydrogenase deficiency (G6PDd) is an X-linked disorder affecting over 400 million people worldwide. Individuals with molecular variants associated with reduced enzymatic activity are susceptible to oxidative stress in red blood cells, thereby increasing the risk of pathophysiological conditions and toxicity to anti-malarial treatments. Globally, the prevalence of G6PDd varies among populations. Accordingly, this study aims to characterize G6PDd distribution within the Ecuadorian population and to describe the spatial distribution of reported malaria cases.

METHODS

Molecular variants associated with G6PDd were genotyped in 581 individuals from Afro-Ecuadorian, Indigenous, Mestizo, and Montubio ethnic groups. Additionally, spatial analysis was conducted to identify significant malaria clusters with high incidence rates across Ecuador, using data collected from 2010 to 2021.

RESULTS

The A- c.202G > A and A- c.968T > C variants underpin the genetic basis of G6PDd in the studied population. The overall prevalence of G6PDd was 4.6% in the entire population. However, this frequency increased to 19.2% among Afro-Ecuadorian people. Spatial analysis revealed 12 malaria clusters, primarily located in the north of the country and its Amazon region, with relative risks of infection of 2.02 to 87.88.

CONCLUSIONS

The findings of this study hold significant implications for public health interventions, treatment strategies, and targeted efforts to mitigate the burden of malaria in Ecuador. The high prevalence of G6PDd among Afro-Ecuadorian groups in the northern endemic areas necessitates the development of comprehensive malaria eradication strategies tailored to this geographical region.

Glucose-6-phosphate dehydrogenase deficiency and long-term risk of immune-related disorders

Introduction

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is an X-linked enzymatic disorder that is particularly prevalent in Africa, Asia, and the Middle East. This study aimed to assess the long-term health risks associated with G6PD deficiency.

Methods

A retrospective cohort study was conducted using data from a national healthcare provider in Israel (Leumit Health Services). A total of 7,473 G6PD-deficient individuals were matched with 29,892 control subjects in a 1:4 ratio, based on age, gender, socioeconomic status, and ethnic groups. The exposure of interest was recorded G6PD diagnosis or positive G6PD diagnostic test. The main outcomes and measures included rates of infectious diseases, allergic conditions, and autoimmune disorders between 2002 and 2022.

Results

Significantly increased rates were observed for autoimmune disorders, infectious diseases, and allergic conditions in G6PD-deficient individuals compared to the control group. Specifically, notable increases were observed for rheumatoid arthritis (odds ratio [OR] 2.41, p<0.001), systemic lupus erythematosus (OR 4.56, p<0.001), scleroderma (OR 6.87, p<0.001), pernicious anemia (OR 18.70, p<0.001), fibromyalgia (OR 1.98, p<0.001), Graves' disease (OR 1.46, p=0.001), and Hashimoto's thyroiditis (OR 1.26, p=0.001). These findings were supported by elevated rates of positive autoimmune serology and higher utilization of medications commonly used to treat autoimmune conditions in the G6PD-deficient group.

Discussion

In conclusion, individuals with G6PD deficiency are at a higher risk of developing autoimmune disorders, infectious diseases, and allergic conditions. This large-scale observational study provides valuable insights into the comprehensive association between G6PD deficiency and infectious and immune-related diseases. The findings emphasize the importance of considering G6PD deficiency as a potential risk factor in clinical practice and further research is warranted to better understand the underlying mechanisms of these associations.

Hunter-gatherer diets and activity as a model for health promotion: Challenges, responses, and confirmations

Beginning in 1985, we and others presented estimates of hunter-gatherer (and ultimately ancestral) diet and physical activity, hoping to provide a model for health promotion. The Hunter-Gatherer Model was designed to offset the apparent mismatch between our genes and the current Western-type lifestyle, a mismatch that arguably affects prevalence of many chronic degenerative diseases. The effort has always been controversial and subject to both scientific and popular critiques. The present article (1) addresses eight such challenges, presenting for each how the model has been modified in response, or how the criticism can be rebutted; (2) reviews new epidemiological and experimental evidence (including especially randomized controlled clinical trials); and (3) shows how official recommendations put forth by governments and health authorities have converged toward the model. Such convergence suggests that evolutionary anthropology can make significant contributions to human health.

Spatial and temporal diversity of positive selection on shared haplotypes at the PSCA locus among worldwide human populations

Selection on standing genetic variation is important for rapid local genetic adaptation when the environment changes. We report that, for the prostate stem cell antigen (PSCA) gene, different populations have different target haplotypes, even though haplotypes are shared among populations. The C-C-A haplotype, whereby the first C is located at rs2294008 of PSCA and is a low risk allele for gastric cancer, has become a target of positive selection in Asia. Conversely, the C-A-G haplotype carrying the same C allele has become a selection target mainly in Africa. However, Asian and African share both haplotypes, consistent with the haplotype divergence time (170 kya) prior to the out-of-Africa dispersal. The frequency of C-C-A/C-A-G is 0.344/0.278 in Asia and 0.209/0.416 in Africa. Two-dimensional site frequency spectrum analysis revealed that the extent of intra-allelic variability of the target haplotype is extremely small in each local population, suggesting that C-C-A or C-A-G is under ongoing hard sweeps in local populations. From the time to the most recent common ancestor (TMRCA) of selected haplotypes, the onset times of positive selection were recent (3-55 kya), concurrently with population subdivision from a common ancestor. Additionally, estimated selection coefficients from ABC analysis were up to ~3%, similar to those at other loci under recent positive selection. Phylogeny of local populations and TMRCA of selected haplotypes revealed that spatial and temporal switching of positive selection targets is a unique and novel feature of ongoing selection at PSCA. This switching may reflect the potential of rapid adaptability to distinct environments.

The adaptive evolution of cancer driver genes

BACKGROUND

Cancer is a life-threatening disease in humans; yet, cancer genes are frequently reported to be under positive selection. This suggests an evolutionary-genetic paradox in which cancer evolves as a secondary product of selection in human beings. However, systematic investigation of the evolution of cancer driver genes is sparse.

RESULTS

Using comparative genomics analysis, population genetics analysis and computational molecular evolutionary analysis, the evolution of 568 cancer driver genes of 66 cancer types were evaluated at two levels, selection on the early evolution of humans (long timescale selection in the human lineage during primate evolution, i.e., millions of years), and recent selection in modern human populations (~ 100,000 years). Results showed that eight cancer genes covering 11 cancer types were under positive selection in the human lineage (long timescale selection). And 35 cancer genes covering 47 cancer types were under positive selection in modern human populations (recent selection). Moreover, SNPs associated with thyroid cancer in three thyroid cancer driver genes (CUX1, HERC2 and RGPD3) were under positive selection in East Asian and European populations, consistent with the high incidence of thyroid cancer in these populations.

CONCLUSIONS

These findings suggest that cancer can be evolved, in part, as a by-product of adaptive changes in humans. Different SNPs at the same locus can be under different selection pressures in different populations, and thus should be under consideration during precision medicine, especially for targeted medicine in specific populations.

Glucose-6-Phosphate Dehydrogenase (G6PD) Deficiency and Long-Term Risk of Immune-Related diseases

BACKGROUND

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is an X-linked recessive enzymatic disorder, particularly prevalent in Africa, Asia and the Middle East. In the US, about 14% of black men are affected. Individuals with G6PD deficiency are often asymptomatic but may develop hemolysis following an infection or upon consumption of specific medications. Despite some evidence that G6PD deficiency affects the immune system, the long- term health risks associated with G6PD deficiency had not been studied in a large population.

METHODS

In this retrospective cohort study, health records from G6PD deficient individuals were compared to matched controls in a national healthcare provider in Israel (Leumit Health Services). Rates of infectious diseases, allergic conditions and autoimmune disorders were compared between groups.

RESULTS

The cohort included 7,473 G6PD deficient subjects (68.7% men) matched with 29,892 control subjects (4:1 ratio) of the same age, gender, socioeconomic status and ethnic group, followed during 14.3±6.2 years.Significantly increased rates for autoimmune disorders, infectious diseases and allergic conditions were observed throughout this period. Notable increases were observed for rheumatoid arthritis (OR 2.41, p<0.001), systemic lupus erythematosus (OR 4.56, p<0.001), scleroderma (OR 6.87, p<0.001), pernicious anemia (OR=18.70, P<0.001), fibromyalgia (OR 1.98, p<0.001), Graves' disease (OR 1.46, P=0.001), and Hashimoto's thyroiditis (OR 1.26, P=0.001). These findings were corroborated with elevated rates of positive autoimmune serology and higher rates of treatment with medications commonly used to treat autoimmune conditions in the G6PD deficient group.

CONCLUSION

G6PD deficient individuals suffer from higher rates of autoimmune disorders, infectious diseases, and allergic conditions.

Using epigenomics to understand cellular responses to environmental influences in diseases

It is a generally accepted model that environmental influences can exert their effects, at least in part, by changing the molecular regulators of transcription that are described as epigenetic. As there is biochemical evidence that some epigenetic regulators of transcription can maintain their states long term and through cell division, an epigenetic model encompasses the idea of maintenance of the effect of an exposure long after it is no longer present. The evidence supporting this model is mostly from the observation of alterations of molecular regulators of transcription following exposures. With the understanding that the interpretation of these associations is more complex than originally recognised, this model may be oversimplistic; therefore, adopting novel perspectives and experimental approaches when examining how environmental exposures are linked to phenotypes may prove worthwhile. In this review, we have chosen to use the example of nonalcoholic fatty liver disease (NAFLD), a common, complex human disease with strong environmental and genetic influences. We describe how epigenomic approaches combined with emerging functional genetic and single-cell genomic techniques are poised to generate new insights into the pathogenesis of environmentally influenced human disease phenotypes exemplified by NAFLD.

Parallel evolution in human populations: A biocultural perspective

Parallel evolution-where different populations evolve similar traits in response to similar environments-has been a topic of growing interest to biologists and biological anthropologists for decades. Parallel evolution occurs in human populations thanks to myriad biological and cultural mechanisms that permit humans to survive and thrive in diverse environments worldwide. Because humans shape and are shaped by their environments, biocultural approaches that emphasize the interconnections between biology and culture are key to understanding parallel evolution in human populations as well as the nuances of human biological variation and adaptation. In this review, we discuss how biocultural theory has been and can be applied to studies of parallel evolution and adaptation more broadly. We illustrate this through four examples of parallel evolution in humans: malaria resistance, lactase persistence, cold tolerance, and high-altitude adaptation.

Diverse selection pressures shaping the genetic architecture of behçet disease susceptibility

Behçet disease (BD) is a polygenic, multifactorial, multisystem inflammatory condition with unknown etiology. Global distribution of BD is geographically structured, highest prevalence observed among East Asian, Middle Eastern, and Mediterranean populations. Although adaptive selection on a few BD susceptibility loci is speculated, a thorough evolutionary analysis on the genetic architecture of BD is lacking. We aimed to understand whether increased BD risk in the human populations with high prevalence is due to past selection on BD associated genes. We performed population genetics analyses with East Asian (high BD prevalence), European (low/very low BD prevalence), and African (very low/no BD prevalence) populations. Comparison of ancestral and derived alleles’ frequencies versus their reported susceptible or protective effect on BD showed both derived and ancestral alleles are associated with increased BD risk. Variants showing higher risk to and more significant association with BD had smaller allele frequency differences, and showed less population differentiation compared to variants that showed smaller risk and less significant association with BD. Results suggest BD alleles are not unique to East Asians but are also found in other world populations at appreciable frequencies, and argue against selection favoring these variants only in populations with high BD prevalence. BD associated gene analyses showed similar evolutionary histories driven by neutral processes for many genes or balancing selection for HLA (Human Leukocyte Antigen) genes in all three populations studied. However, nucleotide diversity in several HLA region genes was much higher in East Asians suggesting selection for high nucleotide and haplotype diversity in East Asians. Recent selective sweep for genes involved in antigen recognition, peptide processing, immune and cellular differentiation regulation was observed only in East Asians. We conclude that the evolutionary processes shaping the genetic diversity in BD risk genes are diverse, and elucidating the underlying specific selection mechanisms is complex. Several of the genes examined in this study are risk factors (such as ERAP1, IL23R, HLA-G) for other inflammatory diseases. Thus, our conclusions are not only limited to BD but may have broader implications for other inflammatory diseases.

Uncovering the extensive trade-off between adaptive evolution and disease susceptibility

Favored mutations in the human genome may make the carriers adapt to changing environments and lifestyles but also susceptible to specific diseases. The scale and details of the trade-off between adaptive evolution and disease susceptibility are unclear because most favored mutations in different populations remain unidentified. As no statistical test can discriminate favored mutations from nearby hitchhiking neutral ones, we report a deep-learning network (DeepFavored) to integrate multiple statistical tests and divide identifying favored mutations into two subtasks. We identify favored mutations in three human populations and analyzed the correlation between favored/hitchhiking mutations and genome-wide association study (GWAS) sites. Both favored and hitchhiking neutral mutations are enriched in GWAS sites with population-specific features, and the enrichment and population specificity are prominent in genes in specific Gene Ontology (GO) terms. These provide evidence for extensive and population-specific trade-offs between adaptive evolution and disease susceptibility. The unveiled scale helps understand and investigate differences and diseases of humans.

Synthetic Red Blood Cell-Specific Glycolytic Intermediate 2,3-Diphosphoglycerate (2,3-DPG) Inhibits Plasmodium falciparum Development In Vitro

Mechanisms of malaria parasite interaction with its host red blood cell may provide potential targets for new antimalarial approaches. Pyruvate kinase deficiency has been associated with resistance to malaria in both experimental models and population studies. Two of the major pyruvate kinase deficient-cell disorders are the decrease in ATP and the increase in 2,3-biphosphoglycerate (2,3-BPG) concentration. High levels of this metabolite, only present in mammalian red blood cell, has an inhibitory effect on glycolysis and we hypothesized that its accumulation may also be harmful to the parasite and be involved in the mechanism of protection provided by that enzymopathy. We examined the effect of a synthetic form, 2,3-DPG, on the Plasmodium falciparum intraerythrocytic developmental cycle in vitro. Results showed an impairment of parasite growth with a direct effect on parasite maturation as significant lower progeny emerged from parasites that were submitted to 2,3-DPG. Further, adding the compound to the culture medium did not result in any effect on the host cell, but instead the metabolic profile of an infected cell became closer to that of a non-infected cell.

Prevalence of glucose-6-phosphate dehydrogenase deficiency (G6PDd), CareStart qualitative rapid diagnostic test performance, and genetic variants in two malaria-endemic areas in Sudan

Glucose-6-phosphate dehydrogenase deficiency (G6PDd) is the most common enzymopathy globally, and deficient individuals may experience severe hemolysis following treatment with 8-aminoquinolines. With increasing evidence of Plasmodium vivax infections throughout sub-Saharan Africa, there is a pressing need for population-level data at on the prevalence of G6PDd. Such evidence-based data will guide the expansion of primaquine and potentially tafenoquine for radical cure of P. vivax infections. This study aimed to quantify G6PDd prevalence in two geographically distinct areas in Sudan, and evaluating the performance of a qualitative CareStart rapid diagnostic test as a point-of-care test. Blood samples were analyzed from 491 unrelated healthy persons in two malaria-endemic sites in eastern and central Sudan. A pre-structured questionnaire was used which included demographic data, risk factors and treatment history. G6PD levels were measured using spectrophotometry (SPINREACT) and first-generation qualitative CareStart rapid tests. G6PD variants (202 G>A; 376 A>G) were determined by PCR/RFLP, with a subset confirmed by Sanger sequencing. The prevalence of G6PDd by spectrophotometry was 5.5% (27/491; at 30% of adjusted male median, AMM); 27.3% (134/491; at 70% of AMM); and 13.1% (64/490) by qualitative CareStart rapid diagnostic test. The first-generation CareStart rapid diagnostic test had an overall sensitivity of 81.5% (95%CI: 61.9 to 93.7) and negative predictive value of 98.8% (97.3 to 99.6). All persons genotyped across both study sites were wild type for the G6PD G202 variant. For G6PD A376G all participants in New Halfa had wild type AA (100%), while in Khartoum the AA polymorphism was found in 90.7%; AG in 2.5%; and GG in 6.8%. Phenotypic G6PD B was detected in 100% of tested participants in New Halfa while in Khartoum, the phenotypes observed were B (96.2%), A (2.8%), and AB (1%). The African A- phenotype was not detected in this study population. Overall, G6PDd prevalence in Sudan is low-to-moderate but highly heterogeneous. Point-of-care testing with the qualitative CareStart rapid diagnostic test demonstrated moderate performance with moderate sensitivity and specificity but high negative predicative value. The two sites harbored primarily the African B phenotype. A country-wide survey is recommended to understand GP6PD deficiencies more comprehensively in Sudan.

Malaria is caused by five species of parasites; of these Plasmodium falciparum and P. vivax cause the majority of global morbidity and mortality. Plasmodium vivax infection is an emerging public health problem in sub-Saharan Africa, including Sudan. Primaquine and other 8-aminoquinolines including tafenoquine are the primary treatments to target the silent liver stage (hypnozoites) in P. vivax infections. However, these regimens can cause severe intravascular hemolysis in patients suffering from glucose-6-phosphate dehydrogenase deficiency (G6PDd). To support safe and efficacious use of primaquine, and potentially tafenoquine in Sudan, this study aimed to estimate the prevalence of G6PDd across two sites in Sudan using spectrophotometry and a qualitative CareStart rapid diagnostic test. Subsequent genetic analysis by PCR/RFLP and sequencing of G6PD genetic variants was performed. This survey found an overall prevalence was 5.5% (27/491; 30% of adjusted male median, AMM), and 27.3% (134/491; 70% of AMM) and 13.1% (64/490) by qualitative CareStart rapid diagnostic test. Important differences in distribution of genetic variants of G6PD were found across the two sites, and the African A- was not observed. In univariate analysis a few parameters showed significant association with G6PD deficiency. In conclusion the prevalence of G6PDd was low to moderate but heterogonous, and the first-generation qualitative CareStart rapid diagnostic test showed moderate performance in both males and females.

Natural selection contributes to the myopia epidemic

The prevalence of myopia, or nearsightedness, has skyrocketed in the past few decades, creating a public health crisis that is commonly attributed to lifestyle changes. Here we report an overall increase in the frequencies of myopia-associated mutant alleles over 25 years among participants of the UK Biobank. Although myopia itself appears to be selected against, many of the mutant alleles are associated with reproductive benefits, suggesting that reproduction-related selection inadvertently contributes to the myopia epidemic. We estimate that, in the UK alone, natural selection adds more than 100 000 myopia cases per generation, and argue that antagonistic pleiotropy be broadly considered in explaining the spreads of apparently disadvantageous phenotypes in humans and beyond.

The Phenotypic Consequences of Genetic Divergence between Admixed Latin American Populations: Antioquia and Chocó, Colombia

Genome-wide association studies have uncovered thousands of genetic variants that are associated with a wide variety of human traits. Knowledge of how trait-associated variants are distributed within and between populations can provide insight into the genetic basis of group-specific phenotypic differences, particularly for health-related traits. We analyzed the genetic divergence levels for 1) individual trait-associated variants and 2) collections of variants that function together to encode polygenic traits, between two neighboring populations in Colombia that have distinct demographic profiles: Antioquia (Mestizo) and Chocó (Afro-Colombian). Genetic ancestry analysis showed 62% European, 32% Native American, and 6% African ancestry for Antioquia compared with 76% African, 10% European, and 14% Native American ancestry for Chocó, consistent with demography and previous results. Ancestry differences can confound cross-population comparison of polygenic risk scores (PRS); however, we did not find any systematic bias in PRS distributions for the two populations studied here, and population-specific differences in PRS were, for the most part, small and symmetrically distributed around zero. Both genetic differentiation at individual trait-associated single nucleotide polymorphisms and population-specific PRS differences between Antioquia and Chocó largely reflected anthropometric phenotypic differences that can be readily observed between the populations along with reported disease prevalence differences. Cases where population-specific differences in genetic risk did not align with observed trait (disease) prevalence point to the importance of environmental contributions to phenotypic variance, for both infectious and complex, common disease. The results reported here are distributed via a web-based platform for searching trait-associated variants and PRS divergence levels at http://map.chocogen.com (last accessed August 12, 2020).

Three Signatures of Adaptive Polymorphism Exemplified by Malaria-Associated Genes

Malaria has been one of the strongest selective pressures on our species. Many of the best-characterized cases of adaptive evolution in humans are in genes tied to malaria resistance. However, the complex evolutionary patterns at these genes are poorly captured by standard scans for nonneutral evolution. Here, we present three new statistical tests for selection based on population genetic patterns that are observed more than once among key malaria resistance loci. We assess these tests using forward-time evolutionary simulations and apply them to global whole-genome sequencing data from humans, and thus we show that they are effective at distinguishing selection from neutrality. Each test captures a distinct evolutionary pattern, here called Divergent Haplotypes, Repeated Shifts, and Arrested Sweeps, associated with a particular period of human prehistory. We clarify the selective signatures at known malaria-relevant genes and identify additional genes showing similar adaptive evolutionary patterns. Among our top outliers, we see a particular enrichment for genes involved in erythropoiesis and for genes previously associated with malaria resistance, consistent with a major role for malaria in shaping these patterns of genetic diversity. Polymorphisms at these genes are likely to impact resistance to malaria infection and contribute to ongoing host-parasite coevolutionary dynamics.

Effect of sodium acetate on serum activity of glucose-6-phosphate dehydrogenase in Plasmodium berghei-infected mice

Malaria is a global health problem with severe morbidity and mortality in Sub-Saharan Africa. Resistance of Plasmodium spp to the current anti-malaria drugs necessitates further search for novel effective drugs. This study, therefore, investigated the effect of sodium acetate on glucose-6-phosphate dehydrogenase in Plasmodium berghei-infected mice. Thirty male Albino mice were randomly distributed into 6 groups, A-F. Animals in Groups B-F were inoculated with P. berghei, intraperitoneally. Subsequently, Group C mice were treated with 20 mg/kg chloroquine, while groups D, E and F received 25, 50 and 100 mg/kg sodium acetate, respectively. All treatments were administered orally for 4 days. At the end of the experiment, animals were sacrificed by cervical dislocation and blood was collected via cardiac puncture for the analyses of serum glucose-6-phosphate dehydrogenase (G6PD), uric acid and lipid profile. Our results showed that Sodium acetate (50 and 100 mg/kg) significantly reduced (p < 0.05) parasitaemia (67.11% and 77.62%, respectively) than chloroquine (61.73%). Besides, body weight and serum G6PD activity in P. berghei infection were improved. Similarly, sodium acetate reduced elevated serum uric acid. Effects of sodium acetate and chloroquine on biochemical parameters were comparable (p > 0.05) but atherogenic lipid ratios were not affected by sodium acetate. These data put together suggested that activity of sodium acetate may be harnessed for development of novel anti-malaria drugs. However, more studies are required to delineate its mechanisms of action.

Epidemiology, evolutionary origin, and malaria-induced positive selection effects of G6PD-deficient alleles in Chinese populations

BACKGROUND

Although glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common inherited disorder in the Chinese population, there is scarce evidence regarding the epidemiology, evolutionary origin, and malaria-induced positive selection effects of G6PD-deficient alleles in various Chinese ethnic populations.

METHODS

We performed a large population-based screening (n = 15,690) to examine the impact of selection on human nucleotide diversity and to infer the evolutionary history of the most common deficiency alleles in Chinese populations.

RESULTS

The frequencies of G6PD deficiency ranged from 0% to 11.6% in 12 Chinese ethnic populations. A frequency map based on geographic information showed that G6PD deficiency was highly correlated with historical malaria prevalence in China and was affected by altitude and latitude. The five most frequently occurring G6PD gene variants were NM_001042351.3:c.1376G>T, NM_001042351.3:c.1388G>A, NM_001042351.3:c.95A>G, NM_001042351.3:c.1311T>C, and NM_001042351.3:c.1024C>T, which were distributed with ethnic features. A pathogenic but rarely reported variant site (NM_001042351.3:c.448G>A) was identified in this study. Bioinformatic analysis revealed a strong and recent positive selection targeting the NM_001042351.3:c.1376G>T allele that originated in the past 3125 to 3750 years and another selection targeting the NM_001042351.3:c.1388G>A allele that originated in the past 5000 to 6000 years. Additionally, both alleles originated from a single ancestor.

CONCLUSION

These results indicate that malaria has had a major impact on the Chinese genome since the introduction of rice agriculture.

Common variants in signaling transcription-factor-binding sites drive phenotypic variability in red blood cell traits

Genome-wide association studies identify genomic variants associated with human traits and diseases. Most trait-associated variants are located within cell-type-specific enhancers, but the molecular mechanisms governing phenotypic variation are less well understood. Here, we show that many enhancer variants associated with red blood cell (RBC) traits map to enhancers that are co-bound by lineage-specific master transcription factors (MTFs) and signaling transcription factors (STFs) responsive to extracellular signals. The majority of enhancer variants reside on STF and not MTF motifs, perturbing DNA binding by various STFs (BMP/TGF-β-directed SMADs or WNT-induced TCFs) and affecting target gene expression. Analyses of engineered human blood cells and expression quantitative trait loci verify that disrupted STF binding leads to altered gene expression. Our results propose that the majority of the RBC-trait-associated variants that reside on transcription-factor-binding sequences fall in STF target sequences, suggesting that the phenotypic variation of RBC traits could stem from altered responsiveness to extracellular stimuli.

Study of Glucose-6-Phosphate Dehydrogenase (G6PD) Deficiency and Genotype Polymorphism of G6PD B and G6PD (A+/A-) in Patients Treated for Plasmodium vivax Malaria in a Tertiary Care Hospital in North East India

Introduction

Glucose-6-phosphate dehydrogenase (G6PD) enzyme deficiency is the most common enzymopathy in humans, and its distribution has been historically described to be closely associated with that of malaria. North East India provides optimal conditions for transmission of malaria and bears a considerable burden of Plasmodium vivax (P. vivax) malaria. Primaquine, a mainstay in the treatment of vivax malaria, may trigger episodes of acute hemolysis in patients with G6PD deficiency. The present study sought to delineate the frequency and genotypes of G6PD deficiency among patients suffering from vivax malaria infections. 

Methods

Blood specimens from 80 individuals diagnosed with vivax malaria underwent enzyme assay for G6PD deficiency. Samples with deficient phenotype underwent isolation of DNA using a genomic DNA isolation kit (Qiagen India Pvt. Ltd., New Delhi, India). The genomic DNA underwent amplification, serial denaturation, annealing, extension, final extension followed by digestion with restriction endonucleases Nla III and Fok I. The digested products were subjected to horizontal agarose electrophoresis for the separation of digested fragments. Samples without nucleotide 376 adenine→guanine (A→G) mutation were classified as G6PD B. Those with the mutation were further classified into G6PD A(+) and G6PD A(-) based on the presence of Nla III site.

Results

Twenty-seven out of 80 individuals (33.75%) with P. vivax malaria were found to have G6PD deficiency, of which a majority (n=24) had G6PD B genotype. Three individuals had Asparagine→Aspartic Acid mutation at position 376 (A→G), of which G6PD A(+) and G6PD A(-) were present in two and one cases, respectively.

Conclusion

G6PD deficiency was noted in about a third of patients with vivax malaria. Since primaquine therapy is contraindicated in this group of patients, there is a rationale for looking into screening patients with vivax malaria from the region prior to primaquine therapy. Further large scale studies may substantiate this and help in better genotypic and geographic characterization of G6PD deficiency in the region. 

Adaptive immunity selects against malaria infection blocking mutations

The mutation responsible for Duffy negativity, which impedes Plasmodium vivax infection, has reached high frequencies in certain human populations. Conversely, mutations capable of blocking the more lethal P. falciparum have not succeeded in malarious zones. Here we present an evolutionary-epidemiological model of malaria which demonstrates that if adaptive immunity against the most virulent effects of malaria is gained rapidly by the host, mutations which prevent infection per se are unlikely to succeed. Our results (i) explain the rarity of strain-transcending P. falciparum infection blocking adaptations in humans; (ii) make the surprising prediction that mutations which block P. falciparum infection are most likely to be found in populations experiencing low or infrequent malaria transmission, and (iii) predict that immunity against some of the virulent effects of P. vivax malaria may be built up over the course of many infections.

Identifying adaptive alleles in the human genome: from selection mapping to functional validation

The suite of phenotypic diversity across geographically distributed human populations is the outcome of genetic drift, gene flow, and natural selection throughout human evolution. Human genetic variation underlying local biological adaptations to selective pressures is incompletely characterized. With the emergence of population genetics modeling of large-scale genomic data derived from diverse populations, scientists are able to map signatures of natural selection in the genome in a process known as selection mapping. Inferred selection signals further can be used to identify candidate functional alleles that underlie putative adaptive phenotypes. Phenotypic association, fine mapping, and functional experiments facilitate the identification of candidate adaptive alleles. Functional investigation of candidate adaptive variation using novel techniques in molecular biology is slowly beginning to unravel how selection signals translate to changes in biology that underlie the phenotypic spectrum of our species. In addition to informing evolutionary hypotheses of adaptation, the discovery and functional annotation of adaptive alleles also may be of clinical significance. While selection mapping efforts in non-European populations are growing, there remains a stark under-representation of diverse human populations in current public genomic databases, of both clinical and non-clinical cohorts. This lack of inclusion limits the study of human biological variation. Identifying and functionally validating candidate adaptive alleles in more global populations is necessary for understanding basic human biology and human disease.

Risk of pneumococcal bacteremia in Kenyan children with glucose-6-phosphate dehydrogenase deficiency

BACKGROUND

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common enzyme deficiency state in humans. The clinical phenotype is variable and includes asymptomatic individuals, episodic hemolysis induced by oxidative stress, and chronic hemolysis. G6PD deficiency is common in malaria-endemic regions, an observation hypothesized to be due to balancing selection at the G6PD locus driven by malaria. G6PD deficiency increases risk of severe malarial anemia, a key determinant of invasive bacterial disease in malaria-endemic settings. The pneumococcus is a leading cause of invasive bacterial infection and death in African children. The effect of G6PD deficiency on risk of pneumococcal disease is undefined. We hypothesized that G6PD deficiency increases pneumococcal disease risk and that this effect is dependent upon malaria.

METHODS

We performed a genetic case-control study of pneumococcal bacteremia in Kenyan children stratified across a period of falling malaria transmission between 1998 and 2010.

RESULTS

Four hundred twenty-nine Kenyan children with pneumococcal bacteremia and 2677 control children were included in the study. Among control children, G6PD deficiency, secondary to the rs1050828 G>A mutation, was common, with 11.2% (n = 301 of 2677) being hemi- or homozygotes and 33.3% (n = 442 of 1329) of girls being heterozygotes. We found that G6PD deficiency increased the risk of pneumococcal bacteremia, but only during a period of high malaria transmission (P = 0.014; OR 2.33, 95% CI 1.19-4.57). We estimate that the population attributable fraction of G6PD deficiency on risk of pneumococcal bacteremia in areas under high malaria transmission is 0.129.

CONCLUSIONS

Our data demonstrate that G6PD deficiency increases risk of pneumococcal bacteremia in a manner dependent on malaria. At the population level, the impact of G6PD deficiency on invasive pneumococcal disease risk in malaria-endemic regions is substantial. Our study highlights the infection-associated morbidity and mortality conferred by G6PD deficiency in malaria-endemic settings and adds to our understanding of the potential indirect health benefits of improved malaria control.

Adaptive reduction of male gamete number in the selfing plant Arabidopsis thaliana

The number of male gametes is critical for reproductive success and varies between and within species. The evolutionary reduction of the number of pollen grains encompassing the male gametes is widespread in selfing plants. Here, we employ genome-wide association study (GWAS) to identify underlying loci and to assess the molecular signatures of selection on pollen number-associated loci in the predominantly selfing plant Arabidopsis thaliana. Regions of strong association with pollen number are enriched for signatures of selection, indicating polygenic selection. We isolate the gene REDUCED POLLEN NUMBER1 (RDP1) at the locus with the strongest association. We validate its effect using a quantitative complementation test with CRISPR/Cas9-generated null mutants in nonstandard wild accessions. In contrast to pleiotropic null mutants, only pollen numbers are significantly affected by natural allelic variants. These data support theoretical predictions that reduced investment in male gametes is advantageous in predominantly selfing species.

Plasmodium vivax Malaria Viewed through the Lens of an Eradicated European Strain

The protozoan Plasmodium vivax is responsible for 42% of all cases of malaria outside Africa. The parasite is currently largely restricted to tropical and subtropical latitudes in Asia, Oceania, and the Americas. Though, it was historically present in most of Europe before being finally eradicated during the second half of the 20th century. The lack of genomic information on the extinct European lineage has prevented a clear understanding of historical population structuring and past migrations of P. vivax. We used medical microscope slides prepared in 1944 from malaria-affected patients from the Ebro Delta in Spain, one of the last footholds of malaria in Europe, to generate a genome of a European P. vivax strain. Population genetics and phylogenetic analyses placed this strain basal to a cluster including samples from the Americas. This genome allowed us to calibrate a genomic mutation rate for P. vivax, and to estimate the mean age of the last common ancestor between European and American strains to the 15th century. This date points to an introduction of the parasite during the European colonization of the Americas. In addition, we found that some known variants for resistance to antimalarial drugs, including Chloroquine and Sulfadoxine, were already present in this European strain, predating their use. Our results shed light on the evolution of an important human pathogen and illustrate the value of antique medical collections as a resource for retrieving genomic information on pathogens from the past.

Evidence of positively selected G6PD A- allele reduces risk of Plasmodium falciparum infection in African population on Bioko Island

BACKGROUND

Glucose-6-phosphate dehydrogenase (G6PD) is an essential enzyme that protects red blood cells from oxidative damage. Although G6PD-deficient alleles appear to confer a protective effect of malaria, the link with clinical protection against Plasmodium infection is conflicting.

METHODS

A case-control study was conducted on Bioko Island, Equatorial Guinea and further genotyping analysis used to detect natural selection of the G6PD A- allele.

RESULTS

Our results showed G6PD A- allele could significantly reduce the risk of Plasmodium falciparum infection in male individuals (adjusted odds ratio [AOR], 0.43; 95% confidence interval [CI], 0.20-0.93; p < .05) and homozygous female individuals (AOR, 0.11; 95% CI, 0.01-0.84; p < .05). Additionally, the parasite densities were significantly different in the individuals with different G6PD A- alleles and individual levels of G6PD enzyme activity. The pattern of linkage disequilibrium and results of the long-range haplotype test revealed a strong selective signature in the region encompassing the G6PD A- allele over the past 6,250 years. The network of inferred haplotypes suggested a single origin of the G6PD A- allele in Africans.

CONCLUSION

Our findings demonstrate that glucose-6-phosphate dehydrogenase (G6PD) A- allele could reduce the risk of P. falciparum infection in the African population and indicate that malaria has a recent positive selection on G6PD A- allele.

THE EXPOSOME IN HUMAN EVOLUTION: FROM DUST TO DIESEL

Global exposures to air pollution and cigarette smoke are novel in human evolutionary history and are associated with about 16 million premature deaths per year. We investigate the history of the human exposome for relationships between novel environmental toxins and genetic changes during human evolution in six phases. Phase I: With increased walking on savannas, early human ancestors inhaled crustal dust, fecal aerosols, and spores; carrion scavenging introduced new infectious pathogens. Phase II: Domestic fire exposed early Homo to novel toxins from smoke and cooking. Phases III and IV: Neolithic to preindustrial Homo sapiens incurred infectious pathogens from domestic animals and dense communities with limited sanitation. Phase V: Industrialization introduced novel toxins from fossil fuels, industrial chemicals, and tobacco at the same time infectious pathogens were diminishing. Thereby, pathogen-driven causes of mortality were replaced by chronic diseases driven by sterile inflammogens, exogenous and endogenous. Phase VI: Considers future health during global warming with increased air pollution and infections. We hypothesize that adaptation to some ancient toxins persists in genetic variations associated with inflammation and longevity.

Prevalence and distribution of G6PD deficiency: implication for the use of primaquine in malaria treatment in Ethiopia

BACKGROUND

G6PD enzyme deficiency is a common enzymatic X-linked disorder. Deficiency of the G6PD enzyme can cause free radical-mediated oxidative damage to red blood cells, leading to premature haemolysis. Treatment of Plasmodium vivax malaria with primaquine poses a potential risk of mild to severe acute haemolytic anaemia in G6PD deficient people. In this study, the prevalence and distribution of G6PD mutations were investigated across broad areas of Ethiopia, and tested the association between G6PD genotype and phenotype with the goal to provide additional information relevant to the use of primaquine in malaria treatment.

METHODS

This study examined G6PD mutations in exons 3-11 for 344 febrile patient samples collected from seven sites across Ethiopia. In addition, the G6PD enzyme level of 400 febrile patient samples from Southwestern Ethiopia was determined by the CareStart™ biosensor. The association between G6PD phenotype and genotype was examined by Fisher exact test on a subset of 184 samples.

RESULTS

Mutations were observed at three positions of the G6PD gene. The most common G6PD mutation across all sites was A376G, which was detected in 21 of 344 (6.1%) febrile patients. Thirteen of them were homozygous and eight were heterozygous for this mutation. The G267+119C/T mutation was found in 4 (1.2%) individuals in South Ethiopia, but absent in other sites. The G1116A mutation was also found in 4 (1.2%) individuals from East and South Ethiopia. For the 400 samples in the south, 17 (4.25%) were shown to be G6PD-deficient. G6PD enzyme level was not significantly different by age or gender. Among a subset of 202 febrile patients who were diagnosed with malaria, 11 (5.45%) were G6PD-deficient. These 11 infected samples were diagnosed with Plasmodium vivax by microscopy. Parasitaemia was not significantly different between the G6PD-deficient and G6PD-normal infections.

CONCLUSIONS

The prevalence of G6PD deficiency is modest among febrile patients in Ethiopia. G6PD deficiency testing is thus recommended before administrating primaquine for radical cure of P. vivax infected patients. The present study did not indicate a significant association between G6PD gene mutations and enzyme levels.

The Future of Genomic Studies Must Be Globally Representative: Perspectives from PAGE

The past decade has seen a technological revolution in human genetics that has empowered population-level investigations into genetic associations with phenotypes. Although these discoveries rely on genetic variation across individuals, association studies have overwhelmingly been performed in populations of European descent. In this review, we describe limitations faced by single-population studies and provide an overview of strategies to improve global representation in existing data sets and future human genomics research via diversity-focused, multiethnic studies. We highlight the successes of individual studies and meta-analysis consortia that have provided unique knowledge. Additionally, we outline the approach taken by the Population Architecture Using Genomics and Epidemiology (PAGE) study to develop best practices for performing genetic epidemiology in multiethnic contexts. Finally, we discuss how limiting investigations to single populations impairs findings in the clinical domain for both rare-variant identification and genetic risk prediction.

Genetic Variation for Ontogenetic Shifts in Metabolism Underlies Physiological Homeostasis in Drosophila

Organismal physiology emerges from metabolic pathways and structures that can vary across development and among individuals. Matoo, Julick, and Montooth found significant variation, both genetic and ontogenetic, in mitochondrial physiology in wild-type and mitochondrial-nuclear...

Organismal physiology emerges from metabolic pathways and subcellular structures like the mitochondria that can vary across development and among individuals. Here, we tested whether genetic variation at one level of physiology can be buffered at higher levels of biological organization during development by the inherent capacity for homeostasis in physiological systems. We found that the fundamental scaling relationship between mass and metabolic rate, as well as the oxidative capacity per mitochondria, changed significantly across development in the fruit fly Drosophila. However, mitochondrial respiration rate was maintained at similar levels across development. Furthermore, larvae clustered into two types—those that switched to aerobic, mitochondrial ATP production before the second instar, and those that relied on anaerobic, glycolytic production of ATP through the second instar. Despite genetic variation for the timing of this metabolic shift, metabolic rate in second-instar larvae was more robust to genetic variation than was the metabolic rate of other instars. We found that larvae with a mitochondrial-nuclear incompatibility that disrupts mitochondrial function had increased aerobic capacity and relied more on anaerobic ATP production throughout development relative to larvae from wild-type strains. By taking advantage of both ways of making ATP, larvae with this mitochondrial–nuclear incompatibility maintained mitochondrial respiratory capacity, but also had higher levels of whole-body reactive oxygen species and decreased mitochondrial membrane potential, potentially as a physiological defense mechanism. Thus, genetic defects in core physiology can be buffered at the organismal level via physiological plasticity, and natural populations may harbor genetic variation for distinct metabolic strategies in development that generate similar organismal outcomes.

Natural Selection and Origin of a Melanistic Allele in North American Gray Wolves

Pigmentation is often used to understand how natural selection affects genetic variation in wild populations since it can have a simple genetic basis, and can affect a variety of fitness-related traits (e.g., camouflage, thermoregulation, and sexual display). In gray wolves, the K locus, a β-defensin gene, causes black coat color via a dominantly inherited KB allele. The allele is derived from dog-wolf hybridization and is at high frequency in North American wolf populations. We designed a DNA capture array to probe the geographic origin, age, and number of introgression events of the KB allele in a panel of 331 wolves and 20 dogs. We found low diversity in KB, but not ancestral ky, wolf haplotypes consistent with a selective sweep of the black haplotype across North America. Further, North American wolf KB haplotypes are monophyletic, suggesting that a single adaptive introgression from dogs to wolves most likely occurred in the Northwest Territories or Yukon. We use a new analytical approach to date the origin of the KB allele in Yukon wolves to between 1,598 and 7,248 years ago, suggesting that introgression with early Native American dogs was the source. Using population genetic simulations, we show that the K locus is undergoing natural selection in four wolf populations. We find evidence for balancing selection, specifically in Yellowstone wolves, which could be a result of selection for enhanced immunity in response to distemper. With these data, we demonstrate how the spread of an adaptive variant may have occurred across a species' geographic range.

Haplotype Loci Under Selection in Canadian Durum Wheat Germplasm Over 60 Years of Breeding: Association With Grain Yield, Quality Traits, Protein Loss, and Plant Height

Durum wheat was introduced in the southern prairies of western Canada in the late nineteenth century. Breeding efforts have mainly focused on improving quality traits to meet the pasta industry demands. For this study, 192 durum wheat lines were genotyped using the Illumina 90K Infinium iSelect assay, and resulted in a total of 14,324 polymorphic SNPs. Genetic diversity changed over time, declining during the first 20 years of breeding in Canada, then increased in the late 1980s and early 1990s. We scanned the genome for signatures of selection, using the total variance Fst-based outlier detection method (Lositan), the hierarchical island model (Arlequin) and the Bayesian genome scan method (BayeScan). A total of 407 outliers were identified and clustered into 84 LD-based haplotype loci, spanning all 14 chromosomes of the durum wheat genome. The association analysis detected 54 haplotype loci, of which 39% contained markers with a complete reversal of allelic state. This tendency to fixation of favorable alleles corroborates the success of the Canadian durum wheat breeding programs over time. Twenty-one haplotype loci were associated with multiple traits. In particular, hap_4B_1 explained 20.6, 17.9 and 16.6% of the phenotypic variance of pigment loss, pasta b∗ and dough extensibility, respectively. The locus hap_2B_9 explained 15.9 and 17.8% of the variation of protein content and protein loss, respectively. All these pleiotropic haplotype loci offer breeders the unique opportunity for further improving multiple traits, facilitating marker-assisted selection in durum wheat, and could help in identifying genes as functional annotations of the wheat genome become available.

Glucose-6-phosphate dehydrogenase deficiency allelic variants and their prevalence in malaria patients in Eritrea

Introduction

Glucose 6-phosphate dehydrogenase (G6PD) deficiency is the most common enzymopathy with a relatively high frequency in malaria-endemic regions. In Eritrea, there is scanty knowledge of G6PD deficiency. The aim of the study was to characterize and determine the prevalence of four common G6PD allelic variants.

Methods

Three hundred and fourteen dried blood spot samples from unrelated microscopically diagnosed malaria patient Eritrean ethnic groups living in five zobas (regions) of Eritrea were analysed by PCR-RFLP method to identify the G6PD B, G6PD A (A376G), G6PD A-(G202A), and G6PD Mediterranean (C563T) variants. To confirm the RFLP results, samples positive for A376G but negative for G202A variants were subjected to Sanger sequencing and a subset of PCR products (exon 5) directly sequenced to identify A376G and other mutations.

Results

For G6PD genotyping, G6PD B was detected in 87.5% and A376G detected in 12.5% of malaria patients, whereas G202A and C563T were absent. Bivariate Statistical analysis showed a statistically significant association between G6PD genotypes and zoba (P < 0.004 < 0.05). Sequencing revealed the expected A376G variant. In exon 5, four common (A376G) mutations, three uncommon mutations rs782669677 (535G→A) and one potentially new mutation (451G→C), relative to the reference, mRNA NM_001042351 were detected. Bioinformatic analysis of these mutations' potential functional impact suggests minimal effect on protein function.

Conclusion

This is the first report indicating that G6PD B and G6PD A genotypes are prevalent in Eritrea. Similar findings were reported in neighboring countries. Further studies including phenotype analysis are needed to corroborate the observed results.

Whole-Genome-Sequence-Based Haplotypes Reveal Single Origin of the Sickle Allele during the Holocene Wet Phase

Five classical designations of sickle haplotypes are made on the basis of the presence or absence of restriction sites and are named after the ethno-linguistic groups or geographic regions from which the individuals with sickle cell anemia originated. Each haplotype is thought to represent an independent occurrence of the sickle mutation rs334 (c.20A>T [p.Glu7Val] in HBB). We investigated the origins of the sickle mutation by using whole-genome-sequence data. We identified 156 carriers from the 1000 Genomes Project, the African Genome Variation Project, and Qatar. We classified haplotypes by using 27 polymorphisms in linkage disequilibrium with rs334. Network analysis revealed a common haplotype that differed from the ancestral haplotype only by the derived sickle mutation at rs334 and correlated collectively with the Central African Republic (CAR), Cameroon, and Arabian/Indian haplotypes. Other haplotypes were derived from this haplotype and fell into two clusters, one composed of Senegal haplotypes and the other composed of Benin and Senegal haplotypes. The near-exclusive presence of the original sickle haplotype in the CAR, Kenya, Uganda, and South Africa is consistent with this haplotype predating the Bantu expansions. Modeling of balancing selection indicated that the heterozygote advantage was 15.2%, an equilibrium frequency of 12.0% was reached after 87 generations, and the selective environment predated the mutation. The posterior distribution of the ancestral recombination graph yielded a sickle mutation age of 259 generations, corresponding to 7,300 years ago during the Holocene Wet Phase. These results clarify the origin of the sickle allele and improve and simplify the classification of sickle haplotypes.

Identifying the favored mutation in a positive selective sweep

Most approaches that capture signatures of selective sweeps in population genomics data do not identify the specific mutation favored by selection. We present iSAFE (for "integrated selection of allele favored by evolution"), a method that enables researchers to accurately pinpoint the favored mutation in a large region (∼5 Mbp) by using a statistic derived solely from population genetics signals. iSAFE does not require knowledge of demography, the phenotype under selection, or functional annotations of mutations.

Prevalence of glucose-6-phosphate dehydrogenase (G6PD) deficiency among malaria patients in Upper Myanmar

Background

Glucose-6-phosphate dehydrogenase (G6PD; EC 1.1.1.49) deficiency is one of the most common X-linked recessive hereditary disorders in the world. Primaquine (PQ) has been used for radical cure of P. vivax to prevent relapse. Recently, it is also used to reduce P. falciparum gametocyte carriage to block transmission. However, PQ metabolites oxidize hemoglobin and generate excessive reactive oxygen species which can trigger acute hemolytic anemia in malaria patients with inherited G6PD deficiency.

Methods

A total of 252 blood samples collected from malaria patients in Myanmar were used in this study. G6PD variant was analysed by a multiplex allele specific PCR kit, DiaPlexC™ G6PD Genotyping Kit [Asian type]. The accuracy of the multiplex allele specific PCR was confirmed by sequencing analysis.

Results

Prevalence and distribution of G6PD variants in 252 malaria patients in Myanmar were analysed. Six different types of G6PD allelic variants were identified in 50 (7 females and 43 males) malaria patients. The predominant variant was Mahidol (68%, 34/50), of which 91.2% (31/34) and 8.8% (3/34) were males and females, respectively. Other G6PD variants including Kaiping (18%, 9/50), Viangchan (6%, 3/50), Mediterranean (4%, 2/50), Union (2%, 1/50) and Canton (2%, 1/50) were also observed.

Conclusions

Results of this study suggest that more concern for proper and safe use of PQ as a radical cure of malaria in Myanmar is needed by combining G6PD deficiency test before PQ prescription. Establishment of a follow-up system to monitor potential PQ toxicity in malaria patients who are given PQ is also required.

Genetic Polymorphisms of Glucose-6-Phosphate Dehydrogenase in Lagos, Nigeria

Glucose-6-phosphate dehydrogenase (G6PD) is an essential enzyme in the pentose phosphate pathway that prevents oxidative damage to cells. This study determined the genotypic and allelic frequencies of G6PD G202A and A376G and also investigated correlation between G6PD polymorphisms and hemoglobin (Hb) phenotypes in children in Lagos, Nigeria. Seventy-eight children [55 with Hb AA (β^Α/β^A) and 23 with Hb AS (β^Α/β^S) trait] and 65 Hb SS (β^S/β^S) (HBB: c.20A>T) subjects in steady state with age range between 5-15 years were recruited for the study. Hemoglobin phenotypes of all study participants were carried out using alkaline electrophoresis and solubility tests. Genomic DNA was extracted from whole blood and restriction fragment length polymorphism-polymerase chain reaction (RFLP-PCR) was used to determine the G202A and the A376G mutations of the G6PD gene. The genotype and allele distributions of G6PD G202A and A376G according to the Hb phenotypes were not statistically significant (p > 0.05). The minor allele frequency 202A was 0.15 (15.0%) and 0.14 (14.0%) in cases and controls, respectively. The overall frequency of 376G allele in the case group was 0.35 (35.0%) and 0.38 (38.0%) in the control group. No statistical significance was observed in the genotype and allele distributions of A376G in both the case and control groups (p > 0.05). The G6PD A- frequency in Hb SS subjects and the control group were 6.2 and 2.6%, respectively. G6PD G202A and A376G polymorphisms were not associated with Hb phenotypes and the allele distributions of 202A and 376G in this study are typical of West African populations.

Complex signatures of natural selection at GYPA

The human MN blood group antigens are isoforms of glycophorin A (GPA) encoded by the gene, GYPA, and are the most abundant erythrocyte sialoglycoproteins. The distribution of MN antigens has been widely studied in human populations yet the evolutionary and/or demographic factors affecting population variation remain elusive. While the primary function of GPA is yet to be discovered, it serves as the major binding site for the 175-kD erythrocyte-binding antigen (EB-175) of the malarial parasite, Plasmodium falciparum, a major selective pressure in recent human history. More specifically, exon two of GYPA encodes the receptor-binding ligand to which P. falciparum binds. Accordingly, there has been keen interest in understanding what impact, if any, natural selection has had on the distribution of variation in GYPA and exon two in particular. To this end, we resequenced GYPA in individuals sampled from both P. falciparum endemic (sub-Saharan Africa and South India) and non-endemic (Europe and East Asia) regions of the world. Observed patterns of variation suggest that GYPA has been subject to balancing selection in populations living in malaria endemic areas and in Europeans, but no such evidence was found in samples from East Asia, Oceania, and the Americas. These results are consistent with malaria acting as a selective pressure on GYPA, but also suggest that another selective force has resulted in a similar pattern of variation in Europeans. Accordingly, GYPA has perhaps a more complex evolutionary history, wherein on a global scale, spatially varying selective pressures have governed its natural history.

Single low-dose primaquine for blocking transmission of Plasmodium falciparum malaria - a proposed model-derived age-based regimen for sub-Saharan Africa

BACKGROUND

In 2012, the World Health Organization recommended blocking the transmission of Plasmodium falciparum with single low-dose primaquine (SLDPQ, target dose 0.25 mg base/kg body weight), without testing for glucose-6-phosphate dehydrogenase deficiency (G6PDd), when treating patients with uncomplicated falciparum malaria. We sought to develop an age-based SLDPQ regimen that would be suitable for sub-Saharan Africa.

METHODS

Using data on the anti-infectivity efficacy and tolerability of primaquine (PQ), the epidemiology of anaemia, and the risks of PQ-induced acute haemolytic anaemia (AHA) and clinically significant anaemia (CSA), we prospectively defined therapeutic-dose ranges of 0.15-0.4 mg PQ base/kg for children aged 1-5 years and 0.15-0.5 mg PQ base/kg for individuals aged ≥6 years (therapeutic indices 2.7 and 3.3, respectively). We chose 1.25 mg PQ base for infants aged 6-11 months because they have the highest rate of baseline anaemia and the highest risks of AHA and CSA. We modelled an anthropometric database of 661,979 African individuals aged ≥6 months (549,127 healthy individuals, 28,466 malaria patients and 84,386 individuals with other infections/illnesses) by the Box-Cox transformation power exponential and tested PQ doses of 1-15 mg base, selecting dosing groups based on calculated mg/kg PQ doses.

RESULTS

From the Box-Cox transformation power exponential model, five age categories were selected: (i) 6-11 months (n = 39,886, 6.03%), (ii) 1-5 years (n = 261,036, 45.46%), (iii) 6-9 years (n = 20,770, 3.14%), (iv) 10-14 years (n = 12,155, 1.84%) and (v) ≥15 years (n = 328,132, 49.57%) to receive 1.25, 2.5, 5, 7.5 and 15 mg PQ base for corresponding median (1st and 99th centiles) mg/kg PQ base of: (i) 0.16 (0.12-0.25), (ii) 0.21 (0.13-0.37), (iii) 0.25 (0.16-0.38), (iv) 0.26 (0.15-0.38) and (v) 0.27 (0.17-0.40). The proportions of individuals predicted to receive optimal therapeutic PQ doses were: 73.2 (29,180/39,886), 93.7 (244,537/261,036), 99.6 (20,690/20,770), 99.4 (12,086/12,155) and 99.8% (327,620/328,132), respectively.

CONCLUSIONS

We plan to test the safety of this age-based dosing regimen in a large randomised placebo-controlled trial (ISRCTN11594437) of uncomplicated falciparum malaria in G6PDd African children aged 0.5 - 11 years. If the regimen is safe and demonstrates adequate pharmacokinetics, it should be used to support malaria elimination.

Glucose-6-Phosphate Dehydrogenase Deficiency Genetic Variants in Malaria Patients in Southwestern Ethiopia

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is an X-linked erythrocyte enzyme disorder with relevance to malaria treatment policy. Treatment with the antimalarial primaquine can result in hemolytic anemia in G6PD-deficient patients. With increased interest in primaquine use, it is important to identify G6PD variants in Ethiopia to inform malaria treatment policy. In the present study, mutations in the G6PD gene are identified in a sample of patients with malaria in Jimma town in southwest Ethiopia. Plasmodium species of infection were confirmed using polymerase chain reaction (PCR) and gel electrophoresis. PCR and Sanger sequencing were performed to observe a portion of the G6PD gene where the common G6PD mutations (A376G, G202A, and C563T) are found. Molecular analysis revealed that most of the samples were single Plasmodium vivax infections (83.7%). For G6PD genotyping, A376G was detected in 23.26% of individuals, whereas G202A and C563T were absent. Three other uncommon mutations were identified: rs782669677 (535G→A), rs370658483, (485 + 37 G→T), and a new mutation at chrX:154535443(C→T). Bioinformatic analysis of these mutations' potential functional impact suggests minimal effect on protein function. The discovery of both common and uncommon G6PD mutations contributes to the discussion on G6PD deficiency and appropriate primaquine treatment in Ethiopia.

Prevalence, genetic variants and clinical implications of G-6-PD deficiency in Burkina Faso: a systematic review

BACKGROUND

It is now well-known that some antimalarials such as primaquine may induce severe hemolytic anemia in people with G-6-PD deficiency. Antimalarial drug prescriptions must, therefore take into account the patient's G-6-PD status in malaria endemic areas such as Burkina Faso, where the prevalence of this genetic abnormality is relatively high. Although great clinical heterogeneity is observed depending on the molecular nature of the deficiency and the residual enzyme activity in the red blood cell, there is very poor data on the prevalence of G-6-PD deficiency and the distribution of involved genetic variants in Burkina Faso. In this systematic review, we present a synthesis of the various studies carried out on the G-6-PD deficiency in Burkina Faso in order to determine its prevalence, probable distribution of the genetic variants involved and their clinical implications for a national systematic screening policy among the groups most vulnerable to malaria.

METHODS

A systematic review was carried out to analyze available published data on the prevalence, phenotypes and mutations responsible for G-6-PD deficiency in Burkina Faso. The key words used were "G-6-PD deficiency AND Burkina Faso" or "Déficit en G-6-PD AND Burkina Faso" in French. To identify the relevant articles, two independent reviewers reviewed the titles, abstracts and the full text of the selected papers.

RESULTS

An average prevalence of 16.6% (183/1100; CI 95%: 0.145-0.190) and 6.5% (69/1066; CI 95%: 0.051-0.081) of G-6-PD deficiency was found respectively in men and women in this systematic review. Although the predominance (99.8% of G-6-PD deficient cases) of 202A/376G G-6-PD A- variant, the Santamaria and Betica Selma variants were identified in Burkina Faso. Independently of the method used, the enzymatic deficiency was significantly higher in males (2.5-20.5%) compared to females (3.3-12.3%).

CONCLUSION

This systematic review suggests that despite the ubiquity of the 202A/376G G-6-PD A- variant in Burkina Faso, it will be necessary to consider the Santamaria and Betica Selma variants although their frequencies remain to be specified. A systematic screening of the G-6-PD deficiency is also needed to prevent the occurrence of iatrogenic hemolytic accidents.

A trade off between catalytic activity and protein stability determines the clinical manifestations of glucose-6-phosphate dehydrogenase (G6PD) deficiency

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common polymorphism and enzymopathy in humans, affecting approximately 400 million people worldwide. It is responsible for various clinical manifestations, including favism, hemolytic anemia, chronic non-spherocytic hemolytic anemia, spontaneous abortion, and neonatal hyperbilirubinemia. Understanding the molecular mechanisms underlying the severity of G6PD deficiency is of great importance but that of many G6PD variants are still unknown. In this study, we report the construction, expression, purification, and biochemical characterization in terms of kinetic properties and stability of five clinical G6PD variants-G6PD Bangkok, G6PD Bangkok noi, G6PD Songklanagarind, G6PD Canton+Bangkok noi, and G6PD Union+Viangchan. G6PD Bangkok and G6PD Canton+Bangkok noi showed a complete loss of catalytic activity and moderate reduction in thermal stability when compared with the native G6PD. G6PD Bangkok noi and G6PD Union+Viangchan showed a significant reduction in catalytic efficiency, whereas G6PD Songklanagarind showed a catalytic activity comparable to the wild-type enzyme. The Union+Viangchan mutation showed a remarkable effect on the global stability of the enzyme. In addition, our results indicate that the location of mutations in G6PD variants affects their catalytic activity, stability, and structure. Hence, our results provide a molecular explanation for clinical manifestations observed in individuals with G6PD deficiency.