Effects of natural selection and gene conversion on the evolution of human glycophorins coding for MNS blood polymorphisms in malaria-endemic African populations

Characterisation of the erythrocyte invasion phenotype of FCB-2: A South American P. falciparum reference strain

The extent of parasite adaptive capability involved in erythrocyte invasion represents a significant challenge for the development of a Plasmodium falciparum vaccine. The parasite's geographical and populational origin may influence such adaptive behaviour; in vitro culture-adapted parasite strains are typically used for such studies. Previous studies have reported invasion phenotypes in strains from Africa and Asia and, to a lesser extent, from Latin America. This study was aimed at expanding the pool of characterised parasite strains from Latin America by describing the invasion phenotype of the P. falciparum Colombia Bogotá 2 (FCB2) strain. The FCB2 genome was sequenced and erythrocyte invasion ligand sequences were analysed and compared to other previously reported ones. RT-PCR was used for assessing Pfeba family erythrocyte invasion ligands and reticulocyte binding homologue (Pfrh) gene transcription. A flow cytometry-based erythrocyte invasion assay (using enzymatically-treated erythrocytes) was used for determining the FCB2 strain's invasion phenotype. The P. falciparum FCB2 genome sequence was analysed, bearing in mind that prolonged in vitro parasite culture may affect its genome sequence and, in some cases, lead to the deletion of certain genes; it was demonstrated that all erythrocyte invasion ligand gene sequences studied here were preserved. Comparative analysis showed that the target genome sequences were conserved whereas transcriptional analysis highlighted Pfebas and Pfrhs gene expression. Erythrocyte invasion analysis demonstrated that the FCB2 strain has a sialic acid-resistant invasion phenotype.

Adaptation in human immune cells residing in tissues at the frontline of infections

Human immune cells are under constant evolutionary pressure, primarily through their role as first line of defence against pathogens. Most studies on immune adaptation are, however, based on protein-coding genes without considering their cellular context. Here, using data from the Human Cell Atlas, we infer the gene adaptation rate of the human immune landscape at cellular resolution. We find abundant cell types, like progenitor cells during development and adult cells in barrier tissues, to harbour significantly increased adaptation rates. We confirm the adaptation of tissue-resident T and NK cells in the adult lung located in compartments directly facing external challenges, such as respiratory pathogens. Analysing human iPSC-derived macrophages responding to various challenges, we find adaptation in early immune responses. Together, our study suggests host benefits to control pathogen spread at early stages of infection, providing a retrospect of forces that shaped the complexity, architecture, and function of the human body.

Genetic variation of glycophorins and infectious disease

Glycophorins are transmembrane proteins of red blood cells (RBCs), heavily glycosylated on their external-facing surface. In humans, there are four glycophorin proteins, glycophorins A, B, C and D. Glycophorins A and B are encoded by two similar genes GYPA and GYPB, and glycophorin C and glycophorin D are encoded by a single gene, GYPC. The exact function of glycophorins remains unclear. However, given their abundance on the surface of RBCs, it is likely that they serve as a substrate for glycosylation, giving the RBC a negatively charged, complex glycan “coat”. GYPB and GYPE (a closely related pseudogene) were generated from GYPA by two duplication events involving a 120-kb genomic segment between 10 and 15 million years ago. Non-allelic homologous recombination between these 120-kb repeats generates a variety of duplication alleles and deletion alleles, which have been systematically catalogued from genomic sequence data. One allele, called DUP4, encodes the Dantu NE blood type and is strongly protective against malaria as it alters the surface tension of the RBC membrane. Glycophorins interact with other infectious pathogens, including viruses, as well as the malarial parasite Plasmodium falciparum, but the role of glycophorin variation in mediating the effects of these pathogens remains underexplored.

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.

Addressing Africa's pandemic puzzle: Perspectives on COVID-19 transmission and mortality in sub-Saharan Africa

To date, SARS-CoV-2 (the virus that causes COVID-19) has spread to almost every region of the world, infecting millions and resulting in the deaths of hundreds of thousands of people. Although it was predicted that Africa would suffer a massive loss of life due to this pandemic, the number of COVID-19 cases has been relatively low across the continent. Researchers have speculated that several factors may be responsible for this outcome in Africa, including the extensive experience that countries have with infectious diseases and the young median age of their populations. However, it is still important for African countries to adopt aggressive and bold approaches against COVID-19, in case the nature of the pandemic changes. This short review will summarize the status of the outbreak in Africa and propose possible reasons for current trends, as well as discuss interventions aimed at preventing a rapid increase in the number of COVID-19 cases in the future.

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.

Two Prevalent ∼100-kb GYPB Deletions Causative of the GPB-Deficient Blood Group MNS Phenotype S-s-U- in Black Africans

The U antigen (MNS5) is one of 49 antigens belonging to the MNS blood group system (ISBT002) carried on glycophorins A (GPA) and B (GPB). U is present on the red blood cells in almost all Europeans and Asians but absent in approximately 1.0% of Black Africans. U negativity coincides with negativity for S (MNS3) and s (MNS4) on GPB, thus be called S-s-U-, and is thought to arise from homozygous deletion of GYPB. Little is known about the molecular background of these deletions. Bioinformatic analysis of the 1000 Genomes Project data revealed several candidate regions with apparent deletions in GYPB. Highly specific Gap-PCRs, only resulting in positive amplification from DNAs with deletions present, allowed for the exact genetic localization of 3 different breakpoints; 110.24- and 103.26-kb deletions were proven to be the most frequent in Black Americans and Africans. Among 157 CEPH DNAs, deletions in 6 out of 8 African ethnicities were present. Allele frequencies of the deletions within African ethnicities varied greatly and reached a cumulative 23.3% among the Mbuti Pygmy people from the Congo. Similar observations were made for U+var alleles, known to cause strongly reduced GPB expression. The 110- and 103-kb deletional GYPB haplotypes were found to represent the most prevalent hereditary factors causative of the MNS blood group phenotype S-s-U-. Respective GYPB deletions are now accessible by molecular detection of homo- and hemizygous transmission.

Genome-wide association studies of severe P. falciparum malaria susceptibility: progress, pitfalls and prospects

BACKGROUND

P. falciparum malaria has been recognized as one of the prominent evolutionary selective forces of human genome that led to the emergence of multiple host protective alleles. A comprehensive understanding of the genetic bases of severe malaria susceptibility and resistance can potentially pave ways to the development of new therapeutics and vaccines. Genome-wide association studies (GWASs) have recently been implemented in malaria endemic areas and identified a number of novel association genetic variants. However, there are several open questions around heritability, epistatic interactions, genetic correlations and associated molecular pathways among others. Here, we assess the progress and pitfalls of severe malaria susceptibility GWASs and discuss the biology of the novel variants.

RESULTS

We obtained all severe malaria susceptibility GWASs published thus far and accessed GWAS dataset of Gambian populations from European Phenome Genome Archive (EGA) through the MalariaGen consortium standard data access protocols. We noticed that, while some of the well-known variants including HbS and ABO blood group were replicated across endemic populations, only few novel variants were convincingly identified and their biological functions remain to be understood. We estimated SNP-heritability of severe malaria at 20.1% in Gambian populations and showed how advanced statistical genetic analytic methods can potentially be implemented in malaria susceptibility studies to provide useful functional insights.

CONCLUSIONS

The ultimate goal of malaria susceptibility study is to discover a novel causal biological pathway that provide protections against severe malaria; a fundamental step towards translational medicine such as development of vaccine and new therapeutics. Beyond singe locus analysis, the future direction of malaria susceptibility requires a paradigm shift from single -omics to multi-stage and multi-dimensional integrative functional studies that combines multiple data types from the human host, the parasite, the mosquitoes and the environment. The current biotechnological and statistical advances may eventually lead to the feasibility of systems biology studies and revolutionize malaria research.

SMIM1 at a glance; discovery, genetic basis, recent progress and perspectives

Recent elucidation of the genetic basis of the Vel blood group system has offered the field of blood transfusion medicine an additional consideration in determining the causes of hemolytic reactions after a patient is transfused. The identification of the SMIM1 gene to be responsible for the Vel blood group allows molecular based tools to be developed to further dissect the function of this antigen. Genetic signatures such as the homozygous 17 bp deletion and the heterozygous 17 bp deletion in combination with other single nucleotide polymorphisms (SNPs) and insertion sequences regulate the expression level of the gene. With this knowledge, it is now possible to study this antigen in-depth.

Nutrition and its role in human evolution

Our understanding of human evolution has improved rapidly over recent decades, facilitated by large-scale cataloguing of genomic variability amongst both modern and archaic humans. It seems clear that the evolution of the ancestors of chimpanzees and hominins separated 7-9 million years ago with some migration out of Africa by the earlier hominins; Homo sapiens slowly emerged as climate change resulted in drier, less forested African conditions. The African populations expanded and evolved in many different conditions with slow mutation and selection rates in the human genome, but with much more rapid mutation occurring in mitochondrial DNA. We now have evidence stretching back 300 000 years of humans in their current form, but there are clearly four very different large African language groups that correlate with population DNA differences. Then, about 50 000-100 000 years ago a small subset of modern humans also migrated out of Africa resulting in a persistent signature of more limited genetic diversity amongst non-African populations. Hybridization with archaic hominins occurred around this time such that all non-African modern humans possess some Neanderthal ancestry and Melanesian populations additionally possess some Denisovan ancestry. Human populations both within and outside Africa also adapted to diverse aspects of their local environment including altitude, climate, UV exposure, diet and pathogens, in some cases leaving clear signatures of patterns of genetic variation. Notable examples include haemoglobin changes conferring resistance to malaria, other immune changes and the skin adaptations favouring the synthesis of vitamin D. As humans migrated across Eurasia, further major mitochondrial changes occurred with some interbreeding with ancient hominins and the development of alcohol intolerance. More recently, an ability to retain lactase persistence into adulthood has evolved rapidly under the environmental stimulus of pastoralism with the ability to husband lactating ruminants. Increased amylase copy numbers seem to relate to the availability of starchy foods, whereas the capacity to desaturase and elongate monounsaturated fatty acids in different societies seems to be influenced by whether there is a lack of supply of readily available dietary sources of long-chain polyunsaturated fatty acids. The process of human evolution includes genetic drift and adaptation to local environments, in part through changes in mitochondrial and nuclear DNA. These genetic changes may underlie susceptibilities to some modern human pathologies including folate-responsive neural tube defects, diabetes, other age-related pathologies and mental health disorders.

Few Plasmodium falciparum merozoite ligand and erythrocyte receptor pairs show evidence of balancing selection

Erythrocyte surface proteins have been identified as receptors of Plasmodium falciparum merozoite proteins. The ligand-receptor interactions enable the parasite to invade human erythrocytes, initiating the clinical symptoms of malaria. These interactions are likely to have had an evolutionary impact on the genes that encode the ligand and receptor proteins. We used sequence data from Kilifi, Kenya to detect departures from neutrality in a paired analysis of P. falciparum merozoite ligands and their erythrocyte receptor genes from the same population. We genotyped parasite and human DNA obtained from 93 individuals with severe malaria. We examined six merozoite ligands EBA175, EBL1, EBA140, MSP1, Rh4 and Rh5, and their corresponding erythrocyte receptors, glycophorin (Gyp) A, GypB, GypC, band 3, complement receptor (CR) 1 and basigin, focusing on the regions involved in the ligand-receptor interactions. Positive Tajima's D values (>1) were observed only in the MSP1 C-terminal region and EBA175 region II, while negative values (<-1) were observed in EBL-1 region II, Rh4, basigin exons 3 and 5, CR1 exon 5, Gyp B exons 2, 3 and 4 and Gyp C exon 2. Additionally, ebl-1 region II and basigin exon 3 showed extreme negative values in all three tests, Tajima's D, Fu & Li D* and F*, ≤ - 2. A large majority of the erythrocyte receptor and merozoite genes have a negative Tajima's D even when compared with previously published whole genome data. Thus, highlighting EBA175 region II and MSP1-33, as outlier genes with a positive Tajima's D (>1). Both these genes contain multiple polymorphisms, which in the case of EBA175 may counteract receptor polymorphisms and/or evade host immune responses and in MSP1 the polymorphisms may primarily evade host immune responses.

Gene tree species tree reconciliation with gene conversion

Gene tree/species tree reconciliation is a recent decisive progress in phylogenetic methods, accounting for the possible differences between gene histories and species histories. Reconciliation consists in explaining these differences by gene-scale events such as duplication, loss, transfer, which translates mathematically into a mapping between gene tree nodes and species tree nodes or branches. Gene conversion is a frequent and important evolutionary event, which results in the replacement of a gene by a copy of another from the same species and in the same gene tree. Including this event in reconciliation models has never been attempted because it introduces a dependency between lineages, and standard algorithms based on dynamic programming become ineffective. We propose here a novel mathematical framework including gene conversion as an evolutionary event in gene tree/species tree reconciliation. We describe a randomized algorithm that finds, in polynomial running time, a reconciliation minimizing the number of duplications, losses and conversions in the case when their weights are equal. We show that the space of optimal reconciliations includes an analog of the last common ancestor reconciliation, but is not limited to it. Our algorithm outputs any optimal reconciliation with a non-null probability. We argue that this study opens a research avenue on including gene conversion in reconciliation, and discuss its possible importance in biology.

Parasite-Related Genetic and Epigenetic Aspects and Host Factors Influencing Plasmodium falciparum Invasion of Erythrocytes

Malaria, a disease caused by Plasmodium parasites, is widespread throughout tropical and sub-tropical regions worldwide; it mostly affects children and pregnant woman. Eradication has stalled despite effective prevention measures and medication being available for this disease; this has mainly been due to the parasite's resistance to medical treatment and the mosquito vector's resistance to insecticides. Tackling such resistance involves using renewed approaches and techniques for accruing a deep understanding of the parasite's biology, and developing new drugs and vaccines. Studying the parasite's invasion of erythrocytes should shed light on its ability to switch between invasion phenotypes related to the expression of gene sets encoding proteins acting as ligands during target cell invasion, thereby conferring mechanisms for evading a particular host's immune response and adapting to changes in target cell surface receptors. This review considers some factors influencing the expression of such phenotypes, such as Plasmodium's genetic, transcriptional and epigenetic characteristics, and explores some host-related aspects which could affect parasite phenotypes, aiming at integrating knowledge regarding this topic and the possible relationship between the parasite's biology and host factors playing a role in erythrocyte invasion.

The Malaria-Protective Human Glycophorin Structural Variant DUP4 Shows Somatic Mosaicism and Association with Hemoglobin Levels

Glycophorin A and glycophorin B are red blood cell surface proteins and are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A-glycophorin B fusion gene and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40%. Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant and reveal somatic variation in copy number of the glycophorin B-glycophorin A fusion gene. By developing a simple, specific, PCR-based assay for DUP4, we show that the DUP4 variant reaches a frequency of 13% in the population of a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation and show that DUP4 is associated with a malarial-related phenotype in a longitudinally followed population.

Genetic and evolutionary determinants of human population variation in immune responses

Humans display remarkable immune response variation when exposed to identical immune challenges. However, our understanding of the genetic, evolutionary, and environmental factors that impact this inter-individual and inter-population immune response heterogeneity is still in its early days. In this review, we discuss three fundamental questions concerning the recent evolution of the human immune system: the degree to which individuals from different populations vary in their innate immune responses, the genetic variants accounting for such differences, and the evolutionary mechanisms that led to the establishment of these variants in modern human populations. We also discuss how past selective events might have contributed to the uneven distribution of immune-related disorders across populations.

Patterns of shared signatures of recent positive selection across human populations

Signatures of recent positive selection often overlap across human populations, but the question of how often these overlaps represent a single ancestral event remains unresolved. If a single selective event spread across many populations, the same sweeping haplotype should appear in each population and the selective pressure could be common across populations and environments. Identifying such shared selective events could identify genomic loci and human traits important in recent history across the globe. In addition, genomic annotations that recently became available could help attach these signatures to a potential gene and molecular phenotype selected across populations. Here, we present a catalogue of selective sweeps in humans, and identify those that overlap and share a sweeping haplotype. We connect these sweep overlaps with potential biological mechanisms at several loci, including potential new sites of adaptive introgression, the glycophorin locus associated with malarial resistance and the alcohol dehydrogenase cluster associated with alcohol dependency.

Population genetic evidence for positive and purifying selection acting at the human IFN-γ locus in Africa

Despite its critical role in the defense against microbial infection and tumor development, little is known about the range of nucleotide and haplotype variation at IFN-γ, or the evolutionary forces that have shaped patterns of diversity at this locus. To address this gap in knowledge, we examined sequence data from the IFN-γ gene in 1461 individuals from 15 worldwide populations. Our analyses uncovered novel patterns of variation in distinct African populations, including an excess of high frequency-derived alleles, unusually long haplotype structure surrounding the IFN-γ gene, and a "star-like" genealogy of African-specific haplotypes carrying variants previously associated with infectious disease. We also inferred a deep time to coalescence of variation at IFN-γ (~ 0.8 million years ago) and ancient ages for common polymorphisms predating the evolution of modern humans. Taken together, these results are congruent with a model of positive selection on standing variation in African populations. Furthermore, we inferred that common variants in intron 3 of IFN-γ are the likely targets of selection. In addition, we observed a paucity of non-synonymous substitutions relative to synonymous changes in the exons of IFN-γ in African and non-African populations, suggestive of strong purifying selection. Therefore, we contend that positive and purifying selection have influenced levels of diversity in different regions of IFN-γ, implying that these distinct genic regions are, or have been, functionally important. Overall, this study provides additional insights into the evolutionary events that have contributed to the frequency and distribution of alleles having a role in human health and disease.

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.

Genetic structure of two erythrocyte binding antigens of Plasmodium falciparum reveals a contrasting pattern of selection

Erythrocyte binding antigens 175 (EBA-175) and 140 (EBA-140) play key roles in erythrocyte invasion by binding to glycophorin A (GPA) and C (GPC) respectively in human malaria. Since antigenic variation in malaria endemic region is a major barrier to development of effective vaccine, we explore the nature and pattern of sequence diversity of these two vaccine candidates in Kolkata, India. Population genetic parameters based on parasite sequences representing region II of Pfeba-175 and Pfeba-140 genes were estimated using DnaSP V.5.10 and MEGA version 6.0. A novel molecular docking approach was implemented to assess the binding affinities of Kolkata Pfeba-175 variants with GPA. P. falciparum Kolkata isolates experienced a recent population expansion as documented by negative Tajima's D, Fu & Li's statistics, unimodal mismatch distribution and star-like median-joining network for both loci. Positive selection seemed to play a major role in shaping the diversity of Pfeba-175 (d_N/d_S=2.45, and McDonald-Kreitman P-value=0.04) with successive accumulation of Q584K/E, E592A and R664S deriving high frequency haplotypes designated here as F2KH3 and F2KH1. In silico molecular docking demonstrated that polypeptides encoded by F2KH1 and F2KH3 were capable of engaging the parasite ligand into energetically favorable interaction with GPA. Our data demonstrated emergence of Pfeba-175 sequences harboring selectively advantageous nonsynonymous substitutions on Pf3D7 sequence background in the Kolkata parasite isolates. A contrasting pattern of Pf3D7-centric expansion of parasite sequences was noted for Pfeba-140. Together, this study provides a firm genetic and biological support favoring a dominant role of EBA-175 in erythrocyte invasion.

High rate of adaptation of mammalian proteins that interact with Plasmodium and related parasites

Plasmodium parasites, along with their Piroplasm relatives, have caused malaria-like illnesses in terrestrial mammals for millions of years. Several Plasmodium-protective alleles have recently evolved in human populations, but little is known about host adaptation to blood parasites over deeper evolutionary timescales. In this work, we analyze mammalian adaptation in ~500 Plasmodium- or Piroplasm- interacting proteins (PPIPs) manually curated from the scientific literature. We show that (i) PPIPs are enriched for both immune functions and pleiotropy with other pathogens, and (ii) the rate of adaptation across mammals is significantly elevated in PPIPs, compared to carefully matched control proteins. PPIPs with high pathogen pleiotropy show the strongest signatures of adaptation, but this pattern is fully explained by their immune enrichment. Several pieces of evidence suggest that blood parasites specifically have imposed selection on PPIPs. First, even non-immune PPIPs that lack interactions with other pathogens have adapted at twice the rate of matched controls. Second, PPIP adaptation is linked to high expression in the liver, a critical organ in the parasite life cycle. Finally, our detailed investigation of alpha-spectrin, a major red blood cell membrane protein, shows that domains with particularly high rates of adaptation are those known to interact specifically with P. falciparum. Overall, we show that host proteins that interact with Plasmodium and Piroplasm parasites have experienced elevated rates of adaptation across mammals, and provide evidence that some of this adaptation has likely been driven by blood parasites.

Glycophorin alleles link to malaria protection

Genome-wide association studies identify variants in the human glycophorin ABE locus

Resistance to malaria through structural variation of red blood cell invasion receptors

The malaria parasite Plasmodium falciparum invades human red blood cells by a series of interactions between host and parasite surface proteins. By analyzing genome sequence data from human populations, including 1269 individuals from sub-Saharan Africa, we identify a diverse array of large copy-number variants affecting the host invasion receptor genes GYPA and GYPB We find that a nearby association with severe malaria is explained by a complex structural rearrangement involving the loss of GYPB and gain of two GYPB-A hybrid genes, which encode a serologically distinct blood group antigen known as Dantu. This variant reduces the risk of severe malaria by 40% and has recently increased in frequency in parts of Kenya, yet it appears to be absent from west Africa. These findings link structural variation of red blood cell invasion receptors with natural resistance to severe malaria.

Inferences of African evolutionary history from genomic data

Africa is the origin of anatomically modern humans and a continent of linguistic, cultural, environmental, phenotypic, and genetic diversity. However, African populations remain underrepresented in genetic studies, which have largely focused on individuals with European and Asian ancestry. The expansion of high-throughput 'omic' technologies to interrogate multiple tissue types across many biomolecules-DNA, proteins, epigenetic modifications, metabolites, and others-has heralded a new era of investigation into African history. In this review, we summarize how some of these recent advances have been applied to contemporary sub-Saharan African populations to inform studies on human origins and adaptation.

Glycophorins, Blood Groups, and Protection from Severe Malaria

In Malaŵi, Malungo alibe odi is a saying that translates as: 'Malaria does not ask permission before coming in'. The recent finding of a new severe malaria resistance locus next to a cluster of glycophorin genes involved in Plasmodium falciparum erythrocyte invasion seems to suggest otherwise: that evolutionary pressure is enabling erythrocytes to lock the door to keep malaria out.

A novel locus of resistance to severe malaria in a region of ancient balancing selection

The high prevalence of sickle haemoglobin in Africa shows that malaria has been a major force for human evolutionary selection, but surprisingly few other polymorphisms have been proven to confer resistance to malaria in large epidemiological studies. To address this problem, we conducted a multi-centre genome-wide association study (GWAS) of life-threatening Plasmodium falciparum infection (severe malaria) in over 11,000 African children, with replication data in a further 14,000 individuals. Here we report a novel malaria resistance locus close to a cluster of genes encoding glycophorins that are receptors for erythrocyte invasion by P. falciparum. We identify a haplotype at this locus that provides 33% protection against severe malaria (odds ratio = 0.67, 95% confidence interval = 0.60-0.76, P value = 9.5 × 10(-11)) and is linked to polymorphisms that have previously been shown to have features of ancient balancing selection, on the basis of haplotype sharing between humans and chimpanzees. Taken together with previous observations on the malaria-protective role of blood group O, these data reveal that two of the strongest GWAS signals for severe malaria lie in or close to genes encoding the glycosylated surface coat of the erythrocyte cell membrane, both within regions of the genome where it appears that evolution has maintained diversity for millions of years. These findings provide new insights into the host-parasite interactions that are critical in determining the outcome of malaria infection.

Positive selection underlies the species-specific binding of Plasmodium falciparum RH5 to human basigin

Plasmodium falciparum, the causative agent of the deadliest form of malaria, is a member of the Laverania subgenus, which includes ape-infecting parasites. P. falciparum is thought to have originated in gorillas, although infection is now restricted to humans. Laverania parasites display remarkable host-specificity, which is partially mediated by the interaction between parasite ligands and host receptors. We analyse the evolution of BSG (basigin) and GYPA (glycophorin A) in primates/hominins, as well as of their Plasmodium-encoded ligands, PfRH5 and PfEBA175. We show that, in primates, positive selection targeted two sites in BSG (F27 and H102), both involved in PfRH5 binding. A population genetics-phylogenetics approach detected the strongest selection for the gorilla lineage: one of the positively selected sites (K191) is a major determinant of PfRH5 binding affinity. Analysis of RH5 genes indicated episodic selection on the P. falciparum branch; the positively selected W447 site is known to stabilize the interaction with human basigin. Conversely, we detect no selection in the receptor-binding region of EBA175 in the P. falciparum lineage. Its host receptor, GYPA, shows evidence of positive selection in all hominid lineages; selected codons include glycosylation sites that modulate PfEBA175 binding affinity. Data herein provide an evolutionary explanation for species-specific binding of the PfRH5-BSG ligand-receptor pair and support the hypothesis that positive selection at these genes drove the host shift leading to the emergence of P. falciparum as a human pathogen.

Blood Groups in Infection and Host Susceptibility

Blood group antigens represent polymorphic traits inherited among individuals and populations. At present, there are 34 recognized human blood groups and hundreds of individual blood group antigens and alleles. Differences in blood group antigen expression can increase or decrease host susceptibility to many infections. Blood groups can play a direct role in infection by serving as receptors and/or coreceptors for microorganisms, parasites, and viruses. In addition, many blood group antigens facilitate intracellular uptake, signal transduction, or adhesion through the organization of membrane microdomains. Several blood groups can modify the innate immune response to infection. Several distinct phenotypes associated with increased host resistance to malaria are overrepresented in populations living in areas where malaria is endemic, as a result of evolutionary pressures. Microorganisms can also stimulate antibodies against blood group antigens, including ABO, T, and Kell. Finally, there is a symbiotic relationship between blood group expression and maturation of the gastrointestinal microbiome.

Biased gene conversion skews allele frequencies in human populations, increasing the disease burden of recessive alleles

Gene conversion results in the nonreciprocal transfer of genetic information between two recombining sequences, and there is evidence that this process is biased toward G and C alleles. However, the strength of GC-biased gene conversion (gBGC) in human populations and its effects on hereditary disease have yet to be assessed on a genomic scale. Using high-coverage whole-genome sequences of African hunter-gatherers, agricultural populations, and primate outgroups, we quantified the effects of GC-biased gene conversion on population genomic data sets. We find that genetic distances (FST and population branch statistics) are modified by gBGC. In addition, the site frequency spectrum is left-shifted when ancestral alleles are favored by gBGC and right-shifted when derived alleles are favored by gBGC. Allele frequency shifts due to gBGC mimic the effects of natural selection. As expected, these effects are strongest in high-recombination regions of the human genome. By comparing the relative rates of fixation of unbiased and biased sites, the strength of gene conversion was estimated to be on the order of Nb ≈ 0.05 to 0.09. We also find that derived alleles favored by gBGC are much more likely to be homozygous than derived alleles at unbiased SNPs (+42.2% to 62.8%). This results in a curse of the converted, whereby gBGC causes substantial increases in hereditary disease risks. Taken together, our findings reveal that GC-biased gene conversion has important population genetic and public health implications.

Natural selection and infectious disease in human populations

Infectious pathogens are among the strongest selective forces that shape the human genome. Migrations and cultural changes in the past 100,000 years exposed populations to dangerous new pathogens.

Host genetics influences susceptibility to infectious disease. Evolutionary adaptations for resistance and symbiosis may underlie common immune-mediated diseases.

Signatures of selection and methods to detect them vary with the age, geographical spread and virulence of the pathogen.

A history of selection on a trait adds power to association studies by driving the emergence of common alleles of strong effect. Combining selection and association metrics can further increase power.

Genome-wide association studies (GWASs) of susceptibility to pathogens that are moderately old (1,000–50,000 years ago), geographically limited in history and exerted strong positive selective pressure will have the most power if GWASs can be done in the historically affected population.

An understanding of host–pathogen interactions can inform the development of new therapies for both infectious diseases and common immune-mediated diseases.

The impact of various infectious agents on human survival and reproduction over thousands of years has exerted selective pressure on numerous regions of the human genome. This Review describes how such signatures of selection can be detected and integrated with data from complementary approaches, such as genome-wide association studies, to provide biological insights into host–pathogen interactions.

The ancient biological 'arms race' between microbial pathogens and humans has shaped genetic variation in modern populations, and this has important implications for the growing field of medical genomics. As humans migrated throughout the world, populations encountered distinct pathogens, and natural selection increased the prevalence of alleles that are advantageous in the new ecosystems in both host and pathogens. This ancient history now influences human infectious disease susceptibility and microbiome homeostasis, and contributes to common diseases that show geographical disparities, such as autoimmune and metabolic disorders. Using new high-throughput technologies, analytical methods and expanding public data resources, the investigation of natural selection is leading to new insights into the function and dysfunction of human biology.

Modulation of glycan recognition by clustered saccharide patches

All cells in nature are covered with a dense and complex array of glycan chains. Specific recognition and binding of glycans is a critical aspect of cellular interactions, both within and between species. Glycan-protein interactions tend to be of low affinity but high specificity, typically utilizing multivalency to generate the affinity required for biologically relevant binding. This review focuses on a higher level of glycan organization, the formation of clustered saccharide patches (CSPs), which can constitute unique ligands for highly specific interactions. Due to technical challenges, this aspect of glycan recognition remains poorly understood. We present a wealth of evidence for CSPs-mediated interactions, and discuss recent advances in experimental tools that are beginning to provide new insights into the composition and organization of CSPs. The examples presented here are likely the tip of the iceberg, and much further work is needed to elucidate fully this higher level of glycan organization.

Using mutagenesis and structural biology to map the binding site for the Plasmodium falciparum merozoite protein PfRh4 on the human immune adherence receptor

To survive and replicate within the human host, malaria parasites must invade erythrocytes. Invasion can be mediated by the P. falciparum reticulocyte-binding homologue protein 4 (PfRh4) on the merozoite surface interacting with complement receptor type 1 (CR1, CD35) on the erythrocyte membrane. The PfRh4 attachment site lies within the three N-terminal complement control protein modules (CCPs 1-3) of CR1, which intriguingly also accommodate binding and regulatory sites for the key complement activation-specific proteolytic products, C3b and C4b. One of these regulatory activities is decay-accelerating activity. Although PfRh4 does not impact C3b/C4b binding, it does inhibit this convertase disassociating capability. Here, we have employed ELISA, co-immunoprecipitation, and surface plasmon resonance to demonstrate that CCP 1 contains all the critical residues for PfRh4 interaction. We fine mapped by homologous substitution mutagenesis the PfRh4-binding site on CCP 1 and visualized it with a solution structure of CCPs 1-3 derived by NMR and small angle x-ray scattering. We cross-validated these results by creating an artificial PfRh4-binding site through substitution of putative PfRh4-interacting residues from CCP 1 into their homologous positions within CCP 8; strikingly, this engineered binding site had an ∼30-fold higher affinity for PfRh4 than the native one in CCP 1. These experiments define a candidate site on CR1 by which P. falciparum merozoites gain access to human erythrocytes in a non-sialic acid-dependent pathway of merozoite invasion.

Patterns of nucleotide and haplotype diversity at ICAM-1 across global human populations with varying levels of malaria exposure

Malaria is one of the strongest selective pressures in recent human evolution. African populations have been and continue to be at risk for malarial infections. However, few studies have re-sequenced malaria susceptibility loci across geographically and genetically diverse groups in Africa. We examined nucleotide diversity at Intercellular adhesion molecule-1 (ICAM-1), a malaria susceptibility candidate locus, in a number of human populations with a specific focus on diverse African ethnic groups. We used tests of neutrality to assess whether natural selection has impacted this locus and tested whether SNP variation at ICAM-1 is correlated with malaria endemicity. We observe differing patterns of nucleotide and haplotype variation in global populations and higher levels of diversity in Africa. Although we do not observe a deviation from neutrality based on the allele frequency distribution, we do observe several alleles at ICAM-1, including the ICAM-1 (Kilifi) allele, that are correlated with malaria endemicity. We show that the ICAM-1 (Kilifi) allele, which is common in Africa and Asia, exists on distinct haplotype backgrounds and is likely to have arisen more recently in Asia. Our results suggest that correlation analyses of allele frequencies and malaria endemicity may be useful for identifying candidate functional variants that play a role in malaria resistance and susceptibility.

Multiple instances of ancient balancing selection shared between humans and chimpanzees

Instances in which natural selection maintains genetic variation in a population over millions of years are thought to be extremely rare. We conducted a genome-wide scan for long-lived balancing selection by looking for combinations of SNPs shared between humans and chimpanzees. In addition to the major histocompatibility complex, we identified 125 regions in which the same haplotypes are segregating in the two species, all but two of which are noncoding. In six cases, there is evidence for an ancestral polymorphism that persisted to the present in humans and chimpanzees. Regions with shared haplotypes are significantly enriched for membrane glycoproteins, and a similar trend is seen among shared coding polymorphisms. These findings indicate that ancient balancing selection has shaped human variation and point to genes involved in host-pathogen interactions as common targets.

Population genomic scan for candidate signatures of balancing selection to guide antigen characterization in malaria parasites

Acquired immunity in vertebrates maintains polymorphisms in endemic pathogens, leading to identifiable signatures of balancing selection. To comprehensively survey for genes under such selection in the human malaria parasite Plasmodium falciparum, we generated paired-end short-read sequences of parasites in clinical isolates from an endemic Gambian population, which were mapped to the 3D7 strain reference genome to yield high-quality genome-wide coding sequence data for 65 isolates. A minority of genes did not map reliably, including the hypervariable var, rifin, and stevor families, but 5,056 genes (90.9% of all in the genome) had >70% sequence coverage with minimum read depth of 5 for at least 50 isolates, of which 2,853 genes contained 3 or more single nucleotide polymorphisms (SNPs) for analysis of polymorphic site frequency spectra. Against an overall background of negatively skewed frequencies, as expected from historical population expansion combined with purifying selection, the outlying minority of genes with signatures indicating exceptionally intermediate frequencies were identified. Comparing genes with different stage-specificity, such signatures were most common in those with peak expression at the merozoite stage that invades erythrocytes. Members of clag, PfMC-2TM, surfin, and msp3-like gene families were highly represented, the strongest signature being in the msp3-like gene PF10_0355. Analysis of msp3-like transcripts in 45 clinical and 11 laboratory adapted isolates grown to merozoite-containing schizont stages revealed surprisingly low expression of PF10_0355. In diverse clonal parasite lines the protein product was expressed in a minority of mature schizonts (<1% in most lines and ∼10% in clone HB3), and eight sub-clones of HB3 cultured separately had an intermediate spectrum of positive frequencies (0.9 to 7.5%), indicating phase variable expression of this polymorphic antigen. This and other identified targets of balancing selection are now prioritized for functional study.

A novel ENU-mutation in ankyrin-1 disrupts malaria parasite maturation in red blood cells of mice

The blood stage of the plasmodium parasite life cycle is responsible for the clinical symptoms of malaria. Epidemiological studies have identified coincidental malarial endemicity and multiple red blood cell (RBC) disorders. Many RBC disorders result from mutations in genes encoding cytoskeletal proteins and these are associated with increased protection against malarial infections. However the mechanisms underpinning these genetic, host responses remain obscure. We have performed an N-ethyl-N-nitrosourea (ENU) mutagenesis screen and have identified a novel dominant (haploinsufficient) mutation in the Ank-1 gene (Ank1^MRI23420) of mice displaying hereditary spherocytosis (HS). Female mice, heterozygous for the Ank-1 mutation showed increased survival to infection by Plasmodium chabaudi adami DS with a concomitant 30% decrease in parasitemia compared to wild-type, isogenic mice (wt). A comparative in vivo red cell invasion and parasite growth assay showed a RBC-autonomous effect characterised by decreased proportion of infected heterozygous RBCs. Within approximately 6–8 hours post-invasion, TUNEL staining of intraerythrocytic parasites, showed a significant increase in dead parasites in heterozygotes. This was especially notable at the ring and trophozoite stages in the blood of infected heterozygous mutant mice compared to wt (p<0.05). We conclude that increased malaria resistance due to ankyrin-1 deficiency is caused by the intraerythrocytic death of P. chabaudi parasites.

Genome-wide comparison of African-ancestry populations from CARe and other cohorts reveals signals of natural selection

The study of recent natural selection in human populations has important applications to human history and medicine. Positive natural selection drives the increase in beneficial alleles and plays a role in explaining diversity across human populations. By discovering traits subject to positive selection, we can better understand the population level response to environmental pressures including infectious disease. Our study examines unusual population differentiation between three large data sets to detect natural selection. The populations examined, African Americans, Nigerians, and Gambians, are genetically close to one another (F(ST) < 0.01 for all pairs), allowing us to detect selection even with moderate changes in allele frequency. We also develop a tree-based method to pinpoint the population in which selection occurred, incorporating information across populations. Our genome-wide significant results corroborate loci previously reported to be under selection in Africans including HBB and CD36. At the HLA locus on chromosome 6, results suggest the existence of multiple, independent targets of population-specific selective pressure. In addition, we report a genome-wide significant (p = 1.36 × 10(-11)) signal of selection in the prostate stem cell antigen (PSCA) gene. The most significantly differentiated marker in our analysis, rs2920283, is highly differentiated in both Africa and East Asia and has prior genome-wide significant associations to bladder and gastric cancers.